Related papers: Constrained Hyperbolic Divergence Cleaning for Smo…
We present a fully Lagrangian conservation form of the general relativistic hydrodynamic equations for perfect fluids with artificial viscosity in a given arbitrary background spacetime. This conservation formulation is achieved by choosing…
We develop a family of stabilized backward differentiation formula (sBDF) schemes of orders one through four for semilinear parabolic equations. The proposed methods are designed to achieve three properties that are rarely available…
Hydrodynamical simulations are the most accurate way to model structure formation in the universe, but they often involve a large number of astrophysical parameters modeling subgrid physics, in addition to cosmological parameters. This…
Approximate Riemann solvers are widely used for solving hyperbolic conservation laws, including those of magnetohydrodynamics (MHD). However, due to the nonlinearity and complexity of MHD, obtaining accurate and robust numerical solutions…
Obtainable computational efficiency is evaluated when using an Adaptive Mesh Refinement (AMR) strategy in time accurate simulations governed by sets of conservation laws. For a variety of 1D, 2D, and 3D hydro- and magnetohydrodynamic…
Central schemes for conservation laws are Riemann solver free methods which are simple and easy to implement. In recent work for Euler equations [Kurganov & Xin, J. Sci. Comput., 96:56, 2023] their accuracy has been enhanced in terms of…
We have written and tested a new general relativistic magnetohydrodynamics (GRMHD) code, capable of evolving MHD fluids in dynamical spacetimes with adaptive-mesh refinement (AMR). Our code solves the Einstein-Maxwell-MHD system of coupled…
We have studied forced turbulence of compressible magnetohydrodynamic (MHD) flows through two-dimensional simulations with different numerical resolutions. First, hydrodynamic turbulence with Mach number $<M_s >_{\rm init} \equiv < v >_{\rm…
A third order shock-capturing numerical scheme for three-dimensional special relativistic magnetohydrodynamics (3-D RMHD) is presented and validated against several numerical tests. The simple and efficient central scheme described in Paper…
This paper describes an energy-preserving and globally time-reversible code for weakly compressible smoothed particle hydrodynamics (SPH). We do not add any additional dynamics to the Monaghan's original SPH scheme at the level of ordinary…
In this paper we propose and analyze a mixed DG method and an HDG method for the stationary Magnetohydrodynamics (MHD) equations with two types of boundary (or constraint) conditions. The mixed DG method is based a recent work proposed by…
The free-boundary compressible 1-D Euler equations with moving physical vacuum boundary are a system of hyperbolic conservation laws which are both characteristic and degenerate. The physical vacuum singularity (or rate-of-degeneracy)…
Motivated by a geometric method employed for the derivation of the Nambu bracket for ideal two-dimensional incompressible hydrodynamics, we reconstruct the reduced magnetohydrodynamic (RMHD) model by a priori imposition of its conservation…
Analysing data from Smoothed Particle Hydrodynamics (SPH) simulations is about understanding global fluid properties rather than individual fluid elements. Therefore, in order to properly understand the outcome of such simulations it is…
I discuss the key features of Smoothed Particle Hydrodynamics (SPH) as a numerical method - in particular the key differences between SPH and more standard grid based approaches - that are important to the practitioner. These include the…
In order to solve the magnetohydrodynamics (MHD) equations with a $\mathbf{\mathcal{H}}(\mathbf{div})$-conforming element, a novel approach is proposed to ensure the exact divergence-free condition on the magnetic field. The idea is to add…
In astrophysical magnetohydrodynamics (MHD) and electrodynamics simulations, numerically enforcing the divB=0 constraint on the magnetic field has been difficult. We observe that for point-based discretization, as used in finite-difference…
A new approach using a hyperbolic-equation system (HES) is proposed to solve for the electron fluids in quasi-neutral plasmas. The HES approach avoids treatments of cross-diffusion terms which cause numerical instabilities in conventional…
An extension of the H-theorem for dissipative particle dynamics (DPD) to the case of a multi-component fluid is made. Detailed balance and an additional H-theorem are proved for an energy-conserving version of the DPD algorithm. The…
In this work, we present a positivity-preserving adaptive filtering approach for discontinuous spectral element approximations of the ideal magnetohydrodynamics equations. This approach combines the entropy filtering method (Dzanic and…