SEARCH · Search NASA

Results for “correlation model”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21

All sources NASA NTRS DOE OSTI

DOE OSTI2024

On the statistical theory of self-gravitating collisionless dark matter flow: Scale and redshift variation of velocity and density distributions

The statistics of velocity and density fields are crucial for cosmic structure formation and evolution. Here, this paper extends our previous work on the two-point second-order statistics for the velocity field [Phys. Fluids 35, 077105 (2023)] to one-point probability distributions for both density and velocity fields. The scale and redshift variation of density and velocity distributions are studied by a halo-based non-projection approach. First, all particles are divided into halo and out-of-halo particles so that the redshift variation can be studied via generalized kurtosis of distributions for halo and out-of-halo particles, respectively. Second, without projecting particle fields onto a structured grid, the scale variation is analyzed by identifying all particle pairs on different scales $r$. We demonstrate that: (i) Delaunay tessellation can be used to reconstruct the density field. The density correlation, spectrum, and dispersion functions were obtained, modeled, and compared with the N-body simulation; (ii) the velocity distributions are symmetric on both small and large scales and are non-symmetric with a negative skewness on intermediate scales due to the inverse energy cascade on small scales with a constant rate $\varepsilon_u$; (iii) On small scales, the even order moments of pairwise velocity $\Delta u_L$ follow a two-thirds law $\propto{(-\varepsilon_ur)}^{2/3}$, while the odd order moments follow a linear scaling $\langle(\Delta u_L)^{2n+1}\rangle=(2n+1)\langle(\Delta u_L)^{2n}\rangle\langle\Delta u_L\rangle\propto{r}$; (iv) The scale variation of the velocity distributions was studied for longitudinal velocities $u_L$ or $u_L^{'}$, pairwise velocity (velocity difference) $\Delta u_L$=$u_L^{'}$-$u_L$ and velocity sum $\Sigma u_L$=$u^{'}_L$+$u_L$. Fully developed velocity fields are never Gaussian on any scale, despite that they can initially be Gaussian; (v) On small scales, $u_L$ and $\Sigma u_L$ can be modeled by a $X$ distribution to maximize the entropy of the system. The distribution of $\Delta u_L$ can be different; (vi) On large scales, $\Delta u_L$ and $\Sigma u_L$ can be modeled by a logistic or a $X$ distribution, while $u_L$ has a different distribution; (vii) the redshift variation of the velocity distributions follows the evolution of the $X$ distribution involving a shape parameter $\alpha(z)$ decreasing with time.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Testing effective field theory with the most general neutron decay correlations

Built on the seminal works by Jackson-Treiman-Wyld and Ebel-Feldman, we derive the most general free neutron differential decay rate where all massive particles (neutron, proton, and electron) are polarized. This introduces 33 new correlations in addition to the 18 existing ones, which overconstrain the coupling constants in the low-energy effective field theory of charged weak interactions, and thus provides stringent tests of the validity of the theory framework itself. We classify the correlation coefficients in terms of their Standard Model limit and discrete symmetries, and study their expansion with respect to the new physics coupling strengths, supplemented by the experiment-independent 𝒪⁡(𝛼) virtual electromagnetic radiative corrections.

Astronomy & Astrophysics↗

DOE OSTI2025

Light and Strange Vector Resonances from Lattice QCD at Physical Quark Masses

We present the first ab initio calculation at physical quark masses of scattering amplitudes describing the lightest pseudoscalar mesons interacting via the strong force in the vector channel. Using lattice quantum chromodynamics, we postdict the defining parameters for two short-lived resonances, the ρ(770) and K*(892) which manifest as complex energy poles in ππ and Kπ scattering amplitudes, respectively. The calculation proceeds by first computing the finite-volume energy spectrum of the two-hadron systems and then determining the amplitudes from the energies using the Lüscher formalism. The error budget includes a data-driven systematic error, obtained by scanning possible fit ranges and fit models to extract the spectrum from Euclidean correlators, as well as the scattering amplitudes from the latter. The final results, obtained by analytically continuing multiple parametrizations into the complex energy plane, are M ρ = 796(5)(50)MeV, Γ ρ = 192(10)(31) MeV, M K* = 893(2)(54) MeV, and Γ K = 51(2)(11) MeV, where the subscript indicates the resonance and M and Γ stand for the mass and width, respectively, and where the first bracket indicates the statistical and the second bracket the systematic uncertainty.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Model Residuals as Shields: A Two-Level Formulation to Defend Smart Grids From Poisoning Attacks

The advancement of smart grids presents both vast opportunities and heightened cybersecurity risks. Data-driven defense mechanisms, though designed as a shield against these threats, can fall prey to poisoning attacks. We delve into regression settings, underscoring the imperative to fortify defenses against a spectrum of poison ratios, notably those above 0.5—an issue scarcely addressed in prior studies. Recognizing the susceptibilities of smart grids and their manipulable sensors, we exploit the very intent of poisoning attacks, compromising model accuracy, as our defense mechanism. Our proposed two-level optimization framework discerns between poisoned and authentic data based on model residuals, outperforming or matching existing methods in 72% to 77% of precision and 75% to 80% of recalls across various poisoning attacks, poison ratios, and datasets. Once the authentic data are identified, the trained model is adaptable for a variety of applications. Comprehensive evaluations on different smart grid datasets, pitted against myriad poisoning schemes, validate our methodology’s edge over existing methods. Here, we also shed light on the implications of model misspecification originating from temporal auto-correlation, a common feature in Internet of Things and smart grid data.

Adversarial machine learning (ML)↗

DOE OSTI2024

Day-to-day reliability of basal heart rate and short-term and ultra short-term heart rate variability assessment by the Equivital eq02+ LifeMonitor in US Army soldiers

Introduction The present study determined the (1) day-to-day reliability of basal heart rate (HR) and HR variability (HRV) measured by the Equivital eq02+ LifeMonitor and (2) agreement of ultra short-term HRV compared with short-term HRV. Methods Twenty-three active-duty US Army Soldiers (5 females, 18 males) completed two experimental visits separated by >48 hours with restrictions consistent with basal monitoring (eg, exercise, dietary), with measurements after supine rest at minutes 20–21 (ultra short-term) and minutes 20–25 (short-term). HRV was assessed as the SD of R–R intervals (SDNN) and the square root of the mean squared differences between consecutive R–R intervals (RMSSD). Results The day-to-day reliability (intraclass correlation coefficient (ICC)) using linear-mixed model approach was good for HR (0.849, 95% CI: 0.689 to 0.933) and RMSSD (ICC: 0.823, 95% CI: 0.623 to 0.920). SDNN had moderate day-to-day reliability with greater variation (ICC: 0.689, 95% CI: 0.428 to 0.858). The reliability of RMSSD was slightly improved when considering the effect of respiration (ICC: 0.821, 95% CI: 0.672 to 0.944). There was no bias for HR measured for 1 min versus 5 min (p=0.511). For 1 min measurements versus 5 min, there was a very modest mean bias of −4 ms for SDNN and −1 ms for RMSSD (p≤0.023). Conclusion When preceded by a 20 min stabilisation period using restrictions consistent with basal monitoring and measuring respiration, military personnel can rely on the eq02+ for basal HR and RMSSD monitoring but should be more cautious using SDNN. These data also support using ultra short-term measurements when following these procedures.

General & Internal Medicine↗

DOE OSTI2025

STARLIGHT

A Python package consisting of standard and novel methods for calculating electronic correlation effects and electrical conductivity in lattice models subject to strong magnetic fields.

Staros, Daniel [Los Alamos National Laboratory]↗

DOE OSTI2024

Combined ambient vibration and surface displacement measurements for improved progressive failure monitoring at a toppling rock slab in Utah, USA

Seismic resonance and surface displacement measurements can be implemented in tandem to improve landslide characterization and progressive failure monitoring. Crack aperture data are frequently used in rock slope stability monitoring and often exhibit recognizable trends prior to failure, such as accelerated crack opening. Alternatively, ambient resonance data offer multiple parameters including modal frequencies, damping, and polarization that can be monitored alongside crack aperture and may respond differently to environmental forcings and complex failure evolution. We analyzed data from continuous ambient vibration monitoring and concomitant crack aperture measurements at the Courthouse Mesa instability, a large toppling sandstone slab in Utah, USA. Three years of data revealed crack aperture increases of 2–4 mm/year with no clearly detectable irreversible changes in modal parameters, including frequency. Annually, frequency and displacement varied by 29% and 19% of the mean, respectively, with average and maximum daily frequency fluctuations of 6.5% and 16%, respectively. These reversible cyclic changes were primarily temperature-driven, but annually, frequency was in-phase with temperature whereas crack aperture lagged temperature changes by ~37 days. Polarization and damping also varied seasonally but were less strongly correlated with temperature. Conceptual 3D finite element modeling demonstrated consistent frequency decreases associated with crack propagation but variable changes in crack aperture measured at a single point; i.e., crack propagation did not always result in increased crack opening but always generated a resonance frequency decrease. Taken together, our data suggest a possible thermal wedging-ratcheting mechanism at the Courthouse Mesa instability, where annual thermoelastic crack closure is impeded by debris infill but the absence of downward crack propagation during the monitoring period is evidenced by no permanent resonance frequency changes. Our study demonstrates that combined seismic resonance and crack aperture data provide an improved description of rock slope instability behavior, supporting refined characterization and monitoring of changes accompanying progressive failure.

58 GEOSCIENCES↗

DOE OSTI2026

Measuring local primordial non-Gaussianity from the clustering of DESI DR1 LRGs and QSOs

We report the first measurement of primordial non-Gaussianity (PNG), parameterized by $f_{\mathrm{NL}}$, in the configuration space two-point correlation function (2pcf). We employ simulation based modeling and a novel approach for the mitigation of imaging systematics. We apply this method to samples of luminous red galaxies (LRG) and quasars (QSO) observed by the Dark Energy Spectroscopic Instrument (DESI) during the first year of its observations (DR1). The observed 68% CL interval on $f_{\mathrm{NL}}$ is $-3^{+22}_{-21}$ using LRGs, and $ 0^{+17}_{-16}$ using QSOs. The joint measurement yields $f_{\mathrm{NL}} = -3^{+12}_{-12}$ at $\ [68\%]$ CL. Our pipeline imposes a Gaussian prior on the value of $p$ (which defines the PNG bias via the Universality relation), with $p_{\rm LRG} = 1.0 \pm 0.1$ and $p_{\rm QSO} = 1.6\pm 0.1$. The observed constraining power of DESI tracers significantly exceeds that of previous large-scale structure (LSS) surveys, and encouragingly, approaches the sensitivity of CMB probes of PNG.

Brown, Z. [Kansas State U.]↗

DOE OSTI2024

Exploring the impact of multivalent substitution on high-temperature electrical conductivity in doped lanthanum chromite perovskites: An experimental and ab-initio study

Doped LaCrO 3 perovskites hold promise as robust materials for electrical interconnects and sensor applications in harsh environments. In this work, we investigated the high-temperature behavior of La 1–x Sr x CrO 3 , La 1–x Ca x CrO 3 and La 0.8 Sr 0.2 Cr 1-x Mn x O3 (0.1 ≤ x ≤ 0.4) through a combination of computational modeling and physical characterization up to 1500 °C. Crystalline structural properties were determined and compared with ab-initio calculations, which demonstrated excellent agreement with experimental findings. High-temperature electrical conductivity measurements were performed under different atmospheres. Calcium and strontium/manganese co-doped lanthanum chromites exhibited typical semiconductor exponential trends and conductivity showed proportional correlation with substitutions levels up to 30 %. The DFT modelling was completed up to 1500 °C, including low and high temperature chromite phases and oxygen vacancies insertion. Calculations were correlated with experimental electrical properties. This work expands the understanding of doped lanthanum chromites and paves the way for the development of materials suitable for demanding high-temperature applications.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2025

Characterization of contaminants in the Lyman-alpha forest auto-correlation with DESI

Baryon Acoustic Oscillations can be measured with sub-percent precision above redshift two with the Lyman-α (Lyα) forest auto-correlation and its cross-correlation with quasar positions. This is one of the key goals of the Dark Energy Spectroscopic Instrument (DESI) which started its main survey in May 2021. We present in this paper a study of the contaminants to the Lyα forest which are mainly caused by correlated signals introduced by the spectroscopic data processing pipeline as well as astrophysical contaminants due to foreground absorption in the intergalactic medium. Notably, an excess signal caused by the sky background subtraction noise is present in the Lyα auto-correlation in the first line-of-sight separation bin. We use synthetic data to isolate this contribution, we also characterize the effect of spectro-photometric calibration noise, and propose a simple model to account for both effects in the analysis of the Lyα forest. We then measure the auto-correlation of the quasar flux transmission fraction of low redshift quasars, where there is no Lyα forest absorption but only its contaminants. We demonstrate that we can interpret the data with a two-component model: data processing noise and triply ionized Silicon and Carbon auto-correlations. This result can be used to improve the modeling of the Lyα auto-correlation function measured with DESI.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

A dynamic solvent chamber propagation estimation framework using RNN for warm solvent injection in heterogeneous reservoirs

Warm solvent injection (WSI), injecting low-temperature solvent into formations to reduce the viscosity of heavy oil, is a clean technology for heavy oil production through reducing greenhouse gas emissions and water usage. The success of WSI operation depends on the uniform development and propagation of solvent chambers in reservoirs. However, reservoir heterogeneity stemming from shale barriers plays a detrimental role in the conformance of solvent chamber development and oil production rate. In this work, we developed a novel recurrent neural network (RNN)-based framework with the capability of efficiently tracking and estimating the solvent chamber positions in heterogeneous reservoirs based on only production time-series data. The developed estimation model utilizes the “sequence-to-sequence" mapping methodology to correlate observed production time-series sequence and solvent chamber edge sequence via a long short-term memory (LSTM) algorithm. The trained RNN models exhibit high accuracy, evidenced by the predicted dynamic solvent chamber locations match the corresponding true locations from numerical simulation, with a high coefficient of determination (R 2 ) and a low mean squared error. Specifically, the achieved R 2 values exceed 0.98 on both the training and testing data. The developed RNN-based workflow was tested via several cases from both regularly- and irregularly-shaped shale barriers, and the results were promising. The predicted solvent chambers showed strong agreement with those obtained from numerical simulations. The major benefits of this workflow include reducing computational time and saving overall monitoring and tracking costs for conventional techniques. In conclusion, the present work would provide a good demonstration of the capability of practical integration of machine learning methods in solving engineering problems.

58 GEOSCIENCES↗

DOE OSTI2025

Characterization of Arsenic and Selenium in Coal Fly Ash to Improve Evaluations for Disposal and Reuse Potential (Final Technical Report)

Coal fly ash is a high volume waste material that is discarded in landfills and surface water impoundments across the U.S. and is also widely recycled for a variety of applications. The leaching of potential of contaminants of concern, such as arsenic (As) and selenium (Se), is often the driver of risk assessments for coal ash disposal and reuse. The extent of leachable As and Se depends on several factors related to environmental conditions and fly ash characteristics. Previous studies employed various methods to delineate the concentration, chemical form, and distribution of As and Se in fly ash materials. However, few studies have attempted to directly correlate these properties to mobilization parameters relevant to disposal and reuse. Instead, the coal residuals industries often rely upon standardized leaching protocols that can be laborious or involve hazardous chemicals. The goals of the project were to: 1) Develop and evaluate a characterization protocol that can be used to screen fly ash samples for leachability of As and Se; 2) Characterize As, Se, and associated constituents of fly ash particles at multiple length scales (nanometer to micrometer) to determine if elemental associations differ as a function of the resolution of characterization; and 3) Establish a predictive model for the chemical composition of coal ash produced annually at major U.S. coal fired power facilities on 50-year national coal supply records. For the first objective, we performed leaching experiments with 52 fly ash samples collected from 15 different U.S. power plants and representing coal feedstocks from the three major domestic coal regions. For this work, we assessed the mobilization potential of As and Se in fly ash based on standardized leaching protocols and performed multivariate and lasso regression analyses to explore correlations of leachable As and Se contents with characteristics such as major element contents, loss on ignition (LOI) and pH. The results of regression models indicated that major elements (Fe, Ca, Al) for a wide range of fly ashes can serve as predictor variables for the leaching potential of As, but not for Se. LOI and pH were not important predictive variables in the models. Both regression approaches resulted in relatively strong fits for leachable As (correlation coefficient R 2 = 0.78 for both models) compared to models for leachable Se (R 2 = 0.49). Overall, these results suggest that correlation models combined with on-site elemental analysis with portable analyzers may enable a screening method for leachable As in coal ash. For the second objective, we utilized nanoscale 2-D imaging (30-50 nm spot size) with the Hard X-ray Nanoprobe (HXN) in combination with microprobe X-ray capabilities (~5 µm resolution) to determine As and Se elemental associations in fly ash particles. Speciation of As and Se was also measured at the nano- to microscale with X-ray absorption spectroscopy. The enhanced resolution of HXN showed As and Se that were diffusely located around or comingled with Ca- and Fe-rich particles. The results also showed nanoparticles of Se attached to the surface of fly ash grains. Overall, a comparison of As and Se species across scales highlights the heterogeneity and complexity of chemical associations for these trace elements of concern in coal fly ash. For the final objective, we developed a predictive model for major element composition of coal ash in reserve at disposal sites of major U.S. coal fired power plants. This model was constructed from coal purchase records of 705 power stations from 1973-2022 and was trained on coal ash composition data showing that coal ash elemental composition is strongly associated with the source of feedstock coal. The model showed regional shifts in the major element contents of ash produced by power plants in the last 50 years, particularly for calcium and iron (expressed as %CaO and %Fe 2 O 3 ), as coal-fired power stations changed their source of coal over this time frame. Our approach enables an estimation of coal ash chemical composition that is stored in waste impoundments at individual power stations. Such information can help delineate the regional market potential for material applications that would utilize coal ash harvested from disposal sites across the U.S.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2026

