pith. sign in

arxiv: 2604.26402 · v1 · submitted 2026-04-29 · 🧮 math.NA · cs.NA

Energy stable auxiliary variable method for Cahn--Hilliard equations

Fei Xie , Nan Lu , Yajuan Sun This is my paper

Pith reviewed 2026-05-07 13:07 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords Cahn-Hilliard equationenergy stabilityauxiliary variable methodquadratic reformulationphase separationnumerical discretizationanisotropic evolution
0
0 comments X

The pith

Quadratic reformulation of rational-like free energies allows exact preservation of the energy dissipation law in Cahn-Hilliard discretizations.

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

The paper establishes a quadratic reformulation theory for rational-like functions and uses it to build the Quadratic Conservation Elevation method. This method pairs the scalar auxiliary variable approach with the implicit midpoint rule so that discretizations of the Cahn-Hilliard equation obey the same energy decrease law as the continuous model. A reader would care because many long-time phase separation simulations lose reliability when numerical energy can grow artificially. The scheme also reproduces the continuous dispersion relation and coarsening rates, and it handles anisotropic free energies to recover missing orientations.

Core claim

For the Cahn-Hilliard equation with rational-like free-energy terms, a quadratic reformulation theory yields the QCE discretization that exactly preserves the original energy dissipation law. The method combines the scalar auxiliary variable technique with the implicit midpoint rule, produces matching discrete dispersion relations and coarsening dynamics, and captures anisotropic evolution under various initial conditions.

What carries the argument

The Quadratic Conservation Elevation (QCE) method, formed by applying a quadratic reformulation to rational-like free energies and advancing the auxiliary variable with the implicit midpoint rule.

If this is right

  • The discrete solutions satisfy an energy inequality identical to the continuous dissipation law at every time step.
  • Dispersion relations and coarsening dynamics of the numerical scheme remain consistent with the continuous equation.
  • Anisotropic free energies can be treated without loss of stability or orientation capture.
  • The approach applies directly to a range of initial conditions that produce phase separation.

Where Pith is reading between the lines

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

  • The same reformulation idea may extend to other phase-field models whose energies are rational-like but not yet covered by the theory.
  • Guaranteed energy stability removes the need for artificial time-step restrictions in long material-science runs.
  • Testing the quadratic assumption on additional classes of rational functions would clarify the scope of exact preservation.

Load-bearing premise

The free-energy density must admit an exact quadratic reformulation whose auxiliary variable can be treated implicitly without breaking the energy identity.

What would settle it

A time-stepping experiment in which the discrete energy increases for one of the rational-like free energies considered, or in which the measured coarsening rate deviates from the continuous prediction, would falsify the preservation claim.

Figures

Figures reproduced from arXiv: 2604.26402 by Fei Xie, Nan Lu, Yajuan Sun.

Figure 1
Figure 1. Figure 1: Three commonly used potential functions, the red view at source ↗
Figure 2
Figure 2. Figure 2: Dispersion relation of the linearized isotropic Cahn–Hilliard equation for view at source ↗
Figure 3
Figure 3. Figure 3: Contours of g(k) in the (kx, ky)-plane for the isotropic case. Top left: exact; top right: explicit Euler; bottom left: implicit Euler; bottom right: implicit midpoint. As shown in view at source ↗
Figure 4
Figure 4. Figure 4: Phase separation and energy decay with different initial conditions. view at source ↗
Figure 5
Figure 5. Figure 5: Coarsening rate of the isotropic CH equation with view at source ↗
Figure 6
Figure 6. Figure 6: Dispersion relation of the linearized anisotropic CH equation with view at source ↗
Figure 7
Figure 7. Figure 7: Discrete dispersion relation of the linearized anisotropic CH equation with view at source ↗
Figure 8
Figure 8. Figure 8: Morphologies and normal-angle polar plots for the twofold anisotropic CH view at source ↗
Figure 9
Figure 9. Figure 9: Morphologies and normal-angle polar plots for the fourfold anisotropic CH view at source ↗
Figure 10
Figure 10. Figure 10: Two-droplet test: phase-field evolution and discrete energy decay. view at source ↗
Figure 11
Figure 11. Figure 11: Temporal convergence order of the proposed scheme. view at source ↗
Figure 12
Figure 12. Figure 12: Single-droplet test: phase-field evolution and discrete energy decay. view at source ↗
Figure 13
Figure 13. Figure 13: Phase-field evolution for the double-droplet test. Left: initial condition at view at source ↗
Figure 14
Figure 14. Figure 14: Evolution of discrete energy for the double-droplet test. view at source ↗
Figure 15
Figure 15. Figure 15: Phase-field evolution for the random test. Left: initial condition at view at source ↗
read the original abstract

In this paper, we propose a quadratic reformulation theory for rational-like functions. Based on this theory, we develop the Quadratic Conservation Elevation (QCE) method, which combines the Scalar Auxiliary Variable (SAV) method with the implicit midpoint rule. We apply this approach to the Cahn-Hilliard (CH) equation with rational-like free-energy terms, obtaining numerical discretizations that preserve the original energy dissipation law. We further derive the discrete dispersion relation and coarsening dynamics, confirming the efficiency and consistency of the method with the continuous counterpart. In addition, we use the proposed method to capture missing orientations for different anisotropic functions. Numerical simulations with various initial conditions illustrate phase separation and anisotropic evolution.

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

1 major / 2 minor

