SEARCH · Engineering Papers

Results for “model data volumes”

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 1,009 records · Page 56

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Big Data For Operation and Maintenance Cost Reduction

The purpose of this research is to develop a first-of-a-kind framework for integrating Big Data capability into the daily activities of our current fleet of nuclear power plants. Big Data is traditionally defined as data sets with high volume, velocity, and heterogeneity, and the existing Big Data analytics capabilities are now widely popular in fields such as finance, weather, e-commerce, healthcare and sports. In the nuclear industry, while the volume and velocity of data may present computational challenges for existing analytics capabilities, data heterogeneity are seen to present the major challenge. This research project mainly focuses on incorporating the wide range of data heterogeneities in nuclear power plants into an integrated Big Data Analytics capability. The primary end-product of this project is a Big Data framework that is capable of dealing with the large volume and heterogeneity of the data found in nuclear power plants to extract timely and valuable information on equipment performance. The framework can generate system insights that are actionable relations between measurable impacts and the corresponding maintenance action plans and enable optimization of plant operation and maintenance based on the extracted information. The developed framework is capable of handling heterogeneous data including both image data and time-series sensor data. Specifically, this developed framework includes the following components. The first component is an overarching maintenance ontology which includes system insights required by maintenance optimization. The maintenance ontology interacts with other components in the developed framework. The second component handles Piping & Instrumentation Diagram (P&ID) data. It can be used to extract system components and their relations automatically from the P&IDs. This extracted information is stored in the first component, i.e., maintenance ontology, and is also used as input to the third component, i.e., a tool for generating the fault tree for the corresponding system. The generated fault tree in turn is stored in the ontology for assessing risk that is used as a criterion in maintenance policy optimization. The fourth component is a tool for inferring the parameters in the Markov degradation model for a nuclear system. It uses basic information from the ontology. The fifth component is a tool for assessing the degradation level using sensor measurement data, for example, pressure, flowrate. This tool can be used for determining corrective maintenance actions. The results obtained from components four and five are returned to the ontology. The sixth component of the framework is a tool for optimizing the maintenance policy for a nuclear system of interest. It takes certain basic information from the ontology, e.g., costs of maintenance actions and system failures, as input, and returns the optimal maintenance policy to the ontology. This tool can be used for determining predictive maintenance actions. A set of experiments have also been conducted to verify the algorithms developed in this project for nuclear system degradation monitoring. The experiments are based on four solenoid valves, similar to the ones used in nuclear power plants. The analyses based on the experimental data using two algorithms, i.e., the Randomized Window Decomposition (RWD) algorithm and the particle filtering algorithm, and the results are introduced in the report. The Big Data framework developed in this project can be used as a support tool in daily activities of plant operation and maintenance and will reduce current costs while maintaining or improving safety levels. Overall, the project will not only benefit existing reactors, however it will open new frontiers to realize the long overdue value of Big Data Analytics in the nuclear sphere.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2021

Quasi-one-dimensional uniform spin- 1 2 Heisenberg antiferromagnet KNaCuP 2 O 7 probed by P 31 and Na 23 NMR

Here we present the structural and magnetic properties of KNaCuP 2 O 7 investigated via x-ray diffraction, magnetization, specific heat, and P 31 and Na 23 NMR measurements and complementary electronic structure calculations. The temperature-dependent magnetic susceptibility and P 31 NMR shift could be modeled very well by the uniform spin- 1 2 Heisenberg antiferromagnetic chain model with a nearest-neighbor interaction J / k B ≃ 58.7 K . The corresponding mapping using first-principles electronic structure calculations leads to J DFT / k B ≃ 59 K with negligibly small interchain couplings, further confirming that the system is indeed a one-dimensional uniform spin- 1 2 Heisenberg antiferromagnet. The diverging trend of NMR spin-lattice relaxation rates ( 1 / T 1 31 and 1 / T 1 23 ) implies the onset of a magnetic long-range ordering at around T N ≃ 1 K . From the value of T N , the average interchain coupling is estimated to be J ' / k B ≃ 0.28 K . Moreover, the NMR spin-lattice relaxation rates show the dominant contributions from uniform ( q = 0 ) and staggered ( q = ± π / a ) spin fluctuations in the high- and low-temperature regimes, respectively, mimicking one-dimensionality of the spin lattice. We have also demonstrated that 1 / T 1 31 in high temperatures varies linearly with 1 / H , reflecting the effect of spin diffusion on the dynamic susceptibility. The temperature-dependent unit cell volume could be described well using the Debye approximation with a Debye temperature of Θ D ≃ 294 K , consistent with the heat capacity data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