Modeling and Experimental Demonstration of Flux Spreading in Light Trapping Planar-Cavity Solar Enclosed Particle Receivers

This study experimentally validates and numerically models the flux-spreading effect in a light-trapping planar-cavity solar receiver) for particle-based concentrating solar power systems. The receiver's shallow cavity with vertical planar walls redistributes concentrated solar flux, reducing peak intensity and achieving uniform heat flux. On-sun tests at National Renewable Energy Laboratory's High-Flux Solar Furnace under flux up to 1500 kW/m2 measured cavity wall temperatures, which were compared with Monte Carlo ray-tracing (SolTrace) and computational numerical simulations. Three angular absorptance models were evaluated: constant absorptance, a Pyromark-based directional model, and a Fresnel-based Cr2O3 model. The Fresnel-derived model showed the best agreement with experiments, achieving high correlation (PC > 0.85), structural similarity (SSIM > 0.98), and signal-to-noise ratios (PSNR > 40 dB), with temperature prediction errors of 1-11%. Results confirm that flux spreading mitigates local overheating and validate the integrated modeling approach, supporting the solar receiver scalability for high-efficiency, high-temperature concentrating solar power applications.

14 SOLAR ENERGY↗

DOE OSTI2025

Water Temperature, Prey Concentration and Salmonid Density Influence Daily Growth of Wild Juvenile Salmonids in Tributaries of the Upper Salmon River, Idaho ( USA )

