Related papers: WHAM: A WENO-based general relativistic numerical …
In this paper, we propose a well-balanced fifth-order finite difference Hermite WENO (HWENO) scheme for the shallow water equations with non-flat bottom topography in pre-balanced form. For achieving the well-balance property, we adopt the…
This paper deals with a new fifth-order weighted essentially non-oscillatory (WENO) scheme improving the WENO-NS and WENO-P methods which are introduced in Ha et al. J. Comput. Phys. (2013) and Kim et al., J. Sci. Comput. (2016)…
The ideal MHD equations are a central model in astrophysics, and their solution relies upon stable numerical schemes. We present an implementation of a new method, which possesses excellent stability properties. Numerical tests demonstrate…
Astrophysical relativistic flow problems require high resolution three-dimensional numerical simulations. In this paper, we describe a new parallel three-dimensional code for simulations of special relativistic hydrodynamics (SRHD) using…
Magnetohydrodynamic (MHD) simulations of subsonic (Mach number~$<1$) turbulence are crucial to our understanding of several processes including oceanic and atmospheric flows, the amplification of magnetic fields in the early universe,…
For more than 40 years the quest to understand how large-scale magnetic fields emerge from turbulent flows in rotating astrophysical systems, such as the Sun, has been a major focus of computational astrophysics research. Using a parameter…
Many problems at the forefront of theoretical astrophysics require the treatment of magnetized fluids in dynamical, strongly curved spacetimes. Such problems include the origin of gamma-ray bursts, magnetic braking of differential rotation…
We present a new version of conservative ADER-WENO finite volume schemes, in which both the high order spatial reconstruction as well as the time evolution of the reconstruction polynomials in the local space-time predictor stage are…
In this paper we analyze the weighted essentially non-oscillatory (WENO) schemes in the finite volume framework by examining the first step of the explicit third-order total variation diminishing Runge-Kutta method. The rationale for the…
Weighted essentially non-oscillatory (WENO) reconstruction schemes are presented that preserve cylindrical symmetry for radial flows on an equal-angle polar mesh. These new WENO schemes are used with a Lagrangian discontinuous Galerkin (DG)…
In this paper, a simple and efficient third-order weighted essentially non-oscillatory (WENO) reconstruction is developed for three-dimensional flows, in which the idea of two-dimensional WENO-AO scheme on unstructured meshes…
Methods to solve the relativistic hydrodynamic equations are a key computational kernel in a large number of astrophysics simulations and are crucial to understanding the electromagnetic signals that originate from the merger of…
KHARMA (an acronym for "Kokkos-based High-Accuracy Relativistic Magnetohydrodynamics with Adaptive mesh refinement") is a new open-source code for conducting general-relativistic magnetohydrodynamic simulations in stationary spacetimes,…
The shock instability problem commonly arises in flow simulations involving strong shocks, particularly when employing high-order schemes, limiting their applications in hypersonic flow simulations. This study focuses on exploring the…
In galaxy clusters, modern radio interferometers observe non-thermal radio sources with unprecedented spatial and spectral resolution. For the first time, the new data allows to infer the structure of the intra-cluster magnetic fields on…
Based on the understandings regarding linear upwind schemes with flux splitting to achieve free-stream preservation (Q. Li, etc. Commun. Comput. Phys., 22 (2017) 64-94), a series of WENO interpolation-based and upwind-biased nonlinear…
Accurate modeling of sea ice dynamics is critical for predicting environmental variables and is important in applications such as navigating ice breaker ships. Research for both modeling and simulating sea ice dynamics is ongoing, with the…
We describe a new hybrid framework to model non-thermal spectral signatures from highly energetic particles embedded in a large-scale classical or relativistic MHD flow. Our method makes use of \textit{Lagrangian} particles moving through…
We report a new finite-difference scheme for two-dimensional magnetohydrodynamics (MHD) simulations, with and without rotation, in unstructured grids with quadrilateral cells. The new scheme is implemented within the code VULCAN/2D, which…
Many systems of current interest in relativistic astrophysics require a knowledge of radiative transfer in a magnetized gas flowing in a strongly-curved, dynamical spacetime. Such systems include coalescing compact binaries containing…