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 1,009 records · Page 56

Carbon stardust: From soot to diamonds

The formation of carbon dust in the outflow from stars and the subsequent evolution of this so called stardust in the interstellar medium is reviewed. The chemical and physical processes that play a role in carbon stardust formation are very similar to those occurring in sooting flames. Based upon extensive laboratory studies of the latter, the structure and physical and chemical properties of carbon soot are reviewed and possible chemical pathways towards carbon stardust are discussed. Grain-grain collisions behind strong interstellar shocks provide the high pressures required to transform graphite and amorphous carbon grains into diamond. This process is examined and the properties of shock-synthesized diamonds are reviewed. Finally, the interrelationship between carbon stardust and carbonaceous meteorites is briefly discussed.

Tielens, Alexander G. G. M.↗

Parabolized Navier-Stokes methods for hypersonic flows

A representative sampling of the techniques used in the integration of the Parabolized Navier-Stokes (PNS) equations is presented. Special atention is given to recent algorithms developed specifically for application to high speed flows, characterized by the presence of strong embedded shock waves and real gas effects. It is shown that PNS solvers are being used in the analysis of sonic boom signatures. Methods for modeling physical effects are discussed, including an overview of commonly used turbulence models and a more detailed discussion of techniques for including equilibrium and finite rate real gas effects.

Lawrence, Scott L.↗

The Space-Time Conservation Element and Solution Element Method: A New High-Resolution and Genuinely Multidimensional Paradigm for Solving Conservation Laws: The Two Dimensional Time Marching Schemes - 1

A new high resolution and genuinely multidimensional numerical method for solving conservation laws is being, developed. It was designed to avoid the limitations of the traditional methods. and was built from round zero with extensive physics considerations. Nevertheless, its foundation is mathmatically simple enough that one can build from it a coherent, robust. efficient and accurate numerical framework. Two basic beliefs that set the new method apart from the established methods are at the core of its development. The first belief is that, in order to capture physics more efficiently and realistically, the modeling, focus should be placed on the original integral form of the physical conservation laws, rather than the differential form. The latter form follows from the integral form under the additional assumption that the physical solution is smooth, an assumption that is difficult to realize numerically in a region of rapid chance. such as a boundary layer or a shock. The second belief is that, with proper modeling of the integral and differential forms themselves, the resulting, numerical solution should automatically be consistent with the properties derived front the integral and differential forms, e.g., the jump conditions across a shock and the properties of characteristics. Therefore a much simpler and more robust method can be developed by not using the above derived properties explicitly.

Chang, Sin-Chung↗

Low Dissipative High Order Numerical Simulations of Supersonic Reactive Flows

The objective of this paper is to evaluate the performance of a newly developed low dissipative sixth-order spatial and fourth-order temporal scheme for viscous reactive flows interacting with shock waves that contain fine scale flow structures. The accuracy and efficiency of the scheme, and to what degree the scheme can capture the correct physical wave speeds of stiff reactive flows will be included.

Sjoegreen, B.↗

Experimental Data from the Benchmark SuperCritical Wing Wind Tunnel Test on an Oscillating Turntable

The Benchmark SuperCritical Wing (BSCW) wind tunnel model served as a semi-blind testcase for the 2012 AIAA Aeroelastic Prediction Workshop (AePW). The BSCW was chosen as a testcase due to its geometric simplicity and flow physics complexity. The data sets examined include unforced system information and forced pitching oscillations. The aerodynamic challenges presented by this AePW testcase include a strong shock that was observed to be unsteady for even the unforced system cases, shock-induced separation and trailing edge separation. The current paper quantifies these characteristics at the AePW test condition and at a suggested benchmarking test condition. General characteristics of the model's behavior are examined for the entire available data set.

Heeg, Jennifer↗

Physics of Interplanetary Dust Collection with Aerogel

