Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial statistics”

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

From bulk to surface: Structure and dynamics of amorphous alumina from deep potential molecular dynamics

Understanding the atomic-scale structure and dynamics of amorphous oxide surfaces is essential for interpreting their chemical reactivity, mechanical stability, and interfacial behavior, yet direct experimental characterization remains challenging. We employ Deep Potential (DP) molecular dynamics to generate large-scale, ab initio -quality models of amorphous Al 2 O 3 bulk glasses and melt-quenched free surfaces, enabling a quantitative analysis of both structure and relaxation dynamics with statistical confidence inaccessible to direct ab initio simulation. The trained DP model reproduces experimental liquid and glass structure, captures the cooling-rate dependence of the bulk glass transition, and corrects systematic biases in the polyhedral populations predicted by widely used classical force fields. At the free surface, mass density recovers to bulk values over ~10 Å, while local coordination requires a slightly wider subsurface region to fully converge. The outermost layer is oxygen-enriched, exhibits altered polyhedral connectivity with contracted Al–O bonds, and hosts a broad population of under-coordinated motifs (notably AlO 3 and OAl 2 ) whose abundances are governed by glass stability. These under-coordinated surface motifs exhibit distinct vibrational signatures and occur as locally paired Lewis acid and Brønsted base sites consistent with bond-valence compensation, yet remain spatially dispersed rather than aggregating into extended clusters. Despite this pronounced structural heterogeneity, surface relaxation and the glass-transition temperature remain comparable to their bulk counterparts, suggesting that the disordered surface is kinetically stable once formed. Together, these results establish a molecular-level picture of amorphous alumina surfaces and demonstrate the capability of machine-learned potentials to resolve structure–property relationships in disordered oxide interfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Hydrological projections over the Upper Indus Basin at 1.5 °C and 2.0 °C temperature increase

We analyse an ensemble of statistically downscaled Global Climate Models (GCMs) to investigate future water availability in the Upper Indus Basin (UIB) of Pakistan for the time horizons when the global and/or regional warming levels cross Paris Agreement (PA) targets. The GCMs data is obtained from the 5th Phase of Coupled Model Inter-Comparison Project under two Representative Concentration Pathways (RCP4.5 and RCP8.5). Based on the five best performing GCMs, we note that global 1.5 °C and 2.0 °C warming thresholds are projected in 2026 and 2047 under RCP4.5 and 2022 and 3036 under RCP8.5 respectively while these thresholds are reached much earlier over Pakistan i.e. 2016 and 2030 under RCP4.5 and 2012 and 2025 under RCP8.5 respectively. Interestingly, the GCMs with the earliest emergence at the global scale are not necessarily the ones with the earliest emergence over Pakistan, highlighting spatial non-linearity in GCMs response. The emergence of 2.0 °C warming at global scale across 5 GCMs ranges from 2031 (CCSM4) to 2049 (NorESM) under RCP8.5. Precipitation generally exhibits a progressive increasing trend with stronger changes at higher warming or radiative forcing levels. Hydrological simulations representing the historical, 1.5 °C and 2.0 °C global and region warming time horizons indicate a robust but seasonally varying increase in the inflows. The highest inflows in the baseline and future are witnessed in July. However, the highest future increase in inflows is projected in October under RCP4.5 (37.99% and 65.11% at 1.5 °C and 2.0 °C) and in April under RCP8.5 (37% and 62.05% at 1.5 °C and 2.0 °C). These hydrological changes are driven by increases in the snow and glacial melt contribution, which are more pronounced at 2.0 °C warming level. These findings should help for effective water management in Pakistan over the coming decades.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Transitions in the morphology and critical stresses of gliding dislocations in multiprincipal element alloys

