Search NASA⌕ Search

SEARCH · Search NASA

Results for “Hydrodynamics simulations”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Wind dynamics in SMC X-1. 1: Hydrodynamic simulation

We present a three-dimensional hydrodynamic simulation of the disrupted stellar wind in the high-mass X-ray binary system SMC X-1. The three dominant processes that determine the geometry of the wind in high X-ray luminosity systems such as SMC X-1 are the X-ray suppression of the stellar wind from the X-ray irradiated face of the primary star, the focusing of the radiatively driven wind in the X-ray shadow by the effects of stellar rotation, and the rapid X-ray heating of gas in the vicinity of the X-ray source, including the X-ray illuminated surface of the primary star. The resulting distribution of circumstellar gas provides a successful explanation for the asymmetric, extended eclipse transitions and the intensity of the deep eclipse X-ray emission in SMC X-1, as well as a possible explanation for the X-ray dips seen near superior conjunction of the X-ray source in Cyg X-1.

Blondin, John M.↗
NASA NTRS2007

Hydrodynamical Simulations of the Jet in the Symbiotic Star MWC 560: Application to X-ray Jets in Symbiotic Stars - 3

In Papers I and II in this series, we presented hydrodynamical simulations of jet models with parameters representative of the symbiotic system MWC 560. These were simulations of a pulsed, initially underdense jet in a high-density ambient medium. Since the pulsed emission of the jet creates internal shocks and since the jet velocity is very high, the jet bow shock and the internal shocks are heated to high temperatures and should therefore emit X-ray radiation. In this paper, we investigate in detail the X-ray properties of the jets in our models. We have focused our study on the total X-ray luminosity and its temporal variability, the resulting spectra, and the spatial distribution of the emission. Temperature and density maps from our hydrodynamical simulations with radiative cooling presented in the second paper are used, together with emissivities calculated with the atomic database ATOMDB. The jets in our models show extended and variable X-ray emission, which can be characterized as a sum of hot and warm components with temperatures that are consistent with observations of CH Cyg and R Aqr. The X-ray spectra of our model jets show emission-line features that correspond to observed features in the spectra of CH Cyg. The innermost parts of our pulsed jets show iron line emission in the 6.4-6.7 keV range, which may explain such emission from the central source in R Aqr. We conclude that MWC 560 should be detectable with Chandra or XMM-Newton, and such X-ray observations will prove crucial for understanding jets in symbiotic stars.

hydrodynamics↗
NASA NTRS2023

Optimizing the Resolution of Hydrodynamic Simulations for MCRaT Radiative Transfer Calculations

Despite their discovery about half a century ago, the Gamma-ray burst (GRB) prompt emission mechanism is still not well understood. Theoretical modeling of the prompt emission has advanced considerably due to new computational tools and techniques. One such tool is the PLUTO hydro dynamics code, which is used to numerically simulate GRB outflows. PLUTO uses Adaptive Mesh Refinement to focus computational efforts on the portion of the grid that contains the simulated jet. Another tool is the Monte Carlo Radiation Transfer (MCRaT) code, which predicts electromagnetic signatures of GRBs by conducting photon scatterings within a jet using PLUTO. The effects of the underlying resolution of a PLUTO simulation with respect to MCRaT post-processing radiative transfer results have not yet been quantified. We analyze an analytic spherical outflow and a hydrodynamically simulated GRB jet with MCRaT at varying spatial and temporal resolutions and quantify how decreasing both resolutions affect the resulting mock observations. We find that changing the spatial resolution changes the hydrodynamic properties of the jet, which directly affect the MCRaT mock observable peak energies. We also find that decreasing the temporal resolution artificially decreases the high energy slope of the mock observed spectrum, which increases both the spectral peak energy and the luminosity. We show that the effects are additive when both spatial and temporal resolutions are modified. Our results allow us to understand how decreased hydrodynamic temporal and spatial resolutions affect the results of post-processing radiative transfer calculations, allowing for the optimization of hydrodynamic simulations for radiative transfer codes.

Jose Arita-Escalante↗
NASA NTRS1996

Comparison of Non-Parabolic Hydrodynamic Simulations for Semiconductor Devices

Parabolic drift-diffusion simulators are common engineering level design tools for semiconductor devices. Hydrodynamic simulators, based on the parabolic band approximation, are becoming more prevalent as device dimensions shrink and energy transport effects begin to dominate device characteristic. However, band structure effects present in state-of-the-art devices necessitate relaxing the parabolic band approximation. This paper presents simulations of ballistic diodes, a benchmark device, of Si and GaAs using two different non-parabolic hydrodynamic formulations. The first formulation uses the Kane dispersion relationship in the derivation of the conservation equations. The second model uses a power law dispersion relation {(hk)(exp 2)/2m = xW(exp Y)}. Current-voltage relations show that for the ballistic diodes considered. the non-parabolic formulations predict less current than the parabolic case. Explanations of this will be provided by examination of velocity and energy profiles. At low bias, the simulations based on the Kane formulation predict greater current flow than the power law formulation. As the bias is increased this trend changes and the power law predicts greater current than the Kane formulation. It will be shown that the non-parabolicity and energy range of the hydrodynamic model based on the Kane dispersion relation are limited due to the binomial approximation which was utilized in the derivation.

