Search NASA⌕ Search

SEARCH · Search NASA

Results for “Critical Velocity”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Impact of Cloud-Base Turbulence on CCN Activation: Single-Size CCN

Abstract This paper examines the impact of cloud-base turbulence on activation of cloud condensation nuclei (CCN). Following our previous studies, we contrast activation within a nonturbulent adiabatic parcel and an adiabatic parcel filled with turbulence. The latter is simulated by applying a forced implicit large-eddy simulation within a triply periodic computational domain of 64 3 m 3 . We consider two monodisperse CCN. Small CCN have a dry radius of 0.01 μ m and a corresponding activation (critical) radius and critical supersaturation of 0.6 μ m and 1.3%, respectively. Large CCN have a dry radius of 0.2 μ m and feature activation radius of 5.4 μ m and critical supersaturation 0.15%. CCN are assumed in 200-cm −3 concentration in all cases. Mean cloud-base updraft velocities of 0.33, 1, and 3 m s −1 are considered. In the nonturbulent parcel, all CCN are activated and lead to a monodisperse droplet size distribution above the cloud base, with practically the same droplet size in all simulations. In contrast, turbulence can lead to activation of only a fraction of all CCN with a nonzero spectral width above the cloud base, of the order of 1 μ m, especially in the case of small CCN and weak mean cloud-base ascent. We compare our results to studies of the turbulent single-size CCN activation in the Pi chamber. Sensitivity simulations that apply a smaller turbulence intensity, smaller computational domain, and modified initial conditions document the impact of specific modeling assumptions. The simulations call for a more realistic high-resolution modeling of turbulent cloud-base activation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Sonic horizons and causality in phase transition dynamics

A system gradually driven through a symmetry-breaking phase transition is subject to the Kibble-Zurek mechanism (KZM). As a consequence of the critical slowing down, its state cannot follow local equilibrium, and its evolution becomes nonadiabatic near the critical point. In the simplest approximation, that stage can be regarded as an “impulse” where the state of the system remains unchanged. It leads to the correct KZM scaling laws. However, such a “freeze-out” might suggest that the coherence length of the nascent order parameter remains unchanged as the critical region is traversed. By contrast, the original causality-based discussion emphasized the role of the sonic horizon: Domains of the broken symmetry phase can expand with a velocity limited by the speed of the relevant sound. This effect was demonstrated in the quantum Ising chain where the dynamical exponent z = 1 and quasiparticles excited by the transition have a fixed speed of sound. To elucidate the role of the sonic horizon, in this paper we study two systems with z > 1 where the speed of sound is no longer fixed, and the fastest excited quasiparticles set the size of the sonic horizon. Their effective speed decays with increasing transition time. In the extreme case, the dynamical exponent z can diverge such as in the Griffiths region of the random Ising chain where the localization of excited quasiparticles freezes the growth of the correlation range when the critical region is traversed. Of particular interest is an example with z < 1 , the long-range extended Ising chain, where there is no upper limit to the velocity of excited quasiparticles with small momenta. Initially, the power-law tail of the correlation function grows adiabatically, but in the nonadiabatic stage it lags behind the adiabatic evolution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

One-Dimensional Multi-Velocity Capabilities for Arbitrary Lagrangian-Eulerian Normal Contact Mechanics

Lagrangian and Arbitrary Lagrangian Eulerian (ALE) hydrodynamics codes such as FLAG form the backbone of many mission-critical multi physics simulations at Los Alamos National Laboratory. Critical to pre forming high fidelity simulations with these codes are Lagrangian and ALE contact algorithms, which allow materials to collide, slide, and sep arate throughout a simulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A Semi-Empirical Correlation for the Onset of Density Wave Oscillations in a Helical Coil Steam Generator

NuScale Power has completed full length helical coil steam generator (HCSG) tests at SIET Laboratories in Piacenza, Italy. The test program included an investigation of the onset and evolution of density wave oscillations (DWO). This paper describes the mechanisms leading to the onset of DWO as inferred from tests performed in the full-height TF-2 counter-flow test facility. This research supports the concept that the onset of DWO is caused by liquid bridging inside the boiling length of the tubes. Using a separate cylinders approximation, a DWO onset correlation is obtained by performing a steady-state momentum balance on the density front created by liquid bridging. This yields a momentum balance equation in terms of a set of dimensionless superficial velocities, a corresponding dimensionless pressure drop term and an average void fraction. Envelope theory is then used to obtain a critical void fraction which is substituted into the dimensionless momentum balance to obtain a DWO onset correlation similar to the Wallis counter-current flooding correlation. A comparison of the DWO onset correlation to the TF-2 DWO experimental data is presented and shows excellent agreement when the average equilibrium vapor quality is greater than 0.2.

