Time-Accurate and highly-Stable Explicit operators for stiff differential equations

被引:16
|
作者
Bassenne, Maxime [1 ,2 ]
Fu, Lin [1 ]
Mani, Ali [1 ]
机构
[1] Stanford Univ, Ctr Turbulence Res, Stanford, CA 94305 USA
[2] Stanford Univ, Lab Artificial Intelligence Med & Biomed Phys, Stanford, CA 94305 USA
基金
芬兰科学院;
关键词
PDEs; Stiffness; Time-marching methods; Runge-Kutta; IMEX; DIRK; RUNGE-KUTTA METHODS; DISCONTINUOUS GALERKIN METHODS; REACTING FLOW; IMPLICIT; STABILITY; EXTRAPOLATION; INTEGRATION; SCHEMES; SYSTEMS;
D O I
10.1016/j.jcp.2020.109847
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
Unconditionally stable implicit time-marching methods are powerful in solving stiff differential equations efficiently. In this work, a novel framework to handle stiff physical terms implicitly is proposed. Both physical and numerical stiffness originating from convection, diffusion and source terms (typically related to reaction) can be handled by a set of predefined Time-Accurate and highly-Stable Explicit (TASE) operators in a unified framework. The proposed TASE operators act as preconditioners on the stiff terms and can be deployed to any existing explicit time-marching methods straightforwardly. The resulting time integration methods remain the original explicit time-marching schemes, yet with nearly unconditional stability. The TASE operators can be designed to be arbitrarily high-order accurate with Richardson extrapolation such that the accuracy order of original explicit time-marching method is preserved. Theoretical analyses and stability diagrams show that the s-stages sth-order explicit Runge-Kutta (RK) methods are unconditionally stable when preconditioned by the TASE operators with order p <= s and p <= 2. On the other hand, the sth-order RK methods preconditioned by the TASE operators with order of p <= s and p > 2 are nearly unconditionally stable. The only free parameter in TASE operators can be determined a priori based on stability arguments. Unlike classical implicit methods, the TASE methodology allows for solving non-linear problems with arbitrary order without requiring solving a nonlinear system of equations. A set of benchmark problems with strong stiffness is simulated to assess the performance of the TASE method. Numerical results suggest that the proposed framework preserves the high-order accuracy of the explicit time-marching methods with very-large time steps for all the considered cases. As an alternative to established implicit strategies, TASE method is promising for the efficient computation of stiff physical problems. (C) 2020 Elsevier Inc. All rights reserved.
引用
收藏
页数:24
相关论文
共 50 条
  • [1] A note on the stability of time-accurate and highly-stable explicit operators for stiff differential equations
    Calvo, M.
    Montijano, J., I
    Randez, L.
    [J]. JOURNAL OF COMPUTATIONAL PHYSICS, 2021, 436
  • [2] On a generalization of time-accurate and highly-stable explicit operators for stiff problems
    Aceto, Lidia
    Conte, Dajana
    Pagano, Giovanni
    [J]. APPLIED NUMERICAL MATHEMATICS, 2024, 200 : 2 - 17
  • [3] Time-accurate and highly-stable explicit peer methods for stiff differential problems
    Conte, Dajana
    Pagano, Giovanni
    Paternoster, Beatrice
    [J]. COMMUNICATIONS IN NONLINEAR SCIENCE AND NUMERICAL SIMULATION, 2023, 119
  • [4] A HIGHLY ACCURATE SOLVER FOR STIFF ORDINARY DIFFERENTIAL EQUATIONS
    Kushnir, Dan
    Rokhlin, Vladimir
    [J]. SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2012, 34 (03): : A1296 - A1315
  • [5] A-STABLE, ACCURATE AVERAGING OF MULTISTEP METHODS FOR STIFF DIFFERENTIAL EQUATIONS
    LINIGER, W
    ODEH, F
    [J]. IBM JOURNAL OF RESEARCH AND DEVELOPMENT, 1972, 16 (04) : 335 - &
  • [6] Explicit Methods for Stiff Stochastic Differential Equations
    Abdulle, Assyr
    [J]. NUMERICAL ANALYSIS OF MULTISCALE COMPUTATIONS, 2012, 82 : 1 - 22
  • [7] Singly TASE Operators for the Numerical Solution of Stiff Differential Equations by Explicit Runge–Kutta Schemes
    Manuel Calvo
    Lin Fu
    Juan I. Montijano
    Luis Rández
    [J]. Journal of Scientific Computing, 2023, 96
  • [8] Explicit methods in solving stiff ordinary differential equations
    Ahmad, RR
    Yaacob, N
    Murid, AHM
    [J]. INTERNATIONAL JOURNAL OF COMPUTER MATHEMATICS, 2004, 81 (11) : 1407 - 1415
  • [9] EXPLICIT STABILIZED MULTIRATE METHOD FOR STIFF DIFFERENTIAL EQUATIONS
    Abdulle, Assyr
    Grote, Marcus J.
    De Souza, Giacomo Rosilho
    [J]. MATHEMATICS OF COMPUTATION, 2022, 91 (338) : 2681 - 2714
  • [10] Composition explicit methods for stiff ordinary differential equations
    Medovikov, AA
    [J]. FINITE DIFFERENCE METHODS: THEORY AND APPLICATIONS, 1999, : 183 - 195