A scalable framework for efficient coupling of thermal and microstructural simulations in additive manufacturing

Predicting microstructure evolution in metal additive manufacturing (AM) is important for process optimization, but spatiotemporal scale disparities between thermal transport and microstructure evolution create significant challenges for efficient data transfer between simulation codes. To address this, we present Stork, a scalable framework for coupling thermal and microstructural simulations. Stork uses a sparse data representation to identify and store active solidification sub-volumes, enabling highly parallel quad-linear interpolation from coarse thermal grids to fine microstructure grids without large intermediate storage. We demonstrate the framework by coupling the semi-analytic heat transfer code 3DThesis with the time-parallel cellular automata code Toucan. This approach achieves over two orders of magnitude reduction in data generation time and file size compared to prior workflows. Numerical studies show that quad-linear interpolation preserves grain morphology and crystallographic texture in laser powder bed fusion (LPBF) simulations for coarsening ratios up to 16. Overall, Stork provides a scalable pathway for high-throughput, component-scale AM simulations on modern high-performance computing systems.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Online peak-aware energy scheduling with untrusted advice

This paper studies the online energy scheduling problem in a hybrid model where the cost of energy is proportional to both the volume and peak usage, and where energy can be either locally generated or drawn from the grid. Inspired by recent advances in online algorithms with Machine Learned (ML) advice, we develop parameterized deterministic and randomized algorithms for this problem such that the level of reliance on the advice can be adjusted by a trust parameter. We then analyze the performance of the proposed algorithms using two performance metrics: robustness that measures the competitive ratio as a function of the trust parameter when the advice is inaccurate, and consistency for competitive ratio when the advice is accurate. Since the competitive ratio is analyzed in two different regimes, we further investigate the Pareto optimality of the proposed algorithms. Our results show that the proposed deterministic algorithm is Pareto-optimal, in the sense that no other online deterministic algorithms can dominate the robustness and consistency of our algorithm. Furthermore, we show that the proposed randomized algorithm dominates the Pareto-optimal deterministic algorithm. Our large-scale empirical evaluations using real traces of energy demand, energy prices, and renewable energy generations highlight that the proposed algorithms outperform worst-case optimized algorithms and fully data-driven algorithms.

Lee, Russell↗

DOE OSTI2020

Geothermal Deep Direct Use for Turbine Inlet Cooling in East Texas

