Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock hydrodynamics”

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

Galactic scale gas flows in colliding galaxies: 3-dimensional, N-body/hydrodynamics experiments

We present some results from three dimensional computer simulations of collisions between models of equal mass galaxies, one of which is a rotating, disk galaxy containing both gas and stars and the other is an elliptical containing stars only. We use fully self consistent models in which the halo mass is 2.5 times that of the disk. In the experiments we have varied the impact parameter between zero (head on) and 0.9R (where R is the radius of the disk), for impacts perpendicular to the disk plane. The calculations were performed on a Cray 2 computer using a combined N-body/smooth particle hydrodynamics (SPH) program. The results show the development of complicated flows and shock structures in the direction perpendicular to the plane of the disk and the propagation outwards of a density wave in both the stars and the gas. The collisional nature of the gas results in a sharper ring than obtained for the star particles, and the development of high volume densities and shocks.

Lamb, Susan A.↗

Colliding winds from early-type stars in binary systems

The dynamics of the wind and shock structure formed by the wind collision in early-type binary systems is examined by means of a 2D hydrodynamics code, which self-consistently accounts for radiative cooling, and represents a significant improvement over previous attempts to model these systems. The X-ray luminosity and spectra of the shock-heated region, accounting for wind attenuation and the influence of different abundances on the resultant level and spectra of X-ray emission are calculated. A variety of dynamical instabilities that are found to dominate the intershock region is examined. These instabilities are found to be particularly important when postshock material is able to cool. These instabilities disrupt the postshock flow and add a time variability of order 10 percent to the X-ray luminosity. The X-ray spectrum of these systems is found to vary with the nuclear abundances of winds. These theoretical models are used to study several massive binary systems, in particular V444 Cyg and HD 193793.

Stevens, Ian R.↗

New Evidence for Efficient Collisionless Heating of Electrons at the Reverse Shock of a Young Supernova Remnant

Although collisionless shocks are ubiquitous in astrophysics, certain key aspects of them are not well understood. In particular, the process known as collisionless electron heating, whereby electrons are rapidly energized at the shock front, is one of the main open issues in shock physics. Here, we present the first clear evidence for efficient collisionless electron heating at the reverse shock of Tycho's supernova remnant (SNR), revealed by Fe K diagnostics using high-quality X-ray data obtained by the Suzaku satellite. We detect K beta (3p yields 1s) fluorescence emission from low-ionization Fe ejecta excited by energetic thermal electrons at the reverse shock front, which peaks at a smaller radius than Fe K alpha (2p yields 1s) emission dominated by a relatively highly ionized component. Comparisons with our hydrodynamical simulations imply instantaneous electron heating to a temperature 1000 times higher than expected from Coulomb collisions alone. The unique environment of the reverse shock, which is propagating with a high Mach number into rarefied ejecta with a low magnetic field strength, puts strong constraints on the physical mechanism responsible for this heating and favors a cross-shock potential created by charge deflection at the shock front. Our sensitive observation also reveals that the reverse shock radius of this SNR is about 10% smaller than the previous measurement using the Fe K alpha morphology from the Chandra observations. Since strong Fe K beta fluorescence is expected only from low-ionization plasma where Fe ions still have many 3p electrons, this feature is key to diagnosing the plasma state and distribution of the immediate postshock ejecta in a young SNR.

Fe K diagnostics↗

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

Non-LTE radiating acoustic shocks and Ca II K2V bright points

We present, for the first time, a self-consistent solution of the time-dependent 1D equations of non-LTE radiation hydrodynamics in solar chromospheric conditions. The vertical propagation of sinusoidal acoustic waves with periods of 30, 180, and 300 s is calculated. We find that departures from LTE and ionization recombination determine the temperature profiles of the shocks that develop. In LTE almost all the thermal energy goes into ionization, so the temperature rise is very small. In non-LTE, the finite transition rates delay the ionization to behind the shock front. The compression thus goes into thermal energy at the shock front leading to a high temperature amplitude. Further behind the shock front, the delayed ionization removes energy from the thermal pool, which reduces the temperature, producing a temperature spike. The 180 s waves reproduce the observed temporal changes in the calcium K line profiles quite well. The observed wing brightening pattern, the violet/red peak asymmetry and the observed line center behavior are all well reproduced. The short-period waves and the 5 minute period waves fail especially in reproducing the observed behavior of the wings.

