Search NASA⌕ Search

SEARCH · Search NASA

Results for “hydrodynamic simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

The shadow wind in high-mass X-ray binaries

We present hydrodynamic simulations of the most X-ray luminous high-mass X-ray binary systems, such as Cen X-3 and SMC X-1. These models illustrate the presence of both a normal radiatively driven wind confined to the X-ray shadow of the primary star -- a shadow wind -- and a thermally driven wind excited by the X-ray heating of the primary's stellar surface -- an X-ray-excited wind. The X-ray flux in these systems is sufficiently intense that any circumstellar gas exposed to the X-ray source will be highly photoionized. These extreme ionization conditions prevent the formation of a normal radiatively driven wind from the irradiated surface of the primary, but such a wind can still form on the shadowed side of the primary. Orbital rotation can then bring this shadow wind into the line of sight toward the X-ray source, enhancing the column density of the wind seen near eclipse egress. Furthermore, such a high X-ray flux can also excite a thermal wind from the irradiated surface of the primary. Again, orbital rotation tends to deflect the wind, this time leading to an enhanced column density near eclipse ingress.

Blondin, John M.↗

Extended hot-gas halos around starburst galaxies

A reanalysis of Einstein Imaging Proportional Counter (IPC) data and new observations from the Ginga Large Area (Proportional) counters (LAC) indicate the presence of extended X-ray emission (10-50 kpc) around the starburst galaxy M82. Here, we discuss our model of this emission, which was obtained by performing numerical hydrodynamic simulations of the starburst event to much later times and larger scales than were previously considered. For our models, we adopted a supernova rate of 0.1/yr, and an extended low-density static halo that is bound to the galaxy. There are three stages to the evolution of the wind-blown bubble and the propagation of the shock front: the bubble expands in an almost uniform density disk gas, with a deceleration of the shock front (t is less than or approximately 3.6 Myr); breakout from the disk and the upward acceleration of the shock front (3.6 Myr is less than or approximately t is less than or approximately 18 Myr); propagation into the halo, leading to a more spherical system and shock deceleration (18 Myr is less than or approximately t). For a halo density of 10(exp -3)/cu cm, the outflow reaches a distance of 40-50 kpc from the center of the starburst galaxy in 50 Myr. We calculated the time evolution of the X-ray luminosity and found that the extended starburst emits 3 x 10(exp 39) erg/s to 10(exp 40) in the Ginga LAC band and approximately 10(exp 41) erg/s in the Einstein band. The degree of the ionization equilibrium in the outflow and its effect on the iron K alpha line emission are discussed.

Tomisaka, Kohji↗

Using X-rays to determine which compact groups are illusory

If the large-scale galaxy distribution is filamentary, as suggested by some observations and recent hydrodynamical simulations, then lengthwise views of filaments will apparently produce compact groups (CGs) that are in reality stretched out along the line of sight. This possibility has been advocated recently by Hernquist, Katz, & Weinberg (1995). Here we propose a test for this hypothesis using X-ray emission from CGs. The observable quantity Q identical with L(sub x)a(sup 3)(sup p)/(L(sup 2)(sup g))T(sup 1/2)(sub x) should be proportional to the axis ratio of the group, a/c, where a and c are the short and long axes of a prolate distribution, a(sub p) is the radius of the group projected onto the sky, L(sub x) is the bolometric X-ray luminosity, L(sub g) is the group blue luminosity, and T(sub x) is the gas temperature. We find that the distribution of Q is consistent with the notion that many CGs with unusually small values of a/c are frauds, i.e., that the values of Q are anomalously small. Absent other information, it is equally possible that CGs are very gas-poor relative to rich clusters; however, this can be tested using the Sunyaev-Zeldovich effect. If the groups have a close to normal ratio of gas to total mass, but are simply stretched out along the line of sight, a Sunyaev-Zeldovich signal should be detectable.

