System Of Nonlinear Equations Research Articles

Comparison of different elastic strain definitions for largely deformed SEI of chemo-mechanically coupled silicon battery particles

Amorphous silicon is a highly promising anode material for next-generation lithium-ion batteries. Large volume changes of the silicon particle have a critical effect on the surrounding solid-electrolyte interphase (SEI) due to repeated fracture and healing during cycling. Based on a thermodynamically consistent chemo-elasto-plastic continuum model we investigate the stress development inside the particle and the SEI. Using the example of a particle with SEI, we apply a higher order finite element method together with a variable-step, variable-order time integration scheme on a nonlinear system of partial differential equations. Starting from a single silicon particle setting, the surrounding SEI is added in a first step with the typically used elastic Green–St-Venant (GSV) strain definition for a purely elastic deformation. For this type of deformation, the definition of the elastic strain is crucial to get reasonable simulation results. In case of the elastic GSV strain, the simulation aborts. We overcome the simulation failure by using the definition of the logarithmic Hencky strain. However, the particle remains unaffected by the elastic strain definitions in the particle domain. Compared to GSV, plastic deformation with the Hencky strain is straightforward to take into account. For the plastic SEI deformation, a rate-independent and a rate-dependent plastic deformation are newly introduced and numerically compared for three half cycles for the example of a radial symmetric particle.

Effect of MHD and radiation on biviscous Bingham fluid flow on Marangoni boundary for heat source/sink with chemical reaction

A significance of Marangoni-driven convective magnetohydrodynamics (MHD) was investigated in a planned study. The present study examines the impact of radiation and chemical reactions on the two dimensional boundary layer flow of bi-viscous Bingham fluid on a thermosolutal Marangoni boundary with magnetic field and heat source/sink. The physical flow problem is mathematically modelled into Navier–Stokes equations. These nonlinear partial differential equations are then mapped into a set of nonlinear ordinary differential equations using similarity transformation. The elegant analytical method is used to obtain analytic approximations for the resulting system of nonlinear differential equations. The analytical solution for the temperature and concentration profiles is expressed in terms of Kummer's confluent hypergeometric function. The features of flow characteristics such as velocity, temperature, and concentration profiles in response to the variations of the emerging parameters are simulated and examined with a physical explanation through graphs. Surface tension is considered to be corresponding to heat and mass. The results reveal that the velocity profile decreases with increasing the magnetic field, while the Marangoni number value increases with increasing the velocity profile. Moreover, the temperature profile rises with rising the radiation parameter and whereas the value chemical reaction parameter upsurges with decline concentration profile. The present problem has great interest due to its applications in industrial applications such as drying of silicon wafers, thin layers of paint, glues, in heat exchangers, crystal growth in space, biological fluids, blood, synovial liquids, plasma, lubricants, many pharmaceutical products, proteins, sewage sludge, china clay, and many more.

High performance computational approach to study model describing reversible two-step enzymatic reaction with time fractional derivative

Enzyme reactions have numerous applications in diverse disciplines of science like chemistry, biology and biomechanics. In this study, we examine the role and act of enzymes in chemical reactions which is considered in the frame of fractional order model. The proposed model includes system of four equations which are studied via Caputo fractional operator. The systems of non-linear equations are evaluated by a semi-analytical approach called q\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$q$$\\end{document}-homotopy analysis transform method. The uniqueness and existence of the solutions has been investigated through fixed point theorem. The solutions of the proposed model are achieved through the considered method and the obtained outcomes are in the form of series which shows rapid convergence. The solutions are computed and graphs are plotted for the obtained results using mathematica software. The achieved results by the proposed method are unique and illustrate the significant dynamics of the considered model via 3D plots and graphs. The results of this study demonstrate the importance and effectiveness of projected derivative and technique in the analysis of time dependent fractional mathematical models. This study also gives an idea to extend the applications of enzymatic reactions in drug development, bio mechanics, and chemical reactions in various cellular metabolisms. Also, enzymatic reactions have a vital role in the fields of the food industry for processing food, in biotechnology for the manufacture of biofuels, and in metabolic engineering to design metabolic pathways.