Carlsson, Mats↗

Formation of standing shocks in stellar winds and related astrophysical flows

Stellar winds and other analogous astrophysical flows can be described, to lowest order, by the familiar one dimensional hydrodynamic equations which, being nonlinear, admit in some instances discontinuous as well as continuous transonic solutions for identical inner boundary conditions. The characteristics of the time dependent differential equations of motion are described to show how a perturbation changes profile in time and, under well defined conditions, develops into a stationary shock discontinuity. The formation of standing shocks in wind type astrophysical flows depends on the fulfillment of appropriate necessary conditions, which are determined by the conservation of mass, momentum and energy across the discontinuity, and certain sufficient conditions, which are determined by the flow's history.

Tsinganos, K.↗

The structure of mass-loading shocks. 1: Comets

A new multifluid model to describe the solar wind-cometary ion plasma in the outer cometary coma is derived. This model is distinguished from previous multifluid models in that wave-particle effects are included explicitly. By considering hydrodynamic timescales, self-consistent dissipative terms (which correspond to the spatial diffusion of cometary ions) are derived, thereby allowing the structure of a cometary bow shock to be resolved. In the hypersonic limit, our diffusive multifluid model reveals an attractive connection to the one-fluid model formulated by Biermann et al. (1967), a connection which serves as an important guide when investigating our more complicated nonhypersonic diffusive model. The detailed model consists of cometary ions diffusing in a warm solar wind which is heated by the dissipation of Alfvenic turbulence (in situ and generated by the initial ring beam distribution). The quasi-parallel cometary shock is found to be at least an order of magnitude thicker than a quasi-perpendicular cometary shock. It is furthermore found that for undisturbed solar wind (sonic) Mach numbers in excess of approximately 6, the cometary bow shock (for both perpendicular and parallel cases) is smoothed completely by the cometary ions and no 'proton' or thermal subshock is necessary. Finally, a detailed comparison of the model with the plasma observations made at Halley by Giotto is undertaken. It is found that the observed shock structure, thickness, location, and plasma parameters compare well with the predictions of the theoretical model, particularly in the case of the quasi-parallel shock.

Zank, G. P.↗

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↗

Instability modes excited by natural screech tones in a supersonic rectangular jet

The evolution of hydrodynamic instability modes self-excited by harmonically related natural screech tones was experimentally investigated. A convergent rectangular nozzle with an aspect ratio of 9.63 was used to produce a supersonic shock containing jet. Measurements in the flow-field were made using standard hot-film probes positioned only in the subsonic (outer) portions of the flow. The hydrodynamic instability mode observed in the shear layer at the screech frequency was observed to be antisymmetric (sinuous) about the smaller dimension of the jet, whereas its harmonic was observed to be symmetric (varicose). In addition, the near-field noise measurements indicated that the radiated screech tone noise was out of phase on either side of the small jet dimension whereas its harmonic was in phase over the same region. To our knowledge such an observation on the nature of the harmonic has thus far gone unreported and therefore is the focus of the present work. The hydrodynamic instability modes occurring at the screech frequency and its harmonic satisfied the conditions for resonance. Detailed measurements of the coherent wave evolution in the streamwise and spanwise directions indicated that strong spanwise variations were present beyond x/h = 8. Details of the screech noise radiated by the coherent instability modes are also presented in this paper.

Raman, Ganesh↗

Hydrodynamic solutions at a sonic-boom focus