Ostriker, Jeremiah P.↗

Hot gas in superclusters and microwave background distortions

Variations in the electron pressure can produce significant distortions of the cosmic microwave background radiation (CMBR) spectrum through the Sunyaev-Zel'dovich effect. In this paper, we compute the amplitude of these distortions in various cosmological scenarios using the results of hydrodynamical simulations combined with analytic theory. For reionized cold dark matter (CDM) with h = 0.5 probed with the OVRO experiment (Theta approximately = 0.1 deg), we find that the amplitude of these fluctuations is (Delta T/T)(sub rms) = 1 x 10(exp -5) (Omega(sub i)/0.06) if the experiment looks at a random part of the sky and 6 x 10(exp -6)(Omega(sub i)/(0.06) if bright areas are avoided. Here Omega is the density of ionized baryons in units of the critical density, rho(sub c) = 2.78 x 10(exp 11)h(exp 2)/cu Mpc. Primary fluctuations would generate (Delta T/T)(sub rms) approximately = 2 x 10(exp -6) for this model. These distortions could be a dominant source of CMBR anisotropics on small and intermediate angular scales. It is important to note that these signals are generated relatively recently, when nonlinear effects are improtant, and do not require early reionization. We also examine the contribution to anisotropies from the Doppler effect, in both the linear and nonlinear regimes. In this case, the bulk of the signal is generated before a redshift of 10, provided that the intercluster medium is ionized back to the surface of last scattering. On the smallest scales thermal effects are bigger, while, on larger scales the Doppler-induced fluctuations are greater with the exact boundary determined by the scenario. However, over a wide range of scales, the two signals are equal to within a factor of a few.

Persi, Fred M.↗

[Kelvin-Helmholtz instability in protostellar jets]

NASA grant NAG 5 2866, funded by the Astrophysics Theory Program, enabled the study the Kelvin-Helmholtz instability in protostellar jets. In collaboration with co-investigator Philip Hardee, the PI derived the analytic dispersion relation for the instability in including a cooling term in the energy equation which was modeled as one of two different power laws. Numerical solutions to this dispersion relation over a wide range of perturbation frequencies, and for a variety of parameter values characterizing the jet (such as Mach number, and density ratio) were found It was found that the growth rates and wavelengths associated with unstable roots of the dispersion relation in cooling jets are significantly different than those associated with adiabatic jets, which have been studied previously. In collaboration with graduate student Jianjun Xu (funded as a research associate under this grant), hydrodynamical simulations were used to follow the growth of the instability into the nonlinear regime. It was found that asymmetric surface waves lead to large amplitude, sinusoidal distortions of the jet, and ultimately to disruption Asymmetric body waves, on the other hand, result in the formation of shocks in the jet beam in the nonlinear regime. In cooling jets, these shocks lead to the formation of dense knots and filaments of gas within the jet. For sufficiently high perturbation frequencies, however, the jet cannot respond and it remains symmetric. Applying these results to observed systems, such as the Herbig-Haro jets HH34, HH111 and HH47 which have been observed with the Hubble Space Telescope, we predicted that some of the asymmetric structures observed in these systems could be attributed to the K-H modes, but that perturbations on timescales associated with the inner disk (about 1 year) would be too rapid to cause disruption. Moreover, it was discovered that weak shock 'spurs' in the ambient gas produced by ripples in the jet surface due to nonlinear, modes of surface and/or body waves could accelerate the ambient gas to low velocity. This latter effect represents a new mechanism by which supersonic jets can accelerate low velocity outflows.

Stone, James↗

The Stability of Radiatively Cooling Jets. 2: Nonlinear Evolution