The National Renewable Energy Laboratory (NREL), the Southern Methodist University Geothermal Laboratory (SMU), Eastman Chemical (Longview, TX), and TAS (Houston, TX) evaluated the feasibility of using geothermal heat to improve the performance of a natural-gas power plant in East Texas. The area of interest is the Eastman Chemical plant in Longview, Texas, which is on the northwestern margin of a geologic region known as the Sabine Uplift. The feasibility study focused on determining the potential for accessing a subsurface hot-water geothermal resource within a 10-km radius of the site to provide thermal energy for absorption chillers. Wells within a 20-km radius are included for broader geological comparison to determine the heat flow, temperature-at-depth, field porosity and permeability. The lithologies of most interest are the Lower Cretaceous Trinity Group and Upper Jurassic Cotton Valley Group. The deeper Cotton Valley formations are hotter (averaging 117 to 130°C), yet permeability and porosity are low. The shallower Trinity Group contains more variability in permeability and porosity and lower temperatures averaging about 98 to 117°C. The shallower formations are considered despite the lower temperature because of increased ability to produce larger volumes of water and extract enough heat before reinjection. The complete SMU analysis is available in the National Geothermal Data System (NGDS). Tapping such deep geothermal sources for direct heating (as opposed to power generation) is known as geothermal deep direct use (DDU). Geothermal DDU has potential across a wide swath of the United States but is underutilized due to challenging project economics associated with developing a deep geothermal resource for what are typically small-scale, variable-demand projects. This project examines the feasibility of geothermal energy integration in a natural-gas combined cycle power station in East Texas. The DDU resource is tapped to drive absorption chillers (24/7) for production of chilled water at 5-10°C (41-50°F). This chilled water is stored until needed, which allows for continuous operation with a relatively small-capacity geothermal/absorption chiller system. When conditions are favorable, the chilled water is dispatched to cool the air entering the compressor stage of a gas combustion turbine. This process, known as turbine inlet cooling (TIC), boosts power production during periods of high temperature and high-power demand. Such systems can enhance grid reliability and reduce the cost for peak-demand power. A simulation model of the power plant was developed in IPSEpro software and validated against operational data from the plant. This model allowed the team to estimate the additional power that could be produced by applying TIC under different operating and ambient conditions. Absorption chiller performance was estimated from vendor sources to determine the production rate of chilled water from the geothermal resource. Geothermal drilling and development costs were estimated using NREL's GEOPHIRES 2.0. The expected lower drilling costs in this region led to an estimated cost of geothermal heat of about $4/MMBtu (1.4 cents/kWh t ). The estimated cost for the absorption chillers and TIC hardware were obtained from literature sources and project partners. Hourly data were obtained for weather, natural gas and electricity prices, and plant operating state for 2017, which served as a representative year. NREL estimated the capital cost, operating cost, and additional electricity production and revenue for different combinations of geothermal capacity, chiller capacity, and water storage-tank size. The analysis drove toward smaller geothermal and chiller systems to reduce equipment cost. A relatively low-cost water storage tank accumulated the near-continuous chilled water output for later use when TIC was most valued.

15 GEOTHERMAL ENERGY↗

DOE OSTI2020

Performing k eff Validation of As-Loaded Criticality Safety Calculations Using UNF-ST&DARDS: Sensitivity Calculations

The general method for performing validation of as loaded criticality safety calculations using UNF ST&DARDS is presented in a paper by Clarity, which includes a description of the UNF-ST&DARDS system. Proof-of-principle analyses were performed in the summer of 2019 for MPC-32 dual purpose canisters (DPCs) containing pressurized water reactor (PWR) fuel assemblies. Summaries of these results are presented in this and a companion paper for this conference. The current paper describes the TSUNAMI-3D calculations performed to generate sensitivity data, and the companion paper discusses the selection of critical experiments applicable for validation of the 11 MPC-32 DPCs considered. The generation of sensitivity data for as-loaded spent nuclear fuel (SNF) DPCs is a challenge given the detailed model of the fuel compositions generated by UNF ST&DARDS. Each fuel assembly is modeled with its own irradiation history in 18 axial nodes, unless the fuel assembly is damaged and thus considered as fresh by design basis. This results in a set of 576 fuel compositions, each of which must be processed separately in a multigroup (MG) calculation. Therefore, a continuous-energy (CE) TSUNAMI-3D method was chosen to alleviate this challenge. Two CE TSUNAMI-3D methods are available in SCALE 6.2.3: the iterated fission probability (IFP) and contribution-linked eigenvalue sensitivity/uncertainty estimation via track-length importance characterization (CLUTCH). Since the IFP method is not feasible because of memory requirements associated with its implementation in SCALE, the CLUTCH method was selected for these calculations. CLUTCH has been implemented in SCALE in parallel, allowing long calculations to be performed in reasonable timeframes. The two primary user inputs necessary for CLUTCH calculations are the F*(r) mesh and the number of latent generations used in determining the F*(r) function. This F*(r) function is used as the importance function for fission chains originating in a given volume element (voxel), and it is calculated using the IFP method in the skipped generations. A large number of skipped generations is thus required to ensure accurate calculation of this importance function. In these calculations, 500 generations were used to calculate the F*(r) function. For more information regarding the calculation of F*(r), see Jones [4]. The remainder of this paper is focused on the selection of the F*(r) mesh and the number of latent generations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

