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 523 records · Page 29

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↗

Cold Fronts in Clusters of Galaxies: Observations and Modeling

Mergers of galaxy clusters -- some of the most energetic events in the Universe -- produce disturbances in hot intracluster medium, such as shocks and cold fronts, that can be used as tools to study the physics of galaxy clusters. Cold fronts may constrain viscosity and the structure and strength of the cluster magnetic fields. Combined with radio data, these observations also shed light on the production of ultrarelativistic particles that are known to coexist with the cluster thermal plasma. This talk will summarize the current X-ray observations of cluster mergers, as well as some recent radio data and high resolution hydrodynamic simulations.

Markevitch, Maxim↗

Shock compression of cemented tungsten carbides to 100 GPa: Structure of shock waves, Hugoniot relations, and strength

Plate impact experiments are conducted on cemented tungsten carbides (WC) with a 3.7 and 6.0 wt. % cobalt binder to better understand their dynamic, high-pressure response to 100 GPa. The measured wave profiles show propagation of steady structured waves. Standard impedance matching procedures are used to determine the Hugoniot relations in the shock velocity–particle velocity (U s –v p ) and Hugoniot stress–specific volume (P–V/V o ) planes. The Hugoniot elastic limit of the samples is controlled by ductility of the Co binder and is determined to be 4.45 ± 0.29 GPa for cemented WC with 3.7 wt. % cobalt and 3.72 ± 0.24 GPa for a 6.0 wt. % cobalt binder. Both grades show a non-linear U s –v p relationship depending on whether the particle velocity is in the strength dominated or hydrodynamic regime. In the strength dominated regime, a non-linear decrease in U s is observed as v p increases from ambient to the material’s hydrodynamic limit. In the hydrodynamic regime, the U s –v p Hugoniot is linear and is determined to be U s = 4.97(±0.006)+1.446(±0.018)v p km/s for WC with 3.7 wt. % Co and U s = 4.93(±0.006)+1.463(±0.017)v p km/s for 6 wt. % Co. Both WC grades indicate shear-stress hardening with mean stress immediately after yield, followed by pressure softening, and then a sharp fall in stress carrying capacity as the mean stress is increased to ≈70 GPa (hydrodynamic limit) and beyond. This behavior is in contrast to pure WC ceramics, which show continued shear-stress hardening with mean stress to ≈80 GPa.

36 MATERIALS SCIENCE↗

Collisions and close encounters between massive main-sequence stars

Smooth-particle hydrodynamics is used to investigate collisions and close encounters between two massive main-sequence stars. Approximate scaling relations and fitting formulas are given for the amount of mass loss and for the critical impact parameters for capture or merging. The more massive stars, which have smaller ratios of specific binding energy to the square of escape velocity, are more easily disrupted in collisions. As the relative velocity increases, the orbital energy of the colliding stars can only be dissipated by shock heating, and the critical capture radius decreases much faster than predicted by linear theory. The average fractional mass loss per collision in a cluster is found not to depend sensitively on the stellar velocity dispersion.

Lai, Dong↗

Pinwheel Experiment Confirmatory Shot 1 Pre-Shot Report

The purpose of this experimental series is to validate the reactive burn model in the PAGOSA hydrodynamic code. It will also help validate the EOS and Strength models for the aluminum and tantalum cylinders. C4 high explosive (HE) is used to drive an aluminum cylinder into an aluminum sandwiched PBX-9501. The follow-on experiment will drive the aluminum cylinder into a tantalum sandwiched PBX-9501. The experiment will provide time of arrival (TOA) data for the burn front in the HE through the metal cylinders using Fiber Bragg Gratings (FBG). Photon Doppler Velocimetry (PDV) will provide shock wave TOA and the velocity of the inward moving cylinder. The PAGOSA code will be used to predict TOA and cylinder wall velocity using PDV and tracers.

42 ENGINEERING↗

The high-latitude cloud MBM 7. I. H I and CO observations

