Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “compression rate”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

At Extreme Strain Rates, Pure Metals Thermally Harden while Alloys Thermally Soften

When materials are deformed at extreme strain rates, >10 6 s -1 , a counterintuitive mechanical response is seen where the strength and hardness of pure metals increases with increasing temperature. This anti-thermal hardening is due to dislocations meeting resistance to their motion from phonons in the crystal lattice. However, here, using optically-driven microballistic impact testing to measure dynamic strength and hardness, we show that when the composition is systematically varied away from high purity, the mechanical response of metals transitions from phonon drag of dislocations back to thermally activated pinning of dislocations, even at the highest strain rates. This boundary from “hotter-is-stronger” to “hotter-is-softer” is observed and mapped for nickel, titanium and gold. Furthermore, the ability to tune between deformation mechanisms with very different temperature dependencies speaks to new directions for alloy design in extreme conditions.

Ballistic transport↗
DOE OSTI2022

Ambient-temperature liquid jet targets for high-repetition-rate HED discovery science

High-power lasers can generate energetic particle beams and astrophysically relevant pressure and temperature states in the high-energy-density (HED) regime. Recently-commissioned high-repetition-rate (HRR) laser drivers are capable of producing these conditions at rates exceeding 1 Hz. However, experimental output from these systems is often limited by the difficulty of designing targets that match these repetition rates. To overcome this challenge, we have developed tungsten microfluidic nozzles, which produce a continuously replenishing jet that operates at flow speeds of approximately 10 m/s and can sustain shot frequencies up to 1 kHz. The ambient-temperature planar liquid jets produced by these nozzles can have thicknesses ranging from hundreds of nanometers to tens of micrometers. In this work, we illustrate the operational principle of the microfluidic nozzle and describe its implementation in a vacuum environment. Further, we provide evidence of successful laser-driven ion acceleration using this target and discuss the prospect of optimizing the ion acceleration performance through an in situ jet thickness scan. Future applications for the jet throughout HED science include shock compression and studies of strongly heated nonequilibrium plasmas. When fielded in concert with HRR-compatible laser, diagnostic, and active feedback technology, this target will facilitate advanced automated studies in HRR HED science, including machine learning-based optimization and high-dimensional statistical analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Dynamic compression and spallation behavior of NbTaTiVZr high-entropy alloy

High-entropy alloys are a new class of materials with promising properties for aerospace and defense applications. The unique behavior of these materials is driven by the complex interactions between dissimilar atoms in a crystal and requires combined theoretical and experimental efforts to unlock their full potential. Here, we evaluate the relationship between the microstructure and the dynamic response of the equiatomic NbTaTiVZr alloy. Specifically, the shock Hugoniot, or equation of state, was measured up to a particle velocity of 0.5 mm μs −1 using gas gun plate impacts. Shock wave profile and incipient spallation experiments were used to characterize wave propagation and damage formation with post-mortem recovery experiments. Molecular dynamics simulations confirm the experimental findings and extend them up to 2.1 mm μs −1 particle velocity. Computational thermodynamic calculation of phase diagram simulations explain details of the material microstructure, which explains the measured strength and experimental damage patterns. Overall, this work provides detailed high-strain-rate characterization of a refractory high-entropy alloy, and more importantly, a framework and demonstration of the utility and necessity of a combined theoretical and experimental approach, outlining the importance of considering processing and manufacturing conditions when evaluating the performance of new materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

On the bulk compaction of brittle granular materials, part I: SeS analysis of axial compression to 4000 MPa

