Related papers: 1+1D implicit disk computations
The planet-forming region of protoplanetary disks is cold, dense, and therefore weakly ionized. For this reason, magnetohydrodynamic (MHD) turbulence is thought to be mostly absent, and another mechanism has to be found to explain gas…
The paper develops high-order physical-constraint-preserving (PCP) methods for general relativistic hydrodynamic (GRHD) equations, equipped with a general equation of state. Here the physical constraints, describing the admissible states of…
Turbulent accretion disks threaded by a large-scale vertical field near equipartition can drive tenuous and fast self-confined jets. Self-similar solutions of these jet-emitting disks (JEDs) have been known for a long time and provide the…
Time-implicit schemes are attractive since they allow numerical time steps that are much larger than those permitted by the Courant-Friedrich-Lewy criterion characterizing time-explicit methods. This advantage comes, however, with a cost:…
Semi-implicit methods are powerful and efficient tools for the three-dimensional modeling of coastal and oceanic processes. A semi-implicit finite difference method for 3D hydrostatic primitive equations is presented in this paper. The…
The local axisymmetric stability of hydrodynamical and magnetized, nearly-Keplerian gaseous accretion disks around non-rotating black holes is examined in the vicinity of the classical marginally-stable orbit (at radii ~ R_ms). An…
We construct both aligned and unaligned (logarithmic spiral) stationary configurations of nonaxisymmetric magnetohydrodynamic (MHD) disks from either a full or a partial razor-thin power-law axisymmetric magnetized singular isothermal disk…
We study 1+1 dimensional relativistic non-resistive magnetohydrodynamics (MHD) with longitudinal boost invariance and shear stress tensor. Several analytical solutions that describe the fluid temperature evolution under the equation of…
One-dimensional integrable and quasi-integrable systems display, on macroscopic scales, a universal form of transport known as Generalized Hydrodynamics (GHD). In its standard Euler-scale formulation, GHD mirrors the equations of a…
We propose a new method to solve the relativistic hydrodynamic equations based on implicit Runge-Kutta methods with a locally optimized fixed-point iterative solver. For numerical demonstration, we implement our idea for ideal hydrodynamics…
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…
A dispersion relation for gravity waves in water covered by disk-like impurities embedded in a viscous matrix is derived. The macroscopic equations are obtained by ensemble-averaging the fluid equations at the disk scale in the asymptotic…
We propose an easy-to-implement iterative method for resolving the implicit (or semi-implicit) schemes arising in solving reaction-diffusion (RD) type equations. We formulate the nonlinear time implicit scheme as a min-max saddle point…
We describe a novel Godunov-type numerical method for solving the equations of resistive relativistic magnetohydrodynamics. In the proposed approach, the spatial components of both magnetic and electric fields are located at zone interfaces…
Context. In Class 0/I and the outskirts of Class II circumstellar discs, the self-gravity of gas significantly affects the disc's vertical hydrostatic equilibrium. The contribution of dust, whose measured mass is still uncertain, could…
We find exact and explicit solutions of the axisymmetric MHD equations of a self-gravitating polytropic gas. These solutions are able to describe a flat (uniform density) subsonic internal core, contracting homologously, of a collapsing…
We present results from three-dimensional, self-gravitating radiation hydrodynamical models of gas accretion by planetary cores. In some cases, the accretion flow is resolved down to the surface of the solid core -- the first time such…
(shortened) We calculate the vertical structure of a local patch of an accretion disk in which heating by dissipation of MRI-driven MHD turbulence is balanced by radiative cooling. Heating, radiative transport, and cooling are computed…
Observations of protoplanetary disks show that some characteristics seem recurrent, even in star formation regions that are physically distant such as surface mass density profiles varying as $r^{-1}$, or aspect ratios about 0.03 to 0.23.…
We present a geometric multigrid solver for the M1 model of radiative transfer without source terms. In radiative hydrodynamics applications, the radiative transfer needs to be solved implicitly because of the fast propagation speed of…