The high-latitude cloud (HLC) MBM 7 has been observed in the 21 cm H I line and the 12CO(1-0) and 13CO(1-0) lines with similar spatial resolutions. The data reveal a total mass approximately 30 M solar for MBM 7 and a complex morphology. The cloud consists of a cold dense core of 5 M solar surrounded by atomic and molecular gas with about 25 M solar, which is embedded in hotter and more diffuse H I gas. We derive a total column density N(H I + 2H2) of 1 x 10(21) cm-2 toward the center and 1 x 10(20) cm-3 toward the envelope of MBM 7. The CO line indicates the existence of dense cores [n(H2) > or = 2000 cm-3] of size (FWHM) approximately 0.5 pc. The morphology suggests shock compression from the southwest direction, which can form molecular cores along the direction perpendicular to the H I distribution. The H I cloud extends to the northeast, and the velocity gradient appears to be about 2.8 km s-1 pc-1 in this direction, which indicates a systematic outward motion which will disrupt the cloud in approximately 10(6) yr. The observed large line widths of approximately 2 km s-1 for CO suggest that turbulent motions exist in the cloud, and hydrodynamical turbulence may dominate the line broadening. Considering the energy and pressure of MBM 7, the dense cores appear not to be bound by gravity, and the whole cloud including the dense cores seem to be expanding. The distance to HLCs suggest that they belong to the galactic plane, since the scale height of the cloud is < or approximately equal to 100 pc. Compared to the more familiar dense dark clouds, HLCs may differ only in their small mass and low density, with their proximity reducing the filling factor and enhancing the contrast of the core and envelope structure.

Non-NASA Center↗

Algorithmic Advancements for High-Order Self-Gravitating Hydrodynamics

Self-gravity plays a key role in the formation and evolution of many astronomical objects. Though gravity is often dominant at large scales, other forces (e.g., gas pressure gradients, radiation, and/or magnetic fields) often compete. It is therefore essential for numerical simulations to evaluate their interplay accurately and robustly. Hanawa & Mullen derived a 4th-order accurate finite volume scheme to solve the equations of self-gravitating hydrodynamics on a uniform Cartesian grid. In this work, we supply improvements to the algorithm that (1) mitigate spurious gravitational circulation and (2) greatly simplify the evaluation of the high order corrections. The proposed algorithm provides the gravitational acceleration (ρg) and the gravitational energy release (ρv · g) as source terms for the hydrodynamic equations, all while preserving conservation of linear momentum. Spurious heating and/or cooling associated with truncation error in the numerical evaluation of the gravitational energy release decreases in proportion to the fourth power of the cell width. We demonstrate fourth order convergence on smooth problems (e.g., 3D inclined sound wave propagation and 3D equilibria). An application test tracks the spherical collapse of a polytrope by an imposed, sudden decrease of the central gas pressure; a bounce and second collapse (associated with a spherical accretion shock) are robustly captured by the high order algorithm.

79 ASTRONOMY AND ASTROPHYSICS↗

The role of high-speed plasma flows in plasmaspheric refilling

A model of time-dependent one-stream interhemispheric plasma flow is used to investigate plasmaspheric refilling. In the model the coupled time-dependent hydrodynamic equations (continuity, momentum and energy) of a two-ion, H(+) and O(+), quasi-neutral, currentless plasma are solved for a closed geomagnetic field line. For the present set of calculations an L = 2 field line was used. A steady state solution was found and used as the initial condition in subsequent simulations of the effect of density depletions. Density depletions were modeled by reducing the densities by an arbitrary factor above 2500 km altitude, while keeping the velocities and temperatures unchanged. Shock structures develop which move up the field line, meeting and reflecting at the equator. The reflected shocks are absorbed by the dense neutral atmosphere when they reach the model flux tube boundaries. After the upwelling streams collide, refilling occurs from the equator downward for 2 to 3 hours; thereafter there is steady refilling from the ionospheres upward.

Guiter, S. M.↗

Dynamics and X-ray emission of a galactic superwind interacting with disk and halo gas