The bulk compaction of granular materials has been studied for decades to interpret and manage responses for soils and powder-based component fabrication, and geophysical, celestial, and ballistic impact. Their bulk or macroscopic compaction response is limited by what occurs at the granular or microstructural scale. Motivation existed to more closely examine that association specific to granular brittle materials (e.g., ceramics and glasses). That examination is offered in a series of three companion papers where Part I describes a new supplemental analysis adopted to bulk compaction response involving relatively high compaction stresses (4000 MPa). Bulk compactions of vitreous silicates and crystalline quartzes were interpreted in three ways, including that of a new analysis that considers the product of void ratio (e) and stress (S) as a function of S, hereafter referred to as “SeS analysis”. In conclusion, the SeS analysis was found to be an informative supplement to conventional bulk compaction analyses because it provides more consistent higher sensitivity for the identification of bulk density rate increase with increasing compaction (softening); a rate increase that arises from the cumulative effect of the onsets and progression of compaction-induced yielding, fracture or comminution, densification, phase change, or combinations thereof occurring at the granular or microstructural scale.

compaction↗
DOE OSTI2021

Constitutive Model Development for Aging Polymer Encapsulants (ASC P&EM FY2021 L2 Milestone 7836)

This SAND report fulfills the completion requirements for the ASC Physics and Engineering Modeling Level 2 Milestone 7836 during Fiscal Year 2021. The Sandia Simplified potential energy clock (SPEC) non-linear viscoelastic constitutive model was developed to predict a whole host of polymer glass physical behaviors in order to provide a tool to assess the effects of stress on these materials over their lifecycle. Polymer glasses are used extensively in applications such as electronics packaging, where encapsulants and adhesives can be critical to device performance. In this work, the focus is on assessing the performance of the model in predicting material evolution associated with long-term physical aging, an area that the model has not been fully vetted in. These predictions are key to utilizing models to help demonstrate electronics packaging component reliability over decades long service lives, a task that is very costly and time consuming to execute experimentally. The initiating hypothesis for the work was that a model calibration process can be defined that enables confidence in physical aging predictions under ND relevant environments and timescales without sacrificing other predictive capabilities. To test the hypothesis, an extensive suite of calibration and aging data was assembled from a combination of prior work and collaborating projects (Aging and Lifetimes as well as the DoD Joint Munitions Program) for two mission relevant epoxy encapsulants, 828DGEBA/DEA and 828DGEBA/T403. Multiple model calibration processes were developed and evaluated against the entire set of data for each material. A qualitative assessment of each calibration's ability to predict the wide range of aging responses was key to ranking the calibrations against each other. During this evaluation, predictions that were identified as non-physical, i.e., demonstrated something that was qualitatively different than known material behavior, were heavily weighted against the calibration performance. Thus, unphysical predictions for one aspect of aging response could generate a lower overall rating for a calibration process even if that process generated better quantitative predictions for another aspect of aging response. This insurance that all predictions are qualitatively correct is important to the overall aim of utilizing the model to predict residual stress evolution, which will depend on the interplay amongst the different material aging responses. The DSC-focused calibration procedure generated the best all-around aging predictions for both materials, demonstrating material models that can qualitatively predict the whole host of different physical aging responses that have been measured. This step forward in predictive capability comes from an unanticipated source, utilization of calorimetry measurements to specify model parameters. The DSC-focused calibration technique performed better than compression-focused techniques that more heavily weigh measurements more closely related to the structural responses to be predicted. Indeed, the DSC-focused calibration procedure was only possible due to recent incorporation of the enthalpy and heat capacity features into SPEC that was newly verified during this L2 milestone. Fundamentally similar aspects of the two material model calibrations as well as parametric studies to assess sensitives of the aging predictions are discussed within the report. A perspective on the next steps to the overall goal of residual stress evolution predictions under stockpile conditions closes the report.

36 MATERIALS SCIENCE↗
DOE OSTI2023

An Investigation of Nanoindentation as Quality Control [Slides]

Nanoindentation Continuous Stiffness Measurement (CSM) mapping, and Strain Rate Jump Testing (SRJT) proved effective to highlight and visualize the extent of property differences between the samples. From the analysis, Lot A has a significant dead zone at its center, especially when compared with the cross section of Lot B. Quick identification of property issues from one of the two feed stocks, both of which fit the specification. Mechanical performance was connected to the microstructure. Thus, using the structure-property results, it is noted that the processing of Lot A may have been insufficient. Nanoindentation allows for the sampling of quality control variables such as hardness and modulus in a short time span at high resolution, which is cheaper and faster than traditional strength experiments (tension and compression tests).

