Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “moving shock”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Symmetry Determining Equations of the Rankine-Hugoniot Equations for Variable Velocity Shock Waves

The “constant velocity piston” problem (Fig. 1), also known as the “piston problem,” is a standard model for a one dimensional, in our case linear, symmetric shock wave moving through an inviscid, perfect gas. The model can be divided into two regions - a perturbed section on the left and an unperturbed section on the right - by a moving shock wave moving left to right. Both the perturbed and unperturbed sections, i.e. the shocked and unshocked regions, respectively, obey the Eulerian conservation equations; however, at the exact location of the shock, there is a mathematical discontinuity not satisfied by the Euler equations. To ensure continuity and conservation of certain quantities when crossing between the unshocked and shocked regions, we evoke a series of equations derived from the Eulerian conservation equations, called the Rankine-Hugoniot equations, or “jump” equations as it is often referred to in the literature on the topic. The classical constant-velocity piston problem assumes the piston features a constant driving velocity (among many other willing suspensions of belief required in the pursuit of a first principles equation model); consequent to this assumption is a constant-velocity shock and a constant-velocity shocked flow state. However, using Lie Group Theory (LGT), also known as symmetry analysis, we can attempt to reinterpret the model with a shock wave of variable velocity in time and space. An extension of the model in this way opens up the possibility for obtaining new analytical solutions to the piston problem for certain shock velocity models. In this report, we use LGT to derive the symmetry determining equations (SDEs), whose solutions are Lie groups, which permit analytical solutions. In the future, we can then use the SDEs to define constraint equations on the shock velocity model and what the successive solutions to the Euler equations might be based off such constraints. This report is structured as follows: Section 2 provides a brief derivation of the Rankine-Hugoniot (“jump”) equations; Section 3 gives an overview of Lie group theory; Section 4 derives the SDEs of the jump equations; Section 5 derives the Euler conservation equations for fluids; and Section 6 presents concluding remarks and opportunities for future studies.

42 ENGINEERING↗

Examining astrophysical gas cloud collapse using an optical depth-scaled, x-ray-irradiated, carbon-foam sphere

When stellar radiation interacts with a molecular cloud, the cloud's fate depends on the strength of the incident radiation and the radiation's mean-free-path within the cloud [F. Bertoldi, Astrophys. J. 346, 735–755 (1989)]. Under the right conditions, the radiation compresses the cloud and a star formation may occur. Where and when the stellar formation occurs in the cloud's collapse are open questions. Direct observation of the complete star–cloud lifecycle is nearly impossible due to the immense timescales and distances over which the interaction occurs. Laboratory astrophysics offers a way to investigate such a system by scaling the important astrophysical parameters to the laboratory. This work describes laboratory experiments to study the radiation-driven implosion of clouds, using x rays from a laser-irradiated, thin, gold foil as a surrogate star and a carbon-foam sphere as a surrogate cloud. An optically thick system, theoretically corresponding to a star-forming regime, was selected by choice of the foam density. Gold foil and sphere motions were imaged by x-ray radiography. Radiographic images show the formation of an interface between rarefied gold and carbon plasmas, a shock moving into the sphere, and a blunting of the initial sphere's shape. Measurements show that the shock moved linearly around 64 μm/ns into the sphere, and the gold–carbon interface formed by 2 ns at the sphere edge remained stationary. The deformation of the sphere was driven by the incident radiation and not by mechanical pressures applied by gold plasma. The blunting of the sphere was likely due to the geometric reduction of flux near the sphere's poles. Higher x-ray flux near the sphere's equator caused high compression and a faster shock, which flattened the sphere. We will discuss the results and implications of our observations.

