Engineering topics

Norman, Michael L.

Publications and source records attributed to Norman, Michael L..

18 records

NASA NTRS1997

Parallel Multiscale Algorithms for Astrophysical Fluid Dynamics Simulations

Our goal is to develop software libraries and applications for astrophysical fluid dynamics simulations in multidimensions that will enable us to resolve the large spatial and temporal variations that inevitably arise due to gravity, fronts and microphysical phenomena. The software must run efficiently on parallel computers and be general enough to allow the incorporation of a wide variety of physics. Cosmological structure formation with realistic gas physics is the primary application driver in this work. Accurate simulations of e.g. galaxy formation require a spatial dynamic range (i.e., ratio of system scale to smallest resolved feature) of 104 or more in three dimensions in arbitrary topologies. We take this as our technical requirement. We have achieved, and in fact, surpassed these goals.

Norman, Michael L.↗

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

Shock interactions with magnetized interstellar clouds. 1: Steady shocks hitting nonradiative clouds

We study the interaction of a steady, planar shock with a nonradiative, spherical, interstellar cloud threaded by a uniform magnetic field. For strong shocks, the sonic Mach number scales out, so two parameters determine the evolution: the ratio of cloud to intercloud density, and the Alfven Mach number. We focus on the case with initial field parallel to the shock velocity, though we also present one model with field perpendicular to the velocity. Even with 100 zones per cloud radius, we find that the magnetic field structure converges only at early times. However, we can draw three conclusions from our work. First, our results suggest that the inclusion of a field in equipartition with the preshock medium can prevent the complete destruction of the cloud found in the field-free case recently considered by Klein, McKee, & Colella. Second, the interaction of the shock with the cloud can amplify the magnetic field in some regions up to equipartition with the post-shock thermal pressure. In the parallel-field case, the shock preferentially amplifies the parallel component of the field, creating a 'flux rope,' a linear structure of concentrated magnetic field. The flux rope dominates the volume of amplified field, so that laminar, rather than turbulent, amplification is dominant in this case. Third, the presence of the cloud enhances the production of X-ray and synchrotron emission. The X-ray emission peaks early, during the initial passage of the shock over the cloud, while the synchrotron emission peaks later, when the flow sweeps magnetic field onto the axis between the cloud and the main shock.

Low, Mordecai-Mark Mac↗

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

The evolution of cosmic-ray-mediated magnetohydrodynamic shocks: A two-fluid approach

We study the shock structure and acceleration efficiency of cosmic-ray mediated Magnetohydrodynamic (MHD) shocks both analytically and numerically by using a two-fluid model. Our model includes the dynamical effect of magnetic fields and cosmic rays on a background thermal fluid. The steady state solution is derived by following the technique of Drury & Voelk (1981) and compared to numerical results. We explore the time evolution of plane-perpendicular, piston-driven shocks. From the results of analytical and numerical studies, we conclude that the mean magnetic field plays an important role in the structure and acceleration efficiency of cosmic-ray mediated MHD shocks. The acceleration of cosmic-ray particles becomes less efficient in the presence of strong magnetic pressure since the field makes the shock less compressive. This feature is more prominent at low Mach numbers than at high Mach numbers.

Jun, Byung-Il↗

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.↗

NASA NTRS1993

Hydrodynamical simulations of star-gas interactions in the interstellar medium with an external gravitational potential

We have calculated 2D models for the global structure of the interstellar medium in a disk galaxy where stars and gas are coupled through star formation, mass loss, stellar heating of the gas, and optically thin radiative cooling. In one set of models, the two dimensions are in the plane of a disk, and for one case, the total size of the region is 1 kpc square. Using a more accurate numerical approach, we confirm the results of Chiang and Bregman (1988) that H I forms into filaments or sheets separated by hotter, more diffuse gas. Although the filaments often have denser subregions, individual clouds of neutral gas are rare. We calculated the vertical distribution of the gas in another set of models where one dimension was in the plane and the other was perpendicular to it. A filament and hot gas network still occurs, although dense neutral complexes are more common and occur near the midplane. The hot gas often cools as it rises, leading to the highest velocity material having both positive and negative velocities. We find other features that are qualitatively similar to observations; for instance, H I 'worms' are produced.

Rosen, Alexander↗

NASA NTRS1993

Nonlinear growth of dynamical overstabilities in blast waves