36 MATERIALS SCIENCE↗
DOE OSTI2022

TEA of a Unique Two-Pathways Process for Post-Combustion CO 2 Capture

A unique two-Pathways process using aqueous sodium glycinate for CO 2 capture from a split flue gas stream emitted from 600 MWe post-combustion coal power plant was developed in Aspen Plus v.10. The split gas flow rate used was 44.75 ton/h and contained 0.0023 mol% SO 2 and 13.33 mol% CO 2 . The process includes a washing unit, a CO 2 absorption unit, a reverse osmosis unit, and a solvent regeneration unit or an ultrafiltration unit. The washing unit uses deionized water to completely remove SO 2 and the CO 2 absorption unit uses SGS to capture at least 90 mol% of the CO 2 in the split flue gas stream. Upon CO 2 and SGS reactions, the resulting liquid products exhibit phase-separation into CO 2 -lean phase and CO 2 -rich phase, allow two distinct pathways. Pathway (i) is to regenerate mostly the CO 2 -rich phase, collect the released CO 2 , and compress it for sequestration purposes. Pathway (ii) is to send the liquid stream from the CO 2 absorption unit to the ultrafiltration unit to separate the solid nanomaterials. The hydraulics and mass transfer characteristics in the washing and CO 2 absorption units were obtained; and techno-economic analysis (TEA) for Pathways (i) and (ii), including Capital Expenditure (CAPEX), Operating Expenditure (OPEX), and Levelized Cost of CO 2 Captured (LCOC), were calculated and compared. The simulation results revealed that the CAPEX, OPEX, and LCOC for Pathway (i) were ($\$12,039,251$), (261 dollar/h), and (54.01 dollar/ton-CO 2 captured), respectively, and those for Pathway (ii) were ($\$5,908,000$), (237.2 dollar/h), and (39.90 dollar/ton-CO 2 captured), respectively. Moreover, in Pathway (ii), 8.19 ton/h of CO 2 were captured to produce 15.62 ton/h NaHCO 3 nanomaterials, which were sold to offset the overall process cost. The LCOC values indicate that Pathway (ii) is more cost-effective than Pathway (i) because LCOC values for Pathway (ii) are much lower than those for Pathway (i).COC values for Pathway (ii) are much lower than those for Pathway (i).

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2025

Numerical analysis of air dehumidification through electrospray-enhanced vortical flow cyclone separator

The dehumidification of air is a highly energy-intensive process, typically requiring vapor compression systems to cool the air significantly below its dew point temperature. Here, this study numerically investigates air dehumidification achieved by electrostatic and dielectrophoresis forces in a converging duct and centrifugal principles in cyclone separators. An electrospray was employed to energize water droplets, which then captured water vapor via dielectrophoresis. A novelty of the present study is the inclusion of a cyclone separator positioned downstream of the electrospray. This combination increased the dehumidification efficiency and controlled the flow of the condensate. By integrating electrospray technology and applying an electric field, this study examines the impact of particle separation efficiency, pressure drops, and cyclone performance on air dehumidification. The control parameters for air dehumidification were investigated. Three cyclones were modelled for airflow rates of 5 cubic feet per minute [cfm] (0.0024 m 3 /s), 38 cfm (0.0179 m 3 /s), and 200 cfm (0.094 m 3 /s). The results indicated air dehumidification was up to 5.7 % under environmental conditions, similar to those used to validate the CFD model. The degree of air dehumidification was primarily influenced by the electric charge acquired by the water droplets, their flight trajectory, and the vorticity intensity within the cyclone. Air dehumidification could increase by up to 13 % by increasing the number of droplets injected into the air stream, the flight time, and the electrical charge. The breakdown electric field strength threshold near the electrospray was a limiting factor for further enhancing dehumidification performance. This study examines the optimization of an electro-assisted dehumidification system that includes cyclone-based droplet separation through numerical analysis. The results provide insights into the design of energy-efficient air treatment technologies.

