Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dense flow”

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 559 records · Page 31

Hall instability driven seeding of helical magneto-Rayleigh–Taylor instabilities in axially premagnetized thin-foil liner Z-pinch implosions

Helical magneto-Rayleigh–Taylor instability (MRTI) structures have been observed in z-pinch-driven liner implosion experiments with a pre-imposed axial magnetic field. We show that the formation of these helical structures can be described by a Hall magnetohydrodynamical (HMHD) model. We used the 3D extended magnetohydrodynamics simulation code PERSEUS (which includes Hall physics) [Seyler and Martin, Phys. Plasmas 18, 012703 (2011)] to study these helical instabilities and show that a Hall interchange instability in low-density coronal plasma immediately surrounding the dense liner is responsible for producing helically oriented effects in the magnetic field and current density within the coronal layer. Here, this seeds the helical pitch angle of the MRTI even when other proposed helical seeding mechanisms are either not present in the experiments or not accounted for in the simulations. For example, this mechanism does not require low-density power-feed plasmas to be swept in from large radius or the development of electrothermal instabilities. The Hall Instability is, thus, a new, independent explanation for the origin of the helical instabilities observed in axially premagnetized liner experiments. Simulation results supporting this mechanism are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Novel Workflow for Coupled Simulation of Hydraulic Stimulation with Simultaneous Injection of Proppant

Enhanced Geothermal Systems (EGS) are located in impermeable host rock. Therefore, to achieve promising circulation, hydraulic stimulation is a key factor. Stimulation involves injecting a combination of fluids and proppant particles into existing fractures under high pressure. The effectiveness of stimulation depends on a variety of different factors, including the pressurization scheme and the optimal choice of carrying fluid and proppant particles to be used. For example, high-density particles can settle out of the injected slurry quickly, whereas lighter particles as well as a dense carrying fluid allow for a long transport into the reservoir. A numerical tool is therefore essential to take all these factors into account and to make robust predictions about the performance of the hydraulic stimulation. We present a novel computational application for coupled numerical simulation of hydraulic stimulation with concurrent injection of proppant in a fractured geothermal reservoir. This application is based on the open-source MOOSE framework that has been enriched with a set of equations to account for the circulation of proppant-laden slurries. It accounts for both the fluid and proppant components of the mixture, as well as particle-driven processes that include gravity settling, particle-particle interactions, and strong density and viscosity contrasts. A particular challenge is the mutual coupling with fracture mechanical processes in a propagating fracture using the FE code Defmod. Aperture and fluid pressure on one side and slurry density and viscosity changes on the other side interact and affect each other. Their basic understanding is therefore essential for the success of the stimulation. We compare our application to common numerical benchmark problems for both slurry flow and hydraulic stimulation.

15 GEOTHERMAL ENERGY↗

Review of the Effects of Polymer Binder Properties on Microstructure and Irreversible Volume Growth of Plastic Bonded Explosives Formulations

The rational design of effective polymeric binders for the formulation of plastic-bonded explosives (PBX) is challenging due to their inherent compositional complexity. The composites comprise irregularly shaped energetic material (EM) powders coated with low weight fractions of polymer via non-equilibrium processes such as slurry coating. Defects can deleteriously affect PBX stability and performance: nano- to micrometer-scale voids can act as loci for hot spots, lowering deflagration and detonation temperatures in unpredictable ways. Furthermore, some nominally desirable polymer properties are at odds with each other: e. g. good flow characteristics are desirable for coating and adhesion, but mechanical stiffness is needed to prevent deformation and cracking of PBX under mechanical stress. Good binder adhesion is critical, but the best means to predict and measure adhesion in PBX is not obvious. Experimental methods of determining binder adhesion on model surfaces may not capture polymer structural configurations relevant to deposition during coating. Molecular dynamics-based computational models have predicted key observables in PBX formulation, suggesting that they may be powerful tools for binder selection. In this review, primarily recent (~2006 and later) literature on polymeric binders for insensitive HE (IHE) is surveyed. We focus on how binder properties influence observable PBX properties as resistance to irreversible volume growth and void formation in PBX formulations mainly (but not exclusively) featuring 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) as the EM (e. g. PBX 9502, LX-17 and others). Conclusions from these studies yield useful guidelines for choosing HE binder candidates, as well as for the general design of highly-filled polymer composite materials. Finally, studies describing challenges in PBX formulation with the newer and more energetically dense high explosive, LLM-105, will be discussed.

