Related papers: Mitigating numerical dissipation in simulations of…
Based on a 14-moment closure for non-resistive (general-) relativistic viscous plasmas, we describe a new numerical scheme that is able to handle all first-order dissipative effects (heat conduction, bulk and shear viscosities), as well the…
This paper aims to numerically verify the large Reynolds number asymptotic theory of magneto-hydrodynamic (MHD) flows proposed in the companion paper Deguchi (2019). To avoid any complexity associated with the chaotic nature of turbulence…
We live in an age in which high-performance computing is transforming the way we do science. Previously intractable problems are now becoming accessible by means of increasingly realistic numerical simulations. One of the most enduring and…
In this work we present a proof of concept of CUDA-capable, resistive, multi-fluid models of relativistic magnetohydrodynamics (RMHD). Resistive and multi-fluid codes for simulating models of RMHD suffer from stiff source terms, so it is…
Observations, mainly of outbursts in dwarf novae, imply that the anomalous viscosity in highly ionized accretion discs is magnetic in origin, and requires that the plasma $\beta \sim 1$. Until now most simulations of the magnetic dynamo in…
Aims. The main goal of the present paper is to provide the first systematic numerical study of the propagation of astrophysical relativistic jets, in the context of high-resolution shock-capturing resistive relativistic magnetohydrodynamics…
Weakly-collisional plasmas, such as the solar wind or the intra-cluster medium (ICM) of galaxy clusters, evolve in the presence of dynamically strong magnetic fields. The turbulent dynamo can amplify magnetic fields to such levels by…
We present a new numerical method of special relativistic resistive magnetohydrodynamics with scalar resistivity that can treat a range of phenomena, from nonrelativistic to relativistic (shock, contact discontinuity, and Alfv\'en wave).…
The performance of high-order gas-kinetic scheme (HGKS) has been investigated for the direct numerical simulation (DNS) of isotropic compressible turbulence up to the supersonic regime. Due to the multi-scale nature and coupled…
Direct numerical simulations (DNS) are an indispensable tool for understanding the fundamental physics of turbulent flows. Because of their steep increase in computational cost with Reynolds number ($R_{\lambda}$), well-resolved DNS are…
In numerical magnetohydrodynamics (MHD), a major challenge is maintaining zero magnetic field-divergence (div-B). Constrained transport (CT) schemes can achieve this at high accuracy, but have generally been restricted to very specific…
We present pseudospectral direct-numerical-simulation (DNS) studies of the three-dimensional magnetohydrodynamic (MHD) equations (3DRFMHD) with a stochastic force that has zero mean and a variance $\sim k^{-3}$, where $k$ is the wavenumber,…
Numerical simulations of magneto-convection have greatly expanded our understanding of stellar interiors and stellar magnetism. Recently, fully compressible hydrodynamical simulations of full-star models have demonstrated the feasibility of…
In this paper, a unified gas-kinetic scheme (UGKS) with simplified multi-scale numerical flux is proposed for the thermodynamic non-equilibrium flow simulation involving the excitation of molecular vibrational degrees of freedom in all flow…
The discovery of dynamical models from data represents a crucial step in advancing our understanding of physical systems. Library-based sparse regression has emerged as a powerful method for inferring governing equations directly from…
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…
The paper presents two contributions in the context of the numerical simulation of magnetized fluid dynamics. First, we show how to extend the ideal magnetohydrodynamics (MHD) equations with an inbuilt magnetic field divergence cleaning…
Many astrophysical small-scale dynamos (SSDs) amplify weak magnetic fields via highly compressible, supersonic turbulence, but established SSD theories have overlooked these compressible effects. To address this, we perform visco-resistive…
We study numerical methods for dissipative particle dynamics (DPD), which is a system of stochastic differential equations and a popular stochastic momentum-conserving thermostat for simulating complex hydrodynamic behavior at mesoscales.…
We present in this work a new reconstruction scheme, so-called MUSCL-THINC-BVD scheme, to solve the five-equation model for interfacial two phase flows. This scheme employs the traditional shock capturing MUSCL (Monotone Upstream-centered…