Abstract
For solving a system of nonlinear algebraic equations (NAEs) of the type: F(x) = 0, or Fi(xj) = 0; i; j = 1; : : : ;n, a Newton-like algorithm has several drawbacks such as local convergence, being sensitive to the initial guess of solution, and the time-penalty involved in finding the inversion of the Jacobian matrix ∂Fi=∂ xj. Based-on an invariant manifold defined in the space of (x; t) in terms of the residual-norm of the vector F(x), we can derive a gradient-flow system of nonlinear ordinary differential equations (ODEs) governing the evolution of x with a fictitious time-like variable t as an independent variable. We can prove that in the present novel Residual-Norm Based Algorithms (RNBAs), the residual-error is automatically decreased to zero along the path of x(t). More importantly, we have derived three iterative algoritms which do not involve the fictitious time and its stepsize Dt. We apply the three RNBAs to several numerical examples, revealing exponential convergences with different slopes and displaying the high efficiencies and accuracies of the present iterative algorithms. All the three presently proposed RNBAs: (i) are easy to implement numerically, (ii) converge much faster than the Newton's method, (iii) do not involve the inversion of the Jacobian ∂Fi=∂ xj, (iv) are suitable for solving a large system of NAEs, and (v) are purely iterative in nature.
Original language | English |
---|---|
Pages (from-to) | 279-304 |
Number of pages | 26 |
Journal | CMES - Computer Modeling in Engineering and Sciences |
Volume | 71 |
Issue number | 3 |
State | Published - 2011 |
Keywords
- Fictitious time integration method (FTIM)
- Iterative algorithm
- Non-linear ordinary differential equations
- Non-linear partial differential equations
- Nonlinear algebraic equations
- Residual-Norm Based Algorithms (RNBAs)