Search NASA⌕ Search

SEARCH · Search NASA

Results for “SPACE MECHANICS”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A spatially-resolved model of neutron-irradiated tungsten coupling stochastic cluster dynamics and finite deformation plasticity

Structural materials used in nuclear reactors face severe degradation in mechanical properties, such as hardening and embrittlement. At the microscopic scale, this occurs due to creation and accumulation of irradiation-induced defects and their interaction with system dislocations. Although techniques exist which can model evolution of irradiation defects, for instance kinetic transport theory-based models, their interaction with mechanical deformation of the bulk material has not been investigated extensively. In this work, we demonstrate a novel spatially-resolved multiscale coupling between microscopic irradiation defect evolution, modeled using Stochastic Cluster Dynamics (SCD) and macroscopic mechanical deformation modeled using a finite-deformation plasticity model. SCD is used to determine the statistically averaged defect cluster spacing, dependent on operating conditions such as irradiation dose and temperature. This acts as an initial condition that governs the critical resolved shear stress of dislocation glide in the macroscopic plasticity model. This framework is used to predict mechanical behavior in post-mortem test of irradiated Tungsten samples, which has found its importance as structural material used in nuclear reactors. The results obtained using the coupled approach are in good agreement with experimental data of uniaxial tension tests. The model is able to capture the effect of temperature and irradiation dose on the material hardening. Two methods are proposed to estimate hardness – using Tabor's Law relating uniaxial yield stress to hardness and from flat-punch simulations. The results are in reasonable agreement with hardness data from micro-indentation experiments of irradiated Tungsten samples. Finally, the model is also able to reveal microstructural details such as spatial variation in defect density and local stress.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

A Unified Interpretation of Variability in Precipitation Isotope Ratios

Abstract Several mechanisms have been proposed to explain why the isotope ratios of precipitation vary in space and time and why they correlate with other climate variables like temperature and precipitation. Here, we argue that this behavior is best understood through the lens of radiative transfer, which treats the depletion of atmospheric vapor transport by precipitation as analogous to the attenuation of light by absorption or scattering. Building on earlier work by Siler et al., we introduce a simple model that uses the equations of radiative transfer to approximate the two-dimensional pattern of the oxygen isotope composition of precipitation ( δ p ) from monthly mean hydrologic variables. The model accurately simulates the spatial and seasonal variability in δ p within a state-of-the-art climate model and permits a simple decomposition of δ p variability into contributions from gradients in evaporation and the length scale of vapor transport. Outside the tropics, δ p is mostly controlled by gradients in evaporation, whose dependence on temperature explains the positive correlation between δ p and temperature (i.e., the temperature effect). At low latitudes, δ p is mostly controlled by gradients in the transport length scale, whose inverse relationship with precipitation explains the negative correlation between δ p and precipitation (i.e., the amount effect). This suggests that the temperature and amount effects are both mostly explained by the variability in upstream rainout, but they reflect distinct mechanisms governing rainout at different latitudes. Significance Statement The isotopic composition of precipitation has long been used to make inferences about past climates based on its observed relationship with precipitation in the tropics and with temperature at higher latitudes. These relationships—known as the “amount effect” and “temperature effect,” respectively—have been attributed to many different mechanisms, most of which are thought to operate at either high or low latitudes but not both. Here, we present a unified framework for interpreting the isotope variability that can explain the latitude dependence of the temperature and amount effects despite making no distinction between high and low latitudes. Although our results are generally consistent with certain interpretations of the amount effect, they suggest that the temperature effect is widely misunderstood.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Superluminal matter waves

The Dirac equation has resided among the greatest successes of modern physics since its emergence as the first quantum mechanical theory fully compatible with special relativity. This compatibility ensures that the expectation value of the velocity is less than the vacuum speed of light. Here, we show that the Dirac equation admits free-particle solutions where the peak amplitude of the wave function can travel at any velocity, including those exceeding the vacuum speed of light, despite having a subluminal velocity expectation value. The solutions are constructed by superposing basis functions with correlations in momentum space. These arbitrary velocity wave functions feature a near-constant profile and may impact quantum mechanical processes that are sensitive to the local value of the probability density as opposed to expectation values.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Inverse design of cellular structures with the targeted nonlinear mechanical response

