Engineering topics

Bryan, Greg L.

Publications and source records attributed to Bryan, Greg L..

14 records

DOE OSTI2019

Numerical and perturbative computations of the fuzzy dark matter model

We investigate nonlinear structure formation in the fuzzy dark matter (FDM) model using both numerical and perturbative techniques. On the numerical side, we examine the virtues and limitations of a Schrödinger-Poisson solver (wave formulation) versus a fluid dynamics solver (Madelung formulation). On the perturbative side, we carry out a computation of the one-loop mass power spectrum, i.e. up to third order in perturbation theory. We find that (1) in many situations, the fluid dynamics solver is capable of producing the expected interference patterns, but it fails in situations where destructive interference causes the density to vanish—a generic occurrence in the nonlinear regime. (2) The Schrödinger-Poisson solver works well in all test cases, but it is demanding in resolution: suppose one is interested in the mass power spectrum on large scales, it is not sufficient to resolve structure on those same scales; one must resolve the relevant de Broglie scale which is often smaller. The fluid formulation does not suffer from this issue. (3) We compare the one-loop mass power spectrum from perturbation theory against the mass power spectrum from the Schrödinger-Poisson solver, and find good agreement in the mildly nonlinear regime. We contrast fluid perturbation theory with wave perturbation theory; the latter has a more limited range of validity. (4) As an application, we compare the Lyman-alpha forest flux power spectrum obtained from the Schrödinger-Poisson solver versus one from an N-body simulation (the latter is often used as an approximate method to make predictions for FDM). At redshift 5, the two, starting from the same initial condition, agree to better than 10% on observationally relevant scales as long as the FDM mass exceeds 2 × 10 –23 eV. Furthermore, we emphasize that the so called quantum pressure is capable of both enhancing and suppressing fluctuations in the nonlinear regime—which dominates depends on the scale and quantity of interest.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Quantifying Supernovae-driven Multiphase Galactic Outflows

Galactic outflows are observed everywhere in star-forming disk galaxies and are critical for galaxy formation. Supernovae (SNe) play the key role in driving the outflows, but there is no consensus as to how much energy, mass, and metal they can launch out of the disk. We perform 3D, high-resolution hydrodynamic simulations to study SNe-driven outflows from stratified media. Assuming the SN rate scales with gas surface density Σ{sub gas} as in the Kennicutt–Schmidt relation, we find that the mass loading factor, η {sub m}, defined as the mass outflow flux divided by the star formation surface density, decreases with increasing Σ{sub gas} as η{sub m}∝Σ{sub gas}{sup −0.61}. Approximately Σ{sub gas} ≲ 50 M {sub ⊙} pc{sup −2} marks when η {sub m} ≳ 1. About 10%–50% of the energy and 40%–80% of the metals produced by SNe end up in the outflows. The tenuous hot phase (T > 3 × 10{sup 5} K), which fills 60%–80% of the volume at the midplane, carries the majority of the energy and metals in the outflows. We discuss how various physical processes, including the vertical distribution of SNe, photoelectric heating, external gravitational field, and SN rate, affect the loading efficiencies. The relative scale height of gas and SNe is a very important factor in determining the loading efficiencies.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

The Impact of Galactic Winds on the Angular Momentum of Disk Galaxies in the Illustris Simulation

Observed galactic disks have specific angular momenta similar to expectations for typical dark matter halos in ΛCDM. Cosmological hydrodynamical simulations have recently reproduced this similarity in large galaxy samples by including strong galactic winds, but the exact mechanism that achieves this is not yet clear. Here we present an analysis of key aspects contributing to this relation: angular momentum selection and evolution of Lagrangian mass elements as they accrete onto dark matter halos, condense into Milky-Way-scale galaxies, and join the z = 0 stellar phase. We contrast this evolution in the Illustris simulation with that in a simulation without galactic winds, where the z = 0 angular momentum is ≈0.6 dex lower. We find that winds induce differences between these simulations in several ways: increasing angular momentum, preventing angular momentum loss, and causing z = 0 stars to sample the accretion-time angular momentum distribution of baryons in a biased way. In both simulations, gas loses on average ≈0.4 dex between accreting onto halos and first accreting onto central galaxies. In Illustris, this is followed by ≈0.2 dex gains in the “galactic wind fountain” and no further net evolution past the final accretion onto the galaxy. Without feedback, further losses of ≈0.2 dex occur in the gas phase inside the galaxies. An additional ≈0.15 dex difference arises from feedback preferentially selecting higher angular momentum gas at accretion by expelling gas that is poorly aligned. These and additional effects of similar magnitude are discussed, suggesting a complex origin of the similarity between the specific angular momenta of galactic disks and typical halos.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

Supernova feedback and the hot gas filling fraction of the interstellar medium

