Related papers: A new class of efficient and robust energy stable …
We present an energy/entropy stable and high order accurate finite difference (FD) method for solving the nonlinear (rotating) shallow water equations (SWEs) in vector invariant form using the newly developed dual-pairing and…
We propose an unconditionally energy-stable, orthonormality-preserving, component-wise splitting iterative scheme for the Kohn-Sham gradient flow based model in the electronic structure calculation. We first study the scheme discretized in…
Modelling interfacial dynamics with soluble surfactants in a multiphase system is a challenging task. Here, we consider the numerical approximation of a phase-field surfactant model with fluid flow. The nonlinearly coupled model consists of…
We develop a stable and efficient numerical scheme for modeling the optical field evolution in a nonlinear dispersive cavity with counter propagating waves and complex, semiconductor physics gain dynamics that are expensive to evaluate. Our…
Recently, the error analysis of BDF$k$ $(1\leqslant k\leqslant5)$ SAV (scalar auxiliary variable) schemes are given in \cite{Huangg:20} for the classical Allen-Cahn equation. However, it remains unavailable for BDF$6$ SAV schemes. In this…
A stable added-mass partitioned (AMP) algorithm is developed for fluid-structure interaction (FSI) problems involving viscous incompressible flow and compressible elastic solids. Deforming composite grids are used to effectively handle the…
In this paper we combine a flexible covariant formulation of the shallow water equations with the semi-implicit numerical scheme developed over the years by Casulli and collaborators. After adopting an orthogonal, but non-orthonormal,…
We carry out a rigorous error analysis of the first-order semi-discrete (in time) consistent splitting scheme coupled with a generalized scalar auxiliary variable (GSAV) approach for the Navier-Stokes equations with no-slip boundary…
This work uses a linear relaxation method to develop efficient numerical schemes for the time-fractional Allen-Cahn and Cahn-Hilliard equations. The L1+-CN formula is used to discretize the fractional derivative, and an auxiliary variable…
We present a novel optimization algorithm, element-wise relaxed scalar auxiliary variable (E-RSAV), that satisfies an unconditional energy dissipation law and exhibits improved alignment between the modified and the original energy. Our…
In this study, a fast and stable machine-learned hybrid algorithm implemented in TensorFlow for the integration of stiff chemical kinetics is introduced. Numerical solutions to differential equations are at the core of computational fluid…
In this paper, we consider numerical approximation of an electrically conductive ferrofluid model, which consists of Navier-Stokes equations, magnetization equation, and magnetic induction equation. To solve this highly coupled, nonlinear,…
We present a second-order strictly length-preserving and unconditionally energy-stable rotational discrete gradient (Rdg) scheme for the numerical approximation of the Oseen-Frank gradient flows with anisotropic elastic energy functional.…
We introduce a class of generic spike-and-slab priors for high-dimensional linear regression with grouped variables and present a Coordinate-ascent Variational Inference (CAVI) algorithm for obtaining an optimal variational Bayes…
Comparing with the classical local gradient flow and phase field models, the nonlocal models such as nonlocal Cahn-Hilliard equations equipped with nonlocal diffusion operator can describe more practical phenomena for modeling phase…
Implicit integration of the viscous term can significantly improve performance in computational fluid dynamics for highly viscous fluids such as lava. We show improvements over our previous proposal for semi-implicit viscous integration in…
The Cahn-Hilliard equation is one of the most common models to describe phase separation processes in mixtures of two materials. For a better description of short-range interactions between the material and the boundary, various dynamic…
Energy stable flux reconstruction (ESFR) is a high-order numerical method used for solving partial differential equations in computational fluid dynamics. This method is designed to preserve the energy stability of the underlying partial…
A stable partitioned algorithm for fluid-structure interaction (FSI) problems that couple viscous incompressible flow with structural shells or beams is described. This added-mass partitioned (AMP) scheme uses Robin (mixed) interface…
We present an efficient and flexible method for solving the non-linear lasing equations of the steady-state ab initio laser theory. Our strategy is to solve the underlying system of partial differential equations directly, without the need…