Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large system simulation”

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 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

New Developments and Verification of the Fusion Blanket Simulation Capabilities in the MOOSE Framework

This poster contains new developments and updates in the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework related to fusion modeling and simulation, and forms the basis for an on-going effort in modeling whole fusion devices at Idaho National Laboratory using MOOSE. In this work, a slice of the midplane of a proposed ceramic breeder blanket design for the Fusion Nuclear Science Facility is simulated using a combination of neutronics, 1D system level thermal hydraulics for helium, 3D heat conduction, and tritium transport. Previous work focused on the engineering scale full assembly of a blanket sector, but the simulations here form the basis for software quality and verification efforts to increase confidence as the large scale simulations continue. This poster was presented at the Idaho National Laboratory Intern Poster Session 2023.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Estimation of the radiation budget during MOSAiC based on ground-based and satellite remote sensing observations

An accurate representation of the radiation budget is essential for investigating the impact of clouds on the climate system, especially in the Arctic, an environment highly sensitive to complex and rapid environmental changes. In this study, we analyse a unique dataset of observations from the central Arctic made during the MOSAiC (Multidisciplinary drifting Observatory for the Study of Arctic Climate) expedition in conjunction with state-of-the-art satellite products from CERES (Clouds and the Earth's Radiant Energy System) to investigate the radiative effect of clouds and radiative closure at the surface and the top of the atmosphere (TOA). We perform a series of radiative transfer simulations using derived cloud macro- and microphysical properties as inputs to the simulations for the entire MOSAiC period, comparing our results to collocated satellite products and ice-floe observations. The radiative closure biases were generally within the instrumental uncertainty, indicating that the simulations are sufficiently accurate to reproduce the radiation budget during MOSAiC. Comparisons of the simulated radiation budget relative to CERES show similar values in the terrestrial flux but relatively large differences in the solar flux, which are attributed to a lower surface albedo and a possible underestimation of atmospheric opacity by CERES. While the simulation results were consistent with the observations, more detailed analyses reveal an overestimation of simulated cloud opacity for cases involving geometrically thick ice clouds. In the annual mean, we found that, during the MOSAiC expedition, the presence of clouds leads to a loss of 5.2 W m −2 of the atmosphere–surface system to space, while the surface gains 25.0 W m −2 and the atmosphere is cooled by 30.2 W m −2 .

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Neutrino emission from binary neutron star mergers: characterising light curves and mean energies

Abstract Neutrinos are copiously emitted by neutron star mergers, due to the high temperatures reached by dense matter during the merger and its aftermath. Neutrinos influence the merger dynamics and shape the properties of the ejecta, including the resulting r -process nucleosynthesis and kilonova emission. In this work, we analyse neutrino emission from a large sample of binary neutron star merger simulations in Numerical Relativity, covering a broad range of initial masses, nuclear equation of state and viscosity treatments. We extract neutrino luminosities and mean energies, and compute quantities of interest such as the peak values, peak broadnesses, time averages and decrease time scales. We provide a systematic description of such quantities, including their dependence on the initial parameters of the system. We find that for equal-mass systems the total neutrino luminosity (several $$10^{53}{\hbox {erg}~{\hbox {s}}^{-1}}$$ 10 53 erg s - 1 ) decreases as the reduced tidal deformability increases, as a consequence of the less violent merger dynamics. Similarly, tidal disruption in asymmetric mergers leads to systematically smaller luminosities. Peak luminosities can be twice as large as the average ones. Electron antineutrino luminosities dominate (initially by a factor of 2-3) over electron neutrino ones, while electron neutrinos and heavy flavour neutrinos have similar luminosities. Mean energies are nearly constant in time and independent on the binary parameters. Their values reflect the different decoupling temperature inside the merger remnant. Despite present uncertainties in neutrino modelling, our results provide a broad and physically grounded characterisation of neutrino emission, and they can serve as a reference point to develop more sophisticated neutrino transport schemes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Quantum simulations in effective model spaces: Hamiltonian-learning variational quantum eigensolver using digital quantum computers and application to the Lipkin-Meshkov-Glick model

Quantum simulations offer the potential to predict the structure and dynamics of nuclear many-body systems that are beyond the capabilities of classical computing. Generally, preparing the ground state of strongly-interacting many-body systems relevant to nuclear physics is however inefficient, even using ideal quantum computers. In addition, currently available NISQ-era quantum devices possess modest numbers of qubits, limiting the size of quantum many-body systems that can be simulated. In this context, a reformulation of the quantum many-body problems using truncated model spaces and Hamiltonians is desirable to make them more amenable to near-term quantum computers. The importance of symmetries in low-energy theories, including effective field theories (EFTs), lattice quantum chromodynamics (QCD), and effective model spaces for nuclear systems, in particular their interplay with the reduction of active Hilbert spaces, is well known. Lesser known is the fact that the non-commutivity of some symmetries and truncations of the model space can be profitably combined with variational calculations to rearrange the entanglement into localized structures and enable more efficient simulations. Here, the goal of the present study is to explore and utilize the non-commutivity of symmetries and model-space truncations of quantum many-body systems important to nuclear physics, particularly in combination with variational algorithms for quantum simulations and effective Hamiltonian learning. We introduce an iterative hybrid classical-quantum algorithm, Hamiltonian learning variational quantum eigensolver (HL-VQE), that simultaneously optimizes an effective Hamiltonian, thereby rearranging entanglement into the effective model space, and the associated ground-state wavefunction. Quantum simulations, using classical computers and IBM's superconducting-qubit quantum computers, are performed to demonstrate the HL-VQE algorithm, in the context of the Lipkin-Meshkov-Glick (LMG) model of interacting fermions, where the Hamiltonian transformation corresponds to an orbital rotation. We use a mapping where the number of qubits scales with the $\log$ of the size of the effective model space, rather than the particle number. HL-VQE is found to provide an exponential improvement in LMG-model calculations of the ground-state energy and wavefunction, compared to naive truncations without Hamiltonian learning, throughout a significant fraction of the Hilbert space. In the context of EFT, this corresponds to counterterms scaling exponentially with the cut-off as opposed to power law. Implementations on IBM's QExperience quantum computers and simulators for 1- and 2-qubit effective model spaces are shown to provide accurate and precise results, reproducing classical predictions. For a range of parameters defining the LMG model, the HL-VQE algorithm is found to have better scaling of quantum resources requirements than previously explored algorithms. In particular, the HL-VQE scales efficiently over a large fraction of the model space, in contrast to VQE alone. This work constitutes a step in the development of entanglement-driven quantum algorithms for descriptions of nuclear many-body systems. This, in part, leverages the potential of noisy intermediate-scale quantum (NISQ) devices. The exponential scaling of counterterms observed in this study suggests the possibility of more general applicability to other non-perturbative EFTs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

An overview of the E3SM version 2 large ensemble and comparison to other E3SM and CESM large ensembles

Abstract. This work assesses a recently produced 21-member climate model large ensemble (LE) based on the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM) version 2 (E3SM2). The ensemble spans the historical era (1850 to 2014) and 21st century (2015 to 2100), using the SSP370 pathway, allowing for an evaluation of the model's forced response. A companion 500-year preindustrial control simulation is used to initialize the ensemble and estimate drift. Characteristics of the LE are documented and compared against other recently produced ensembles using the E3SM version 1 (E3SM1) and Community Earth System Model (CESM) versions 1 and 2. Simulation drift is found to be smaller, and model agreement with observations is higher in versions 2 of E3SM and CESM versus their version 1 counterparts. Shortcomings in E3SM2 include a lack of warming from the mid to late 20th century, likely due to excessive cooling influence of anthropogenic sulfate aerosols, an issue also evident in E3SM1. Associated impacts on the water cycle and energy budgets are also identified. Considerable model dependence in the response to both aerosols and greenhouse gases is documented and E3SM2's sensitivity to variable prescribed biomass burning emissions is demonstrated. Various E3SM2 and CESM2 model benchmarks are found to be on par with the highest-performing recent generation of climate models, establishing the E3SM2 LE as an important resource for estimating climate variability and responses, though with various caveats as discussed herein. As an illustration of the usefulness of LEs in estimating the potential influence of internal variability, the observed CERES-era trend in net top-of-atmosphere flux is compared to simulated trends and found to be much larger than the forced response in all LEs, with only a few members exhibiting trends as large as observed, thus motivating further study.