Density Wave Oscillations, Helical Coil Steam Gene↗
DOE OSTI2020

A Satellite-Based Estimate of Convective Vertical Velocity and Convective Mass Flux: Global Survey and Comparison with Radar Wind Profiler Observations

Convective vertical velocity (w c ) and convective mass flux (M c ) lie at the heart of GCM cumulus parameterizations, but few observations of these critical parameters are available. In this paper, we develop and evaluate a novel, satellite-based method for estimating profiles of w c and M c . Here, comparisons with collocated ground-based radar wind profiler (RWP) observations show that satellite estimated median w c is slightly greater than the RWP estimates, but they show solid agreement when compared at the 95th percentiles (intense updrafts). RWP-derived and satellite estimated M c are broadly comparable in the lower and middle troposphere, with some differences in the upper troposphere due to differences in convective core sampling. A k-means cluster analysis of multiple years of w c data shows that convective characteristics are distinctly different among extratropical convection, tropical land convection, and tropical oceanic convection. Tropical land convection is significantly more intense and more variable than the oceanic counterpart.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Mapping the magnetic field in the Taurus/B211 filamentary cloud with SOFIA HAWC + and comparing with simulation

ABSTRACT Optical and infrared polarization mapping and recent Planck observations of the filametary cloud L1495 in Taurus show that the large-scale magnetic field is approximately perpendicular to the long axis of the cloud. We use the HAWC + polarimeter on SOFIA to probe the complex magnetic field in the B211 part of the cloud. Our results reveal a dispersion of polarization angles of 36°, about five times that measured on a larger scale by Planck. Applying the Davis–Chandrasekhar–Fermi (DCF) method with velocity information obtained from Institut de Radioastronomie Millimétrique 30 m C18O(1-0) observations, we find two distinct sub-regions with magnetic field strengths differing by more than a factor 3. The quieter sub-region is magnetically critical and sub-Alfv$\acute{\rm e}$nic; the field is comparable to the average field measured in molecular clumps based on Zeeman observations. The more chaotic, super-Alfv$\acute{\rm e}$nic sub-region shows at least three velocity components, indicating interaction among multiple substructures. Its field is much less than the average Zeeman field in molecular clumps, suggesting that the DCF value of the field there may be an underestimate. Numerical simulation of filamentary cloud formation shows that filamentary substructures can strongly perturb the magnetic field. DCF and true field values in the simulation are compared. Pre-stellar cores are observed in B211 and are seen in our simulation. The appendices give a derivation of the standard DCF method that allows for a dispersion in polarization angles that is not small, present an alternate derivation of the structure function version of the DCF method, and treat fragmentation of filaments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Measurement of Thomson-scattering spectra with continuous angular resolution (invited)

Here, a novel Thomson-scattering diagnostic with continuous angular resolution over a span of 120° was developed for the characterization of plasmas produced at the Omega Laser Facility. Spectrally resolving light scattered from electron plasma wave features as a function of emission angle provides a means to efficiently probe a large range of plasma frequencies and k vectors. Together, these spectra contain critical constraints on the plasma-physics models used to interpret the data and enable experimental measurements of the electron-velocity distribution function over several orders of magnitude without assumptions about its mathematical form. Major components of the instrument include (1) a reflective collection objective that gathers light over a range of 120° × 12°; (2) a spatial-filter image relay for measurement localization; (3) cylindrical optics for producing a line image of the collection aperture; (4) a transmission grating spectrometer; and (5) a time-gated, image-intensified camera. Thomson-scattered light collected from an ~(50 – μM) 3 volume of plasma is recorded with 0.8-nm spectral and 1° angular resolution. Initial experiments examined the properties of the electron-velocity distribution in gas-jet-produced plasmas in the presence of heating via inverse bremsstrahlung absorption.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Role of Inlet Boundary Conditions on Fuel-Air Mixing at Supercritical Conditions