Summary. The manuscript proposes a quadratic reformulation theory for rational-like free-energy functions and develops the Quadratic Conservation Elevation (QCE) method by combining the Scalar Auxiliary Variable (SAV) approach with the implicit midpoint rule. This is applied to the Cahn-Hilliard equation to produce discretizations claimed to preserve the exact original energy dissipation law. The authors further derive the discrete dispersion relation and coarsening dynamics, demonstrate capture of missing orientations in anisotropic cases, and present numerical simulations of phase separation under various initial conditions.

Significance. If the quadratic reformulation is shown to be algebraically exact and the energy identity is preserved exactly under discretization, the work would extend energy-stable SAV-type schemes to a useful class of non-polynomial free energies in phase-field modeling. The additional dispersion and coarsening analysis, together with the anisotropic examples, would strengthen its applicability for long-time interfacial simulations.

major comments (1)
  1. [Quadratic reformulation theory and QCE discretization] The central claim that the QCE scheme preserves the original (not surrogate) energy dissipation law requires that the quadratic identity F(ϕ) = Q(r, ϕ) holds exactly for the rational-like densities considered and that this identity remains invariant under the implicit-midpoint averaging. The manuscript must supply the explicit algebraic construction of the auxiliary variable and quadratic form for each free-energy example, together with a direct verification that the midpoint discretization produces a telescoping sum identical to the continuous dissipation law without remainder terms.
minor comments (2)
  1. The abstract and introduction would benefit from a concise statement of the precise class of rational-like functions for which the quadratic reformulation is guaranteed to exist.
  2. Numerical tables or plots comparing the discrete energy decay against the continuous law for at least one non-trivial initial condition would make the preservation claim easier to inspect.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading and constructive comments. We address the major comment below and will incorporate the requested clarifications in a revised manuscript.

read point-by-point responses

  1. Referee: [Quadratic reformulation theory and QCE discretization] The central claim that the QCE scheme preserves the original (not surrogate) energy dissipation law requires that the quadratic identity F(ϕ) = Q(r, ϕ) holds exactly for the rational-like densities considered and that this identity remains invariant under the implicit-midpoint averaging. The manuscript must supply the explicit algebraic construction of the auxiliary variable and quadratic form for each free-energy example, together with a direct verification that the midpoint discretization produces a telescoping sum identical to the continuous dissipation law without remainder terms.

    Authors: We agree that explicit constructions strengthen the verification of exact energy preservation. Our quadratic reformulation theory defines the auxiliary variable r such that the identity F(ϕ) = Q(r, ϕ) holds algebraically and pointwise for the rational-like densities by direct substitution. Because the implicit midpoint rule is applied uniformly to the quadratic terms in both the evolution equation and the auxiliary update, the discrete inner product yields an exact telescoping sum identical to the continuous dissipation law, with no remainder. In the revised manuscript we will add, for each free-energy example, the explicit algebraic form of r and Q together with the step-by-step verification that the midpoint averaging produces the required telescoping identity without remainder terms. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation extends SAV framework via independent quadratic reformulation theory.

full rationale

The paper introduces a quadratic reformulation theory for rational-like functions as a new contribution, then combines it with the existing SAV approach and implicit midpoint discretization to obtain a scheme that preserves the original energy dissipation law for the Cahn-Hilliard equation. This preservation follows from applying the discrete scheme to the reformulated system rather than defining the target law in terms of the discrete solution itself. No load-bearing step reduces by construction to a fitted parameter, self-citation chain, or ansatz smuggled from prior work; the central claim rests on the proposed theory and its algebraic exactness for the functions considered, which is presented as verifiable independently of the numerical results. The derivation chain is therefore self-contained against external benchmarks such as the continuous energy law and standard SAV literature.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Insufficient information; abstract mentions a quadratic reformulation theory and rational-like functions but supplies no explicit free parameters, axioms, or new entities.

pith-pipeline@v0.9.0 · 5412 in / 1065 out tokens · 56508 ms · 2026-05-07T13:07:06.703374+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

3 extracted references · 3 canonical work pages

  1. [1]

    [1]W. J. Boettinger, J. A. W arren, C. Beckermann, and A. Karma,Phase-field simulation of solidification, Annual Review of Materials Research, 32 (2002), pp. 163–194, https: //doi.org/10.1146/annurev.matsci.32.101901.155803. FEI XIE, NAN LU, YAJUAN SUN25 [2]J. W. Cahn and J. E. Hilliard,Free energy of a nonuniform system. i. interfacial free energy, The J...

    work page doi:10.1146/annurev.matsci.32.101901.155803 2002

  2. [2]

    average vector field

    [3]J. W. Cahn and D. W. Hoffman,A vector thermodynamics for anisotropic surfaces ii. curved and faceted surfaces, Acta Metallurgica, 22 (1974), pp. 1205–1214, https://doi.org/10.101 6/0001-6160(74)90134-5. [4]E. Celledoni, V. Grimm, R. I. McLachlan, D. I. McLaren, D. O’Neale, B. Owren, and G. R. W. Quispel,Preserving energy resp. dissipation in numerical ...

    work page doi:10.1016/j.jcp.2012.06.022 1974

  3. [3]

    CRC Press, 2010.doi:10.1201/b10387

    [8]D. Furihata and T. Matsuo,Discrete Variational Derivative Method: A Structure-Preserving Numerical Method for Partial Differential Equations, Chapman and Hall/CRC, 2010, ht tps://doi.org/10.1201/b10387. [9]C. Herring,Surface tension as a motivation for sintering, in Fundamental contributions to the continuum theory of evolving phase interfaces in solid...

    work page doi:10.1201/b10387 2010