A linear, second-order, unconditionally energy-stable scalar-auxiliary-variable scheme for the CHQ viscoelastic phase-separation model
Abstract
The CHQ model describes viscoelastic phase separation in a polymer--solvent mixture: a Cahn--Hilliard equation for the polymer volume fraction ϕ is two-way coupled, through a phase-dependent cross-diffusion flux, to a damped relaxation equation for a scalar bulk elastic stress q. We propose a scalar-auxiliary-variable (SAV) time discretisation of this system, combining the SAV reformulation of the double-well potential with a second-order backward differentiation formula (BDF2) and a second-order extrapolation of all phase-dependent mobilities: a linear, computationally cheaper complement to the reference scheme's fully nonlinear Newton solve. The resulting method is (i) linear at every time step, with a coupling structure we prove has exactly one solution at every step, for every choice of time-step size; (ii) mass conserving, exactly, as an algebraic identity of the discrete equations themselves (the only drift ever seen numerically traces to the finite tolerance of the iterative linear solver, not a failure of the conservation property); and (iii) unconditionally stable for a modified discrete energy, via an exact one-step energy--dissipation identity. We show the energy identity extends to this filtered variant as well, with two explicit extra non-negative terms. Hence, the scheme actually producing the physics results is itself provably stable for the same energy. A Fourier--Galerkin implementation verifies that the energy identity holds at every time step to within floating-point round-off. It also confirms second-order temporal accuracy, including in the steep-interface regime, and agreement with an independently implemented nonlinear Newton--BDF2 solver. Compared with that independent nonlinear Newton--BDF2 solver, at equal grid resolution, the SAV method uses $2.5−−3.0\times$ fewer linear-solver iterations and $1.7−−2.0\times$ less wall-clock time.
Related
- Convergence Analysis of Variable-Separation-Based Model Reduction for Linear Parametric Partial Differential Equations
- A Generalized Scalar Auxiliary Variable Method for Structure-Preserving and Efficient Integration of Nonlinear port-Hamiltonian DAEs
- Dissipative solutions for viscoelastic phase separation with degenerate mobility: a dynamic transport approach
- Blow-up at arbitrary energy levels for a strongly damped viscoelastic wave equation with a variable-exponent logarithmic source
- A pullback-corrected scalar auxiliary variable optimizer with momentum and adaptive mobility
BibTeX
@misc{bhola2026linear,
title = {{A linear, second-order, unconditionally energy-stable scalar-auxiliary-variable scheme for the CHQ viscoelastic phase-separation model}},
author = {Ritesh Bhola and Michael te Vrugt},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.30077},
eprint = {2610.30077},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.30077},
}
RIS
TY - UNPB
TI - A linear, second-order, unconditionally energy-stable scalar-auxiliary-variable scheme for the CHQ viscoelastic phase-separation model
AU - Ritesh Bhola
AU - Michael te Vrugt
PY - 2026
DA - 2026/10/21
PB - arXiv
DO - 10.48550/arxiv.2610.30077
UR - https://arxiv.org/abs/2610.30077
AB - The CHQ model describes viscoelastic phase separation in a polymer--solvent mixture: a Cahn--Hilliard equation for the polymer volume fraction ϕ is two-way coupled, through a phase-dependent cross-diffusion flux, to a damped relaxation equation for a scalar bulk elastic stress q. We propose a scalar-auxiliary-variable (SAV) time discretisation of this system, combining the SAV reformulation of the double-well potential with a second-order backward differentiation formula (BDF2) and a second-order extrapolation of all phase-dependent mobilities: a linear, computationally cheaper complement to the reference scheme's fully nonlinear Newton solve. The resulting method is (i) linear at every time step, with a coupling structure we prove has exactly one solution at every step, for every choice of time-step size; (ii) mass conserving, exactly, as an algebraic identity of the discrete equations themselves (the only drift ever seen numerically traces to the finite tolerance of the iterative linear solver, not a failure of the conservation property); and (iii) unconditionally stable for a modified discrete energy, via an exact one-step energy--dissipation identity. We show the energy identity extends to this filtered variant as well, with two explicit extra non-negative terms. Hence, the scheme actually producing the physics results is itself provably stable for the same energy. A Fourier--Galerkin implementation verifies that the energy identity holds at every time step to within floating-point round-off. It also confirms second-order temporal accuracy, including in the steep-interface regime, and agreement with an independently implemented nonlinear Newton--BDF2 solver. Compared with that independent nonlinear Newton--BDF2 solver, at equal grid resolution, the SAV method uses $2.5−−3.0\times$ fewer linear-solver iterations and $1.7−−2.0\times$ less wall-clock time.
ER -