Advanced engine design and alternative fuels present the possibility of fuel injection at purely supercritical conditions in diesel engines and gas turbines. The complex interactions that govern this phenomenon still need significant research, particularly the boundary conditions for fuel injection are critical for accurate simulation. However, the flow inside the injector itself is often omitted to reduce the computational efforts, and thus, velocity, mass flux, or total pressure is specified at the injector exit (or domain inlet), often with simplified velocity profiles and turbulence levels. This simplified inlet boundary treatment has minimal effects on results for conventional fuel injection conditions, however, the validity of this approach at supercritical conditions has not been assessed. Comprehensive real-gas and binary fluid mixing models have been implemented for computational fluid dynamics (CFD) analysis of fuel-air mixing at supercritical conditions. The model is verified using prior CFD results from the literature. The model is used to investigate the effects of the shape of axial velocity and mass fraction profiles at the inlet boundary with the goal to improve the comparison of predictions to experimental data. Results show that the boundary conditions have a significant effect on the predictions, and none of the cases match precisely with experimental data. Furthermore, the study reveals that the physical location of the inlet boundary might be difficult to infer correctly from the experiments and highlights the need for high-quality, repeatable measurements at supercritical conditions to support the development of relevant high-fidelity models for fuel-air mixing.

02 PETROLEUM↗
DOE OSTI2025

Investigation into the instantaneous centre of rotation for enhanced design of floating offshore wind turbines

The dynamic behaviour of floating offshore wind turbines (FOWTs) involves complex interactions of multivariate loads from wind, waves, and currents, which result in complex motion characteristics. Although methods for analysing global motion responses are well-established, the time- and location-dependent kinematics remain underexplored. This paper investigates the instantaneous centre of rotation (ICR), a point of zero velocity at a time instance of general plane motion. Understanding and strategically positioning the ICR can reduce the dynamic motion in critical structural locations, enhancing the performance and structural robustness of FOWTs. The paper presents a method for computing the ICR using time-domain simulation results and proposes a statistical analysis approach suitable for design studies. Building on prior research, it examines the sensitivity of the ICR to external loading and design features, providing insights into how these factors influence motion response and how the motion response influences the statistics of the ICR, structural loads, and other performance metrics of interest. The study explores two FOWT configurations, a spar and a semisubmersible, identifying design variables that most effectively control the ICR statistics and identifying the ICR statistics most correlated with the responses of interest. Finally, through two case studies, we demonstrate how to apply these new insights in a practical design scenario. By adjusting the design variables most correlated with the ICR (fairlead vertical position and centre of mass for the spar and mooring line length and offset column diameter for the semisubmersible), we successfully modified the designs of the floating support structures to reduce the loads in the mooring lines, tower base, and blade roots, improving the ultimate strength and fatigue characteristics compared to the original designs.

17 WIND ENERGY↗
DOE OSTI2025

Stokes-dependent droplet collection efficiency on a NACA 0012 airfoil from droplet-informed simulations with statistical overloading

Accurate modelling of ice accretion on aircraft wings requires analysing droplet impingement on the surface to optimize the design of ice-protection systems. We perform Euler–Lagrange simulations of a droplet-laden flow impinging on a NACA 0012 airfoil. Our study includes water droplets with eight discrete sizes ranging from 1 to 160 microns. We vary the free-stream velocity of the incoming airflow in the range 60 ≤ U ≤ 240 m s −1 and the chord length of the airfoil in the range 0.5 ≤ c ≤ 2 m. Due to the dilute nature of supercooled clouds, one-way coupling is used in the simulations. The effects of droplet breakup and collision are also neglected. To reduce the computational cost, we employ statistical overloading of droplets, allowing us to simulate millions of impinging droplets in a time span on the order of milliseconds. Our results show that the droplet collection efficiency, which measures the likelihood of droplet impingement on the airfoil surface, increases with droplet size and free-stream velocity but decreases with airfoil size. We demonstrate that collection efficiency, impingement velocity and impingement angle are primarily dictated by a single non-dimensional parameter, the droplet Stokes number. We also identify a critical stagnation-streamline Stokes number below which impingements do not occur and use it to estimate the minimum droplet size for impingement. In addition, we observe droplet behaviour to become Stokes number independent at large values of the Stokes number. This article is part of the theme issue ‘Heat and mass transfer in frost and ice’.

