Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Ramp compression”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

SPE-4-prime Calculations on Sympathetic detonation and Preliminary Confinement [Slides]

Peak pressure is very low compared with input pressures that lead to detonation in most explosives (> 100 MPa). As an illustration, Gibbs and Popolato data for PETN at density 1.4 g/cc shows a 500 MPa shock transits to detonation only after a run distance of > 10 mm. The above data are consistent with the RP-83 design - the PETN in the detonator gets a real whack from the EBW to achieve prompt initiation. 500 MPa is 5 kilobar, which is a low pressure in the shock initiation arena. The 4 MPa input pressure is 0.04 kilobar, which is very low. After traveling through the intervening geologic media, the pressure will arrive as a ramp rather than a shock, which is much less effective at initiating explosives. In fact, a compression wave may desensitize the explosive. The most sensitive HE, the PETN in the detonator, is shielded by the mass of ANFO and the pentolite boosters, and so the pressure input to the detonator will be significantly less than that experienced by the ANFO. An additional factor is that the abandoned charge was thermally compromised, and the effect of the somewhat elevated temperature on the PETN and the SHANFO is most likely to desensitize them. This has been observed at LLNL with PETN, and the manufacturer of the SHANFO is the source of the information on its response to temperature. The only possibility posited by one of us is failure of the geologic media resulting in a puncture / perforation of the SHANFO charge. If sufficient shear energy were put into this charge, it could conceivably initiate. This is difficult to quantify and seems like a remote possibility, and at any rate is outside the expected behavior of the geologic media as I understand the problem. Information on the effect of temperature on SHANFO is in the enclosed report that analyzes the failure of SPE 4. I think there are additional possibilities that were not evaluated in the report, but it does provide some useful data on SHANFO sensitivity at elevated temperatures. I don't know if you have seen it, so I enclose a copy. PETN (Jon L. Maienschein, PhD Director, LLNL Energetic Materials Center)

58 GEOSCIENCES↗
DOE OSTI2017

Generation of Ramped Current Profiles in Relativistic Electron Beams Using Wakefields in Dielectric Structures

We show that temporal pulse tailoring of charged-particle beams is essential to optimize efficiency in collinear wakefield acceleration schemes. In this Letter, we demonstrate a novel phase space manipulation method that employs a beam wakefield interaction in a dielectric structure, followed by bunch compression in a permanent magnet chicane, to longitudinally tailor the pulse shape of an electron beam. This compact, passive, approach was used to generate a nearly linearly ramped current profile in a relativistic electron beam experiment carried out at the Brookhaven National Laboratory Accelerator Test Facility. Here, we report on these experimental results including beam and wakefield diagnostics and pulse profile reconstruction techniques.

43 PARTICLE ACCELERATORS↗
DOE OSTI2017

Axial magnetic field injection in magnetized liner inertial fusion

MagLIF is a fusion concept using a Z-pinch implosion to reach thermonuclear fusion. In current experiments, the implosion is driven by the Z-machine using 19 MA of electrical current with a rise time of 100 ns. MagLIF requires an initial axial magnetic field of 30 T to reduce heat losses to the liner wall during compression and to confine alpha particles during fusion burn. This field is generated well before the current ramp starts and needs to penetrate the transmission lines of the pulsed-power generator, as well as the liner itself. Consequently, the axial field rise time must exceed hundreds of microseconds. Any coil capable of being submitted to such a field for that length of time is inevitably bulky. The space required to fit the coil near the liner, increases the inductance of the load. In turn, the total current delivered to the load decreases since the voltage is limited by driver design. Yet, the large amount of current provided by the Z-machine can be used to produce the required 30 T field by tilting the return current posts surrounding the liner, eliminating the need for a separate coil. However, the problem now is the field penetration time, across the liner wall. This paper discusses why skin effect arguments do not hold in the presence of resistivity gradients. Numerical simulations show that fields larger than 30 T can diffuse across the liner wall in less than 60 ns, demonstrating that external coils can be replaced by return current posts with optimal helicity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Magnesium oxide at extreme temperatures and pressures studied with first-principles simulations