Refractory multiprincipal element alloys (MPEAs) are promising material candidates for high-temperature, high-strength applications. However, their deformation mechanisms, particularly at the dislocation scale, are unusual compared to those of conventional alloys, making it challenging to understand the origin of their high strength. Here, using an atomistically informed phase-field dislocation dynamics model, we study transitions in the morphology and critical stresses of long gliding screw dislocations over extended distances in a refractory MPEA. The model MPEA crystal accounts for the atomic-scale fluctuations in chemical composition across the glide planes via spatially correlated lattice energies and for the differences in glide resistances between screw and edge dislocations of unit length. Here, we show that the dislocation moves in a stop-start motion, alternating between wavy morphology in free flight and nearly recovered straight screw orientation in full arrest. The periods of wavy glide are due to variable kink-pair formation and migration rates along the length, where portions with higher rates glide more quickly. The critical stress to initiate motion corresponds to the stress required to form and migrate a kink pair at the weakest region along the length of the dislocation. Heterogeneity in lattice energy leads to variability in the local stress-strain response and to a strain hardening-like response, in which the critical stress to reactivate glide increases with glide distance. Statistical assessment of hundreds of realizations of dislocations indicates that the amount of hardening directly scales with the dispersion in underlying lattice energy.

3-dimensional systems↗
DOE OSTI2022

Imaging Phase Segregation in Nanoscale Li x CoO 2 Single Particles

Li x CoO 2 (LCO) is a common battery cathode material that has recently emerged as a promising material for other applications including electrocatalysis and as electrochemical random access memory (ECRAM). During charge– discharge cycling LCO exhibits phase transformations that are significantly complicated by electron correlation. While the bulk phase diagram for an ensemble of battery particles has been studied extensively, it remains unclear how these phases scale to nanometer dimensions and the effects of strain and diffusional anisotropy at the single-particle scale. Understanding these effects is critical to modeling battery performance and for predicting the scalability and performance of electrocatalysts and ECRAM. Here we investigate isolated, epitaxial LiCoO 2 islands grown by pulsed laser deposition. After electrochemical cycling of the islands, conductive atomic force microscopy (c-AFM) is used to image the spatial distribution of conductive and insulating phases. Above 20 nm island thicknesses, we observe a kinetically arrested state in which the phase boundary is perpendicular to the Li-planes; we propose a model and present image analysis results that show smaller LCO islands have a higher conductive fraction than larger area islands, and the overall conductive fraction is consistent with the lithiation state. Thinner islands (14 nm), with a larger surface to volume ratio, are found to exhibit a striping pattern, which suggests surface energy can dominate below a critical dimension. When increasing force is applied through the AFM tip to strain the LCO islands, significant shifts in current flow are observed, and underlying mechanisms for this behavior are discussed. The c-AFM images are compared with photoemission electron microscopy images, which are used to acquire statistics across hundreds of particles. Finally, the results indicate that strain and morphology become more critical to electrochemical performance as particles approach nanometer dimensions.

25 ENERGY STORAGE↗
DOE OSTI2021

Examining the Impacts of Great Lakes Temperature Perturbations on Simulated Precipitation in the Northeastern United States

Most inland water bodies are not resolved by general circulation models, requiring that lake surface temperatures be estimated. Given the large spatial and temporal variability of the surface temperatures of the North American Great Lakes, such estimations can introduce errors when used as lower boundary conditions for dynamical downscaling. Lake surface temperatures (LSTs) influence moisture and heat fluxes, thus impacting precipitation within the immediate region and potentially in regions downwind of the lakes. For this study, the Advanced Research version of the Weather Research and Forecasting Model (WRF-ARW) was used to simulate precipitation over the six New England states during a 5-yr historical period. The model simulation was repeated with perturbed LSTs, ranging from 10°C below to 10°C above baseline values obtained from reanalysis data, to determine whether the inclusion of erroneous LST values has an impact on simulated precipitation and synoptic-scale features. Results show that simulated precipitation in New England is statistically correlated with LST perturbations, but this region falls on a wet–dry line of a larger bimodal distribution. Wetter conditions occur to the north and drier conditions occur to the south with increasing LSTs, particularly during the warm season. Additionally, the precipitation differences coincide with large-scale anomalous temperature, pressure, and moisture patterns. Care must therefore be taken to ensure reasonably accurate Great Lakes surface temperatures when simulating precipitation, especially in southeastern Canada, Maine, and the mid-Atlantic region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Scaling Analysis of Two–Phase Flow in Fractal Permeability Fields