Science & Technology - Other Topics↗
DOE OSTI2023

Scrambling Transition in a Radiative Random Unitary Circuit

Here we study quantum information scrambling in a random unitary circuit that exchanges qubits with an environment at a rate p. As a result, initially localized quantum information not only spreads within the system, but also spills into the environment. Using the out-of-time-order correlator (OTOC) to characterize scrambling, we find a nonequilibrium phase transition in the directed percolation universality class at a critical swap rate p c : for p < p c the ensemble-averaged OTOC exhibits ballistic growth with a tunable light cone velocity, while for p > p c the OTOC fails to percolate within the system and vanishes uniformly within a finite timescale, indicating that all local operators are rapidly swapped into the environment. To elucidate its information-theoretic consequences, we demonstrate that the transition in operator spreading coincides with a transition in an observer’s ability to decode the system’s initial quantum information from the swapped-out, or “radiated,” qubits. We present a simple decoding scheme which recovers the system’s initial information with perfect fidelity in the nonpercolating phase, and with continuously decreasing fidelity with decreasing swap rate in the percolating phase. Depending on the initial state of the swapped-in qubits, we further observe a corresponding entanglement transition in the coherent information from the system into the radiated qubits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Dominant Controls on Preferential Flow and Their Implications for Future Soil Water Fluxes

Abstract Soil water flow, particularly preferential flow (PF), is a critical control on hydrological and biogeochemical processes, including groundwater recharge, contaminant transport, and carbon cycling. However, it remains challenging to predict PF occurrence across large environmental gradients. Here, we developed a deep learning (DL) model to estimate event‐scale soil water flow velocity and the probability of PF occurrence using high‐frequency soil moisture and precipitation data from 33 sites across the National Ecological Observatory Network. The model demonstrated high skill in predicting the binary occurrence of PF (91% F1‐score; 85% accuracy) but the performance was limited in predicting soil water velocity ( R 2 = 0.31). We found that precipitation characteristics (duration, volume, and intensity) were the most important predictors for soil water velocity. Among the non‐precipitation event variables, sand content showed relatively high predictive skill, though differences among non‐event climate variables were generally modest. Lower sand content was associated with increased predicted soil water velocity, a finding that highlights the role of soil structure in producing more non‐uniform flow, which contrasts with traditional uniform flow models. Projecting a reduced DL model under both moderate and high‐emissions future climate scenarios (2060–2099 Representative Concentration Pathways 4.5 and 8.5), we found ∼7.3% increase under RCP4.5 and ∼15% under RCP8.5 of soil water velocities compared to the historical simulation, while modeled likelihood of PF changed little. These findings suggest climate change is not making PF more frequent, but it is making existing PF pathways more efficient with important consequences for associated nutrient and contaminant transport under climate change. Plain Language Summary Water movement in soil is critical for water quality. While often modeled as a uniform flow process, in reality water moves rapidly through cracks and burrows in what is called “preferential flow” (PF), which limits natural filtration and can transport pollutants. We developed a deep learning model, trained on data from 33 U.S. sites, to predict when and how fast this PF occurs based on precipitation, soil, and climate data. The model showed that precipitation characteristics (duration, intensity, volume) were the most important predictors of PF. Lower soil sand content/higher clay content was associated with faster water flow, likely due to clay soils forming aggregates and cracks that water moves through rather than infiltrating uniformly. Further analyses based on climate projections suggest that the speed at which PF occurs will become more rapid under future climate scenarios compared to historical simulation. This highlights the need to represent PF in soil water models when assessing future water quality. Key Points The effect of precipitation peak intensity on soil water velocities declined with increasing precipitation intensity Antecedent soil moisture failed to predict preferential flow (PF), contrasting the high predictive power of sand content Climate predictions suggest that soil water velocities through PF paths will increase ∼15% by 2099

Li, Bonan↗
DOE OSTI2023

Low-Velocity Impact Performances of Healed Polymer Fiber Reinforced Plastics