We combine two first-principles computer simulation techniques, path integral Monte Carlo and density functional theory molecular dynamics, to determine the equation of state of magnesium oxide in the regime of warm dense matter, with densities ranging from 0.35 to 71 g cm -3 and temperatures ranging from 10 000 K to 5 × 10 8 K. These conditions are relevant for the interiors of giant planets and stars as well as for shock wave compression measurements and inertial confinement fusion experiments. We study the electronic structure of MgO and the ionization mechanisms as a function of density and temperature. We show that the L-shell orbitals of magnesium and oxygen hybridize at high density. This results in a gradual ionization of the L-shell with increasing density and temperature. In this regard, MgO behaves differently from pure oxygen, which is reflected in the shape of the MgO principal shock Hugoniot curve. The curve of oxygen shows two compression maxima, while that of MgO shows only one. We predict a maximum compression ratio of 4.66 to occur for a temperature of 6.73 × 10 7 K. Lastly, we research how multiple shocks and ramp waves can be used to cover a large range of densities and temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2011

Test Problem: Tilted Rayleigh-Taylor for 2-D Mixing Studies

The 'tilted-rig' test problem originates from a series of experiments (Smeeton & Youngs, 1987, Youngs, 1989) performed at AWE in the late 1980's, that followed from the 'rocket-rig' experiments (Burrows et al., 1984; Read & Youngs, 1983), and exploratory experiments performed at Imperial College (Andrews, 1986; Andrews and Spalding, 1990). A schematic of the experiment is shown in Figure 1, and comprises a tank filled with light fluid above heavy, and then 'tilted' on one side of the apparatus, thus causing an 'angled interface' to the acceleration history due to rockets. Details of the configuration given in the next chapter include: fluids, dimensions, and other necessary details to simulate the experiment. Figure 2 shows results from two experiments, Case 110 (which is the source for this test problem) that has an Atwood number of 0.5, and Case 115 (a secondary source described in Appendix B), with Atwood of 0.9 Inspection of the photograph in Figure 2 (the main experimental diagnostic) for Case 110. reveals two main areas for mix development; 1) a large-scale overturning motion that produces a rising plume (spike) on the left, and falling plume (bubble) on the right, that are almost symmetric; and 2) a Rayleigh-Taylor driven mixing central mixing region that has a large-scale rotation associated with the rising and falling plumes, and also experiences lateral strain due to stretching of the interface by the plumes, and shear across the interface due to upper fluid moving downward and to the right, and lower fluid moving upward and to the left. Case 115 is similar but differs by a much larger Atwood of 0.9 that drives a strong asymmetry between a left side heavy spike penetration and a right side light bubble penetration. Case 110 is chosen as the source for the present test problem as the fluids have low surface tension (unlike Case 115) due the addition of a surfactant, the asymmetry small (no need to have fine grids for the spike), and there is extensive reasonable quality photographic data. The photographs in Figure 2 also reveal the appearance of a boundary layer at the left and right walls; this boundary layer has not been included in the test problem as preliminary calculations suggested it had a negligible effect on plume penetration and RT mixing. The significance of this test problem is that, unlike planar RT experiments such as the Rocket-Rig (Youngs, 1984), Linear Electric Motor - LEM (Dimonte, 1990), or the Water Tunnel (Andrews, 1992), the Tilted-Rig is a unique two-dimensional RT mixing experiment that has experimental data and now (in this TP) Direct Numerical Simulation data from Livescu and Wei. The availability of DNS data for the tilted-rig has made this TP viable as it provides detailed results for comparison purposes. The purpose of the test problem is to provide 3D simulation results, validated by comparison with experiment, which can be used for the development and validation of 2D RANS models. When such models are applied to 2D flows, various physics issues are raised such as double counting, combined buoyancy and shear, and 2-D strain, which have not yet been adequately addressed. The current objective of the test problem is to compare key results, which are needed for RANS model validation, obtained from high-Reynolds number DNS, high-resolution ILES or LES with explicit sub-grid-scale models. The experiment is incompressible and so is directly suitable for algorithms that are designed for incompressible flows (e.g. pressure correction algorithms with multi-grid); however, we have extended the TP so that compressible algorithms, run at low Mach number, may also be used if careful consideration is given to initial pressure fields. Thus, this TP serves as a useful tool for incompressible and compressible simulation codes, and mathematical models. In the remainder of this TP we provide a detailed specification; the next section provides the underlying assumptions for the TP, fluids, geometry details, boundary conditions (and alternative set-ups), initial conditions, and acceleration history (and ways to treat the acceleration ramp at the start of the experiment). This is followed by a section that defines data to be collected from the simulations, with results from the experiments and DNS from Livescu using the CFDNS code, and ILES simulations from Youngs using the compressible TURMOIL code and Andrews using the incompressible RTI3D code. We close the TP with concluding remarks, and Appendices that includes details of the sister Case 115, initial condition specifications for density and pressure fields. The Tilted-Rig Test Problem is intended to serve as a validation problem for RANS models, and as such we have provided ILES and DNS simulations in support of the test problem definition. The generally good agreement between experiment, ILES and DNS supports our assertion that the Tilted-Rig is useful, and the only 2-D TP that can be used to validate RANS models.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Structure of Shock Waves and Inelasticity in Shock-Compressed Cemented Tungsten Carbides