58 GEOSCIENCES↗
DOE OSTI2020

Impact of particle arrays on phase separation composition patterns

We examine the symmetry-breaking effect of fixed constellations of particles on the surface-directed spinodal decomposition of binary blends in the presence of particles whose surfaces have a preferential affinity for one of the components. Our phase-field simulations indicate that the phase separation morphology in the presence of particle arrays can be tuned to have a continuous, droplet, lamellar, or hybrid morphology depending on the interparticle spacing, blend composition, and time. In particular, when the interparticle spacing is large compared to the spinodal wavelength, a transient target pattern composed of alternate rings of preferred and non-preferred phases emerges at early times, tending to adopt the symmetry of the particle configuration. We reveal that such target patterns stabilize for certain characteristic length, time, and composition scales characteristic of the pure phase-separating mixture. To illustrate the general range of phenomena exhibited by mixture-particle systems, we simulate the effects of single-particle, multi-particle, and cluster–particle systems having multiple geometrical configurations of the particle characteristic of pattern substrates on phase separation. Our simulations show that tailoring the particle configuration, or substrate pattern configuration, a relative fluid-particle composition should allow the desirable control of the phase separation morphology as in block copolymer materials, but where the scales accessible to this approach of organizing phase-separated fluids usually are significantly larger. Limited experiments confirm the trends observed in our simulations, which should provide some guidance in engineering patterned blend and other mixtures of technological interest.

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