ABSTRACT Theory, experiments and field studies indicate that the somatic growth rate of freshwater consumers is shaped by the individual, additive and multiplicative effects of multiple factors, including consumer size and condition, temperature, prey resources and biotic interactions. While our understanding of how these factors affect wild populations of freshwater consumers is improving, the topic remains poorly studied, especially with respect to mobile species. Here, we report on an 8‐year, seven‐stream ( n = 49 stream‐year combinations) observational study examining the individual and interactive effects of invertebrate prey concentration (F, mg/m 3 ), mean daily water temperature (T, °C) and juvenile Chinook salmon ( Oncorhynchus tshawytscha ) density (D, fish/100 m 2 ) on summer daily growth rates (%/d) of mobile, anadromous, juvenile Chinook salmon (age‐0+, n = 382) and sub‐yearling (age‐0+, n = 61) and yearling (age‐1+) steelhead trout ( O. mykiss , n = 70) rearing in cold (mean daily summer: 12.1°C, range: 4.2°C–16.7°C) mountain tributaries of the Salmon River basin in central Idaho (USA). AIC c model selection indicated that daily juvenile salmonid growth positively correlated with water temperature, prey biomass concentration, local juvenile Chinook density and the interaction between water temperature and food but with species and age‐specific differences. Water temperature was a covariate in all top‐ranked models, with daily growth (%/day) rate increasing (0.05%–0.23%/d) linearly with mean daily summer water temperature. In addition to a direct positive relationship with daily growth rate, there was evidence that prey concentration positively interacted with water temperature to accelerate daily growth (F × T). The positive relationship between juvenile salmonid daily growth rate and juvenile Chinook density is difficult to explain and could result from confounding factors. The individual success observed in these streams may contribute to population‐level benefits for the focal consumers, as prey‐rich, warm summers may result in larger individuals with higher energy reserves at the end of the summer/autumn growing season, contributing to improved overwinter survival. Our results, taken in combination with evidence from models, experiments and observational studies, have climate change implications. Current and predicted increases in water temperature will necessitate higher rates of prey consumption by aquatic ectothermic consumers to offset accelerated metabolic demands. Thus, to improve the resilience of mobile freshwater consumers in a warming climate, we suggest that natural resource managers not only consider physical and chemical habitat conditions but also biotic conditions, including the spatiotemporal quantity and quality of prey resources.

