pith. machine review for the scientific record. sign in

arxiv: 2605.05645 · v1 · submitted 2026-05-07 · 🧮 math.NA · cs.NA

Recognition: unknown

Long-time stability of implicit-explicit Runge-Kutta methods for two-dimensional incompressible flows

Cao Wen, Hong-lin Liao, Xiaoming Wang, Xuping Wang

Pith reviewed 2026-05-08 06:55 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords implicit-explicit Runge-Kuttalong-time stabilityNavier-Stokes equationsvorticity-stream formulationGrönwall inequalityHölder inequalitymathematical inductionadaptive time-stepping
0
0 comments X

The pith

High-order implicit-explicit Runge-Kutta methods achieve long-time stability for two-dimensional incompressible Navier-Stokes flows.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper constructs a family of implicit-explicit Runge-Kutta methods up to fourth-order accuracy for the two-dimensional incompressible Navier-Stokes equations written in vorticity-stream function form. It proves that these schemes remain stable over arbitrarily long times in both the L2 and H1 norms by applying a convolution-type Hölder inequality together with a damping-type multistage Grönwall inequality. A mathematical induction step bounds the vorticity stage by stage in a fractional Sobolev norm. This framework allows the methods to use larger time steps in adaptive simulations without encountering instability, addressing a gap in the analysis of high-order schemes for fluid flows.

Core claim

The authors establish a unified analytical framework that proves long-time stability in both the L² and H¹ norms for a family of implicit-explicit Runge-Kutta methods applied to the two-dimensional incompressible Navier-Stokes equations. The schemes employ stiffly accurate diagonally implicit Runge-Kutta approximations for the linear diffusive term together with explicit Runge-Kutta approximations for the nonlinear advection term. By exploiting the model structure they reduce the order conditions to 5 for third-order and 11 for fourth-order methods. The proof combines a convolution-type Hölder inequality with a damping-type multistage Grönwall inequality and relies on a mathematical-indation

What carries the argument

Convolution-type Hölder inequality combined with damping-type multistage Grönwall inequality, supported by mathematical induction for stage-wise boundedness of vorticity in the H^δ norm.

Load-bearing premise

The mathematical-induction argument ensures stage-wise boundedness of the vorticity in the H^δ norm for the specific IERK schemes and the Navier-Stokes vorticity-stream formulation.

What would settle it

A concrete numerical simulation or analytical counterexample in which the vorticity fails to remain bounded in the H^δ norm for some initial data or time-step sequence, causing the L2 or H1 norms to grow without bound over long times.

Figures

Figures reproduced from arXiv: 2605.05645 by Cao Wen, Hong-lin Liao, Xiaoming Wang, Xuping Wang.