Fluid mixing in permeable media is essential in many practical applications. The mixing process is a consequence of velocity fluctuations owing to geological heterogeneities and mobility contrast of fluids. Heterogeneities in natural rocks are often spatially correlated, and their properties, such as permeability, may be described using fractal distributions. This work models the fractal characteristics of such permeability fields in which the covariance function is expressed as a power-law function. A generalized scaling relation is derived relating various fractal permeability fields using the magnitude of their fluctuations. Here, this relation reveals the self-similar behavior of two-phase flow in such permeable media. To that end, a recently developed, high-resolution numerical simulator is employed to validate the analytically derived scaling relations. Two flow problems are considered in which flow is governed by 1) a linear, and 2) a nonlinear transport equation. Due to the probabilistic representation of the fractal permeability fields, a sensitivity study is conducted for each flow scenario to determine the number of realizations required for statistical convergence. Scaling analysis is performed using ensemble averages of simulated saturation profiles and their mixing lengths. Results support the validity of the developed scaling relation across the range of investigated flow conditions at intermediate times. The dynamics of linear flow in the asymptotic regime is affected by the correlation structure of heterogeneity. In nonlinear flow, scaling behavior appears to be dominated by the degree of nonlinearity.

58 GEOSCIENCES↗
DOE OSTI2025

Videos, photos, and AI-derived grain size data associated with “High-throughput AI Video Surveys Enable Reproducible Multiscale Sediment Size Mapping, with Implications for Hydrobiogeochemical Parameterization”

NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the manuscript “High-throughput AI Video Surveys Enable Reproducible Multiscale Sediment Size Mapping, with Implications for Hydrobiogeochemical Parameterization” under review. This data package includes five data types: 1) raw photos and videos from drone survey and walking smartphone surveys; 2) images derived from raw videos; 3) manual labeling of reference scales; 4) metadata for all images and photo resolution derived from artificial intelligence (AI) models or manual labels, 5) grain size data obtained from AI models for all photos, 6) metadata and grain size data after quality control, 7) summaries of sample efficiency for all data, and 8) computational fluid dynamics (CFD) data used to support hydro-biogeochemical (HBGC) parameter estimation. Such data is used to 1) demonstrate significant improvements in accuracy, efficiency, and quality control for grain size data collection with the help of AI models, 2) study the spatial heterogeneity of grain size and observation reproducibility based on tens of thousands of data points generated by the AI models, and 3) evaluate the impacts of grain size heterogeneity on key HBGC parameters across sediment-to-reach and hourly-to-yearly scales. In particular, the data package contains 116 folders and 179696 files. The files include 41 videos in .mov format, 64047 photos in .jpg format, 13541 video-derived photos in .png format, 12747 segmentation mask data in .tif format, 12747 segmentation data in .json format, 24771 .csv files that with metadata and grain size for each individual photo as well as water depth and velocity data from CFD and observation, 51791 .txt files of raw AI predicted labels, and 11 flight record data in .srt format. The summary for all metadata and grain size statistics information is included in “Scales_V3_NG.csv” and “Statistics_V3_NG.csv”. The summary for data that pass data quality control (QC) level 0-2 is included in “QCStatistics_V3_NG.csv”. The QC level 0 represents photos whose photo resolution is positive, excluding photos that miss reference scale. The QC level 1 means reference scale circularity uncertainty is less than 5% for smartphone images while representing photo resolution is larger than 0.44 mm/pixel for drone images. The QC level 2 means excluding photos whose grain number is less than 100, a minimum number of grains recommended by classic literature. The summary for each video’s name, length, frame rates, survey area, grain number, survey efficiency, etc. can be found in “QCSummary_V3_NG.csv”. The summary for site name, GPS coordinates, and number of images at each site can be found in “SitesSummary_V3_*.csv” files. Overall computational efficiency summary is reported in Table 4 of accompanying manuscript. Additionally, the nitrate concentration data used in this work was downloaded from an existing dataset published on ESS-DIVE (Boat-Dragged Sensor Hanford Reach.csv; Conner A. et al., 2020). We thank the United States Forest Service, Washington Department of Fish and Wildlife, Washington Department of Natural Resources, Cowiche Canyon Conservatory, Port of Benton, and the Confederated Tribes and Bands of the Yakama Nation for access to field locations where the data were collected. We also thank the Yakama Nation Tribal Council and Yakama Nation Fisheries for working with us to facilitate data collection and optimization of data usage according to their values and worldview.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

