Structure-Preserving Augmented Lagrangian Finite Element Methods for Cavitation in Reynolds and Stokes Flows
Abstract
In this paper we propose augmented Lagrangian finite element methods for cavitation in the Reynolds and Stokes models of lubrication, for which the discrete complementarity conditions hold exactly and the augmentation parameter drops out of the method. For the Reynolds equation this follows from nodal quadrature in a mixed piecewise linear method, for which we prove the classical first-order error estimate. We also determine the computable stability threshold of a multiplier-free stabilised alternative. For Stokes flow we show that the constrained scalar is the mechanical pressure when the deviatoric stress, with zero bulk viscosity, is used, but not with the customary incompressible stress. We discretise with a jump-stabilised Crouzeix-Raviart element and piecewise constant pressure, and prove well-posedness, stability in two and three dimensions, and a first-order error estimate. We present numerical examples which verify these properties and show that close pressure profiles in the two models do not imply close cavity predictions, the latter also being sensitive to the end conditions of the computational domain.
Related
- Augmented Lagrangian methods with telescopic subproblems for MPCCs and cardinality constraints with general constraints: constraint-preserving swaps and B-stationarity
- A Structure-Preserving and Thermodynamically Compatible Cell-Centered Lagrangian Finite Volume Scheme for Continuum Mechanics
- A Structure-Preserving Semi-Lagrangian based bi-fidelity method for the Vlasov--Poisson system with uncertainties
- Memory-Augmented Phase-Field Modeling: Onsager Formulation, Internal-Variable Lifting, and Structure-Preserving Computation
- The Augmented Lagrangian Methods: Overview and Recent Advances
BibTeX
@misc{hansbo2026structurepreserving,
title = {{Structure-Preserving Augmented Lagrangian Finite Element Methods for Cavitation in Reynolds and Stokes Flows}},
author = {Peter Hansbo and Mats G. Larson},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.15223},
eprint = {2609.15223},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.15223},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Structure-Preserving Augmented Lagrangian Finite Element Methods for Cavitation in Reynolds and Stokes Flows
AU - Peter Hansbo
AU - Mats G. Larson
PY - 2026
DA - 2026/09/14
PB - arXiv
DO - 10.48550/arxiv.2609.15223
UR - https://arxiv.org/abs/2609.15223
AB - In this paper we propose augmented Lagrangian finite element methods for cavitation in the Reynolds and Stokes models of lubrication, for which the discrete complementarity conditions hold exactly and the augmentation parameter drops out of the method. For the Reynolds equation this follows from nodal quadrature in a mixed piecewise linear method, for which we prove the classical first-order error estimate. We also determine the computable stability threshold of a multiplier-free stabilised alternative. For Stokes flow we show that the constrained scalar is the mechanical pressure when the deviatoric stress, with zero bulk viscosity, is used, but not with the customary incompressible stress. We discretise with a jump-stabilised Crouzeix-Raviart element and piecewise constant pressure, and prove well-posedness, stability in two and three dimensions, and a first-order error estimate. We present numerical examples which verify these properties and show that close pressure profiles in the two models do not imply close cavity predictions, the latter also being sensitive to the end conditions of the computational domain.
ER -