Modeling Hydrodynamic Instabilities, Shocks, and Radiation Waves in High Energy Density Experiments [Slides]
Abstract not provided.
SEARCH · Engineering Papers
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.
Abstract not provided.
Mechanophores (MPs), stimuli-responsive molecules that respond chromatically to mechanochemical reactions, are important for understanding the coupling between mechanics and chemistry as well as in engineering applications. However, the atomic-level understanding of their activation originates from gas phase studies or under simple linear elongation forces directly on molecules or polymer chains containing MPs. The effect of many-body distortions, pervasive in condensed-phase applications, is not understood. Therefore, we performed large-scale molecular dynamics (MD) simulations of a poly(methyl methacrylate)-spiropyran copolymer under dynamic mechanical loading and studied the activation of the MP under various conditions from dynamical compression to tension during unloading. Detailed analysis of the all-atom MD trajectories shows that the MP blocks experience significant many-body intramolecular distortion that can significantly decrease the activation barrier as compared with when deformation rates are slow relative to molecular relaxation time scales. Here we find that the reactivity of MPs under material compression states is governed by many-body effects of intramolecular torsions, whereas under tension, the reactions are governed by tensile stresses.
Surface modification methods have been applied to metals and alloys to change the surface integrity, obtain superior mechanical properties, and improve service life irrespective of the field of application. In this review paper, current state-of-the-art of peening techniques are demonstrated. More specifically, classical and advanced shot peening (SP), ultrasonic impact peening (UIP), and laser shock peening (LSP) have been discussed. The effect of these techniques on mechanical properties, such as hardness, wear resistance, fatigue life, surface roughness, and corrosion resistance of various metals and alloys, are discussed. This study also reports the comparisons, advantages, challenges, and potential applications of these processes.
An axisymmetric full Navier-Stokes computational fluid dynamics (CFD) study was conducted to examine nozzle exhaust jet plume effects on the sonic boom signature of a supersonic aircraft. A simplified axisymmetric nozzle geometry, representative of the nozzle on the NASA Dryden NF-15B Lift and Nozzle Change Effects on Tail Shock (LaNCETS) research airplane, was considered. The highly underexpanded nozzle flow is found to provide significantly more reduction in the tail shock strength in the sonic boom N-wave pressure signature than perfectly expanded and overexpanded nozzle flows. A tail shock train in the sonic boom signature, similar to what was observed in the LaNCETS flight data, is observed for the highly underexpanded nozzle flow. The CFD results provide a detailed description of the nozzle flow physics involved in the LaNCETS nozzle at different nozzle expansion conditions and help in interpreting LaNCETS flight data as well as in the eventual CFD analysis of a full LaNCETS aircraft. The current study also provided important information on proper modeling of the LaNCETS aircraft nozzle. The primary objective of the current CFD research effort was to support the LaNCETS flight research data analysis effort by studying the detailed nozzle exhaust jet plume s imperfect expansion effects on the sonic boom signature of a supersonic aircraft. Figure 1 illustrates the primary flow physics present in the interaction between the exhaust jet plume shock and the sonic boom coming off of an axisymmetric body in supersonic flight. The steeper tail shock from highly expanded jet plume reduces the dip of the sonic boom N-wave signature. A structured finite-volume compressible full Navier-Stokes CFD code was used in the current study. This approach is not limited by the simplifying assumptions inherent in previous sonic boom analysis efforts. Also, this study was the first known jet plume sonic boom CFD study in which the full viscous nozzle flow field was modeled, without coupling to a sonic boom propagation analysis code, from the stagnation chamber of the nozzle to the far field external flow, taking into account all nonisentropic effects in the shocks, boundary layers, and free shear layers, and their interactions at distances up to 30 times the nozzle exit diameter from the jet centerline. A CFD solution is shown in Figure 2. The flow field is very complicated and multi-dimensional, with shock-shock and shockplume interactions. At the time of this reporting, a full three-dimensional CFD study was being conducted to evaluate the effects of nozzle vectoring on the aircraft tail shock strength.
Mixing of materials in porous media can cause a significant impact on fusion yield as previously demonstrated by the National Ignition Facility (NIF) MARBLE Campaign. Initially, the reactants are separated, with deuterium in the lattice struts and a tritium gas fill in the voids. Lattice parameters such as the strut thickness and relative pitch, provide a control for the mix parameters in the experiment. Los Alamos National Laboratory’s (LANL) BOSQUE project looks to better understand how the mix of the reactants and shell materials impact the fusion burn and resultant yield on various laser platforms. xRAGE’s Eulerian hydrodynamics and adaptive mesh refinement (AMR) provide the unique ability to study the impacts of multiscale features of complex lattice structures. This modeling provides the ability to measure shock front variations as the wave progress’ through a given media. Initial conditions of the lattice are essential to accurately model mix and burn measured by experiment. By varying the initial orientation and densities of these lattice regions the early time dynamics of how the shock is launched into the system is changed and advocates for the study of resulting effects. In this work, we will study the sensitivities of shock effects in varying 3D printed geometric systems and how these shocks alter the structure and mix in the lattice. We will discuss both preliminary experimental results and simulations to help plan and constrain future experiments where we will study the impact of different lattice geometries and lattice bulk densities. This work concludes with the relative impacts of lattice geometries on shock speeds at different bulk densities and the resultant mix widths due to those shock interactions. We see agreement with theory at the higher end of our bulk density study, however, as we approach lower bulk density systems the dynamics of these interactions begin to change.
Abstract Inner heliosheath (IHS) ions are expected to be exposed to various waves, shocks, and turbulence, which can affect ion distributions and thus their charge-exchange rates with interstellar neutral atoms. This work addresses the potential significance of electromagnetic ion cyclotron (EMIC) waves under expected IHS conditions. From a kinetic dispersion relation, we find the possibility of frequent triggering of EMIC instability in the IHS. The threshold anisotropy of proton temperatures required for the instability is small, ( T ⊥ – T ∣∣ )/ T ∣∣ ≈ 0.1 or less, mainly due to high plasma β (ratio of the plasma to magnetic pressures). Numerical calculations on the scattering of ions (protons, He + , and He 2+ with energy of 0.01–50 keV) based on two models for EMIC waves with a moderate intensity indicate significant scattering in the pitch angle (mostly a few tens of degrees) and energy (mostly a few tens of percent) although details depend on the energy and pitch angle of each species and adopted EMIC wave models. This occurs on a short timescale (<100 times the gyro-period of each ion species). Resonant scattering in a few to a few tens of keV (corresponding to the expected pickup ion energy in the IHS) is easily expected unless the wavenumber is too large. The scattering effect is distinguished among different species such that for lower gyrofrequency ions (He + versus He 2+ and He 2+ versus protons), the main scattering effect moves toward a lower energy domain. All these results imply continuous disturbance of ion distributions by possibly prevailing EMIC waves in the IHS.
Despite making up 5-20 wt.% of Earth’s predominantly iron core, the melting properties of elemental nickel at core conditions remain poorly understood, due largely to a dearth of experimental data. We present here an in situ X-ray diffraction study performed on laser shock-compressed samples of bulk nickel, reaching pressures up to ~ 500 GPa. Hugoniot states of nickel were targeted using a flat-top laser drive, with in situ X-ray diffraction data collected using the Linac Coherent Light Source. Rietveld methods were used to determine the densities of the shocked states from the measured diffraction data, while peak pressures were determined using a combination of measured particle velocities, shock transit times, hydrodynamic simulations, and laser intensity calibrations. We observed solid compressed face-centered cubic (fcc) Ni up to at least 332 ± 30 GPa along the Hugoniot—significantly higher than expected from the majority of melt lines that have been proposed for nickel. We also bracket the partial melting onset to between 377 ± 38 GPa and 486 ± 35 GPa.
A physical model is proposed for the estimation of the screech amplitude in underexpanded supersonic jets. The model is based on the hypothesis that the interaction of a plane acoustic wave with stationary shock waves provides amplification of the transmitted acoustic wave upon traversing the shock. Powell's discrete source model for screech incorporating a stationary array of acoustic monopoles is extended to accommodate variable source strength. The proposed model reveals that the acoustic sources are of increasing strength with downstream distance. It is shown that the screech amplitude increases with the fuiiy expanded jet Mach number. Comparisons of predicted screech amplitude with available test data show satisfactory agreement. The effect of variable source strength on directivity of the fundamental (first harmonic, lowest frequency mode) and the second harmonic (overtone) is found to be unimportant with regard to the principal lobe (main or major lobe) of considerable relative strength, and is appreciable only in the secondary or minor lobes (of relatively weaker strength
Engineering applications for aircraft noise prediction contain models for physical phenomenon that enable solutions to be computed quickly. These models contain parameters that have an uncertainty not accounted for in the solution. To include uncertainty in the solution, nondeterministic computational methods are applied. Using prediction models for supersonic jet broadband shock-associated noise, fixed model parameters are replaced by probability distributions to illustrate one of these methods. The results show the impact of using nondeterministic parameters both on estimating the model output uncertainty and on the model spectral level prediction. In addition, a global sensitivity analysis is used to determine the influence of the model parameters on the output, and to identify the parameters with the least influence on model output.
Presented is a comprehensive, up-to-date review of the shock-wave boundary-layer interaction problem. A detailed physical description of the phenomena for transonic and supersonic speed regimes is given based on experimental observations, correlations, and theoretical concepts. Approaches for solving the problem are then reviewed in depth. Specifically, these include: global methods developed to predict sudden changes in boundary-layer properties; integral or finite-difference methods developed to predict the continuous evolution of a boundary-layer encountering a pressure field induced by a shock wave; coupling methods to predict entire flow fields; analytical methods such as multi-deck techniques; and finite-difference methods for solving the time-dependent Reynolds-averaged Navier-Stokes equations used to predict the development of entire flow fields. Examples are presented to illustrate the status of the various methods and some discussion is devoted to delineating their advantages and shortcomings. Reference citations for the wide variety of subject material are provided for readers interested in further study.
The effects of different initial density distributions on the evolution of buoyancy-driven homogeneous variable-density turbulence (HVDT) at low (0.05) and high (0.75) Atwood numbers are studied by using high-resolution direct numerical simulations. HVDT aims to mimic the acceleration-driven Rayleigh–Taylor and shock-driven Richtmyer–Meshkov instabilities and reveals new physics that arise from variable-density effects on the turbulent mixing. Here, the initial amounts of pure light and pure heavy flows are altered primarily to mimic the variable-density turbulence at the different locations of the Rayleigh–Taylor and Richtmyer–Meshkov instabilities’ mixing layers where the amounts of the mixing fluids are not equal. It is found that for the low Atwood number cases, the asymmetric initial density distribution has limited effects on both global and local flow evolution for HVDT. However, at high Atwood number, both global flow evolution and the local flow structures are strongly affected by the initial composition ratio. The flow composed of more light fluid reaches higher turbulent levels and the local statistics reach their fully-developed behavior earlier in the time evolution. In contrast, during the late time decay, where most of the flow is well-mixed, all parameters become independent of the initial composition ratio for both low and high Atwood number cases.
We have carried out high-resolution near-infrared spectroscopic observations toward 16 Galactic supernova remnants (SNRs) showing strong H 2 emission features. A dozen bright H 2 emission lines are clearly detected for individual SNRs, and we have measured their central velocities, line widths, and fluxes. For all SNRs except one (G9.9–0.8), the H 2 line ratios are well consistent with that of thermal excitation at T ~ 2000 K, indicating that the H 2 emission lines are most likely from shock-excited gas and therefore that they are physically associated with the remnants. The kinematic distances to the 15 SNRs are derived from the central velocities of the H 2 lines using a Galactic rotation model. We derive for the first time the kinematic distances to four SNRs: G13.5+0.2, G16.0–0.5, G32.1–0.9, and G33.2–0.6. Among the remaining 11 SNRs, the central velocities of the H 2 emission lines for six SNRs are well consistent (±5 km s -1 ) with those obtained in previous radio observations, while for the other five SNRs (G18.1–0.1, G18.9–1.1, Kes 69, 3C 396, W49B) they are significantly different. We discuss the velocity discrepancies in these five SNRs. In G9.9–0.8, the H 2 emission shows nonthermal line ratios and narrow line width (~4 km s -1 ), and we discuss its origin.
The prospect that density waves and galactic shock waves are present on the large-scale in disk-shaped galaxies has received support in recent years from both theoretical and observational studies. Large-scale galactic shock waves in the interstellar gas are suggested to play an important, governing role in star formation, molecule formation, and the degree of development of spiral structure. Through the dynamics of the interstellar gas and the galactic shock wave phenomenon, a new insight into the physical basis underlying the morphological classification system of galaxies is suggested.
The prospect that density waves and galactic shock waves are present on the large scale in disk shaped galaxies has received support in recent years from both theoretical and observational studies. Large-scale galactic shock waves in the interstellar gas are suggested to play an important governing role in star formation, molecule formation, and the degree of development of spiral structure. Through the dynamics of the interstellar gas and the galactic shock-wave phenomenon, a new insight into the physical basis underlying the morphological classification system of galaxies is suggested.
Numerical computations are presented which illustrate and test various effects pertinent to the amplification and generation of turbulence in shock wave turbulent boundary layer interactions. Several fundamental physical mechanisms are identified. Idealizations of these processes are examined by nonlinear numerical calculations. The results enable some limits to be placed on the range of validity of existing linear theories.
The dark fine grained matrix of gas-rich ordinary chondrites replicates many of the physical, morphological, and spectral characteristics of the highly shocked and optically altered black chondrites. Spectral mixture modeling shows that the darkening and spectral attenuation seen in the dark matrix can be simulated with realistic mass fractions of light host material and black chondritic material. All these factors point to the conclusion that the dark matrix of gas-rich ordinary chondrites is dark due to the same processes that darkens black chondrites, shock distributed small particle size FeNi metal and troilite. Because the darkening is not seen in any of the non-gas-rich light portions and is only seen in the gas rich grains of the meteorite, the shock darkening would have to occur as part of the matrix's exposure to regolith processes. Since all gas-rich grains are darkened, it follows that darkening is not only common, but pervasive in asteroidal regoliths. These results imply that the upper, optically active layer of an ordinary chondrite parent body should have the spectral characteristics of a black chondrite, which are a dark, relatively featureless spectrum with modest red slope in the infrared. These are the characteristics of special type C asteroids.
A version of the dynamical adjustment procedure is generalized in order to make the procedure applicable to a diabatic model. In this new initialization procedure, backward adiabatic model integration is followed by forward diabatic model integration, with a high-frequency (low-pass) time filter in the form of the Euler backward time-differencing scheme applied throughout the integration. The new procedure was applied to the Goddard Laboratory for Atmospheres 4D data-assimilation system. It was found that, right from the beginning of the forecast integration, the forecast tendencies (and fields) were free of any noise due to imbalance in initial conditions, and the shocks related to an initial imbalance between model physics and dynamics, or the initial spinup effect, were practically removed.
The Active Magnetospheric Particle Tracer Explorer/Charge Composition Explorer (AMPTE/CCE) spacecraft, with an apogee of approximately 8.8R(sub E) and an inclination of approximately 4.3 deg, sampled the outer dayside equatorial magnetosphere for extended time periods and often crossed into the magnetosheath whenever the solar wind pressure was sufficiently high to compress the magnetopause to less than 8.8R(sub E). We have analyzed ion distributions on both sides of the magnetopause in order to investigate any local time and energy dependence, giving information about physical processes at the magnetopause and the bow shock. Particle measurements are from the Charge Energy Mass (CHEM)(1.5 to 300 keV/e) and Medium-Energy Particle Analyzer (MEPA)(10 keV to 2 MeV) instruments. The wide total energy range (1.5 keV to 2 MeV) covered describes the magnetospheric distributions quite well, and for the purpose of this study, adequately describes the high-energy part of the shocked solar wind. Thus both solar wind and magnetospheric components can be recognized in a mixed particle distribution. Case studies of representative magnetopause crossings at dawn, noon, and dusk, as well as a survey of several other crossings, indicated: (1) a local time and energy dependence of magnetosheath spectra at energies greater than or equal to 50 KeV; spectra were harder at the duskside than at the dawnside and also correlated with magnetospheric activity, (2) constantly much higher intensities in the magnetosphere than in the magnetosheath at energies greater than 10 KeV and an earthward gradient in the subsolar magnetosheath. In addition to the steady state magnetosheath population there exists a burst-type component indicative of a magnetospheric source, and most of the time this is recognized as a flux transfer event. Overall, the results about the origin of the greater than or equal to 50 KeV magnetosheath ions are consistent with the continuous leakage of magnetospheric particles across a tangential discontinuity magnetopause, locally distributed according to magnetospheric drift paths. Magnetic reconnection, although present, should not be a dominant source on average, because it is not continuous in time. Fermi acceleration should not be dominant because it predicts the opposite local time asymmetry, and shock drift acceleration should be a minor contributor at E greater than or equal to 50 keV because of upper-energy cutoff limitations. Our observations also indicate a significant magnetospheric contribution to energies as low as approximately 10 KeV, where the magnetosphere-magnetosheath intensity gradient reverses. However, in orderto examine the relative strength and local time distribution of all possible sources at these energies, a detailed analysis is required.