ELM2.1-XGBfire1.0: improving wildfire prediction by integrating a machine learning fire model in a land surface model

Wildfires have shown increasing trends in both frequency and severity across the contiguous United States (CONUS). However, process-based fire models have difficulties in accurately simulating the burned area over the CONUS due to a simplification of the physical process and cannot capture the interplay among fire, ignition, climate, and human activities. The deficiency of burned area simulation deteriorates the description of fire impact on energy balance, water budget, and carbon fluxes in the Earth system models (ESMs). Alternatively, fire models based on machine learning (ML), which capture statistical relationships between the burned area and environmental factors, have shown promising burned area predictions and corresponding fire impact simulation. We develop a hybrid framework (ELM2.1-XGBFire1.0) that integrates an eXtreme Gradient Boosting (XGBoost) wildfire model with the Energy Exascale Earth System Model (E3SM) land model (ELM) version 2.1. A Fortran–C–Python deep learning bridge is adapted to support online communication between ELM and the ML fire model. Specifically, the burned area predicted by the ML-based wildfire model is directly passed to ELM to adjust the carbon pool and vegetation dynamics after disturbance, which are then used as predictors in the ML-based fire model in the next time step. Evaluated against the historical burned area from Global Fire Emissions Database 5 from 2001–2019, the ELM2.1-XGBFire1.0 outperforms process-based fire models in terms of spatial distribution and seasonal variations. The ELM2.1-XGBFire1.0 has proven to be a new tool for studying vegetation–fire interactions and, more importantly, enables seamless exploration of climate–fire feedback, working as an active component of E3SM.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Simulated wildfire burned area over the CONUS during 2001-2020

Wildfires have shown increasing trends in both frequency and severity across the Contiguous United States (CONUS). However, process-based fire models have difficulties in accurately simulating the burned area over the CONUS due to a simplification of the physical process and cannot capture the interplay among fire, ignition, climate, and human activities. The deficiency of burned area simulation deteriorates the description of fire impact on energy balance, water budget, and carbon fluxes in the Earth System Models (ESMs). Alternatively, machine learning (ML) based fire models, which capture statistical relationships between the burned area and environmental factors, have shown promising burned area predictions and corresponding fire impact simulation. We develop a hybrid framework (ML4Fire-XGB) that integrates a pretrained eXtreme Gradient Boosting (XGBoost) wildfire model with the Energy Exascale Earth System Model (E3SM) land model (ELM). A Fortran-C-Python deep learning bridge is adapted to support online communication between ELM and the ML fire model. Specifically, the burned area predicted by the ML-based wildfire model is directly passed to ELM to adjust the carbon pool and vegetation dynamics after disturbance, which are then used as predictors in the ML-based fire model in the next time step. Evaluated against the historical burned area from Global Fire Emissions Database 5 from 2001-2020, the ML4Fire-XGB model outperforms process-based fire models in terms of spatial distribution and seasonal variations. Sensitivity analysis confirms that the ML4Fire-XGB well captures the responses of the burned area to rising temperatures. The ML4Fire-XGB model has proved to be a new tool for studying vegetation-fire interactions, and more importantly, enables seamless exploration of climate-fire feedback, working as an active component in E3SM.

