Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock initiation”

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 739 records · Page 41

A Wind-driven Origin for the Firework Morphology of the Supernova Remnant Pa 30

Pa 30—the likely remnant of the Galactic Type Iax supernova of 1181 AD—displays an unusual, firework-like morphology, consisting of radial filaments extending from a common center, where a white dwarf (WD) currently drives a very fast wind (speed ≳10 4 km s −1 ). We propose the filaments arose from the Rayleigh–Taylor-unstable nature of the interface between the circumstellar medium (CSM) and the shocked wind launched by the natal WD; the filaments then elongated intact due to the Kelvin–Helmholtz-stable nature of the large initial density contrast between the wind and CSM, supplemented by the slowly declining wind density profile (relative to homologously expanding ejecta). To support this interpretation, we present two-dimensional hydrodynamical simulations and derive the filament properties, including their speed, density, and temperature, all of which are consistent with observations. We suggest the filaments elongate until the wind and CSM densities become comparable at the contact discontinuity, which occurs within 1–10 yr, and then truncate because the Rayleigh–Taylor instability (RTI) halts. The subsequent Kelvin–Helmholtz instability growth timescale across the current width of the filaments is longer than the age of Pa 30, so they remain intact. The filament-less central region in Pa 30 is therefore more likely a consequence of the finite timescale over which the RTI operates, rather than a wind termination shock. In general, firework-like filaments may form in other systems, provided there is a sufficiently large density contrast between the ejecta and its surroundings.

79 ASTRONOMY AND ASTROPHYSICS↗

Measurement and Computation of Supersonic Flow in a Lobed Diffuser-Mixer for Trapped Vortex Combustors

The trapped vortex combustor (TVC) pioneered by Air Force Research Laboratories (AFRL) is under consideration as an alternative to conventional gas turbine combustors. The TVC has demonstrated excellent operational characteristics such as high combustion efficiency, low NO(x) emissions, effective flame stabilization, excellent high-altitude relight capability, and operation in the lean-burn or rich burn-quick quench-lean burn (RQL) modes of combustion. It also has excellent potential for lowering the engine combustor weight. This performance at low to moderate combustor mach numbers has stimulated interest in its ability to operate at higher combustion mach number, and for aerospace, this implies potentially higher flight mach numbers. To this end, a lobed diffuser-mixer that enhances the fuel-air mixing in the TVC combustor core was designed and evaluated, with special attention paid to the potential shock system entering the combustor core. For the present investigation, the lobed diffuser-mixer combustor rig is in a full annular configuration featuring sixfold symmetry among the lobes, symmetry within each lobe, and plain parallel, symmetric incident flow. During hardware cold-flow testing, significant discrepancies were found between computed and measured values for the pitot-probe-averaged static pressure profiles at the lobe exit plane. Computational fluid dynamics (CFD) simulations were initiated to determine whether the static pressure probe was causing high local flow-field disturbances in the supersonic flow exiting the diffuser-mixer and whether shock wave impingement on the pitot probe tip, pressure ports, or surface was the cause of the discrepancies. Simulations were performed with and without the pitot probe present in the modeling. A comparison of static pressure profiles without the probe showed that static pressure was off by nearly a factor of 2 over much of the radial profile, even when taking into account potential axial displacement of the probe by up to 0.25 in. (0.64 cm). Including the pitot probe in the CFD modeling and data interpretation lead to good agreement between measurement and prediction. Graphical inspection of the results showed that the shock waves impinging on the probe surface were highly nonuniform, with static pressure varying circumferentially among the pressure ports by over 10 percent in some cases. As part of the measurement methodology, such measurements should be routinely supplemented with CFD analyses that include the pitot probe as part of the flow-path geometry.

Brankovic, Andreja↗

Physics consistent machine learning framework for inverse modeling with applications to ICF capsule implosions