Signal and image reconstruction with a double parameter Hager–Zhang‐type conjugate gradient method for system of nonlinear equations

AbstractThe one parameter conjugate gradient method by Hager and Zhang (Pac J Optim, 2(1):35–58, 2006) represents a family of descent iterative methods for solving large‐scale minimization problems. The nonnegative parameter of the scheme determines the weight of conjugacy and descent, and by extension, the numerical performance of the method. The scheme, however, does not converge globally for general nonlinear functions, and when the parameter approaches 0, the scheme reduces to the conjugate gradient method by Hestenes and Stiefel (J Res Nat Bur Stand, 49:409–436, 1952), which in practical sense does not perform well due to the jamming phenomenon. By carrying out eigenvalue analysis of an adaptive two parameter Hager–Zhang type method, a new scheme is presented for system of monotone nonlinear equations with its application in compressed sensing. The proposed scheme was inspired by nice attributes of the Hager–Zhang method and the various schemes designed with double parameters. The scheme is also applicable to nonsmooth nonlinear problems. Using fundamental assumptions, analysis of the global convergence of the scheme is conducted and preliminary report of numerical experiments carried out with the scheme and some recent methods indicate that the scheme is promising.

Influence of Rotation and Viscosity on Parallel Rolls of Electrically Conducting Fluid

Rayleigh–Bénard convection is a fundamental fluid dynamics phenomenon that significantly influences heat transfer in various natural and industrial processes, such as geophysical dynamics in the Earth’s liquid core and the performance of heat exchangers. Understanding the behavior of conductive fluids under the influence of heating, rotation, and magnetic fields is critical for improving thermal management systems. Utilizing the Boussinesq approximation, this study theoretically examines the nonlinear convection of a planar layer of conductive liquid that is heated from below and subjected to rotation about a vertical axis in the presence of a magnetic field. We focus on the onset of stationary convection as the temperature difference applied across the planar layer increases. Our theoretical approach investigates the formation of parallel rolls aligned with the magnetic field under free–free boundary conditions. To analyze the system of nonlinear equations, we expand the dependent variables in a series of orthogonal functions and express the coefficients of these functions as power series in a parameter ϵ. A solution for this nonlinear problem is derived through Fourier analysis of perturbations, extending to O(ϵ8), which allows for a detailed visualization of the parallel rolls. Graphical results are presented to explore the dependence of the Nusselt number on the Rayleigh number (R) and Ekman number (E). We observe that both the local Nusselt number and average Nusselt number increase as the Ekman number decreases. Furthermore, the flow appears to become more deformed as E decreases, suggesting an increased influence of external factors such as rotation. This deformation may enhance mixing within the fluid, thereby improving heat transfer between different regions.

Tripled Fixed Point Methodologies for Finding a Solution to a System of Nonlinear Integral Equations on Orthogonal Metric Spaces

Tripled Fixed Point Methodologies for Finding a Solution to a System of Nonlinear Integral Equations on Orthogonal Metric Spaces

Multimodal vortex-induced vibration mitigation and design approach of bistable nonlinear energy sink inerter on bridge structure