Smith, A. W.↗
NASA NTRS1988

Analysis of the Space Shuttle Orbiter skin panels under simulated hydrodynamic loads

The Space Shuttle orbiter skin panels were analyzed under pressure loads simulating hydrodynamic loads to determine their capability to sustain a potential ditching and to determine pressures that typically would produce failures. Two Dynamic Crash Analysis of Structures (DYCAST) finite element models were used. One model was used to represent skin panels (bays) in the center body, while a second model was used to analyze a fuselage bay in the wing region of the orbiter. From an assessment of the DYCAST nonlinear computer results, it is concluded that the probability is extremely high that most, if not all, of the lower skin panels would rupture under ditching conditions. Extremely high pressure loads which are produced under hydrodynamic planning conditions far exceed the very low predicted failure pressures for the skin panels. Consequently, a ditching of the orbiter is not considered to have a high probability of success and should not be considered a means of emergency landing unless no other option exists.

Carden, Huey D.↗
NASA NTRS1994

Two-dimensional hydrodynamical simulations of wind-compressed disks around rapidly rotating B stars

We use a two-dimensional piecewise parabolic method (PPM) code to simulate numerically the hydrodynamics of a radiation-driven stellar wind from a rapidly rotating Be star. The results generally confirm predictions of the semianalytic 'wind-compressed disk' model recently proposed by Bjorkman and Cassinelli to explain the circumstellar disks inferred observationally to exist around such rapidly rotating stars. However, this numerical simulation is able to incorporate several important effects not accounted for in the simple model, including a dynamical treatment of the outward radiative driving and gas pressure, as well as a rotationally distorted, oblate stellar surface. This enables us to model quantitatively the compressed wind and shock that forms the equatorial disk. The simulation results thus do differ in several important details from the simple method, showing, for example, an inner disk inflow not possible in the heuristic approach of assuming a fixed outward velocity law. There is also no evidence for the predicted detachment of the disk that arises in the fixed outflow picture. The peak equatorward velocity in the dynamical models is furthermore about a factor of 2 smaller than the lytically predicted value of approximately 50% of the stellar equatorial rotation speed. As a result, the dynamical disks are somewhat weaker than predicted, with a wider opening angle, lower disk/pole density ratio, and smaller shock velocity jump. The principal cause of these latter differences appears to be an artificially strong equatorward drift of the subsonic outflow in the original analytic model. Much better agreement with the dynamical results can be obtained, however, from a slightly modified, analytic wind-compression model with a more detailed specification of the fixed wind outflow and a lower boundary set to the sonic radius along a rotationally oblate stellar surface. Hence, despite these detailed differences, the general predicted effect of disk formation by wind compression toward the equator is substantially confirmed.

Owocki, Stanley P.↗
NASA NTRS1991

Radiative-hydrodynamical simulations of accretion disk coronae

An algorithm has been developed that will carry out detailed, two-dimensional, fully self-consistent radiative-hydrodynamical simulations of accretion disk coronae in X-ray-luminous compact sources. The calculation reported here, for an accreting neutron star radiating at 0.5 times the Eddington luminosity, reveals several striking features. (1) The corona is comprised of two main regions - an inner highly dynamic portion whose vertical structure varies cyclically on a dynamical time scale, and an outer, more stable zone in which the evaporated plasma rises to form a 'sheath' that gradually merges into a wind at large radii. (2) The coronal structure shows a density inversion that contrasts sharply with the Gaussian profiles of earlier hydrostatic models. (3) Interestingly, flow velocities as high as a few billion cm/s are not uncommon in portions of the corona. The relevance of this point to the large measured width of the emission features seen in low-mass X-ray binaries and the Galactic black hole candidates is discussed.

Melia, Fulvio↗
NASA NTRS1979

Hydrodynamic simulations of flare/surge events

A one-dimensional, hydrodynamic, time-dependent model which simulates the major observed dynamics of flare associated surges was developed. The thermodynamics, the surge mass, the time scales, the physical dimensions, and the velocities of typical surge events are reproduced. The surge is created by a sudden pressure increase at the top of the chromosphere, and this pressure pulse produces a disturbance which is followed, with a time-dependent numerical solution, as it propagates upward through the transition region and into the corona. The leading edge of the disturbance is a weak shock which has only a slight effect on the original transition region and coronal thermodynamics. The major effect occurs behind the shock where the temperature is decreased and the density is increased, and this cool, dense region moves upward then falls downward. After the material begins returning to the chromosphere, a second shock is formed which propagates upward, brings the infalling material to rest, and returns the atmosphere to hydrostatic equilibrium.