A robust spectral element implementation of the $k - τ$ RANS model in Nek5000/NekRS

The $k - ω$ Reynolds Averaged Navier Stokes (RANS) model is one of the industry standard approaches for modeling of turbulent flows. It performs better than the $k - ϵ$ model for low Reynolds number flows and is also more suitable for boundary layers with adverse pressure gradients. Major drawback of the model, however, is that the asymptotic value of $ω$ at the walls is singular, necessitating the use of a contrived “sufficiently” large value for $ω$ as the boundary condition for its transport equation. Here, this invariably leads to the solution being sensitive to near wall grid spacing. While an acceptable solution for low order (finite volume) methods, the excessive near wall gradients lead to persistent numerical stability issues in high order codes. To alleviate the problem, specifically in the context of the high order spectral element code Nek5000, a regularized $k - ω$ approach was formulated in our prior work (Tomboulides et al., 2018). The formulation, however, relies on the use of wall distance and its gradients for modeling the closure terms and can pose problems for simulations in complex geometries. This work presents a novel implementation of the $k - τ$ RANS model in Nek5000, where $τ = 1/ω$, eliminating the need for regularization, owing to the asymptotically bounded behavior of the source terms in the $τ$ transport equation, and also eliminating dependence on wall distance. Robustness and stability of the $k - τ$ model is ensured through implicit treatment of the source terms and their careful numerical implementation and demonstrated through several cases aimed at verification and validation. Studies include both canonical and engineering relevant problems, viz., turbulent channel flow, pipe flow, backward facing step, flow over NACA 0012 airfoil and flow in a T-junction. Results from the $k - τ$ model are shown to be consistent with regularized $k - ω$ model and also with the $k - ω$ SST model in OpenFOAM (for select studies). Comparison with experimental data is also shown, where available, to bolster validation efforts for the $k - τ$ model implementation through prediction of key turbulent quantities of interest.

Nek5000↗

DOE OSTI2024

CarbonSAFE Phase II: Optimizing Alabama’s CO 2 Storage in Shelby County, Alabama (Project OASIS), Milestone M3 - Site Specific Drilling Report

The Phase II Storage Complex Feasibility project, entitled "Optimizing Alabama's CO 2 Storage in Shelby County, Alabama (Project OASIS)," is a part of the DOE/NETL's CarbonSAFE initiative. The Project is managed by the Southern States Energy Board (SSEB), and includes participation from Advanced Resources International, Inc. (ARI), Crescent Resource Innovation, Southern Company, Alabama A&M University, Auburn University, and Oklahoma State University. Project OASIS is working to establish the foundation for a commercial-scale geologic storage complex for CO 2 captured from Plant Gaston (home of the National Carbon Capture Center) and surrounding industrial sources of CO 2 located in Shelby County, Alabama. The Project objectives are: • Demonstrate that the subsurface saline formations at the storage complex can store commercial volumes of CO 2 safely and permanently. • Develop a comprehensive Community Benefits Plan. • Develop the infrastructure framework for a CO 2 storage hub. • Develop a rigorous risk registry and to conduct a comprehensive risk assessment. • Develop a monitoring plan. • Develop a comprehensive site characterization plan to support an Underground Injection Control Class VI Permit in a future Phase III program. • Evaluate the commercial viability of the project. Project OASIS is about 30 miles southeast of Birmingham, Alabama within a geologic province called the Valley and Ridge (Figure 1.1). The Valley and Ridge Province comprises a sequence of Paleozoic carbonate and clastic rocks that underwent structural deformation during the Alleghanian Orogeny. Storage prospects occur in relatively flat lying structural panels located between thrust faults. Available geologic studies related to hydrocarbon exploration suggest that Cambro-Ordovician carbonates and Cambrian clastic units offer multiple potential storage intervals, and that regional confining systems are present, such as the tectonically thickened Floyd-Parkwood Shale. The surface property is owned by a timber and land stewardship company, The Westervelt Company, Inc., who worked with the Project Team to select and prepare adequate sites for geologic assessment. The purpose of drilling the Westover Stratigraphic Test Well #2 was to collect geologic data to model the feasibility of commercial scale CO 2 injection and storage. This includes geological and geophysical evaluations, reservoir engineering analyses, and risk assessments.

