Related papers: A novel directly energy-preserving method for char…
In this work, we develop a novel numerical scheme to solve the classical Keller--Segel (KS) model which simultaneously preserves its intrinsic mathematical structure and achieves optimal accuracy. The model is reformulated into a gradient…
We present a compact discontinuous Galerkin (CDG) method for an elliptic model problem. The problem is first cast as a system of first order equations by introducing the gradient of the primal unknown, or flux, as an additional variable. A…
In this paper, we improve the original Lagrange multiplier approach \cite{ChSh22,ChSh_II22} and introduce a new energy correction approach to construct a class of robust, positivity/bound-preserving, mass conserving and energy dissipative…
In this paper, we present a shock capturing discontinuous Galerkin (SC-DG) method for nonlinear systems of conservation laws in several space dimensions and analyze its stability and convergence. The scheme is realized as a space-time…
Semi-implicit multilevel spectral deferred correction (SI-MLSDC) methods provide a promising approach for high-order time integration for nonlinear evolution equations including conservation laws. However, existing methods lack robustness…
In this paper, we introduce and analyze a class of numerical schemes that demonstrate remarkable superiority in terms of efficiency, the preservation of positivity, energy stability, and high-order precision to solve the time-dependent…
Many practical problems can be described by second-order system $\ddot{q}=-M\nabla U(q)$, in which people give special emphasis to some invariants with explicit physical meaning, such as energy, momentum, angular momentum, etc. However,…
A hybrid method for the incompressible Navier--Stokes equations is presented. The method inherits the attractive stabilizing mechanism of upwinded discontinuous Galerkin methods when momentum advection becomes significant, equal-order…
This article presents a new numerical scheme for the discretization of dissipative particle dynamics with conserved energy. The key idea is to reduce elementary pairwise stochastic dynamics (either fluctuation/dissipation or thermal…
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…
We report a new particle in cell (PIC) method based on the semi-implicit approach. The novelty of the new method is that unlike any of its semi-implicit predecessors at the same time retains the explicit computational cycle and conserves…
This paper proposes a numerical method, termed PosDiv-CDG, that provably preserves both positivity and the globally divergence-free (DF) condition at arbitrarily high order in multiple dimensions. It resolves the fundamental structural…
We develop a Discrete Element Method (DEM) for elastodynamics using polyhedral elements. We show that for a given choice of forces and torques, we recover the equations of linear elastodynamics in small deformations. Furthermore, the…
The Energy Conserving Semi-Implicit Method (ECSIM) introduced by Lapenta (2017) has many advantageous properties compared to the classical semi-implicit and explicit PIC methods. Most importantly, energy conservation eliminates the growth…
We develop an efficient algorithm, coordinate descent FCI (CDFCI), for the electronic structure ground state calculation in the configuration interaction framework. CDFCI solves an unconstrained non-convex optimization problem, which is a…
In this paper, we derive a theoretical analysis of an interior penalty discontinuous Galerkin methods for solving the Cahn-Hilliard-Navier-Stokes model problem. We prove unconditional unique solvability of the discrete system, obtain…
In this work, a nonlinear momentum method is introduced to enhance the convergence performance of momentum-based gradient optimization algorithms. Classical momentum methods, such as the Heavy Ball method, can be viewed as a dynamical…
Many modern discontinuous Galerkin (DG) methods for conservation laws make use of summation by parts operators and flux differencing to achieve kinetic energy preservation or entropy stability. While these techniques increase the robustness…
A method for constructing first integral preserving numerical schemes for time-dependent partial differential equations on non-uniform grids is presented. The method can be used with both finite difference and partition of unity approaches,…
We introduce a new hybridized discontinuous Galerkin method for the incompressible magnetohydrodynamics equations. If particular velocity, pressure, magnetic field, and magnetic pressure spaces are employed for both element and trace…