Structure-preserving upwind Lagrange multiplier schemes for solid-state dewetting with a logarithmic Flory--Huggins potential
Abstract
Phase-field simulations of solid-state dewetting based on polynomial potentials may exhibit spurious bulk-diffusion coarsening, inconsistent with the surface-diffusion-dominated physics. To address this issue, we formulate a degenerate Cahn--Hilliard model with the logarithmic Flory--Huggins potential and dynamic contact line boundary conditions for the film--substrate--vapor triple junction. The logarithmic barrier confines the phase variable to its physical range and suppresses spurious coarsening at the continuum level. We develop a fully discrete, structure-preserving scheme that combines a Lagrange multiplier approach with an upwind finite-volume discretization and rigorously guarantees pointwise boundedness, mass conservation, and energy dissipation without artificial cut-offs or projections. A dimensional-splitting strategy further reduces computational cost while preserving these properties in each one-dimensional sweep. We also derive an explicit estimate of the equilibrium radius contraction caused by spontaneous film shrinking, showing that the logarithmic potential produces weaker spurious shrinkage than the polynomial potential. Numerical experiments confirm the analysis and demonstrate that, for small temperature parameters, the proposed scheme suppresses spurious coarsening and pinch-off while accurately reproducing surface-diffusion-dominated dynamics and relaxation toward equilibrium island structures.
Related
- A Second Order Positivity-Preserving Scheme for the Nonlocal Cahn-Hilliard System with Variable Mobility and Logarithmic Flory-Huggins Potential
- Uniqueness and structure-preserving numerical simulations for an Ultrahyperbolic Schrödinger equation
- A structure preserving finite volume scheme for the multidimensional model in cardiac electrophysiology
- A Structure-Preserving and Thermodynamically Compatible Cell-Centered Lagrangian Finite Volume Scheme for Continuum Mechanics
- A dual-field structure-preserving mixed finite element discretization for incompressible Hall MHD equations
BibTeX
@misc{huang2026structurepreserving,
title = {{Structure-preserving upwind Lagrange multiplier schemes for solid-state dewetting with a logarithmic Flory--Huggins potential}},
author = {Qiong-Ao Huang and Ying-wei Wang and Cheng Yuan},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.36719},
eprint = {2609.36719},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.36719},
note = {CC BY-NC-ND 4.0},
}
RIS
TY - UNPB
TI - Structure-preserving upwind Lagrange multiplier schemes for solid-state dewetting with a logarithmic Flory--Huggins potential
AU - Qiong-Ao Huang
AU - Ying-wei Wang
AU - Cheng Yuan
PY - 2026
DA - 2026/09/29
PB - arXiv
DO - 10.48550/arxiv.2609.36719
UR - https://arxiv.org/abs/2609.36719
AB - Phase-field simulations of solid-state dewetting based on polynomial potentials may exhibit spurious bulk-diffusion coarsening, inconsistent with the surface-diffusion-dominated physics. To address this issue, we formulate a degenerate Cahn--Hilliard model with the logarithmic Flory--Huggins potential and dynamic contact line boundary conditions for the film--substrate--vapor triple junction. The logarithmic barrier confines the phase variable to its physical range and suppresses spurious coarsening at the continuum level. We develop a fully discrete, structure-preserving scheme that combines a Lagrange multiplier approach with an upwind finite-volume discretization and rigorously guarantees pointwise boundedness, mass conservation, and energy dissipation without artificial cut-offs or projections. A dimensional-splitting strategy further reduces computational cost while preserving these properties in each one-dimensional sweep. We also derive an explicit estimate of the equilibrium radius contraction caused by spontaneous film shrinking, showing that the logarithmic potential produces weaker spurious shrinkage than the polynomial potential. Numerical experiments confirm the analysis and demonstrate that, for small temperature parameters, the proposed scheme suppresses spurious coarsening and pinch-off while accurately reproducing surface-diffusion-dominated dynamics and relaxation toward equilibrium island structures.
ER -