Uranium recovery from seawater using amidoxime-based braided polymers synthesized from acrylic fibers

Global demand for nuclear energy is expected to rise in the coming decades. To meet these growing needs, new uranium resources must be explored. One of the potential alternatives to traditional uranium mining is oceanic uranium. The capture and recovery of uranium from the ocean has been under investigation for some time, with many recent studies focused on amidoxime based adsorbents. These adsorbents while able to achieve high uranium recovery capacities are, nevertheless, expensive to produce and adsorb a significant amount of hard to remove vanadium. The purpose of this study is to evaluate the adsorption performance of amidoxime based polymer braids synthesized from acrylic fibers which are designed to significantly cut polymer synthesis and conditioning costs. Adsorption experiments were performed in a recirculating raceway flume system at environmental conditions, approximately 10.8 oC, with 40-micron prefiltered seawater over 28 days for small and large sized LCW polymer braids. For both braid sizes, the adsorption of vanadium was far lower on the LCW materials than on Oak Ridge National Laboratory’s amidoxime based AF1 and AI8 small braids, with the LCW small braids outperforming both materials with respect to uranium adsorption under similar conditions. Adsorption modeling was used to simulate the performance at higher temperatures based on 20oC experiments previously performed with similar materials. Simulation results indicated that the small LCW braids would have a less significant advantage with respect to uranium adsorption over both AF1 and AI8 at 20oC and adsorb slightly less uranium than AF1 at 31oC. Vanadium adsorption by LCW small braids was less than one third of vanadium adsorption by either AF1 or AI8 for all temperatures. This behavior indicates that the LCW materials are able to achieve comparable uranium adsorption to other amidoxime adsorbents previously developed, while being cheaper to produce and adsorbing significantly less vanadium.

Wiechert, Alexander I.↗
DOE OSTI2026

From Rules to Reasoning: A Survey of Large Language Model-Based Approaches to Scientific Hypothesis and Idea Generation

Scientific hypothesis generation represents a fundamental challenge in contemporary research due to exponentially expanding literature volumes and increasing disciplinary specialization. Large language models (LLMs) have emerged as transformative tools for automated scientific discovery, moving beyond traditional rule-based and literature-mining approaches. Four paradigmatic approaches define current LLM-driven hypothesis generation: direct prompting and fine-tuning methods, knowledge-enhanced frameworks integrating retrieval-augmented generation (RAG), multi-agent collaborative systems simulating research teams, and reasoning-focused approaches implementing cognitive architectures. Domain-specific applications demonstrate statistical equivalence to human expert performance in social psychology, experimental validation in biomedical research, and near-expert quality in astronomy. Evaluation methodologies encompass human expert assessment, LLM-as-judge frameworks, and comprehensive benchmarking systems. Technical challenges include hallucination management, knowledge integration limitations, and balancing novelty with feasibility. Future directions emphasize hybrid neural-symbolic architectures and sophisticated human-AI collaboration models for responsible scientific discovery acceleration.