In high energy density physics (HEDP) and inertial confinement fusion (ICF), predictive modeling is complicated by uncertainty in parameters that characterize various aspects of the modeled system, such as those characterizing material properties, equation of state (EOS), opacities, and initial conditions. Typically, however, these parameters are not directly observable. What is observed instead is a time sequence of radiographic projections using X-rays. In this work, we define a set of sparse hydrodynamic features derived from the outgoing shock profile and outer material edge, which can be obtained from radiographic measurements, to directly infer such parameters. Our machine learning (ML)-based methodology involves a pipeline of two architectures, a radiograph-to-features network (R2FNet) and a features-to-parameters network (F2PNet), that are trained independently and later combined to approximate a posterior distribution for the parameters from radiographs. We show that the machine learning architectures are able to accurately infer initial conditions and EOS parameters, and that the estimated parameters can be used in a hydrodynamics code to obtain density fields, shocks, and material interfaces that satisfy thermodynamic and hydrodynamic consistency. Finally, we demonstrate that features resulting from an unknown EOS model can be successfully mapped onto parameters of a chosen analytical EOS model, implying that network predictions are learning physics, with a degree of invariance to the underlying choice of EOS model. To the best of our knowledge, our framework is the first demonstration of recovering both thermodynamic and hydrodynamic consistent density fields from noisy radiographs.

97 MATHEMATICS AND COMPUTING↗

Synopsis of SURF and SURFplus

Developed by Menikoff and Shaw, the SURF reactive burn model builds on the Ignition and Growth concept by incorporating the lead shock pressure directly into the volumetric hot-spot burn rate. The plus-extension augments the model with a late time surface burn rate due to carbon clustering. Together, SURF and SURFplus allow for robust modeling of both conventional and insensitive high explosives. Practically, the SURF burn models require shock detection and the advection of the lead shock pressure as an additional material field. Calibrated parameters for given equations of state and thermodynamic closure are dependent on initial temperature and density. An implementation in Python is given.

42 ENGINEERING↗

Stability and convergence of nuclear detonations in white dwarf collisions

We investigate the numerical stability of thermonuclear detonations in 1D accelerated reactive shocks and 2D binary collisions of equal-mass magnetized and unmagnetized white dwarf stars. To achieve high resolution at initiation sites, we devised geometric gridding and mesh velocity strategies specially adapted to the unique requirements of head-on collisional geometries, scenarios in which one expects maximum production of iron-group products. We study the effects of grid resolution and the limiting of temperature, energy, and reactants for different stellar masses, separations, magnetic fields, inigenerationtial compositions, detonation mechanisms, and limiter parameters across a range of cell sizes from 1 to 100 km. Our results set bounds on the parameter space of limiter amplitudes for which both temperature- and energy-limiting procedures yield consistent and monotonically convergent solutions. Within these bounds, we find that grid resolutions of 5 km or better are necessary for uncertainties in total released energy and iron-group products to drop below 10%. Intermediate-mass products (e.g., calcium) exhibit similar convergence trends but with somewhat greater uncertainty. These conclusions apply equally to pure C/O white dwarfs, multispecies compositions (including helium shells), magnetized and unmagnetized cores, and either single or multiple detonation scenarios.

79 ASTRONOMY AND ASTROPHYSICS↗

Anorthite - Thermal equation of state to high pressures

New shock-wave data are presented for anorthite from which a full high-temperature, high-pressure equation of state is derived. Whereas anorthite has relatively low values of thermal expansion and Grueneisen parameter at zero pressure, it is found that these attain relatively high values in the high density state corresponding to the high-pressure phase Hugoniot but decrease upon compression as expected. It is noted that higher order anharmonic contributions decrease more rapidly with pressure and that the thermal expansion therefore saturates to a high temperature value at pressures above about 100 GPa. Reduction of the Hugoniot data permits shock temperatures to be calculated; it also yields a principal adiabat for the high pressure branch of the Hugoniot. The initial bulk modulus of this adiabat is essentially identical to that of anorthite, whereas the initial density is about 3.40 Mg/cu m.

Jeanloz, R.↗

Dynamic compression and volatile release of carbonates

Particle velocity profiles upon shock compression and isentropic releases were measured for polycrystalline calcite. The Solenhofen limestone release paths lie, close to the Hugoniot. Calcite 3 to 2 transition, upon release, was observed, but rarefaction shocks were not detected. The equation of state is used to predict the fraction of material devolatilized upon isentropic release as a function of shock pressure. The effect of ambient partial pressure of CO2 on the calculations is demonstrated and considered in models of atmospheric evolution by impact induced mineral devolatilization. The radiative characteristics of shocked calcite indicate that localization of thermal energy occurs under shock compression. Shock entropy calculations result in a minimum estimate of 90% devolatilization upon complete release from 10 GPa. Isentropic release paths from calculated continuum Hugoniot temperatures cross into the CaO (solid) + CO2 (vapor) field at improbably low pressures. It is found that release paths from measured shock temperatures cross into the melt plus vapor field at pressures greater than .5 GPa, which suggests that devolatilization is initiated at the shear banding sites.