Advanced additive manufacturing capabilities have enabled a transformational ability to create sophisticated cellular structures using diverse materials. By altering the topology of the unit cell, the mechanical behavior, such as the stress-strain response during compression, can be modulated. Nevertheless, identifying a printable topology within an enormous design space that would precisely deliver the targeted nonlinear material response is challenging. We propose a data-driven generative framework based on a conditional variational autoencoder (cVAE) architecture that can inverse design the cellular structure based on the intended nonlinear stress-strain response. Trained on a dataset of structure-property pairs, the cVAE learns a compact and expressive latent space that enables efficient mapping from targets to feasible geometries. Two inference modes are explored: (1) decoder-only generation, which enables the exploration of diverse designs conditioned solely on the desired mechanical response, and (2) encoder-decoder generation, which further allows for the incorporation of desired topologies, ensuring the generated structure conforms to both mechanical properties and to desired-topology constraints. The results demonstrate that the model can generate structurally plausible and mechanically accurate designs, with the predicted stress-strain curves closely matching the targets. Even under joint conditioning, the model effectively balances geometric fidelity and functional performance.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Ultrasonic-Assisted Extrusion Processing for Enhancing Physical Properties of High-Density Polyethylene by Flow-Induced Crystallization

The evolution of crystallinity resulting from stress imposed on a melt, known as flow-induced crystallinity, can strongly influence the mechanical and physical properties of semicrystalline polymers. This study investigates shear-induced crystallization by applying an ultrasonic field to the melt flow as it passes through dies with various geometries. A custom-built sonication die is employed for controlling the dynamic temperature and shear environment, resulting in molecular alignment and potential for flow-induced crystallization. Application of both conventional and ultrasonic shear rates at the equilibrium melt temperature of high-density polyethylene (HDPE) was investigated to accelerate crystallinity and manipulate the crystal morphology across the film in pursuit of improved mechanical and gas barrier properties without the need for additives or other polymer layers. The relationships among ultrasonic-assisted extrusion processing, polymer structure, and performance were analyzed using wide- and small-angle X-ray scattering (WAXS and SAXS), tensile testing, and oxygen transmission rate (OTR) analysis. Multiple linear regression models were implemented to predict the correlation among HDPE structure, process, and properties. Structural analysis revealed that both conventional and ultrasonic shear rates had the most significant influence on lamellar spacing and redistribution of rigid and soft amorphous fractions within the crystalline domains, ultimately dictating the mechanical and physical properties of the films. The goal is to explore the potential of the ultrasonic-assisted high crystallinity monolayer that can replace some of the functionality of complex, heterogeneous multilayer packaging with a single-material film having enhanced oxygen barrier properties.

crystallinity↗
DOE OSTI2026

Spaced out: An economic framework to explore the impacts of PV panel spacing on large-scale farming in Colorado

