Related papers: A Multi-Fluid Dust Module in Athena++: Algorithms …
A new code and methodology are introduced for solving the general relativistic magnetohydrodynamic (GRMHD) equations in fixed background spacetimes using time-explicit, finite-volume discretization. The code has options for solving the…
Forward modeling is often used to interpret substructures observed in protoplanetary disks. To ensure the robustness and consistency of the current forward modeling approach from the community, we conducted a systematic comparison of…
We present a detailed study of localised magnetohydrodynamical (MHD) instabilities occuring in two--dimensional magnetized accretion disks. We model axisymmetric MHD disk tori, and solve the equations governing a two--dimensional magnetized…
In this chapter, we aim at presenting the basic techniques necessary to go beyond the widely accepted paradigm of second-order numerics. We specifically focus on finite-volume schemes for hyperbolic conservation laws occuring in fluid…
Simulating and predicting dynamics of quantum many-body systems is extremely challenging, even for state-of-the-art computational methods, due to the spread of entanglement across the system. However, in the long-wavelength limit, quantum…
Magnetohydrodynamics (MHD) describes the interaction between electrically conducting fluids and electromagnetic fields. We propose and analyze a symplectic, second-order algorithm for the evolutionary MHD system in Els\"asser variables. We…
We present numerical simulations and explore scalings and anisotropy of compressible magnetohydrodynamic (MHD) turbulence. Our study covers both gas pressure dominated (high beta) and magnetically dominated (low beta) plasmas at different…
This work presents the publicly available moving-mesh magnetohydrodynamics code DISCO. DISCO is efficient and accurate at evolving orbital fluid motion in two and three dimensions, especially at high Mach number. DISCO employs a moving-mesh…
A previously-developed hybrid particle-continuum method [J. B. Bell, A. Garcia and S. A. Williams, SIAM Multiscale Modeling and Simulation, 6:1256-1280, 2008] is generalized to dense fluids and two and three dimensional flows. The scheme…
This paper describes a numerical scheme for multi-fluid hydrodynamics in the limit of small mass densities of the charged particles. The inertia of the charged particles can then be neglected, which makes it possible to write an evolution…
Numerical simulations with self-similar initial and boundary conditions provide a link between theoretical and numerical investigations of jet dynamics. We perform axisymmetric resistive magnetohydrodynamic (MHD) simulations for a…
(abridged) We apply a second-order Godunov code, Athena, to studies of the magnetorotational instability using unstratified shearing box simulations with a uniform net vertical field and a sinusoidally varying zero net vertical field. The…
We present a novel combination of numerical techniques to improve the efficiency, accuracy, and robustness of multi-component compressible flow simulations. At the core of our approach is an Entropy-Stable formulation that preserves kinetic…
Multiphase gas, with hot ($\sim10^6$K) and cold ($\sim10^4$K) gas, is ubiquitous in astrophysical media across a wide range of scales. However, simulating multiphase gas has been a long-standing challenge, due to the large separation…
Circumstellar discs, from which planetary systems are formed, consist of gas, dust and solids. Simulations of self-consistent dynamics of gas, dust and solids in circumstellar discs is a challenging problem. In the paper we present fast…
We introduce TIGRESS, a novel framework for multi-physics numerical simulations of the star-forming interstellar medium (ISM) implemented in the Athena MHD code. The algorithms of TIGRESS are designed to spatially and temporally resolve key…
The present work is devoted to introduce the backward Euler based modular time filter method for MHD flow. The proposed method improves the accuracy of the solution without a significant change in the complexity of the system. Since time…
Numerical solutions of the incompressible magnetohydrodynamic (MHD) equations are reported for the interior of a rotating, perfectly-conducting, rigid spherical shell that is insulator-coated on the inside. A previously-reported spectral…
The Kelvin-Helmholtz instability (KHI) is a standard test of hydrodynamic and magnetohydrodynamic (MHD) simulation codes and finds many applications in astrophysics. The classic linear theory considers a discontinuity in density and…
This paper introduces the Sheffield Magnetohydrodynamics Algorithm Using GPUs (SMAUG+), an advanced numerical code for solving magnetohydrodynamic (MHD) problems, using multi-GPU systems. Multi-GPU systems facilitate the development of…