Large-scale structures, e.g., long-span bridge structures, are prone to induce multi-modal vibrations due to their densely spaced low modal frequencies. Due to the limited frequency bandwidth of linear dynamic absorbers, they are incapable of effectively mitigating vibrations across multiple modes. To this end, the bistable nonlinear energy sink inerter (BNESI) is used to mitigate the multimodal vortex-induced vibration (VIV) of the beam structure. The highly nonlinear equilibrium differential equations of the beam-BNESI system are numerically solved, and the simulated annealing (SA) algorithm is employed to determine the optimal VIV reduction ratio and BNESI parameters. In comparison to the cubic-type nonlinear energy sink inerter (CNESI), BNESI is found to possess more stable equilibrium positions, smaller stiffness coefficients, and higher VIV mitigation efficiency. The selection of design modes has been found to influence the efficiency of multimodal VIV mitigation, with the use of the intermediate modal order as the design mode resulting in the highest efficiency for multimodal VIV mitigation. The performance-based multimodal VIV mitigation design can be realized with three parameters, i.e., inertance ratio, damping coefficient, and stiffness coefficient. Moreover, the performance-based multimodal VIV mitigation approach and models proposed in this study demonstrate a high level of precision.

Stability Analysis of Dual Solutions for Nonlinear Radiative Magnetohydrodynamic Flow of [formula omitted] Hybrid Nanofluid over a Nonlinearly Shrinking Surface

The present work explores the existence and temporal stability of dual solutions with the consequence of nonlinear thermal radiation on hybrid nanofluid flow (including Silver and Titanium dioxide nanoparticles in water) past a permeable shrinking sheet. The analysis also considers the united outcomes of suction/injection, Joule heating, viscous dissipation, and magnetohydrodynamics. Employing realistic assumptions and suitable similarity transformations, the governing nonlinear partial differential equations are formulated and altered into a nonlinear ordinary differential equations system. These equations are then solved numerically using the Lobatto IIIA technique. It is observed that dual solutions can occur within a precise range of the shrinking surface parameter. Analyzing the temporal stability of these dual solutions under small disturbances shows that the upper solution branch is stable and represents a physically realistic solution to the problem. In contrast, the lower solution branch is unstable and not physically realizable. No solution exists beyond the critical λc=−1.3915 for the shrinking parameter. The critical values of shrinking parameter λc=−1.2565,−1.3188, and −1.4161 correspond to volume fractions of hybrid nanofluid with ϕAg−TiO2 of 2%,4%, and 8%, respectively. No solution exists beyond this critical point. The graphical representation and quantitative explanation demonstrate how various emergent characteristics affect momentum and thermal boundary layer profiles, Nusselt number, and skin friction. The velocity and temperature profile of hybrid nanofluid can be improved by adding a small volume of nanoparticle volume fraction. Thermal boundary layer thickness is amplified by an accumulation in the shrinking parameter, power index of velocity component, thermal radiation parameter, and temperature difference ratio. Sheering stress and heat transfer coefficient improve with adding more nanoparticles in the base fluid for the first solution braches, but opposite effects are found in the second solutions. The contributions of this research are significant both theoretically, in terms of advancing the mathematical modeling of hybrid nanofluid flow with heat transfer in engineering systems, and practically, in terms of applications to engineering cooling systems.

Enhanced heat transfer in novel star-shaped enclosure with hybrid nanofluids: A neural network-assisted study

The current research is a numerical investigation of the thermo-fluidic transport trend within a new star-shaped enclosure filled with Fe3O4−MWCNT (Multiwall Carbon Nanotubes) suspended hybrid nanoparticles in a base fluid (water), induced by a horizontal magnetizing field. A rectangular rod inside the cavity positioned at center is kept at a high temperature, as the low temperature is chosen for outer corrugated boundary. Formulation yelled in mathematical principles, along with the boundary conditions and some physical assumptions, are formulated as dimensionless PDEs. The finite element method (FEM) assisted by “COMSOL” software, employed for numerical simulations, and the PARDISO solver is used for solving nonlinear system of equations. Also, an artificial neural network (ANN) model is assembled to analyze the effect of these parameters on thermal and flow transport properties of Hybrid nanofluid. The training of this ANN model is done by a dataset generated from numerical simulations, suggesting predictions regarding heat transport under varying parameters. This approach does not only reduce the effort, but also reduces computing time when exploring the thermal behaviors across various datasets. The visualization of velocity and temperature fields through streamlines and isotherms, along with the evaluation of the Nusselt number, illustrate the enhanced heat transfer facilitated by the incorporation of Fe3O4−MWCNT nanoparticles. The results show that increased magnetic field strength reduces fluid velocity due to the generated Lorentz force, which counteracts convective flow. Percentage analysis indicates that hybrid nanofluids (Fe3O4−MWCNT) significantly improve thermal distribution compared to conventional fluids. This study demonstrates the effectiveness of ANN models in forecasting the influence of diverse factors on the flow and heat transport characteristics of hybrid nanofluids. These insights are essential to the design and optimization of heat transfer systems.

