Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock physics”

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 361 records · Page 20

Data-scarce surrogate modeling of shock-induced pore collapse process

Understanding the mechanisms of shock-induced pore collapse is of great interest in various disciplines in sciences and engineering, including materials science, biological sciences, and geophysics. However, numerical modeling of the complex pore collapse processes can be costly. To this end, a strong need exists to develop surrogate models for generating economic predictions of pore collapse processes. Here, in this work, we study the use of a data-driven reduced-order model, namely dynamic mode decomposition, and a deep generative model, namely conditional generative adversarial networks, to resemble the numerical simulations of the pore collapse process at representative training shock pressures. Since the simulations are expensive, the training data are scarce, which makes training an accurate surrogate model challenging. To overcome the difficulties posed by the complex physics phenomena, we make several crucial treatments to the plain original form of the methods to increase the capability of approximating and predicting the dynamics. In particular, physics information is used as indicators or conditional inputs to guide the prediction. In realizing these methods, the training of each dynamic mode composition model takes only around 30 s on CPU. In contrast, training a generative adversarial network model takes 8 h on GPU. Moreover, using dynamic mode decomposition, the final-time relative error is around 0.3% in the reproductive cases. We also demonstrate the predictive power of the methods at unseen testing shock pressures, where the error ranges from 1.3 to 5% in the interpolatory cases and 8 to 9% in extrapolatory cases.

97 MATHEMATICS AND COMPUTING↗

The Spectroscopic Evolution of the Symbiotic-like Recurrent Nova V407 Cygni During Its 2010 Outburst. 2. The Circumstellar Environment and the Aftermath

The nova outburst of V407 Cyg in 2010 Mar. 10 was the first observed for this star but its close resemblance to the well known symbiotic-like recurrent nova RS Oph suggests that it is also a member of this rare type of Galactic novae. The nova was the first detected at gamma-ray energies and is the first known nova explosion for this system. The extensive multiwavelength coverage of this outburst makes it an ideal comparison with the few other outbursts known for similar systems. We extend our previous analysis of the Mira and the expanding shock from the explosion to detail the time development of the photoionized Mira wind, circumstellar medium, and shocked circumstellar environment to derive their physical parameters and how they relate to large scale structure of the environment, extending the previous coverage to more than 500 days after outburst. We use optical spectra obtained at high resolution with the Nordic Optical Telescope (NOT) (R approx. =.45000 to 65000) and medium resolution Ondrejov Observatory (R approx. = 12000) data and compare the line variations with publicly available archival measurements at 30 GHz OVNR and at X-rays with Swift during the first four months of the outburst, through the end of the epoch of strong XR emission. We use nebular diagnostics and high resolution profile variations to derive the densities and locations of the extended emission. We find that the higher the ionization and/or the higher the excitation energy, the more closely the profiles resemble the He II/Ca V-type high velocity shock profile discussed in Paper I. This also accounts for the comparative development of the [N II] and [O III] isoelectronic transitions: the [O III] 4363A profile does not show the low velocity peaks while the excited [N II] 5754A does. If nitrogen is mainly N(+3) or higher in the shock, the upper state of the [N II] nebular lines will contribute but if the oxygen is O(+2) then this line is formed by recombination, masking the nebular contributor, and the lower states are collisionally quenched but emit from the low density surroundings. Absorption lines of Fe-peak ions formed in the Mira wind were visible as P Cyg profiles at low velocity before Day 69, around the time of the X-ray peak and we identified many absorption transitions without accompanying emission for metal lines. The H Balmer lines showed strong P Cyg absorption troughs that weakened during the 2010 observing period, through Day 128. The Fe-peak line profiles and flux variations were different for permitted and forbidden transitions: the E1 transitions were not visible after Day 128 but had shown a narrow peak superimposed on an extended (200 km/s) blue wing, while the M1 and E2 transitions persisted to Day 529, the last observation, and showed extended redshifted wings up of the same velocity. We distinguish the components from the shock, the photoionized environment, and the chromosphere and inner Mira wind using spectra taken more than one year after outburst. The multiple shells and radiative excitation phenomenology are similar to those recently cited for GRBs and SNIa .

Shore, S. N.↗

Collisionless ion-electron energy exchange in magnetized shocks