This report presents the results of research undertaken to study various problems associated with hypervelocity capture of dust particles in aerogel. The primary topics investigated were the properties of shocked aerogel and the requirements for reliable capture of particles on the STARDUST mission. In particular, the viscosity of shocked aerogel has been an open question. The results presented here suggest that the viscosity of aerogel at high impact velocities is negligible, although there remains some uncertainty about lower velocities. The model adopted for viscosity treats the mixture of polymeric silica and decomposition products and finds that, for particle velocities of 6-7 km/s, the viscosity is similar to that typical of light gasses at STP. Expressions for the Hugoniot of aerogel as a function of density were also obtained from the available data. All aerogels of interest for cosmic dust collectors have very similar shock velocity-particle velocity Hugoniot curves. The strength behavior of aerogel for low-speed penetration was measured, but further work is needed to study the proper way to apply this to the issue of terminal deceleration of a dust particle. Preliminary calculations designed to maximize the penetration depths were performed to determine the required density of aerogel to reliably stop a particle in a 3 cm thickness of aerogel (the path length expected for a normal impact into the STARDUST collector). In order to stop a particle of density rho(sub p) and diameter d(sub p), the mean density of the aerogel collector should be no less than that given by the expression bar rho(sub 0) = 1.085 X 10(exp -4 )rho(sub p)d(sub p), for densities measured in g/ cu cm and the particle diameter measured in micrometers.

Anderson, William W.↗

Inflationary stimulated Raman scattering in shock-ignition plasmas

In the shock-ignition inertial confinement fusion scheme, high-intensity lasers propagate through an inhomogeneous coronal plasma, driving a shock designed to cause fuel ignition. During the high-intensity ignitor laser pulse, in the long scale length coronal plasma, back-scattered stimulated Raman scattering (SRS) is likely to be in the kinetic regime. In this work, we use one-dimensional particle-in-cell simulations to show that there is a non-linear frequency shift caused by kinetic effects, resulting in the growth of SRS in an inhomogeneous plasma far exceeding the predictions of fluid theory, the so-called inflationary SRS or iSRS. We find that iSRS occurs over a wide range of density scale lengths relevant to shock-ignition and other directly-driven inertial confinement fusion schemes. The presence of iSRS in shock-ignition plasmas has implications for the theoretical gains from shock-ignition inertial confinement fusion. Here, we quantify the intensity threshold for the onset of iSRS for shock-ignition relevant parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Kinetic simulation study of magnetized collisionless shock formation on a terawatt laser system

Perpendicular, magnetized, collisionless shocks in hydrogen and neon plasmas are studied with 2-D particle-in-cell simulations for parameters accessible to experiments on OMEGA EP. The simulations are performed with realistic ion-electron mass ratios by which the relative importance of different micro-instabilities can be accurately captured. The dispersion relation of the modified two-stream instability (MTSI), the main dissipation mechanism for these shocks, is used to find suitable parameters for upcoming experiments. Simulations show that magnetized collisionless shocks can be readily formed within a few tenths of an ion gyro-period in both hydrogen and neon gases, with a background magnetic field of 50 T, achievable using the magneto-inertial fusion electrical discharge system (MIFEDS). Here, a portion of the ions are reflected to the upstream region and accelerated in both shock normal and tangential directions, indicating the formation of a supercritical shock. Shock front reformation is seen in longer time 1-D simulations. The results show that the formation time and width of these shocks are determined by MTSI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Shock compression of crystalline TeO 2 to the high-pressure fluid regime: Insights from ab initio molecular dynamics simulations

The shock response of fully-dense and porous crystalline tellurium dioxide (TeO 2 ⁠) to the high-pressure and high-temperature fluid regime was investigated within the framework of density functional theory with Mermin’s generalization to finite temperatures. The principal and porous shock Hugoniot curves were predicted from canonical ab initio molecular dynamics (AIMD) simulations, with the phase space sampled along isotherms up to 80 000 K, for densities ranging from ρ = 3 to 17 g/cm 3 . The polymorphs investigated are α-TeO 2 paratellurite (⁠P4 1 2 1 2), TeO 2 cotunnite (⁠Pnma⁠), and TeO 2 post-cotunnite (⁠P2 1 /m⁠). Based on the discontinuity found in the calculated U s – u p slope of TeO 2 post-cotunnite at a shock velocity of U s ≃ 8.35km/s and a particle velocity of u p ≃ 3.64 km/s, the shock melting temperature and pressure are predicted to be ≃ 6500 K and ≃ 170 GPa. Results from the AIMD simulations are in line with the static compression data of TeO 2 paratellurite and cotunnite, and with the recent shock Hugoniot data for single-crystal α- TeO 2 for pressures up to 85 GPa, obtained using the inclined-mirror method and the velocity interferometer system for any reflector combined with powder gun and two-stage light-gas gun.