On a priori estimate of periodic solutions of a system of nonlinear ordinary differential equations of the second order

The question of a priori estimation of periodic solutions for a system of non-linear ordinary differential equations of the second order with a distinguished main positively homogeneous part is investigated. In this question, the well-known methods for deriving an a priori estimation of periodic solutions for similar systems of first-order ordinary differential equations are not directly applicable.Combining these methods with the idea of qualitative research of singularly perturbed ordinary differential equations, conditions are found that provide an a priori estimation of periodic solutions for the system of second-order equations.The conditions for the a priori estimation are formulated in terms of the properties of the main positively homogeneous part of the system of equations. The existence of periodic solutions is proved to be invariant under a continuous change of the main positively homogeneous part and the conditions of a priori estimation are preserved. Based on the results obtained, in the future, it is possible to investigate the existence of periodic solutions.

A new shifted generalized Chebyshev approach for multi-dimensional sinh-Gordon equation

The numerical treatment of multi-dimensional non-linear sinh-Gordon equations is the focus of this paper. We numerically solve the (1 + 1) and (2 + 1) sinh-Gordon equations using two collocation algorithms. We select the set of basis functions as a set of generalized Chebyshev polynomials (CPs), which we express as orthogonal combinations of CPs. We develop and utilize some formulas related to these polynomials to propose our numerical algorithms. Specific values for the high-order derivatives of the basis functions serve in the derivation of the two presented algorithms. Additionally, we provide an estimation of the basis functions used in the convergence analysis study. We follow the two collocation algorithms to transform the sinh-Gordon equations into non-linear equation systems, which any suitable solver can handle. We provide some examples and comparisons to confirm the effectiveness of our presented algorithms.

An intriguing numerical strategy for Zakharov–Kuznetsov equation through graph-theoretic polynomials

This paper explores graph-theoretic polynomials to find the approximate solution of the (2+1)D Time-fractional Zakharov-Kuznetsov(TF-Z-K) equation. The Zakharov-Kuznetsov equations govern the behavior of nonlinear acoustic waves in the plasma of hot isothermal electrons and cold ions in the presence of a homogeneous magnetic field. Independence polynomials of the Ladder-Rung graph serve as the polynomial approximation for the suggested Independence Polynomial Collocation Method (IPCM). The Caputo fractional derivatives are adopted to determine the fractional derivatives in the TF-Z-K equation. The TF-Z-K equation is converted into a system of nonlinear algebraic equations using the collocation points in IPCM. The Newton-Raphson approach yields the solution of the suggested method by solving the resulting system. We’ve compared a few scenarios with the tangible outcomes to validate the procedure. Quantitative outcomes match the current findings and validate the exactness of IPCM compared t o the recent numerical and semi-analytical approaches in the literature.

A compact combination of second-kind Chebyshev polynomials for Robin boundary value problems and Bratu-type equations

AbstractThis paper presents a new algorithm for resolving linear and non-linear second-order Robin boundary value problems (BVPS) and the Bratu-type equations in one and two dimensions using spectral approaches. Basis functions according to second-kind shifted and modified shifted Chebyshev polynomials that comply with the Robin conditions are created. It has produced operational matrices for its derivatives. The provided solutions are the result of applying the collocation and tau approaches. These methods convert the problem dictated by its boundary conditions into a system of linear or non-linear algebraic equations that may be solved using any suitable numerical solver. Convergence analysis has been provided and it accords with the numerical results. Six numerical problems are provided to investigate and demonstrate the practical utility of the suggested method. The current results show that our method outperforms the previous methods in terms of accuracy which are presented in tables and figures.