42 ENGINEERING↗
DOE OSTI2024

First-principles investigation of the thermophysical properties of NaCl, PuCl 3 , and NaCl-PuCl 3 Molten salts

Molten salts have a variety of applications that span the nuclear and solar industries, and which involve thermal storage and heat transfer. There is a present knowledge gap in the thermophysical properties of molten salts, which limits the readiness level of molten salt applications. This is especially pertinent for molten salt reactors where the fissile material is dissolved within the molten salt. Ab initio molecular dynamics is a common method employed to investigate the structural and thermophysical properties at elevated temperatures. NaCl-PuCl 3 is a candidate fuel salt in molten salt reactors. The scope of this study is to investigate the seven unique compositions and a range of temperatures of NaCl-PuCl 3 and calculate the density, heat capacity, compressibility, enthalpy of mixing, and coefficient of thermal expansion. Within this work, the van der Waals (vdW) interactions are primarily handled with the vdW-DF2 functional but spot-checked with both the dDsC and DFT-D3 methods. The calculated densities are in good agreement with the NaCl and the eutectic compositions, and no experimental values exist for PuCl 3 at these temperatures. The densities are fit to a first-order Redlich-Kister expansion as a function of both composition and temperature. The heat capacity in the literature is scattered but the calculated values agree well with one of the experimental results. Here, the heat capacity increases as a linear function with respect to concentration at a rate of about 7.5 J/mol-K per 10 mol% PuCl 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2025

Century: Zap Energy’s 100-kW-Scale Repetitive Sheared-Flow-Stabilized Z -Pinch System with Liquid Metal Cooling

Zap Energy is developing the sheared-flow-stabilized (SFS) Z-pinch concept for commercial applications. The SFS Z pinch relies on plasma self-organization, in the sense that plasma dynamics play a critical role in confinement. Using plasma axial current for confinement and compression eliminates the need for external confinement or heating technologies. This compact magnetic confinement technology could, in turn, provide the basis for a cost-effective deuterium-tritium fusion power plant. In addition to a robust experimental program pushing plasma performance towards breakeven conditions, Zap Energy has parallel programs developing power handling systems suitable for future power plants. Technologies under development include high average-power repetitive pulsed power, high duty-cycle cathodes, and liquid metal wall systems. Century is the name of Zap Energy’s first effort to integrate these three components into an operational system capable of firing non-reacting hydrogen SFS Z-pinch plasmas into a liquid-metal-lined container at sustained repetition rates on the order of 0.1 Hz. Here, the pulsed power driver and liquid metal heat exchanger are both designed to sustain input powers of 100 kW. Construction and initial operations with an interim ~10 kW liquid metal heat exchanger are described.

Century↗
DOE OSTI2023

A numerical study of an atomizing jet in a resonant acoustic field

By applying state-of-the-art, high-fidelity methods for compressible multiphase flows, this work features a parametric study of a turbulent liquid jet atomizing under the influence of standing acoustic waves with five different amplitude values. We perform these simulations in a semi-periodic domain to isolate interactions between acoustic and atomization processes and to emulate a pure liquid fuel spray placed at the pressure node of a resonant acoustic field. Characterization of the flow includes measurements of the mean dimensions, droplet number, and surface area of the liquid distribution, alongside qualitative depictions of the interface evolution, breakup events, and hydrodynamic fields. As acoustic forcing is provided to the flow, the added energy provokes atomization when shear would otherwise be insufficient for breakup. The magnitude of the forcing determines the rate of droplet production, which accelerates as the evolving liquid surface modifies the induced motion of the gas phase. At the nodal plane, the acoustic field produces an organizing effect on the liquid structures.

42 ENGINEERING↗
DOE OSTI2022