The numerical gasdynamics code ZEUS-2D is used to directly model the dynamical overstabilities in blast waves. The linear analysis is confirmed by perturbing a blast wave with a low-amplitude eigenfunction of the overstability. The amplitude of the perturbations is increased in order to determine the nonlinear behavior of the overstabilities. The overstability is found to saturate due to weak transverse shocks in the shell. Transverse velocities in the dense shell reach the postshock sound speed, and high-density regions with sizes of the order of the shell thickness form. Transverse oscillations continue even after saturation. This confirms and explains the damping of the overstability experimentally discovered by Grun et al. (1991).

Mac Low, Mordecai-Mark↗

NASA NTRS1993

The hydrodynamics of galaxy formation on Kiloparsec scales

Two dimensional numerical simulations of Zeldovich pancake fragmentation in a dark matter dominated universe were carried out to study the hydrodynamical and gravitational effects on the formation of structures such as protogalaxies. Preliminary results were given in Yuan, Centrella and, Norman (1991). Here we report a more exhaustive study to determine the sensitivity of protogalaxies to input parameters. The numerical code we used for the simulations combines the hydrodynamical code ZEUS-2D (Stone and Norman, 1992) which was modified to include the expansion of the universe and radiative cooling of the gas with a particle-mesh code which follows the motion of dark matter particles. The resulting hybrid code is able to handle highly nonuniform grids which we utilized to obtain a high resolution (much greater than 1 kpc) in the dense region of the pancake.

Norman, Michael L.↗

NASA NTRS1993

Numerical modeling of the interstellar medium in galactic disks

We have been developing detailed hydrodynamic models of the global interstellar medium in the hope of understanding the mass and volume occupied by various phases, as well as their structure and kinematics. In our model, the gas is modeled by one fluid while representative Pop 1 stars are modeled by a second fluid. The two fluids are coupled in that the gas forms into stars at a rate given by a Schmidt law while stellar mass loss returns matter into the gas phase (on a time scale of 100 Myr). Also, the stars heat the gas through stellar winds and the gas cools through optically thin radiation. The time behavior of these two fluids is studied in two spatial dimensions with the Eulerian finite difference numerical hydrodynamic code Zen. The two spatial dimensions are along the plane of a disk (x, total length of 2 kpc) and perpendicular to the disk (z, total height of +/- 15 kpc) and a galactic gravitational field in the z direction, typical of that at the solar circle, is imposed upon the simulation; self-gravity and rotation are absent. For the boundary conditions, outflow is permitted at the top and bottom of the grid (z = +/- 15 kpc) while periodic boundary conditions are imposed upon left and right sides of the grid. As initial conditions, we assumed a gaseous distribution like that seen for the H1 by earlier researchers, although the results are insensitive to the initial conditions. We have run simulations in which the heating due to stars, parameterized as a stellar wind velocity, a, is varied from low (a = 150 km/s), to intermediate (a = 300 km/s), to high (a = 600 km/s). Since the intermediate case is roughly equivalent to the Galactic energy injection rate from supernovae, this summary will concentrate on results from this simulation.

Rosen, A.↗

NASA NTRS1992

A test suite for magnetohydrodynamical simulations

A collection is presented of MHD problems that will provide researchers who have interests in modeling MHD phenomena with a battery of tests (a 'test suite') for calibrating their numerical algorithms. Use of these tests will provide a common reference for comparison of different numerical MHD algorithms. The test suite includes both 1D and 2D problems. Taken together, these problems test the abilities of a numerical method to propagate accurately all of the MHD wave families in moving and stationary media, to capture MHD shocks and contact discontinuities, and to model accurately Lorentz force terms in multiple dimensions. An example solution of each test problem is presented. Diagnostic convergence-testing procedures that provide a quantitative evaluation of MHD algorithms are demonstrated.

Stone, James M.↗

NASA NTRS1992

Instabilities in astrophysical jets. II - Numerical simulations of slab jets

Numerical simulations of an unstable supersonic slab-symmetric jet are described. The instabilities within the jet are characterized by growing internal body waves and their coupled surface waves that are also predicted in linear perturbation theory. The characteristic theory of fluid dynamics is used to help interpret the wave morphologies. It is demonstrated that these waves can be excited by imposing an arbitrary disturbance. From the numerical simulations, it is found that the sound waves propagating against the flow slow down as they propagate outward, and they grow in amplitude. These waves eventually disrupt the jet at a certain length. This disruption length is related to the jet Mach number and the perturbation intensity. Thus, the Mach number of a jet observed with a radio telescope can be estimated by measuring the disruption length and estimating the perturbation intensity. The jet Mach numbers in radio tailed sources determined in this way agree quite well with estimates from ram pressure bending arguments. The wiggles and flares observed in many extragalactic jets, especially in tailed radio sources, appear to be intimately related to instabilities and the jet disruption process.