Figure 1
Figure 1. Figure 1: Errors of IERK(2,3;0.35), IERK(3,5;1.2) and IERK(4,7;-0.8) schemes. view at source ↗
Figure 2
Figure 2. Figure 2: Errors of IERK(2,3;c2), IERK(3,5;a55) and IERK(4,7;ˆa43). We employ the IERK(2,3;0.35), IERK(3,5;1.2) and IERK(4,7;-0.8) as the repre￾sentative second-, third-, and fourth-order methods, respectively. Numerical solutions are computed using time steps τ = 2−k/10 for 0 ≤ k ≤ 7, and the errors are measured until the final time T = 1 view at source ↗
Figure 3
Figure 3. Figure 3: Numerical simulation of Example 4.2 using the first-order scheme (1.3) specifically, for the first periodic region and a positive integer lt, we define (4.1) g2(t) := ( 0, t ∈ [0, T1], sin2 2πlt(t−T1) T −T1  , t ∈ [T1, T]. This forcing exhibits two distinct temporal regimes: a quiescent interval on [0, T1] and a highly oscillatory interval on [T1, T] when lt ≫ 1. In the regime of large viscosity (i.e., sm… view at source ↗
Figure 4
Figure 4. Figure 4: Simulation of Example 4.2 using the first-order scheme (1.3) with τmax = 0.1. Further numerical results for the case τmax = 0.1 are presented in view at source ↗
Figure 5
Figure 5. Figure 5: Numerical simulations of Example 4.2 using IERK(2,3;c2) schemes with ATS-LDLB strategy for τmax = 1 (left), τmax = 0.5 (middle), and τmax = 0.1 (right). Further inspections indicate that for the cases c2 = 1.4 and c2 = 4, the numerical errors near t = T1 exceed 1.0 regardless of the values of τmax (for clarity, the errors for c2 = 1.4 and c2 = 4 are only plotted for the case τmax = 1). From the perspective… view at source ↗
Figure 6
Figure 6. Figure 6: The mixed errors ∥e n mix∥∞ of IERK(3,5;a55) (top) and IERK(4,7;ˆa43) (bottom) schemes with the ATS-LDLB strategy for Example 4.2 view at source ↗
Figure 7
Figure 7. Figure 7: Performance of IERK(4,7;2) method with ATS and ATS-LDLB strategies view at source ↗
Figure 8
Figure 8. Figure 8: Performance of IERK(4,7;2) with ATS-LDLB strategy ( view at source ↗
read the original abstract

High-order adaptive time-stepping algorithms are of significant practical value and theoretical interest for accelerating long-time fluid-flow simulations and resolving complex dynamical behaviors. While several high-order implicit-explicit schemes have been proposed in the literature, their long-time stability properties remain largely unexplored. We develop a family of long-time stable implicit-explicit Runge-Kutta (IERK) methods, up to fourth-order temporal accuracy, for the two-dimensional incompressible Navier-Stokes equations in vorticity-stream function formulation. By combining a convolution-type H\"{o}lder inequality with a damping-type multistage Gr\"{o}nwall inequality, we establish a unified analytical framework that proves long-time stability in both the $L^2$ and $H^1$ norms. A key component of the analysis is a mathematical-induction argument that ensures stage-wise boundedness of the vorticity in the $H^\delta$ norm for some $\delta>0$. To the best of our knowledge, this is the first work to establish large-time stability results for high-order IERK algorithms for the two-dimensional incompressible Navier-Stokes equations. Our IERK schemes employ stiffly accurate diagonally implicit Runge-Kutta approximations for the linear diffusive term together with explicit Runge-Kutta approximations for the nonlinear advection term. By exploiting the specific structure of the Navier-Stokes model, we derive a reduced set of order conditions-requiring only 5 and 11 conditions for the third- and fourth-order methods, respectively, in contrast to the classical 6 and 18-allowing the construction of a parameterized family of efficient schemes. These IERK methods are particularly well suited for adaptive time-stepping, as they permit significantly enlarged step sizes in actual computations.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript develops a family of implicit-explicit Runge-Kutta (IERK) methods, up to fourth-order temporal accuracy, for the two-dimensional incompressible Navier-Stokes equations in vorticity-stream function formulation. These schemes employ stiffly accurate diagonally implicit Runge-Kutta approximations for the linear diffusive term and explicit Runge-Kutta approximations for the nonlinear advection term. The central claim is a unified analytical framework proving long-time stability in both L² and H¹ norms, obtained by combining a convolution-type Hölder inequality with a damping-type multistage Grönwall inequality, together with a mathematical-induction argument that ensures stage-wise boundedness of the vorticity in the H^δ norm for some δ>0. The methods are constructed via a reduced set of order conditions (5 for order 3, 11 for order 4) that exploit the structure of the Navier-Stokes model, and the work asserts this is the first such large-time stability result for high-order IERK algorithms on this problem.

Significance. If the stability analysis is rigorous, the result would be significant for numerical analysis of fluid flows: it supplies the first long-time L²/H¹ stability guarantees for high-order IERK schemes on 2D incompressible NS, directly supporting adaptive time-stepping with enlarged step sizes in long-time simulations. The framework relies on standard inequalities applied after an induction step rather than on machine-checked proofs or fully parameter-free derivations, but the reduced order conditions and explicit exploitation of the vorticity-stream structure are positive features that could be leveraged in follow-up work.

major comments (2)
  1. [Abstract / stability analysis] Abstract and stability analysis section: the mathematical-induction argument that keeps every internal stage value ω_i^n bounded in H^δ (δ>0) is load-bearing for the long-time claim, yet its closure is not verified in detail; specifically, the explicit advection term evaluated at the intermediate stages of the stiffly accurate DIRK/ERK pair must absorb commutator and embedding constants into the damping term without producing time-dependent growth that would prevent the multistage Grönwall factor from remaining uniform.
  2. [Order conditions and stability analysis] The reduced order conditions (5 for third-order, 11 for fourth-order) are stated to suffice for the schemes, but it is not shown that these conditions automatically guarantee the stage-order or stability-function properties needed for the nonlinear H^δ estimate to close under the induction; any failure here would undermine the subsequent application of the convolution Hölder and damping Grönwall inequalities.
minor comments (2)
  1. [Introduction] The abstract's phrasing 'to the best of our knowledge' should be supported by a concise literature comparison in the introduction that distinguishes the present long-time result from existing short-time or low-order IERK analyses for NS.
  2. [Preliminaries / notation] Notation for the stage values and the precise definition of the H^δ norm (including the value of δ) should be introduced earlier and used consistently when stating the induction hypothesis.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. The two major points raised concern the level of detail in the induction argument and the explicit connection between the reduced order conditions and the required stability properties. We address both below and will revise the manuscript to provide the requested clarifications and additional lemmas.

read point-by-point responses

  1. Referee: [Abstract / stability analysis] Abstract and stability analysis section: the mathematical-induction argument that keeps every internal stage value ω_i^n bounded in H^δ (δ>0) is load-bearing for the long-time claim, yet its closure is not verified in detail; specifically, the explicit advection term evaluated at the intermediate stages of the stiffly accurate DIRK/ERK pair must absorb commutator and embedding constants into the damping term without producing time-dependent growth that would prevent the multistage Grönwall factor from remaining uniform.

    Authors: We agree that the closure of the induction step requires more explicit verification. In the revised manuscript we will expand the stability analysis section with a dedicated paragraph (or short subsection) that applies the induction hypothesis stage by stage. We will show that the stiffly accurate property of the DIRK part supplies a uniform damping factor that absorbs the commutator terms arising from the vorticity-stream formulation together with the embedding constants from the Hölder inequality. Because the damping coefficient is independent of the time step and of the stage index, the multistage Grönwall factor remains bounded uniformly in time, exactly as claimed. A remark on the admissible range for δ > 0 will also be added. revision: yes

  2. Referee: [Order conditions and stability analysis] The reduced order conditions (5 for third-order, 11 for fourth-order) are stated to suffice for the schemes, but it is not shown that these conditions automatically guarantee the stage-order or stability-function properties needed for the nonlinear H^δ estimate to close under the induction; any failure here would undermine the subsequent application of the convolution Hölder and damping Grönwall inequalities.

    Authors: The referee is correct that an explicit link between the reduced order conditions and the stage-order/stability-function properties used in the H^δ estimate was not spelled out. While the conditions were obtained by exploiting the divergence-free structure and the specific form of the advection term, we acknowledge that a direct verification is needed. In the revision we will insert a short lemma immediately after the statement of the order conditions. The lemma will verify that the reduced Butcher tableaux satisfy the necessary stage-order conditions and that the stability functions remain bounded in a manner compatible with the induction hypothesis, thereby allowing the convolution Hölder and multistage Grönwall arguments to close without additional growth factors. revision: yes

Circularity Check

0 steps flagged

Stability proof self-contained with external inequalities and induction

full rationale

The derivation chain rests on applying a convolution-type Hölder inequality and a damping-type multistage Grönwall inequality after a preliminary mathematical induction that secures stage-wise H^δ vorticity bounds. These are standard external analytic tools whose validity does not presuppose the target long-time L²/H¹ stability result; the reduced order conditions (5 for order 3, 11 for order 4) are obtained directly from the stiffly accurate DIRK/ERK structure of the NS vorticity-stream system rather than from the stability claim itself. No self-definitional loops, fitted parameters renamed as predictions, or load-bearing self-citations appear in the provided derivation steps.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on two standard inequalities and an induction argument whose closure depends on the specific IERK coefficients and the NS vorticity equation structure.

axioms (2)
  • standard math Convolution-type Hölder inequality
    Invoked to bound nonlinear terms in the stability analysis.
  • standard math Damping-type multistage Grönwall inequality
    Used to obtain long-time bounds from the multistage formulation.

pith-pipeline@v0.9.0 · 5620 in / 1280 out tokens · 39464 ms · 2026-05-08T06:55:06.940078+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

48 extracted references · 1 canonical work pages

  1. [1]

    Archilla and J

    B. Archilla and J. Novo , Robust error bounds for the Navier-Stokes equations using implicit-explicit second order BDF method with variable steps. IMA J. Numer. Anal. , 43:5 (2023), pp. 2892--2933

    2023

  2. [2]

    Ascher, S

    U. Ascher, S. Ruuth and R.J. Spiteri , Implicit-explicit Runge-Kutta methods for time dependent partial differential equations, Appl. Numer. Math. , 25:2 3 (1997), pp. 151--167

    1997

  3. [3]

    Boscarino, L

    S. Boscarino, L. Pareschi and G. Russo , A unified IMEX Runge-Kutta approach for hyperbolic systems with multiscale relaxation, SIAM J. Numer. Anal. , 55:4 (2017), pp. 2085--2109

    2017

  4. [4]

    Boisneault, S

    A. Boisneault, S. Dubuis and M. Picasso , An adaptive space-time algorithm for the incompressible Navier-Stokes equations, J. Comput. Phys. , 493 (2023), 112457

    2023

  5. [5]

    Boscarino and G

    S. Boscarino and G. Russo , On a class of uniformly accurate IMEX Runge-Kutta schemes and applications to hyperbolic systems with relaxation, SIAM J. Sci. Comput. , 31:3 (2009), pp. 1926--1945

    2009

  6. [6]

    Calvo, J

    M.P. Calvo, J. De Frutos and J. Novo , Linearly implicit Runge-Kutta methods for advection-reaction-diffusion equations, Appl. Numer. Math. , 37:4 (2001), pp. 535--549

    2001

  7. [7]

    Cardone, Z

    A. Cardone, Z. Jackiewicz, A. Sandu and H. Zhang , Extrapolated implicit-explicit Runge-Kutta methods, Math. Model. Anal. , 19:1 (2014), pp. 18--43

    2014

  8. [8]

    W. Chen, M. Gunzburger, D. Sun, X.-M. Wang , Efficient and long-time accurate second-order methods for the Stokes-Darcy system, SIAM J. Numer. Anal. , 51:5 (2013), pp. 2563--2584

    2013

  9. [9]

    W. Chen, M. Gunzburger, D. Sun, X.-M. Wang , An efficient and long-time accurate third-order algorithm for the Stokes--Darcy system , Numer. Math. , 134:4 (2016), pp. 857--879

    2016

  10. [10]

    Y. Chen, Y. Huang, H.-L. Liao and N. Yi , Energy behaviors of implicit-explicit Runge-Kutta methods in adaptive finite element simulations for the Cahn-Hilliard model, Commun. Comput. Phys. , 2026, to appear

    2026

  11. [11]

    Cheng and X.-M

    W. Cheng and X.-M. Wang , A semi-implicit scheme for stationary statistical properties of the infinite Prandtl number model, SIAM J. Numer. Anal. , 47:1 (2008), pp. 250--270

    2008

  12. [12]

    Cheng, C

    K. Cheng, C. Wang , Long time stability of high order multistep numerical schemes for two-dimensional incompressible Navier-Stokes equations, SIAM J. Numer. Anal. , 54:5 (2016), pp. 3123--3144

    2016

  13. [13]

    Cheng, C

    K. Cheng, C. Wang, S. Wise and X. Yue , A second-order, weakly energy-stable pseudo-spectral scheme for the Cahn-Hilliard equation and its solution by the homogeneous linear iteration method, J. Sci. Comput. , 69:3 (2016), pp. 1083--1114

    2016

  14. [14]

    T. Chu, J. Wang, N. Wang and Z. Zhang , Optimal-order convergence of a two-step BDF method for the Navier-Stokes equations with H^1 initial data, J. Sci. Comput. , 96:2 (2023), 62

    2023

  15. [15]

    Constantin and C

    P. Constantin and C. Foias , Navier-Stokes equations , University of Chicago Press, New York, 1988

    1988

  16. [16]

    DeCaria and M

    V. DeCaria and M. Schneier , An embedded variable step IMEX scheme for the incompressible Navier-Stokes equations, Comput. Methods Appl. Mech. Eng. , 376 (2021), 113661

    2021

  17. [17]

    Dimarco and L

    G. Dimarco and L. Pareschi , Asymptotic preserving implicit-explicit Runge-Kutta methods for nonlinear kinetic equations, SIAM J. Numer. Anal. , 51:2 (2013), pp. 1064--1087

    2013

  18. [18]

    Foias, O

    C. Foias, O. Manley, R. Rosa, and R. Temam , Navier-Stokes equations and turbulence , Cambridge University Press, Cambridge, 2001

    2001

  19. [19]

    Foias and R

    C. Foias and R. Temam , Gevrey class regularity for the solutions of the Navier-Stokes equations, J. Funct. Anal. , 87:2 (1989), pp. 359--369

    1989

  20. [20]

    Z. Fu, T. Tang and J. Yang , Energy diminishing implicit-explicit Runge-Kutta methods for gradient flows, Math. Comput. , 93:350 (2024), pp. 2745--2767

    2024

  21. [21]

    Gottlieb, F

    S. Gottlieb, F. Tone, C. Wang, X.-M. Wang, D. Wirosoetisno , Long time stability of a classical efficient scheme for two-dimensional Navier--Stokes equations, SIAM J. Numer. Anal. , 50:1 (2012), pp. 126--150

    2012

  22. [22]

    Gottlieb and C

    S. Gottlieb and C. Wang , Stability and convergence analysis of fully discrete Fourier collocation spectral method for 3-D viscous Burgers' equation, J. Sci. Comput. , 53:2 (2012), pp. 102--128

    2012

  23. [23]

    Guesmi, M

    M. Guesmi, M. Grotteschi and J. Stiller , Assessment of high-order IMEX methods for incompressible flow, Int. J. Numer. Meth. Fluids , 95:6 (2023), pp. 954--978

    2023

  24. [24]

    Izzo and Z

    G. Izzo and Z. Jackiewicz , Highly stable implicit-explicit Runge-Kutta methods, Appl. Numer. Math. , 113 (2017), pp. 71-92

    2017

  25. [25]

    Kanevsky, M.H

    A. Kanevsky, M.H. Carpenter, D. Gottlieb and J.S. Hesthaven , Application of implicit-explicit high order Runge-Kutta methods to discontinuous-Galerkin schemes, J. Comput. Phys. , 225:2 (2007), pp. 1753--1781

    2007

  26. [26]

    A. T. Hill and E. S\" u li , Approximation of the global attractor for the incompressible Navier-Stokes equations, IMA J. Numer. Anal. , 20:4 (2000), pp. 633--667

    2000

  27. [27]

    Horn and C.R

    R.A. Horn and C.R. Johnson , Matrix Analysis, Second Edition, Cambridge University Press, New York, 2013

    2013

  28. [28]

    Kennedy and M.H

    C.A. Kennedy and M.H. Carpenter , Additive Runge-Kutta schemes for convection-diffusion-reaction equations, Appl. Numer. Math. , 44:1 2 (2003), pp. 139--181

    2003

  29. [29]

    Li and H.-L

    Z. Li and H.-L. Liao , Stability of variable-step BDF2 and BDF3 methods , SIAM J. Numer. Anal. , 60:4 (2022), pp. 2253--2272

    2022

  30. [30]

    H.-L. Liao, T. Tang, X.-P. Wang and T. Zhou , A class of refined implicit-explicit Runge-Kutta methods with robust time adaptability and unconditional convergence for the Cahn-Hilliard model , Math. Comput. , 2025, doi: 10.1090/mcom/4090

    work page doi:10.1090/mcom/4090 2025

  31. [31]

    Liao, X.-P

    H.-L. Liao, X.-P. Wang and C. Wen , Original energy dissipation preserving corrections of integrating factor Runge-Kutta methods for gradient flow problems , J. Comput. Phys. , 2024, 519: 113456

    2024

  32. [32]

    Liao, X.-P

    H.-L. Liao, X.-P. Wang and C. Wen , Average energy dissipation rates of additive implicit-explicit Runge-Kutta methods for gradient flow problems, CSIAM Trans. Appl. Math. , 2026, to appear

    2026

  33. [33]

    Liao and Z

    H.-L. Liao and Z. Zhang , Analysis of adaptive BDF2 scheme for diffusion equations , Math. Comput. , 90:329 (2021), pp. 1207--1226

    2021

  34. [34]

    Liu and J

    H. Liu and J. Zou , Some new additive Runge-Kutta methods and their applications , J. Comput. Appl. Math. , 190:1 2 (2006), pp. 74--98

    2006

  35. [35]

    Liu and C.-W

    J.-G. Liu and C.-W. Shu , A high-order discontinuous Galerkin method for 2D incompressible flows, J. Comput. Phys. , 160:2 (2000), pp. 577--596

    2000

  36. [36]

    N. C. Nguyen, J. Peraire and B. Cockburn , An implicit high-order hybridizable discontinuous Galerkin method for the incompressible Navier-Stokes equations, J. Comput. Phys. , 230:4 (2011), pp. 1147--1170

    2011

  37. [37]

    Pareschi and G

    L. Pareschi and G. Russo , Implicit-explicit Runge-Kutta schemes and applications to hyperbolic systems with relaxation, J. Sci. Comput. , 25:1 2 (2005), pp. 129--155

    2005

  38. [38]

    J. Shen, T. Tang and L. Wang , Spectral methods: Algorithms, analysis and applications, Springer-Verlag, Berlin Heidelberg, 2011

    2011

  39. [39]

    A. M. Stuart and A. R. Humphries , Dynamical systems and numerical analysis , Cambridge University Press, New York, 1998

    1998

  40. [40]

    Tachim-Medjo , Navier-Stokes equations in the vorticity-velocity formulation: The two-dimensional case, Appl

    T. Tachim-Medjo , Navier-Stokes equations in the vorticity-velocity formulation: The two-dimensional case, Appl. Numer. Math. , 21:2 (1996), pp. 185--206

    1996

  41. [41]

    Temam , Infinite-dimensional dynamical systems in mechanics and physics , 2nd ed., Springer-Verlag, New York, 1997

    R. Temam , Infinite-dimensional dynamical systems in mechanics and physics , 2nd ed., Springer-Verlag, New York, 1997

    1997

  42. [42]

    Tone, X.-M

    F. Tone, X.-M. Wang, and D. Wirosoetisno , Long-time dynamics of 2d double-diffusive convection: analysis and/of numerics. Numer. Math. , 130:3 (2015), pp. 541--566

    2015

  43. [43]

    Tone and D

    F. Tone and D. Wirosoetisno , On the long-time stability of the implicit Euler scheme for the two-dimensional Navier-Stokes equations, SIAM J. Numer. Anal. , 44:1 (2006), pp. 29--40

    2006

  44. [44]

    X.-M. Wang, Z. Wang and M. Mao , Linearly implicit variable step-size BDF schemes with Fourier pseudospectral approximation for incompressible Navier-Stokes equations, Appl. Numer. Math. , 172 (2022), pp. 393--412

    2022

  45. [45]

    Wang , Approximation of stationary statistical properties of dissipative dynamical systems: time discretization, Math

    X.-M. Wang , Approximation of stationary statistical properties of dissipative dynamical systems: time discretization, Math. Comput. , 79:269 (2010), pp. 259--280

    2010

  46. [46]

    Wang , An efficient second order in time scheme for approximating long time statistical properties of the two dimensional Navier-Stokes equations, Numer

    X.-M. Wang , An efficient second order in time scheme for approximating long time statistical properties of the two dimensional Navier-Stokes equations, Numer. Math. , 121:4 (2012), pp. 753--779

    2012

  47. [47]

    Wang , Numerical algorithms for stationary statistical properties of dissipative dynamical systems, Dis

    X.-M. Wang , Numerical algorithms for stationary statistical properties of dissipative dynamical systems, Dis. Contin. Dyn. Sys. , 36:8 (2016), pp. 4599--4618

    2016

  48. [48]

    X.-P. Wang, X. Zhao and H.-L. Liao , A unified framework on the original energy laws of three effective classes of Runge-Kutta methods for phase field crystal type models, SIAM J. Numer. Anal. , 63:4 (2025), pp. 1808--1832

    2025