اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدTesting Hydrodynamics Schemes in Galaxy Disc Simulations
43
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
We examine how three fundamentally different numerical hydrodynamics codes follow the evolution of an isothermal galactic disc with an external spiral potential. We compare an adaptive mesh refinement code (RAMSES), a smoothed particle hydrodynamics code (sphNG), and a volume-discretised meshless code (GIZMO). Using standard refinement criteria, we find that RAMSES produces a disc that is less vertically concentrated and does not reach such high densities as the sphNG or GIZMO runs. The gas surface density in the spiral arms increases at a lower rate for the RAMSES simulations compared to the other codes. There is also a greater degree of substructure in the sphNG and GIZMO runs and secondary spiral arms are more pronounced. By resolving the Jeans length with a greater number of grid cells we achieve more similar results to the Lagrangian codes used in this study. Other alterations to the refinement scheme (adding extra levels of refinement and refining based on local density gradients) are less successful in reducing the disparity between RAMSES and sphNG/GIZMO. Although more similar, sphNG displays different density distributions and vertical mass profiles to all modes of GIZMO (including the smoothed particle hydrodynamics version). This suggests differences also arise which are not intrinsic to the particular method but rather due to its implementation. The discrepancies between codes (in particular, the densities reached in the spiral arms) could potentially result in differences in the locations and timescales for gravitational collapse, and therefore impact star formation activity in more complex galaxy disc simulations.
قيم البحث
اقرأ أيضاً
There is strong evidence that the diffuse ionized gas (DIG) in disc galaxies is photoionized by radiation from UV luminous O and B stars in the galactic disc, both from observations and detailed numerical models. However, it is still not clear what m
echanism is responsible for providing the necessary pressure support for a diffuse gas layer at kpc-scale above the disc. In this work we investigate if the pressure increase caused by photoionization can provide this support. We run self-consistent radiation hydrodynamics models of a gaseous disc in an external potential. We find that photoionization feedback can drive low levels of turbulence in the dense galactic disc, and that it provides pressure support for an extended diffuse gas layer. Our results show that there is a natural fine-tuning between the total ionizing radiation budget of the sources in the galaxy and the amount of gas in the different ionization phases of the ISM, and provide the first fully consistent radiation hydrodynamics model of the DIG.
Growth of the black holes (BHs) from the seeds to supermassive BHs (SMBHs, $sim!10^9,M_odot$) is not understood, but the mass accretion must have played an important role. We performed two-dimensional radiation hydrodynamics simulations of line-drive
n disc winds considering the metallicity dependence in a wide range of the BH mass, and investigated the reduction of the mass accretion rate due to the wind mass loss. Our results show that denser and faster disc winds appear at higher metallicities and larger BH masses. The accretion rate is suppressed to $sim! 0.4$--$0.6$ times the mass supply rate to the disc for the BH mass of $M_{rm BH}gtrsim 10^5,M_{odot}$ in high-metallicity environments of $Zgtrsim Z_odot$, while the wind mass loss is negligible when the metallicity is sub-solar ($sim 0.1Z_odot$). By developing a semi-analytical model, we found that the metallicity dependence of the line force and the BH mass dependence of the surface area of the wind launch region are the cause of the metallicity dependence ($propto! Z^{2/3}$) and BH mass dependencies ($propto! M_{rm BH}^{4/3}$ for $M_{rm BH}leq 10^6,M_odot$ and $propto! M_{rm BH}$ for $M_{rm BH}geq 10^6,M_odot$) of the mass-loss rate. Our model suggests that the growth of BHs by the gas accretion effectively slows down in the regime $gtrsim 10^{5}M_odot$ in metal-enriched environments $gtrsim Z_odot$. This means that the line-driven disc winds may have an impact on late evolution of SMBHs.
We present results from a subset of simulations from the Evolution and Assembly of GaLaxies and their Environments (EAGLE) suite in which the formulation of the hydrodynamics scheme is varied. We compare simulations that use the same subgrid models w
ithout re-calibration of the parameters but employing the standard GADGET flavour of smoothed particle hydrodynamics (SPH) instead of the more recent state-of-the-art ANARCHY formulation of SPH that was used in the fiducial EAGLE runs. We find that the properties of most galaxies, including their masses and sizes, are not significantly affected by the details of the hydrodynamics solver. However, the star formation rates of the most massive objects are affected by the lack of phase mixing due to spurious surface tension in the simulation using standard SPH. This affects the efficiency with which AGN activity can quench star formation in these galaxies and it also leads to differences in the intragroup medium that affect the X-ray emission from these objects. The differences that can be attributed to the hydrodynamics solver are, however, likely to be less important at lower resolution. We also find that the use of a time step limiter is important for achieving the feedback efficiency required to match observations of the low-mass end of the galaxy stellar mass function.
Smoothed Particle Hydrodynamics (SPH) is a ubiquitous numerical method for solving the fluid equations, and is prized for its conservation properties, natural adaptivity, and simplicity. We introduce the Sphenix SPH scheme, which was designed with th
ree key goals in mind: to work well with sub-grid physics modules that inject energy, be highly computationally efficient (both in terms of compute and memory), and to be Lagrangian. Sphenix uses a Density-Energy equation of motion, along with variable artificial viscosity and conduction, including limiters designed to work with common sub-grid models of galaxy formation. In particular, we present and test a novel limiter that prevents conduction across shocks, preventing spurious radiative losses in feedback events. Sphenix is shown to solve many difficult test problems for traditional SPH, including fluid mixing and vorticity conservation, and it is shown to produce convergent behaviour in all tests where this is appropriate. Crucially, we use the same parameters within Sphenix for the various switches throughout, to demonstrate the performance of the scheme as it would be used in production simulations.
We have performed Smoothed Particle Hydrodynamics (SPH) simulations to study the time evolution of one and two protoplanets embedded in a protoplanetary accretion disc. We investigate accretion and migration rates of a single protoplanet depending on
several parameters of the protoplanetary disc, mainly viscosity and scale height. Additionally, we consider the influence of a second protoplanet in a long time simulation and examine the migration of the two planets in the disc, especially the growth of eccentricity and chaotic behaviour. One aim of this work is to establish the feasibility of SPH for such calculations considering that usually only grid-based methods are adopted. To resolve shocks and to prevent particle penetration, we introduce a new approach for an artificial viscosity, which consists of an additional artificial bulk viscosity term in the SPH-representation of the Navier-Stokes equation. This allows for an accurate treatment of the physical kinematic viscosity to describe the shear, without the use of an artificial shear viscosity.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