Liu, Ye↗
DOE OSTI2021

DSGAN

This study develops a neural network-based approach for emulating high-resolution modeled precipitation data with comparable statistical properties but at greatly reduced computational cost. The key idea is to use combination of low-and high- resolution simulations to train a neural network to map from the former to the latter. Specifically, we define two types of CNNs, one that stacks variables directly and one that encodes each variable before stacking, and we train each CNN type both with a conventional loss function, such as mean square error (MSE), and with a conditional generative adversarial network (CGAN), for a total of four CNN variants. We compare the four new CNN-derived high-resolution precipitation results with precipitation generated from original high resolution simulations, a bilinear interpolater and the state-of-the-art CNN-based super-resolution (SR) technique. Results show that the SR technique produces results similar to those of the bilinear interpolator with smoother spatial and temporal distributions and smaller data variabilities and extremes than the original high resolution simulations. While the new CNNs trained by MSE generate better results over some regions than the interpolator and SR technique do, their predictions are still not as close as the original high resolution simulations. The CNNs trained by CGAN generate more realistic and physically reasonable results, better capturing not only data variability in time and space but also extremes such as intense and long-lasting storms. The new proposed CNN-based downscaling approach can downscale precipitation from 50~km to 12~km in 14~min for 30~years

LIU, ZHENGCHUN↗
DOE OSTI2022

Performance of borated scintillator screens for high-resolution neutron imaging

The most commonly used screens for neutron imaging consist of 6 LiF + ZnS. This type of screen yields the highest light output per detected neutron. For high resolution, gadolinium oxysulfide (GOS, Gadox) screens are employed, which have a much higher detection efficiency, but a light output so much lower than LiF + ZnS that measurements are often limited by photon statistics. Historically, screens using boron as a neutron-sensitive material have not been very successful. However, a new preparation method was introduced recently that produces light output higher than Gadox with detection efficiency greater than LiF + ZnS. Measurements of these new borated screens were performed at the NeXT facility at ILL, Grenoble, in comparison to a high resolution Gadox screen.

10B↗
DOE OSTI2020

Thermal transport exceeding bulk heat conduction due to nonthermal micro/nanoscale phonon populations

While classical size effects usually lead to a reduced effective thermal conductivity, we report here that nonthermal phonon populations produced by a micro/nanoscale heat source can lead to enhanced heat conduction, exceeding the prediction from Fourier's law. We study nondiffusive thermal transport by phonons at small distances within the framework of the Boltzmann transport equation (BTE) and demonstrate that the transport is significantly affected by the distribution of phonons emitted by the source. We discuss analytical solutions of the steady-state BTE for a source with a sinusoidal spatial profile, as well as for a three-dimensional Gaussian “hot spot,” and provide numerical results for single crystal silicon at room temperature. If a micro/nanoscale heat source produces a thermal phonon distribution, it gets hotter than that predicted by the heat diffusion equation; however, if the source predominantly produces low-frequency acoustic phonons with long mean free paths, it may get significantly cooler than that predicted by the heat equation, yielding an enhanced heat transport beyond bulk heat conduction.

42 ENGINEERING↗
DOE OSTI2020

Grid-based diffusion Monte Carlo for fermions without the fixed-node approximation

A diffusion Monte Carlo algorithm is introduced that can determine the correct nodal structure of the wave function of a few-fermion system and its ground-state energy without an uncontrolled bias. This is achieved by confining signed random walkers to the points of a uniform infinite spatial grid, allowing them to meet and annihilate one another to establish the nodal structure without the fixed-node approximation. An imaginary-time propagator is derived rigorously from a discretized Hamiltonian, governing a non-Gaussian, sign-flipping, branching, and mutually annihilating random walk of particles. The accuracy of the resulting stochastic representations of a fermion wave function is limited only by the grid and imaginary-time resolutions and can be improved in a controlled manner. Here, the method is tested for a series of model problems including fermions in a harmonic trap as well as the He atom in its singlet or triplet ground state. For the latter case, the energies approach from above with increasing grid resolution and converge within 0.015 E h of the exact basis-set-limit value for the grid spacing of 0.08 a.u. with a statistical uncertainty of 10 –5 E h without an importance sampling or Jastrow factor.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Improving North American Wildfire Prediction by Integrating a Machine-Learning Fire Model in a Land Surface Model