AI-driven discovery↗
DOE OSTI2025

Risk assessment of wellbore leakage during underground hydrogen storage

The expansion of renewable energy sources would require large-scale energy storage options to overcome the intermittent nature of these sources. Underground hydrogen storage (UHS) in depleted hydrocarbon reservoirs offers a scalable and practical energy storage solution. These reservoirs are chosen for their availability and large capacity, but the unique properties of hydrogen raise concerns about potential leakage pathways, particularly through wellbores. In this study, we develop and apply, for the first time, reduced-order models (ROMs) specifically designed for efficient leakage risk prediction in UHS systems operating in depleted hydrocarbon reservoirs. Using 3,000 high-fidelity simulation scenarios, we examine the influence of 11 key parameters, including reservoir and aquifer depths, wellbore permeability and porosity, initial saturations of water, oil and gas fractions (hydrogen, light, intermediate, and heavy hydrocarbons), reservoir pressure multiplier, and the aquifer-to-reservoir volume ratio, to simulate leakage behavior over a 1,000-year timescale. We train ROMs using a two-step classification-regression approach, achieving R 2 values exceeding 99 % across all targets. These ROMs effectively capture the leakage evolution and identify critical controls of leakage, guiding the design of mitigation strategies. Results indicate that gas leakage occurs in about 27 % of scenarios as early as five years post-operation, reaching volumes of up to 106 ft3. Oil leakage is less frequent (~17 %) and typically begins decades later. Our findings also show that hydrogen often migrates first, owing to its smaller molecular size and higher buoyancy, followed by heavier hydrocarbons. Over time, these heavier components contribute significantly to the total leaked volume, reinforcing the need for targeted monitoring and remediation strategies. Our analysis highlights that deeper storage reservoirs, shallower aquifers, and low-permeability wellbores significantly reduce leakage risks. In conclusion, this work offers a robust framework for risk-informed UHS deployment, supporting energy security through reliable large-scale hydrogen storage while safeguarding environmental integrity.

08 HYDROGEN↗
DOE OSTI2022

Large-Area Lasing in Nanoscale Complex Media: The Critical Role of Local Dielectric Environment

Controlling how electromagnetic waves interact with complex media is critical for applications in imaging and focusing. Such lightwave interactions with complex media can lead to dramatic optical effects like lasing. While much work in random lasing focus on understanding how gain and scattering co-operatively generate lasing, little work has focused on how to manipulate the lasing threshold without modifying the structural disorder. Here, a simple, mostly unexplored, strategy is demonstrated that employs atomic layer deposition (ALD) to tune the local near-field environment while preserving the underpinning disorder—controlling lasing in a nanoscale complex medium on a large scale (>cm 2 ). The nanoscale complex medium is a quasi-2D system of coupled zinc oxide nanospheres with overall thickness deep in the sub-wavelength regime (≈λ/4). Near-ultraviolet femtosecond spectroscopy probes the broadband response of the gain nanomaterial, details how ALD process fundamentally modifies the fast-picosecond and slow-nanosecond carrier dynamics, and informs on the relevant timescales critical for lasing. Full-field electromagnetic simulations provide critical insights about how near-field dielectric environment modifies the nanostructure's scattering cross-section, which ultimately results in enhanced lasing. Importantly, these results highlight a simple path to control how electromagnetic waves interact in a complex medium, a key step toward large-scale implementation of complex lasers.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Thermal and electric multidomain dynamic model for integration of power grid distribution with behind-the-meter devices