Energy partition between ions and electrons in collisionless shocks has been a long-standing unsolved fundamental physical question. Here, we show that kinetic simulations of moderate Alfv´enic Mach number, magnetized, collisionless shocks reveal rapid, faster-than-Coulomb, energy exchange between ions and electrons when the plasma is sufficiently magnetized. Using kinetic and multi-fluid models with counter-streaming ions, we identify resonances between electron whistler and ion magnetohydrodynamic waves that account for this rapid energy exchange.

High-energy-density plasmas↗

Numerical study of strong plasma shock waves produced in an electromagnetic shock tube.

A set of two-fluid Navier-Stokes equations with classical physical transport coefficients is used to compute the evolution and structure of collisional plasma shock waves in an electromagnetic shock tube. Shock speeds up to 200 cm/microsec and shocked-plasma temperatures of the order of kilovolts are studied. A strong transverse bias magnetic field is employed, which significantly alters the size and shape of the shock profiles, when compared with the zero bias field case. The wave structure is different from that of a two-fluid gasdynamic shock without transverse magnetic field. Especially significant is the effect of the small ion Larmor radius in reducing the gasdynamic shock thickness by at least an order of magnitude in the transverse shock case, which permits collisional shocks to have thicknesses much smaller than the post-shock mean free path. These collisional shock waves produce very hot ions, the ion temperature increasing with shock speed, although the formation time and distances also increase substantially with shock speed.

Schneider, S. H.↗

RaDIATE Collaboration Thermal shock studies

As next-generation accelerator target facilities (High Energy Physics, Spallation Sources, ... ) become increasingly more powerful and intense, high power target systems face key technical challenges, such as radiation damage and thermal shock. Those ultimately degrade the performance and lifetime of targets and have been identified as the leading cross-cutting challenges of high-power target facilities. In order to operate reliable beam-intercepting devices in the framework of energy and intensity increase for next generation accelerators, the RaDIATE Collaboration (Radiation Damage In Accelerator Target Environment), established in 2012 and managed by Fermilab, brings together existing expertise in nuclear material and accelerator targets from 20 international institutions, including CERN, to execute a coordinated strategy for high power targetry R&D. In this context, several thermal shock studies were performed at CERN's HiRadMat (High-Radiation to Materials) facility, that took a key step towards improving our knowledge on target damage tolerance. The first experiment HRMT-24, supported by EURCARD2 and completed in 2015, successfully validated the Johnson-Cook strength model developed at SwRI on beryllium S200FH (used for beam window material), providing a better confidence in simulating the thermal shock response of current and future S200FH beryllium components. HRMT-43, supported by ARIES and completed in 2018, tested various materials (Be, C, SiC, Si, Ti and ceramic nanofiber). It was a first and unique test which included pre-irradiated specimens from high energy proton beam irradiation to identify thermal shock response differences between non-irradiated and previously irradiated materials. Real-time measurement of dynamic thermomechanical response of graphite helped to benchmark numerical simulations.

43 PARTICLE ACCELERATORS↗

An investigation of bleed configurations and their effect on shock wave/boundary layer interactions

The design of high efficiency supersonic inlets is a complex task involving the optimization of a number of performance parameters such as pressure recovery, spillage, drag, and exit distortion profile, over the flight Mach number range. Computational techniques must be capable of accurately simulating the physics of shock/boundary layer interactions, secondary corner flows, flow separation, and bleed if they are to be useful in the design. In particular, bleed and flow separation, play an important role in inlet unstart, and the associated pressure oscillations. Numerical simulations were conducted to investigate some of the basic physical phenomena associated with bleed in oblique shock wave boundary layer interactions that affect the inlet performance.

Hamed, Awatef↗

Optimizing time integration for accurate recovery of shockwave interface location in radiography

We present simulations and experiments of time integrated radiographic imaging of a moving 1D shock wave front and a quantitative method for determining the statistical error in locating the shock front as a function of integration time and noise in the radiograph. We discuss the trade-off between increasing motion blur, which leads to decreased shock front location certainty, and increasing signal-to-noise, which leads to improved image quality with increasing integration time. We find an optimum integration time between a short integration time, where noise limits the error, and a long integration time, where motion blurring limits the error. This methodology can be used to tune experimental configurations to obtain the highest quality radiograph for a given experimental configuration.