Wildfires have shown increasing trends in both frequency and severity across the Contiguous United States (CONUS). However, process-based fire models have difficulties in accurately simulating the burned area over the CONUS due to a simplification of the physical process and cannot capture the interplay among fire, ignition, climate, and human activities. The deficiency of burned area simulation deteriorates the description of fire impact on energy balance, water budget, and carbon fluxes in the Earth System Models (ESMs). Alternatively, machine learning (ML) based fire models, which capture statistical relationships between the burned area and environmental factors, have shown promising burned area predictions and corresponding fire impact simulation. We develop a hybrid framework (ML4Fire-XGB) that integrates a pretrained eXtreme Gradient Boosting (XGBoost) wildfire model with the Energy Exascale Earth System Model (E3SM) land model (ELM) version 2.1. A Fortran-C-Python deep learning bridge is adapted to support online communication between ELM and the ML fire model. Specifically, the burned area predicted by the ML-based wildfire model is directly passed to ELM to adjust the carbon pool and vegetation dynamics after disturbance, which are then used as predictors in the ML-based fire model in the next time step. Evaluated against the historical burned area from Global Fire Emissions Database 5 from 2001-2020, the ML4Fire-XGB model outperforms process-based fire models in terms of spatial distribution and seasonal variations. Sensitivity analysis confirms that the ML4Fire-XGB well captures the responses of the burned area to rising temperatures. The ML4Fire-XGB model has proved to be a new tool for studying vegetation-fire interactions, and more importantly, enables seamless exploration of climate-fire feedback, working as an active component in E3SM.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Role of Chemical Disorder in High Temperature Dislocation Glide in Refractory Multi-principal Element Alloys

Refractory multi-principal element alloys (RMPEAs) combine a chemically disordered lattice with a structurally ordered, single‐phase body‐centered cubic (bcc) crystal structure. Chemical fluctuations in these alloys give rise to significant energy barriers that impede dislocation motion. In this study, we use phase field dislocation dynamics to examine how spatial temperature fluctuations compete with randomness in energy barriers to affect the motion of long screw dislocations in three equi-atomic MoNbTa‐based RMPEAs: MoNbTa, MoNbTaW, and MoNbTaVW. Over a wide range of homologous temperatures (T h ≈ 0–0.6), we determined a screw dislocation ‘flow stress’ as the minimum applied stress to sustain continuous motion over a long excursion distance within a fixed timeframe. All three RMPEAs exhibit temperature dependent flow stress with three characteristic glide regimes: at low homologous temperatures, flow stress drops sharply, and screw glide remains planar and rectilinear; at intermediate homologous temperatures, the flow stress levels off as glide becomes planar but wavy; and at high homologous temperatures, flow stress plateaus as screws exhibit nonplanar, three‐dimensional motion. Glide kinetics and transition temperatures are controlled by the chemically induced fluctuations in the energy landscape. The rectilinear to wavy transition temperature is controlled by statistically weakest local barriers in the glide plane, whereas the wavy to 3D transition temperature is governed by statistically strongest local barriers. At low and intermediate homologous temperatures, the relative spread in energy barriers governs glide behavior by controlling the local kinetics of kink pair formation and kink pinning. At high homologous temperatures, the average barrier height governs the glide behavior by controlling the number of out-of-plane excursions during 3D glide. These findings reveal how random chemical fluctuations determine screw‐driven plasticity in RMPEAs, providing critical insight for the design of high temperature structural alloys.

Defects↗
DOE OSTI2023

Continental Scale Hydrostratigraphy: Comparing Geologically Informed Data Products to Analytical Solutions