In this study, shock wave experiments are conducted on General Carbide cemented tungsten carbide with 3.7wt.% cobalt binder to determine its shock-induced compression behavior up to 100 GPa. The measured wave profiles indicate the cemented tungsten carbide to undergo elastic-plastic deformation during shock compression. A three-stage particle velocity profile is observed in the experiments -- an initial elastic-rise to the Hugoniot Elastic Limit (HEL), an elastic-plastic ramp indicating substantial post-yield hardening, and finally a rise to the peak shocked Hugoniot state. The results of the experiments are used to determine the HEL, the shock velocity (U s ) vs. particle velocity (u p ) Hugoniot relation, and the longitudinal stress (σ x ) vs. specific volume (V) curve for the samples. The HEL of the material was determined to lie between 4.41 GPa and 4.58 GPa. The U s - u p relation was determined to be U s = 4.97 + 1. 457u p for particle velocities greater than 0.75 km/s. The measured plastic shock velocities for particle velocities less than 0.7 km/s were found to be larger than those predicted using the linear U s - u p Hugoniot relationship, indicating the cemented WC samples to preserve substantial shear strength in the post-yield deformation region. No phase transformation was observed up to 100 GPa.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Shockless compression and release behavior of beryllium to 110 GPa

A magnetohydrodynamic loading technique was used to shocklessly compress beryllium to peak longitudinal stresses of 19–110 GPa and, subsequently, unload in order to determine both the compressive response and also the shear stress supported upon release. Loading strain rates were on the order of 10{sup 6 }s{sup −1}, while the unloading rates were nearly constant at 3 × 10{sup 5 }s{sup −1}. Velocimetry was used to monitor the ramp and release behavior of a beryllium/lithium fluoride window interface. After applying window corrections to infer in situ beryllium velocities, a Lagrangian analysis was employed to determine the material response. The Lagrangian wavespeed-particle velocity response is integrated to generate the stress-strain path, average change in shear stress over the elastic unloading, and estimates of the shear modulus at peak compression. These data are used to infer the pressure dependence of the flow strength at the unloading rate. Comparisons to several strength models reveal good agreement to 45 GPa, but the data indicate 20%–30% higher strength near 100 GPa.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Compression-induced solidification of shock-melted cerium