74 ATOMIC AND MOLECULAR PHYSICS↗

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↗

A simple model of a strong shock driven by a spherical or cylindrical piston

A simple model of piston-driven spherical and cylindrical shocks is suggested. The model is based on a consistent use of two factors: (a) an almost uniform pressure across the shocked layer and, (b) continuous geometrical stretching of the surface elements of the expanding piston. It turns out that for a uniform pre-shock medium the gas between the piston and the shock behaves essentially as an incompressible fluid. An algebraic equation for the shock vs piston position is obtained. Detailed evaluation of the accuracy of the proposed solution shows that its accuracy is a few percent for the adiabatic index γ = 5/3 (as in ideal plasma). A closed-form solution describing enhancement of a weak ambient magnetic field by the shock is presented. The proposed model of piston-driven shocks goes beyond the classical self-similar solutions in that it: (1) naturally covers an early, non-asymptotic dynamics and its transition to asymptotic regime; (2) allows for smooth radial density variation of an ambient gas of the form of bumps, dimples or ramps between two constant values; and (3) allows for smooth temporal variation of piston velocity of the form of bumps, dimples, or ramps. Furthermore, this simple and versatile model provides some new insights into a classical hydrodynamical problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Instability seeding mechanisms due to internal defects in inertial confinement fusion targets

Performance degradation in laser-driven inertial confinement fusion (ICF) implosions is caused by several effects, one of which is Rayleigh–Taylor instability growth. Defects in ICF targets, such as internal voids and surface roughness, create instability seeds in the shell as shocks propagate through the target. A comprehensive understanding of seeding mechanisms is essential to characterize the impact of target defects on inflight shell integrity and mass injection into the central, lower-density vapor region. An analysis of early-time behavior of both single-mode shell mass modulations and isolated voids is performed by examining the evolution of the acoustic waves launched by these target imperfections. A systematic study of localized perturbation growth as a function of defect placement and size is presented. Furthermore, the use of low-density ablator materials (such as foams) is suggested as a potential mitigation strategy to improve target robustness against the impact of defect-initiated growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Radiative cooling effects on reverse shocks formed by magnetized supersonic plasma flows

We study the structure of reverse shocks formed by the collision of supersonic, magnetized plasma flows driven by an inverse (or exploding) wire array with a planar conducting obstacle. We observe that the structure of these reverse shocks varies dramatically with wire material, despite the similar upstream flow velocities and mass densities. For aluminum wire arrays, the shock is sharp and well-defined, consistent with magneto-hydrodynamic theory. In contrast, we do not observe a well-defined shock using tungsten wires, and instead we see a broad region dominated by density fluctuations on a wide range of spatial scales. We diagnose these two very different interactions using interferometry, Thomson scattering, shadowgraphy, and a newly developed imaging refractometer that is sensitive to small deflections of the probing laser corresponding to small-scale density perturbations. We conclude that the differences in shock structure are most likely due to radiative cooling instabilities, which create small-scale density perturbations elongated along magnetic field lines in the tungsten plasma. These instabilities grow more slowly and are smoothed by thermal conduction in the aluminum plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Amplification of turbulence through multiple planar shocks

