Related papers: A Shock-Patching Code for Ultra-Relativistic Fluid…
We investigate asymptotic convergence in the~$\Delta x \!\rightarrow\! 0$ limit as a tool for determining whether numerical computations involving shocks are accurate. We use one-dimensional operator-split finite-difference schemes for…
The simulation of large open water surface is challenging using a uniform volumetric discretization of the Navier-Stokes equations. Simulating water splashes near moving objects, which height field methods for water waves cannot capture,…
Simulating inhomogeneous flows with different characteristic scales in different coordinate directions using the collide-and-stream based lattice Boltzmann methods (LBM) can be accomplished efficiently using rectangular lattice grids. We…
In an attempt to investigate the structures of ultra-relativistic jets injected into the intracluster medium (ICM) and the associated flow dynamics, such as shocks, velocity shear, and turbulence, we have developed a new special…
A new multidimensional simulation code for relativistic two-fluid electrodynamics (RTFED) is described. The basic equations consist of the full set of Maxwell's equations coupled with relativistic hydrodynamic equations for separate two…
High Mach number shocks are ubiquitous in interstellar turbulence. The Pencil Code is particularly well suited to the study of magnetohydrodynamics in weakly compressible turbulence and the numerical investigation of dynamos because of its…
Astrophysical plasmas in relativistic spacetimes, such as black hole accretion flows, are often weakly collisional and require kinetic modeling to capture non-local transport and particle acceleration. However, the extreme scale separation…
Analysing two-dimensional shallow water equations with idealised bottom topographies have many applications in the atmospheric and oceanic sciences; however, restrictive flow pattern assumptions have been made to achieve explicit solutions.…
Relativistic shocks are present in all high-energy astrophysical processes involving relativistic plasma outflows interacting with their ambient medium. While a well understood process in the context of relativistic hydrodynamics and ideal…
We propose a novel numerical method for solving multi-dimensional, special relativistic Boltzmann equations for neutrinos coupled to hydrodynamics equations. It is meant to be applied to simulations of core-collapse supernovae. We handle…
We present a new approximate Riemann solver for the augmented system of equations of resistive relativistic magnetohydrodynamics (RRMHD) that belongs to the family of Harten-Lax-van Leer contact wave (HLLC) solvers. In HLLC solvers, the…
We present a new special relativistic hydrodynamics (SRHD) code capable of handling coexisting ultra-relativistically hot and non-relativistically cold gases. We achieve this by designing a new algorithm for conversion between primitive and…
We survey a number of moment hierarchies and test their performances in computing one-dimensional shock structures. It is found that for high Mach numbers, the moment hierarchies are either computationally expensive or hard to converge,…
High-speed turbulent flows are encountered in most space-related applications (including exploration, tourism and defense fields) and represent a subject of growing interest in the last decades. A major challenge in performing high-fidelity…
Some of the most interesting scenarios that can be studied in astrophysics, contain fluids and plasma moving under the influence of strong gravitational fields. To study these problems it is required to implement numerical algorithms robust…
The lattice Boltzmann method (LBM) is routinely employed in the simulation of complex multiphase flows comprising bulk phases separated by non-ideal interfaces. LBM is intrinsically mesoscale with an hydro-dynamic equivalence popularly set…
This paper presents numerical cross-comparisons and benchmark results for two different kinetic numerical methods, capable of describing relativistic dissipative fluid dynamics in a wide range of kinematic regimes, typical of relevant…
We establish the hydrodynamic limit of the one-dimensional Boltzmann equation with hard-sphere collisions toward Riemann solutions of the compressible Euler system. The Riemann solutions covered by our result include generic superpositions…
Accurately resolving steady electrohydrodynamic (EHD) flows presents a formidable computational challenge due to the strong nonlinear coupling between charged-particle density, velocity fields, and electric potential. These interactions…
A D2Q9 Hybrid Lattice Boltzmann Method (HLBM) is proposed for the simulation of both compressible subsonic and supersonic flows. This HLBM is an extension of the model of Feng et al: [12], which has been found, via different test cases, to…