Steinolfson, R. S.↗
NASA NTRS2013

Constraints on Decreases in Eta Carinae's Mass-loss from 3D Hydrodynamic Simulations of Its Binary Colliding Winds

Recent work suggests that the mass-loss rate of the primary star Eta-A in the massive colliding wind binary Eta Carinae dropped by a factor of 2-3 between 1999 and 2010. We present result from large- (+/- 1545 au) and small- (+/- 155 au) domain, 3D smoothed particle hydrodynamics (SPH) simulations of Eta Car's colliding winds for three Eta-A mass-loss rates ( (dot-M(sub Eta-A) = 2.4, 4.8 and 8.5 × 10(exp −4) M(solar)/ yr), investigating the effects on the dynamics of the binary wind-wind collision (WWC). These simulations include orbital motion, optically thin radiative cooling and radiative forces. We find that dot-M Eta-A greatly affects the time-dependent hydrodynamics at all spatial scales investigated. The simulations also show that the post-shock wind of the companion star Eta-B switches from the adiabatic to the radiative-cooling regime during periastron passage (Phi approx.= 0.985-1.02). This switchover starts later and ends earlier the lower the value of dot-M Eta-A and is caused by the encroachment of the wind of Eta-A into the acceleration zone of Eta-B's wind, plus radiative inhibition of Eta-B's wind by Eta-A. The SPH simulations together with 1D radiative transfer models of Eta-A's spectra reveal that a factor of 2 or more drop in dot-M EtaA should lead to substantial changes in numerous multiwavelength observables. Recent observations are not fully consistent with the model predictions, indicating that any drop in dot- M Eta-A was likely by a factor of approx. < 2 and occurred after 2004. We speculate that most of the recent observed changes in Eta Car are due to a small increase in the WWC opening angle that produces significant effects because our line of sight to the system lies close to the dense walls of the WWC zone. A modest decrease in dot-M Eta-A may be responsible, but changes in the wind/stellar parameter of Eta-B, while less likely, cannot yet be fully ruled out. We suggest observations during Eta-Car's next periastron in 2014 to further test for decreases in dot-M Eta-A. If dot-M Eta-A is declining and continues to do so, the 2014 X-ray minimum should be even shorter than that of 2009.

Eta Carinae↗
NASA NTRS1994

X-ray clusters in a cold dark matter + lambda universe: A direct, large-scale, high-resolution, hydrodynamic simulation

A new, three-dimensional, shock-capturing, hydrodynamic code is utilized to determine the distribution of hot gas in a cold dark matter (CDM) + lambda model universe. Periodic boundary conditions are assumed: a box with size 85/h Mpc, having cell size 0.31/h Mpc, is followed in a simulation with 270(exp 3) = 10(exp 7.3) cells. We adopt omega = 0.45, lambda = 0.55, h identically equal to H/100 km/s/Mpc = 0.6, and then, from the cosmic background explorer (COBE) and light element nucleosynthesis, sigma(sub 8) = 0.77, omega(sub b) = 0.043. We identify the X-ray emitting clusters in the simulation box, compute the luminosity function at several wavelength bands, the temperature function and estimated sizes, as well as the evolution of these quantities with redshift. This open model succeeds in matching local observations of clusters in contrast to the standard omega = 1, CDM model, which fails. It predicts an order of magnitude decline in the number density of bright (h nu = 2-10 keV) clusters from z = 0 to z = 2 in contrast to a slight increase in the number density for standard omega = 1, CDM model. This COBE-normalized CDM + lambda model produces approximately the same number of X-ray clusters having L(sub x) greater than 10(exp 43) erg/s as observed. The background radiation field at 1 keV due to clusters is approximately the observed background which, after correction for numerical effects, again indicates that the model is consistent with observations.

Cen, Renyue↗
NASA NTRS1992

Radiative-hydrodynamical simulations of X-ray burst-induced accretion disk coronae

A fully self-consistent radiative hydrodynamical algorithm developed earlier is used to study the 'steady' state coronae in order to simulate the response of the evaporated outflow to a time-dependent irradiation of the disk's surface. It is found that the burst definition is altered at both low and high inclination angles i. At i greater than or approximately equal to 45 deg, the detectable flux is enhanced (as much as 50 percent) by X-rays scattered into the line of sight from the radiation propagation through the corona, and (greater than or approximately equal to 30 percent) by X-rays reflected off the corona when i is less than or approximately equal to 45 deg.

Mella, Fulvio↗
NASA NTRS1992

