pith. sign in

arxiv: 1907.04254 · v1 · pith:HAG32RN6new · submitted 2019-07-09 · 🧮 math.NA · cs.NA

Arbitrarily High-order Unconditionally Energy Stable Schemes for Gradient Flow Models Using the Scalar Auxiliary Variable Approach

Yuezheng Gong , Jia Zhao , Qi Wang This is my paper

Pith reviewed 2026-05-25 00:06 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords gradient flowscalar auxiliary variableenergy stable schemeshigh-order methodsunconditional stabilitynumerical schemesFourier pseudospectral method
0
0 comments X

The pith

High-order scalar auxiliary variable schemes achieve arbitrary temporal accuracy for gradient flow models while remaining unconditionally energy stable.

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

The paper introduces the high-order scalar auxiliary variable (HSAV) method to construct numerical schemes for gradient flow models that attain arbitrarily high order in time. These schemes preserve the energy dissipation law with no restriction on the size of the time step. The approach is presented as general, independent of the specific form of the free energy, so it applies across a class of thermodynamically consistent models. Spatial discretization via the Fourier pseudospectral method yields fully discrete schemes that inherit the same unconditional stability. Numerical experiments on standard models confirm the expected accuracy orders and demonstrate that larger time steps suffice compared with earlier SAV schemes.

Core claim

The HSAV schemes reach arbitrarily high order in time while preserving the energy dissipation law without any restriction on the time step size, and the strategy applies to a class of thermodynamically consistent gradient flow models to produce semi-discrete high-order energy-stable schemes; the fully discrete versions with Fourier pseudospectral spatial discretization are likewise unconditionally energy stable.

What carries the argument

The high-order scalar auxiliary variable (HSAV) reformulation, which introduces a scalar auxiliary variable to rewrite the energy so that high-order time discretizations can be built while inheriting unconditional energy stability from the underlying SAV framework.

If this is right

  • The method yields semi-discrete schemes that remain unconditionally energy stable for any time step while attaining any prescribed temporal order.
  • The same unconditional stability carries over to the fully discrete schemes once the Fourier pseudospectral method is used in space.
  • Because the construction does not depend on the explicit form of the free energy, the same HSAV procedure produces high-order stable schemes for any gradient flow model that admits an SAV reformulation.
  • Numerical experiments show that the schemes attain their design order and permit substantially larger time steps than standard SAV schemes to reach a given accuracy.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The unconditional stability for arbitrary order could allow long-time simulations of phase separation or crystal growth at time steps set by accuracy rather than stability constraints.
  • If the SAV reformulation can be extended to other spatial discretizations, the HSAV approach might produce high-order stable schemes on unstructured meshes or finite-element spaces.
  • The generality claim suggests the method could be applied to coupled or nonlocal gradient flows once an appropriate scalar auxiliary variable is identified.

Load-bearing premise

Gradient flow models admit a scalar auxiliary variable reformulation of the energy that permits construction of high-order time discretizations while retaining the unconditional energy dissipation property.

What would settle it

A concrete numerical test on one of the standard gradient flow models (Allen-Cahn, Cahn-Hilliard, or similar) in which an HSAV scheme of claimed order p produces an energy increase for some time step size larger than the stability limit of lower-order methods.

read the original abstract

In this paper, we propose a novel family of high-order numerical schemes for the gradient flow models based on the scalar auxiliary variable (SAV) approach, which is named the high-order scalar auxiliary variable (HSAV) method. The newly proposed schemes could be shown to reach arbitrarily high order in time while preserving the energy dissipation law without any restriction on the time step size (i.e., unconditionally energy stable). The HSAV strategy is rather general that it does not depend on the specific expression of the effective free energy, such that it applies to a class of thermodynamically consistent gradient flow models arriving at semi-discrete high-order energy-stable schemes. We then employ the Fourier pseudospectral method for spatial discretization. The fully discrete schemes are also shown to be unconditionally energy stable. Furthermore, we present several numerical experiments on several widely-used gradient flow models, to demonstrate the accuracy, efficiency and unconditionally energy stability of the HSAV schemes. The numerical results verify that the HSAV schemes can reach the expected order of accuracy, and it allows a much larger time step size to reach the same accuracy than the standard SAV schemes.

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

