TY - JOUR
T1 - Design of order-preserving algorithms for transient first-order systems with controllable numerical dissipation
AU - Masuri, S. U.
AU - Sellier, M.
AU - Zhou, X.
AU - Tamma, K. K.
PY - 2011/12/30
Y1 - 2011/12/30
N2 - Using a new design procedure termed as Algorithms by Design, which we have successfully introduced in our previous efforts for second-order systems, alternatively, we advance in this exposition, the design and development of a computational framework that permits order-preserving second-order time accurate, unconditionally stable, zero-order overshooting behavior, and features with controllable numerical dissipation and dispersion via a family of algorithms for effectively solving transient first-order systems. The key feature is the incorporation of a spurious root to introduce controllable numerical dissipation while preserving second-order accuracy (order-preserving feature) resulting in a two-root system, namely, the principal root (ρ1∞) and a spurious root (ρ2∞). In contrast to the classical Trapezoidal family of algorithms which are the most popular, the present framework has the same order of computational complexity, but a higher payoff that is a significant advance to the field for tackling a wide class of applications dealing with first-order transient systems. We also present the special case with selection of ρ1∞ = 1 and any ρ2∞ leading to the design of a family of generalized single-step single-solve [GS4-1] algorithms recovering the Crank-Nicolson method at one end (ρ2∞ = 1) and the Midpoint Rule at the other end (ρ2∞ = 0) and anything in between, all of which have spectral radius features resembling that of the Crank-Nicolson method. More interestingly, with the particular choice of ρ1∞ = ρ2∞ = 0, the developed framework additionally inherits L-stable features. We illustrate the successful design of the developed GS4-1 framework using two simple illustrative numerical examples.
AB - Using a new design procedure termed as Algorithms by Design, which we have successfully introduced in our previous efforts for second-order systems, alternatively, we advance in this exposition, the design and development of a computational framework that permits order-preserving second-order time accurate, unconditionally stable, zero-order overshooting behavior, and features with controllable numerical dissipation and dispersion via a family of algorithms for effectively solving transient first-order systems. The key feature is the incorporation of a spurious root to introduce controllable numerical dissipation while preserving second-order accuracy (order-preserving feature) resulting in a two-root system, namely, the principal root (ρ1∞) and a spurious root (ρ2∞). In contrast to the classical Trapezoidal family of algorithms which are the most popular, the present framework has the same order of computational complexity, but a higher payoff that is a significant advance to the field for tackling a wide class of applications dealing with first-order transient systems. We also present the special case with selection of ρ1∞ = 1 and any ρ2∞ leading to the design of a family of generalized single-step single-solve [GS4-1] algorithms recovering the Crank-Nicolson method at one end (ρ2∞ = 1) and the Midpoint Rule at the other end (ρ2∞ = 0) and anything in between, all of which have spectral radius features resembling that of the Crank-Nicolson method. More interestingly, with the particular choice of ρ1∞ = ρ2∞ = 0, the developed framework additionally inherits L-stable features. We illustrate the successful design of the developed GS4-1 framework using two simple illustrative numerical examples.
KW - Controllable numerical dissipation
KW - First-order systems
KW - Order preserving
KW - Time integration
UR - http://www.scopus.com/inward/record.url?scp=81955160854&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=81955160854&partnerID=8YFLogxK
U2 - 10.1002/nme.3228
DO - 10.1002/nme.3228
M3 - Article
AN - SCOPUS:81955160854
SN - 0029-5981
VL - 88
SP - 1411
EP - 1448
JO - International Journal for Numerical Methods in Engineering
JF - International Journal for Numerical Methods in Engineering
IS - 13
ER -