Related papers: Monotonicity-preserving finite element schemes bas…
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…
This paper develops an efficient and robust solution technique for the steady Boussinesq model of non-isothermal flow using Anderson acceleration applied to a Picard iteration. After analyzing the fixed point operator associated with the…
This paper addresses the analysis and numerical assessment of a computational method for solving the Cahn--Hilliard equation defined on a surface. The proposed approach combines the stabilized trace finite element method for spatial…
We present an energy-stable scheme for numerically approximating the governing equations for incompressible two-phase flows with different densities and dynamic viscosities for the two fluids. The proposed scheme employs a scalar-valued…
A second order accurate, linear numerical method is analyzed for the Landau-Lifshitz equation with large damping parameters. This equation describes the dynamics of magnetization, with a non-convexity constraint of unit length of the…
We introduce an explicit adaptive Milstein method for stochastic differential equations (SDEs) with no commutativity condition. The drift and diffusion are separately locally Lipschitz and together satisfy a monotone condition. This method…
This work extends the concepts of algebraic flux correction and convex limiting to continuous high-order Bernstein finite element discretizations of scalar hyperbolic problems. Using an array of adjustable diffusive fluxes, the standard…
In the case of hyperbolic conservation laws, high-order methods, such as the classical DG method, experience the phenomenon of unwanted high-frequency oscillations in the vicinity of a shock. Shock-capturing methods such as artificial…
We revise the structure-preserving finite element method in [K. Hu, Y. MA and J. Xu. (2017) Stable finite element methods preserving $\nabla \cdot \mathbf{B}=0$ exactly for MHD models. Numer. Math.,135, 371-396]. The revised method is…
The Allen-Cahn equation (ACE) inherently possesses two crucial properties: the maximum principle and the energy dissipation law. Preserving these two properties at the discrete level is also necessary in the numerical methods for the ACE.…
In this paper we construct a third order method for solving additively split autonomous stiff systems of ordinary differential equations. The constructed additive method is L-stable with respect to the implicit part and allows to use an…
Under non-global Lipschitz condition, Euler Explicit method fails to converge strongly to the exact solution, while Euler implicit method converges but requires much computational efforts. Tamed scheme was first introduced in [2] to…
We develop local discontinuous Galerkin (LDG) methods for conservation laws with heterogeneous stochastic fluxes, where the Stratonovich-driven transport terms may be linear or nonlinear. Such equations arise, for example, in simplified…
We construct a finite element discretization and time-stepping scheme for the incompressible Euler equations with variable density that exactly preserves total mass, total squared density, total energy, and pointwise incompressibility. The…
This study aims to construct a stable, high-order compact finite difference method for solving Sobolev-type equations with Dirichlet boundary conditions in one-space dimension. Approximation of higher-order mixed derivatives in some…
The Residual Smooting Scheme (RSS) have been introduced in \cite{AverbuchCohenIsraeli} as a backward Euler's method with a simplified implicit part for the solution of parabolic problems. RSS have stability properties comparable to those of…
We consider the nonlinear Schr{\"o}dinger equation (NLSE) in 1+1 dimension with scalar-scalar self interaction $\frac{g^2}{\kappa+1} (\psi^\star \psi)^{\kappa+1}$ in the presence of the external forcing terms of the form $r e^{-i(kx +…
In this paper, we introduce both monotone and nonmonotone variants of LiBCoD, a \textbf{Li}nearized \textbf{B}lock \textbf{Co}ordinate \textbf{D}escent method for solving composite optimization problems. At each iteration, a random block is…
The numerical solution of problems in nonlinear magnetostatics is typically based on a variational formulation in terms of magnetic potentials, the discretization by finite elements, and iterative solvers like the Newton method. The vector…
We propose a novel formulation for parametric finite element methods to simulate surface diffusion of closed curves, which is also called as the curve diffusion. Several high-order temporal discretizations are proposed based on this new…