Numerical simulations of downward convective overshooting in giants

An attempt at understanding downward overshooting in the convective envelopes of post-main-sequence stars has been made on the basis of three-dimensional large-eddy simulations, using artificially modified OPAL opacity and taking into account radiation and ionization in the equation of state. Two types of star, an intermediate-mass star and a massive star, were considered. To avoid a long thermal relaxation time of the intermediate-mass star, we increased the stellar energy flux artificially while trying to maintain a structure close to the one given by a 1D stellar model. A parametric study of the flux factor was performed. For the massive star, no such process was necessary. Numerical results were analysed when the system reached the statistical steady state. It was shown that the penetration distance in pressure scaleheights is of the order of unity. The scaling relations between penetration distance, input flux and vertical velocity fluctuations studied by Singh et al. were checked. The anisotropy of the turbulent convection and the diffusion models of the third-order moments representing the non-local transport were also investigated. These models are dramatically affected by the velocity fields and no universal constant parameters seem to exist. The limitations of the numerical results were also discussed.     

Hydrodynamical adaptive mesh refinement simulations of turbulent flows – II. Cosmological simulations of galaxy clusters

The development of turbulent gas flows in the intra-cluster medium and in the core of a galaxy cluster is studied by means of adaptive mesh refinement (AMR) cosmological simulations. A series of six runs was performed, employing identical simulation parameters but different criteria for triggering the mesh refinement. In particular, two different AMR strategies were followed, based on the regional variability of control variables of the flow and on the overdensity of subclumps, respectively. We show that both approaches, albeit with different results, are useful to get an improved resolution of the turbulent flow in the ICM. The vorticity is used as a diagnostic for turbulence, showing that the turbulent flow is not highly volume filling but has a large area-covering factor, in agreement with previous theoretical expectations. The measured turbulent velocity in the cluster core is larger than 200 km s⁻¹, and the level of turbulent pressure contribution to the cluster hydrostatic equilibrium is increased by using the improved AMR criteria.     

Hydrodynamical adaptive mesh refinement simulations of turbulent flows – I. Substructure in a wind

The problem of the resolution of turbulent flows in adaptive mesh refinement (AMR) simulations is investigated by means of three-dimensional (3D) hydrodynamical simulations in an idealized setup, representing a moving subcluster during a merger event. AMR simulations performed with the usual refinement criteria based on local gradients of selected variables do not properly resolve the production of turbulence downstream of the cluster. Therefore, we apply novel AMR criteria which are optimised to follow the evolution of a turbulent flow. We demonstrate that these criteria provide a better resolution of the flow past the subcluster, allowing us to follow the onset of the shear instability, the evolution of the turbulent wake and the subsequent back-reaction on the subcluster core morphology. We discuss some implications for the modelling of cluster cold fronts.     

In-shock cooling in numerical simulations

We model a one-dimensional shock-tube using smoothed particle hydrodynamics and investigate the consequences of having finite shock-width in numerical simulations caused by finite resolution of the codes. We investigate the cooling of gas during passage through the shock for three different cooling regimes.
For a theoretical shock temperature of 10⁵ K, the maximum temperature of the gas is much reduced. When the ratio of the cooling time to shock-crossing time was 8, we found a reduction of 25 per cent in the maximum temperature reached by the gas. When the ratio was reduced to 1.2, the maximum temperature reached dropped to 50 per cent of the theoretical value. In both cases the cooling time was reduced by a factor of 2.
At lower temperatures, we are especially interested in the production of molecular hydrogen, and so we follow the ionization level and H₂ abundance across the shock. The effect of in-shock cooling is substantial: the maximum temperature the gas reaches compared with the theoretical temperature is found to vary between 0.15 and 0.81, depending upon the shock strength and mass resolution. The downstream ionization level is reduced from the theoretical level by a factor of between 2.4 and 12.5, and the resulting H₂ abundance by a factor of 1.35 to 2.22.
At temperatures above 10⁵ K, radiative shocks are unstable and will oscillate. We find that the shock jump temperature varies by a factor of 20 because of these oscillations.
We conclude that extreme caution must be exercised when interpreting the results of simulations of galaxy formation.     

Parameter tests within cosmological simulations of galaxy formation