Abstract This study synthesizes two different methods for estimating hydraulic conductivity (K) at large scales. We derive analytical approaches that estimate K and apply them to the contiguous United States. We then compare these analytical approaches to three‐dimensional, national gridded K data products and three transmissivity (T) data products developed from publicly available sources. We evaluate these data products using multiple approaches: comparing their statistics qualitatively and quantitatively and with hydrologic model simulations. Some of these datasets were used as inputs for an integrated hydrologic model of the Upper Colorado River Basin and the comparison of the results with observations was used to further evaluate the K data products. Simulated average daily streamflow was compared to daily flow data from 10 USGS stream gages in the domain, and annually averaged simulated groundwater depths are compared to observations from nearly 2000 monitoring wells. We find streamflow predictions from analytically informed simulations to be similar in relative bias and Spearman's rho to the geologically informed simulations. R ‐squared values for groundwater depth predictions are close between the best performing analytically and geologically informed simulations at 0.68 and 0.70 respectively, with RMSE values under 10 m. We also show that the analytical approach derived by this study produces estimates of K that are similar in spatial distribution, standard deviation, mean value, and modeling performance to geologically‐informed estimates. The results of this work are used to inform a follow‐on study that tests additional data‐driven approaches in multiple basins within the contiguous United States.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Simulating the filamentation of smoothed laser beams with three-dimensional nonlinear dynamics

In a plasma, the ponderomotive force of an inhomogeneous electromagnetic field expels plasma from regions of high intensity. When a laser propagates through a plasma, this force creates density wells that subsequently modify the index of refraction. The beam refracts and focuses into these wells and may filament. In extreme cases, the laser beam will spray due to increasing angular divergence of the beam. The threshold for ponderomotive self-focusing is well established for isolated laser hotspots or speckles. Here, we define a practical threshold for characterizing the filamentation of thousands of speckles that are found in the focal plane of high-power laser beams spatially smoothed with random phase plates as used at high energy and power laser facilities studying inertial confinement fusion. This threshold is tested against three-dimensional simulations of speckled laser light propagating through plasma. Four metrics are applied to assess filamentation: the fraction of power above five times the average intensity, an effective f-number, the mean-squared perpendicular wavenumber, and the fraction of rarefied density with deviation from the initial density exceeding | δn/ n| = 0.1. The speckled beams studied are generated by random phase plates, both with and without additional polarization smoothing, in a parameter regime of relevance to indirect drive experiments. While filamentation has been discussed extensively in the literature, we believe this to be the first published simulation study with three-dimensional nonlinear hydrodynamics that addresses the onset threshold of ponderomotive filamentation and establishes the lengths and time scales necessary to reach a statistical steady state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Measuring cosmological distances using cluster edges as a standard ruler

ABSTRACT The line-of-sight velocity dispersion profile of galaxy clusters exhibits a ‘kink’ corresponding to the spatial extent of orbiting galaxies. Because the spatial extent of a cluster is correlated with the amplitude of the velocity dispersion profile, we can utilize this feature as a gravity-calibrated standard ruler. Specifically, the amplitude of the velocity dispersion data allows us to infer the physical cluster size. Consequently, observations of the angular scale of the ‘kink’ in the profile can be translated into a distance measurement to the cluster. Assuming the relation between cluster radius and cluster velocity dispersion can be calibrated from simulations, we forecast that with existing data from the Sloan Digital Sky Survey we will be able to measure the Hubble constant with 3.0 per cent precision. Implementing our method with data from the Dark Energy Spectroscopic Instrument (DESI) will result in a 1.3 per cent measurement of the Hubble constant. Adding cosmological supernova data improves the uncertainty of the DESI measurement to 0.7 per cent. While these error estimates are statistical only, they provide strong motivation for pursuing the necessary simulation program required to characterize and calibrate the systematic uncertainties impacting our proposed measurement. Whether or not our proposed measurement can in fact result in competitive H0 constraints will depend on what the eventual systematics floor for this method is.

79 ASTRONOMY AND ASTROPHYSICS↗