Modeling Sequences of Earthquakes and Aseismic Slip (SEAS) in Elasto-Plastic Fault Zones With a Hybrid Finite Element Spectral Boundary Integral Scheme

We present a coupled finite element spectral boundary integral framework for modeling sequences of earthquakes and aseismic slip on a 2-D planar rate-and-state fault with off-fault visco-plastic response in the plane strain approximation. The model resolves both slow aseismic deformation and inertia effects during rapid slip. As an application, we perform two sets of simulations with different choices of cohesion to explore the co-evolution of fault slip, bulk plasticity and local stress fields. The first set implements a relatively large value of the cohesion parameter, which results in limiting inelastic strain accumulation to dynamic rupture phases. The second set implements a smaller cohesion, allowing for plastic strain to accumulate in both seismic and aseismic phases. For the first model, our results indicate that the extent and distribution of plastic strain depend on the angle of maximum compressive principal stress. At larger angles, inelastic strain accumulates on the extensional side of a dynamically propagating rupture. At smaller angles, the extent of plasticity is limited to the compressional side of the domain. At smaller cohesion values, off-fault plasticity may occur during aseismic slip, which alters the nucleation characteristics and earthquake sequence pattern. Furthermore, our results at lower cohesion values indicate that plastic strain accumulation may occur in both the extensional and compressional sides of the off-fault bulk even at higher angles of maximum compression. In conclusion, this produces damage patterns that deviate from the traditional off-fault fan-like distribution observed in dynamic rupture simulations and emphasizes the significance of long-term deformation in interpreting observations.

58 GEOSCIENCES↗
DOE OSTI2020

Geometric and Material Variability Influences Stress States Relevant to Coastal Permafrost Bluff Failure

Scientific knowledge and engineering tools for predicting coastal erosion are largely confined to temperate climate zones that are dominated by non-cohesive sediments. The pattern of erosion exhibited by the ice-bonded permafrost bluffs in Arctic Alaska, however, is not well-explained by these tools. Investigation of the oceanographic, thermal, and mechanical processes that are relevant to permafrost bluff failure along Arctic coastlines is needed. We conducted physics-based numerical simulations of mechanical response that focus on the impact of geometric and material variability on permafrost bluff stress states for a coastal setting in Arctic Alaska that is prone to toppling mode block failure. Our three-dimensional geomechanical boundary-value problems output static realizations of compressive and tensile stresses. We use these results to quantify variability in the loci of potential instability. We observe that niche dimension affects the location and magnitude of the simulated maximum tensile stress more strongly than the bluff height, ice wedge polygon size, ice wedge geometry, bulk density, Young's Modulus, and Poisson's Ratio. Our simulations indicate that variations in niche dimension can produce radically different potential failure areas and that even relatively shallow vertical cracks can concentrate displacement within ice-bonded permafrost bluffs. These findings suggest that stability assessment approaches, for which the geometry of the failure plane is delineated a priori, may not be ideal for coastlines similar to our study area and could hamper predictions of erosion rates and nearshore sediment/biogeochemical loading.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Constitutive behavior and processing maps of a new wrought magnesium alloy ZE20 (Mg-2Zn-0.2Ce)

ZE20 (Mg-2Zn-0.2Ce) 2 is a new wrought magnesium alloy with improved extrudability and mechanical properties [1]. To understand the constitutive behavior and workability of this new alloy, Gleeble thermomechanical testing has been carried out in this study. The flow stress behavior of ZE20 was investigated between 250 °C–450 °C and 10 –3 s –1 –1.0 s –1 in isothermal compression. Constitutive descriptions of the flow stress are provided. A new general approach at application of the extended Ludwik equation is demonstrated and was found to be more accurate than the hyperbolic sine Arrhenius model while having a similar number of model constants. Processing maps were developed based on the experimental results and are verified with microstructural investigation. A region of safe processing with non-basal texture and high activity of dynamic recrystallization (DRX) was found between 375 °C and 450 °C, from 10 –1 s –1 to 10–2.5 s –1 . A region of potentially safe processing with annealing that is associated with shear band nucleation of non-basal grains was identified for temperatures as low as 300 °C and rates as high as 10 –1 s –1 .

