Related papers: Numerical overcooling in shocks
To make relevant predictions about observable emission, hydrodynamical simulation codes must employ schemes that account for radiative losses, but the large dimensionality of accurate radiative transfer schemes is often prohibitive.…
We present analysis of the evolution of dark matter halos in dense environments of groups and clusters in dissipationless cosmological simulations. The premature destruction of halos in such environments, known as the overmerging, reduces…
Smoothed particle hydrodynamics (SPH) employs an artificial viscosity to properly capture hydrodynamical shock waves. In its original formulation, the resulting numerical viscosity is large enough to suppress structure in the velocity field…
This work provides analytical solutions describing the post-shock structure of radiative shocks growing in astrophysics and in laboratory. The equations including a cooling function $\Lambda \propto \rho^{\epsilon} P^{\zeta} x^{\theta}$ are…
N-body simulations that follow only a collisionless dark matter component have failed to produce galaxy halos or substructure within dense environments. We investigate the `over-merging' problem analytically and with numerical simulations,…
Recent simulations of self-gravitating accretion discs, carried out using a three-dimensional Smoothed Particle Hydrodynamics (SPH) code by Meru and Bate, have been interpreted as implying that three-dimensional global discs fragment much…
We test and improve the numerical schemes in our smoothed particle hydrodynamics (SPH) code for cosmological simulations, including the pressure-entropy formulation (PESPH), a time-dependent artificial viscosity, a refined timestep…
We present a suite of cosmological radiation-hydrodynamical simulations of the assembly of galaxies driving the reionization of the intergalactic medium (IGM) at z >~ 6. The simulations account for the hydrodynamical feedback from…
(Abridged) We here continue our effort to model the behaviour of matter when orbiting or accreting onto a generic black hole by developing a new numerical code employing advanced techniques geared solve the equations of in…
We compare the results of numerical simulations of thin and quasi-spherical (thick) accretion flows with existing analytical solutions. We use a Lagrangian code based on the Smooth Particle Hydrodynamics (SPH) scheme and an Eulerian finite…
We present a carefully designed, systematic study of the angular resolution dependence of simulations with the Prometheus-Vertex neutrino-hydrodynamics code. Employing a simplified neutrino heating-cooling scheme in the Prometheus…
We carry out simulations of gravitationally unstable discs using a Smoothed Particle Hydrodynamics (SPH) code and a grid-based hydrodynamics code, FARGO, to understand the previous non-convergent results reported by Meru & Bate (2011a). We…
We present the newly-incorporated gray radiation hydrodynamics capabilities of the FLASH code based on a radiation flux-limiter aware hydrodynamics numerical implementation designed specifically for applications in astrophysical problems.…
We analyze heating and cooling processes in an idealized simulation of a cool-core cluster, where momentum-driven AGN feedback balances radiative cooling in a time-averaged sense. We find that, on average, energy dissipation via shock waves…
We describe an accurate, one-dimensional, spherically symmetric, Lagrangian hydrodynamics/gravity code, designed to study the effects of radiative cooling and photo-ionization on the formation of protogalaxies. The code can treat an…
We investigate the nature of gas accretion onto haloes and galaxies at z=2 using cosmological hydrodynamic simulations run with the moving mesh code AREPO. Implementing a Monte Carlo tracer particle scheme to determine the origin and…
We present a series of 2-d ($r,\phi$) hydrodynamic simulations of marginally self gravitating disks around protostars using an SPH code. We implement simple dynamical heating and we cool each location as a black body, using a photosphere…
We carry out high-resolution adaptive mesh refinement simulations of a cool core cluster, resolving the flow from Mpc scales down to pc scales. We do not (yet) include any AGN heating, focusing instead on cooling in order to understand how…
Numerical effects are known to plague adaptive mesh refinement (AMR) codes when treating massive particles, e.g. representing massive black holes (MBHs). In an evolving background, they can experience strong, spurious perturbations and then…
Numerical shock instability is a complexity which may occur in supersonic simulations. Riemann solver is usually the crucial factor that affects both the computation accuracy and numerical shock stability. In this paper, several classical…