adhesion↗

A pushered capsule implosion as an alternate approach to the ignition regime for inertial confinement fusion

We report in inertial confinement fusion, the threshold for ignition is a highly dynamic quantity as the sources and sinks of power in the hot spot can vary rapidly. In this article, we consider the ignition condition as a race between heating and disassembly rates and make use of a prior solution to the fusion hot-spot thermodynamics to develop a Lawson-like ignition criteria for pressure × confinement time (p-τ) vs temperature. Low-Z capsule designs reach the temperature for this threshold using as much of the shell as feasible as ablator but then are limited in τ by low stagnated mass. An alternate approach, the pushered single shell (PSS) design [D. D.-M. Ho, S. MacLaren, and Y. Wang, “High-yield implosions via radiation trapping and high rho-R,” paper presented at the 60th Annual Meeting of the APS Division of Plasma Physics, 2018], introduces a dense inner layer of Mo-Be alloy that is smoothly graded outward to pure Be, increasing the confinement time at stagnation and lowering the temperature requirement at the ignition threshold. Here, we describe a PSS ignition design for the National Ignition Facility and use the theory as well as simulations to compare it with the low-Z capsule approach. Additionally, we show how an adjustment to the design is used to anticipate the effects of mixing at the fuel–ablator interface.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Synthesized spectra of turbulent clouds

A comparison has been made of the velocity structure of non-star-forming dense molecular clouds and the output of a 512(exp 3) hydrodynamic code for compressible gas displaying features of turbulence. A remarkable similarity exists between the interstellar line profiles and the velocity profiles found by integrating, with an optically thin assumption, through subsections of the simulation. This similarity is in terms of both the structure of individual spectra, showing skewness, multiple peaks, and wings, and also the variations of the spectra from point to point. It is possible to see the origin of the various velocity substructures and how they correlate with regions of high density and velocity divergence, or regions of high vorticity, by deprojecting individual spectra of the simulation. It is found that the interstellar medium spectra are best represented by the epoch of the simulation in which the energy spectrum is essentially Kolmogorov-like and initially developed shocks have been largely converted into vorticity structures. Further, it is found that the vorticity peaks ascribed to intermittency contribute to the non-Gaussian features and wings of the line profile, but at an intensity level that is weak compared to the contribution of the bulk of the flow.

Falgarone, E.↗

Plasma in Saturn's Nightside Magnetosphere and the Implications for Global Circulation

We present a bulk ion flow map from the nightside, equatorial region of Saturn's magnetosphere derived from the Cassini CAPS ion mass spectrometer data. The map clearly demonstrates the dominance of corotation flow over radial flow and suggests that the flux tubes sampled are still closed and attached to the planet up to distances of 50RS. The plasma characteristics in the near-midnight region are described and indicate a transition between the region of the magnetosphere containing plasma on closed drift paths and that containing flux tubes which may not complete a full rotation around the planet. Data from the electron spectrometer reveal two plasma states of high and low density. These are attributed either to the sampling of mass-loaded and depleted flux tubes, respectively, or to the latitudinal structure of the plasma sheet. Depleted, returning flux tubes are not, in general, directly observed in the ions, although the electron observations suggest that such a process must take place in order to produce the low-density population. Flux-tube content is conserved below a limit defined by the mass-loading and magnetic field strength and indicates that the flux tubes sampled may survive their passage through the tail. The conditions for mass-release are evaluated using measured densities, angular velocities and magnetic field strength. The results suggest that for the relatively dense ion populations detectable by the ion mass spectrometer (IMS), the condition for flux-tube breakage has not yet been exceeded. However, the low-density regimes observed in the electron data suggest that loaded flux tubes at greater distances do exceed the threshold for mass-loss and subsequently return to the inner magnetosphere significantly depleted of plasma.

Saturn↗

CO2 lidar observations of Mount Pinatubo debris: FIRE 2 and longer-term measurements

