pith. sign in

arxiv: 1907.04142 · v2 · pith:XTU5K6TVnew · submitted 2019-07-01 · 🧮 math.NA · cs.NA· math.AP

Energy stable schemes for gradient flows based on novel auxiliary variable with energy bounded above

Zhengguang Liu This is my paper

Pith reviewed 2026-05-25 12:11 UTC · model grok-4.3

classification 🧮 math.NA cs.NAmath.AP
keywords gradient flowsenergy stable schemesauxiliary variable methodNAEVSAVIEQunconditional stabilitynumerical simulation
0
0 comments X

The pith

A novel auxiliary energy variable method for gradient flows guarantees the square root term stays positive by assuming the energy is bounded above.

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

The paper develops energy stable schemes for gradient flows using a new auxiliary variable approach. It defines the variable based on an upper bound of the energy functional instead of the usual lower bound. This change guarantees that expressions under the square root remain positive without additional restrictions. The author rigorously proves unconditional energy stability for the resulting semi-discrete schemes. Numerical comparisons with the SAV method and 2D simulations confirm the approach's accuracy and efficiency.

Core claim

The NAEV method introduces an auxiliary variable constructed from the assumption that the energy functional is bounded from above. This construction ensures the computed values for the functional inside the square root are guaranteed to be positive. Unlike IEQ and SAV approaches, the method does not require any bounded-below restrictions on the energy. The paper proves unconditional energy stability for all semi-discrete schemes and demonstrates accuracy through comparisons and simulations.

What carries the argument

The novel auxiliary energy variable (NAEV) defined using an upper bound on the energy functional to keep the radicand positive by construction.

If this is right

  • Unconditional energy stability holds for all the semi-discrete schemes.
  • The values inside the square root are guaranteed positive without extra conditions.
  • The bounded-below assumption required by SAV and IEQ is no longer needed.
  • Comparative tests show accuracy and efficiency similar to the classical SAV approach.
  • Two-dimensional numerical simulations confirm stability and accuracy for the schemes.

Where Pith is reading between the lines

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

  • This formulation could apply directly to gradient flows whose energy decreases without a lower limit but respects a known upper limit.
  • Choosing the upper bound constant may influence conditioning or accuracy in long simulations.
  • The positivity guarantee could extend the method to other dissipative systems where standard quadratization encounters sign issues.
  • Similar upper-bound reformulations might be tested on higher-order time discretizations while preserving the stability proof structure.

Load-bearing premise

The energy functional of the gradient flow problem is bounded from above.

What would settle it

Apply the NAEV scheme to a gradient flow whose energy functional has no upper bound and check whether the square root expression becomes negative or the discrete energy fails to decrease.

Figures

Figures reproduced from arXiv: 1907.04142 by Zhengguang Liu.