We use two-dimensional time-dependent hydrodynamical simulations to follow the growth of the Kelvin-Helmholtz (K-H) instability in cooling jets into the nonlinear regime. We focus primarily on asymmetric modes that give rise to transverse displacements of the jet beam. A variety of Mach numbers and two different cooling curves are studied. The growth rates of waves in the linear regime measured from the numerical simulations are in excellent agreement with the predictions of the linear stability analysis presented in the first paper in this series. In the nonlinear regime, the simulations show that asymmetric modes of the K-H instability can affect the structure and evolution of cooling jets in a number of ways. We find that jets in which the growth rate of the sinusoidal surface wave has a maximum at a so-called resonant frequency can be dominated by large-amplitude sinusoidal oscillations near this frequency. Eventually, growth of this wave can disrupt the jet. On the other hand, nonlinear body waves tend to produce low-amplitude wiggles in the shape of the jet but can result in strong shocks in the jet beam. In cooling jets, these shocks can produce dense knots and filaments of cooling gas within the jet. Ripples in the surface of the jet beam caused by both surface and body waves generate oblique shock "spurs" driven into the ambient gas. Our simulations show these shock "spurs" can accelerate ambient gas at large distances from the jet beam to low velocities, which represents a new mechanism by which low-velocity bipolar outflows may be driven by high-velocity jets. Rapid entrainment and acceleration of ambient gas may also occur if the jet is disrupted. For parameters typical of protostellar jets, the frequency at which K-H growth is a maximum (or highest frequency to which the entire jet can respond dynamically) will be associated with perturbations with a period of - 200 yr. Higher frequency (shorter period) perturbations excite waves associated with body modes that produce internal shocks and only small-amplitude wiggles within the jet. The fact that most observed systems show no evidence for large-amplitude sinusoidal oscillation leading to disruption is indicative that the perturbation frequencies are generally large, consistent with the suggestion that pro- tostellar jets arise from the inner regions (r less than 1 AU) of accretion disks.

Stone, James M.↗

Damped Ly alpha absorbers at high redshift: Large disks or galactic building blocks?

The nature of the physical structures giving rise to damped Lyman alpha absorption systems (DLAS) at high redshifts is investigated. The proposal that rapidly rotating large disks are the only viable explanation for the observed asymmetric profiles of low ionization absorption lines is examined. Using hydrodynamic simulations of galaxy formation, it is demonstated that irregular protogalactic clumps can reproduce the observed velocity width distribution and asymmetries of the absorption profiles equally well. The velocity broadening in the simulated clumps is due to a mixture of rotation, random motions, infall and merging. The observed velocity width correlates with the virial velocity for the dark matter halo of the forming protogalactic clump. The typical virial velocity of the halos required to lead to the DLAS population is approximately 100 km/s. It is concluded that the evidence that DLAS at high redshift are related to large, rapidly rotating disks, is not compelling.

Haehnelt, Martin G.↗

Statistical Properties of Line Centroid Velocity Increments in the rho Ophiuchi Cloud

We present a comparison of histograms of CO (2-1) line centroid velocity increments in the rho Ophiuchi molecular cloud with those computed for spectra synthesized from a three-dimensional, compressible, but non-starforming and non-gravitating hydrodynamic simulation. Histograms of centroid velocity increments in the rho Ophiuchi cloud show clearly non-Gaussian wings, similar to those found in histograms of velocity increments and derivatives in experimental studies of laboratory and atmospheric flows, as well as numerical simulations of turbulence. The magnitude of these wings increases monotonically with decreasing separation, down to the angular resolution of the data. This behavior is consistent with that found in the phase of the simulation which has most of the properties of incompressible turbulence. The time evolution of the magnitude of the non-Gaussian wings in the histograms of centroid velocity increments in the simulation is consistent with the evolution of the vorticity in the flow. However, we cannot exclude the possibility that the wings are associated with the shock interaction regions. Moreover, in an active starforming region like the rho Ophiuchi cloud, the effects of shocks may be more important than in the simulation. However, being able to identify shock interaction regions in the interstellar medium is also important, since numerical simulations show that vorticity is generated in shock interactions.

