Divergence-free and mass-conservative virtual element methods for the Navier-Stokes-Cahn-Hilliard system

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Silgado, Alberth, Vacca, Giuseppe
Format: Preprint
Published: 2026
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866915756448415744
author Silgado, Alberth
Vacca, Giuseppe
author_facet Silgado, Alberth
Vacca, Giuseppe
contents In this work, we design and analyze semi/fully-discrete virtual element approximations for the time-dependent Navier--Stokes-Cahn--Hilliard equations, modeling the dynamics of two-phase incompressible fluid flows with diffuse interfaces. A new variational formulation is derived involving solely the velocity, pressure, and phase field, together with corresponding a priori energy estimates. The spatial discretization is based on the coupling divergence-free and $C^1$-conforming elements of high-order, while the time discretization employs a classical backward Euler scheme. By introducing a novel skew-symmetric trilinear form to discretize the convective term in the Cahn--Hilliard equation, we propose discrete schemes that satisfy mass conservation and energy bounds. Moreover, optimal error estimates are provided for both formulations. Finally, two numerical experiments are presented to support our theoretical findings and to illustrate the good performance of the proposed schemes for different polynomial degrees and polygonal meshes.
format Preprint
id arxiv_https___arxiv_org_abs_2601_18758
institution arXiv
publishDate 2026
record_format arxiv
spellingShingle Divergence-free and mass-conservative virtual element methods for the Navier-Stokes-Cahn-Hilliard system
Silgado, Alberth
Vacca, Giuseppe
Numerical Analysis
65M60, 35K55, 76D05
In this work, we design and analyze semi/fully-discrete virtual element approximations for the time-dependent Navier--Stokes-Cahn--Hilliard equations, modeling the dynamics of two-phase incompressible fluid flows with diffuse interfaces. A new variational formulation is derived involving solely the velocity, pressure, and phase field, together with corresponding a priori energy estimates. The spatial discretization is based on the coupling divergence-free and $C^1$-conforming elements of high-order, while the time discretization employs a classical backward Euler scheme. By introducing a novel skew-symmetric trilinear form to discretize the convective term in the Cahn--Hilliard equation, we propose discrete schemes that satisfy mass conservation and energy bounds. Moreover, optimal error estimates are provided for both formulations. Finally, two numerical experiments are presented to support our theoretical findings and to illustrate the good performance of the proposed schemes for different polynomial degrees and polygonal meshes.
title Divergence-free and mass-conservative virtual element methods for the Navier-Stokes-Cahn-Hilliard system
topic Numerical Analysis
65M60, 35K55, 76D05
url https://arxiv.org/abs/2601.18758