The volcanic debris in the stratosphere from the June 1991 eruption of Mt. Pinatubo first appeared over the NOAA Wave Propagation Laboratory (WPL) field site near Boulder, Colorado (40.15 N, 105.23 W), in July of 1991. The presence of the Pinatubo cloud has allowed us to characterize both the tropospheric and stratospheric aerosol backscatter using the NOAA/WPL CO2 Doppler lidar. The lidar has measured vertical backscatter profiles at lambda = 10.59 mu m for over a decade. Analysis of this dense set of profiles reveals the effects of atmospheric and microphysical processes during the buildup and decay of Mt. Pinatubo's clouds. Further information on the NOAA lidar, specifically calibrations using a hard target, can be found in Post and Cupp (1990). We present results of those measurements for June 15, 1991, through December 31, 1992. During that period of longer-term measurements, WPL took part in FIRE II (First ISCCP (International Satellite Cloud Climatology Project) Regional Experiment II), from November 12 through December 8, 1991, measuring vertical backscatter profiles almost daily. One of the mechanisms for purging stratospheric aerosols is tropopause folding, which occurs in cold-core extratropical cyclones. Tropospheric mass loading occurs during folding events which can substantially increase the amount of ice nuclei in the upper troposphere, and may affect the formation of cirrus in that region. Spring and fall are prominent times for tropopause folding events because of the migration of the subtropical and polar jet streams during the transition seasons. Sassen has suggested that the volcanic aerosols from Pinatubo played a role in the formation of cirrus during FIRE II, particularly during a period of moist subtropical flow on December 5-6, 1991.

Levinson, David H.↗

Origin of water masses in Floridan Aquifer System revealed by 81 Kr

Groundwater in coastal aquifers serves as an essential resource in densely populated areas throughout the world. However, hydrological connection to the ocean leaves it at risk of salinization via changes in climate and hydrological cycles. Therefore, an accurate hydrological characterization of coastal aquifers is of the utmost importance. Although localized salinization of shallow aquifers has commonly been studied, geochemical constraints on the time scale of regional freshwater and seawater in coastal aquifers are limited especially in carbonate aquifers. Such information is not only crucial for water resource management, but also for understanding solute flux over the land-ocean boundary and for interpreting the geochemical signatures that each water mass carries. Toward this goal, we report tracer-based, subsurface residence times of groundwater in the Floridan Aquifer System in southern Florida, which is one of the most productive aquifers on Earth. The first application of Kr 81-an emerging tracer of old groundwater-to this aquifer has identified freshwater recharge during the last glacial period as anticipated from previous studies using other geochemical tools including radiocarbon, stable isotopes and noble gas concentrations. Moreover, freshwater of Holocene age was also detected at 100 km-distance from the recharge area, suggesting the possibility of an active flow system in the upgradient region. A contribution of fossil seawater that predates the last glacial maximum was locally identified, suggesting slow seawater circulation and limited but significant solute flux from the ocean into the aquifer for dolomitization. This also implies that the deep saline coastal groundwater can potentially serve as a paleo-seawater archive.

58 GEOSCIENCES↗

Behavior of a Polydisperse Cluster of Drops Evaporating in a Vortex

A model has been developed to describe the dynamics and evaporation of polydisperse collections of liquid drops in an axisymmetric, infinite, cylindrical vortex. This formulation is valid both in the dense regime, where interactions between particles are important, and in the dilute regime, where interactions between particles are not important. In contrast to the standard way of discretizing polydisperse collections of particles, where the initial size distribution is partitioned into a finite number of bins and these bins are fixed during the calculation, here it is only the initial size distribution which is partitioned into a finite number of sizes. Each initial size class thus defined is followed dynamically and thermodynamically in its own system of coordinates which moves with the drops. Therefore, each initial size class develops a continuum of sizes as the drops are centrifuged towards the hotter ambient and evaporate. The gas phase is followed in its own system of coordinates.

liquid evaporation liquid dispersal particle flow↗

Software Testbed for Developing and Evaluating Integrated Autonomous Subsystems