Tyburczy, J. A.↗

Dynamic compression and volatile release of carbonates

Particle velocity profiles upon shock compression and isentropic releases were measured for polycrystalline calcite. The Solenhofen limestone release paths lie close to the Hugoniot. Calcite 3 to 2 transition, upon release was observed, but rarefaction shocks were not detected. The equation of state is used to predict the fraction of material devolatilized upon isentropic release as a function of shock pressure. The effect of ambient partial pressure of CO2 on the calculations is demonstrated and considered in models of atmospheric evolution by impact induced mineral devolatilization. The radiative characteristics of shocked calcite indicate that localization of thermal energy occurs under shock compression. Shock entropy calculations result in a minimum estimate of 90 percent devolatilization upon complete release from 10 GPa. Isentropic release paths from calculated continuum Hugoniot temperatures cross into the CaO (solid) + CO2 (vapor) field at improbably low pressures. It is found that release paths from measured shock temperatures cross into the melt plus vapor field at pressures greater than .5 GPa, which suggests that devolatilization is initiated at the shear banding sites.

Tyburczy, J. A.↗

Calculating Sonic-Boom Propagation

Nonlinear effects included, enabling more-realistic modeling. Modified Method of Characteristics Sonic Boom Extrapolation Program (MMOC) is computer program for sonic-boom propagation that includes shock coalescence and incorporates effects of asymmetries due to volume and lift. Numerically integrates nonlinear governing equations using data on initial data line approximately one body length from aircraft and yields sonic-boom pressure at ground as function of time or of position at given time. MMOC written in FORTRAN IV.

Darden, C. M.↗

Pick-up ions at Comet P/Halley's bow shock - Observations with the IIS spectrometer on Giotto

Gaseous material expanding from the nucleus of Comet Halley into space forms the neutral coma around the comet. Ionization in the solar UV radiation removes particles from the coma and injects them into the solar wind plasma. These freshly created ions are accelerated by the interplanetary electric field on cycloidal trajectories with gyrocenters moving with the speed of the magnetic field lines. In the solar wind frame of reference these particles move along the magnetic field lines with a fixed pitch angle. Pitch-angle scattering and energy diffusion reduce quickly the initial energy anisotropy which is associated with the narrow pick-up structures. First observations of heavy cometary pick-up ions (water group ions) at the bow shock are presented. The evolution of the distribution function in the vicinity of the shock and radial density profiles are discussed.

Wilken, B.↗

Impact craters on Venus - Initial analysis from Magellan

The general features of impact craters are described emphasizing two aspects: the effect of the atmosphere on crater and ejecta morphology and the implications of the distribution and appearance of the craters for the volcanic and tectonic resurfacing history of Venus. Magellan radar images reveal 135 craters about 15 km in diameter containing central peaks, multiple central peaks, and peak rings. Craters smaller than 15 km exhibit multiple floors or appear in clusters. Surface flows of material initially entrained in the atmosphere are characterized. Zones of low radar albedo originated from deformation of the surface by the shock or pressure wave associated with the incoming meteoroid surround many craters. A spectrum of surface ages on Venus ranging from 0 to 800 million years indicates that Venus must be a geologically active planet.

Phillips, Roger J.↗

The making of an Alfvenic fluctuation: The resolution of a second-order analysis