Study of the focusing of sonic-boom N waves using a numerical shock-following method. The computer code used in the calculations employs a moving mesh so that weak N wave shocks can be followed for large distances without numerical dissipation. In several sample calculations, overpressures at the focus ranging from 2 to 20 times the nominal overpressure have been computed. It is also found that weak shocks tend to cusp at the axis in accord with the geometric-acoustics description, whereas strong shocks straighten out without cusping.

Parker, L. W.↗

Computational Relativistic Astrophysics Using the Flow Field-Dependent Variation Theory

We present our method for solving general relativistic nonideal hydrodynamics. Relativistic effects become pronounced in such cases as jet formation from black hole magnetized accretion disks which may lead to the study of gamma-ray bursts. Nonideal flows are present where radiation, magnetic forces, viscosities, and turbulence play an important role. Our concern in this paper is to reexamine existing numerical simulation tools as to the accuracy and efficiency of computations and introduce a new approach known as the flow field-dependent variation (FDV) method. The main feature of the FDV method consists of accommodating discontinuities of shock waves and high gradients of flow variables such as occur in turbulence and unstable motions. In this paper, the physics involved in the solution of relativistic hydrodynamics and solution strategies of the FDV theory are elaborated. The general relativistic astrophysical flow and shock solver (GRAFSS) is introduced, and some simple example problems for computational relativistic astrophysics (CRA) are demonstrated.

Richardson, G. A.↗

Interacting winds and the shaping of planetary nebulae

Using two-dimensional hydrodynamics, the problem of shaping of planetary nebulae by the interacting winds model is studied. The shock structure is obtained and it is shown that interacting winds are capable of producing morphologies similar to the ones observed in planetary nebulae, provided that a density contrast exists between the equatorial and polar directions. This confirms recent suggestions by Balick. Mechanisms that can produce the required density contrast are discussed and it is shown that binary central stars provide the most natural mechanism, especially via common envelope evolution. It is not yet entirely clear whether stellar rotation is sufficient to generate the required contrast in the case of single central stars.

Soker, Noam↗

Numerical calculation of shock-induced initiation of detonations

Results of numerical calculations of the impact of steel cylinders and spheres on the plastic bonded high explosive PBX 9501 are described. The calculations were carried out by a reactive, multicomponent, two dimensional, Eulerian hydrodynamic computer code, 2DE. The 2DE computer code is a finite difference code that uses the donor acceptor cell method to compute mixed cell fluxes. The parameters in the Forest Fire burn model are developed from experiments where the induced shock approximates a plane wave and are applied, in this case, to a situation where the induced shock is a divergent wave with curvature that depends on the size and shape of the projectile. The calculated results are compared with results from experiments involving instrumented mock and live high explosives, with projectiles of varying size, shapes, and velocities.

Cort, G. E.↗

Interpretation of OAO-2 ultraviolet light curves of beta Doradus

Middle-ultraviolet light curves of beta Doradus, obtained by OAO-2, are presented along with other evidence indicating that the small additional bumps observed on the rising branches of these curves have their origin in shock-wave phenomena in the upper atmosphere of this classical Cepheid. A simple piston-driven spherical hydrodynamic model of the atmosphere is developed to explain the bumps, and the calculations are compared with observations. The model is found to be consistent with the shapes of the light curves as well as with measurements of the H-alpha radial velocities.

Hutchinson, J. L.↗

Lunar and Planetary Science XXXV: Origin of Planetary Systems

The session"Origin of Planetary Systems" presented the following reports:Planetary Formation in the the Gamma Cephei System by Core:Accretion; Accretion and Heating of Particles by Supersonic Planetesimals; Planetesimal Accretion in Close Binary Systems; Collisions, Gas Flow, and the Formation of Planetesimals; Hydrodynamic Escape of a Proto-Atmosphere Just After a Giant Impact; X-Ray Flare Induced Shock Waves and Chondrule Formation in Upper Solar Nebula; and Making Water Worlds: The Role of 26Al.

Source record↗

Post-shock temperatures in minerals