As renewable energy sources like solar and wind power become more integrated into the grid, coordinated control of behind-the-meter devices is crucial for enhancing grid flexibility and reliability and for meeting cost targets, with standardized models being developed to support this transition. The increasing flexibility and uncertainty of integrated renewable energy grids, along with interactions between various subsystems, make traditional steady-state modeling insufficient to capture transient and dynamic behaviors. Current models (e.g., composite load and battery equivalent models) focus on thermodynamic or electrical characteristics but overlook critical electromechanical interactions. This limits the ability to share performance information for grid services and hampers fast dynamic simulations. In addition, motor stalling is usually triggered by a fault event and attributed to the characteristics of the mechanical torque of the motor, resulting in absorption of a large amount of reactive power during the stalling period. Further, this significant withdrawal of reactive power will deteriorate the dynamic voltage stability of power grids and cause delayed voltage recovery. Therefore, an in-depth modeling of the thermodynamics or mechanical torque is essential to study the impacts of the realistic torque characteristics of those behind-the-meter devices on power system voltage stability. This study developed a dynamic multidomain model for building HVAC systems, such as air-source heat pumps, to simulate their thermal and electrical responses to grid transients. The model can accurately predict power metrics with a mean absolute percentage error of 10%, by validating against with power system computer-aided design performance data. Case studies demonstrate the model capability of capturing the transient response to sudden voltage changes, rapid load fluctuations, and system shutdowns respectively. During a sudden voltage drop (30% for 0.1s), a fully loaded heat pump’s motor speed dropped, continued declining, and shut down after 3.6s, with severe power oscillations and a torque spike. A partially loaded unit experienced temporary oscillations but stabilized. Under higher building loads, compressor speed increased from 64% to 100%, with power and torque rising before stabilizing. In safety-triggered shutdowns, power decreased after minor fluctuations, and torque briefly spiked before dropping to zero.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Cathodic Protection Modeling for Hanford Underground Double-Shell Tank Farms

Hanford stores millions of gallons of radioactive and chemically hazardous waste from the production of weapon materials in tank farms consisting of underground carbon-steel storage tanks surrounded by reinforced concrete. Six of these Hanford tank farms use double-shell storage tanks (DSTs). The DST farms were constructed from 1968 to 1986 with a planned 40–50 year design life, so some are already operating beyond their initial life expectancy. Ultrasonic testing (UT) has indicated significant thinning on the bottom of the secondary (outer) liner of these tanks, believed to arise from groundwater intrusion driving concrete side corrosion. There is no direct access to the steel/concrete interface between the tank and the concrete pad, making it difficult to apply a chemical-based mitigation strategy or to conduct repairs, but cathodic protection (CP) is a possible method to inhibit further concrete-side corrosion. Hanford already uses CP to protect below grade steel piping within the tank farms and connected to the tanks, but this system was not designed to protect the tank bottoms. CP design must account for the structures surrounding the DSTs, including the steel reinforcing bars (rebar) within the concrete pad and vault, various process lines, and the existing CP system. In this study, finite element analysis (FEA) modeling was carried out to simulate CP protection of 1) a single tank and CP anode to develop options for modeling the rebar and to compare to a simpler circuit model and 2) the entire Hanford AN tank farm as a representative example consisting of seven tanks, associated piping, and both existing and new CP anodes. Both circuit and FEA models predict that significant protective current could be delivered to the bottoms of the tanks with the addition of tank-protection anodes below the depth of the tanks. Simulations with only the existing pipe-protection anodes active confirmed that only a very small current to the tank bottoms is predicted under present conditions. Multiple simplified representations of the dome and wall rebar were tested to reduce the computational complexity of the tank-farm simulations, resulting in modeling the rebar as edge elements with a prescribed effective circumference that matches the real rebar surface area. The geometry of the rebar is also simplified into horizontal hoops around the tank walls and radial rebar over the dome with increased effective circumference to retain the target surface area. This simplification was found to greatly reduce the complexity and solution time of the models without large changes in current distributions, especially to the tank bottom. A range of values were tested for model parameters such as soil and concrete resistivities and polarization resistance to investigate their impact on the current and electric potential distributions. Depending on the parameters used, FEA simulations predict some risk of overprotection, particularly on the piping system; since overprotection can also lead to surface damage associated with hydrogen gas generation at the interface (e.g. hydrogen embrittlement or damage to coatings), this needs to be considered when refining the design of the new CP system. Comparison between the FEA models and the circuit model representation demonstrated that the circuit model could not match the predicted FEA current distribution, even when using the exact same surface areas. This discrepancy appeared to be at least partly attributable to the impact of the relative positions of the tank components and anodes to each other and to the ground surface. The FEA model accounts for the relative positions since it solves the governing equations in three dimensions, but the circuit model cannot account for the positioning. In particular, the circuit model underpredicts the current to the tank bottom and overpredicts the current to the dome compared to FEA for the baseline geometry. The FEA models omitted the electrically isolated rebar in the bottom concrete slab. However, a circuit based stray current model estimated that only 2.1% of the total current through the slab would stray into the rebar, corresponding to ~0.21 A for a target current density of 2 mA/ft2 to the tank bottom. The estimated corrosion driven by this amount of stray current is predicted to yield a lifetime of >400 years for the minimum rebar diameter, assuming an acceptable cross-section area loss of 10%.