Lis, D. C.↗

Accretion Onto Magnetic Degenerate Stars

While the original objectives of this research program included the study of radiative processes in cataclysmic variables and the evolution of neutron star magnetic fields, the scope of the reported research expanded to other related topics as this project developed. This final report therefore describes the results of our research in the following areas: 1) Irradiation-driven mass transfer cycles in cataclysmic variables and low-mass X-ray binaries; 2) Propeller effect and magnetic field decay in isolated old neutron stars; 3) Decay of surface magnetic fields in accreting neutron stars and pulsars; 4) Finite-Difference Hydrodynamic simulations of mass transfer in binary stars.

Frank, Juhan↗

Comet Shoemaker-Levy 9, Jupiter, and Impact Shock Chemistry

Four years ago this month, a hitherto unknown comet in loose orbit around Jupiter passed so near the giant planet that it was torn apart into 20 fragments by tides. One orbit later, two years ago this month, the fragments of doomed comet P/Shoemaker Levy (SL9) fell into Jupiter. The enormous energies of these impacts (the largest fragments were nearly 1 km across and, hitting at 60 km/s, released some 2-4 x 10(exp 27) ergs) produced enormous explosions. Several of the ejecta plumes were imaged towering 3000 km above Jupiter's limb. The heat released when the plumes fell was considerable and easily observed on Earth. The impacts produced strong shocks, both promptly at the impact site and again, later, and over thousands of kilometers, when the ejecta plume reentered the atmosphere. The focus of this talk will be to discuss what the SL9 impacts taught us about impact shock chemistry - the processes, the ingredients, the results - and what inferences we may draw for impacts on early Earth. Shock chemistry generates a suite of molecules not usually seen on Jupiter. The most surprising report was of a huge amount of diatomic sulfur S2 at the site of the G impact. Other unusual products include CS, CS2, OCS, H2S, SO2, HCN, CO, and H2O; although H2S and H2O are doubtless abundant below the visible clouds. Hot or enhanced CH4 and NH3 were also detected. A general rule of shock chemistry is that CO forms until either C or O is exhausted. If O greater than C, the other products are oxidized, and excess O goes to H2O. If C greater than O, the other products are reduced, and excess C goes to HCN, C2H2, and a wide variety of more complicated organics. Ultimately, given time, the carbon would react all the way to graphite, but in practice the reactions are incomplete. The dark ejecta debris were probably composed in part of carbonaceous particles generated by the shocks. In a sense, the SL9 impacts performed the famous Miller-Urey experiment on a grand scale, with one result being the production of a lot of complex brown organic solids (called "tholins"). We use, a straightforward chemical kinetics model for the H, N, C, O, S system to follow the nonequilibrium chemistry behind the shocks. The model traces the evolving chemical composition of a parcel of gas by directly integrating the web of chemical reactions. Pressure and temperature histories of the parcels are patterned after those calculated by numerical hydrodynamic simulations of the ejecta plume. A given plume parcel is generally shocked twice; t.e a parcel shocked near the impact site is ejected at high velocity and is shocked again when it reenters the atmosphere. The final state of the gas depends mostly on the second shock, provided that the latter is hot enough. The chemical evidence is ambiguous, but most indications are that C greater than O in the shocked, reacting gas. Telltale signatures of abundant oxygen - SO2, SO, CO2, O2 - were not seen, while signatures of abundant carbon - CS, CS2, and HCN - were. On the other hand, abundant H2O would appear to require O greater than C, and two other observed sulfur species, S2 and OCS, appear to form more easily in a somewhat oxidized gas, presumable vaporized from the comet itself. Since on general principles one expects the -comet to have had a more-or-less cosmic composition, i.e. O greater than C, the production of CS, CS2, and HCN probably requires C greater than O in the shocked jovian air. This in turn implies that even the largest fragments released the bulk of their energy above the jovian water table, in all likelihood above 5 bars . There is no evidence in favor of the proposition that a significant amount of wet jovian air was shocked strongly enough to coax water to react; i.e. wet jovian air saw only temperatures significantly below 2000 K.