Figure 3
Figure 3. Figure 3 [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 3
Figure 3. Figure 3 [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 3
Figure 3. Figure 3: shows the evolution of [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figure 3
Figure 3. Figure 3 [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figure 3
Figure 3. Figure 3 [PITH_FULL_IMAGE:figures/full_fig_p012_3.png] view at source ↗
read the original abstract

In this paper, we consider a novel auxiliary variable method to obtain energy stable schemes for gradient flows. The auxiliary variable based on energy bounded above does not limited to the hypothetical conditions adopted in previous approaches. We proved the unconditional energy stability for all the semi-discrete schemes carefully and rigorously. The novelty of the proposed schemes is that the computed values for the functional in square root are guaranteed to be positive. This method, termed novel auxiliary energy variable (NAEV) method does not consider any bounded below restrictions any longer. However, these restrictions are necessary in invariant energy quadratization (IEQ) and scalar auxiliary variable (SAV) approaches which are very popular methods recently. This property of guaranteed positivity is not available in previous approaches. A comparative study of classical SAV and NAEV approaches is considered to show the accuracy and efficiency. Finally, we present various 2D numerical simulations to demonstrate the stability and accuracy.

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 manuscript proposes the novel auxiliary energy variable (NAEV) method for energy-stable schemes for gradient flows. It defines an auxiliary variable using an upper bound on the energy (so that the square-root term remains positive) rather than requiring a lower bound as in SAV and IEQ. The authors state that they prove unconditional energy stability for the semi-discrete schemes, present a comparative study with classical SAV, and show 2D numerical simulations demonstrating stability and accuracy. The key claimed novelty is the guaranteed positivity of the auxiliary variable without bounded-below restrictions.

Significance. If the stability proofs are rigorous as claimed, NAEV removes a structural limitation of SAV/IEQ and guarantees positivity of the auxiliary quantity by construction. This could extend auxiliary-variable techniques to a wider class of dissipative systems. The comparative study and simulations provide concrete evidence of practical performance.

minor comments (3)
  1. [Abstract] Abstract: the claim of 'rigorous' unconditional stability proofs is stated without any derivation outline or key estimate; the manuscript should include at least a one-paragraph sketch of the stability argument in the abstract or introduction.
  2. The comparative study section should report quantitative metrics (e.g., L2 errors, CPU times, or convergence rates) rather than qualitative statements to support the accuracy and efficiency claims.
  3. Notation for the upper bound C and the auxiliary variable r should be introduced with an explicit equation number at first use to improve readability.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive summary, significance assessment, and recommendation of minor revision. No major comments were provided in the report, so we have no point-by-point responses to address. We are prepared to make any minor changes requested by the editor or in a subsequent round.

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained

full rationale

The paper introduces the NAEV auxiliary variable by explicitly assuming an a priori upper bound on the energy functional (to ensure the square-root term remains positive) and then proves unconditional energy stability for the resulting semi-discrete schemes. This construction is presented as a direct reformulation that removes the lower-bound hypothesis required by SAV/IEQ; the stability proof is stated to be rigorous and independent of fitted parameters or prior self-citations. No load-bearing step reduces by definition or by self-referential fitting to the input data or to an unverified uniqueness claim. The central guarantee (positivity of the auxiliary expression) follows immediately from the stated upper-bound assumption rather than from any circular redefinition of the energy itself.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 1 invented entities

The central claim rests on redefining the auxiliary variable via an upper bound on the energy functional; no free parameters, additional axioms, or invented physical entities are mentioned in the abstract.

axioms (1)
  • domain assumption The energy functional is bounded from above
    Invoked to construct the auxiliary variable that guarantees positivity of the square-root expression.
invented entities (1)
  • Novel auxiliary energy variable (NAEV) no independent evidence
    purpose: To replace the SAV/IEQ auxiliary variable and remove the bounded-below restriction while preserving unconditional stability
    Newly introduced construction whose only stated support is the positivity guarantee and stability proof claimed in the abstract.

pith-pipeline@v0.9.0 · 5686 in / 1249 out tokens · 39849 ms · 2026-05-25T12:11:20.808150+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

22 extracted references · 22 canonical work pages · 3 internal anchors

  1. [1]

    Ainsworth and Z

    M. Ainsworth and Z. Mao , Analysis and approximation of a fractional cahn–hilliard equation , SIAM Journal on Nu- merical Analysis, 55 (2017), pp. 1689–1718

    work page 2017

  2. [2]

    Ambati, T

    M. Ambati, T. Gerasimov, and L. De Lorenzis , A review on phase-field models of brittle fracture and a new fast hybrid formulation, Computational Mechanics, 55 (2015), pp. 383–405

    work page 2015

  3. [3]

    L. Chen, J. Zhao, W. Cao, H. W ang, and J. Zhang , An accurate and efficient algorithm for the time-fractional molecular beam epitaxy model with slope selection , arXiv preprint arXiv:1803.01963, (2018)

    work page arXiv 2018

  4. [4]

    L. Chen, J. Zhao, and H. W ang , On power law scaling dynamics for time-fractional phase field models during coarsening, arXiv preprint arXiv:1803.05128, (2018)

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  5. [5]

    Q. Du, L. Ju, X. Li, and Z. Qiao , Stabilized linear semi-implicit schemes for the nonlocal cahn–hilliard equation, Journal of Computational Physics, 363 (2018), pp. 39–54

    work page 2018

  6. [6]

    D. J. Eyre , Unconditionally gradient stable time marching the cahn-hilliard equation , MRS Online Proceedings Library Archive, 529 (1998)

    work page 1998

  7. [7]

    Guo and P

    Z. Guo and P. Lin , A thermodynamically consistent phase-field model for two-phase flows with thermocapillary effects , Journal of Fluid Mechanics, 766 (2015), pp. 226–271

    work page 2015

  8. [8]

    H. G. Lee , An energy stable method for the swift–hohenberg equation with quadratic–cubic nonlinearity , Computer Meth- ods in Applied Mechanics and Engineering, 343 (2019), pp. 40–51

    work page 2019

  9. [9]

    L. Lin, Z. Yang, and S. Dong , Numerical approximation of incompressible navier-stokes equations based on an auxiliary energy variable, Journal of Computational Physics, 388 (2019), pp. 1–22

    work page 2019

  10. [10]

    Liu and X

    Z. Liu and X. Li , Efficient modified techniques of invariant energy quadratization approach for gradient flows , Applied Mathematics Letters, 98 (2019), pp. 206–214

    work page 2019

  11. [11]

    Marth, S

    W. Marth, S. Aland, and A. Voigt , Margination of white blood cells: a computational approach by a hydrodynamic phase field model , Journal of Fluid Mechanics, 790 (2016), pp. 389–406

    work page 2016

  12. [12]

    Miehe, M

    C. Miehe, M. Hofacker, and F. Welschinger , A phase field model for rate-independent crack propagation: Robust algorithmic implementation based on operator splits , Computer Methods in Applied Mechanics and Engineering, 199 (2010), pp. 2765–2778

    work page 2010

  13. [13]

    J. Shen, J. Xu, and J. Yang , A new class of efficient and robust energy stable schemes for gradient flows , arXiv preprint arXiv:1710.01331, (2017)

    work page internal anchor Pith review Pith/arXiv arXiv 2017

  14. [14]

    J. Shen, J. Xu, and J. Yang , The scalar auxiliary variable (SAV) approach for gradient flows , Journal of Computational Physics, 353 (2018), pp. 407–416

    work page 2018

  15. [15]

    Shen and X

    J. Shen and X. Yang , Numerical approximations of Allen-Cahn and Cahn-Hilliard equations , Discrete Contin. Dyn. Syst, 28 (2010), pp. 1669–1691

    work page 2010

  16. [16]

    J. Shen, X. Yang, and H. Yu , Efficient energy stable numerical schemes for a phase field moving contact line model , Journal of Computational Physics, 284 (2015), pp. 617–630

    work page 2015

  17. [17]

    J. Shin, H. G. Lee, and J.-Y. Lee , First and second order numerical methods based on a new convex splitting for phase-field crystal equation , Journal of Computational Physics, 327 (2016), pp. 519–542

    work page 2016

  18. [18]

    A. A. Wheeler, W. J. Boettinger, and G. B. McF adden , Phase-field model for isothermal phase transitions in binary alloys, Physical Review A, 45 (1992), p. 7424

    work page 1992

  19. [19]

    A. A. Wheeler, B. T. Murray, and R. J. Schaefer , Computation of dendrites using a phase field model , Physica D: 12 ZHENGGUANG LIU (a) t=2 (b) t=10 (c) t=30 (d) t=40 (e) t=60 (f) t=120 Fig. 3.5 . Snapshots of the phase variable φ are taken at t=2, 10, 30, 40, 60, 120 for example 6. Nonlinear Phenomena, 66 (1993), pp. 243–262

    work page 1993

  20. [20]

    X. Yang , Linear, first and second-order, unconditionally energy stable numerical schemes for the phase field model of homopolymer blends, Journal of Computational Physics, 327 (2016), pp. 294–316

    work page 2016

  21. [21]

    Yang and D

    X. Yang and D. Han , Linearly first-and second-order, unconditionally energy stable schemes for the phase field crystal model, Journal of Computational Physics, 330 (2017), pp. 1116–1134

    work page 2017

  22. [22]

    Numerical approximations of the Cahn-Hilliard and Allen-Cahn Equations with general nonlinear potential using the Invariant Energy Quadratization approach

    X. Yang and G. Zhang , Numerical approximations of the Cahn-Hilliard and Allen-Cahn equations with general nonlinear potential using the Invariant Energy Quadratization approach , arXiv preprint arXiv:1712.02760, (2017)

    work page internal anchor Pith review Pith/arXiv arXiv 2017