VanDervort, R. W. [University of Michigan 1 , Ann ↗

Forward sensitivity approach for estimating eddy viscosity closures in nonlinear model reduction

In this paper, we propose a variational approach to estimate eddy viscosity using forward sensitivity method (FSM) for closure modeling in nonlinear reduced order models. FSM is a data assimilation technique that blends model's predictions with noisy observations to correct initial state and/or model parameters. We apply this approach on a projection based reduced order model (ROM) of the one-dimensional viscous Burgers equation with a square wave defining a moving shock, and the two-dimensional vorticity transport equation formulating a decay of Kraichnan turbulence. We investigate the capability of the approach to approximate an optimal value for eddy viscosity with different measurement configurations. Specifically, we show that our approach can sufficiently assimilate information either through full field or sparse noisy measurements to estimate eddy viscosity closure to cure standard Galerkin reduced order model (GROM) predictions. Therefore, our approach provides a modular framework to correct forecasting error from a sparse observational network on a latent space. We highlight that the proposed GROM-FSM framework is promising for emerging digital twin applications, where real-time sensor measurements can be used to update and optimize surrogate model's parameters.

42 ENGINEERING↗

Shock-induced kinetics and cellular structures of liquid nitromethane detonation

Using a combination of high-speed diagnostics, optical pyrometry, velocimetry and video photography, we examine the spatiotemporal reaction kinetics of a prototypical high-energy explosive, liquid nitromethane (NM), as we drive this high explosive into detonation. The detonations were initiated by powerful shock waves whose durations (4 ns) were shorter than the characteristic time associated with the reaction (7 ns). Simple optical spectroscopy alone cannot characterize the kinetics with high time resolution because the reactants and products are flowing at a high velocity of ~6 km/s (6 µm/ns). Additionally, every detection volume behind the shock front contains material initiated over a range of times. High spatial resolution can be obtained by probing interfaces where the shock enters or breaks out of the NM. In addition, we obtain cellular patterns imprinted on the luminous shock front by the two-stage explosion in NM. These spatiotemporal patterns arise naturally as a result of the asymmetry produced by the moving shock front. In this way the shock front serves as a thin moving intrinsic optical gauge that reports and characterizes the two-stage NM explosion behind the front.

42 ENGINEERING↗

Reduced order models for Lagrangian hydrodynamics

It is reported, as a mathematical model of high-speed flow and shock wave propagation in a complex multimaterial setting, Lagrangian hydrodynamics is characterized by moving meshes, advection-dominated solutions, and moving shock fronts with sharp gradients. These challenges hinder the existing projection-based model reduction schemes from being practical. We develop several variations of projection-based reduced order model techniques for Lagrangian hydrodynamics by introducing three different reduced bases for position, velocity, and energy fields. A time-windowing approach is also developed to address the challenge imposed by the advection-dominated solutions. Lagrangian hydrodynamics is formulated as a nonlinear problem, which requires a proper hyper-reduction technique. Therefore, we apply the over-sampling DEIM and SNS approaches to reduce the complexity due to the nonlinear terms. Finally, we also present both a posteriori and a priori error bounds associated with our reduced order model. We compare the performance of the spatial and time-windowing reduced order modeling approaches in terms of accuracy and speed-up with respect to the corresponding full order model for several numerical examples, namely Sedov blast, Gresho vortices, Taylor-Green vortices, and triple-point problems.

97 MATHEMATICS AND COMPUTING↗

Dynamics of Three-Dimensional Shock-Wave/Boundary-Layer Interactions

Advances in measuring and understanding separated, nominally two-dimensional (2D) shock-wave/turbulent-boundary-layer interactions (STBLI) have triggered recent campaigns focused on three-dimensional (3D) STBLI, which display far greater configuration diversity. Nonetheless, unifying properties emerge for semi-infinite interactions, taking the form of conical asymptotic behavior where shock-generator specifics become insignificant. The contrast between 2D and 3D separation is substantial; the skewed vortical structure of 3D STBLI reflects the essentially 2D influence of the boundary layer on the 3D character of the swept shock. As with 2D STBLI, conical interactions engender prominent spectral content below that of the turbulent boundary layer. However, the uniform separation length scale, which is crucial to normalizing the lowest-frequency dynamics in 2D STBLI, is absent. Comparatively, the spectra of 3D STBLI are more representative of the mid-frequency, convective, shear-layer dynamics in 2D, while phenomena associated with 2D separation-shock breathing are muted. Asymptotic behavior breaks down in many regions important to 3D-STBLI dynamics, occurring in a configuration-dependent manner. Aspects of inceptive regions near shock generators and symmetry planes are reviewed. Focused efforts toward 3D modal and nonmodalanalyses, moving-shock/boundary-layer interactions, fluid/structure interactions, and flow control are suggested as directions for future work.

42 ENGINEERING↗

Effects of Coronal Magnetic Field Configuration on Particle Acceleration and Release during the Ground Level Enhancement Events in Solar Cycle 24

Ground level enhancements (GLEs) are extreme solar energetic particle (SEP) events that are of particular importance in space weather. In solar cycle 24, two GLEs were recorded on 2012 May 17 (GLE 71) and 2017 September 10 (GLE 72), respectively, using a range of advanced modern instruments. Here we conduct a comparative analysis of the two events by focusing on the effects of large-scale magnetic field configuration near active regions on particle acceleration and release. Although the active regions are both located near the western limb, temporal variations of SEP intensities and energy spectra measured in situ display different behaviors at early stages. By combining a potential field model, we find the coronal mass ejection (CME) in GLE 71 originated below the streamer belt, while in GLE 72 it originated near the edge of the streamer belt. We reconstruct the CME shock fronts with an ellipsoid model based on nearly simultaneous coronagraph images from multiple viewpoints and further derive the 3D shock geometry at the GLE onset. The highest-energy particles are primarily accelerated in the shock–streamer interaction regions, i.e., likely at the nose of the shock in GLE 71 and the eastern flank in GLE 72, due to quasi-perpendicular shock geometry and confinement of closed fields. Subsequently, they are released to the field lines connecting to near-Earth spacecraft when the shocks move through the streamer cusp region. This suggests that magnetic structures in the corona, especially shock–streamer interactions, may have played an important role in the acceleration and release of the highest-energy particles in the two events.

79 ASTRONOMY AND ASTROPHYSICS↗

Weak shock compaction on granular salt

This study conducted integrated experiments and computational modeling to investigate the speeds of a developing shock within granular salt and analyzed the effect of various impact velocities up to 245 m/s. Experiments were conducted on table salt utilizing a novel setup with a considerable bore length for the sample, enabling visualization of a moving shock wave. Experimental analysis using particle image velocimetry enabled the characterization of shock velocity and particle velocity histories. Mesoscale simulations further enabled advanced analysis of the shock wave’s substructure. In simulations, the shock front’s precursor was shown to have a heterogeneous nature, which is usually modeled as uniform in continuum analyses. The presence of force chains results in a spread out of the shock precursor over a greater ramp distance. With increasing impact velocity, the shock front thickness reduces, and the precursor of the shock front becomes less heterogeneous. Furthermore, mesoscale modeling suggests the formation of force chains behind the shock front, even under the conditions of weak shock. This study presents novel mesoscale simulation results on salt corroborated with data from experiments, thereby characterizing the compaction front speeds in the weak shock regime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Richtmyer–Meshkov instability when a shock is reflected for fluids with arbitrary equation of state

First predicted by Richtmyer in 1960 and experimentally confirmed by Meshkov in 1969, the Richtmyer–Meshkov instability (RMI) is crucial in fields such as physics, astrophysics, inertial confinement fusion and high-energy-density physics. These disciplines often deal with strong shocks moving through condensed materials or high-pressure plasmas that exhibit non-ideal equations of state (EoS), thus requiring theoretical models with realistic fluid EoS for accurate RMI simulations. Approximate formulae for asymptotic growth rates, like those proposed by Richtmyer, are helpful but rely on heuristic prescriptions for compressible materials. These prescriptions can sometimes approximate the RMI growth rate well, but their accuracy remains uncertain without exact solutions, as the fully compressible RMI growth rate is influenced by both vorticity deposited during shock refraction and multiple sonic wave refractions. This study advances previous work by presenting an analytic, fully compressible theory of RMI for reflected shocks with arbitrary EoS. It compares theoretical predictions with heuristic prescriptions using ideal gas, van der Waals gas and three-term constitutive equations for simple metals, the latter being analysed with detailed and simplified ideal-gas-like EoS. We additionally offer an alternative explicit approximate formula for the asymptotic growth rate. The comprehensive model also incorporates the effects of constant-amplitude acoustic waves at the interface, associated with the D'yakov–Kontorovich instability in shocks.

Napieralski, Mario (ORCID:0009000692344901)↗

Detection of visible-wavelength aurora on Mars

Mars hosts various auroral processes despite the planet’s tenuous atmosphere and lack of a global magnetic field. To date, all aurora observations have been at ultraviolet wavelengths from orbit. We describe the discovery of green visible-wavelength aurora, originating from the atomic oxygen line at 557.7 nanometers, detected with the SuperCam and Mastcam-Z instruments on the Mars 2020 Perseverance rover. Near–real-time simulations of a Mars-directed coronal mass ejection (CME) provided sufficient lead-time to schedule an observation with the rover. The emission was observed 3 days after the CME eruption, suggesting that the aurora was induced by particles accelerated by the moving shock front. To our knowledge, detection of aurora from a planetary surface other than Earth has never been reported, nor has visible aurora been observed at Mars. This detection demonstrates that auroral forecasting at Mars is possible, and that during events with higher particle precipitation, or under less dusty atmospheric conditions, aurorae will be visible to future astronauts.

Science & Technology - Other Topics↗

A particle resolved simulation approach for studying shock interactions with moving, colliding solid particles

This work applies a new combination of techniques for the fully resolved simulation of compressible, gas–particle multiphase flows. Here, the adaptive wavelet collocation method is used to dynamically, and efficiently, adapt the computational grid to localized flow features and the particles. A characteristic-based volume penalization method that imposes arbitrary Dirichlet, Neumann, or Robin-type immersed boundary conditions, is used to enforce the no-slip condition at particle surfaces. A hard-sphere collision model is applied to capture the particle–particle collisions. Proof of concept test cases are presented, showcasing the dynamic grid adaptation and fully resolved two-way coupling between the phases that is possible with this approach. Results for a shock-driven single cylinder under viscous and inviscid conditions are presented along with a demonstration of a shock interacting with a cloud of randomly distributed cylinders and spheres.

42 ENGINEERING↗

Broadband multiwavelength properties of the archetypal blazar 3C 279 during the 2017 Event Horizon Telescope campaign

The archetypal blazar 3C 279 hosts a prominent relativistic jet and exhibits strong broadband variability across the electromagnetic spectrum. In April 2017, the Event Horizon Telescope (EHT) observed 3C 279, alongside one of the most extensive quasi-simultaneous multiwavelength (MWL) campaigns ever conducted. With the aim of investigating the physical processes governing 3C 279, we analyzed individual observations and multiband light curves, and constructed a new quasi-simultaneous MWL spectrum. We also performed phenomenological modeling using the turbulent extreme multi-zone (TEMZ) model to constrain the fundamental physical properties of the source. The EHT observations reveal a clear flux increase in the innermost core between April 5 and 11, 2017. Over a broader timescale, radio measurements at longer wavelengths show concurrent enhancements in core flux and polarization around mid-April, coinciding with the ejection of a superluminal knot. Record UV-optical flares with strong polarization variability occurred in late March, followed by gamma-ray activity that declined before the end of the EHT observing period. During this interval, the source remained in a low X-ray state and showed no detectable VHE emission. The TEMZ modeling suggests that the broadband spectrum and variability of 3C 279 can be explained within a jet scenario in which turbulent plasma cells are compressed by a stationary conical shock. However, alternative interpretations, such as magnetic reconnection or a moving shock-in-jet event, remain plausible. This coordinated MWL campaign advances our understanding of the origin of jet and gamma-ray emission in 3C 279, while also providing a comprehensive publicly available dataset that will serve as a valuable reference for future studies.

Principe, G. [Trieste U.; INFN, Trieste; Bologna, ↗

Prediction of Detonation-Induced Disturbances Propagating Upstream into Inlets of Rotating Detonation Combustors

Disturbances caused by the detonation wave in a rotating detonation combustor (RDC) propagate upstream through the inlet, and can potentially affect and couple to upstream components, such as turbomachinery or isolators. These disturbances can potentially also affect the operation of the RDC itself. By drawing from the analogy of a detonation wave bounded by an inert gas, the pressure disturbances observed upstream of the inlet are explained as the consequence of the passage of an upstream propagating oblique shock. In this study, the pressure rise in the plenum from the oblique shock is measured in an axial air inlet RDC. The speed of the upstream propagating wave is estimated to be moderately above the acoustic speed of the oxidizer in the plenum. The wave propagates into the plenum despite local regions of choking in the inlet. It is estimated that the time it takes a fluid particle to transit from the plenum to the detonation channel through the inlet is much larger than the rotational time of the detonation wave. This implies that a fluid particle experiences multiple shocks prior to entering the detonation channel. The oblique shock propagating upstream through the inlet area change is modeled by leveraging an analogy with a quasi-1D shock wave moving in a variable area duct with mean (incoming) flow. Due to flow expansion along the area change, fluid particles are found to experience stronger shocks in the inlet than in the plenum, thereby creating different thermodynamic states within the fill region as the oxidizer emerges from the inlet.

Feleo, Alexander↗

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↗

Simulation to a Newborn Supernova Remnant from a Low-mass Iron Core Star

Supernova remnant observations show a high degree of asymmetry, mixing, and inhomogeneity. These asymmetries are seeded during the early seconds of the explosion and are further enhanced and modified as the shock and ejecta move through the stellar progenitor and into the circumstellar medium. We present simulations of a 9.6 M⊙ zero-metallicity progenitor initialized after shock revival and evolved for several years when the ejecta is in the circumstellar medium. A suite of 1D and 2D simulations examines the effects of neutron-star wind and radioactive decay heating. In 1D, decay heating forms a low-density bubble that suppresses the reverse shock. While in 2D, the heating is localized to metal-rich pockets, inflating them and compressing the surrounding material into dense shells. In 3D, the neutron-star wind and decay heating modify the plume morphology, producing more large-scale structures. The extended plume morphology leads to an asymmetrical shock breakout. After breakout, the leading plumes cannot keep up with the shock front, resulting in deceleration and fragmentation by the reverse shock while retaining the large-scale asymmetry. The projected ejecta morphology and velocities are strongly viewing angle dependent. The relatively uniform metal-rich distribution does not resemble the strongly inhomogeneous ejecta structure of Cas A. The 160-isotope decay network shows that 24.4% of the radioactive heating comes from decay chains other than the canonical 56Ni chain. The low explosion energy, low 56Ni yield, and Ni/Fe ratio greater than unity suggest an observational signature similar to an electron capture supernova.

Neopane, Sudarshan [University of Tennessee (UT)]↗

Dynamical Modulation of Solar Flare Electron Acceleration due to Plasmoid-shock Interactions in the Looptop Region

A fast-mode shock can form in the front of reconnection outflows and has been suggested as a promising site for particle acceleration in solar flares. Recent developments in the study of magnetic reconnection have shown that numerous plasmoids can be produced in a large-scale current layer. Here we investigate the dynamical modulation of electron acceleration in the looptop region when plasmoids intermittently arrive at the shock by combining magnetohydrodynamics simulations with a particle kinetic model. As plasmoids interact with the shock, the looptop region exhibits various compressible structures that modulate the production of energetic electrons. The energetic electron population varies rapidly in both time and space. The number of 5–10 keV electrons correlates well with the compression area, while that of >50 keV electrons shows good correlation with the strong compression area but only moderate correlation with shock parameters. We further examine the impacts of the first plasmoid, which marks the transition from a quasi-steady shock front to a distorted and dynamical shock. The number of energetic electrons is reduced by ~20% at 15–25 keV and nearly 40% for 25–50 keV, while the number of 5–10 keV electrons increases. In addition, the electron energy spectrum above 10 keV evolves softer with time. We also find that double or even multiple distinct sources can develop in the looptop region when the plasmoids move across the shock. Our simulations have strong implications to the interpretation of nonthermal looptop sources, as well as the commonly observed fast temporal variations in flare emissions, including the quasi-periodic pulsations.

79 ASTRONOMY AND ASTROPHYSICS↗

Planetesimal Impact Vapor Plumes and Nebular Shocks Form Chondritic Mixtures

The origin of chondrules and the chondritic sedimentary rocks that dominate the meteoritic record is a long-standing problem in planetary science. Here, we develop a physical model for the formation of chondritic mixtures as an outcome of vaporizing collisions between planetesimals that were dynamically excited by the growth and migration of planets. We present calculations of nebular shock waves generated by impact vapor plumes and focus on aspects of the plume interaction with the nebular gas and dust that have been neglected in previous studies of impact ejecta. We find that when water dominates the vapor, the plumes are relatively cool. However, the plume expansion is supersonic and can drive strong shock waves in the dusty nebular gas. Portions of these nebular shock fronts initially melt nebular dust, forming chondrules that are coupled to the moving front. As the shock front expands and cools, the chondrules solidify, while the shock front entrains additional dust. Eventually, the plume expansion stalls and then hydrodynamically collapses, turbulently mixing variably processed dust and size-sorted chondrules. For probable impact parameters and nebular conditions during giant planet growth and migration, the impact-generated mixtures have characteristics that span the range observed in chondritic meteorites, providing an environment for rapid formation of chondritic assemblages after chondrule formation. Our impact vapor and nebular shocks model, named the IVANS model, links chondrule formation to the overall context of planet formation and provides a framework for interpreting the detailed chronological and geochemical record contained in chondritic meteorites.

79 ASTRONOMY AND ASTROPHYSICS↗

Laser-driven collisionless shock acceleration of ions from near-critical plasmas

Here, this paper overviews experimental and numerical results on acceleration of narrow energy spread ion beams by an electrostatic collisionless shockwave driven by 1 μm (Omega EP) and 10 μm (UCLA Neptune Laboratory) lasers in near critical density CH and He plasmas, respectively. Shock waves in CH targets produced high-energy ~50 MeV protons (ΔE/E of ≤30%) and 314 MeV C 6+ ions (ΔE/E of ≤10%). Observation of acceleration of both protons and carbon ions to similar velocities is consistent with reflection of particles off the moving potential of a shock front. For shocks driven by CO 2 laser in a gas jet, ~30 MeV peak in He ion spectrum was detected. Particle-in-cell simulations indicate that regardless of the target further control over its density profile is needed for optimization of accelerated ion beams in part of energy spread, yield and maximum kinetic energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