Zhao, Jun-Hui↗

NASA NTRS1991

Numerical models of extragalactic radio sources

When supercomputer-implemented numerical simulations analyzing the nonlinear physics inherent in the hydrodynamic and MHD equations are applied to extragalactic radio sources, many of the complex structures observed on telescopic images are reproduced. Attention is presently given to recently obtained results from 2D and 3D numerical simulations of the formation and evolution of extended radio morphologies; these numerical models allow the exploration of such physical phenomena as the role of magnetic fields in the dynamics and emissivity of extended radio galaxies, intermittent outflow from the cores of active galaxies, fluid-jet instabilities, and the bending of collimated outflows by motion through the intergalactic medium.

Burns, Jack O.↗

NASA NTRS1991

Computational astrophysical fluid dynamics

The field of astrophysical fluid dynamics (AFD) is described as an emerging discipline which derives historically from both the theory of stellar evolution and space plasma physics. The fundamental physical assumption behind AFD is that fluid equations of motion accurately describe the evolution of plasmas on scales that are large in comparison with particle interaction length scales. Particular attention is given to purely fluid models of large-scale astrophysical plasmas. The role of computer simulation in AFD research is also highlighted and a suite of general-purpose application codes for AFD research is discussed. The codes are called ZEUS-2D and ZEUS-3D and solve the equations of AFD in two and three dimensions, respectively, in several coordinate geometries for general initial and boundary conditions. The topics of bipolar outflows from protostars, galactic superbubbles and supershells, and extragalactic radio sources are addressed.

Norman, Michael L.↗

NASA NTRS1990

Magnetic collimation of protostellar winds into bipolar outflows

Researchers describe self-consistent 2-D magnetohydrodynamic (MHD) simulations of the collimation of an isotropic protostellar wind into bipolar outflows by magnetic stresses in the ambient medium. A variety of ambient field strengths, wind luminosities, and density profiles were studied. Collimation occurs when the energy of the magnetic field swept up by the expanding bubble approaches the bubble thermal energy. Measured axial and radial expansion rates are in good agreement with the analytical predictions of Konigl (1982).

Norman, Michael L.↗

NASA NTRS1990

Expanding hydrodynamical jets crossing a galactic halo/intergalactic medium interface

Parameters within ranges that are plausible for radio sources are presently used to perform two-dimensional hydrodynamical calculations of axisymmetric, initially conical, jets whose initial propagation is through isothermal galactic halos with power-law density distributions; these emerge across a pressure-matched interface into a hotter, but less dense medium whose parameters are typical of an intracluster or intergalactic gas. Upon crossing this interface, the jets accelerate and focused toward cylindrical shapes having long, narrow cocoons.

Wiita, Paul J.↗

NASA NTRS1989

Low-temperature transonic cooling flows in galaxy clusters

Calculations are presented which demonstrate that cooling flow models with large sonic radii may be consistent with observed cluster gas properties. It is found that plausible cluster parameters and cooling flow mass accretion rates can produce sonic radii of 10-20 kpc for sonic point temperatures of 1-3 x 10 to the 6th K. The numerical calculations match these cooling flows to hydrostatic atmosphere solutions for the cluster gas beyond the cooling flow region. The cooling flows produce no appreciable 'holes' in the surface brightness toward the cluster center, and the model can be made to match the observed X-ray surface brightness of three clusters in which cooling flows had been believed to be absent. It is suggested that clusters with low velocity dispersion may be the natural location for such 'cool' cooling flows, and fits of these models to the X-ray surface brightness profiles for three clusters are presented.

Sulkanen, Martin E.↗

NASA NTRS1989

Superbubble blowout dynamics

Multiple supernovae and stellar winds from OB associations carve large holes filled with hot gas in the galactic disk. These superbubbles sweep up H I into cold, thin, dense shells and eventually grow large enough to blow completely out of the galactic H I disk. When superbubbles blow out of the disk, they vent hot gas and supernova energy into the galactic corona. In this paper ZEUS, a two-dimensional hydrodynamics code, is used to model the blowout of a superbubble from exponential and Gaussian models for the vertical density stratification. The results are compared to those from the Kompaneets (thin-shell) approximation. It is found that this approximation works very well, and that most of the mass of the shell remains in the plane, with 5 percent of it accelerating upward. The venting of the hot gas and the stability of the shell depends strongly on the model of the density distribution. It is suggested that the low galactic halo actually consists of a froth of merged superbubbles.

Mac Low, Mordecai-Mark↗