Ulysses observations of the high speed polar streams show that they are largely occupied by very large amplitude Alfvenic fluctuations accompanied by many rotational discontinuities. These fluctuations have a nearly constant magnetic intensity or amplitude, and the magnetic field direction per wave cycle sweeps only through a limited arc, much as a car wiperblade would do. Barnes and Hollweg (JGR, 79, 2302, 1974) suggested that this unusual waveform could arise from an obliquely propagating and linearly polarized Alfven wave of finite amplitude. From a second-order analysis, they showed that the existence of a particular solution with a constant amplitude but could not resolve the outcome of the homogeneous solution which consisted of fast waves. They suggested that Landau damping of these fast waves may be needed to get the observed waveform. We present a 1 1/2 D hybrid simulation which is fully nonlinear and correctly describes the ion kinetics for an initially monochromatic and linearly polarized Alfven wave propagating obliquely to the background magnetic field. The wave has a large amplitude and a wavelength so long that it can be considered dispersionless for simulation times. At early times, the second harmonic in density and in magnetic field transverse to the initial wave magnetic field are generated and have more power than other harmonics. Steepening is observed with a weak fast shock emerging, but no rotational discontinuity is left behind, and instead a constant amplitude and an arc-shaped waveform is made. The compressional component which develops after the shocks have dissipated is to zeroth order better described as a pure acoustic wave than as a fast wave. This might be explained by the relaxing of the Alfven wave to a state where its ponderomotive force vanishes so that the compressional component can travel almost independently of it.

Vasquez, Bernard J.↗

Collisionless shock structures of Earth and other planets

This report summarizes the closing segment of our multi-spacecraft, multi-instrument study of collisionless shock structure. In this last year of our study, we have necessarily concentrated on subjects that limited time and remaining resources could be expected to bring to reasonable stopping points, if not full conclusions. Our attention has been focused therefore on matters that were either well underway when the year began or that could be expected to yield rapidly completed reports publishable quickly in abbreviated versions. Contemporary publication delays prevent any new initiatives from reaching the literature within the year in the best of circumstances. The topics that fell into these categories were detailed plasma wave (pw) phenomenology in slow shocks in the Earth's distant geomagnetic tail, instantaneous orientations of theta(sub Bn) in quasiparallel (Q(sub parallel)) shock structure, and a comprehensive overview of the relationship between structural ULF waves in the Qll shock environment and waves in the magnetosphere, i.e. geomagnetic ULF pulsations. The remainder of this report describes our freshly completed results, discusses two related investigations of pw waves in the foreshock and magnetosheath, and appends the abstracts of published papers and the texts of papers in press.

Greenstadt, Eugene W.↗

Silicate Melting and Vaporization During Rocky Planet Formation

Collisions that induce melting and vaporization can have a substantial effect on the thermal and geochemical evolution of planets. However, the thermodynamics of major minerals are not well known at the extreme conditions attained during planet formation.We obtained new data at the Sandia Z Machine and use published thermodynamic data for the major mineral forsterite (Mg 2 SiO 4 ) to calculate the specific entropy in the liquid region of the principal Hugoniot.We use our calculated specific entropy of shocked forsterite, and revised entropies for shocked silica, to determine the critical impact velocities for melting or vaporization upon decompression from the shocked state to 1 bar and the triple points, which are near the pressures of the solar nebula.We also demonstrate the importance of the initial temperature on the criteria for vaporization. Applying these results to N-body simulations of terrestrial planet formation, we find that up to 20% to 40% of the total system mass is processed through collisions with velocities that exceed the criteria for incipient vaporization at the triple point. Vaporizing collisions between small bodies are an important component of terrestrial planet formation.

58 GEOSCIENCES↗

Comet Shoemaker-Levy 9: Impact on Jupiter and plume evolution

The impact of fragments of Comet Shoemaker-Levy 9 on Jupiter and the resulting vapor plume expansion are investigated by conducting three-dimensional numerical simulations using the smoothed particle hydrodynamics (SPH) method. An icy body, representing the cometary fragments, with a velocity of 60 km/sec and a diameter of 2 km can penetrate to 350 km below the 1-bar pressure level in the atmosphere. Most of the initial kinetic energy of the fragment is transferred to the atmosphere between 50 km and 300 km below the 1-bar pressure level. The shock-heated atmospheric gas in the wake is totally dissociated and partially ionized. Scaling our SPH results to other sizes indicates that fragments larger than approximately 100 m in diameter can penetrate to below the visible cloud decks. The energy deposited in the atmosphere is explosively released in the upward expansion of the resulting plume. The plume preferentially expands upward rather than horizontally due to the density gradient of the ambient atmosphere. It rises greater than or equal to 10(exp 2) km in approximately 10(exp 2) sec. Eventually the total atmospheric mass ejected to above 1 bar is greater than or equal to 40 times the initial mass of the impactor. The plume temperature at a radius approximately 10(exp 3) km is greater than 10(exp 3) K for 10(exp 3) sec for a 2-km fragment. We predict that impact-induced plumes will be observable with the remote sensing instruments of the Galileo spacecraft. As the impact site rotates into the view of Earth some 20 min after the impact, the plume expansion will be observable using the Hubble Space Telescope (HST) and from visible and infrared instruments on groundbased telescopes. The rising plume reaches approximately 3000 km altitude in approximately 10 min and will be visible from Earth.