20 FOSSIL-FUELED POWER PLANTS↗

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 OSTI2021

Risk-based area of review estimation in overpressured reservoirs to support injection well storage facility permit requirements for CO 2 storage projects

This paper by the Energy & Environmental Research Center presents a workflow and modeling approach for delineating a risk-based area of review (AOR) to support a U.S. Environmental Protection Agency (EPA) Class VI permit for a carbon dioxide (CO 2 ) storage project. The approach combines semianalytical solutions for estimating formation fluid leakage through a hypothetical leaky wellbore with the results of numerical reservoir simulations to define the AOR. The modeling utilizes 1) semianalytical solutions from the peer-reviewed literature for formation fluid leakage through abandoned wellbores by Raven (1990) and Avci (1994), 2) a FORTRAN model compiled and described in Cihan et al. (2011, 2012) called ASLMA (Analytical Solution for Leakage in Multilayered Aquifers), and 3) a computational framework for estimating a risk-based AOR first proposed by Oldenburg et al. (2014, 2016). Therefore, the approach builds upon well-established research and underlying hydrogeological principles that have been upheld for nearly three decades. Moreover, the ASLMA model has been broadly applied to an array of storage projects. The work presented herein extends these earlier works using a custom wrapper written in the software environment, R (R Core Team, 2020), which was developed to perform multiple runs of the ASLMA model using given ranges for one or more input parameters. In addition, the current work simulates the pressure buildup within the storage reservoir in response to CO 2 injection using a compositional simulator to better accommodate the temporospatial evolution of pressure buildup within the storage reservoir that is more accurately modeled using a heterogeneous geologic model and a compositional simulator that accounts for the multiphase interactions. The workflow is demonstrated using a case study for a 180,000-metric-ton-per-year storage project located in the PCOR (Plains CO 2 Reduction) Partnership region. For the storage project evaluated here, under the scenario where the leaky wellbore is open to a saline aquifer (thief zone) between the overlying seal (cap rock) and the underground sources of drinking water (USDW), the risk-based AOR essentially collapses to the areal extent of the CO 2 plume in the storage reservoir because the pressure buildup in the storage reservoir beyond the CO 2 plume is insufficient to drive formation fluids up a hypothetical leaky wellbore into the USDW. However, even under the conservative assumption that the leaky wellbore is not open to a thief zone, beyond the areal extent of the CO 2 plume, the incremental leakage is less than 400 m 3 over 20 years, which represents ~0.0001% or less of the total volume of water contained within the USDW rock volume. As discussed in the text, the threshold criterion for defining the risk-based AOR is site-specific and should be informed by the results of the sensitivity analysis and available site characterization data. The approach outlined in this paper is designed to be protective of USDWs and, therefore, comply with the Safe Drinking Water Act requirements and provisions for the U.S. EPA Class VI Underground Injection Control (UIC) Program (Class VI Rule) and North Dakota Administrative Code Chapter 43-05-01.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Shifting groundwater fluxes in bedrock fractures: Evidence from stream water radon and water isotopes