CONTEXT Agrivoltaic systems co-locate solar technologies with agricultural operations on an integrated plot of land and potentially provide benefits to both energy and agricultural systems. To date, large-scale (>5-MW) agrivoltaic projects in the United States have been limited to grazing and ecovoltaic applications, raising questions about the impact and scalability of agrivoltaic crop systems. Many agrivoltaic designs raise the height of the solar panels to accommodate agricultural practices while keeping energy density high. However, raising the panels results in increased photovoltaic (PV) development costs, which often are higher than the economic returns of crop production underneath the panels. This leads to unfavorable project economics and the need for other agrivoltaic solutions than raising panels. OBJECTIVE To explore other solutions, we perform an initial feasibility analysis for an agrivoltaic solution that can integrate with large-scale farming practices by increasing the row spacing in between panels. Increased PV row spacing is a low-cost approach for scaling agrivoltaics to accommodate crop production and this spacing can be tailored to required crop equipment for different regions. Increasing row spacing will reduce the power density (PV installed per acre), but in areas that are not land limited, these agrivoltaic designs could be economically feasible. Our analysis establishes a framework for a feasibility analysis for where and with what crops spaced out panel agrivoltaic solutions might be economical. METHODS Using a case study for large-scale agriculture crops in Colorado, we establish a framework for wide-row agrivoltaic economic feasibility analysis. We utilized the System Advisor Model to calculate technoeconomic metrics to compare different row spacing solutions and capture tradeoffs of these system designs. RESULTS AND CONCLUSIONS We find that, in some circumstances, wider row agrivoltaic solutions that allow for continued mechanized crop production can provide economic benefits over a traditional utility-scale PV system. For most crops examined in this analysis, roughly $\$$200/acre in agricultural profit justified spacing out the panels to at least 31.7 ft. to accommodate agrivoltaic configurations versus PV only configurations. Additionally, opportunities for increased agricultural revenue with agrivoltaic systems allow PV project economics to tolerate a larger range of CAPEX variability while remaining economically viable relative to the PV only configurations. SIGNIFICANCE This framework can be adapted for a wide variety of crops and regions and allows for examination of economically favorable sites for future agrivoltaic systems that utilize different configuration and expand opportunities for agrivoltaics.

14 SOLAR ENERGY↗
DOE OSTI2026

Theoretical assessment of the transition between electron emission mechanisms in a collisional crossed-field gap

Understanding the transition to the limiting current under vacuum and collisional conditions for thermal-field emission in a crossed-field gap is critical for characterizing the operating regime for numerous applications, such as high-power microwaves and Hall thrusters. By deriving the exact solution for the current density as a function of cathode temperature (and, by extension, initial velocity), cathode electric field, applied voltage, collision frequency, and magnetic field, we determine the transition from the canonical current source (thermal-field emission) to the vacuum and collisional limiting currents under different operating conditions. Here, we further show that the exact solution recovers the asymptotic solutions for the canonical and limiting currents for vacuum and collisional crossed-field gaps under appropriate limits. The exact solution demonstrates that the critical current decreases with both increasing magnetic field strength and increasing collisionality. Nexus phase space plots constructed by matching the individual asymptotic solutions provide the conditions for the transitions between mechanisms, elucidating when a single mechanistic equation would suffice and when theories combining multiple mechanisms are necessary.