To implement fault tolerant autonomy in future space systems, it will be necessary to integrate planning, adaptive control, and state estimation subsystems. However, integrating these subsystems is difficult, time-consuming, and error-prone. This paper describes Intelliface/ADAPT, a software testbed that helps researchers develop and test alternative strategies for integrating planning, execution, and diagnosis subsystems more quickly and easily. The testbed's architecture, graphical data displays, and implementations of the integrated subsystems support easy plug and play of alternate components to support research and development in fault-tolerant control of autonomous vehicles and operations support systems. Intelliface/ADAPT controls NASA's Advanced Diagnostics and Prognostics Testbed (ADAPT), which comprises batteries, electrical loads (fans, pumps, and lights), relays, circuit breakers, invertors, and sensors. During plan execution, an experimentor can inject faults into the ADAPT testbed by tripping circuit breakers, changing fan speed settings, and closing valves to restrict fluid flow. The diagnostic subsystem, based on NASA's Hybrid Diagnosis Engine (HyDE), detects and isolates these faults to determine the new state of the plant, ADAPT. Intelliface/ADAPT then updates its model of the ADAPT system's resources and determines whether the current plan can be executed using the reduced resources. If not, the planning subsystem generates a new plan that reschedules tasks, reconfigures ADAPT, and reassigns the use of ADAPT resources as needed to work around the fault. The resource model, planning domain model, and planning goals are expressed using NASA's Action Notation Modeling Language (ANML). Parts of the ANML model are generated automatically, and other parts are constructed by hand using the Planning Model Integrated Development Environment, a visual Eclipse-based IDE that accelerates ANML model development. Because native ANML planners are currently under development and not yet sufficiently capable, the ANML model is translated into the New Domain Definition Language (NDDL) and sent to NASA's EUROPA planning system for plan generation. The adaptive controller executes the new plan, using augmented, hierarchical finite state machines to select and sequence actions based on the state of the ADAPT system. Real-time sensor data, commands, and plans are displayed in information-dense arrays of timelines and graphs that zoom and scroll in unison. A dynamic schematic display uses color to show the real-time fault state and utilization of the system components and resources. An execution manager coordinates the activities of the other subsystems. The subsystems are integrated using the Internet Communications Engine (ICE). an object-oriented toolkit for building distributed applications.

Testbed↗

Triphenylphosphine Oxide-Derived Anolyte for Application in Nonaqueous Redox Flow Battery

Recent advances in redox flow batteries have made them a viable option for grid-scale energy storage, however they exhibit low energy density. One way to boost energy density is by increasing the cell potential using a nonaqueous system. Molecular engineering has proven to be an effective strategy to develop redox-active compounds with extreme potentials but these are usually challenged by resource sustainability of the newly developed redox materials. Here, we investigate the utility of phosphine oxides as anolytes with extremely negative potentials. Specifically, we found that cyclic triphenylphosphine oxide (CPO), has a highly negative potential (–2.4 V vs Fc/Fc + ). Importantly, CPO is synthesized from triphenylphosphine oxide, a common industrial chemical waste with no commercial value. Structural and electrochemical characterization of the reduced radical anion showed that enhanced stability is due to cyclization or extended pi-conjugation. Importantly, mechanistic investigation into the decomposition of CPO under various solvents and electrochemical conditions allowed us to utilize an acetonitrile/DMF binary solvent system to enable a long-lived anolyte which exhibited no fade over 350 cycles. In conclusion, this work led to the development of a waste-derived cyclic phosphine oxide that exhibits excellent cycling stability making it an ideal anolyte toward the development of energy-dense RFBs.

25 ENERGY STORAGE↗

Bipolar Membranes with an Electrospun 3D Junction

Freestanding bipolar membranes (BPMs) with an entirely new and transformative morphology were fabricated and characterized. The key disruptive design element was a water-splitting/water-generation junction layer of finite thickness composed of interwoven and interlocking electrospun anion-exchange polymer (AEP) fibers and cation-exchange polymer (CEP) fibers. Methods were developed to sandwich the 3D junction layer between two pre-formed dense AEP and CEP films., with the closure of all interfiber voids, where the processing steps are easily incorporated into a commercial roll-to-roll membrane manufacturing scheme. A series of membranes were made using different polymers and junction layers to identify the optimum morphology and composition for water-splitting and water-generation applications. For water splitting, the 3D junction BPMs worked remarkably well with operating current densities at/above 1.0 A/cm2 (10X greater than commercial BPMs) at a transmembrane voltage drop of only 1.1 V. In water generation mode, the 3D junction BPMs operated at 0.5 A/cm2, a world-record current density. The membranes were found to exhibit outstanding durability and can be manufactured at scale for low cost. The high operating current densities and stable long-term operation are due to the 3D junction layer design, where there is a high interfacial area for water splitting or water generation reactions and where the interlocking fibers prevent delamination of the outer films. These new BPMs are ideal candidate materials for both existing water-splitting electrodialysis separations and new electrochemical processes, such as large-scale direct air capture of CO 2 , reactors for CO 2 reduction, self-hydrating fuel cells, and redox flow batteries.