Compression induced solidification has been observed in cerium on nanosecond timescales. For this study, a series of experiments were conducted on Sandia National Laboratories' Z Facility in which cerium was shock melted and subsequently shocklessly, or ramp, loaded across the melt line inducing solidification. The signature of solidification manifested in the recovery of material strength and the propagation of waves at the local elastic sound velocity. Density functional theory simulations of cerium along the experimental phase-space path exhibit spontaneous freezing to a tetragonal phase at the same pressure and closely predict the observed physical properties of solid and liquid cerium near melt.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Shock response of cobalt-based cemented tungsten carbides at pressures up to 100 GPa

In the present study, shock wave experiments are conducted on General Carbide cemented tungsten carbides with 3.7 wt.% and 6.0 wt.% cobalt binder to determine their shock compression response up to 100 GPa. A three-stage particle velocity profile is observed in the experiments -- an initial elastic-rise to the Hugoniot elastic limit (HEL), an elastic-plastic ramp indicating post-yield hardening, and finally a rise to the peak shocked Hugoniot state. The results of the experiments are used to determine the HEL, shock velocity (U s ) vs. particle velocity (u p ) Hugoniot relationship, and the longitudinal stress vs. specific volume Hugoniot relationship for the two samples. Here, the HEL for the WC 3.7wt% Co and WC 6.0wt% Co samples was determined to be ~4.45 GPa and 3.72 GPa, respectively. The U s – u p relation was determined to be U s = 4.97(0.006)+1. 446 (0.018) u p for WC 3.7 wt.% Co and U s = 4.93(0.006)+1. 454(0.017)u p for the WC 6.0 wt.% Co sample at peak particle velocities>0.75 km/s. For both WC grades, in the particle velocity regime less than ~0.75 km/s, the measured shock wave velocities were found to be larger than those predicted by the linear U s –u p Hugoniot relationship, indicating the two WC samples to preserve substantial strength in the post-yield deformation regime. Both WC grades show a catastrophic drop in shear stress carrying capacity when shocked to longitudinal stresses greater than ~70 GPa.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Outer shell symmetry for double shell capsules with aluminum ablators

Double shell targets are a promising potential avenue to obtain robust neutron yield at current laser facilities. Similar to single shell designs, double shells require the symmetric implosion of an ablator in order to uniformly compress and heat a fuel volume, with the goal of achieving thermonuclear burn. Significant differences between double and single shells include the usage of an aluminum ablator as well as a reverse ramp laser pulse. In addition, double shells require a different convergence than single shells for fuel ignition. Numerical implosion studies at various energies with comparisons to experimental outcomes are required to gain confidence that simulations can capture the ablator shape from subscale to full scale. The current work builds on previous implosion experiments conducted at 1-MJ laser energy to confirm achieved ablator symmetry at 1.25 and 1.5 MJ. Average ablator P2 and P4 shapes measured in these experiments are within 5% of the simulated shape, which merits the platforms for further experimental studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Outer shell symmetry for double shell capsules with aluminum ablators

Double shell targets are a promising potential avenue to obtain robust neutron yield at current laser facilities. Similar to single shell designs, double shells require the symmetric implosion of an ablator in order to uniformly compress and heat a fuel volume, with the goal of achieving thermonuclear burn. Significant differences between double and single shells include the usage of an aluminum ablator as well as a reverse ramp laser pulse. In addition, double shells require a different convergence than single shells for fuel ignition. Numerical implosion studies at various energies with comparisons to experimental outcomes are required to gain confidence that simulations can capture the ablator shape from subscale to full scale. The current work builds on previous implosion experiments conducted at 1-MJ laser energy to con firm achieved ablator symmetry at 1.25 and 1.5 MJ. Average ablator P2 and P4 shapes measured in these experiments are within 5% of the simulated shape, which merits the platforms for further experimental studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Modeling of Thermal Decomposition of TATB-Based Explosive for Safety Analysis