Breen, Lorin I. [Purdue Univ., West Lafayette, IN ↗
DOE OSTI2024

Adaptive Scalpel Scanning Probe Microscopy for Enhanced Volumetric Sensing in Tomographic Analysis

Controlling nanoscale tip‐induced material removal is crucial for achieving atomic‐level precision in tomographic sensing with atomic force microscopy (AFM). While advances have enabled volumetric probing of conductive features with nanometer accuracy in solid‐state devices, materials, and photovoltaics, limitations in spatial resolution and volumetric sensitivity persist. This work identifies and addresses in‐plane and vertical tip‐sample junction leakage as sources of parasitic contrast in tomographic AFM, hindering real‐space 3D reconstructions. Novel strategies are proposed to overcome these limitations. First, the contrast mechanisms analyzing nanosized conductive features are explored when confining current collection purely to in‐plane transport, thus allowing reconstruction with a reduction in the overestimation of the lateral dimensions. Furthermore, an adaptive tip‐sample biasing scheme is demonstrated for the mitigation of a class of artefacts induced by the high electric field inside the thin oxide when volumetrically reduced. This significantly enhances vertical sensitivity by approaching the intrinsic limits set by quantum tunneling processes, allowing detailed depth analysis in thin dielectrics. The effectiveness of these methods is showcased in tomographic reconstructions of conductive filaments in valence change memory, highlighting the potential for application in nanoelectronics devices and bulk materials and unlocking new limits for tomographic AFM.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Structural Analysis of a Serpentine Superconducting Magnet for the Interaction Region of the Electron Ion Collider

The electron ion collider (EIC) is under design to collide the high energy and highly polarized hadron with electron beams with luminosities up to 10 34 cm −2 s −1 . The superconducting magnet designs at the interaction region (IR) are quite challenging due to the close proximity of the hadron and electron beams. Several serpentine types of superconducting magnets have been designed due to the space restrictions at the IR. Here, to validate the designs in terms of the mechanical strength, structural integrity, as well as the allowed deformations, a detailed structural mechanical analysis for a 12-layer superconducting quadrupole serpentine magnet (Q1BpR) in the rear side hadron beams has been carried out using the finite element method (FEM) considering the pretension stress, shrinking due to cool down, and the electromagnetic force on the conductors. Anisotropy material properties and individual roving tension have been considered. The electro-magnetic (EM) simulations were performed using COMSOL, confirmed with RAT, the structural mechanical analysis were performed using COMSOL by considering the contact elements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Optimizing Metadata Exchange: Leveraging DAOS for ADIOS Metadata I/O

In HPC I/O middleware like the Adaptable I/O System (ADIOS) often mediates data transfers between applications. The metadata I/O generated by such systems often presents significant scaling and performance limitations. This work seeks improvement opportunities for metadata I/O by leveraging the DAOS storage systems, a recent storage system solution deployed on high-end systems such as the Aurora supercomputer. We investigate the tradeoffs and the design space for integrating I/O engines for the ADIOS middleware based on the different storage mechanisms supported by DAOS. We present a new DAOS-Array-ChunkSize-aligned engine which provides up to 2.3× improved performance than when using the existing DAOS-POSIX interface, without requiring any application modifications.

Venkatesh, Ranjan Sarpangala↗
DOE OSTI2026

Strong-Drive Limits in Josephson Circuits: From Chaos to an Unbound-Resonance Threshold

Strong microwave drives enable fast measurement and parametric control in superconducting circuits but can induce transitions out of the intended low-energy manifold. We develop a unified description of strong-drive limits in flux- and charge-driven Josephson circuits across drive frequency and dc flux bias. Using classical phase-space analysis and Floquet--Markov simulations, we identify distinct low- and high-frequency mechanisms. At low frequency, we characterize bound-state resonances and separatrix chaos and find that the flux-drive chaos threshold depends strongly on dc flux bias. At high frequency, these mechanisms are suppressed, and the dissipative steady state transfers from the central bound-state sector to outer resonances formed by above-barrier running trajectories. The resulting unbound-resonance threshold is nearly independent of drive frequency and circuit parameters over the regime studied and is controlled primarily by dc flux bias. Coherent simulations show that parametric operation persists beyond this threshold, but at a reduced rate, setting an effective upper bound on the achievable operation speed. We derive analytical criteria for both thresholds, validate them numerically, and experimentally confirm the predicted dc-bias dependence of the low-frequency threshold in a flux-driven SQUID. We also determine the timescales for transfer into the unbound-resonance regime and relaxation back to the bound-state manifold after the drive is removed. Finally, we relate the stability limits to a complementary picture based on the junction critical current and extend the framework to multitone drives and inductively shunted circuits. Together, these results identify the mechanisms limiting strong driving and suggest routes to extend the stable operating range of Josephson circuits.

You, Xinyuan [Fermilab] (ORCID:0000000291789419)↗
DOE OSTI2024

Design And Control of Thermal Storage for Ventilative Cooling in Multifamily Buildings (Final Report)

The concept of natural ventilation is to provide a heat sink for cooling the building during the occupancy periods at night and improve indoor thermal comfort during the daytime without or with minimal need for mechanical cooling. Commonly, natural ventilation through the window or mechanical circulation of the air through the ventilation ducts are two key methods used to achieve night ventilation. Due to the increasing electric need for space heating and cooling, renewable energies combined with energy storage systems are receiving attention worldwide nowadays. When designed and controlled properly, this strategy can reduce energy use for cooling, reduce the size of mechanical cooling equipment, reduce peak electrical demand for cooling, and better enable demand flexibility for mechanical cooling equipment. In some cases, ventilative cooling can eliminate the need for mechanical cooling altogether. For these reasons, California’s Building Energy Efficiency Standards have recently added prescriptive requirements that all new single-family residences (in most California climates) must include ventilative cooling systems (aka: “whole house fans”). However, the current standards do not address multifamily buildings because market-available ventilative cooling products are designed for single-family residences and are physically incompatible with many multifamily building archetypes. The results from the study show that the combination of thermal energy storage on building walls along with nighttime ventilation cooling assisted with nighttime ventilation can save total building energy usage by up to 9% and peak load by 16% in moderate climate zones. The controlled nighttime ventilation cooling (activate ventilation only when outdoor conditions are favorable) can also provide total electricity savings of up to 5% and 7% in hot and cold climates, respectively. In the case of moderate climates, kitchen exhaust-assisted nighttime ventilation cooling can eliminate the need for cooling during overnight hours. This can also significantly assist in reducing carbon emissions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Robust finite-temperature many-body scarring on a quantum computer

Mechanisms for suppressing thermalization in disorder-free many-body systems, such as Hilbert space fragmentation and quantum many-body scars, have recently attracted much interest in foundations of quantum statistical physics and potential quantum information processing applications. However, their sensitivity to realistic effects such as finite temperature remains largely unexplored. Here, we have utilized IBM's Kolkata quantum processor to demonstrate an unexpected robustness of quantum many-body scars at finite temperatures when the system is prepared in a thermal Gibbs ensemble. We identify such robustness in the PXP model, which describes quantum many-body scars in experimental systems of Rydberg atom arrays and ultracold atoms in tilted Bose-Hubbard optical lattices. By contrast, other theoretical models which host exact quantum many-body scars are found to lack such robustness and their scarring properties quickly decay with temperature. Our study sheds light on the important differences between scarred models in terms of their algebraic structures, which impacts their resilience to finite temperature. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Pathways for Nucleation and Growth in Confined Spaces and at Interfaces

Mineral crystallization is central to myriad natural processes from the formation of snowflakes to stalagmites, but the molecular-scale mechanisms are often far more complex than models reflect. Feedbacks between the hydro-, bio-, and geo-spheres drive complex crystallization processes that challenge our ability to observe and quantify them, motivating an expansion of crystallization theories. Here, in this article, we discuss how the driving forces and timescales of nucleation are influenced by factors ranging from simple geometric confinement to distinct interfacial solution structures involving solvent organization, electrical double layers, and surface charging effects. Taken together, these ubiquitous natural phenomena can preserve metastable intermediates, drive precipitation of undersaturated phases, and modulate crystallization in time and space.

biosphere↗
DOE OSTI2025

Unveiling Atomistic Mechanisms Governing Additive Manufacturing Processability and Mechanical Behavior of a Refractory Complex Concentrated Alloy

Extending the concept of complex concentrated alloys (CCAs) to the refractory alloys (solidus temperature over 2000 °C) space potentially facilitates the design of lightweight structural alloys with service temperatures that exceed those of Ni and Co‐based alloys. However, the room and elevated temperature tensile properties of the current refractory‐CCAs (R‐CCAs) are inferior to those of the Ni/Co‐based alloys. Furthermore, the manufacturing scalability of R‐CCAs remains challenging, in that cracks are prevalent in all R‐CCAs when processed using near‐net shape manufacturing processes, such as fusion‐based additive manufacturing (F‐BAM). Still, mechanisms governing the poor F‐BAM processability of R‐CCAs remain unexplored. Here, to this end, this work unveils the atomistic mechanisms underlying F‐BAM process‐induced cracking in a NbTiTaMoHfZrC R‐CCA. The implications of light elements’ presence for intrinsic ductility and grain boundary cohesion, and subsequently for F‐BAM processability and mechanical behavior, are revealed. Leveraging the insights, we accomplish what is, to the best of the knowledge, the first instance of crack‐free F‐BAM processing of any R‐CCA. Additionally, the R‐CCA exhibits over 20% tensile ductility and ≈160 MPa tensile yield strength at 1200 °C. In addition to facilitating the design of lightweight R‐CCAs, findings enable scalable manufacturing of these ultra‐high temperature alloys for structural applications.

Refractory alloys↗
DOE OSTI2025

Microstructure‐Dependent Sodium Storage Mechanisms in Hard Carbon Anodes

Sustainable energy storage is essential to support the transition to renewables and meet the increasing demand for energy. Sodium‐ion batteries (NIBs) are attractive for grid‐scale energy storage due to the abundance and low cost of sodium, sustainability of other battery components, and electrochemical performance. Hard carbon (HC) is a leading anode material for NIBs, but its complex microstructure complicates the understanding of sodium storage mechanisms. Using X‐ray total scattering and density functional theory calculations, this study clarifies how HC's microstructural variations influence sodium storage across the slope (high potential) and plateau (low potential) regions of the potential capacity curve. In the slope region, sodium initially adsorbs at high‐binding energy defect sites and subsequently intercalates between graphene layers, adsorbing at low‐binding energy defect sites, correlating with different slopes observed during initial sodiation. Initial irreversibility arises from sodium trapping at surface defects and solid electrolyte interface formation. In the plateau region, sodium simultaneously intercalates and fills pores, influenced by pore size, interlayer spacing, and defect concentration. HCs with larger pore sizes form larger sodium clusters. In conclusion, the proposed mechanism underscores the role of microstructure engineering in enhancing HC performance and advancing NIBs for grid‐scale energy storage.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Delving into the depths of NGC 3783 with XRISM

We present our study of the X-Ray Imaging and Spectroscopy Mission (XRISM) observation of the Seyfert-1 galaxy NGC 3783. XRISM’s Resolve microcalorimeter has enabled, for the first time, a detailed characterization of the highly ionized outflows in this active galactic nucleus. Our analysis constrains their outflow and turbulent velocities, along with their ionization parameter (ξ) and column density (N H ). The high-resolution Resolve spectrum reveals a distinct series of Fe absorption lines between 6.4 and 7.8 keV, ranging from Fe XVIII to Fe XXVI. At lower energies (1.8−3.3 keV), absorption features from Si, S, and Ar are also detected. Our spectroscopy and photoionization modeling of the time-averaged Resolve spectrum uncovers six outflow components, five of which exhibit relatively narrow absorption lines with outflow velocities ranging from 560 to 1170 km s −1 . In addition, a broad absorption feature is detected, which is consistent with Fe XXVI outflowing at 14 300 km s −1 (0.05 c). The kinetic luminosity of this component is 0.8−3% of the bolometric luminosity. Our analysis of the Resolve spectrum shows that more highly ionized absorption lines are intrinsically broader than those of lower-ionization species, indicating that the turbulent velocity of the six outflow components (ranging from 0 to 3500 km s −1 ) increases with ξ. Furthermore, we find that the column density (N H ) of the outflows generally declines with the ionization parameter up to log ξ = 3.2 but rises beyond this point, suggesting a complex ionization structure. The absorption profile of the Fe XXV resonance line is intriguingly similar to UV absorption lines (Lyα and C IV) observed by the Hubble Space Telescope, from which we infer that the outflows are clumpy in nature. Our XRISM/Resolve results from lower- and higher-ionization regimes support a “hybrid wind” scenario in which the observed outflows have multiple origins and driving mechanisms. We explore various interpretations of our findings within active galactic nucleus wind models.

X-rays: galaxies↗
DOE OSTI2024

Air-coupled tsunamis generated from impacts and airbursts: Our understanding before Hunga-Tonga Hunga-Ha'apai

The effort to prevent or mitigate the effects of an impact on Earth is known as planetary defense. A significant component of planetary defense research involves risk assessment. Much of our understanding of the risk from near-Earth objects comes from the geologic record in the form of impact craters, but not all asteroid impacts are crater-forming events. Small asteroids explode before reaching the surface, generating an airburst, and most impacts into the ocean do not penetrate the water to form a crater in the sea floor. The risk from these non-crater-forming ocean impacts and airbursts is difficult to quantify and represents a significant uncertainty in our assessment of the overall threat. We are currently working to better understand impact scenarios that can generate dangerous tsunamis. One of the suggested mechanisms for the production of asteroid–generated tsunamis is by direct coupling of the pressure wave to the water, analogous to the means by which a moving weather front can generate a meteotsunami. To test this hypothesis, we ran a series of airburst simulations and provided time-resolved pressure and wind profiles to use as source functions for tsunami models. We used the CTH hydrocode to model the various airburst scenarios to compare to the results of other simulations and provide time dependent boundary conditions as input to shallow-water wave propagation codes. The strongest and most destructive meteotsunamis are generated by atmospheric pressure oscillations with amplitudes of only a few hPa1 (mbar), corresponding to changes in sea level of a few cm. The resulting wave is strongest when there is a resonance between the ocean and the atmospheric forcing. A Proudman resonance takes place when the atmospheric disturbance’s translational speed (U) equals the longwave phase speed $\sqrt{gh}$ of shallow water wave. Coupling is strongest when the Froude number (Fr=U/c) is unity. A weather front propagates much slower than the speed of sound, so meteotsunamis are most common and dangerous in shallow bodies of water such as the Mediterranean Sea or Lake Michigan. By contrast, the blast wave from an airburst or crater-forming impact propagates at a speed faster than a tsunami in the deepest ocean, and a Proudman resonance cannot be achieved even though the overpressures are orders of magnitude greater. However, blast wave profiles are N-waves in which a sharp shock wave leading to overpressure is followed by a more gradual rarefaction to a much longer-duration underpressure phase. Even though the blast outruns the water wave it is forcing, the tsunami should continue to be driven by the out-of-resonance gradient associated with the suction phase, which may depend strongly on the details of the airburst or impact scenario. The open question is whether there are any conditions under which such an airburst-driven tsunami can be dangerous enough to contribute to the overall impact risk. We have also identified other potential mechanisms for airburst-generated tsunamis: 1) reaction force at the surface from the plume ejected into space, which carries significant momentum, 2) expanding toroidal vortices at the surface, which travel more slowly than the shock wave and can generate a Proudman resonance in relatively shallow ocean (such as continental shelf), and 3) steam explosion from seawater ablation by a “Type II” (Libyan Desert Glass-type) airburst in which the hot vapor jet descends to the surface. On January 15, 2022, the Hunga-Tonga Hunga-Ha’apai volcano, located approximately 60 km north of Tongatapu, the main island of Tonga, violently erupted with a powerful explosion, culminating the period of volcanic activity that started in December of 2021. This event and resulting tsunamis provided an existence proof for the air pressure wave coupling mechanism we proposed. It also suggests that it can be stronger and more significant over much greater distances than we contemplated, leading to global tsunamis associated with impact events on land as well as in the water. Large atmospheric explosions generate global Lamb waves with larger amplitudes, longer periods, and slower speeds than the local and regional blast waves we modeled prior to that event. This paper reviews our analysis and modeling of airburst-driven tsunamis prior to the 2022 Hunga-Tonga Hunga-Ha’apai tsunami, which was the subject of two presentations at the 2023 Planetary Defense Conference and is the subject of another paper currently in preparation.

54 ENVIRONMENTAL SCIENCES↗