Geologic features (e.g., fractures and alluvial fans) can play an important role in the locations and volumes of groundwater discharge and degree of groundwater-surface water (GW-SW) interactions. However, the role of these features in controlling GW-SW dynamics and streamflow generation processes are not well constrained. GW-SW interactions and streamflow generation processes are further complicated by variability in precipitation inputs from summer and fall monsoon rains, as well as declines in snowpack and changing melt dynamics driven by warming temperatures. Using high spatial and temporal resolution radon and water stable isotope sampling and a 1D groundwater flux model, we evaluated how groundwater contributions and GW-SW interactions varied along a stream reach impacted by fractures (fractured-zone) and below the fractured hillslope (non-fractured zone) in Coal Creek, a Colorado River headwater stream affected by summer monsoons. During early summer, groundwater contributions from the fractured zone dominated, but declined throughout the summer. Groundwater contributions from the non-fractured zone were constant throughout the summer and became proportionally more important later in the summer. We hypothesize that groundwater in the non-fractured zone is dominantly sourced from a high-storage alluvial fan at the base of a tributary that is connected to Coal Creek throughout the summer and provides consistent groundwater influx. Water isotope data revealed that Coal Creek responds quickly to incoming precipitation early in the summer, and summer precipitation becomes more important for streamflow generation later in the summer. We quantified the change in catchment dynamic storage and found it negatively related to stream water isotope values, and positively related to modeled groundwater discharge and the ratio of fractured zone to non-fractured zone groundwater. Importantly, we interpret these relationships as declining hydrologic connectivity throughout the summer leading to late summer streamflow supported predominantly by shallow flow paths, with variable response to drying from geologic features based on their storage. As groundwater becomes more important for sustaining summer flows, quantifying local geologic controls on groundwater inputs and their response to variable moisture conditions may become critical for accurate predictions of streamflow.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Thermoelastic properties of tungsten at simultaneous high pressure and temperature

In this work, the compressional ( P ) and shear wave velocities ( S ) and unit cell volumes (densities) of polycrystalline tungsten (W) have been measured simultaneously up to 10.5 GPa and 1073 K using ultrasonic interferometry in conjunction with x-ray diffraction and x-radiography techniques. Thermoelastic properties of W were derived using different methods. We obtained the isothermal bulk modulus K T 0 = 310.3(1.5) GPa, its pressure derivative K ' T 0 = 4.4(3), its temperature derivative at constant pressure ( ∂ K T / ∂ T ) P = – 0.0138 ( 1 ) GPa K – 1 and at constant volume ( ∂ K T / ∂ T ) V = – 0.0050 GPa K – 1 , the thermal expansion α ( 0, T ) = 1.02(27) × 10 –5 + 7.39(3.2) × 10 –9 T (K –1 ), as well as the pressure derivative of thermal expansion ( ∂ α / ∂ P ) T = – 1.44 ( 1 ) × 10 – 7 K – 1 GP a – 1 based on the high-temperature Birch–Murnaghan equation of state (EOS), the Vinet EOS, and thermal pressure approach. Finite strain analysis allowed us to derive the elastic properties and their pressure/temperature derivatives independent of the choice of pressure scale. A least-squares fitting yielded K S 0 = 314.5(2.5) GPa, K S 0 ' = 4.45(9), (∂ K S /∂ T ) P = – 0.0076(6) GPa K –1 , G 0 = 162.4(9) GPa, G 0 ' = 1.8(1), (∂ G /∂ T ) P = – 0.0175(9) GPa K –1 , and α 298 K = 1.23 × 10 – 5 K – 1 . Fitting current data to the Mie–Grüneisen–Debye EOS with derived θ 0 = 383.4 K yielded γ 0 = 1.81 ( 6 ) and q = 0.3 . The thermoelastic parameters obtained from various approaches are consistent with one another and comparable with previous results within uncertainties. Our current study provides a complete and self-consistent dataset for the thermoelastic properties of tungsten at high P – T conditions, which is important to improve the theoretical modeling of these materials under dynamic conditions.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Optimal gas cooling with minimal aero-opticaldistortion for next generation high-average-power lasers (Progress report, Year 2)

