Related papers: A new Lagrange multiplier approach for constructin…
This paper deals with time stepping schemes for the Cahn--Hilliard equation with three different types of dynamic boundary conditions. The proposed schemes of first and second order are mass-conservative and energy-dissipative and -- as…
Simulating infiltration in porous media using Richards' equation remains computationally challenging due to its parabolic structure and nonlinear coefficients. While a wide range of numerical methods for differential equations have been…
Multidimensional consistency has emerged as a key integrability property for partial difference equations (P$\Delta$Es) defined on the "space-time" lattice. It has led, among other major insights, to a classification of scalar affine-linear…
This work focuses on the development of a non-conforming domain decomposition method for the approximation of PDEs based on weakly imposed transmission conditions: the continuity of the global solution is enforced by a discrete number of…
Modeling the chemical, electric, and thermal transport as well as phase transitions and the accompanying mesoscale microstructure evolution within a material in an electronic device setting involves the solution of partial differential…
We develop a general framework for numerically solving differential equations while preserving invariants. As in standard projection methods, we project an arbitrary base integrator onto an invariant-preserving manifold, however, our method…
For decades, the computational multiphase flow community has grappled with mass loss in the level set method. Numerous solutions have been proposed, from fixing the reinitialization step to combining the level set method with other…
In this work, we focus on the numerical approximation of the shallow water equations in two space dimensions. Our aim is to propose a well-balanced, all-regime and positive scheme. By well-balanced, it is meant that the scheme is able to…
We consider the semiclassical limit for the nonlinear Schrodinger equation. We introduce a phase/amplitude representation given by a system similar to the hydrodynamical formulation, whose novelty consists in including some asymptotically…
We propose and analyse a boundary-preserving numerical scheme for the weak approximation for some stochastic partial differential equations (SPDEs) with bounded state-space. We impose regularity assumptions on the drift and diffusion…
This paper aims to construct structure-preserving numerical schemes for multi-dimensional space fractional Klein-Gordon-Schr\"{o}dinger equation, which are based on the newly developed partitioned averaged vector field methods. First, we…
In this article we propose a novel strategy for choosing the Lagrange multipliers in the Levenberg-Marquardt method for solving ill-posed problems modeled by nonlinear operators acting between Hilbert spaces. Convergence analysis results…
We develop several efficient numerical schemes which preserve exactly the global constraints for constrained gradient flows. Our schemes are based on the SAV approach combined with the Lagrangian multiplier approach. They are as efficient…
We propose a boundary value correction approach for cases when curved boundaries are approximated by straight lines (planes) and Lagrange multipliers are used to enforce Dirichlet boundary conditions. The approach allows for optimal order…
It is well known that the classic Allen-Cahn equation satisfies the maximum bound principle (MBP), that is, the absolute value of its solution is uniformly bounded for all time by certain constant under suitable initial and boundary…
A high-order numerical method is developed for solving the Cahn-Hilliard-Navier-Stokes equations with the Flory-Huggins potential. The scheme is based on the $Q_k$ finite element with mass lumping on rectangular grids, the second-order…
The proximal Galerkin finite element method is a high-order, low-iteration complexity, nonlinear numerical method that preserves the geometric and algebraic structure of point-wise bound constraints in infinite-dimensional function spaces.…
An iterative optimization approach that simultaneously minimizes the energy and optimizes the Lagrange multipliers enforcing desired constraints is presented. The method is tested on previously established benchmark systems and it is proved…
We propose a model for the coupling between free fluid and a linearized poro-hyperelastic body. In this model, the Brinkman equation is employed for fluid flow in the porous medium, incorporating inertial effects into the fluid dynamics. A…
In this work, we present a novel family of high order accurate numerical schemes for the solution of hyperbolic partial differential equations (PDEs) which combines several geometrical and physical structure preserving properties. First, we…