Numerical simulations of galaxy formation require a number of parameters. Some of these are intrinsic to the numerical integration scheme (e.g., the time-step), while others describe the physical model (e.g., the gas metallicity). In this paper we present results of a systematic exploration of the effects of varying a subset of these parameters on simulations of galaxy formation. We use N -body and 'Smoothed Particle Hydrodynamics' techniques to follow the evolution of cold dark matter and gas in a small volume. We compare a fiducial model with 24 different simulations, in which one parameter at a time is varied, focusing on properties such as the relative fraction of hot and cold gas, and the abundance and masses of galaxies. We find that for reasonable choices of numerical values, many parameters have relatively little effect on the galaxies, with the notable exception of the parameters that control the resolution of the simulation and the efficiency with which gas cools.     

Hydrodynamic simulations of merging clusters of galaxies

We present the results of high-resolution AP³M+SPH simulations of merging clusters of galaxies. We find that the compression and shocking of the core gas during a merger can lead to large increases in bolometric X-ray luminosities and emission-weighted temperatures of clusters. Cooling flows are completely disrupted during equal-mass mergers, with the mass deposition rate dropping to zero as the cores of the clusters collide. The large increase in the cooling time of the core gas strongly suggests that cooling flows will not recover from such a merger within a Hubble time. Mergers with subclumps having one eighth of the mass of the main cluster are also found to disrupt a cooling flow if the merger is head-on. However, in this case the entropy injected into the core gas is rapidly radiated away and the cooling flow restarts within a few Gyr of the merger. Mergers in which the subcluster has an impact parameter of 500 kpc do not disrupt the cooling flow, although the mass deposition rate is reduced by ∼30 per cent. Finally, we find that equal mass, off-centre mergers can effectively mix gas in the cores of clusters, while head on mergers lead to very little mixing. Gas stripped from the outer layers of subclumps results in parts of the outer layers of the main cluster being well mixed, although they have little effect on the gas in the core of the cluster. None of the mergers examined here resulted in the intracluster medium being well mixed globally.     

A numerical and experimental study of the time-evolution of a low Mach number jet

We present a comparison between a plasma-generated 'starting jet' experiment and an axisymmetric numerical simulation of the flow. The experimental flow and the numerical simulation give results that agree both qualitatively and quantitatively, showing that the complex vortical structures arising in the flow are surprisingly well reproduced by the numerical model. This result inspires confidence in the accuracy of astrophysical jet numerical simulations. Also, even though the Mach number of our laboratory jet is somewhat low ( M ∼0.5), the dimensionless parameters of this jet are not very far from those expected for Faranoff–Riley class I radio jets.     

Figure rotation of dark haloes in cold dark matter simulations

We investigate the figure rotation of dark matter haloes identified in Λ cold dark matter (CDM) simulations. We find that when strict criteria are used to select suitable haloes for study, five of the 222 haloes identified in our   z = 0  simulation output undergo coherent figure rotation over a  5 h ⁻¹ Gyr  period. We discuss the effects of varying the selection criteria and find that pattern speeds for a much larger fraction of the haloes can be measured when the criteria are relaxed. Pattern speeds measured over a  1 h ⁻¹ Gyr  period follow a lognormal distribution, centred at  Ω_p= 0.2 h rad Gyr⁻¹  with a maximum value of 0.94 h rad Gyr⁻¹. Over a  5 h ⁻¹ Gyr  period, the average pattern speed of a halo is about  0.1 h rad Gyr⁻¹  and the largest pattern speed found is  0.24 h rad Gyr⁻¹  . Less than half of the selected haloes showed alignment between their figure rotation axis and minor axis, the exact fraction being somewhat dependent on how one defines a halo. While the pattern speeds observed are lower than those generally thought capable of causing spiral structure, we note that coherent figure rotation is found over very long periods and argue that further simulations would be required before strong conclusions about spiral structure in all galaxies could be drawn. We find no correlation between halo properties such as total mass and the pattern speed.     

Formation of accretion centers in simulations of colliding uniform density H2 cores

We test here the first stage of a route of modifications to be applied to the public GADGET2 code for dynamically identifying accretion centers during the collision process of two adjacent and identical gas cores. Each colliding core has a uniform density profile and rigid body rotation; its mass and size have been chosen to represent the observed core L1544; for the thermal and rotational energy ratios with respect to the potential energy, we assume the values α = 0.3 and β = 0.1, respectively. These values favor the gravitational collapse of the core. We here study cases of both head‐on and off‐center collisions, in which the pre‐collision velocity increases the initial sound speed of the barotropic gas by up to several times. In a simulation the accretion centers are formed by the highest density particles, so we here report their location and properties in order to realize the collision effects on the collapsing and colliding cores. In one of the models, we observe a roughly spherical distribution of accretion centers located at the front wave of the collision. In a forthcoming publication we will apply the full modified GADGET code to study the collision of turbulent cores. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Hydrodynamic simulations of propagating warps and bending waves in accretion discs