Next-generation lasers that combine high-peak-power ultrashort pulses with high repetition rates have the potential to efficiently deliver the high-average-power outputs required for inertial fusion energy generation, laser-based particle acceleration, and to generate compact secondary x- and 𝛾-ray radiation sources, among others. The combination of long gain length with high surface-to-volume ratio to remove residual heat makes gas-cooled multi-slab amplifier heads an inherently aperture- and average power scalable technology that enables this new class of solid-state lasers, and their thermal management becomes the most pressing technical challenge. As the index-of refraction depends on the gas density, heat removal from the slab surfaces unavoidably entails scattering of propagated light from regions of high temperature fluctuations. The requirement to maximize convective heat extraction from the gain medium under the constraint of minimum optical aberration of the beam renders amplifier slab cooling a technological challenge at the intersection of aero-optics, fluid-mechanics and heat transfer. The goal of this research is to leverage computational fluid dynamics and aero optics simulations of different levels of fidelity to devise a numerical modeling and optimization framework for gas-cooled multi-slab amplifiers. The project will be conducted in close collaboration with the Advanced Photon Technologies Group at Lawrence Livermore National Laboratory, where a gas-cooled test stand for laser materials and thermal management structures will be developed and installed. The final goal of the collaborative effort is to manufacture and benchmark a physical prototype of an numerically optimized amplifier head design. Measurement data on solid-state and thermo-optical distortions will be used for the validation of the numerical models. The key technical challenge in this multi-objective design problem is that of balancing optimal heat removal with minimal optical distortion. To facilitate computationally efficient design optimization and, at the same time, deepen our basic understanding of the interaction of light with turbulence, a complementary study that combines high and low-fidelity flow simulations and aero-optical models is proposed. Along the lines of these two objectives and solution strategies, the project is organized into two thrusts.

42 ENGINEERING↗

DOE OSTI2020

Compression behavior of dense H 2 -He mixtures up to 160 GPa

We have studied the compression behavior of H 2 -He mixtures in comparison with pure H 2 and He using powder synchrotron x-ray diffraction and present the pressure-volume (PV) compression data of H 2 -He mixtures to 160 GPa. The results indicate that both H 2 and He in H 2 -He mixtures remain in hcp to the maximum pressure studied, yet develop a substantial level of lattice distortion in the (100) plane, most profound in He-rich solids and below 66 GPa. The measured PV data also indicate softening of He (or H 2 )-rich lattice upon increasing the level of the guest H 2 (or He) concentration. We suggest that the observed softening and lattice distortion are due to a substitutional incorporation of H2 (guest) molecules into the basal plane of hcp-He (host) lattice and, thereby, reflect the miscibility between H 2 and He in H 2 -He mixtures. Interestingly, solid He exhibits a lesser degree of preferred orientation in H2-He mixtures than in pure He, likely due to the presence of solid H 2 disturbing the crystalline ordering of He-rich solids. Finally, the present PV compression data of H 2 -rich and He-rich solids to 160 GPa deviate from those of pure H 2 and pure He above ~70 and 45 GPa respectively, providing new constraints for development of the EOS for H 2 -He mixtures for planetary models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Enhancing high-fidelity nonlinear solver with reduced order model

Abstract We propose the use of reduced order modeling (ROM) to reduce the computational cost and improve the convergence rate of nonlinear solvers of full order models (FOM) for solving partial differential equations. In this study, a novel ROM-assisted approach is developed to improve the computational efficiency of FOM nonlinear solvers by using ROM’s prediction as an initial guess. We hypothesize that the nonlinear solver will take fewer steps to the converged solutions with an initial guess that is closer to the real solutions. To evaluate our approach, four physical problems with varying degrees of nonlinearity in flow and mechanics have been tested: Richards’ equation of water flow in heterogeneous porous media, a contact problem in a hyperelastic material, two-phase flow in layered porous media, and fracture propagation in a homogeneous material. Overall, our approach maintains the FOM’s accuracy while speeding up nonlinear solver by 18–73% (through suitable ROM-assisted FOMs). More importantly, the proximity of ROM’s prediction to the solution space leads to the improved convergence of FOMs that would have otherwise diverged with default initial guesses. We demonstrate that the ROM’s accuracy can impact the computational efficiency with more accurate ROM solutions, resulting in a better cost reduction. We also illustrate that this approach could be used in many FOM discretizations (e.g., finite volume, finite element, or a combination of those). Since our ROMs are data-driven and non-intrusive, the proposed procedure can easily lend itself to any nonlinear physics-based problem.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Dynamic mode decomposition with core sketch

