SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
The physics governing the response of materials to shock loading has many applications in various fields, ranging from Earth and planetary science to engineering and ballistics. Elucidating the behavior of materials experiencing the passage of shock waves motivates active research worldwide, including experimental developments with in situ and time-resolved measurements, modern computational approaches, and multi-scale simulations.
High-energy-density laser facilities and advances in dynamic compression techniques have expanded access to material states in the Terapascal regime relevant to inertial confinement fusion, planetary science, and geophysics. However, experimentally determining the material temperature in these extreme conditions has remained a difficult challenge. Extended X-ray Absorption Fine Structure (EXAFS), referring to the modulations in x-ray absorption above an absorption edge from photoelectrons’ interactions with neighboring atoms, has proven to be a versatile and robust technique for probing material temperature and density for mid-to-high Z elements under dynamic compression. The current platform at the National Ignition Facility has developed six configurations for EXAFS measurements between 7 and 18 keV for different absorption edges (Fe K, Co K, Cu K, Ta L3, Pb L3, and Zr K) using a curved-crystal spectrometer and a bright, continuum foil x-ray source. Here, in this work, we describe the platform geometry, x-ray source performance, spectrometer resolution and throughput, design considerations, and data in ambient and dynamic-compression conditions.
Shock-compressing material samples that are precompressed with diamond-anvil cells (DACs) enables benchmarking equation of state measurements over a wide range of high-pressure conditions, in particular, warm dense matter relevant to planetary science, astrophysics, and inertial confinement fusion. These targets provide optical access to the shocked sample, so optical velocity interferometry and pyrometry are typical diagnostics. We present a scaling for the coupling of the laser drive to the shock pressure in DACs based on over 150 prior experiments that can be used to design future experiments. For initial static pressures ≤ 5 GPa, using a sapphire window mitigates the effect of photoionization compared to diamond windows and allows the measurement of shock pressures up to 1 TPa with optical diagnostics.
Understanding the complex plasma dynamics in ultra-intense relativistic laser–solid interactions is of fundamental importance for applications of laser–plasma-based particle accelerators, the creation of high-energy-density matter, understanding planetary science, and laser-driven fusion energy. However, experimental efforts in this regime have been limited by the lack of accessibility of over-critical densities and the poor spatiotemporal resolution of conventional diagnostics. Over the last decade, the advent of femtosecond brilliant hard X-ray free-electron lasers (XFELs) has opened new horizons to overcome these limitations. Here, for the first time, we present full-scale spatiotemporal measurements of solid-density plasma dynamics, including preplasma generation with tens of nanometer scale length driven by the leading edge of a relativistic laser pulse, ultrafast heating and ionization at the main pulse arrival, the laser-driven blast wave, and transient surface return current-induced compression dynamics up to hundreds of picoseconds after interaction. These observations are enabled by utilizing a novel combination of advanced X-ray diagnostics including small-angle X-ray scattering, resonant X-ray emission spectroscopy, and propagation-based X-ray phase-contrast imaging simultaneously at the European XFEL-HED beamline station.
Understanding longer timescale material reactions under dynamic stress loading is critical for applications in materials engineering, shock physics, and planetary science. Prompt impulse experiments generate lower pressures since the ablator—the material first removed by the laser—is thicker and farther from the diagnostic plane, capturing aggregate material responses from the initial shock wave, rarefaction waves, and later time effects. This complexity demands thorough material characterization and simulation support. Since traditional sample construction is specific to supported shock experiments, designing prompt impulse experiments requires reconsideration around target design and sample engineering. Here, we present sample preparation techniques, experimental investigations, and theoretical simulations to investigate glue layer impacts, aiming to standardize samples for consistent data at lower laser fluences. We find that glue layers <30 μm have a minimal impact on peak velocity and pulse shape. The peak velocity scales linearly with glue layer thickness until a glue layer of 75 μm. For glue layers >75 μm, the peak velocity no longer scales with thickness; however, the pulse shape continues to degrade as described by simulations.
Significance Understanding the history of volatile species such as water in the Earth–Moon system is a major objective of planetary science. In this work, we use the moderately volatile element Rb, which has a long-lived isotope ( 87 Rb) that decays to 87 Sr, to show that lunar volatile element depletion was not caused by the Moon-forming impact. The Rb–Sr systematics of lunar rocks mandate that the bodies involved in the impact that formed the Earth–Moon system were depleted in volatile elements relative to the bulk solar system prior to the impact. As such, Earth’s relatively small proportion of water is either primarily indigenous or was added after the Giant Impact from a source that contained essentially no moderately volatile elements.
The radiogenic isotopic compositions of basaltic Martian meteorites (shergottites) and clinopyroxene/olivine cumulate meteorites (nakhlite/chassignites) are used to define the global evolution of Mars. However, the two main groups of meteorites demonstrate that their sources underwent divergent styles of magmatic evolution. The shergottites portray a planet that differentiated ~4.52 billion years ago via solidification of a magma ocean, producing incompatible element-depleted and -enriched reservoirs that remained isolated until melt production. In contrast, the reservoir from which the nakhlite/chassignites derive may have formed earlier, produced melts that fractionated Sm/Nd and Hf/W differently, was compositionally less variable, and experienced a significantly more complex history following primordial differentiation than the shergottite sources. The disparate histories recorded by these two groups of meteorites elucidate important questions that could be addressed by acquiring additional samples. Obtaining samples that shared the isotopic systematics of the shergottites would provide confidence that extrapolating the primordial differentiation history of Mars from shergottite radiogenic isotope systematics is reasonable. Returned samples from Mars will also constrain the physical locations of the meteorite source regions, providing insights into the general structure of the Martian mantle. In addition, they will help constrain the phases present in the martian mantle during melting and the conditions under which they are stable. Finally, identifying an evolved lithology that satisfies the geochemical and isotopic constraints placed on the incompatible element-enriched endmember observed in the shergottites would define the nature of magmatic evolution on Mars and whether it is more akin to processes on the Earth or the Moon.
SUMMARY Meteorite impacts have proved to be a significant source of seismic signal on the Moon, and have now been recorded on Mars by InSight seismometers. Understanding how impacts produce seismic signal is key to the interpretation of this unique data, and to improve their identification in continuous seismic records. Here, we use the seismic Representation Theorem, and particularly the stress glut theory, to model the seismic motion resulting from impact cratering. The source is described by equivalent forces, some resulting from the impactor momentum transfer, and others from the stress glut, which represents the mechanical effect of plasticity and non linear processes in the source region. We condense these equivalent forces into a point-source with a time-varying single force and nine-component moment tensor. This analytical representation bridges the gap between the complex dynamics of crater formation, and the linear point-source representation classically used in seismology. Using the multiphysics modelling software HOSS, we develop a method to compute the stress glut of an impact, and the associated point-source from hypervelocity impact simulations. For a vertical and an oblique impact at 1000 m s−1, we show that the moment tensor presents a significant deviatoric component. Hence, the source is not an ideal isotropic explosion contrary to previous assumptions, and draws closer to a double couple for the oblique impact. The contribution of the point force to the seismic signal appears negligible. We verify this model by comparing two signals: (1) HOSS is coupled to SPECFEM3D to propagate the near-source signal elastically to remote seismic stations; (2) the point-source model derived from the stress-glut theory is used to generate displacements at the same distance. The comparison shows that the point-source model is accurately simulating the low-frequency impact seismic waveform, and its seismic moment is in trend with Lunar and Martian impact data. High-frequencies discrepancies exist, which are partly related to finite-source effects, but might be further explained by the difference in mathematical framework between classical seismology and HOSS’ numerical modelling.
The thermodynamic properties of hexagonal-close-packed iron (ε–Fe) are essential for investigating the internal structure and dynamic properties of planetary cores. Despite their importance to planetary sciences, experimental investigations of ε–Fe at relevant conditions are still challenging. Therefore, ab initio calculations are crucial to elucidating the thermodynamic properties of this system. Here, we use a free energy calculation scheme based on the phonon gas model compatible with temperature-dependent phonon frequencies. We investigate the effects of electronic thermal excitations, which introduces a temperature dependence on phonon frequencies, and the implication for the thermodynamic properties of ε–Fe at extreme pressure (P) and temperature (T) conditions. We disregard phonon-phonon interactions, i.e., anharmonicity and their effect on phonon frequencies. Nevertheless, the current scheme is also applicable to T -dependent anharmonic frequencies. We conclude that the impact of thermal electronic excitations on vibrational properties is not significant up to ~4000 K at 200 GPa but should not be ignored at higher temperatures or pressures. However, the static free energy F st must always include the effect of thermal excitation fully in a continuum of T. Furthermore, our results for isentropic equations of state show good agreement with data from recent ramp compression experiments up to 1400 GPa conducted at the National Ignition Facility.
Experimental progress finally reached the metallic solid hydrogen phase, which was predicted by Wigner and Huntington over 80 years ago. However, the different structures in the phase diagram are still being debated due to the difficulty of diffraction experiments for high-pressured hydrogen. The determination of crystal structures under extreme condition is both of the basic condensed matter physics, and in planetary science: the behavior of giant gaseous planets (e.g., Jupiter and Saturn) strongly depends on the properties of inner high-pressured hydrogen. This work describes possible structures appearing under high pressures of 400~600 GPa. We applied a structural search using particle swarm optimization with density functional theory (DFT) to propose several candidate structures. For these structures, we performed fixed-node diffusion Monte Carlo simulations combined with DFT zero-point energy corrections to confirm their relative stability. We found P2 1 /c-8 as a promising candidate structure for the H 2 -PRE phase. P2 1 /c-8 is predicted the most stable at 400 and 500 GPa. P2 1 /c-8 reproduces qualitatively the IR spectrum peaks observed in the H 2 -PRE phase.
X-ray diffraction measurements under laser-driven dynamic compression now allow us to investigate the atomic structure of matter at TPa pressures and thousands of degree temperatures with broad implications for condensed matter physics, planetary science and astronomy. Furthermore, pressure determination in these experiments often relies on velocimetry measurements coupled with modeling that requires accurate knowledge of the optical and thermo-mechanical properties of a window material, resulting in significant systematic uncertainty. Here we demonstrate different approach applicable to X-ray diffraction experiments under quasi-isentropic ramp-compression based on the use of in-situ pressure calibrants, similar to the methods often adopted in static-compression experiments with diamond anvil cells. Focusing on experiments using a diamond window, we discuss challenges and mitigation strategies for the novel approach. Our study, in addition to providing new structural information of 5 metals up to hundreds of GPa, provides validation to the currently used methods based on time-resolved measurement of the diamond free-surface velocity, and reveals that the use of in-situ calibrants enables a factor of four reduction in the pressure uncertainty in these experiments.
The high-pressure, temperature phase diagram and transport properties of materials are of broad interest to planetary sciences and high-energy-density sciences and applications. Alumina (Al 2 O 3 ) or its various forms (e.g., solutions with other oxides or silicates) are important constituents in Earth’s and super-Earths’ mantles, common window materials in dynamic compression experiments, and standard pressure calibrators in diamond-anvil-cell experiments. Its structures and transport properties are of particular importance but have not been well studied at above 100 GPa pressures. Based on extensive first-principles molecular dynamics calculations, we obtain atomic level insights on structural differences and accurate results on electrical conductivity (σ dc ), thermal conductivity (κ), and reflectivity of Al 2 O 3 in several phases at pressures of ~130–1300 GPa and temperatures of 4000–20,000 K. We find the solid-to-liquid state changes of Al 2 O 3 are accompanied by an insulator-to-semimetal transition with enhanced σ dc and κ, which is similar to other silicates and oxides and can facilitate magnetic field generation in deep interiors of rocky planets. We also find the solid-to-solid transition (from Rh 2 O 3 (II) to CaIrO 3 ) is accompanied by negligible changes in σ dc and κ at 160 GPa, contrary to previous expectations, whereas another higher pressure phase transition (CaIrO 3 to U 2 S 3 at 450 GPa) is associated with increases in σ dc and κ. Furthermore, we show the transport properties can be enhanced by oxygen vacancy in Al 2 O 3 , but are not significantly affected by forming solutions with MgSiO 3 , particularly when in the solid states.
Identification of neutron sources is central to nuclear physics and its applications, from planetary science to nuclear security, yet direct source discrimination from measured neutron spectra remains fundamentally elusive. Here, we introduce a Bayesian protocol that directly infers source ensembles from measured neutron spectra by combining full-spectrum template matching with probabilistic evidence evaluation. Applying this protocol to recoil and time-of-flight spectroscopy, we recover single- and two-source configurations with strong statistical significance (beyond 4𝜎) at event counts as low as ∼10 3 . These results demonstrate that neutron spectral signatures can be leveraged for robust source identification, opening a new observational window for both fundamental research and operationally driven applications.
Existing COTS inorganic scintillators all have the characteristic of being very good at certain desirable properties, but not sufficient at other desirable properties for HEP. The demand for suitable scintillators (with regards to both scintillation detector properties and suitable pricing), to be used for modern intensities frontier (Mu2e-II), energy frontier (High luminosity large hadron collider) and future e+e- collider projects (aimed as Higgs bosons factory, such as the International Linear Collider (ILC) and the Future Circular Collider (FCC) are putting even higher challenges on crystal scintillators.In this work, we report the development of a novel high-performance scintillators that can address the issues associated with existing scintillators, the activator doped Hg2Br2. Initial results are very encouraging on the detection of gamma and alpha particles using a solid-state photomultiplier (SSPM). The responses have been stable and repeatable. Hg2Br2 also has many advantages over existing COTS scintillators such as: high density, bright, fast, good energy resolution, no intrinsic radiation, radiation hard and cost-effectiveness. Here, we present here why Hg2Br2 is the next generation scintillator for high energy physics experiments as well as other scientific and imaging applications such as planetary science and medical imaging.
In the event of a potentially catastrophic asteroid impact, with sufficient warning time, deploying a nuclear device remains a powerful option for planetary defense if a kinetic impactor proves insufficient. Predicting the effectiveness of a potential nuclear deflection or disruption mission depends on accurate multiphysics simulations of the device's x-ray energy deposition into the asteroid and the resulting material ablation. These simulations span many orders of magnitude, require a variety of different complex physics packages, and are computationally expensive. Having an efficient and accurate way of modeling this system is necessary for exploring a mission's sensitivity to the asteroid's range of physical properties. To expedite future simulations, we present a completed x-ray energy deposition model developed using the radiation-hydrodynamics code Kull which can be used to initiate a nuclear mitigation mission hydrocode. The model spans a wide variety of possible mission initial conditions: four different asteroid-like materials (Silicon Dioxide, Forsterite, Iron, and Ice), two different source spectra (1 and 2 keV blackbodies), and then a broad range of radiation fluences (0.0001 to 1 kt per square meter), source durations (10 to 100 ns), and asteroid porosities (0 to 80 percent). Using blowoff momentum as the primary metric, the modelinitiated simulation results match the full radiation-hydrodynamics results to within 10 percent. Please reference the journal article: Burkey et al., X-Ray Energy Deposition Model for Simulating Asteroid Response to a Nuclear Planetary Defense Mitigation Mission, Planetary Science Journal, (2023) for more information.
Apatite is a mineral of widespread importance in Earth and planetary science. Here we examine the behavior of a natural fluorapatite (FAp) crystal from Durango (Mexico) under compression to 61 GPa. Single-crystal X-ray diffraction experiments were carried out in a diamond-anvil cell using a synchrotron source. The apatite structure persists up to 32.4 GPa. Birch-Murnaghan equation of state parameters were fit to the pressure-volume data for fluorapatite for two cases: fixing V 0 at its measured ambient value resulted in a bulk modulus, K 0T , of 97.0(8) GPa and a pressure derivative of the bulk modulus, K' 0T , of 3.3(1), while fixing V 0 and K 0T at its ambient value 90.5 GPa (derived from ultrasonically measured elastic constants) resulted in a K' 0T value of 4.1(1). At 35.6 GPa, fluorapatite transforms to a triclinic phase (P$\bar{1}$, Z = 4), designated here as fluorapatite II (FAp-II). This phase persists up to at least 61 GPa. The major structural differences between FAp and FAp-II involve the buckling of the Ca polyhedra along the c-axis and changes in the number and coordination of the Ca sites. Furthermore, our study extends the pressure range over which fluorapatite has been examined by more than a factor of three, providing new insights into its structural response to high-pressure conditions.
Determining the internal structure of Uranus is a key objective for planetary science. Knowledge of Uranus’s bulk composition and the distribution of elements is crucial to understanding its origin and evolutionary path. In addition, Uranus represents a poorly understood class of “intermediate-mass planets" (intermediate in size between the relatively well studied terrestrial and gas giant planets), which appear to be very common in the Galaxy. As a result, a better characterization of Uranus will also help us to better understand exoplanets in this mass and size regime.
Abstract. Depressions – inwardly draining regions – are common to many landscapes. When there is sufficient moisture, depressions take the form of lakes and wetlands; otherwise, they may be dry. Hydrological flow models used in geomorphology, hydrology, planetary science, soil and water conservation, and other fields often eliminate depressions through filling or breaching; however, this can produce unrealistic results. Models that retain depressions, on the other hand, are often undesirably expensive to run. In previous work we began to address this by developing a depression hierarchy data structure to capture the full topographic complexity of depressions in a region. Here, we extend this work by presenting the Fill–Spill–Merge algorithm that utilizes our depression hierarchy data structure to rapidly process and distribute runoff. Runoff fills depressions, which then overflow and spill into their neighbors. If both a depression and its neighbor fill, they merge. We provide a detailed explanation of the algorithm and results from two sample study areas. In these case studies, the algorithm runs 90–2600 times faster (with a reduction in compute time of 2000–63 000 times) than the commonly used Jacobi iteration and produces a more accurate output. Complete, well-commented, open-source code with 97 % test coverage is available on GitHub and Zenodo.