d'Entremont, Anna [Savannah River National Laborat↗
DOE OSTI2025

Designing reinforcement learning algorithms for building HVAC control: From experimental observation to simulation comparisons

Advanced supervisory-level control with reinforcement learning (RL) is regarded as a promising solution for HVAC systems to minimize energy consumption while maintaining thermal comfort and indoor air quality. However, most RL applications were conducted in the simulation environment rather than real-world HVAC systems. This paper developed a value-based RL controller termed Deep Q-Network (DQN) for a typical central HVAC system and evaluated its performance in a building test facility. By comparing DQN with a rule-based controller, the study not only demonstrated the cases where DQN could properly maintain indoor comfort but also discussed possible reasons why DQN failed in some other situations. Recognizing the limitations of value-based RL algorithms from the experimental tests, a simulation study was conducted to compare DQN with an alternative RL approach, an actor–critic algorithm termed Deep Deterministic Policy Gradient (DDPG). In scenarios with a relatively large action space, DDPG outperformed DQN by requiring fewer computational resources and achieving better thermal comfort, lower energy consumption, and more stable control actions. The findings suggest that the ability of DDPG to handle continuous control variables more effectively allows for faster convergence in training and more precise control in practice, which enhances the overall efficiency and reliability of the HVAC system.

Guo, Fangzhou↗
DOE OSTI2026

Optimizing district energy systems by integrating Borehole Thermal Energy Storage Using a Mixed-Integer Linear Programming g-function framework with a Multi-Timescale Rolling Horizon method

Shallow geothermal has gained increasing attention in recent years; however, a reliable framework for its accurate incorporation into large-scale energy system optimization remains lacking. This study proposes a Mixed-Integer Linear Programming (MILP) framework combined with the g-function approach to integrate Borehole Thermal Energy Storage (BTES) technology into energy system optimization. Validation against a Modelica-based reservoir network simulation demonstrates that the proposed framework effectively captures the ground thermal response under varying energy loads and accurately estimates the borefield energy supply. To enhance scalability, a Rolling Horizon with Multi-Timescale (RH-MTS) method is further introduced, reducing computational time by 73 % for the 1-year optimization model with only minor loss of optimality. The framework is demonstrated through the case study of the UC Berkeley campus. Results indicate that BTES is a cost-effective and low-carbon solution: two borefields comprising 382 boreholes can meet 8.0 % and 6.6 % of the total campus heating and cooling demand, respectively, at an average energy rate of 0.70–0.77 USD/kWh and carbon intensity of 0.54 kg-CO2/kWh. Short-term analysis reveals a 35%–65% decline in BTES energy flow after 3–6 months of continuous heating/cooling operation, while long-term simulation shows that annual energy production of BTES can vary by up to 12.0 % after four years before stabilizing. Overall, this study develops a novel optimization framework that couples physics-based g-function method with MILP optimization framework, thereby advancing methodological development for shallow-geothermal integration and providing actionable guidance for BTES deployment in district-energy systems.

Yang, Jiahui↗
DOE OSTI2024

Global peak water limit of future groundwater withdrawals

Over the past 50 years, humans have extracted the Earth’s groundwater stocks at a steep rate, largely to fuel global agro-economic development. Given society’s growing reliance on groundwater, we explore ‘peak water limits’ to investigate whether, when and where humanity might reach peak groundwater extraction. Using an integrated global model of the coupled human–Earth system, we simulate groundwater withdrawals across 235 water basins under 900 future scenarios of global change over the twenty-first century. Here we find that global non-renewable groundwater withdrawals exhibit a distinct peak-and-decline signature, comparable to historical observations of other depletable resources (for example, minerals), in nearly all (98%) scenarios, peaking on average at 625 km 3 yr -1 around mid-century, followed by a decline through 2100. The peak and decline occur in about one-third (82) of basins, including 21 that may have already peaked, exposing about half (44%) of the global population to groundwater stress. Most of these basins are in countries with the highest current extraction rates, including the United States, Mexico, Pakistan, India, China, Saudi Arabia and Iran. Furthermore, these groundwater-dependent basins will probably face increasing costs of groundwater and food production, suggesting important implications for global agricultural trade and a diminished role for groundwater in meeting global water demands during the twenty-first century.

