Improved iteration scheme for partial equilibrium flow (Q1083974): Difference between revisions
From MaRDI portal
Set profile property. |
ReferenceBot (talk | contribs) Changed an Item |
||
Property / cites work | |||
Property / cites work: Partial chemical equilibrium in fluid dynamics / rank | |||
Normal rank | |||
Property / cites work | |||
Property / cites work: Numerical method for partial equilibrium flow / rank | |||
Normal rank | |||
Property / cites work | |||
Property / cites work: Q5652137 / rank | |||
Normal rank |
Latest revision as of 16:34, 17 June 2024
scientific article
Language | Label | Description | Also known as |
---|---|---|---|
English | Improved iteration scheme for partial equilibrium flow |
scientific article |
Statements
Improved iteration scheme for partial equilibrium flow (English)
0 references
1985
0 references
A numerical method for partial equilibrium flow was recently described [the first author and \textit{L. D. Cloutman}, ibid. 39, 405-417 (1981; Zbl 0469.76068)]. The heart of the method is a quadratic iteration scheme in which the progress increment for equilibrium reaction s is obtained from the equilibrium constraint for reaction s alone. The scheme therefore converges best when the equilibrium reactions are weakly coupled, which can frequently be achieved by a judicious choice of the independent equilibrium reactions. The new iteration scheme consists of the following ingredients: (a) preconditioning of the equilibrium constraints conditions to make them more nearly linear in the progress variables, (b) application of a one-step Gauss-Seidel-Newton iteration to the preconditioned system, followed by (c) switching to a full Newton-Raphson iteration if the simpler Gauss-Seidel-Newton iteration fails to converge in a specified number of steps. If the equilibrium reactions are weakly coupled, convergence usually occurs before the full Newton-Raphson iteration is called into play. However, if the reactions are not weakly coupled the interaction between them is properly accounted for by the matrix inversion in the Newton-Raphson procedure.
0 references
weakly coupled equilibrium
0 references
partial equilibrium flow
0 references
quadratic iteration scheme
0 references
equilibrium reaction
0 references
equilibrium constraint
0 references
preconditioning of the equilibrium constraints conditions
0 references
one-step Gauss-Seidel-Newton iteration
0 references
preconditioned system
0 references
full Newton-Raphson iteration
0 references
convergence
0 references
matrix inversion
0 references