There is a general agreement that the conspicuous extranuclear X-ray, optical-line, and radio-contiuum emission of starbursts is associated with powerful galactic superwinds blowing from their centers. However, despite the significant advances in observational studies of superwinds, there is no consensus on the nature of the emitting material and even on the emission mechanisms themselves. This is to a great extent a consequence of a poor understanding of dynamical processes in the starburst superwind regions. To address this issue, we have conducted two-dimensional hydrodynamical simulations of galactic superwinds. While previous similar studies have used a single (disk) component to represent the ISM of the starburst galaxy, we analyze the interaction of the wind with a two-component disk-halo ambient interstellar medium and argue that this two-component representation is crucial for adequate modeling of starbursts. The emphasis of this study is on the geometry and structure of the wind region and the X-ray emission arising in the wind material and the shocked gas in the disk and the halo of the galaxy. The simulation results have shown that a clear-cut bipolar wind can easily develop under a range of very different conditions. On the other hand, a complex 'filamentary' structure associated with the entrained dense disk material is found to arise within the hot bubble blown out by the wind. The flow pattern within the bubble is dominated equally by the central biconic outflow and a system of whirling motions r elated to the origin and development of the 'filaments'. The filament parameters make them a good candidate for optical-emission-line filamentary gas observed in starburst halos. We find that the history of mass and energy deposition in the starburst region of the galaxy is crucial for wind dynamics. A 'mild' early wind, which arises as a result of the cumulative effect of stellar winds from massive stars, produces a bipolar vertical cavity in the disk and halo gas without strongly affecting the gaseous disk, thus creating conditions for virtually free vertical escape of the hot gas at the later, much more violent supernova-dominated phases of the starburst. We calculate the luminosity, mass, and effective temperature of the X-ray emitting gas in the 'soft' (0.1 to 0.7 keV, 0.7 to 2.2 keV, and 0.1 to 2.2 keV) and 'hard' (1.6 to 8.3 keV) energy bands and estimate the contribution of different gaseous components to the X-ray flux in these bands. Analysis of these parameters enables us to make conclusions regarding the nature of the X-ray-emitting material. We have inferred that the bulk of the soft thermal X-ray emission from starbursts arises in the wind-shocked material of the disk and halo gas rather than in the wind material itself. This enables us to predict that the integrated soft X-ray spectra of starbursts need not show an overabundance of heavy elements which are believed to be produced copiously in the centers of starbursts. Unlike soft X-ray emission, the hard component of thermal X-ray emission is found to originate in the wind material ejected from the starburst region. However, the derived ratio of hard-to-soft X-ray luminosities is too small compared to that observed in starbursts. We conclude therefore that the observed hard X-ray emission of starbursts is probably not associated with the thermal emission of hot wind or ambient shocked gas. Typical temperatures of the bulk of the soft X-ray-emitting material in our very different models have been found to agree well with the ones estimated on the basis of the ROSAT data for the soft component of X-ray emission of nearby starbursts. We predict that temperatures of the extranuclear soft X-ray-emitting gas in starburst galaxies with heavy element abundances near solar should be close to T(sub Xs = 2 to 5 x 10(exp 6)K.

Suchkov, Anatoly A.↗

Measurements of the imploding plasma sheath in triple-nozzle gas-puff z pinches

Gas-puff z-pinch implosions are characterized by the formation of a dense annular plasma shell, the sheath, that is driven to the axis by magnetic forces and therefore subject to the magneto-Rayleigh–Taylor instability. Here, the conditions within these sheaths are measured on the 1-MA COBRA generator at Cornell University for various gas species and initial fill densities. The gas-puff loads are initialized by a 7 cm diameter triple-nozzle gas valve assembly with concentric outer and inner annular nozzles and a central gas jet. Thomson scattering and laser interferometry provide spatially resolved flow, temperature, and electron density profiles midway through the implosion, while extreme ultraviolet pinhole cameras record the evolution of the plasma column and photoconducting diodes measure x-ray emission. Analysis of the scattering spectra includes a means of discriminating between thermal and non-thermal broadening to test for the presence of hydrodynamic turbulence. Two types of sheath profiles are observed, those with sharp discontinuities at the leading edge and those with smooth gradients. In both cases, non-thermal broadening is generally peaked at the front of the sheath and exhibits a characteristic decay length that roughly scales with the sheath ion mean free path. Here, we demonstrate that this non-thermal broadening term is inconsistent with laminar velocity gradients and is more consistent with dissipative turbulence driven by unstable plasma waves in a collisionless shock. The resulting differences in sheath profile are then set by the sheath ion collisionality in a manner consistent with recent 1D kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Design and modeling of indirectly driven magnetized implosions on the NIF

The use of magnetic fields to improve the performance of hohlraum-driven implosions on the National Ignition Facility (NIF) is discussed. The focus is on magnetically insulated inertial confinement fusion, where the primary field effect is to reduce electron-thermal and alpha-particle loss from the compressed hotspot (magnetic pressure is of secondary importance). We summarize the requirements to achieve this state. The design of recent NIF magnetized hohlraum experiments is presented. These are close to earlier shots in the three-shock, high-adiabat (BigFoot) campaign, subject to the constraints that magnetized NIF targets must be fielded at room-temperature, and use ≲1 MJ of laser energy to avoid the risk of optics damage from stimulated Brillouin scattering. We present results from the original magnetized hohlraum platform, as well as a later variant that gives a higher hotspot temperature. In both platforms, imposed fields (at the capsule center) of up to 28 T increase the fusion yield and hotspot temperature. Integrated radiation-magneto-hydrodynamic modeling with the Lasnex code of these shots is shown, where laser power multipliers and a saturation clamp on cross-beam energy transfer are developed to match the time of peak capsule emission and the P2 Legendre moment of the hotspot x-ray image. The resulting fusion yield and ion temperature agree decently with the measured relative effects of the field, although the absolute simulated yields are higher than the data by 2.0−2.7×. The tuned parameters and yield discrepancy are comparable for experiments with and without an imposed field, indicating the model adequately captures the field effects. Self-generated and imposed fields are added sequentially to simulations of one BigFoot NIF shot to understand how they alter target dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Nature of the Vela Supernova Remnant as Revealed by O VI and C IV Absorption

Highly ionized gas, in particular C IV and O VI, is produced in the interstellar medium in regions with hot (T approx. 10(exp 6) K) X-ray emitting gas and at the boundaries where hot gas and cooler (T approx. 10(exp 4) K) gas interact. Supernova remnant shocks produce most of the hot gas in the ISM and, if they are in the correct range of speeds, should produce observable quantities of C IV and O VI absorption. In turn, the column densities of these ions are potentially powerful diagnostics of the shock speed and interstellar environment in which the SNR is evolving. With the advent of FUSE, the power of this diagnostic technique is now available. We have FUSE data toward 8 stars behind the Vela SNR, and have developed a data reduction and analysis method that produces reasonably reliable O VI column densities, in spite of the complexities of the FUSE spectra in this region. In order to gain insight into the observational results, the Vela SNR evolution was modelled using Piecewise Parabolic Method numerical hydrodynamics code. The code is 1-D and incorporates non-equilibrium ionization, radiative cooling, thermal conduction and magnetic pressure.

Lines, Nichols J.↗

Benchmarking boron carbide equation of state using computation and experiment

Boron carbide ( B 4 C ) is of both fundamental scientific and practical interest due to its structural complexity and how it changes upon compression, as well as its many industrial uses and potential for use in inertial confinement fusion (ICF) and high-energy density physics experiments. Here, we report the results of a comprehensive computational study of the equation of state (EOS) of B 4 C in the liquid, warm dense matter, and plasma phases. Our calculations are cross-validated by comparisons with Hugoniot measurements up to 61 megabar from planar shock experiments performed at the National Ignition Facility (NIF). Our computational methods include path integral Monte Carlo, activity expansion, as well as all-electron Green's function Korringa-Kohn-Rostoker and molecular dynamics that are both based on density functional theory. We calculate the pressure-internal energy EOS of B 4 C over a broad range of temperatures ( ~ 6 × 10 3 – 5 × 10 8 K) and densities (0.025–50 g / cm 3 ). We assess that the largest discrepancies between theoretical predictions are ≲ 5 % near the compression maximum at 1– 2 × 10 6 K. This is the warm-dense state in which the K shell significantly ionizes and has posed grand challenges to theory and experiment. By comparing with different EOS models, we find a Purgatorio model (LEOS 2122) that agrees with our calculations. The maximum discrepancies in pressure between our first-principles predictions and LEOS 2122 are ~ 18 % and occur at temperatures between 6 × 10 3 – 2 × 10 5 K, which we believe originate from differences in the ion thermal term and the cold curve that are modeled in LEOS 2122 in comparison with our first-principles calculations. To account for potential differences in the ion thermal term, we have developed three new equation-of-state models that are consistent with theoretical calculations and experiment. We apply these new models to 1D hydrodynamic simulations of a polar direct-drive NIF implosion, demonstrating that these new models are now available for future ICF design studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The density distribution and physical origins of intermittency in supersonic, highly magnetized turbulence with diverse modes of driving

The probability density function (PDF) of the logarithmic density contrast, s = ln (ρ/ρ 0 ), with gas density ρ and mean density ρ 0 , for hydrodynamical supersonic turbulence is well known to have significant non-Gaussian (intermittent) features that monotonically increase with the turbulent Mach number, M. By studying the mass- and volume-weighted s-PDF for an ensemble of 36 sub-to-trans-Alfvenic mean-field, supersonic, isothermal turbulence simulations with different modes of driving, relevant to molecular gas in the cool interstellar medium, we show that a more intricate picture emerges for the non-Gaussian nature of s. Using four independent measures of the non-Gaussian components, we find hydrodynamical-like structure in the highly magnetized plasma for M≲4. However, for M≳4, the non-Gaussian signatures disappear, leaving approximately Gaussian s-statistics – exactly the opposite of hydrodynamical turbulence in the high-M limit. We also find that the non-Gaussian components of the PDF increase monotonically with more compressive driving modes. To understand the M≲4 non-Gaussian features, we use one-dimensional pencil beams to explore the dynamics along and across the large-scale magnetic field, B 0 . Here we discuss kinetic, density, and magnetic field fluctuations from the pencil beams, and identify physical sources of non-Gaussian components to the PDF as single, strong shocks coupled to fast magnetosonic compressions that form along B 0 . We discuss the Gaussianization of the M≳4 s-fields through the lens of two phenomenologies: the self-similarity of the s-field and homogenization of the dynamical time-scales between the over- and underdense regions in the compressible gas.

79 ASTRONOMY AND ASTROPHYSICS↗

The Modified Void Nucleation and Growth Model (MNAG) for Damage Evolution in BCC Ta

A void coalescence term was proposed as an addition to the original void nucleation and growth (NAG) model to accurately describe void evolution under dynamic loading. The new model, termed as modified void nucleation and growth model (MNAG model), incorporated analytic equations to explicitly account for the evolution of the void number density and the void volume fraction (damage) during void nucleation, growth, as well as the coalescence stage. The parameters in the MNAG model were fitted to molecular dynamics (MD) shock data for single-crystal and nanocrystalline Ta, and the corresponding nucleation, growth, and coalescence rates were extracted. The results suggested that void nucleation, growth, and coalescence rates were dependent on the orientation as well as grain size. Compared to other models, such as NAG, Cocks–Ashby, Tepla, and Tonks, which were only able to reproduce early or later stage damage evolution, the MNAG model was able to reproduce all stages associated with nucleation, growth, and coalescence. The MNAG model could provide the basis for hydrodynamic simulations to improve the fidelity of the damage nucleation and evolution in 3-D microstructures.

36 MATERIALS SCIENCE↗

Oxgen-burning hydrodynamics. 1: Steady shell burning

With new hydrodynamic techniques, the relatively fast evolutionary stages of a star prior to core collapse may be explicitly computed in two spatial dimensions, with a treatment of the microphysics (e.g., nuclear reactions, equation of state, neutrino cooling) which is comparable to typical one-dimensional simulations. The nature of shell oxygen burning in a massive star, prior to core collapse, is used as a first example; it is of particular interest because it is (1) the region in which Ni-56 will be produced by the supernova shock, (2) the region of the 'mass cut', which will separate the collapsed core from the ejected mantle, (3) the site of much of the explosive nucleosynthesis, and (4) a suggested source of symmetry breaking to drive mixing instabilities which were observed in SN 1987A. The nature of the shell burning affects the size of the core which will collapse. The method is illustrated on this test case, and the character of the convection is examined.

Arnett, David↗

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

Deciphering the nature of the pulsar wind nebula CTB 87 with XMM-Newton

CTB 87 (G74.9+1.2) is an evolved supernova remnant (SNR) which hosts a peculiar pulsar wind nebula (PWN). The X-ray peak is offset from that observed in radio and lies towards the edge of the radio nebula. The putative pulsar, CXOU J201609.2+371110, was first resolved with Chandra and is surrounded by a compact and a more extended X-ray nebula. Here, we use a deep XMM–Newton observation to examine the morphology and evolutionary stage of the PWN and to search for thermal emission expected from a supernova shell or reverse shock interaction with supernova ejecta. We do not find evidence of thermal X-ray emission from the SNR and place an upper limit on the electron density of 0.05 per cu.cm for a plasma temperature kT ∼ 0.8 keV. The morphology and spectral properties are consistent with a ∼20-kyr-old relic PWN expanding into a stellar wind-blown bubble. We also present the first X-ray spectral index map from the PWN and show that we can reproduce its morphology by means of 2D axisymmetric relativistic hydrodynamical simulations.

B Guest↗