An experimental technique for the measurement of post-shock temperatures in minerals is presented, and silicate post-shock temperature measurements are compared with theoretical calculations. The technique involves the use of infrared detectors to determine the brightness temperatures at 4.5 to 5.75 microns and 7 to 14 microns of samples shocked to between 5 and 30 GPa. Results obtained for aluminum 2024 and stainless steel 304, as well as for Bamble bronzite and synthetic crystal forsterite are found at low pressures to be considerably in excess of the temperatures predicted assuming a hydrodynamic rheology and isentropic release parallel to the Hugoniot. The results are shown to be in better agreement, however, with values calculated assuming elastoplastic behavior, and the post-shock temperatures of crystalline quartz are found to be in good agreement with those calculated by Mashimo et al. (1979) from release adiabat data.

Raikes, S. A.↗

A High-Resolution X-Ray and Optical Study of SN1006: Asymmetric Expansion and Small-Scale Structure in a Type Ia Supernova Remnant

We introduce a deep (670 ks) X-ray survey of the entire SN 1006 remnant from the Chandra X-Ray Observatory, together with a deep Hα image of SN 1006 from the 4 m Blanco telescope at CTIO. Comparison with Chandra images from 2003 gives the first measurement of the X-ray proper motions around the entire periphery, carried out over a 9 yr baseline. We find that the expansion velocity varies significantly with azimuth. The highest velocity of approx.7400 km/s (almost 2.5 times that in the northwest (NW)) is found along the southeast (SE) periphery, where both the kinematics and the spectra indicate that most of the X-ray emission stems from ejecta that have been decelerated little, if at all. Asymmetries in the distribution of ejecta are seen on a variety of spatial scales. Si-rich ejecta are especially prominent in the SE quadrant, while O and Mg are more uniformly distributed, indicating large-scale asymmetries arising from the explosion itself. Neon emission is strongest in a sharp filament just behind the primary shock along the NWrim, where the pre-shock density is highest. Here the Ne is likely interstellar, while Ne within the shell may include a contribution from ejecta. Within the interior of the projected shell we find a few isolated "bullets" of what appear to be supernova ejecta that are immediately preceded by bowshocks seen in Hα, features that we interpret as ejecta knots that have reached relatively dense regions of the surrounding interstellar medium, but that appear in the interior in projection. Recent three-dimensional hydrodynamic models for Type Ia supernovae display small-scale features that strongly resemble the ones seen in X-rays in SN 1006; an origin in the explosion itself or from subsequent hydrodynamic instabilities both remain viable options. We have expanded the search for precursor X-ray emission ahead of a synchrotron-dominated shock front, as expected from diffusive shock acceleration theory, to numerous regions along both the northeast and southwest rims of the shell. Our data require that a precursor be thinner than about 3, and fainter than about 5% of the post-shock peak. These limits suggest that the magnetic field is amplified by a factor of seven or more in a narrow precursor region, promoting diffusive particle acceleration.

galaxies: active↗

On the Determination of Ejecta Structure and Explosion Asymmetry from the X-ray Knots of Cassiopeia A

We present a detailed analysis of Chandra X-ray spectra from individual ejecta knots in the supernova remnant Cassiopeia A. The spectra are fitted to give the electron temperature T(sub e), and (single) ionization age n(sub e)t. These quantities are compared with the predictions of self similar hydrodynamic models incorporating time dependent ionization and radiation losses, and Coulomb electron-ion equilibration behind the reverse shock, for a variety of different ejecta density profiles described by a uniform density core and a power law envelope. We find that the ejecta close to the 'jet' region in the NE, but not actually in the jet itself, have a systematically shallower outer envelope than ejecta elsewhere in the remnant, and we interpret this as being due to more energy of the initial explosion being directed in this polar direction as opposed to equatorially. The degree of asymmetry we infer is at the low end of that generally modelled in asymmetric core-collapse simulations, and may be used to rule out highly asymmetric explosion models.

Laming, J. Martin↗