With the increase in collected data volumes, either from experimental measurements or high fidelity simulations, there is an ever-growing need to develop computationally efficient tools to process, analyze, and interpret these datasets. Modal analysis techniques have gained great interest due to their ability to identify patterns in the data and extract valuable information about the system being considered. Dynamic mode decomposition (DMD) relies on elements of the Koopman approximation theory to compute a set of modes, each associated with a fixed oscillation frequency and a decay/growth rate. Extracting these details from large datasets can be computationally expensive due to the need to implement singular value decomposition of the input data matrix. Sketching algorithms have become popular in numerical linear algebra where statistical theoretic approaches are utilized to reduce the cost of major operations. A sketch of a matrix is another matrix, which is significantly smaller, but still sufficiently approximates the original system. We put forth an efficient DMD framework, SketchyDMD, based on a core sketching algorithm that captures information about the range and corange (their mutual relationship) of input data. The proposed sketching-based framework can accelerate various portions of the DMD routines, compared to classical methods that operate directly on the raw input data. We conduct numerical experiments using the spherical shallow water equations as a prototypical model in the context of geophysical flows. In conclusion, we show that the proposed SketchyDMD is superior to existing randomized DMD methods that are based on capturing only the range of the input data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Lowering the barrier to access information-rich transient kinetic data for machine learning methods

Transient kinetic data contain a wealth of information about intrinsic features of a catalyst as well as the reaction mechanism. Currently, high volume transient data is underutilized, and data science methods could both increase the value of information that can be extracted from this data, integrate experimental with theoretical data sources, and accelerate the pace of catalyst technology advancement. Transient kinetic characterizations with simple probe molecules exhibiting reversible adsorption, irreversible adsorption and bulk-surface diffusion are presented as training components for similar experiments with more complex surface reactions. In conclusion, by increasing the availability and accessibility of transient kinetic data through details of its structure and acquisition, we aim to decrease the barrier for data scientists to apply machine learning methods to this valuable data source.

Catalysis↗

DOE OSTI2022

Early-time measure in eternal inflation

In a situation like eternal inflation, where our data is replicated at infinitely-many other space-time events, it is necessary to make a prior assumption about our location to extract predictions. The principle of mediocrity entails that we live at asymptotic late times, when the occupational probabilities of vacua has settled to a near-equilibrium distribution. In this paper we further develop the idea that we instead exist during the approach to equilibrium, much earlier than the exponentially-long mixing time. In this case we are most likely to reside in vacua that are easily accessed dynamically. Using first-passage statistics, we prove that vacua that maximize their space-time volume at early times have: 1. maximal ever-hitting probability; 2. minimal mean first-passage time; and 3. minimal decay rate. These requirements are succinctly captured by an early-time measure. The idea that we live at early times is a predictive guiding principle, with many phenomenological implications. First, our vacuum should lie deep in a funneled region, akin to folding energy landscapes of proteins. Second, optimal landscape regions are characterized by relatively short-lived vacua, with lifetime of order the de Sitter Page time. For our vacuum, this lifetime is ~ 10 130 years, which is consistent with the Standard Model estimate due to Higgs metastability. Third, the measure favors vacua with small, positive vacuum energy. This can address the cosmological constant problem, provided there are sufficiently many vacua in the entire ensemble of funnels. As a concrete example, we study the Bousso-Polchinski lattice of flux vacua, and find that the early-time measure favors lattices with the fewest number of flux dimensions. This favors compactifications with a large hierarchy between the lightest modulus and all other Kähler and complex structure moduli.

79 ASTRONOMY AND ASTROPHYSICS↗