36 MATERIALS SCIENCE↗
DOE OSTI2022

Time-dependent THMC properties and microstructural evolution of damaged rocks in excavation damage zone

Modeling coupled thermo-hydro-mechanical-chemical (THMC) processes in host rocks near high-level nuclear waste (HLW) repositories at various time scales is an extremely challenging task. The current study integrates experimental, theoretical, and numerical methods in assessing the evolution of excavation damage zone (EDZ) over time and its implication on the long-term migration of hazardous species. Argillite and rock salt and are the focus of this study. The first part of the report presents a novel time-dependent directional microcrack damage theory for generic brittle rocks. It features detailed statistical description of the microcracks within a damaged solid, permitting a direct upscaling of microscale processes such as crack growth kinetics, crack closure/opening, sliding friction to explain the macroscopic creep, nonlinear elasticity, shear dilation, and loading-unloading hysteresis. This provides a basic platform for describing the anisotropic mechanical and transport properties of damaged rocks during excavation and subsequent THMC loadings. The model is validated and numerically implemented to Finite Element (FE) package ABAQUS through the user-defined material (UMAT) interface and have demonstrated great potential in resolving the time-dependent and anisotropic evolution of damage in the EDZ. The second part of the report presents a multi-scale experimental effort in characterizing the thermal, hydraulic, and mechanical properties of Mancos shale and Avery Island salt. For the Mancos shale, triaxial compression tests are performed at different confining pressures and temperatures to probe its thermomechanical properties relevant to HLW repositories. The obtained stress-strain data are interpreted using the proposed directional damage theory. Post-test specimens are subjected to gas permeability tests to reveal the correlation between permeability and the degree of microcracking. At microscale, temperature-controlled nanoindentation tests are performed and found a linear correlation between fracture toughness and elastic modulus from 25°C to 300°C. For Avery Island salt, we have designed and manufactured a novel relative-humidity controlled uniaxial creep device. Long-term creep tests at low stresses (< 5 MPa) are performed at different levels of relative humidity (RH). Besides confirming the much higher creep rates as one would expect through extrapolating the high-stress creep data, the results reveal that the steady-state creep rate of rock salt is strongly dependent on the ambient RH, an aspect that is often neglected in the literature. Both behaviors can be attributed to the pressure-solution creep mechanism which dominates at low stress and high RH levels. The third part of the report explores a set of numerical strategies in modeling the THMC behavior of porous geomaterials. A fully implicit, monolithic FE solution that can flexibly interface with different material models and coupling mechanisms for THMC problems is developed and verified through the ABAQUS user-defined element (UEL) interface. The scheme is used to study the THM response of a hypothetical HLW storage site with reference to an existing in-situ heater test. Strategies for integrating the UEL and the microcrack UMAT are suggested. Another numerical endeavor of this study is to implement a higher-order asymptotic homogenization method to account for the heterogeneous porous structures. The same method is then extended to perform microstructure-informed thermo-mechanical modeling of generalized continua. The above outcomes of this project provide a strong thrust towards enhancing the fundamental understanding and modeling capability of the long-term evolution of host rocks in EDZ, thus helping achieve the design goal of 1-million-year isolation of high-level nuclear wastes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Residual fracture energy of natural and recycled aggregate concrete after exposure to high temperatures