Three-dimensional hydrodynamical simulations of colliding stars . III - Collisions and tidal captures of unequal-mass main-sequence stars

A 3D smooth-particle hydrodynamic code with 7000 particles is used to simulate collisions between lower main sequence (MS) stars whose masses differ by a factor of five. Collisions between point-mass intruders and MS stars are simulated to study how the physics of the encounters depends on the finite radius of the smaller MS stars. It is found that the two MS stars become gravitationally bound in a physical collision if their relative velocity at infinity is less than a critical velocity Vd. Vd decreases from 1000 km/s in a head-on collision to 150 km/s in a grazing one. If the less massive star is replaced by a point mass, Vd remains the same for grazing collision and tidal encounters but drops to about 600 km/s in head-on collisions. The collision and coalescence of two MS stars resets the nuclear clock of the coalesced star, so it contracts to the MS. Tidal capture of binary stars in globular clusters can only occur in encounters in which the closest approach of the two stars to their center of mass is less than 2.0 times the sum of their radii.

Benz, W.↗
NASA NTRS1998

Three-Dimensional Hydrodynamic Simulations of the Postimpact Proto-Earth

Astrophysical fluid configurations are susceptible to a variety of nonaxisymmetric instabilities under the combined effects of rotation, self-gravity, and thermal pressure. When strong enough, they can induce rapid transport of mass and angular momentum. Our own previous studies of nonaxisymmetric instabilities in model protostars and protostellar disks show that significant transport can occur on orbital timescales and that material can be ejected to large distances. In this contribution, we present three-dimensional simulations of the circumterrestrial debris belt that may have resulted from a giant impact. Our three-dimensional hydrodynamics code with self-gravity and artificial viscosity is fully second-order in space and time; the equations of hydrodynamics and the Poisson equation are solved on an Eulerian cylindrical grid. In the preliminary calculations presented here, we use a simplified EOS where the central proto-Earth is treated as an n = 1/2 polytropic fluid, surrounded by a more compressible, rapidly rotating, fluid disk that represents silicate vapor. Our initial disk parameters are generated from the endstate data of recent smoothed particle hydrodynamics giant-impact calculations. Ultimately, we wish to detennine under what conditions nonaxisymmetric instabilities grow in the postimpact disk and whether they facilitate the transport of material outside the proto-Earth's Roche Limit, leading to the formation of the Moon. In future work, we hope to include a more realistic EOS and the consequences of heating, cooling, and phase transitions.

Pickett, B. K.↗
NASA NTRS1987

Three-dimensional hydrodynamical simulations of stellar collisions. I - Equal-mass main-sequence stars

Two distinct mass-loss mechanisms are noted in the present, fully three-dimensional calculations of collisions between identical stars. While strong shocks in nearly head-on collisions lead to high-velocity jets perpendicular to the collision axis, with increasing mass loss as impact velocity at infinity increases from zero to 2.3 times the escape velocity from the stellar surface, low velocity encounters lead to a sharp increase in mass loss at impact parameters that correspond to nearly-grazing collisions in a two-stage process. In the first stage, the two stars become gravitationally bound due to the encounter's energy dissipation; these binary components then violently coalesce during subsequent periastron passage.

Benz, Willy↗
NASA NTRS1990

Hydrodynamic simulations of stellar wind disruption by a compact X-ray source

This paper presents two-dimensional numerical simulations of the gas flow in the orbital plane of a massive X-ray binary system, in which the mass accretion is fueled by a radiation-driven wind from an early-type companion star. These simulations are used to examine the role of the compact object (either a neutron star or a black hole) in disturbing the radiatively accelerating wind of the OB companion, with an emphasis on understanding the origin of the observed soft X-ray photoelectric absorption seen at late orbital phases in these systems. On the basis of these simulations, it is suggested that the phase-dependent photoelectric absorption seen in several of these systems can be explained by dense filaments of compressend gas formed in the nonsteady accreation bow shock and wake of the compact object.

Blondin, John M.↗
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 NTRS2003

Impact Tsunami Calculations: Hydrodynamical Simulations vs. Linear Theory

Tsunamis generated by the impacts of asteroids and comets into the Earth oceans are widely recognized as a potential catastrophic hazard to the Earth s population. Our general conclusion is that linear theory is a reasonably accurate guide to behavior of tsunamis generated by impactors of moderate size, where the initial transient impact cavity is of moderate depth compared to the ocean depth. This is particularly the case for long wavelength waves that propagate fastest and would reach coastlines first. Such tsunamis would be generated in the open ocean by impactors of 300 meters in diameter, which might be expected to strike the Earth once every few thousand years, on the average. Larger impactors produce cavities deep enough to reach the ocean floor; even here, linear theory is applicable if the starting point is chosen at a later phase in the calculation when the impact crater has slumped back to produce a cavity of moderate depth and slope.

Korycansky, E.↗