30 DIRECT ENERGY CONVERSION↗

Thermal and chemical convection in planetary mantles

Melting of the upper mantle and extraction of melt result in the formation of a less dense depleted mantle. This paper describes series of two-dimensional models that investigate the effects of chemical buoyancy induced by these density variations. A tracer particles method has been set up to follow as closely as possible the chemical state of the mantle and to model the chemical buoyant force at each grid point. Each series of models provides the evolution with time of magma production, crustal thickness, surface heat flux, and thermal and chemical state of the mantle. First, models that do not take into account the displacement of plates at the surface of Earth demonstrate that chemical buoyancy has an important effect on the geometry of convection. Then models include horizontal motion of plates 5000 km wide. Recycling of crust is taken into account. For a sufficiently high plate velocity which depends on the thermal Rayleigh number, the cell's size is strongly coupled with the plate's size. Plate motion forces chemically buoyant material to sink into the mantle. Then the positive chemical buoyancy yields upwelling as depleted mantle reaches the interface between the upper and the lower mantle. This process is very efficient in mixing the depleted and undepleted mantle at the scale of the grid spacing since these zones of upwelling disrupt the large convective flow. At low spreading rates, zones of upwelling develop quickly, melting occurs, and the model predicts intraplate volcanism by melting of subducted crust. At fast spreading rates, depleted mantle also favors the formation of these zones of upwelling, but they are not strong enough to yield partial melting. Their rapid displacement toward the ridge contributes to faster large-scale homogenization.

Dupeyrat, L.↗

Gaussian process hydrodynamics

Abstract We present a Gaussian process (GP) approach, called Gaussian process hydrodynamics (GPH) for approximating the solution to the Euler and Navier-Stokes (NS) equations. Similar to smoothed particle hydrodynamics (SPH), GPH is a Lagrangian particle-based approach that involves the tracking of a finite number of particles transported by a flow. However, these particles do not represent mollified particles of matter but carry discrete/partial information about the continuous flow. Closure is achieved by placing a divergence-free GP prior ξ on the velocity field and conditioning it on the vorticity at the particle locations. Known physics (e.g., the Richardson cascade and velocity increment power laws) is incorporated into the GP prior by using physics-informed additive kernels. This is equivalent to expressing ξ as a sum of independent GPs ξ l , which we call modes, acting at different scales (each mode ξ l self-activates to represent the formation of eddies at the corresponding scales). This approach enables a quantitative analysis of the Richardson cascade through the analysis of the activation of these modes, and enables us to analyze coarse-grain turbulence statistically rather than deterministically. Because GPH is formulated by using the vorticity equations, it does not require solving a pressure equation. By enforcing incompressibility and fluid-structure boundary conditions through the selection of a kernel, GPH requires significantly fewer particles than SPH. Because GPH has a natural probabilistic interpretation, the numerical results come with uncertainty estimates, enabling their incorporation into an uncertainty quantification (UQ) pipeline and adding/removing particles (quanta of information) in an adapted manner. The proposed approach is suitable for analysis because it inherits the complexity of state-of-the-art solvers for dense kernel matrices and results in a natural definition of turbulence as information loss. Numerical experiments support the importance of selecting physics-informed kernels and illustrate the major impact of such kernels on the accuracy and stability. Because the proposed approach uses a Bayesian interpretation, it naturally enables data assimilation and predictions and estimations by mixing simulation data and experimental data.

Mathematics↗

Ancient lakes on Mars?

