Related papers: A direct Eulerian GRP scheme for spherically symme…
This paper presents a fully discrete numerical scheme for one-dimensional nonlocal wave equations and provides a rigorous theoretical analysis. To facilitate the spatial discretization, we introduce an auxiliary variable analogous to the…
We construct centered rarefaction wave solutions given background solutions to the compressible Euler equations. The flow considered in this article is the homentropic flow of perfect gas governed by compressible Euler equations and the…
We consider the relativistic Euler equations governing spherically symmetric, perfect fluid flows on the outer domain of communication of Schwarzschild spacetime, and we introduce a version of the finite volume method which is formulated…
We present the new open-source spherically-symmetric general-relativistic (GR) hydrodynamics code GR1D. It is based on the Eulerian formulation of GR hydrodynamics (GRHD) put forth by Romero-Ibanez-Gourgoulhon and employs radial-gauge,…
We consider controlled differential equations and give new estimates for higher order Euler schemes. Our proofs are inspired by recent work of A. M. Davie who considers first and second order schemes. In order to implement the general case…
An efficient algorithm for solving Poisson's equation in two and three spatial dimensions is discussed. The algorithm, which is described in detail, is based on the integral form of Poisson's equation and utilizes spherical coordinates and…
We present the derivation of hydrodynamical equations for a perfect fluid in General Relativity, within the 3+1 decomposition of spacetime framework, using only primitive variables. Primitive variables are opposed to conserved variables, as…
We present a Gaussian Process (GP) approach (Gaussian Process Hydrodynamics, GPH) for approximating the solution of the Euler and Navier-Stokes equations. As in Smoothed Particle Hydrodynamics (SPH), GPH is a Lagrangian particle-based…
This paper presents entropy symmetrization and high-order accurate entropy stable schemes for the relativistic magnetohydrodynamic (RMHD) equations. It is shown that the conservative RMHD equations are not symmetrizable and do not admit a…
The theory of Wasserstein gradient flows in the space of probability measures has made an enormous progress over the last twenty years. It constitutes a unified and powerful framework in the study of dissipative partial differential…
This paper studies high-order accurate entropy stable nodal discontinuous Galerkin (DG) schemes for the ideal special relativistic magnetohydrodynamics (RMHD). It is built on the modified RMHD equations with a particular source term, which…
In a recent paper we presented a new ultra efficient numerical method for solving kinetic equations of the Boltzmann type (G. Dimarco, R. Loubere, Towards an ultra efficient kinetic scheme. Part I: basics on the 689 BGK equation, J. Comp.…
This paper is concerned with the numerical approximation of the isothermal Euler equations for charged particles subject to the Lorentz force. When the magnetic field is large, the so-called drift-fluid approximation is obtained. In this…
In the computation of compressible fluid flows, numerical boundary conditions are always necessary for all physical variables at computational boundaries while just partial physical variables are often prescribed as physical boundary…
In this article we develop a Primitive Variable Recovery Scheme (PVRS) to solve any system of coupled differential conservative equations. This method obtains directly the primitive variables applying the chain rule to the time term of the…
High order algorithms have emerged in numerical astrophysics as a promising avenue to reduce truncation error (proportional to a power of the linear resolution $\Delta x$) with only a moderate increase to computational expense. Significant…
By a semi-Lagrangian change of coordinates, the hydrostatic Euler equations describing free-surface sheared flows is rewritten as a system of quasilinear equations, where stability conditions can be determined by the analysis of its…
We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply-connected fluid domain that includes…
There has been increasing interest in statistical analysis of data lying in manifolds. This paper generalizes a smoothing spline fitting method to Riemannian manifold data based on the technique of unrolling and unwrapping originally…
We present a second-order strictly length-preserving and unconditionally energy-stable rotational discrete gradient (Rdg) scheme for the numerical approximation of the Oseen-Frank gradient flows with anisotropic elastic energy functional.…