Related papers: A second order accurate scalar auxiliary variable …
This paper develops a two-level fourth-order scheme for solving time-fractional convection-diffusion-reaction equation with variable coefficients subjected to suitable initial and boundary conditions. The basis properties of the new…
We introduce novel entropy-dissipative numerical schemes for a class of kinetic equations, leveraging the recently introduced scalar auxiliary variable (SAV) approach. Both first and second order schemes are constructed. Since the…
It is well-known that the Allen-Cahn equation not only satisfies the energy dissipation law but also possesses the maximum bound principle (MBP) in the sense that the absolute value of its solution is pointwise bounded for all time by some…
This paper proposes a new class of arbitrarily high-order conservative numerical schemes for the generalized Korteweg-de Vries (KdV) equation. This approach is based on the scalar auxiliary variable (SAV) method. The equation is…
The study of polycrystalline materials requires theoretical and computational techniques enabling multiscale investigations. The amplitude expansion of the phase field crystal model (APFC) allows for describing crystal lattice properties on…
The classic second-order average vector field (AVF) method can exactly preserve the energy for Hamiltonian ordinary differential equations and partial differential equations. However, the AVF method inevitably leads to fully-implicit…
The Smectic-A (SmA) phase is modeled by a modified Landau-de Gennes (mLdG) model proposed by Xia et al. [Phys. Rev. Lett., 126 (2021), 177801], in which a tensor order parameter $\mathbf{Q}$ for the orientational order is coupled with a…
The energy dissipation law and the maximum bound principle (MBP) are two important physical features of the well-known Allen-Cahn equation. While some commonly-used first-order time stepping schemes have turned out to preserve…
We consider solid state phase transformations that are caused by free energy densities with domains of non-convexity in strain-composition space; we refer to the non-convex domains as mechano-chemical spinodals. The non-convexity with…
In this paper, we present energy-stable numerical schemes for a Smectic-A liquid crystal model. This model involve the hydrodynamic velocity-pressure macroscopic variables $({\bf u},p)$ and the microscopic order parameter of Smectic-A…
We present a numerical scheme for approximating the incompressible Navier-Stokes equations based on an auxiliary variable associated with the total system energy. By introducing a dynamic equation for the auxiliary variable and…
In this paper, we propose and analyze a first-order and a second-order time-stepping schemes for the anisotropic phase-field dendritic crystal growth model. The proposed schemes are based on an auxiliary variable approach for the Allen-Cahn…
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 consider the numerical approximations for a phase field model consisting of incompressible Navier--Stokes equations with a generalized Navier boundary condition, and the Cahn-Hilliard equation with a dynamic moving contact line boundary…
This report extends the mathematical support of a subgrid artificial viscosity (SAV) method to simulate the incompressible Navier-Stokes equations to better performing a linearly extrapolated BDF2 (BDF2LE) time discretization. The method…
Crystalline defects critically influence material properties, necessitating accurate simulation methods. Existing approaches, from atomic-scale configurations to continuum elasticity, face inherent limitations in modeling…
How to develop efficient numerical schemes while preserving the energy stability at the discrete level is a challenging issue for the three component Cahn-Hilliard phase-field model. In this paper, we develop first and second order temporal…
We propose a novel second order in time numerical scheme for Cahn-Hilliard-Navier- Stokes phase field model with matched density. The scheme is based on second order convex-splitting for the Cahn-Hilliard equation and pressure-projection…
Self-consistent field theory (SCFT) is one of the most widely-used framework in studying the equilibrium phase behaviors of inhomogenous polymers. For liquid crystalline polymeric systems, the main numerical challenges of solving SCFT…
In this paper, we propose and analyze an efficient numerical method for the anisotropic phase field dendritic crystal growth model, which is challenging because we are facing the nonlinear coupling and anisotropic coefficient in the model.…