Takata, Toshiko↗

The effect of jet–ejecta interaction on the viewing angle dependence of kilonova light curves

ABSTRACT The merger of two neutron stars produces an outflow of radioactive heavy nuclei. Within a second of merger, the central remnant is expected to also launch a relativistic jet, which shock-heats and disrupts a portion of the radioactive ejecta. Within a few hours, emission from the radioactive material gives rise to an ultraviolet, optical, and infrared transient (a kilonova). We use the endstates of a suite of 2D relativistic hydrodynamic simulations of jet–ejecta interaction as initial conditions for multidimensional Monte Carlo radiation transport simulations of the resulting viewing angle-dependent light curves and spectra starting at $1.5\, \mathrm{h}$ after merger. We find that on this time-scale, jet shock heating does not affect the kilonova emission for the jet parameters we survey. However, the jet disruption to the density structure of the ejecta does change the light curves. The jet carves a channel into the otherwise spheroidal ejecta, revealing the hot, inner regions. As seen from near (≲30°) the jet axis, the kilonova is brighter by a factor of a few and bluer. The strength of this effect depends on the jet parameters, since the light curves of more heavily disrupted ejecta are more strongly affected. The light curves and spectra are also more heavily modified in the ultraviolet than in the optical.

Klion, Hannah↗

Initial Response of the Aurora to the January 10, 1997 Magnetic Cloud

On January 10th, 1997, a magnetic cloud originating at the Sun was incident on the Earth. The initial disturbance to the magnetosphere, as reflected in the intensification of the aurora, was measured by the Ultraviolet Imager on the Polar Spacecraft. The first activation of the aurora at local noon occurred within minutes of the arrival of the shock. The subsequent evolution of the aurora over the next 18 minutes shows that the magnetic disturbance proceeds from local noon, symmetrically around the dawn and dusk flanks to local midnight. The substorm onset was observed to occur 174 minutes after the initial brightening of the aurora and 78 minutes after the southward turning of the IMF (Interplanetary Magnetic Field). During the intervening time, significant polar cap precipitation is observed. The polar cap precipitation begins at the poleward edge of the oval in the post midnight region and develops to form several complex transpolar structures. The polar cap precipitation subsides and quiet conditions are observed for 40 minutes prior to the onset of the substorm. During this event we have observed several unusual unique auroral forms develop that are different from the standard substorm models. We will present interpretation of the development of the pre-substorm events in light of the interplanetary conditions.

Spann, James F., Jr.↗

Impact and collisional processes in the solar system

As impact cratered terrains have been successively recognized on certain planets and planetary satellites, it has become clear that impact processes are important to the understanding of the accretion and evolution of all solid planets. The noble gases in the normalized atmospheric inventories of the planets and the normalized gas content of meteorites are grossly similar, but demonstrate differences from each other which are not understood. In order to study shock devolatilization of the candidate carrier phases which are principally thought to be carbonaceous or hydrocarbons in planetesimals, experiments were conducted on noble gase implantation in various carbons: carbon black, activated charcoal, graphite, and carbon glass. These were candidate starting materials for impact devolatilization experiments. Initial experiments were conducted on vitreous amorphous carbon samples which were synthesized under vapor saturated conditions using argon as the pressurizing medium. An amino acid and surface analysis by laser ionization analyses were performed on three samples of shocked Murchison meteorite. A first study was completed in which a series of shock loading experiments on a porous limestone and on a non-porous gabbro in one and three dimensions were performed. Also a series of recovery experiments were conducted in which shocked molten basalt a 1700 C is encapsulated in molybdenum containers and shock recovered from up to 6 GPa pressures.

Ahrens, Thomas J.↗