We study the amplification of isotropic, incompressible turbulence through multiple planar, collisional shocks, using analytical linear theory. There are two limiting cases we explore. The first assumes shocks occur rapidly in time such that the turbulence does not evolve between shocks. Whereas the second case allows enough time for turbulence to isotropize between each shock. For the latter case, through a quasi-equation-of-state, we show that the weak multishock limit is agnostic to the distinction between thermal and vortical turbulent pressures, like an isotropic volumetric compression. When turbulence does not return to isotropy between shocks, the generated anisotropy—itself a function of shock strength—can feedback on amplification by further shocks, altering choices for maximal or minimal amplification. In addition for this case, we find that amplification is sensitive to the shock ordering. As a result, we map how choices of shock strength can impact these amplification differences due to ordering, finding, for example, shock pairs which lead to identical mean postshock fields (density, temperature, pressure) but maximally distinct turbulent amplification.

Astrophysical fluid dynamics↗

Molecular dynamics simulations of inelastic x-ray scattering from shocked copper

By taking the spatial and temporal Fourier transforms of the coordinates of the atoms in molecular dynamics simulations conducted using an embedded-atom-method potential, we calculate the inelastic scattering of x rays from copper single crystals shocked along [001] to pressures of up to 70 GPa. Above the Hugoniot elastic limit, we find that the copious stacking faults generated at the shock front introduce strong quasi-elastic scattering (QES) that competes with the inelastic scattering signal, which remains discernible within the first Brillouin zone; for specific directions in reciprocal space outside the first zone, the QES dominates the inelastic signal overwhelmingly. The synthetic scattering spectra we generate from our Fourier transforms suggest that energy resolutions of order 10 meV would be required to distinguish inelastic from quasi-elastic scattering within the first Brillouin zone of shock-loaded copper. We further note that high-resolution inelastic scattering also affords the possibility of directly measuring particle velocities via the Doppler shift. These simulations are of relevance to future planned inelastic scattering experiments at x-ray Free Electron Laser facilities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Equation-of-state, sound speed, and reshock of shock-compressed fluid carbon dioxide

Mechanical equation-of-state data of initially liquid and solid CO2 shock-compressed to terapascal conditions are reported. Diamond-sapphire anvil cells were used to vary the initial density and state of CO2 samples that were then further compressed with laser-driven shock waves, resulting in a data set from which precise derivative quantities, including Grüneisen parameter and sound speed, are determined. Reshock states are measured to 800 GPa and map the same pressure-density conditions as the single shock using different thermodynamic paths. The compressibility data reported here do not support current density-functional-theory calculations, but are better represented by tabular equation-of-state models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Observation of breakdown wave mechanism in avalanche ionization produced atmospheric plasma generated by a picosecond CO 2 laser

Understanding the formation and long-timescale evolution of atmospheric plasmas produced by ultrashort, long-wavelength infrared (LWIR) pulses is an important but partially understood problem. Of particular interest are plasmas produced in air with a peak laser intensity ~10 12 W/cm 2 , the so-called clamping intensity observed in LWIR atmospheric guiding experiments where tunneling and multi-photon ionization operative at near-IR or shorter wavelengths are inoperative. We find that avalanche breakdown on the surface of aerosol (dust) particles can act to seed the breakdown of air observed above the 200 GW/cm 2 threshold when a train of 3 ps 10.6 μm laser pulses separated by 18 ps is used. The breakdown first appears at the best focus but propagates backward toward the focusing optic as the plasma density approaches critical density and makes forward propagation impossible. The velocity of the backward propagating breakdown can be as high as 10 9 cm/s, an order of magnitude greater than measured with ns pulse-produced breakdown, and can be explained rather well by the so-called breakdown wave mechanism. Transverse plasma expansion with a similar velocity is assisted by UV photoionization and is observed as a secondary longitudinal breakdown mechanism in roughly 10% of the shots. When a cm-size, TW power beam is propagated, interception of aerosol particles is guaranteed and several (40 cm -3 ) breakdown sites appear, each initially producing a near-critical density plasma. On a 10 ns–1 μs timescale, shockwaves from each site expand radially and coalesce to produce a large hot gas channel. The radial velocity of the expansion agrees well with the prediction of the blast wave theory developed for ultrafast atmospheric detonations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