Related papers: A Multistate Low-dissipation Advection Upstream Sp…
We present the FARGO3D code, recently publicly released. It is a magnetohydrodynamics code developed with special emphasis on protoplanetary disks physics and planet-disk interactions, and parallelized with MPI. The hydrodynamics algorithms…
We present a scalable, high-order implicit large-eddy simulation (ILES) approach for incompressible transitional flows. This method employs the mass-conserving mixed stress (MCS) method for discretizing the Navier-Stokes equations. The MCS…
Modern astrophysical simulations aim to accurately model an ever-growing array of physical processes, including the interaction of fluids with magnetic fields, under increasingly stringent performance and scalability requirements driven by…
Liou-Steffen splitting (AUSM) schemes are popular for low Mach number simulations, however, like many numerical schemes for compressible flow they require careful modification to accurately resolve convective features in this regime.…
We propose two sets of initial conditions for magnetohydrodynamics (MHD) in which both the velocity and the magnetic fields have spatial symmetries that are preserved by the dynamical equations as the system evolves. When implemented…
Entropy stable schemes can be constructed with a specific choice of the numerical flux function. First, an entropy conserving flux is constructed. Secondly, an entropy stable dissipation term is added to this flux to guarantee dissipation…
Robust and accurate fully implicit finite-volume schemes applied to Darcy-scale multiphase flow and transport in porous media are highly desirable. Recently, a smooth approximation of the saturation-dependent flux coefficients based on…
Accurate simulations of flows in stellar interiors are crucial to improving our understanding of stellar structure and evolution. Because the typically slow flows are merely tiny perturbations on top of a close balance between gravity and…
In this paper, we consider numerical approximations for solving the nonlinear magneto-hydrodynamical system, that couples the Navier-Stokes equations and Maxwell equations together. A challenging issue to solve this model numerically is…
This study proposes an algorithm for modeling compressible flows in spherical shells in nearly incompressible and weakly compressible regimes based on an implicit direction splitting approach. The method retains theoretically expected…
The moment-of-fluid method (MOF) is known as an extension of the volume-of-fluid method with piecewise linear interface construction (VOF-PLIC). In this study, several directional splitting advection algorithms are extended from the…
We present THC: a new high-order flux-vector-splitting code for Newtonian and special-relativistic hydrodynamics designed for direct numerical simulations of turbulent flows. Our code implements a variety of different reconstruction…
The majority of available numerical algorithms for interfacial two-phase flows either treat both fluid phases as incompressible (constant density) or treat both phases as compressible (variable density). This presents a limitation for the…
This paper presents an extension of the hybrid scheme proposed by Wang et al. (J. Comput. Phys. 229 (2010) 169-180) for numerical simulation of compressible isotropic turbulence to flows with higher turbulent Mach numbers. The scheme still…
A general procedure to construct a class of simple and efficient high resolution Total Variation Diminishing (TVD) schemes for non-linear hyperbolic conservation laws by introducing anti-diffusive terms with the flux limiters is presented.…
We consider the simulation of isentropic flow in pipelines and pipe networks. Standard operating conditions in pipe networks suggest an emphasis to simulate low Mach and high friction regimes -- however, the system is stiff in these regimes…
We describe a three-dimensional, Godunov-type numerical magnetohydrodynamics (MHD) method designed for studying disk accretion to a rotating magnetized star in the general case where the star's rotation axis, its magnetic moment, and the…
The HLLEM scheme is a popular contact and shear preserving approximate Riemann solver for cheap and accurate computation of high speed gasdynamical flows. Unfortunately this scheme is known to be plagued by various forms of numerical shock…
We present a numerical study of an idealized magnetohydrodynamic (MHD) configuration consisting of a planar wake flow embedded into a three-dimensional (3D) sheared magnetic field. Our simulations investigate the possibility for in-situ…
In order to prevent velocity, pressure, and temperature spikes at material discontinuities occurring when the interface-capturing schemes inconsistently simulate compressible multi-material flows(when the specific heats ratio is…