Premium
An effective numerical integration method for typical stiff systems
Author(s) -
Aiken Richard C.,
Lapidus Leon
Publication year - 1974
Publication title -
aiche journal
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.958
H-Index - 167
eISSN - 1547-5905
pISSN - 0001-1541
DOI - 10.1002/aic.690200225
Subject(s) - singular perturbation , ordinary differential equation , method of matched asymptotic expansions , perturbation (astronomy) , mathematics , stiff equation , boundary layer , equivalence (formal languages) , mathematical analysis , stability (learning theory) , boundary (topology) , feature (linguistics) , numerical analysis , boundary value problem , differential equation , control theory (sociology) , mechanics , computer science , physics , control (management) , linguistics , philosophy , discrete mathematics , quantum mechanics , machine learning , artificial intelligence
Abstract There is an equivalence between stiff and singularly perturbed systems of ordinary differential equations. This feature is exploited in this paper by numerically employing recent singular perturbation methods to attack troublesome boundary layer stage of the solution in which some variables have very short response times. The numerical method affords a means of essentially determining the thickness of this boundary layer. The algorithm is capable of high stability and accuracy for the commonly occurring stiff system, whether or not it is in singularly perturbed form. Application to a singularly perturbed reaction system and a highly stiff reactor system not in singularly perturbed form demonstrate the effectiveness and utility of this approach.