Abstract Fracture energy ( G F ) represents the average total energy consumed in the whole process of crack propagation. Using one parameter of G F , it is not sufficient to describe the fracture behavior during crack stable and unstable propagation periods. In this study, the G F , stable fracture energy ( G FS ) and unstable fracture energy ( G FU ) of natural and recycled aggregate concrete (RAC) were investigated by using three‐point bending tests. Totally, 128 cubic specimens (100 mm × 100 mm × 100 mm) and 64 notched beams with size of 510 mm × 100 mm × 100 mm were manufactured to test the compressive strength ( f cu ), splitting tensile strength ( f t ), load–deflection (P‐δ), and load‐cracking mouth opening displacement (P‐CMOD) curves, considering different replacement levels (0%, 30%, 70%, and 100%) of recycled coarse aggregate (RA) and high temperatures exposure (20, 100, 200, and 300°C). The effects of RA substitutions and temperatures on basic mechanical properties, G F , G FS , and G FU were analyzed. The results demonstrated that the f cu , f t , and elastic modulus ( E ) decreased with the incremental temperature. The G FS of concrete with different RA substitutions first decreased and then increased, and the minimum value appeared at 200°C. The G F and G FU of concrete with different RA substitutions increased with the elevated temperature, and the increase rate decreased with RA substitutions. When the temperature is less than 100°C, the G F and G FU of concrete remained mostly unchanged with the increase of RA substitutions, while the temperature reached 200 and 300°C, the G F and G FU of concrete increased and decreased with the increase of RA substitutions, respectively.

Chen, Wei↗
DOE OSTI2020

Dispersion of cryogenic hydrogen through high-aspect ratio nozzles

Liquid hydrogen is increasingly being used as a delivery and storage medium for stations that provide compressed gaseous hydrogen for fuel cell electric vehicles. In efforts to provide scientific justification for separation distances for liquid hydrogen infrastructure in fire codes, the dispersion characteristics of cryogenic hydrogen jets (50–64 K) from high aspect ratio nozzles have been measured at 3 and 5 bar abs stagnation pressures. These nozzles are more characteristic of unintended leaks, which would be expected to be cracks, rather than conventional round nozzles. Spontaneous Raman scattering was used to measure the concentration and temperature field along the major and minor axes. Within the field of interrogation, the axis-switching phenomena was not observed, but rather a self-similar Gaussian-profile flow regime similar to room temperature or cryogenic hydrogen releases through round nozzles. The concentration decay rate and half-widths for the planar cryogenic jets were found to be nominally equivalent to that of round nozzle cryogenic hydrogen jets indicating a similar flammable envelope. The results from these experiments will be used to validate models for cryogenic hydrogen dispersion that will be used for simulations of alternative scenarios and quantitative risk assessment.

08 HYDROGEN↗
DOE OSTI2023

Molecular Rotations, Multiscale Order, Hyperuniformity, and Signatures of Metastability during the Compression/Decompression Cycles of Amorphous Ices

We model, via large-scale molecular dynamics simulations, the isothermal compression of low-density amorphous ice (LDA) to generate high-density amorphous ice (HDA) and the corresponding decompression extending to negative pressures to recover the low-density amorphous phase (LDA HDA ). Both LDA and HDA are nearly hyperuniform and are characterized by a dynamical HBN, showing that amorphous ices are nonstatic materials and implying that nearly hyperuniformity can be accommodated in dynamical networks. In correspondence with both the LDA-to-HDA and the HDA-to-LDA HDA phase transitions, the (partial) activation of rotational degrees of freedom activates a cascade effect that induces a drastic change in the connectivity and a pervasive reorganization of the HBN topology which, ultimately, break the samples’ hyperuniform character. Key to this effect is the rapid rate at which changes occur, and not their magnitude. The inspection of structural properties from the short- to the long-range shows that signatures of metastability are present at all length-scales, hence providing further solid evidence in support of the liquid–liquid critical point scenario. LDA and LDA HDA differ in terms of HBN and structural properties, implying that they are distinct low-density glasses. Our work unveils the role of molecular rotations in the phase transitions between amorphous ices and shows how the unfreezing of rotational degrees of freedom generates a cascade effect that propagates over multiple length-scales. Our findings greatly improve our basic understanding of water and amorphous ices and can potentially impact the field of molecular network-forming materials at large.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