Supernovae (SNe), the most energetic stellar feedback mechanism, are crucial for regulating the interstellar medium (ISM) and launching galactic winds. We explore how supernova remnants (SNRs) create a multiphase medium by performing three-dimentional hydrodynamical simulations at various SN rates, S, and ISM average densities, n-bar . The evolution of an SNR in a self-consistently generated three-phase ISM is qualitatively different from that in a uniform or a two-phase warm/cold medium. By traveling faster and further in the low-density hot phase, the domain of an SNR increases by >10{sup 2.5}. Varying n-bar and S, we find that a steady state can only be achieved when the hot gas volume fraction f{sub V,hot}≲0.6±0.1. Above that level, overlapping SNRs render connecting topology of the hot gas, and the ISM is subjected to thermal runaway. Photoelectric heating (PEH) has a surprisingly strong impact on f{sub V,hot}. For n-bar ≳3 cm{sup −3}, a reasonable PEH rate is able to suppress the thermal runaway. Overall, we determine the critical SN rate for the onset of thermal runaway to be S{sub crit}=200( n-bar /1 cm{sup −3}){sup k}(E{sub SN}/10{sup 51} erg){sup −1} kpc{sup −3} Myr{sup −1}, where k = (1.2, 2.7) for n-bar ⩽1 and >1 cm{sup −3}, respectively. We present a fitting formula of the ISM pressure P( n-bar ,S), which can be used as an effective equation of state in cosmological simulations. Despite the five orders of magnitude span of ( n-bar ,S), the average Mach number varies little: M ≈ 0.5 ± 0.2, 1.2 ± 0.3, and 2.3 ± 0.9 for the hot, warm, and cold phases, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

Cooling, AGN feedback, and star formation in simulated cool-core galaxy clusters

Numerical simulations of active galactic nuclei (AGNs) feedback in cool-core galaxy clusters have successfully avoided classical cooling flows, but often produce too much cold gas. We perform adaptive mesh simulations that include momentum-driven AGN feedback, self-gravity, star formation, and stellar feedback, focusing on the interplay between cooling, AGN heating, and star formation in an isolated cool-core cluster. Cold clumps triggered by AGN jets and turbulence form filamentary structures tens of kpc long. This cold gas feeds both star formation and the supermassive black hole (SMBH), triggering an AGN outburst that increases the entropy of the intracluster medium (ICM) and reduces its cooling rate. Within 1–2 Gyr, star formation completely consumes the cold gas, leading to a brief shutoff of the AGN. The ICM quickly cools and redevelops multiphase gas, followed by another cycle of star formation/AGN outburst. Within 6.5 Gyr, we observe three such cycles. There is good agreement between our simulated cluster and the observations of cool-core clusters. ICM cooling is dynamically balanced by AGN heating, and a cool-core appearance is preserved. The minimum cooling time to free-fall time ratio typically varies between a few and ≳20. The star formation rate (SFR) covers a wide range, from 0 to a few hundred M{sub ⊙} yr{sup −1}, with an average of ∼40 M{sub ⊙} yr{sup −1}. The instantaneous SMBH accretion rate shows large variations on short timescales, but the average value correlates well with the SFR. Simulations without stellar feedback or self-gravity produce qualitatively similar results, but a lower SMBH feedback efficiency (0.1% compared to 1%) results in too many stars.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

KINETIC ENERGY FROM SUPERNOVA FEEDBACK IN HIGH-RESOLUTION GALAXY SIMULATIONS

We describe a new method for adding a prescribed amount of kinetic energy to simulated gas modeled on a cartesian grid by directly altering grid cells’ mass and velocity in a distributed fashion. The method is explored in the context of supernova (SN) feedback in high-resolution (∼10 pc) hydrodynamic simulations of galaxy formation. Resolution dependence is a primary consideration in our application of the method, and simulations of isolated explosions (performed at different resolutions) motivate a resolution-dependent scaling for the injected fraction of kinetic energy that we apply in cosmological simulations of a 10{sup 9} M{sub ⊙} dwarf halo. We find that in high-density media (≳50 cm{sup −3}) with coarse resolution (≳4 pc per cell), results are sensitive to the initial kinetic energy fraction due to early and rapid cooling. In our galaxy simulations, the deposition of small amounts of SN energy in kinetic form (as little as 1%) has a dramatic impact on the evolution of the system, resulting in an order-of-magnitude suppression of stellar mass. The overall behavior of the galaxy in the two highest resolution simulations we perform appears to converge. We discuss the resulting distribution of stellar metallicities, an observable sensitive to galactic wind properties, and find that while the new method demonstrates increased agreement with observed systems, significant discrepancies remain, likely due to simplistic assumptions that neglect contributions from SNe Ia and stellar winds.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

COSMOLOGICAL SIMULATIONS OF GALAXY FORMATION WITH COSMIC RAYS