0 major / 3 minor

Summary. The paper proposes the high-order scalar auxiliary variable (HSAV) method as an extension of the standard SAV reformulation to construct arbitrarily high-order time discretizations for gradient flow models. These schemes are claimed to be unconditionally energy stable (preserving the exact energy dissipation law for any time step size), independent of the specific form of the free energy, and applicable to a class of thermodynamically consistent models. The approach is combined with Fourier pseudospectral spatial discretization to yield fully discrete unconditionally stable schemes, with numerical experiments on standard models (e.g., Allen-Cahn, Cahn-Hilliard) verifying the expected temporal orders and improved efficiency over standard SAV.

Significance. If the central claims on arbitrary-order accuracy and unconditional stability hold with the stated generality, the work provides a broadly applicable framework for efficient, structure-preserving simulations of gradient flows. The independence from the concrete free-energy expression and the extension to both semi-discrete and fully discrete settings are potentially valuable for long-time integration of phase-field and related models.

minor comments (3)
  1. [Abstract] Abstract: the phrasing 'could be shown to reach arbitrarily high order' is tentative; revise to reflect that the manuscript contains the proofs and constructions.
  2. The distinction between the standard SAV energy identity and the new HSAV reformulation that enables arbitrary order should be made explicit with a dedicated equation or subsection early in the paper.
  3. Numerical experiments section: include a direct comparison table of maximum stable time steps between HSAV and standard SAV for the same accuracy tolerance to quantify the efficiency gain.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of the HSAV approach, its claimed generality, and the recommendation of minor revision. No specific major comments appear in the provided report, so we have no individual points to address.

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper constructs the HSAV schemes explicitly to enforce the energy dissipation law via the auxiliary variable reformulation and time discretization choices, which is the intended design rather than a reduction of a claimed prediction to its inputs. The abstract states the generality (no dependence on specific free-energy form) and unconditional stability as properties of the construction, with numerical experiments serving as verification rather than fitted inputs renamed as predictions. No self-citation chains, uniqueness theorems imported from prior author work, or ansatzes smuggled via citation are load-bearing in the provided description; the derivation remains self-contained against the stated SAV framework without internal redefinition or statistical forcing.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The approach relies on the standard SAV reformulation and domain assumptions typical for gradient flow models; no free parameters, new entities, or ad-hoc axioms are introduced beyond the existing framework.

axioms (2)
  • domain assumption Gradient flow models are thermodynamically consistent with a free energy functional bounded from below.
    Invoked to ensure the energy dissipation law holds and to support the SAV reformulation.
  • domain assumption The Fourier pseudospectral spatial discretization preserves the discrete energy stability property.
    Stated for the fully discrete schemes in the abstract.

pith-pipeline@v0.9.0 · 5733 in / 1288 out tokens · 26442 ms · 2026-05-25T00:06:36.588431+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

  • IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel echoes
    ?
    echoes

    ECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.

    The HSAV strategy is rather general that it does not depend on the specific expression of the effective free energy... preserving the energy dissipation law without any restriction on the time step size (i.e., unconditionally energy stable).

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