Kiffney, Peter. M. [Fish Ecology, Northwest Fisher↗

DOE OSTI2025

Root genetics in the field to understand drought adaptation and carbon sequestration (Final Scientific/Technical Report)

For all crop plants, roots play a critical role in growth. Roots anchor the plants, and are the primary site of nutrient and water uptake. Roots are also the main source of C to soil in the form of root tissues and exudates, and thus greatly influence SOM stocks. To perform these functions, primary roots extend into soil, producing a network of branching roots of characteristic form, known as its root system architecture (RSA). RSA varies among species, and among varieties within a species that are adapted to different environments. Root traits are major targets for the second green revolution because of their potential to improve crop productivity, increase drought tolerance and nutrient acquisition, and increase C capture of soil. Improving the quality of roots in maize will be particularly valuable, since this crop is planted on over 92 million acres annually in the US. The future sustainability of agricultural systems relies on their ability to enhance soil organic matter (SOM) storage and reduce GHG emissions, while maintaining or enhancing productivity. This program had two components, Sensors and Models. For the first component, we designed and built a high-throughput phenotyping platform for root pulling of maize plants. This eliminated the physical labor of manually pulling up plants and reduced the number of personnel required down to one. The standardized pulling mechanism allowed recording force curves during the pulling process, providing additional information. We validated that the maximum force for pulling the root system was well-correlated with the root system mass and provided root crowns for further RSA analysis. These root crowns identified significant correlations with 2D root area and root depth, along with 3D root volume, total root length and number of root tips. We then used this system for field-based studies in maize on the genetics of root system architecture and its relation to nitrogen-use efficiency (NUE), including using lines relevant to the Corteva breeding program. Varieties were also evaluated at Corteva sites in the cornbelt and Danforth farm in Missouri, to establish responses across sites. From these studies we have identified genetic loci associated with root traits and created mutant lines for these loci and correlations of root traits with NUE. For the Models component, we worked to incorporate root and soil characteristics into the MEMS 2.0 soil and ecosystem biogeochemical model. Existing soil C models, such as Century, are unable to represent specific root trait interactions with the soil environment and therefore to accurately forecast the potential C sequestration benefits of root breeding under different climatic and soil type conditions. We have developed the MEMS 2.0 ecosystem biogeochemical model to improve quantification of farm-scale soil carbon and greenhouse gas emissions. The new knowledge and large datasets produced by this project will be used to develop and drive an innovative model capable of forecasting the impacts on soil C stocks and nutrient dynamics. An innovation was to use the empirical data from the field studies (in 1, above) to model genetic variation in nitrogen use efficiencies and soil C input. Our work demonstrated that maize root-derived C rapidly replaces existing soil C and after 3 years of continuous maize, up to 20% of soil organic C in the topsoil (0-15cm) and 3% in the subsoil (15-30cm) was contributed by maize. However, this contribution did not entirely represent a net increase. Root C contribution to soil was affected by maize genetics. We have analyzed soils derived from the CSU field trials for C and N stocks, in the different soil physical fractions represented by the MEMS model, using both physical fractionation with elemental analyses, and Fourier transformed infrared spectroscopy. Data will be used to link crop nitrogen use efficiencies with soil C sequestration and provide data to bridge the field trials with the model development, for verification of model predictions. The project had a number of successful outcomes: we have used the new phenotyping platform to identify new genetic loci that can enhance root phenotypes; we have partnered with multiple maize seed companies phenotype varieties in their breeding programs; we have developed the MEMS model that can help inform industry on the potential for carbon sequestration in the agricultural sector, and which is now available at the CSU Soil Carbon Solutions Center for use.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Advanced Computational Modeling of High-Level Waste Vitrification at the Hanford Site

The U.S. Department of Energy (DOE) has selected vitrification for stabilizing legacy tank waste at the Hanford site, where radioactive waste from plutonium production was historically stored in underground tanks. This waste will be separated into low-activity waste (LAW) and high-level waste (HLW) fractions and processed at the Waste Treatment and Immobilization Plant (WTP). At WTP, glass melters are used for the vitrification of radioactive tank waste, transforming it into a stable borosilicate glass form for safe long-term storage. The melter vessel is constructed from highly durable and heat-resistant materials, where the vitrification process occurs. The main regions that are modeled are the melt pool, plenum, cold cap, riser/discharge chamber, and surrounding structure with insulation layers. Forced convection induced by air bubblers at the base of the melter ensure uniform temperature distribution and provide heat to the cold cap layer. The cold cap is a region of reacting batch feed that floats on top of the molten glass and is where the batch-to-glass reactions occur. Joule heating provided by electrodes mounted along the vertical walls of the melter and immersed directly in the glass, generates the necessary heat for the net endothermic conversion processes that occur in the cold cap. The high temperatures, radioactivity, and opaque nature of the glass prevent direct observation inside the melters. Therefore, computational models are essential for providing insight into factors that affect melter throughput. Thermocouples in the plenum provide operators with plenum temperature measurements. Operational adjustments include bubbling rate, voltage supplied to the electrodes, feed adjustments, and glass removal rate. Different computational fluid dynamics (CFD) models have been developed, each serving a specific purpose. There are CFD models of different scale melters, as well as models that capture the two-phase flow interfaces of rising bubbles in the molten glass or models with a simplified molten glass region so that the surrounding structure and plenum can be feasibly incorporated. Pilot-scale melter models have been developed to serve as validation of the methods employed in the simulation of the full-scale WTP melters. Models incorporating resolved bubbling are used to develop momentum source terms to implement into a single phase, multi-region, steady-state flow model that is being validated by measured process parameters such as glass production rate, voltage, input power, plenum temperatures, etc. The resolved bubbling model uses the multiphase volume of fluid approach to model the system with a high-resolution interface capturing scheme to maintain sharp interfaces between the molten glass and the air phase. The suite of CFD models is continually being improved to incorporate more realistic physics and achieve faster turnaround time. For example, an incremental controller is implemented to automatically adjust electrode voltage within the simulation to a molten glass set point temperature of 1150°C. Newer models feature improved meshes to ensure conformal meshes between regions and eliminate unnecessary mesh refinement in areas that are not of interest (such as boundary layers in offgas ports). Instead of explicitly modeling the structural, refractory, and insulation layers of the melter, a thermal resistance approach is used with published correlations used for boundary conditions. The development of robust and efficient CFD models will be instrumental in enabling the WTP to successfully fulfill its mission of safely stabilizing legacy nuclear waste.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗

DOE OSTI2025

Multipoint Correlations in Poisson Media

Multipoint correlations in randomly heterogeneous composite media are determined by the probability that a set of points belong to specific phases. They determine a wide range of macroscopic transport properties such as conductivity, dielectric constant, diffusion coefficient, and transmittance. The Poisson model—a random tesselation of space by hyperplanes—provides realistic descriptions of heterogeneous media in, e.g., radiation-transport applications; yet, until now, it has lacked closed-form expressions for its multipoint correlations. We resolve this problem by presenting an exact solution for the multipoint correlations in the Poisson model. By comparing it to Monte Carlo simulations of four-point correlations in three dimensions, we demonstrate the accuracy of our solution. In conclusion, we visualize the multipoint correlations and discuss their features.

Amorphous materials↗

DOE OSTI2024

Sparsity-Independent Lyapunov Exponent in the Sachdev-Ye-Kitaev Model

The saturation of a recently proposed universal bound on the Lyapunov exponent has been conjectured to signal the existence of a gravity dual. This saturation occurs in the low-temperature limit of the dense Sachdev-Ye-Kitaev (SYK) model, N Majorana fermions with q body ( q > 2 ) infinite-range interactions. We calculate certain out-of-time-order correlators (OTOCs) for N ≤ 64 fermions for a highly sparse SYK model and find no significant dependence of the Lyapunov exponent on sparsity up to near the percolation limit where the Hamiltonian breaks up into blocks. This provides strong support to the saturation of the Lyapunov exponent in the low-temperature limit of the sparse SYK. A key ingredient to reaching N = 64 is the development of a novel quantum spin model simulation library that implements highly optimized matrix-free Krylov subspace methods on graphical processing units. This leads to a significantly lower simulation time as well as vastly reduced memory usage over previous approaches, while using modest computational resources. Strong sparsity-driven statistical fluctuations require both the use of a much larger number of disorder realizations with respect to the dense limit and a careful finite size scaling analysis. The saturation of the bound in the sparse SYK points to the existence of a gravity analog that would enlarge substantially the number of field theories with this feature. Published by the American Physical Society 2024

Physics↗