We investigate and model the cook-off behavior of LX-17 to understand the response of explosive systems in abnormal thermal environments. Decomposition has been explored via conventional ODTX (One-Dimensional Time-to-eXplosion), PODTX (ODTX with pressure-measurement), TGA (Thermo-Gravimetric Analysis), and DSC (Differential Scanning Calorimetry) experiments under isothermal and ramped temperature profiles. The data were used to fit reaction rate parameters for proposed schemes in an ALE3D computational model. This model includes chemical reactions, thermo- and hydro-dynamics, and material properties, including thermal expansion, compressibility, and strength. These parameterizations were carried out utilizing a Python evolutionary optimization method on LLNL’s high-performance computing clusters. Additional experiments are being developed to further characterize and monitor decomposition intermediates to improve the model. Once experimentally validated, this model will be scalable to several applications involving LX-17. Furthermore, the optimization methodology developed herein should be applicable to other high explosive materials.

Chemistry - Chemical explosives↗
DOE OSTI2012

Control of post-disruption runaway electron beams in DIII-D

Recent experiments in the DIII-D tokamak have demonstrated real-time control and dissipation of post-disruption runaway electron (RE) beams. In the event that disruption avoidance, control, and mitigation schemes fail to avoid or suppress RE generation, active control of the RE beam may be an important line of defense to prevent the rapid, localized deposition of RE beam energy onto vulnerable vessel sections. During and immediately after the current quench, excessive radial compression of the runaway beams is avoided by a combination of techniques, improving the likelihood of the beams surviving this dynamic period without a fast termination. Once stabilized, the runaway beams are held in a steady state (out to the ohmic flux limit) with the application of active plasma current and position controls. Beam interaction with the vessel wall is minimized by avoiding distinct thresholds for enhanced wall interaction at small and large radii, corresponding to inner wall and outer limiter interaction, respectively. Staying within the 'safe zone' between those radial thresholds allows for the sustainment of long-lived, quiescent runaway beams. The total beam energy and runaway electron population are then dissipated gradually by a controlled ramp-down of the runaway current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

IOTA experiment for proton pulse compression at extreme space-charge

A gammaT scheme may be required for the PIP-II era performance or ACE-MIRT era performance of the Booster. PIP-II era operations of the Fermilab proton complex will require the Fermilab Booster to increase beam intensity from 4.5e12 to 6.5e12 protons, while also increasing its ramp from 15 Hz to 20 Hz. These changes pose particular challenges for transition-crossing in the Booster, where longitudinal beam quality must be controlled in order to facilitate slip-stacking in the Recycler Ring later in the Main Injector cycle. Two novel gammaT jump schemes are proposed, termed “double gammaT jump” and “partial gammaT jump”, which optimize the magnitude of the gammaT jump within optics and power supply constraints.

Eldred, Jeffrey [Fermilab]↗
DOE OSTI2016

Short-pulse, compressed ion beams at the Neutralized Drift Compression Experiment

We have commenced experiments with intense short pulses of ion beams on the Neutralized Drift Compression Experiment (NDCX-II) at Lawrence Berkeley National Laboratory, with 1-mm beam spot size within 2.5 ns full-width at half maximum. The ion kinetic energy is 1.2 MeV. To enable the short pulse duration and mm-scale focal spot radius, the beam is neutralized in a 1.5-meter-long drift compression section following the last accelerator cell. A short-focal-length solenoid focuses the beam in the presence of the volumetric plasma that is near the target. In the accelerator, the line-charge density increases due to the velocity ramp imparted on the beam bunch. The scientific topics to be explored are warm dense matter, the dynamics of radiation damage in materials, and intense beam and beam-plasma physics including select topics of relevance to the development of heavy-ion drivers for inertial fusion energy. Below the transition to melting, the short beam pulses offer an opportunity to study the multi-scale dynamics of radiation-induced damage in materials with pump-probe experiments, and to stabilize novel metastable phases of materials when short-pulse heating is followed by rapid quenching. First experiments used a lithium ion source; a new plasma-based helium ion source shows much greater charge delivered to the target.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Evaluating microstructure evolution in an SOFC electrode using digital volume correlation