Bremsstrahlung↗

Formation of collisionless shocks driven by strongly magnetized relativistic electrons in the laboratory

In experiments performed with the OMEGA EP laser system, proton deflectometry captured magnetic field dynamics consistent with collisionless shock formation driven by strongly magnetized relativistic electrons. During laser-foil interactions, shocks can form as relativistic electrons and strong surface magnetic fields generated by a short-pulse laser impinge on a cooler plasma produced by a longer-pulse laser. Three-dimensional particle-in-cell simulations reproduce the magnetic draping and fast formation speeds measured in the experiment and reveal that this relativistic-electron-driven shock forms at an interface that is unstable to shear and streaming instabilities. The simulation results provide insight into the microphysics that may influence high-energy shocks observed in extreme astrophysical environments. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Nonlinear simulation of under-resolved flows with shocks

Here, we consider the numerical simulation of advection-dominated flows whose wide range of physical length scales exceed the memory capacity of finite computers. Simulating flows with shocks and turbulence presented challenges for the earliest computers that were quickly overcome by the development of new numerical methodology. Principal among those new ideas were artificial viscosity and finite volume methods, concepts that remain in common use today. We begin by describing the history of those methods, the innovators and their motivations. We then describe the development of finite scale theory, a reformulation of Navier–Stokes theory that exposes the physical principles on which artificial viscosity is based. We discuss the essential properties of the finite scale equations, the observer, unresolved kinetic energy and inviscid energy dissipation. We briefly consider the implementation of the finite scale equations on the computer from the point of view of Gisin’s conjectures about finite information.

97 MATHEMATICS AND COMPUTING↗

Narrow Radiative Recombination Continua: A Signature of Ions Crossing the Contact Discontinuity of Astrophysical Shocks

X-rays from planetary nebulae (PNs) are believed to originate from a shock driven into the fast stellar wind (v ~ 1000 kilometers per second) as it collides with an earlier circumstellar slow wind (v ~ 10 kilometers per second). In theory, the shocked fast wind (hot hubble) and the ambient cold nebula can remain separated by magnetic fields along a surface referred to as the contact discontinuity (CD) that inhibits diffusion and heat conduction. The CD region is extremely difficult to probe directly owing to its small size and faint emission. This has largely left the study of CDs, stellar-shocks, and the associated micro-physics in the realm of theory. This paper presents spectroscopic evidence for ions from the hot bubble (kT approximately equal to 100 eV) crossing the CD and penetrating the cold nebular gas (kT approximately equal to 1 eV). Specifically, a narrow radiative recombination continuum (RRC) emission feature is identified in the high resolution X-ray spectrum of the PN BD+30degree3639 indicating bare C VII ions are recombining with cool electrons at kT(sub e) = 1.7 plus or minus 1.3 eV. An upper limit to the flux of the narrow RRC of H-like C VI is obtained as well. The RRCs are interpreted as due to C ions from the hot bubble of BD+30degree3639 crossing the CD into the cold nebula, where they ultimately recombine with its cool electrons. The RRC flux ratio of C VII to C VI constrains the temperature jump across the CD to deltakT greater than 80 eV, providing for the first time direct evidence for the stark temperature disparity between the two sides of an astrophysical CD, and constraining the role of magnetic fields and heat conduction accordingly. Two colliding-wind binaries are noted to have similar RRCs suggesting a temperature jump and CD crossing by ions may be common feature of stellar wind shocks.

Behar, Ehud↗

Refractive Imaging of Air Shock Above Microscale Defects in Pentaerythritol Tetranitrate (PETN) Films

Physical vapor deposition (PVD) of high explosives can produce energetic samples with unique microstructure and morphology compared to traditional powder processing techniques, but challenges may exist in fabricating explosive films without defects. Deposition conditions and substrate material may promote microcracking and other defects in the explosive films. Here, we investigate effects of engineered microscale defects (gaps) on detonation propagation and failure for pentaerythritol tetranitrate (PETN) films using ultra-high-speed refractive imaging and hydrocode modelling. Observations of the air shock above the gap reveal significant instabilities during gap crossing and re-ignition.