58 GEOSCIENCES↗
DOE OSTI2026

Quantum effects on the dynamics and properties of soft materials

The quantum effects of nuclear and electronic motion play an important role in the structure, dynamics, and function of soft materials, yet they are difficult to capture with conventional classical simulations or static electronic–structure methods. In this work several complementary approaches for treating quantum effects in polymeric and soft–matter systems are demonstrated, with a focus being on the hydrogen-bonded networks, ion and charge transport, and photoactive chromophores. The proton transfer, tunneling, and isotope effects are captured within the reduced-dimensionality models by implementing grid-based nuclear quantum dynamics in terms of the discrete variable and Fourier bases. The nuclear quantum dynamics is extended to larger systems by employing the quantum trajectories and quantum–thermal bath schemes combined with on-the-fly electronic structure, enabling the description of high-dimensional polymeric environments at feasible cost. The dynamics in the electronic degrees of freedom, simulating the optical response in large chromophores such as chlorophylls, is performed using the real-time time-dependent density functional theory implemented in the real-space multigrid (RMG) code. These approaches are demonstrated on case studies of the proton and hydroxide transport in hydrated polymer membranes, charge transfer in conjugated polymers, and the optical spectra of chlorophyll chromophores relevant to polymerized chlorophyll materials and chlorophyll–polymer hybrids. The reviewed methods and applications highlight practical routes of including quantum effects in simulations of soft functional materials.

Garashchuk, Sophya [Univ. of South Carolina, Colum↗
DOE OSTI2021

Avalanche transport of energetic-ions in magnetic confinement plasmas: nonlinear multiple wave-number simulation

Large burst activity, identified as toroidal Alfvén eigenmode (TAE) avalanche, occurs frequently in neutral-beam heated plasmas in National Spherical Torus Experiment (NSTX). Based on the typical experimental observation of TAE avalanche on NSTX, a self-consistent nonlinear multiple wave-number (k ∥ ≃ n/R, where n toroidal mode-number and R major radius) simulation associated with TAE avalanches is performed using the experimental parameters and profiles before the occurrence of TAE avalanche as the M3D-K input. The wave–wave nonlinear coupling among different modes and the resonant interaction between different modes and energetic-ions during TAE avalanches are identified in the nonlinear multiple wave-number simulations. The resonance overlap during the TAE avalanche is clearly observed in the simulation. It is found that the effective wave–wave coupling and a sufficiently strong drive are two important ingredients for the onset of TAE avalanches. TAE avalanche is considered to be a strongly nonlinear process and it is always accompanied by the simultaneous rapid frequency-chirping and large amplitude bursting of multiple modes and significant energetic-ion losses. The experimental phenomenon is observed on NSTX and is qualitatively reproduced by the simulation results in this work. These findings indicate that the onset of avalanche is triggered by nonlinearity of the system, and are also conducive to understanding the underlying mechanism of avalanche transport of energetic particles in the future burning plasmas, such as International Thermonuclear Experiment Reactor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Microscopic Theory of Fluctuating Hydrodynamics in Nonlinear Lattices

The theory of fluctuating hydrodynamics has been an important tool for analyzing macroscopic behavior in nonlinear lattices. However, despite its practical success, its microscopic derivation is still incomplete. In this work, we provide the microscopic derivation of fluctuating hydrodynamics, using the coarse-graining and projection technique; the equivalence of ensembles turns out to be critical. The Green-Kubo (GK)-like formula for the bare transport coefficients are presented in a numerically computable form. Our numerical simulations show that the bare transport coefficients exist for a sufficiently large but finite coarse-graining length in the infinite lattice within the framework of the GK-like formula. Finally, this demonstrates that the bare transport coefficients uniquely exist for each physical system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