The valley systems in Mars' ancient cratered terrain provide strong evidence for a warmer and wetter climate very early in planetary history. The valley systems in some instances debouch into closed depressions that could have acted as local ponding basins for the flow. A survey of the Martian equatorial region shows that numerous local depressions at the confluence of valley systems exist. These depressions (approximately 100 km) typically are characterized by many valleys flowing into them and few or none flowing out. If ponding did take place, these basin would have contained lakes for some period during Mars' early warmer epoch. Although the collection basins are numerous, location of ones that have not suffered significant subsequent geologic modification is difficult. Some morphologic features suggest that volcanic lavas may have filled them subsequent to any early fluvial activity. Two detailed maps of valley systems and local ponding basins in USGC 1:2,000,000 subquadrangles were completed and a third is in progress. The completed regions are in Mare Tyrrhenum (MC-22 SW) and Margarifter Sinus (MC-19 SE), and the region in progress is in Iapygia (MC-21 NW). On the maps, the valley systems and interpreted margins of ponding basins are indicated. The depressions are of interest for two reasons. First, the depressions were surely the sites in which the materials eroded from the valleys were deposited. Such sediments could preserve important information about the physical conditions at the time of deposition. Second, the sediments could preserve evidence of water-atmosphere interactions during the early period of the Martian climate. Atmospheric carbon dioxide would dissolve in water, and solid carbonate minerals would tend to precipitate out to form carbonate sedimentary deposits. Formation of carbonates in this manner might account for some of the CO2 lost from the early more dense atmosphere.

Goldspiel, J. M.↗

Plasmasphere dynamics in the duskside bulge region: A new look at old topic

Data acquired during several multiday periods in 1982 at ground stations Siple, Halley, and Kerguelen and on satellites Dynamics Explorer 1, International Sun Earth Explorer 1, and GEOS 2 have been used to investigate thermal plasma structure and dynamics in the duskside plasmasphere bulge region of the Earth. The distribution of thermal plasma in the dusk bulge sector is difficult to describe realistically, in part because of the time integral manner in which the thermal plasma distribution depends upon on the effects of bulk cross-B flow and interchange plasma flows along B. While relatively simple MHD models can be useful for qualitatively predicting certain effects of enhanced convection on a quiet plasmasphere, such as an initial sunward entrainment of the outer regions, they are of limited value in predicting the duskside thermal plasma structures that are observed. Furthermore, use of such models can be misleading if one fails to realize that they do not address the question of the formation of the steep plasmapause profile or provide for a possible role of instabilities or other irreversible processes in plasmapause formation. Our specific findings, which are based both upon the present case studies and upon earlier work, include the following: (1) during active periods the plasmasphere appears to become divided into two entities, a main plasmasphere and a duskside bulge region. (2) in the aftermath of an increase in convection activity, the main plasmasphere tends (from a statistical point of view) to become roughly circular in equatorial cross section, with only a slight bulge at dusk; (3) the abrupt westward edge of the duskside bulge observed from whistlers represents a state in the evolution of sunward extending streamers; (4) in the aftermath of a weak magnetic storm, 10 to 30% of the plasma 'removed' from the outer plasmasphere appears to remain in the afternoon-dusk sector beyond the main plasmasphere. (5) outlying dense plasma structures may circulate in the outer duskside magetosphere for many days following an increase in convection, unless there is extremely deep quieting; (6) a day-night plasmatrough boundary may be identified in equatorial satellite data; (7) factor-of-2-to-10 density irregularities appear near the plasmatrough from the ionosphere at L = 4.6, predominantly bidirectional field aligned and equatorially trapped light ion pitch angle distributions give away to a predominantly isotropic distribution (as seen by DE 1) when the plasma density reaches a level a factor of about 3 below the satured plasmasphere level; (9) some outlying dense plasma structures are effectively detached from the main plasmasphere, while others appear to be connected to that body.

Carpenter, D. L.↗

Techniques for studying materials under extreme states of high energy density compression

