Related papers: A positivity-preserving adaptive-order finite-diff…
Conservative numerical schemes for general relativistic magnetohydrodynamics (GRMHD) require a method for transforming between ``conserved'' variables such as momentum and energy density and ``primitive'' variables such as rest-mass…
We propose an adaptive stencil construction for high order accurate finite volume schemes aposteriori stabilized devoted to solve one-dimensional steady-state hyperbolic equations. High-accuracy (up to the sixth-order presently) is achieved…
We present and compare third- as well as fifth-order accurate finite difference schemes for the numerical solution of the compressible ideal MHD equations in multiple spatial dimensions. The selected methods lean on four different…
We consider a class of finite element approximations for fourth-order parabolic equations that can be written as a system of second-order equations by introducing an auxiliary variable. In our approach, we first solve a variational problem…
The paper develops high-order accurate physical-constraints-preserving finite difference WENO schemes for special relativistic hydrodynamical (RHD) equations, built on the local Lax-Friedrich splitting, the WENO reconstruction, the…
Non-local systems of conservation laws play a crucial role in modeling flow mechanisms across various scenarios. The well-posedness of such problems is typically established by demonstrating the convergence of robust first-order schemes.…
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…
In this paper, we propose a novel approximation strategy for time-dependent hyperbolic systems of conservation laws for the Euler system of gas dynamics that aims to represent the dynamics of strong interacting discontinuities. The goal of…
The compact finite difference method is a powerful tool for discretizing conservation laws, owing to its inherent flexibility in developing high-resolution and highly stable schemes. In this paper, we propose a framework for the design of…
In this paper we propose the first better than second order accurate method in space and time for the numerical solution of the resistive relativistic magnetohydrodynamics (RRMHD) equations on unstructured meshes in multiple space…
A new code and methodology are introduced for solving the general relativistic magnetohydrodynamic (GRMHD) equations in fixed background spacetimes using time-explicit, finite-volume discretization. The code has options for solving the…
The paper describes a new upwind conservative numerical scheme for special relativistic resistive magnetohydrodynamics with scalar resistivity. The magnetic field is kept approximately divergence free and the divergence of the electric…
Model order reduction provides low-complexity high-fidelity surrogate models that allow rapid and accurate solutions of parametric differential equations. The development of reduced order models for parametric \emph{nonlinear} Hamiltonian…
We construct high order symmetric volume-preserving methods for the relativistic dynamics of a charged particle by the splitting technique with processing. Via expanding the phase space to include time $t$, we give a more general…
We present a high-order conservative, positivity-preserving, and non-oscillatory scheme for solving the Vlasov equation. The scheme attains formal fifth-order accuracy through a convex combination of positive and non-oscillatory polynomials…
We present the implementation of a three-dimensional, second order accurate Godunov-type algorithm for magneto-hydrodynamic (MHD), in the adaptive-mesh-refinement (AMR) cosmological code {\tt CHARM}. The algorithm is based on the full…
We describe a high-order numerical magnetohydrodynamics (MHD) solver built upon a novel non-linear entropy stable numerical flux function that supports eight travelling wave solutions. By construction the solver conserves mass, momentum,…
We propose and analyze a novel approach to construct structure preserving approximations for the Poisson-Nernst-Planck equations, focusing on the positivity preserving and mass conservation properties. The strategy consists of a standard…
We propose a conservative algorithm to test the geometrical validity of simplicial (triangles, tetrahedra), tensor product (quadrilaterals, hexahedra), and mixed (prisms) elements of arbitrary polynomial order as they deform over a…
For high-order accurate schemes such as discontinuous Galerkin (DG) methods solving Fokker-Planck equations, it is desired to efficiently enforce positivity without losing conservation and high-order accuracy, especially for implicit time…