Related papers: An Adaptive Grid, Implicit Code for Spherically Sy…
In this paper, we propose an adaptive high-order method for hyperbolic systems of conservation laws. The proposed method is based on a dual formulation approach: Two numerical solutions, corresponding to conservative and nonconservative…
Recent years have seen a significant progress in the development of general relativistic codes for the numerical solution of the equations of magnetohydrodynamics in spacetimes with high and dynamical curvature. These codes are valuable…
In this review we briefly summarize the so-called effective fluid approach, which is a compact framework that can be used to describe a plethora of different modified gravity models as general relativity (GR) and a dark energy (DE) fluid.…
Hot, diffuse, relativistic plasmas such as sub-Eddington black hole accretion flows are expected to be collisionless, yet are commonly modeled as a fluid using ideal general relativistic magnetohydrodynamics (GRMHD). Dissipative effects…
We present an implicit numerical method to solve the time-dependent equations of radiation hydrodynamics (RHD) in axial symmetry assuming hydrostatic equilibrium perpendicular to the equatorial plane (1+1D) of a gaseous disk. The equations…
The Fully Constrained Formulation (FCF) of General Relativity is a novel framework introduced as an alternative to the hyperbolic formulations traditionally used in numerical relativity. The FCF equations form a hybrid elliptic-hyperbolic…
In this work, we report the development of a spatially fourth order temporally second order compact scheme for incompressible Navier-Stokes (N-S) equations in time-varying domain. Sen [J. Comput. Phys. 251 (2013) 251-271] put forward an…
Numerical radiation-hydrodynamics (RHD) for non-relativistic flows is a challenging problem because it encompasses processes acting over a very broad range of timescales, and where the relative importance of these processes often varies by…
We consider models of relativistic matter containing sharp interfaces across which the matter model changes. These models will be relevant for neutron stars with crusts, phase transitions, or for viscous boundaries where the length scale is…
We develop a semi-implicit algorithm for time-accurate simulation of the compressible Navier-Stokes equations, with special reference to wall-bounded flows. The method is based on linearization of the partial convective fluxes associated…
Dynamical equations in generalized hydrodynamics (GHD), a hydrodynamic theory for integrable quantum systems at the Euler scale, take a rather simple form, even though an infinite number of conserved charges are taken into account. We show…
Hydrodynamics of the non-relativistic compressible fluid in the curved spacetime is derived using the generalized framework of the stochastic variational method (SVM) for continuum medium. The fluid-stress tensor of the resultant equation…
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…
A detailed derivation of the Lattice Boltzmann (LB) scheme for relativistic fluids recently proposed in Ref. [1], is presented. The method is numerically validated and applied to the case of two quite different relativistic fluid dynamic…
We describe an extension of the cosmological hydrodynamics code ENZO to include the self-consistent transport of ionizing radiation modeled in the flux-limited diffusion approximation. A novel feature of our algorithm is a coupled implicit…
In this paper, we describe a new hydrodynamics code for 1D and 2D astrophysical simulations, BETHE-hydro, that uses time-dependent, arbitrary, unstructured grids. The core of the hydrodynamics algorithm is an arbitrary Lagrangian-Eulerian…
We present a novel numerical implementation of radiative transfer in the cosmological smoothed particle hydrodynamics (SPH) simulation code {\small GADGET}. It is based on a fast, robust and photon-conserving integration scheme where the…
We present a new numerical code designed to solve the Einstein field equations for axisymmetric spacetimes. The long term goal of this project is to construct a code that will be capable of studying many problems of interest in axisymmetry,…
We present a hybrid particle/grid approach for simulating incompressible fluids on collocated velocity grids. We interchangeably use particle and grid representations of transported quantities to balance efficiency and accuracy. A novel…
In an attempt to investigate the structures of ultra-relativistic jets injected into the intracluster medium (ICM) and the associated flow dynamics, such as shocks, velocity shear, and turbulence, we have developed a new special…