A least-change secant algorithm for solving generalized complementarity problem

AbstractIn this paper, we propose a least-change secant algorithm to solve the generalized complementarity problem indirectly trough its reformulation as a nonsmooth system of nonlinear equations using a one-parametric family of complementarity functions. We present local and superlinear convergence results of new algorithm and analyze its numerical performance.

A numerical study of heat and mass transfer characteristic of three-dimensional thermally radiated bi-directional slip flow over a permeable stretching surface

Within fluid mechanics, the flow of hybrid nanofluids over a stretching surface has been extensively researched due to their influence on the flow and heat transfer properties. Expanding on this concept by introducing porous media, the current study explore the flow and heat and mass transport characteristics of hybrid nanofluid. This investigation includes the effect of magnetohydrodynamic (MHD) with chemical reaction, thermal radiation, and slip effects. The nanoparticles, copper, and alumina are combined with water for the formation of a hybrid nanofluid. Using the self-similar method for the reduction of Partial differential equations (PDEs) to the system of Ordinary differential equations (ODEs). These nonlinear equation systems are solved numerically using the bvp4c (boundary value solver) technique. The effect of the different physical non-dimensional flow parameters on different flow profiles such as velocity, temperature, concentration, skin friction, Nusselt and mass transfer rate are depicted through graphs and tables. The velocity profiles diminish with the effect of magnetic and slip parameters. The temperature and concentration slip parameters reduce the temperature and concentration profile respectively. The higher values of magnetic factor lessened the skin friction coefficient for both slip and no-slip conditions. An elevation in the thermal slip parameter reduced the boundary layer thickness and the heat transfer from the surface to the fluid. The Nusselt number amplified with the climbing values of the radiation parameter. The mass transfer rate depressed with the solutal slip parameter. Comparison is made with the published work in the literature and there is excellent agreement between them.

On solvability of a two-dimensional symmetric nonlinear system of difference equations

We show that the system of difference equations xn+k=xn+lyn−efxn+l+yn−e−f,yn+k=yn+lxn−efyn+l+xn−e−f,n∈N0,\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ x_{n+k}=\\frac{x_{n+l}y_{n}-ef}{x_{n+l}+y_{n}-e-f},\\quad y_{n+k}= \\frac{y_{n+l}x_{n}-ef}{y_{n+l}+x_{n}-e-f},\\quad n\\in {\\mathbb{N}}_{0}, $$\\end{document} where k∈N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$k\\in {\\mathbb{N}}$\\end{document}, l∈N0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$l\\in {\\mathbb{N}}_{0}$\\end{document}, l<k\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$l< k$\\end{document}, e,f∈C\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$e, f\\in {\\mathbb{C}}$\\end{document}, and xj,yj∈C\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$x_{j}, y_{j}\\in {\\mathbb{C}}$\\end{document}, j=0,k−1‾\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$j=\\overline{0,k-1}$\\end{document}, is theoretically solvable and present some cases of the system when the general solutions can be found in a closed form.

A new solution to force analysis including Coulomb friction in mechanism joints

The objective of this paper is to propose a novel approximate solution for determining the reactions of joints in mechanical systems, which involves the influence of Coulomb friction. It is widely acknowledged that if Coulomb friction is used in determining an exact solution to the equilibrium equations for a mechanism, then it will involve nonlinear systems of equations. With the abundance of computer tools now available, tasks of this kind, especially numerical computations, are not particularly challenging. However, new analytically tractable approximate methods are still valuable as a straightforward way of solving the problem. Several studies have been carried out in this field to find a simple solution. In this paper, a new approach based on friction circle concept and Babylonian algorithm is developed for various mechanisms, which is exceptionally well-suited for calculating and streamlining the solution process for mechanical systems by eliminating the necessity for iterative steps at each stage of the force analysis.