The properties of materials under extreme conditions of pressure and density are of key interest to a number of fields, including planetary geophysics, materials science, and inertial confinement fusion. In geophysics, the equations of state of planetary materials, such as hydrogen and iron, under ultrahigh pressure and density provide a better understanding of their formation and interior structure [Celliers et al., “Insulator-metal transition in dense fluid deuterium,” Science 361, 677–682 (2018) and Smith et al., “Equation of state of iron under core conditions of large rocky exoplanets,” Nat. Astron. 2, 591–682 (2018)]. The processes of interest in these fields occur under conditions of high pressure (100 GPa–100 TPa), high temperature (>3000 K), and sometimes at high strain rates (>103 s−1) depending on the process. With the advent of high energy density (HED) facilities, such as the National Ignition Facility (NIF), Linear Coherent Light Source, Omega Laser Facility, and Z, these conditions are reachable and numerous experimental platforms have been developed. To measure compression under ultrahigh pressure, stepped targets are ramp-compressed and the sound velocity, measured by the velocity interferometer system for any reflector diagnostic technique, from which the stress-density of relevant materials is deduced at pulsed power [M. D. Knudson and M. P. Desjarlais, “High-precision shock wave measurements of deuterium: Evaluation of exchange-correlation functionals at the molecular-to-atomic transition,” Phys. Rev. Lett. 118, 035501 (2017)] and laser [Smith et al., “Equation of state of iron under core conditions of large rocky exoplanets,” Nat. Astron. 2, 591–682 (2018)] facilities. To measure strength under high pressure and strain rates, experimenters measure the growth of Rayleigh–Taylor instabilities using face-on radiography [Park et al., “Grain-size-independent plastic flow at ultrahigh pressures and strain rates,” Phys. Rev. Lett. 114, 065502 (2015)]. The crystal structure of materials under high compression is measured by dynamic x-ray diffraction [Rygg et al., “X-ray diffraction at the national ignition facility,” Rev. Sci. Instrum. 91, 043902 (2020) and McBride et al., “Phase transition lowering in dynamically compressed silicon,” Nat. Phys. 15, 89–94 (2019)]. Medium range material temperatures (a few thousand degrees) can be measured by extended x-ray absorption fine structure techniques, Yaakobi et al., “Extended x-ray absorption fine structure measurements of laser-shocked V and Ti and crystal phase transformation in Ti,” Phys. Rev. Lett. 92, 095504 (2004) and Ping et al., “Solid iron compressed up to 560 GPa,” Phys. Rev. Lett. 111, 065501 (2013), whereas more extreme temperatures are measured using x-ray Thomson scattering or pyrometry. This manuscript will review the scientific motivations, experimental techniques, and the regimes that can be probed for the study of materials under extreme HED conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hydrogen Refueling Reference Station Lot Size Analysis for Urban Sites

Hydrogen Fueling Infrastructure Research and Station Technology (H2FIRST) is a project initiated by the DOE in 2015 and executed by Sandia National Laboratories and the National Renewable Energy Laboratory to address R&D barriers to the deployment of hydrogen fueling infrastructure. One key barrier to the deployment of fueling stations is the land area they require (i.e. "footprint"). Space is particularly a constraint in dense urban areas where hydrogen demand is high but space for fueling stations is limited. This work presents current fire code requirements that inform station footprint, then identifies and quantifies opportunities to reduce footprint without altering the safety profile of fueling stations. Opportunities analyzed include potential new methods of hydrogen delivery, as well as alternative placements of station technologies (i.e. rooftop/underground fuel storage). As interest in heavy-duty fueling stations and other markets for hydrogen grows, this study can inform techniques to reduce the footprint of heavy-duty stations as well. This work characterizes generic designs for stations with a capacity of 600 kg/day hydrogen dispensed and 4 dispenser hoses. Three base case designs (delivered gas, delivered liquid, and on-site electrolysis production) have been modified in 5 different ways to study the impacts of recently released fire code changes, colocation with gasoline refueling, alternate delivery assumptions, underground storage of hydrogen, and rooftop storage of hydrogen, resulting in a total of 32 different station designs. The footprints of the base case stations range from 13,000 to 21,000 ft 2 . A significant focus of this study is the NFPA 2 requirements, especially the prescribed setback distances for bulk gaseous or liquid hydrogen storage. While the prescribed distances are large in some cases, these setback distances are found to have a nuanced impact on station lot size; considerations of the delivery truck path, traffic flow, parking, and convenience store location are also important. Station designs that utilize underground and rooftop storage can reduce footprint but may not be practical or economical. For example, burying hydrogen storage tanks underground can reduce footprint, but the cost savings they enable depend on the cost of burial and the cost land. Siting and economic analysis of station lot sizes illustrate the benefit of smaller station footprints in the flexibility and cost savings they can provide. This study can be used as a reference that provides examples of the key design differences that fueling stations can incorporate, the approximate sizes of generic station lots, and considerations that might be unique to particular designs.

08 HYDROGEN↗