Extending the lifecycle of traditional carbon or glass fiber-reinforced plastics is a complicated problem. The lack of sustainability limits the applications of the traditional composite materials in the vehicle industries where recycling and repurposing are critical issues. Alternatives for the low-stressed structural components are polymer fiber-reinforced plastics (PFRPs). In PFRPs, both the fibers and matrix are composed of thermoplastic polymers (e.g. polypropylene or polyethylene). They are lightweight, easy to manufacture, and cost-effective. Additionally, recycling and repurposing thermoplastic polymers are well understood. Therefore, the PFRPs have strong advantages compared to the traditional fiber-reinforced composites in low-stressed structural applications. In this study, we investigated the low-velocity impact (LVI) performances of the PFRPs and compared them with carbon fiber-reinforced plastics (CFRPs). A semi-spherical impactor was dropped to flat, square panels, and the absorbed impact energy was measured. The damage mechanisms were examined using a Xray µCT scan. The PFRPs outperformed the CFRPs in terms of perforation energy normalized by plate thickness and density. After the perforation, we healed the fractured plates by leveraging the recyclability of the thermoplastic polymers. The healing process of the panels was identical to the initial panel manufacturing process. No additional materials were added during the healing process. The healed PFRP panels were impacted again and substantially recovered energy absorption capability. We also conducted the repeated-impact test with several different impact energies. Unlike the CFRPs where the impact peak load decreased as the impacts repeated, the PFRPs showed an increasing trend. Such a unique mechanism was due to the strain-hardening behavior of the polymer fibers and matrix. As a result, the repeated-impact life of the PFRPs was significantly enhanced. These results are particularly interesting in the automotive or aerospace industries where repeated LVI is frequently observed.

Ko, Seunghyun↗
DOE OSTI2020

Machine learning coupled multi-scale modeling for redox flow batteries

The reaction distribution in macro or device-scale has been studied for redox flow batteries. The reaction distribution on electrode pore-scale structure however is not well understood, lacking especially on how the reaction distribution on the pore-scale may impact the overall performance of a flow battery. This study introduces for the first time a framework of a multi-scale model that provides understanding of the relationship between the pore-scale electrode structure reaction and the device-scale electrochemical reaction uniformity within the flow battery. A reduced order model is constructed based on 128 pore-scale simulations, which provide a quantitative relationship between the battery operation conditions (inlet velocity, current density, inlet concentration) and the surface reaction uniformity for the pore-scale sample. The multi-scale framework upscales this pore-scale surface reaction uniformity to device-scale combined uniformity. Based on the multi-scale model, a time-varying optimization of the inlet velocity is established, leading to significant reduction on pump power consumption with targeted surface reaction uniformity. The multi-scale model establishes the critical link between the micro-structure of a flow battery component and its performance at the macro-scale, therefore providing rationale for further operational or material optimization.

flow batteries, machine learning, multi-scale mode↗
DOE OSTI2023

Inversion-based correction of Double-Torsion (DT) subcritical crack growth tests for crack profile geometry

Because of its simplicity and the ability to produce a stable, slow-propagating crack, the Double-Torsion (DT) method has been used widely for investigating the critical and subcritical propagation of a slow-propagating tensile (mode-I) crack. However, to determine the complex relationship between the crack velocity $\mathcal{v_c}$ vs. the strain energy release rate $\mathscr{G}$ (or the stress intensity factor K) from laboratory measurements, several corrections must be made to account for the impact of sample and crack geometry. Particularly, DT test typically produces a crack with a curved edge profile instead of a straight line, causing the local $\mathcal{v_c}$ and $\mathscr{G}$ vary along the crack front. The experimentally measured $\mathcal{v_c}$ and $\mathscr{G}$ data merely reflect collective, averaged behavior of the crack. This makes inversion for the intrinsic, “true” crack growth kinetics necessary, based upon the knowledge of the crack geometry. Simple and effective correction methods have been proposed and validated for the slow, chemical-reaction-controlled part (Region I) of the $\mathcal{v_c}$-$\mathscr{G}$ curve. However, reliable methods for the highly nonlinear, transport-dominated part (Region II) and its sudden transition to the dynamic propagation part (Region III) are still lacking. Here we propose a method for determining the intrinsic $\mathcal{v_c}$-$\mathscr{G}$ relationship across all three Regions based upon DT test data, using a simple model function and its numerical inversion. The performance of this approach is examined and demonstrated using both synthetic and laboratory data for subcritical crack growth in soda lime glass.

42 ENGINEERING↗
DOE OSTI2022

Dissipation and Bathymetric Sensitivities in an Unstructured Mesh Global Tidal Model