Exact cosmological solutions for a Chaplygin Gas in anisotropic petrov type D spacetimes in Eddington-inspired-Born-Infeld gravity: Dark Energy model

We consider a Petrov Type D physical metric g, an auxiliary metric q and a Chaplygin Gas of pressure P in Eddington-inspired-Born-Infeld theory. From the Eddington-inspired-Born-Infeld-Chaplygin Gas equations, we first derive a system of second order nonlinear ordinary differential equations. Then, by a suitable change of variables, we arrive at a system of first order linear ordinary differential equations for the non-vanishing components of the pressure P, the physical metric g and the auxiliary metric q. Thanks to the superposition method, we collect an analytical solution for the nonlinear system obtained, which allows to retrieve new exact cosmological solutions for the model considered. By studying the Kretschmann invariant, we see that a singularity exists at the origin of the cosmic time. By the Kruskal-like coordinates, we conclude that this solution is the counterpart of the Friedman-Lemaître-Robertson-Walker spacetime in the Eddington-inspired-Born-Infeld theory. The Hubble and deceleration parameters in both directions of the physical metric g and the auxiliary metric q, as well as their behaviours over time, are also studied. The thermodynamic behaviour of the Chaplygin Gas model is investigated and, as a result, we show that the third-law of thermodynamics is verified. This means that the value of the entropy of the Chaplygin Gas in the perfect crystal state is zero at a temperature of zero Kelvin, which yields a determined value of the entropy and not an additive constant. Finally, we show that the solutions change asymptotically to the isotropic regime of expansion of Dark Energy. With this, we infer that the Chaplygin Gas can show a unified picture of Dark Energy and Dark Matter cooling during the expansion of the Universe.

Vibration suppression of a platform by a fractional type electromagnetic damper and inerter-based nonlinear energy sink

The theory of linear and nonlinear dynamic vibration absorbers is a well-established topic for many years. However, many recent contributions paid attention to the nonlinear vibration absorbers and different practical realizations of corresponding devices. Here, we propose a mechanical system constituted of the inerter-based nonlinear energy sink attached to the main body that is resting on an elastic foundation and is grounded through the fractional type electromagnetic damper. The two-degree-of-freedom system is described via two coupled differential equations with one of them having a fractional-order derivative term and the other one containing cubic stiffness nonlinearity. The incremental harmonic balance (IHB) method is employed to solve the equations and studies the strongly nonlinear periodic responses of the system. Applied approximated solution methodology is validated by the numerical Newmark method adapted to deal with the system of nonlinear fractional-order differential equations. The appropriate and necessary number of harmonics used in the IHB solution is commented and validated. This study can be a first step in understanding the dynamics and giving directions for the future design of vibration-isolating platforms based on inerter-based nonlinear vibration absorbers and electromagnetic dampers.

A spectral collocation scheme for solving nonlinear delay distributed-order fractional equations

The basic purpose of this paper is to implement a variant of the Jacobi–Gauss collocation scheme for solving nonlinear delay (in difference format) fractional differential equations of distributed-order type. At the first stage of the solving procedure, the Legendre–Gauss quadrature rule is implemented to approximate the main problem by a difference delay fractional differential equation of multi-term type. So, one can apply the Jacobi–Gauss collocation approach for localizing the resulting multi-term fractional differential equation with the difference delay factor. Moreover, some special cases of Jacobi–Gauss quadrature rules can be used for approximating the involved integrals in the aforementioned problem. By considering the Lagrange interpolating polynomials, as the numerical solutions, the solution of the main problem can be approximated by solving the associated system of nonlinear algebraic equations in terms of unknown Lagrange multipliers. Convergence analysis of the problem is investigated rigorously under some mild conditions at a spectral rate. Extensive test problems are considered to justify the theoretical analysis.