We present the results of a study of propagating warp or bending waves in accretion discs. Three-dimensional hydrodynamic simulations were performed using smoothed particle hydrodynamics (SPH), and the results are compared with calculations based on the linear theory of warped discs.
We examine the response of a gaseous disc to an initially imposed warping disturbance under a variety of physical conditions. We consider primarily the physical regime in which the dimensionless viscosity parameter α < H r , where H r is the disc aspect ratio, so that bending waves are expected to propagate. We also performed calculations for disc models in which α > H r , where the warps are expected to evolve diffusively. Small-amplitude (linear) perturbations are studied in both Keplerian and slightly non-Keplerian discs, and we find that the results of the SPH calculations can be reasonably well fitted by those of the linear theory. The main results of these calculations are: (i) the warp in Keplerian discs when α < H r propagates with little dispersion, and damps at a rate expected from estimates of the code viscosity; (ii) warps evolve diffusively when α > H r ; (iii) the slightly non-Keplerian discs lead to a substantially more dispersive behaviour of the warps, which damp at a similar rate to the Keplerian case, when α < H r .
Initially imposed higher amplitude, non-linear warping disturbances were studied in Keplerian discs. The results indicate that non-linear warps can lead to the formation of shocks, and that the evolution of the warp becomes less wave-like and more diffusive in character.
This work is relevant to the study of the warped accretion discs that may occur around Kerr black holes or in misaligned binary systems, and is mainly concerned with discs in which α < H r . The results indicate that SPH can model the hydrodynamics of warped discs, even when using rather modest numbers of particles.     

Star formation and chemical evolution in smoothed particle hydrodynamics simulations: a statistical approach

In smoothed particle hydrodynamics (SPH) codes with a large number of particles, star formation as well as gas and metal restitution from dying stars can be treated statistically. This approach allows one to include detailed chemical evolution and gas re-ejection with minor computational effort. Here we report on a new statistical algorithm for star formation and chemical evolution, especially conceived for SPH simulations with large numbers of particles, and for parallel SPH codes.
For the sake of illustration, we also present two astrophysical simulations obtained with this algorithm, implemented into the Tree-SPH code by Lia & Carraro .
In the first simulation, we follow the formation of an individual disc-like galaxy, predict the final structure and metallicity evolution, and test resolution effects. In the second simulation we simulate the formation and evolution of a cluster of galaxies, to demonstrate the capabilities of the algorithm in investigating the chemo-dynamical evolution of galaxies and of the intergalactic medium in a cosmological context.     

Simulations of multiphase turbulence in jet cocoons

The interaction of optically emitting clouds with warm X-ray gas and hot, tenuous radio plasma in radio jet cocoons is modelled by 2D compressible hydrodynamic simulations. The initial setup is the Kelvin–Helmholtz instability at a contact surface of density contrast 10⁴. The denser medium contains clouds of higher density. Optically thin radiation is realized via a cooling source term. The cool phase effectively extracts energy from the other gas which is both, radiated away and used for acceleration of the cold phase. This increases the system's cooling rate substantially and leads to a massively amplified cold mass dropout. We show that it is feasible, given small seed clouds of the order of  100 M_⊙  , that all of the optically emitting gas in a radio jet cocoon may be produced by this mechanism on the propagation time-scale of the jet. The mass is generally distributed as   T ^−1/2  with temperature, with a prominent peak at 14 000 K. This peak is likely to be related to the counteracting effects of shock heating and a strong rise in the cooling function. The volume filling factor of cold gas in this peak is of the order of  10⁻⁵–10⁻³  and generally increases during the simulation time.
The simulations tend towards an isotropic scale-free Kolmogorov-type energy spectrum over the simulation time-scale. We find the same Mach-number density relation as Kritsuk & Norman and show that this relation may explain the velocity widths of emission lines associated with high-redshift radio galaxies, if the environmental temperature is lower, or the jet-ambient density ratio is less extreme than in their low-redshift counterparts.