Zahnle, Kevin↗

The Soft X-Ray Emission in a Large Sample of Galaxy Clusters with ROSAT PSPC

The study of soft X-ray emission of 38 X-ray selected galaxy clusters observed by ROSAT PSPC indicates that the soft excess phenomenon may be a common occurrence in galaxy clusters. Excess soft X-ray radiation, above the contribution from the hot intra-cluster medium, is evident in a large fraction of sources, and is clearly detected with large statistical significance in the deepest observations. The investigation relies on new, high resolution 21 cm HI observations. The sample selection also features analysis of infrared images, to further ensure reliability of results with respect to the characteristics of Galactic absorption. The possibility of background or calibration effects as cause of the excess emission is likewise investigated; a detailed analysis of the distribution of the excess emission with respect to detector position and Galactic HI column density shows that the excess emission is a genuine celestial phenomenon. We find evidence for a preferential distribution of the soft excess emission at distances larger than approx. 150-200 kpc from the centers of clusters; this behavior may be naturally explained in the context of a non-thermal Inverse-Compton scenario. Alternatively, we propose that the phenomenon maybe caused by thermal emission of very large-scale 'warm' filaments seen in recent hydrodynamic simulations. This new interpretation relieves the very demanding requirements of either the traditional intra-cluster 'warm' gas and the non-thermal scenarios. We also investigate the possibility of the soft excess originating from unresolved, X-ray faint cluster galaxies.

Bonamente, Massimiliano↗

STIS Spectroscopy of the Lyman-Alpha Forest Toward 3C 273

We present results on the low-redshift Lyman-alpha forest as based on high-resolution (7 km/s) STIS spectra of 3C 273. A total of 121 intergalactic Lyman-alpha-absorbing systems were detected, of which 60 are above the 3.5sigma completeness limit, logNu(sub HI) approximately equals 12.3. The median line-width parameter, b = 27 km/s, is similar to that seen at high redshift. However the distribution of HI column densities has a steeper slope, beta = 2.02+/-0.21, than is seen at high redshift. Overall, the observed Nu(sub HI)-b distribution is consistent with that derived from a ACDM hydrodynamic simulation. We have used NED to compile a list of 300 galaxies (91 from SDSS, 98 from APM) within 1 Mpc of the line of sight to 3C 273 and are working to find line-of-sight velocity correlations between the galaxies and detected Lyman-alpha absorbers.

Heap, Sara R.↗

FUSE Observations of QSOs behind Galaxy Clusters and of Galactic O VI Emission

We observed five quasars located behind clusters of galaxies with the goal of detecting resonant O VI absorption from warm (T approx. 10(exp 6) K) gas in the clusters' intracluster medium. The presence of such warm gas is predicted by cosmological hydrodynamic simulations and is observationally supported by the detection of "soft excess" emission in several galaxy clusters. A second goal was the detection of diffuse O VI emission from warm gas in our own Galaxy.

Dixon, William V.↗

Gas Density Discontinuities in Merging Clusters

Chandra has discovered a new phenomenon in galaxy clusters, the sharp gas density edges. Depending on the sign of the temperature jump across the edge, these features may either be bow shocks or cold fronts. While bow shocks obviously are driven by merging sub-clusters, what causes cold fronts is not entirely clear, as they are observed both in mergers and in relaxed clusters. The purpose of the XMM study of A3376, an interesting cluster with density edges, is to understand the origin of cold fronts and to look for possible shocks. The XMM data for A3376 have been mostly analyzed (the X-ray edge turned out to be a cold front). Preliminary results have been shown at a conference and a paper is in preparation. We also have Chandra data for this cluster, and are comparing and combining the two datasets. In the course of analyzing the X-ray data for this cluster as well as several others, it has become apparent that we need the help of hydrodynamic simulations to study the precise mechanism by which cold fronts are formed, the main goal of the present project. A postdoc (Yago Ascasibar) is currently running SPH simulations of an idealized sub- cluster merger. These advanced simulations are nearing completion and two papers with their results are in preparation.