36 MATERIALS SCIENCE↗

Solar Sources of Earth-affecting Energetic Particles

Particle radiation from the Sun is one of the most important sources of hazardous space weather in the vicinity of Earth. Detailed studies of the origin of the so-called large solar energetic particle (SEP) events became possible only during the solar cycle 23, thanks to the availability of nearly continuous observation of the solar sources ofthese events. In particular, coronal mass ejections (CMEs), which are found to be a key requirement for the occurrence of an SEP event, have been recorded continuously only starting in the 1990s. The physical connection between CMEs and SEPs is that the CMEs drive a fast-mode MHD shock, which accelerates SEPs in the corona and interplanetary medium. The earliest indication of a shock is the occurrence of a type II radio burst at frequencies anywhere from more than a hundred MHz to a few MHz. Recent investigations using STEREO observations have revealed that the shock forms very close to the Sun - a mere 100,000 km above the surface. The shock formation depends not only on the CME properties, but also on the physical conditions in the ambient medium that supports shock propagation. This paper considers extreme cases of SEP events and the associated CMEs and type II radio bursts to illustrate the variability observed in SEP event properties. Comparison will be made between the events of solar cycles 23 and 24.

Gopalswamy, Nat↗

Fast imaging of laser ablation of cerium: Dependence of plasma plume dynamics and structure on native oxide growth

We have imaged the evolution of plasma following nanosecond pulsed laser ablation of cerium in air with sub-10 ns resolution. We find the ablated plasma plume expansion and shape to be strongly dependent on the time after mechanical polishing of the sample surface. After an initial shock front common to both new and aged samples (velocities up to 16 km s ₋1 ), we observe a second front that is relatively localized to the surface vicinity of newly polished samples, but in the case of older samples it moves outward at velocities of up to 4 km s ₋1 . We attribute this behavior to the formation and growth of a native oxide layer on the order of hours after polishing. These results demonstrate that plasma imaging can be used to diagnose the presence of, and evaluate the extent of, thin surface oxides or other compounds that form soon after exposure to some atmosphere. These findings also highlight the need to explicitly consider the presence of such overlayers when modeling laser/metal interaction, when making measurements of vapor phase chemistry, and when analyzing the composition, phase, and morphology of solid particulates formed after the ablation of reactive metals.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigating growth to detonation in vapor-deposited hexanitrostilbene and pentaerythritol tetranitrate films using high-throughput methods

In this work, a high-throughput experimental setup was used to characterize initiation threshold and growth to detonation in the explosives hexanitrostilbene (HNS) and pentaerythritol tetranitrate (PETN). The experiment sequentially launched an array of laser-driven flyers to shock samples arranged in a 96-well microplate geometry, with photonic Doppler velocimetry diagnostics to characterize flyer velocity and particle velocity at the explosive–substrate interface. Vapor-deposited films of HNS and PETN were used to provide numerous samples with various thicknesses, enabling characterization of the evolution of growth to detonation. One-dimensional hydrocode simulations were performed with reactions disabled to illustrate where the experimental data deviate from the predicted inert response. Prompt initiation was observed in 144 μm thick HNS films at flyer velocities near 3000 m/s and in 125 μm thick PETN films at flyer velocities near 2400 m/s. This experimental setup enables rapid quantification of the growth of reactions in explosive materials that can reach detonation at sub-millimeter length scales. These data can subsequently be used for parameterizing reactive burn models in hydrocode simulations, as discussed in Paper II [D. E. Kittell, R. Knepper, and A. S. Tappan, J. Appl. Phys. 131, 154902 (2022)].

36 MATERIALS SCIENCE↗

Chandra X-ray Observation of a Mature Cloud-Shock Interaction in the Bright Eastern Knot of Puppis A

