Related papers: A positivity-preserving adaptive-order finite-diff…
Numerical schemes provably preserving the positivity of density and pressure are highly desirable for MHD, but the rigorous positivity-preserving (PP) analysis remains challenging. The difficulties mainly arise from the intrinsic complexity…
This article considers a new discretization scheme for conservation laws. The discretization setting is based on a discontinuous Galerkin scheme in combination with an approximation space that contains high-order polynomial modes as well as…
The relativistic magnetohydrodynamics (RMHD) set of equations has recently seen increased use in astrophysical computations. Even so, RMHD codes remain fragile. The reconstruction can sometimes yield superluminal velocities in certain parts…
We present a new computational method for smoothly matching general relativistic ideal magnetohydrodynamics (MHD) to its force-free limit. The method is based on a flux-conservative formalism for MHD and its force-free limit, and a vector…
In this paper, we design, analyze, and numerically validate positive and energy-dissipating schemes for solving the time-dependent multi-dimensional system of Poisson-Nernst-Planck (PNP) equations, which has found much use in the modeling…
Numerically solving magnetohydrodynamic (MHD) equations faces many challenges: avoiding divergence error, maintaining positivity, and satisfying entropy conditions. Among discontinuous Galerkin (DG) schemes, there has been a modal version…
This paper develops a family of fast, structure-preserving numerical algorithms for the nonlinear Maxwell-Ampere Nernst-Planck equations. For the first-order scheme, the Slotboom transformation rewrites the Nernst-Planck equation to enable…
We present a new numerical tool for solving the special relativistic ideal MHD equations that is based on the combination of the following three key features: (i) a one-step ADER discontinuous Galerkin (DG) scheme that allows for an…
Solving high-dimensional Fokker-Planck (FP) equations is a challenge in computational physics and stochastic dynamics, due to the curse of dimensionality (CoD) and unbounded domains. Existing deep learning approaches, such as…
In a recent work [10], we have introduced a pressure-robust Hybrid High-Order method for the numerical solution of the incompressible Navier-Stokes equations on matching simplicial meshes. Pressure-robust methods are characterized by error…
In ordinary turbulence research it has been a long standing tradition to solve the equations in spectral space giving the best possible accuracy. This is indeed a natural choice for incompressible problems with periodic boundaries, but it…
In this work we construct a low-order nonconforming approximation method for linear elasticity problems supporting general meshes and valid in two and three space dimensions. The method is obtained by hacking the Hybrid High-Order method,…
In this work, a new hybrid predictive Reduced Order Model (ROM) is proposed to solve reacting flow problems. This algorithm is based on a dimensionality reduction using Proper Orthogonal Decomposition (POD) combined with deep learning…
In this paper, we develop and present an arbitrary high order well-balanced finite volume WENO method combined with the modified Patankar Deferred Correction (mPDeC) time integration method for the shallow water equations. Due to the…
A second order accurate numerical scheme is proposed and analyzed for the periodic three-component Macromolecular Microsphere Composite(MMC) hydrogels system, a ternary Cahn-Hilliard system with a Flory-Huggins-deGennes free energy…
A novel approach for selecting appropriate reconstructions is implemented to the hyperbolic conservation laws in the high-order local polynomial-based framework, e.g., the discontinuous Galerkin (DG) and flux reconstruction (FR) schemes.…
We develop new numerical schemes for Vlasov--Poisson equations with high-order accuracy. Our methods are based on a spatially monotonicity-preserving (MP) scheme and are modified suitably so that positivity of the distribution function is…
A robust second order, shock-capturing numerical scheme for multi-dimensional special relativistic magnetohydrodynamics on computational domains with adaptive mesh refinement is presented. The base solver is a total variation diminishing…
A general method for constructing high order upwind schemes for multidimensional magnetohydrodynamics (MHD), having as a main built-in condition the divergence-free constraint $\divb=0$ for the magnetic field vector $\bb$, is proposed. The…
High-order reconstruction schemes for the solution of hyperbolic conservation laws in orthogonal curvilinear coordinates are revised in the finite volume approach. The formulation employs a piecewise polynomial approximation to the…