Mushotzky, Richard↗

The X-ray halo of an extremely luminous LSB disk galaxy

We are continuing to refine our upper limit on emission from halo gas in Malin 2. The upper limit is, of course, below the detected flux, but is made more difficult to quantify by the disk and possible AGN sources. We are also exploring spectral and spatial-size constraints to help separate the sources of emission. On the theory side, more recent work on the X-ray halo luminosity from halo gas leftover from galaxy formation has lowered the prediction for disk galaxies (e.g. Toft et al. 2002, MNRAS, 335, 799). While our upper limit is well below the original prediction, refinements in model have moved the theoretical goalposts, so that the observation may be consistent with newer models. A recent theoretical development, which our observations of Malin 2 appear to support, is that a substantial amount of mass can be accreted onto galaxies without being heated at a virial shock. The previous standard theory was that gas accreting into a halo hits a virial shock and is heated to high temperatures, which could produce X-ray halos in massive galaxies. Recent models show that "smooth accretion" of matter bypasses the virial shocking (Murali e t al. 2002, ApJ, 571, 1; Birnboim & Dekel 2003, MNRAS, 345, 349). Additionally, new hydrodynamical simulations of galaxy mergers by UCSC graduate student T. J. Cox show that hot gas halos can be created by gas blown out from the merger, taking up orbital energy of the merging galaxies (Cox et al. 2004, ApJ, 607, L87). If mergers rather than virial shocking are the origin of hot gas halos, the existence of an X-ray halo should depend more on past merger activity than halo mass. Then it makes sense that elliptical galaxies and poor groups with ellipticals, which are probably formed in mergers, have X-ray gas halos; while a giant, quiescent LSB disk galaxy like Malin 2, which has never suffered a major merger, does not have an X-ray halo. While both the observational expectations and theoretical models have changed since we began this project, which has forced us to re-evaluate the goals, we are pressing forward to firm up the observations and put them in context of the current models for X-ray halos of massive galaxies.

Weiner, Benjamin J.↗

Velocity Characteristics of Evaporated Plasma using Hinode/EIS

This paper presents a detailed study of chromospheric evaporation using the EUV Imaging Spectrometer (EIS) onboard Hinode in conjunction with HXR observat,ions from RHESSI. The advanced capabilities of EIS were used to measure Doppler shifts in 15 emission lines covering the temperature range T=0.05-16 MK during the impulsive phase of a C-class flare on 2007 December 14. Blueshifts indicative of the evaporated material were observed in six emission lines from Fe XIV-XXIV (2-16 MK). Upflow velocity was found to scale with temperature as v(sub up) (kilometers per second) approximately equal to 5-17 T (MK). Although the hottest emission lines, Fe XXIII and Fe XXIV, exhibited upflows of greater than 200 kilometers per second, their line profiles were found to be dominated by a stationary component in stark contrast to the predictions of the standard flare model. Emission from O VI-Fe XIII lines (0.5-1.5 MK) was found to be redshifted by v(sub down) (kilometers per second) approximately equal to 60-17 T (MK) and was interpreted as the downward-moving 'plug' characteristic of explosive evaporation. These downflows occur at temperatures significantly higher than previously expected. Both upflows and downflows were spatially and temporally correlated with HXR emission observed by RHESSI that provided the properties of the electron beam deemed to be the driver of the evaporation. The energy contained in the electron beam was found to be greater than or equal to 10(sup 11) ergs per square centimeter per second consistent with the value required to drive explosive chromospheric evaporation from hydrodynamic simulations.

Milligan, Ryan O.↗