We investigate the dynamical impact of cosmic rays (CR) in cosmological simulations of galaxy formation using adaptive-mesh refinement simulations of a 10{sup 12} M {sub ☉} halo. In agreement with previous work, a run with only our standard thermal energy feedback model results in a massive spheroid and unrealistically peaked rotation curves. However, the addition of a simple two-fluid model for CRs drastically changes the morphology of the forming disk. We include an isotropic diffusive term and a source term tied to star formation due to (unresolved) supernova-driven shocks. Over a wide range of diffusion coefficients, the CRs generate thin, extended disks with a significantly more realistic (although still not flat) rotation curve. We find that the diffusion of CRs is key to this process, as they escape dense star-forming clumps and drive outflows within the more diffuse interstellar medium.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

AN ADAPTIVE PARTICLE-MESH GRAVITY SOLVER FOR ENZO

We describe and implement an adaptive particle-mesh algorithm to solve the Poisson equation for grid-based hydrodynamics codes with nested grids. The algorithm is implemented and extensively tested within the astrophysical code Enzo against the multigrid solver available by default. We find that while both algorithms show similar accuracy for smooth mass distributions, the adaptive particle-mesh algorithm is more accurate for the case of point masses, and is generally less noisy. We also demonstrate that the two-body problem can be solved accurately in a configuration with nested grids. In addition, we discuss the effect of subcycling, and demonstrate that evolving all the levels with the same timestep yields even greater precision.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Modeling active galactic nucleus feedback in cool-core clusters: The formation of cold clumps

We perform high-resolution (15-30 pc) adaptive mesh simulations to study the impact of momentum-driven active galactic nucleus (AGN) feedback in cool-core clusters, focusing in this paper on the formation of cold clumps. The feedback is jet-driven with an energy determined by the amount of cold gas within 500 pc of the super-massive black hole. When the intracluster medium in the core of the cluster becomes marginally stable to radiative cooling, with the thermal instability to the free-fall timescale ratio t{sub TI}/t{sub ff} < 3-10, cold clumps of gas start to form along the propagation direction of the AGN jets. By tracing the particles in the simulations, we find that these cold clumps originate from low entropy (but still hot) gas that is accelerated by the jet to outward radial velocities of a few hundred km s{sup –1}. This gas is out of hydrostatic equilibrium and so can cool. The clumps then grow larger as they decelerate and fall toward the center of the cluster, eventually being accreted onto the super-massive black hole. The general morphology, spatial distribution, and estimated Hα morphology of the clumps are in reasonable agreement with observations, although we do not fully replicate the filamentary morphology of the clumps seen in the observations, probably due to missing physics.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Modeling active galactic nucleus feedback in cool-core clusters: The balance between heating and cooling

We study the long-term evolution of an idealized cool-core galaxy cluster under the influence of momentum-driven active galactic nucleus (AGN) feedback using three-dimensional high-resolution (60 pc) adaptive mesh refinement simulations. The feedback is modeled with a pair of precessing jets whose power is calculated based on the accretion rate of the cold gas surrounding the supermassive black hole (SMBH). The intracluster medium first cools into clumps along the propagation direction of the jets. As the jet power increases, gas condensation occurs isotropically, forming spatially extended structures that resemble the observed Hα filaments in Perseus and many other cool-core clusters. Jet heating elevates the gas entropy, halting clump formation. The cold gas that is not accreted onto the SMBH settles into a rotating disk of ∼10{sup 11} M {sub ☉}. The hot gas cools directly onto the disk while the SMBH accretes from its innermost region, powering the AGN that maintains a thermally balanced state for a few Gyr. The mass cooling rate averaged over 7 Gyr is ∼30 M {sub ☉} yr{sup –1}, an order of magnitude lower than the classic cooling flow value. Medium resolution simulations produce similar results, while in low resolution runs, the cluster experiences cycles of gas condensation and AGN outbursts. Owing to its self-regulating mechanism, AGN feedback can successfully balance cooling with a wide range of model parameters. Our model also produces cold structures in early stages that are in good agreement with the observations. However, the long-lived massive cold disk is unrealistic, suggesting that additional physical processes are still needed.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

ENZO: AN ADAPTIVE MESH REFINEMENT CODE FOR ASTROPHYSICS

This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in one, two, and three dimensions, and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gravity of fluids and particles), primordial gas chemistry, optically thin radiative cooling of primordial and metal-enriched plasmas (as well as some optically-thick cooling models), radiation transport, cosmological expansion, and models for star formation and feedback in a cosmological context. In addition to explaining the algorithms implemented, we present solutions for a wide range of test problems, demonstrate the code's parallel performance, and discuss the Enzo collaboration's code development methodology.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1994

The evolution of X-ray clusters in a cold plus hot dark matter universe