51 extracted references · 51 canonical work pages · 2 internal anchors

  1. [1]

    S. M. Allen and J. W. Cahn. Ground state structures in ordered binary alloys with second neighbor interactions. Acta Metall., 20(423), 1972

    work page 1972

  2. [2]

    V. E. Badalassi, H. D. Ceniceros, and S. Banerjee. Computation of multiphase systems with phase field models. Journal of Computational Physics, 190(2):371– 397, 2003

    work page 2003

  3. [3]

    Boyer and C

    F. Boyer and C. Lapuerta. Study of a three component cahn-hilliard flow model. ESAIM: Mathematical Modeling and Numerical Analysis, 40:653–687, 2006

    work page 2006

  4. [4]

    J. W. Cahn and J. E. Hilliard. Free energy of a nonuniform system. I. interfacial free energy. Journal of Chemical Physics, 28:258–267, 1958

    work page 1958

  5. [5]

    Chen and M

    J. Chen and M. Qin. Multi-symplectic Fourier pseudospectral method for the nonlinear Schr¨ odinger equation. Electronic Transactions on Numerical Analysis, 12:193–204, 2001

    work page 2001

  6. [6]

    L. Chen, J. Zhao, and X. Yang. Regularized linear schemes for the molecular beam epitaxy model with slope selection. Applied Numerical Mathematics, 128:138–156, 2018

    work page 2018

  7. [7]

    L. Q. Chen. Phase-field models for microstructure evolution. Annual Review of Material Research, 32(1):113–140, 2002

    work page 2002

  8. [8]

    W. Chen, S. Conde, C. Wang, X. Wang, and S. Wise. A linear energy stable scheme for a thin film model without slope selection. Journal of Scientific Computing, 52:546–562, 2012

    work page 2012

  9. [9]

    W. Chen, C. Wang, X. Wang, and S. M. Wise. A linear iteration algorithm for a second-order energy stable scheme for a thin film model without slope selection. Journal of Scientific Computing, 59(3):574–601, 2014

    work page 2014

  10. [10]

    D. J. Eyre. Unconditionally gradient stable time marching the Cahn-Hilliard equa- tion. Mrs Online Proceedings Library Archive, 529:39–46, 1998

    work page 1998

  11. [11]

    W. Feng, C. Wang, and S. Wise. Linearly preconditioned nonlinear conjugate gradeint solvers for the epitaxial thin film equation with slope selection. arXiv, pages 1–15, 2017

    work page 2017

  12. [12]

    W. Feng, C. Wang, S. Wise, and Z. Zhang. A second-order energy stable back- ward differentiation formula method for the epitaxial thin film equation with slope selection. ArXiv, page 1706.01943, 2017. 19

    work page internal anchor Pith review Pith/arXiv arXiv 2017

  13. [13]

    Gomez and T

    H. Gomez and T. J. R. Hughes. Provably unconditionally stable, second-order time- accurate, mixed variational methods for phase-field models. Journal of Computa- tional Physics, 230(13):5310–5327, 2011

    work page 2011

  14. [14]

    Y. Gong, J. Cai, and Y. Wang. Multi-symplectic Fourier pseudospectral method for the Kawahara equation. Communications in Computational Physics, 16(1):35–55, 2014

    work page 2014

  15. [15]

    Gong and J

    Y. Gong and J. Zhao. Energy-stable runge–kutta schemes for gradient flow models using the energy quadratization approach. Applied Mathematics Letters, 94:224– 231, 2019

    work page 2019

  16. [16]

    Y. Gong, J. Zhao, and Q. Wang. Linear second order in time energy stable schemes for hydrodynamic models of binary mixtures based on a spatially pseudospectral approximation. Advances in Computational Mathematics, 44:1573–1600, 2018

    work page 2018

  17. [17]

    Y. Gong, J. Zhao, and Q. Wang. Second order fully discrete energy stable methods on staggered grids for hydrodynamic phase field models of binary viscous fluids. SIAM J. Sci. Comput., 40(2):B528–B553, 2018

    work page 2018

  18. [18]

    Y. Gong, J. Zhao, X. Yang, and Q. Wang. Second-order linear schemes for hydro- dynamic phase field models of viscous fluid flows with variable densities. SIAM Journal on Scientific Computing, 4(1):B138–B167, 2018

    work page 2018

  19. [19]

    Gonzalez-Ferreiro, H

    B. Gonzalez-Ferreiro, H. Gomez, and I. Romero. A thermodynamically consistent numerical method for a phase field model of solidification. Commun Nonlinear Sci Numer Simulat, 19:2309–2323, 2014

    work page 2014

  20. [20]

    Guillen-Gonzailez and G

    F. Guillen-Gonzailez and G. Tierra. On linear schemes for a Cahn-Hilliard diffuse interface model. Journal of Computational Physics, 234:140–171, 2013

    work page 2013

  21. [21]

    Guo and Y

    R. Guo and Y. Xu. Semi-implicit spectral deferred correction method based on the invariant energy quadratization approach for phase field problems. Communication in Computational Physics, 26(1):87–113, 2019

    work page 2019

  22. [22]

    Z. Guo, P. Lin, J. Lowengrub, and S. Wise. Mass conservative and energy stable finite difference methods for the quasi-incompressible navier-stokes-cahn-hilliard system: primitive variable and projection-type schemes. Computer Methods in Applied Mechanics and Engineering, 326:144–174, 2017

    work page 2017

  23. [23]

    Hairer, C

    E. Hairer, C. Lubich, and G. Wanner. Geometric Numerical Integration: structure- preserving algorithms for ordinary differential equations, volume 31. Springer, 2006

    work page 2006

  24. [24]

    Han, A Brylev, X

    D. Han, A Brylev, X. Yang, and Z. Tan. Numerical analysis of second order, fully discrete energy stable schemes for phase field models of two phase incompressible flows. Journal of Scientific Computing, 70(3):965–989, 2017

    work page 2017

  25. [25]

    Hu and M

    W. Hu and M. lai. Unconditionally energy stable immersed boundary method with application to vesicle dyanmics. East Asian Journal on Applied Mathematics, 3(3):247–262, 2013

    work page 2013

  26. [26]

    L. Ju, X. Li, Z. Qiao, and H. Zhang. Energy stability and error estimates of ex- ponential time differencing schemes for the epitaxial growth model without slope selection. Mathematics of Computation, 87(312):1859–1885, 2018

    work page 2018

  27. [27]

    Phase-field models for multi-component fluid flows

    Junseok Kim. Phase-field models for multi-component fluid flows. Communication in Computational Physics, 12(3):613–661, 2012. 20

    work page 2012

  28. [28]

    R. V. Kohn and X. Yan. Upper bound on the coarsening rate for an epitaxial growth model. Comm. Pure Appl. Math., 56:1549–1564, 2003

    work page 2003

  29. [29]

    Li and J

    B. Li and J. Liu. Epitaxial growth without slope selection: energetics, coarsening and dynamic scaling. Journal of Nonlinear Science, 14:429–451, 2004

    work page 2004

  30. [30]

    Li and J

    B. Li and J. G. Liu. Thin film epitaxy with or without slope selection. European Journal of Applied Mathematics, 14:713–743, 2003

    work page 2003

  31. [31]

    Y. Li, J. Kim, and N. Wang. An unconditionally energy-stable second-order time- accurate scheme for the cahn-hilliard equation on surfaces. Commun Nonlinear Sci Numer Simulat, 53:213–227, 2017

    work page 2017

  32. [32]

    L. Onsager. Reciprocal relations in irreversible processes I. Physical Review, 37:405–426, 1931

    work page 1931

  33. [33]

    L. Onsager. Reciprocal relations in irreversible processes II. Physical Review, 38:2265–2279, 1931

    work page 1931

  34. [34]

    Z. Qiao, Z. Sun, and Z. Zhang. Stability and convergence of second-order schemes for the nonlinear epitaxial growth model without slope selection. Math. Comp., 84:653–674, 2015

    work page 2015

  35. [35]

    Z. Qiao, C. Wang, S. Wise, and Z. Zhang. Error analysis of a finite difference scheme for the epitaxial thin film model with slope selection with an improved convergence constant. International Journal of Numerical Analysis and Modeling, 14(2):1–23, 2017

    work page 2017

  36. [36]

    I. Romero. Thermodynamically consistent time-stepping algorithms for non-linear thermomechanical systems. International Journal for Numerical Methods in Engi- neering, 79:706–732, 2009

    work page 2009

  37. [37]

    J. Shen, C. Wang, X. Wang, and S. Wise. Second-order convex splitting schemes for gradient flows with Ehrlich-Schwoebel type energy: application to thin film epitaxy. SIAM Journal of Numerical Analysis, 50(1):105–125, 2012

    work page 2012

  38. [38]

    J. Shen, J. Xu, and J. Yang. A new class of efficient and robust energy stable schemes for gradient flows. ArXiv, page 1710.01331, 2017

    work page internal anchor Pith review Pith/arXiv arXiv 2017

  39. [39]

    J. Shen, J. Xu, and J. Yang. The scalar auxiliary variable (sav) approach for gra- dient flows. Journal of Computational Physics, 353:407–416, 2018

    work page 2018

  40. [40]

    Shen and X

    J. Shen and X. Yang. Numerical approximations of Allen-Cahn and Cahn-Hilliard equations. Disc. Conti. Dyn. Sys.-A, 28:1669–1691, 2010

    work page 2010

  41. [41]

    J. Shin, H. Lee, and J. Lee. Unconditionally stable methods for gradient flow using convex splitting runge-kutta scheme. Journal of Computational Physics, 347:367– 381, 2017

    work page 2017

  42. [42]

    C. Wang, X. Wang, and S. Wise. Unconditionally stable schemes for equations of thin film epitaxy. Discrete and Continuous Dynamic Systems, 28(1):405–423, 2010

    work page 2010

  43. [43]

    Wang and S

    C. Wang and S. Wise. An energy stable and convergent finite-difference scheme for the modified phase field crystal equation. SIAM Journal on Numerical Analysis, 49(3):945–969, 2011. 21

    work page 2011

  44. [44]

    S. Wang, R. Sekerka, A. Wheeler, B. T. Murray, S. R. Coriell, R. J. Braun, and G. B. McFadden. Thermodynamically-consistent phase-field models for solidifica- tion. Physica D, 69:189–200, 1993

    work page 1993

  45. [45]

    Xu and T

    C. Xu and T. Tang. Stability analysis of large time-stepping methods for epitaxial growth models. SIAM Journal of Numerical Analysis, 44(4):1759–1779, 2006

    work page 2006

  46. [46]

    X. Yang, J. Li, G. Forest, and Q. Wang. Hydrodynamic theories for flows of active liquid crystals and the generalized onsager principle. Entropy, 18(6):202, 2016

    work page 2016

  47. [47]

    Yang and J

    X. Yang and J. Zhao. On linear and unconditionally energy stable algorithms for variable mobility Cahn-Hilliard type equation with logarithmic Flory-Huggins po- tential. Communications in Computational Physics, 25:703–728, 2019

    work page 2019

  48. [48]

    X. Yang, J. Zhao, and X. He. Linear, second order and unconditionally energy stable schemes for the viscous cahn-hilliard equation with hyperbolic relaxation. Journal of Computational and Applied Mathematics, 343:80–97, 2018

    work page 2018

  49. [49]

    X. Yang, J. Zhao, and Q. Wang. Numerical approximations for the molecular beam epitaxial growth model based on the invariant energy quadratization method. Jour- nal of Computational Physics, 333:102–127, 2017

    work page 2017

  50. [50]

    X. Yang, J. Zhao, Q. Wang, and J. Shen. Numerical approximations for a three components cahn-hilliard phase-field model based on the invariant energy quadra- tization method. Mathematical Models and Methods in Applied Sciences, 27:1993– 2023, 2017

    work page 1993

  51. [51]

    J. Zhao, X. Yang, Y. Gong, X. Zhao, X. Yang, J. Li, and Q. Wang. A general strategy for numerical approximations of non-equilibrium models–part i thermo- dynamical systems. International Journal of Numerical Analysis and Modeling, 15(6):884–918, 2018. 22

    work page 2018