Towards polyalgorithmic linear system solvers for nonlinear elliptic problems

Journal Article · Sun May 01 04:00:00 EDT 1994 · SIAM Journal on Scientific and Statistical Computing (Society for Industrial and Applied Mathematics); (United States)

DOI:https://doi.org/10.1137/0915044· OSTI ID:7112077

Ern, A ^[1]; Giovangigli, V ^[2]; Keyes, D E ^[3]; Smooke, M D ^[4]

---

1. Yale Univ., New Haven, CT (United States). Dept. of Mechanical Engineering Ecole Nationale des Ponts et Chausses, Paris (France)
2. Ecole Polytechnique, Palaiseau (France)
3. Institute for Computer Applications in Science and Engineering, Hampton, VA (United States)
4. Yale Univ., New Haven, CT (United States)

The authors investigate the performance of several preconditioned conjugate gradient-like algorithms and a standard stationary iterative method (block-line successive overrelaxation (SOR)) on linear systems of equations that arise from a nonlinear elliptic flame sheet problem simulation. The nonlinearity forces a pseudotransient continuation process that makes the problem parabolic and thus compacts the spectrum of the Jacobian matrix so that simple relaxation methods are viable in the initial stages of the solution process. However, because of the transition from parabolic to elliptic character as the timestep is increased in pursuit of the steady-state solution, the performance of the candidate linear solvers spreads as the domain of convergence of Newton's method is approached. In numerical experiments over the course of a full nonlinear solution trajectory, short recurrence or optimal Krylov algorithms combined with a Gauss-Seidel (GS) preconditioning yield better execution times with respect to the standard block-line SOR techniques, but SOR performs competitively at a smaller storage cost until the final stages. Block-incomplete factorization preconditioned methods, on the other hand, require nearly a factor of two more storage than SOR and are uniformly less effective during the pseudotransient stages. The advantage of GS preconditioning is partly attributable to the exploitation of a dominant convection direction in the examples; nevertheless, a multidomain version of GS with streamwise coupling lagged at rows between adjacent subdomains incurs only a modest penalty.

• DOE Contract Number: FG02-88ER13966
• OSTI ID: 7112077
• Journal Information: SIAM Journal on Scientific and Statistical Computing (Society for Industrial and Applied Mathematics); (United States), Journal Name: SIAM Journal on Scientific and Statistical Computing (Society for Industrial and Applied Mathematics); (United States) Vol. 15:3; ISSN 0196-5204; ISSN SIJCD4
• Country of Publication: United States
• Language: English

Related Subjects

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CHEMISTRY
400800 -- Combustion
Pyrolysis
& High-Temperature Chemistry
99 GENERAL AND MISCELLANEOUS
990200* -- Mathematics & Computers
ALGORITHMS
BOUNDARY-VALUE PROBLEMS
CALCULATION METHODS
CHEMICAL REACTIONS
COMBUSTION
COMPUTERIZED SIMULATION
DIFFERENTIAL EQUATIONS
EQUATIONS
FLAME PROPAGATION
ITERATIVE METHODS
MATHEMATICAL LOGIC
NONLINEAR PROBLEMS
OXIDATION
PARTIAL DIFFERENTIAL EQUATIONS
PERFORMANCE
SIMULATION
THERMOCHEMICAL PROCESSES