Degradation mechanisms within solid oxide fuel cells (SOFC) during thermal cycling limit operational start-up times and cell lifetime, and must therefore be better understood and mitigated. This work explores such mechanisms using digital volume correlation (DVC) techniques applied to lab-based X-ray tomograms where the microstructural evolution is evaluated during the operational cycling of a Ni–YSZ/YSZ cell. To emulate reduced start-up times, five tomograms were collected over four operational thermal cycles to 750 °C at various ramp-rates: 3, 10, 20 and 30 °C min –1 . Two key features are observed in the 3D strain profiles. Firstly, during low ramp-rates the material produces microscopic channels towards cluster-points. This is thought to be caused by the ceramic skeleton inhibiting the effects of sintering. Secondly, previously unseen macroscopic ‘waves’ developed after high-rate cycling, consisting of linear regions of compression and tension throughout the sample. These wave features decay away from a heterogeneous defect which is thought to be responsible for the non-uniform strain profile. This work demonstrates the first use of sub-micron DVC computations applied to an SOFC exposed to operationally relevant temperatures. Lastly, these findings will assist in the development of new electrode materials from their fabrication to operation, ultimately supporting commercial viability of SOFCs.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

Adaptive mesh based combustion simulations of direct fuel injection effects in a supersonic cavity flame-holder

Here, we present high-fidelity reacting simulations of a supersonic cavity flame-holder configuration. The focus of this work is on flame stabilization brought about by varying the location of fuel injection in a cavity stabilized supersonic flow of air. Central to our approach is a compressible multi-species reacting flow solver that uses adaptive-mesh-refinement (AMR), enabling the resolution of flame, shock-waves, boundary-layers, and small-scale structures in the computational domain. Our analysis indicates that fuel injection closer to the ramp at the aft end of the cavity allows for greater mixing and lower peak temperatures compared to fuel injection upstream that is closer to the backward facing step of the cavity. This difference is mainly due to greater turbulent fluctuations generated from the shear-layer towards the cavity ramp, thereby enhancing the mixing of fuel and air. A low frequency oscillatory behaviour in heat-release and pressure was also observed for the upstream injection case while a much higher-frequency phenomena was observed in the near-ramp injection case. By identifying the important physical determinants of the combustion processes, this study illustrates a promising pathway to design and optimize direct fuel injection strategies in supersonic cavity flame-holders that can improve flame stability, combustion efficiency, and reduce emissions.

42 ENGINEERING↗
DOE OSTI2019

Measurements of the Strain Dependence of Critical Current of Commercial REBCO Tapes at 15 T Between 4.2 and 40 K for High Field Magnets

Interest for high magnetic fields (>16 T) for applications in high energy physics (HEP) and fusion machines, requires the development of high current cables capable to withstand the large forces, mechanical and electromagnetic, experienced during manufacturing and operations. The critical current (I c ) of REBCO tapes depends on strain, magnetic fields, and operational temperatures. Understanding how these parameters affect the I c of the conductor will be critical to develop robust high-current REBCO cables. However, there are limited reports on the strain dependence of I c , in particular at high fields and elevated temperatures relevant for future high-field compact fusion reactor magnets. We present here I c of commercial REBCO tapes measured as a function of compressive and tensile strain (between -0.6% and +0.65%) at high magnetic fields (12 and 15 T) and different temperatures (within 4.2-40 K). Results at 4.2 and 20 K showed less than 5% reduction in the normalized I c at high strain, while a stronger strain dependence was observed at 40 K. Samples tested at 12 T and 4.2 K showed similar strain dependence as 15 T and 4.2 K. In all tested conditions, the tape experienced reversible I c reduction in both tension and compression. Finite element analysis was used to predict the residual thermal strain accumulated in the REBCO layer prior of testing to account for the effect of the cooldown. A method was also developed to account for the current sharing observed between the sample and the sample holder during the ramp of the current. Our results provide useful input for the development of high-field fusion and HEP magnets using REBCO conductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