Large Eddy Simulations of a Mach 2 Expansion Compression Flow
Contains results on simulations of a Mach 2 expansion compression ramp which mirrors work done in the Tri Sonic Wind tunnel.
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.
Contains results on simulations of a Mach 2 expansion compression ramp which mirrors work done in the Tri Sonic Wind tunnel.
Alumina (Al 2 O 3 ) is an important ceramic material notable for its compressive strength and hardness. It represents one of the major oxide components of the Earth's mantle. Static compression experiments have reported evidence for phase transformations from the trigonal α-corundum phase to the orthorhombic Rh 2 O 3 (II)-type structure at ~90 GPa, and then to the post-perovskite structure at ~130 GPa, but these phases have yet to be directly observed under shock compression. In this work, we describe laser-driven shock compression experiments on polycrystalline alumina conducted at the Matter in Extreme Conditions endstation of the Linac Coherent Light Source. Ultrafast x-ray pulses (50 fs, 10 12 photons/pulse) were used to probe the atomic-level response at different times during shock propagation and subsequent pressure release. At 107 ± 8 GPa on the Hugoniot, we observe diffraction peaks that match the orthorhombic Rh 2 O 3 (II) phase with a density of 5.16 ± 0.03 g cm –3 . Upon unloading, the material transforms back to the α-corundum structure. Upon release to ambient pressure, densities are lower than predicted assuming isentropic release, indicating additional lattice expansion due to plastic work heating. Furthermore, using temperature values calculated from density measurements, we provide an estimate of alumina's strength on release from shock compression.
In this study, an analysis of toroidal plasma compression by a collapsing flux conserver reveals linearly stable scenarios of operation to high compression ratios. The resistive and ideal MHD stability is calculated in full toroidal geometry using the asymptotic matching method in realistic conditions. A time dependent MHD simulation of the compression is conducted to confirm the conservation principles used to calculate the equilibrium states for the analysis. The near edge current profile, controlled by toroidal field ramping during compression, is shown to be critical to stability due to coupling between poloidal components of the least stable mode. Two cases with slight differences in edge current are examined, one with a stable corridor to high compression and one without.
In this study, a toroidal plasma compressed by a collapsing flux conserver is analyzed to reveal stable scenarios of operation to high compression ratios. The resistive and ideal MHD stability is calculated in full toroidal geometry, using the asymptotic matching method in realistic conditions, and comparing with nonlinear simulations. The near edge current profile, controlled by toroidal field ramping during compression, is shown to be critical to stability due to coupling between poloidal components of the least stable mode. The extension of a length of shaft on axis is also found to be critical at high compression, as the resulting good curvature region in magnetic field stabilizes pressure driven modes that would otherwise be unstable. This work extends from previous studies, which initially showed the existence of a stable scenario, to include findings of more extensive stable zones, detailed effects of geometry, and nonlinear simulations of the instabilities. The nonlinear simulations of the compression are consistent with the linear analyses, confirming both the conservation and stability properties.
The goal of this proposal was to develop the capability to subject solids to ultrahigh compression using ramp loading techniques at NIF. Specifically, we proposed to examine the equation of state, strength, and phase transitions in iron and diamond to 30 Mbar pressures. Diamond and iron are fundamental materials of profound interest in condensed matter physics, materials science, and planetary science. The behavior of these solids under ultrahigh compression is essential to understanding planetary interiors and for probing the limits to our understanding of solid state behavior at extreme conditions.
We have conducted molecular dynamics (MD) simulations of quasi-isentropic ramp-wave compression to very high pressures over a range of strain rates from 10 11 down to 10 8 1/s. Using scaling methods, we collapse wave profiles from various strain rates to a master profile curve, which shows deviations when material response is strain-rate dependent. Thus, we can show with precision where, and how, strain-rate dependence affects the ramp wave. We find that strain rate affects the stress-strain material response most dramatically at strains below 20%, and that above 30% strain the material response is largely independent of strain rate. We show good overall agreement with experimental stress-strain curves up to approximately 30% strain, above which simulated response is somewhat too stiff. We postulate that this could be due to our interatomic potential or to differences in grain structure and/or size between simulation and experiment. Strength is directly measured from per-atom stress tensor and shows significantly enhanced elastic response at the highest strain rates. As a result, this enhanced elastic response is less pronounced at higher pressures and at lower strain rates.
Diamond is, by virtue of the covalent bonding between atoms and the very strong carbon to carbon bonds, the hardest natural material. It has been a fascinating material since its discovery, first as a decorative gem and more recently, for its numerous industrial uses because of its extreme hardness, elastic modulus, and optical transparency. In recent years, it has become a preferred ablator for laser shock experiments, and this has led to its choice as the capsule material for fusion experiments at the National Ignition Facility. Further, this review covers both experimental and computational (including machine learning) advancements in research on diamond subjected extreme conditions of temperature and pressure. The synergy between shock and ramp loading experiments and atomic level simulations is proving to be powerful in advancing our understanding of diamond under extremes.
Laser-driven ramp wave compression experiments have been used to investigate the strength (flow stress) of tantalum and other metals at high pressures and high strain rates. Recently this kind of experiment has been used to assess the dependence of the strength on the average grain size of the material, finding no detectable variation with grain size. The insensitivity to grain size has been understood theoretically to result from the dominant effect of the high dislocation density generated at the extremely high strain rates of the experiment. Here we review the experiments and describe in detail the multiscale strength model used to simulate them. The multiscale strength model has been extended to include the effect of geometrically necessary dislocations generated at the grain boundaries during compatible plastic flow in the polycrystalline metal. Lastly, we use the extended model to make predictions of the threshold strain rates and grain sizes below which grain size strengthening would be observed in the laser-driven Rayleigh-Taylor experiments.
Abstract Magnetic compression has long been proposed a promising method for plasma heating in a field reversed configuration (FRC). However, it remains a challenge to fully understand the physical mechanisms underlying the compression process, due to its highly dynamic nature beyond the one-dimensional (1D) adiabatic theory model (Spencer et al 1983 Phys. Fluids 26 1564). In this work, magnetohydrodynamics simulations on the magnetic compression of FRCs using the NIMROD code (Sovinec et al 2004 J. Comput. Phys. 195 355) and their comparisons with the 1D theory have been performed. The effects of the assumptions of the theory on the compression process have been explored, and the detailed profiles of the FRC during compression have been investigated. The pressure evolution agrees with the theoretical prediction under various initial conditions. The axial contraction of the FRC can be affected by the initial density profile and the ramping rate of the compression magnetic field, but the theoretical predictions on the FRC’s length in general and the relation r s = 2 r o in particular hold approximately well during the whole compression process, where r s is the major radius of FRC separatrix and r o is that of the magnetic axis. The evolutions of the density and temperature can be affected significantly by the initial equilibrium profile and the ramping rate of the compression magnetic field. During the compression, the major radius of the FRC is another parameter that is susceptible to the ramping rate of the compression field. Basically, for the same magnetic compression ratio, the peak density is higher and the FRC’s radius r s is smaller than the theoretical predictions.
Andra is currently implementing the detailed engineering phase of the step-wise design development of Cigeo (the French HLW and IL-LLW Deep Geological Repository). Various technical issues are at stake, including the requirements linked to 'Reversibility / Retrievability', at different stages of the repository life cycle (from its construction start-up to its final closure and monitoring). These reversibility requirements have a significant impact on the way Cigeo is designed, will be progressively developed, operated, later closed and monitored. On a shorter term, one must demonstrate that the waste container retrieval operations are technically safe and feasible. In that sense, the structural dimensioning of the ramps and mains, as well as that of the IL-LLW disposal vaults, is critical, since the concrete liner supporting the rock walls must last over a period of some 100 to 150 years (lapse of time during which retrieval operations are likely to take place). The present paper starts with a description of the geomechanical conditions prevailing in the Callovo-Oxfordian (COX) clay-stone formation (aka the argillite). The structural dimensioning of the IL-LLW disposal vault liner, as calculated at the end of the basic engineering phase is indicated. The thickness values obtained so far have led to considering an alternative design in which a compressible material is positioned between the concrete liner intrados and the rock walls: in such situation, the 'creeping effect' of the rock is mitigated, for a time, by the compressible material (nearly a perfect elastoplastic material), hence reducing the stress level exerted inside the concrete liner. This situation enables to significantly reduce the liner thickness and provides potentially significant saving in the construction of the works concerned (preliminary sensitivity results are given). The following part of the paper is then focused on a practical case story: a technological demonstration (still ongoing) of a tunnel boring experiment, implemented underground (in a dedicated drift of the Meuse Haute Marne URL infrastructures, at scale 1:2 of the Cigeo drifts/vaults). In this excavation and construction experiment, the tunnel boring machine (TBM) lays (via an erector) 'classical' concrete arch segments (forming a ring) in a first portion of the drift and 'compressible' ones (the compressible material is positioned on their extrados) in a second portion of the drift. The technical demonstration objective is two-fold: the capacity to lay compressible arch segments with a TBM is checked, their structural behavior (with time) is compared with that of 'classical' ones. The way the compressible material is produced and then fixed on the arch concrete segment is also documented, c/w a description of the complementary testing campaign to be carried out on a test bench (on surface) of a concrete ring (at a scale 1:2 of Cigeo) formed of compressible arch segments. This paper concludes with a timeline of the experiment. Should the results obtained provide some perspectives of technical and economic improvement (i.e. saving vis-a-vis the 'classical' solution) a decision could be taken by Andra to incorporate this alternative solution in the 'design freeze' of Cigeo underground infrastructures. The positive outcomes of this technical test campaign would also pave the way for specifying the Cigeo TBM characteristics (dimensions and handling equipment), scheduled for manufacturing order around 2020. (authors)
Presentation on use of grooved graded-density impactor to produce one-dimensional ramp waves for simplified isentropic compression experiments.
This study examines the qualitative features of an anelasticity model associated with the bowing of dislocations in the presence of phase transition. A simple physically plausible mechanism is introduced to describe the interaction of anelasticity and the transformation. Varying the anelastic parameters results in strong differences in the deviatoric stress response. The model is applied to study the behavior of tin (Sn) and compared to data from ramp driven compression-release experiments. Tin exhibits a complex phase diagram within a relatively accessible range of temperature and pressures and the characterization of its phases is considered an open problem with significant scientific merit. The coupling between anelasticity, plasticity, and phase transformation contributes to release wave features traditionally associated with the phase transition effect alone suggesting the importance of accounting for the effects jointly. Posterior distributions of the plastic and anelastic parameters are computed using Bayesian-inference-based methods, further highlighting the importance of anelasticity in this regime.
Recent X-ray observations of merger shocks in galaxy clusters have shown that the postshock plasma has two temperatures, with the protons hotter than the electrons. By means of two-dimensional particle-in-cell simulations, we study the physics of electron irreversible heating in low-Mach-number perpendicular shocks, for a representative case with sonic Mach number of 3 and plasma beta of 16. We find that two basic ingredients are needed for electron entropy production: (1) an electron temperature anisotropy, induced by field amplification coupled to adiabatic invariance; and (2) a mechanism to break the electron adiabatic invariance itself. In shocks, field amplification occurs at two major sites: at the shock ramp, where density compression leads to an increase of the frozen-in field; and farther downstream, where the shock-driven proton temperature anisotropy generates strong proton cyclotron and mirror modes. The electron temperature anisotropy induced by field amplification exceeds the threshold of the electron whistler instability. The growth of whistler waves breaks the electron adiabatic invariance and allows for efficient entropy production. For our reference run, the postshock electron temperature exceeds the adiabatic expectation by $\simeq 15 \% $, resulting in an electron-to-proton temperature ratio of $\simeq 0.45$. We find that the electron heating efficiency displays only a weak dependence on mass ratio (less than $\simeq 30 \% $ drop, as we increase the mass ratio from ${m}_{i}/{m}_{e}=49$ up to ${m}_{i}/{m}_{e}=1600$). We develop an analytical model of electron irreversible heating and show that it is in excellent agreement with our simulation results.
Evolutions of elastic-plastic waves have been recorded in three series of plate impact experiments with annealed vanadium samples under conditions of shockless and combined ramp and shock dynamic compression. The shaping of incident wave profiles was realized using intermediate base plates made of different silicate glasses through which the compression waves were entered into the samples. Measurements of the free surface velocity histories revealed an apparent growth of the Hugoniot elastic limit with decreasing average rate of compression. The growth was explained by “freezing” of the elastic precursor decay in the area of interaction of the incident and reflected waves. A set of obtained data show that the current value of the Hugoniot elastic limit and plastic strain rate is rather associated with the rate of the elastic precursor decay than with the local rate of compression. The study has revealed the contributions of dislocation multiplications in elastic waves. It has been shown that independently of the compression history the material arrives at the minimum point between the elastic and plastic waves with the same density of mobile dislocations.
Here, we present equations of state relevant to conditions encountered in ramp and multiple-shock compression experiments of water. These experiments compress water from ambient conditions to pressures as high as about 14 GPa and temperatures of up to several hundreds of Kelvin. Water may freeze into ice VII during this process. Although there are several studies on the thermodynamic properties of ice VII, an accurate and analytic free energy model from which all other properties may be derived in a thermodynamically consistent manner has not been previously determined. We have developed such a free energy model for ice VII that is calibrated with pressure-volume-temperature measurements and melt curve data. Furthermore, we show that liquid water in the pressure and temperature range of interest is well-represented by a simple Mie-Grüneisen equation of state. Our liquid water and ice VII equations of state are validated by comparing to sound speed and Hugoniot data. Although they are targeted towards ramp and multiple-shock compression experiments, we demonstrate that our equations of state also behave reasonably well at pressures and temperatures that lie somewhat beyond those found in the experiments.
Recent advances in the chemical process industry have allowed for the intensification of reactors and unit operations, by enhancing heat and mass transfer or combining multiple unit operations. Process intensification can enable chemical plants to be constructed in a more compact, modular fashion, offering improvements over conventional on-site approaches to capital construction, operations and maintenance. This modular chemical process intensification (MCPI) offers several benefits over conventional stick-built (CSB) plant construction in terms of reduced footprint, reduced energy consumption, lower cost, less waste and improved safety and quality control. However, acceptance of MCPI over CSB construction practices within the chemical industry can be impeded by the uncertain risks associated with investing in new technology. Furthermore, the work here documents a case study during the development of a modular chemical plant for capturing distributed energy resources within chemical production. This MCPI approach to plant construction is contrasted with a CSB approach for producing the same chemical product. Data collection tools were developed, based on a literature review, and data was collected from the technology developer to understand the process technology. Sensitivity analysis was then conducted to analyze the business rationale for the application of MCPI over CSB across several market scenarios. It was found that MCPI would be better suited for capacities up to 150 000 metric tons per year, but that improvement in payback period was needed. Additionally, for MCPI approach to achieve acceptable payback periods, efforts are needed to reduce the cost of capital equipment and compress the schedule for ramping up modular production.
In this paper, we report on the quasi-isentropic compression of an iron sample using ramp shaped laser irradiation. This technique allows us to quasi-isentropically compress iron up to 700 GPa and 8500 K. To our knowledge, these data are the highest pressures reached on iron in off-Hugoniot conditions and the closest to the thermodynamic states thought to exist in Earth-like planetary cores. The experiment was performed on the Ligne d'Intégration laser facility at CESTA, Bordeaux, France.