We present Chandra X-ray images and spectra of the most prominent cloud-shock interaction region in the Puppis A supernova remnant. The Bright Eastern Knot (BEK) has two main morphological components: (1) a bright compact knot that lies directly behind the apex of an indentation in the eastern X-ray boundary and (2) lying 1 westward behind the shock, a curved vertical structure (bar) that is separated from a smaller bright cloud (cap) by faint diffuse emission. Based on hardness images and spectra, we identify the bar and cap as a single shocked interstellar cloud. Its morphology strongly resembles the "voided sphere" structures seen at late times in Klein et al. experimental simulat.ions of cloud-shock interactions, when the crushing of the cloud by shear instabilities is well underway. We infer an intera.ction time of roughly cloud-crushing timescales, which translates to 2000-4000 years, based on the X-ray temperature, physical size, and estimated expansion of the shocked cloud. This is the first X-ray identified example of a cloud-shock interaction in this advanced phase. Closer t o the shock front, the X-ray emission of the compact knot in the eastern part of the BEK region implies a recent interaction with relatively denser gas, some of which lies in front of the remnant. The complex spatial relationship of the X-ray emission of the compact knot to optical [O III] emission suggests that there are multiple cloud interactions occurring along the line of sight.

Hwang, Una↗

Aeroacoustic Flow Phenomena Accurately Captured by New Computational Fluid Dynamics Method

One of the challenges in the computational fluid dynamics area is the accurate calculation of aeroacoustic phenomena, especially in the presence of shock waves. One such phenomenon is "transonic resonance," where an unsteady shock wave at the throat of a convergent-divergent nozzle results in the emission of acoustic tones. The space-time Conservation-Element and Solution-Element (CE/SE) method developed at the NASA Glenn Research Center can faithfully capture the shock waves, their unsteady motion, and the generated acoustic tones. The CE/SE method is a revolutionary new approach to the numerical modeling of physical phenomena where features with steep gradients (e.g., shock waves, phase transition, etc.) must coexist with those having weaker variations. The CE/SE method does not require the complex interpolation procedures (that allow for the possibility of a shock between grid cells) used by many other methods to transfer information between grid cells. These interpolation procedures can add too much numerical dissipation to the solution process. Thus, while shocks are resolved, weaker waves, such as acoustic waves, are washed out.

Blech, Richard A.↗

High strain-rate strength response of single crystal tantalum through in-situ hole closure imaging experiments

The properties of crystalline materials often depend on directionality and operating conditions. Specifically, the strength of materials can depend anisotropically on crystal direction and the loading condition. To probe these effects, a preliminary series of high strain-rate (> 105/s) strength plate-impact hole closure experiments were performed on high purity single crystal Tantalum cubes. The orientation of the single crystals with respect to impact/loading were varied to provide data to inform crystal plasticity modeling efforts. The experiments consist of in-situ high-resolution X-ray radiographic imaging of the hole collapse under dynamic compression conditions to infer the material strength via its resistance to closure at increasing levels of plastic strain. The experiments are compared against hydrocode simulation predictions. Here, a comparison with simple elastic perfectly plastic strength model predictions is presented to elucidate the response of the different crystal orientations at high strain-rate and large plastic strains.

36 MATERIALS SCIENCE↗

The Effects of Climate Change on Interregional Electricity Market Dynamics on the U.S. West Coast

The United States (U.S.) West Coast power system is strongly influenced by variability and extremes in air temperatures (which drive electricity demand) and streamflows (which constrain hydropower production). As hydroclimate changes across the West Coast, a combination of forces may work in tandem to make its bulk power system more vulnerable to physical reliability issues and market price shocks. In particular, a warmer climate is expected to increase summer cooling (electricity) demands and shift the average timing of peak streamflow (hydropower production) away from summer to the spring and winter, depriving power systems of hydropower when it is needed the most. Here, we investigate how climate change could alter interregional electricity market dynamics on the West Coast, including the potential for hydroclimatic changes in one region (e.g. Pacific Northwest (PNW)) to “spill over” and cause price and reliability risks in another (e.g. California). We find that the most salient hydroclimatic risks for the PNW power system are changes in streamflow, while risks for the California system are driven primarily by changes in summer air temperatures, especially extreme heat events that increase peak system demand. Altered timing and amounts of hydropower production in the PNW do alter summer power deliveries into California but show relatively modest potential to impact prices and reliability there. Instead, it is future climate conditions (extreme heat) in California that may exert a stronger significant influence on prices and reliability in the PNW, especially if California continues to rely on its northern neighbor for imported power to meet higher summer demands.

24 POWER TRANSMISSION AND DISTRIBUTION↗