Related papers: The energy technique for the six-step BDF method
This paper studies multistep methods for the integration of reversible dynamical systems, with particular emphasis on the planar Kepler problem. It has previously been shown by Cano & Sanz-Serna that reversible linear multisteps for…
We study positive solutions to the steady state reaction diffusion systems of the form: \begin{equation} \left\{\begin{array}{ll} -\Delta u = \lambda f(v)+\mu h(u), & \Omega,\\ -\Delta v = \lambda g(u)+\mu q(v),& \Omega,\\ \frac{\partial…
We propose some finite element schemes to solve a class of fourth-order nonlinear PDEs, which include the vector-valued Landau--Lifshitz--Baryakhtar equation, the Swift--Hohenberg equation, and various Cahn--Hilliard-type equations with…
We present several first-order and second-order numerical schemes for the Cahn-Hilliard equation with discrete unconditional energy stability. These schemes stem from the generalized Positive Auxiliary Variable (gPAV) idea, and require only…
The two-phase Navier-Stokes Cahn-Hilliard (NSCH) mixture model is a key framework for simulating multiphase flows with non-matching densities. Developing fully discrete, energy-stable schemes for this model remains challenging, due to the…
Two primary scalar auxiliary variable (SAV) approaches are widely applied for simulating gradient flow systems, i.e., the nonlinear energy-based approach and the Lagrange multiplier approach. The former guarantees unconditional energy…
In this paper, we propose and analyze a linear, structure-preserving scalar auxiliary variable (SAV) method for solving the Allen--Cahn equation based on the second-order backward differentiation formula (BDF2) with variable time steps. To…
In earlier work [H. Liu and Z. Wang, J. Comput. Phys., 328(2017)], an arbitrary high-order conservative and energy-dissipative direct discontinuous Galerkin (DDG) scheme was developed. Although this scheme enforced solution positivity using…
{\bf Abstract} \,\, We consider the following nonlinear Schr\"{o}dinger equation on exterior domain. \begin{equation} \begin{cases} iu_t+\Delta_g u + ia(x)u - |u|^{p-1}u = 0 \qquad (x,t) \in \Omega\times (0,+\infty), \qquad (1)\cr…
A challenging difficulty in solving the radial Dirac eigenvalue problem numerically is the presence of spurious (unphysical) eigenvalues among the correct ones that are neither related to mathematical interpretations nor to physical…
The fully discrete problem for convection-diffusion equation is considered. It comprises compact approximations for spatial discretization, and Crank-Nicolson scheme for temporal discretization. The expressions for the entries of inverse of…
Let $L:= -a(x) (-\Delta)^{\alpha/2}+ (b(x), \nabla)$, where $\alpha\in (0,2)$, and $a:\rd\to (0,\infty)$, $b: \rd\to \rd$. Under certain regularity assumptions on the coefficients $a$ and $b$, we associate with the $C_\infty(\rd)$-closure…
We consider the nonlinear Gross-Pitaevskii equation at positive density, that is, for a bounded solution not tending to 0 at infinity. We focus on infinite ground states, which are by definition minimizers of the energy under local…
We obtain eigenvalues and eigenfunctions of the Schr\"{o}dinger equation with a hyperbolic double-well potential. We consider exact polynomial solutions for some particular values of the potential-strength parameter and also numerical…
A general method using multipliers for finding the conserved integrals for any system of partial differential equations (PDEs) is reviewed and further developed in several ways. Multipliers are expressions whose (summed) product with a PDE…
We present a diagrammatic approach to construct self-energy approximations within many-body perturbation theory with positive spectral properties. The method cures the problem of negative spectral functions which arises from a…
The equation arising from Prandtl boundary layer theory is considered. The existence of the entropy solution can be proved by BV estimate method. The interesting problem is that, since a may be degenerate on the boundary, the usual boundary…
A new method to remove the stiffness of partial differential equations is presented. Two terms are added to the right-hand-side of the PDE : the first is a damping term and is treated implicitly, the second is of the opposite sign and is…
We introduce and analyze a discontinuous Petrov-Galerkin method with optimal test functions for the heat equation. The scheme is based on the backward Euler time stepping and uses an ultra-weak variational formulation at each time step. We…
A high-order, well-balanced, positivity-preserving quasi-Lagrange moving mesh DG method is presented for the shallow water equations with non-flat bottom topography. The well-balance property is crucial to the ability of a scheme to…