We present the first self-consistently computed results on the evolution of X-ray properties of galaxy clusters in a cold + hot dark matter (CHDM) model. We have performed a hydrodynamic plus N-body simulation for the COBE-compatible CHDM model with standard mass components: Omega(sub hot) = 0.3, Omega (sub cold) = 0.6 and Omega(sub baryon) = 0.1 (h = 0.5). In contrast with the CDM model, which fails to reproduce the observed temperature distribution function dN/dT (Bryan et al. 1994b), the CHDM model fits the observational dN/dT quite well. Our results on X-ray luminosity are less firm but even more intriguing. We find that the resulting X-ray luminosity functions at redshifts z = 0.0, 0.2, 0.4, 0.7 are well fit by observations, where they overlap. The fact that both temperatures and luminosities provide a reasonable fit to the available observational data indicates that, unless we are missing some essential physics, there is neither room nor need for a large fraction of gas in rich clusters: 10% (or less) in baryons is sufficient to explain their X-ray properties. We also see a tight correlation between X-ray luminosity and gas temperature.

Bryan, Greg L.↗

NASA NTRS1994

A comparison of cosmological hydrodynamic codes

We present a detailed comparison of the simulation results of various hydrodynamic codes. Starting with identical initial conditions based on the cold dark matter scenario for the growth of structure, with parameters h = 0.5 Omega = Omega(sub b) = 1, and sigma(sub 8) = 1, we integrate from redshift z = 20 to z = O to determine the physical state within a representative volume of size L(exp 3) where L = 64 h(exp -1) Mpc. Five indenpendent codes are compared: three of them Eulerian mesh-based and two variants of the smooth particle hydrodynamics 'SPH' Lagrangian approach. The Eulerian codes were run at N(exp 3) = (32(exp 3), 64(exp 3), 128(exp 3), and 256(exp 3)) cells, the SPH codes at N(exp 3) = 32(exp 3) and 64(exp 3) particles. Results were then rebinned to a 16(exp 3) grid with the exception that the rebinned data should converge, by all techniques, to a common and correct result as N approaches infinity. We find that global averages of various physical quantities do, as expected, tend to converge in the rebinned model, but that uncertainites in even primitive quantities such as (T), (rho(exp 2))(exp 1/2) persists at the 3%-17% level achieve comparable and satisfactory accuracy for comparable computer time in their treatment of the high-density, high-temeprature regions as measured in the rebinned data; the variance among the five codes (at highest resolution) for the mean temperature (as weighted by rho(exp 2) is only 4.5%. Examined at high resolution we suspect that the density resolution is better in the SPH codes and the thermal accuracy in low-density regions better in the Eulerian codes. In the low-density, low-temperature regions the SPH codes have poor accuracy due to statiscal effects, and the Jameson code gives the temperatures which are too high, due to overuse of artificial viscosity in these high Mach number regions. Overall the comparison allows us to better estimate errors; it points to ways of improving this current generation ofhydrodynamic codes and of suiting their use to problems which exploit their best individual features.

Kang, Hyesung↗

NASA NTRS1994

X-ray clusters from a high-resolution hydrodynamic PPM simulation of the cold dark matter universe

A new three-dimensional hydrodynamic code based on the piecewise parabolic method (PPM) is utilized to compute the distribution of hot gas in the standard Cosmic Background Explorer (COBE)-normalized cold dark matter (CDM) universe. Utilizing periodic boundary conditions, a box with size 85 h(exp-1) Mpc, having cell size 0.31 h(exp-1) Mpc, is followed in a simulation with 270(exp 3)=10(exp 7.3) cells. Adopting standard parameters determined from COBE and light-element nucleosynthesis, Sigma(sub 8)=1.05, Omega(sub b)=0.06, we find the X-ray-emitting clusters, compute the luminosity function at several wavelengths, the temperature distribution, and estimated sizes, as well as the evolution of these quantities with redshift. The results, which are compared with those obtained in the preceding paper (Kang et al. 1994a), may be used in conjuction with ROSAT and other observational data sets. Overall, the results of the two computations are qualitatively very similar with regard to the trends of cluster properties, i.e., how the number density, radius, and temeprature depend on luminosity and redshift. The total luminosity from clusters is approximately a factor of 2 higher using the PPM code (as compared to the 'total variation diminishing' (TVD) code used in the previous paper) with the number of bright clusters higher by a similar factor. The primary conclusions of the prior paper, with regard to the power spectrum of the primeval density perturbations, are strengthened: the standard CDM model, normalized to the COBE microwave detection, predicts too many bright X-ray emitting clusters, by a factor probably in excess of 5. The comparison between observations and theoretical predictions for the evolution of cluster properties, luminosity functions, and size and temperature distributions should provide an important discriminator among competing scenarios for the development of structure in the universe.

Bryan, Greg L.↗