Abstract The mechanisms and geographic distribution of global tidal dissipation in barotropic tidal models are examined using a high resolution unstructured mesh finite element model. Mesh resolution varies between 2 and 25 km and is especially focused on inner shelves and steep bathymetric gradients. Tidal response sensitivities to bathymetric changes are examined to put into context response sensitivities to frictional processes. We confirm that the Ronne Ice Shelf dramatically affects Atlantic tides but also find that bathymetry in the Hudson Bay system is a critical control. We follow a sequential frictional parameter optimization process and use TPXO9 data‐assimilated tidal elevations as a reference solution. From simulated velocities and depths, dissipation within the global model is estimated and allows us to pinpoint dissipation at high resolution. Boundary layer dissipation is extremely focused with 1.4% of the ocean accounting for 90% of the total. Internal tide friction is much more distributed with 16.7% of the ocean accounting for 90% of the total. Often highly regional dissipation can impact basin‐scale and even ocean wide tides. Optimized boundary layer friction parameters correlate very well with the physical characteristics of the locality with high friction factors associated with energetic tidal regions, deep ocean island chains, and ice covered areas. Global complex M 2 tide errors are 1.94 cm in deep waters. Total global boundary layer and internal tide dissipation are estimated, respectively, at 1.83 and 1.49 TW. This continues the trend in the literature toward attributing more dissipation to internal tides.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Monitoring Terrestrial Water Storage, Drought and Seasonal Changes in Central Oklahoma With Ambient Seismic Noise

Abstract Significant imbalances in terrestrial water storage (TWS) and severe drought have been observed around the world as a consequence of climate changes. Improving our ability to monitor TWS and drought is critical for water‐resource management and water‐deficit estimation. We use continuous seismic ambient noise to monitor temporal evolution of near‐surface seismic velocity, dv / v , in central Oklahoma from 2013 to 2022. The derived dv / v is found to be negatively correlated with gravitational measurements and groundwater depths, showing the impact of groundwater storage on seismic velocities. The hydrological effects involving droughts and recharge of groundwater occur on a multi‐year time scale and dominate the overall derived velocity changes. The thermoelastic response to atmospheric temperature variations occurs primarily on a yearly timescale and dominates the superposed seasonal velocity changes in this study. The occurrences of droughts appear simultaneously with local peaks of dv / v , demonstrating the sensitivity of near‐surface seismic velocities to droughts.

58 GEOSCIENCES↗
DOE OSTI2020

Advanced High-Performance Computational Modeling of the Seismic Response of High-Hazard and/or Nuclear Facilities and Critical Infrastructure at the NNSS

New methods for predicting the amplitude and variability of ground shaking from earthquakes (and explosions) are needed for seismic hazard analysis for buildings, nuclear power plants, and critical infrastructure at the NNSS. We are comparing existing 1-D and new 3-D geophysical methods for estimating the shear-wave velocity structure in the upper 30 meters of the ground surface (Vs30), which plays a major role in ground motion amplification and seismic response of buildings. We evaluate the performance of these methodologies at the U1a Complex at the NNSS and develop simple 1-D and high-resolution 3-D Vs30 models. We then emplace these high-resolution models into a background seismic velocity model. We will collaborate with Lawrence Livermore National Laboratory (LLNL) to conduct numerical modeling of the ground shaking at the NNSS using their high-performance computing technology and state-of-the-art ground motion simulation methodology. The primary work that was completed in FY 2019 was to acquire the seismic systems and familiarize staff at the NNSS with their use. We also worked on developing a collection plan with the Device Assembly Facility (DAF) at the NNSS, but due to time constraints and other ongoing projects at the DAF, we had to use U1a as a backup. We were able to coordinate the seismic survey, and we will complete the collection of seismic data in FY 2020. Additionally, during FY 2019, we completed the geologic framework model (GFM) for the U1a Complex and modeled the Yucca fault. LLNL worked with us through FY 2019 to prepare the data files for their modeling software and tested the software for reliability. In FY 2020 we will develop the end-to-end capability so that any facility could easily be modeled and the expected shaking from a local earthquake understood. The work in FY 2020 will include building fault models from the GFM and finalizing the velocity model analysis. The final simulations will be run for multiple rupture models, and final assessments will demonstrate the seismic hazard at the U1a Complex.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