Related papers: Orbital Advection by Interpolation: A Fast and Acc…
In this paper, we propose high order numerical methods to solve a 2D advection diffusion equation, in the highly oscillatory regime. We use an integrator strategy that allows the construction of arbitrary high-order schemes {leading} to an…
The concept of an advection-dominated accretion flow (ADAF), with or without wind loss, was used to describe the spectra of the galactic center source Sgr A*, low-luminosity AGN and nuclei of elliptical galaxies (including M87). The…
We present a high-order radial basis function finite difference (RBF-FD) framework for the solution of advection-diffusion equations on time-varying domains. Our framework is based on a generalization of the recently developed Overlapped…
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…
Using only physical mechanisms, i.e., 3D MHD with no phenomenological viscosity, we have simulated the dynamics of a moderately thin accretion disk subject to torques whose radial scaling mimics those produced by lowest-order post-Newtonian…
The evolution of a large-scale poloidal magnetic field in accretion discs is an important problem because of its role in the launching of jets and winds and in determining the intensity of turbulence. In this paper, we develop a formalism…
We present a forward semi-Lagrangian numerical method for systems of transport equations able to advect smooth and discontinuous fields with high-order accuracy. The numerical scheme is composed of an integration of the transport equations…
We have developed a three-dimensional numerical model and applied it to simulate plasma flows in semi-detached binary systems whose accretor possesses a strong intrinsic magnetic field. The model is based on the assumption that the plasma…
We propose a novel numerical method for the solution of the shallow water equations in different regimes of the Froude number making use of general polygonal meshes. The fluxes of the governing equations are split such that advection and…
An algorithm for simulating self-gravitating cosmological astrophysical fluids is presented. The advantages include a large dynamic range, parallelizability, high resolution per grid element and fast execution speed. The code is based on a…
Due to the prevalence of magnetic fields in astrophysical environments, magnetohydrodynamic (MHD) simulation has become a basic tool for studying astrophysical fluid dynamics. To further advance the precision of MHD simulations, we have…
Modelling the flow in a thin accretion disc like a dynamical system, we analyse the nature of the critical points of the steady solutions of the flow. For the simple inviscid disc there are two critical points, with the outer one being a…
We analyze the mixed frame equations of radiation hydrodynamics under the approximations of flux-limited diffusion and a thermal radiation field, and derive the minimal set of evolution equations that includes all terms that are of leading…
A new moving mesh scheme based on the Lagrange-Galerkin method for the approximation of the one-dimensional convection-diffusion equation is studied. The mesh movement, which is prescribed by a discretized dynamical system for the nodal…
Vortices, turbulence, and unsteady non-laminar flows are likely both prominent and dynamically important features of astrophysical disks. Such strongly nonlinear phenomena are often difficult, however, to simulate accurately, and are…
The migration of growing protoplanets depends on the thermodynamics of the ambient disc. Standard modelling, using locally isothermal discs, indicate in the low planet mass regime an inward (type-I) migration. Taking into account…
This article presents a high order conservative flux optimization (CFO) finite element method for the elliptic diffusion equations. The numerical scheme is based on the classical Galerkin finite element method enhanced by a flux…
We present an efficient dimension-by-dimension finite-volume method which solves the adiabatic magnetohydrodynamics equations at high discretization order, using the constrained-transport approach on Cartesian grids. Results are presented…
We study supersonic flow past a convex corner which is surrounded by quiescent gas. When the pressure of the upstream supersonic flow is larger than that of the quiescent gas, there appears a strong rarefaction wave to rarefy the supersonic…
A Cartesian grid method combined with a simplified gas kinetic scheme is presented for subsonic and supersonic viscous flow simulation on complex geometries. Under the Cartesian mesh, the computational grid points are classified into four…