Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “process coupling”

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

Model Data Archive for Manuscript Titled "Evaluation of a Coupled Surface–Subsurface Hydrologic Model Using Dense Water‑Level Sensors in a Mixed Urban–Rural Watershed"

This archive provides scripts, input files, and datasets used for the implementation and evaluation of a fully coupled surface–subsurface hydrologic model in the Neches River Basin, southeast Texas. The study uses the Advanced Terrestrial Simulator (ATS) to simulate coupled surface–subsurface hydrologic processes over a mixed urban–rural watershed and evaluates model performance using a dense network of 136 in situ water-level sensors, nine U.S. Geological Survey (USGS) stream gauges, and SSEBop-derived evapotranspiration estimates during the period October 2014–June 2024. The workflow is implemented primarily in Python 3 using the Watershed Workflow package. The Jupyter notebooks can be executed using open-source software such as Anaconda JupyterLab or Visual Studio Code. Other data files include TXT, CSV, XML, SHP, TIF, NetCDF, HDF5, and ExodusII files, which can be processed using the provided Python scripts. ATS input files are provided in XML format and can be edited using any commonly used text editor. This archive contains: *Scripts and input files used to generate the ATS model setup, including watershed discretization, mesh generation, parameter mapping, and model configuration. *Jupyter notebooks used for preprocessing observational data, evaluating streamflow, water levels, and evapotranspiration, computing performance metrics, and generating the figures presented in the manuscript. *ATS simulation outputs and processed observational datasets, including OneRain and DD6 water-level sensors, USGS streamflow observations, GIS data, and supporting spatial datasets used throughout the study.

Dense water-level sensor network↗
DOE OSTI2023

Reactive Transport Modeling of Anthropogenic Carbon Mineralization in Stacked Columbia River Basalt Reservoirs

Numerical simulation of CO2 storage in basalts and related reactive lithologies requires modeling complex, coupled hydrologic and chemical processes, including multi-phase flow and transport, partitioning of CO2 into the aqueous phase, and chemical interactions with aqueous fluids and rock minerals. We conducted reactive transport simulations of the Wallula pilot-scale CO2 injection into the flow tops of the Grande Ronde Basalt using our PNNL STOMP-CO2 simulator with the ECKEChem reactive module. Our mineralization simulation of the ~1,000 tons of injected CO2 into the interflow zones was based on the hydrologic transport model we previously developed. For this work, the simulations considered geochemical reactions involving the basalt components, precipitates, formation brine, and injected CO2. In our benchmark case, carbonate minerals precipitated, resulting in ~20% of the CO2 being mineralized in 10 years. Increasing the reaction rate of a single primary mineral phase (clinopyroxene) by an order of magnitude resulted in a carbon mineralization reaction extent of ~90% over the same time interval. Based on these initial sensitivity analysis results, it is clear that a thorough understanding of primary mineral dissolution rates is required for accurately predicting long-term fate and transport of injected CO2 into basalt formations. Our reactive transport numerical simulations will be key components of commercial-scale CO2 storage operation permitting, de-risking, and optimization in mafic and ultramafic reservoirs.

Cao, Ruoshi↗
DOE OSTI2023

Understanding and Predictability of Integrated Mountain Hydroclimate (Workshop Report)

Mountainous systems cover approximately 23% of Earth’s land and are distributed across all continents. They can capture and store atmospheric moisture that is then cycled through the terrestrial surface and subsurface system, released to downstream communities, and cycled back to the atmosphere. Mountain hydroclimate—characterized by steep gradients, geological, ecological, and biogeochemical diversity—is influenced by topographic forcing and elevated warming and susceptible to large subseasonal to multidecadal variability and rapid changes. Terrestrial hydrological and biogeochemical cycles also experience cascading effects from global warming impacts, such as multidecadal declines in mountain snowpack, longer growing seasons, and increased frequency and severity of extreme events like droughts and wildfires. However, little is known about the effects of these impacts and their feedbacks on climate systems and surface-subsurface compartments. Also unknown are the full implications of changing hydroclimate and extreme events on hydro biogeochemical cycles across atmosphere, terrestrial, and human systems in mountain regions and beyond. This knowledge gap is critical, given human reliance on mountain systems for stable water supply and quality. Mountain systems’ increasing vulnerability to climate change and human perturbations motivates the need to improve understanding of integrated mountain hydroclimate (IMHC) systems and their feedbacks and impacts on humans across scales. However, due to large heterogeneity and strong gradients, coupled natural-human processes in mountain regions present significant challenges for observations, modeling, predictions, and projections. Motivated by gaps in mountain hydroclimate understanding, observations, and modeling and the need for credible projections of future changes, the U.S. Department of Energy’s (DOE) Biological and Environmental Research (BER) program organized a virtual workshop on “Understanding and Predictability of Integrated Mountain Hydroclimate.” Sponsored by BER’s Earth and Environmental Systems Sciences Division (EESSD), the workshop aimed to inform and catalyze EESSD’s growing interests in enhancing predictive understanding of IMHC. Organizers structured the workshop to identify (1) knowledge gaps, (2) observational and modeling challenges, (3) short-term (1 to 3 years) and long-term (10 years and beyond) research opportunities, and (4) strategies for fostering collaboration and coordination. To address the outstanding challenges of IMHC, the workshop included two sessions organized by disciplinary, cross-disciplinary, and crosscutting science topics. The disciplinary and cross-disciplinary topics focused on essential IMHC elements: atmosphere, terrestrial, and human systems and their interactions. Breakout sessions on disciplinary and cross-disciplinary topics facilitated identification of crosscutting topics and central emerging themes. Session 1 focused on connecting existing DOE investments to accelerate progress related to scientific challenges in understanding mountain hydroclimate. In Session 2, participants further explored key Session 1 takeaways through the lens of multiagency collaborations and coordination.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

An Innovative Zero-Liquid Discharge Intermediate-Cold-Liquid Eutectic-Freeze Desalination System

A novel freeze desalination technology is presented for purification of water produced from various industrial processes, including oil and gas extraction. The main innovative idea is introducing a water-immiscible intermediate-cold-liquid (ICL) medium which absorbs the cold from a refrigeration cycle and transfers it to a freezing chamber where it is mixed with brine. The cold ICL enters the freezing chamber at temperatures between -10°C to -30°C. Precooled brine in injected into the freezing chamber and pure water freezes out of the brine upon mixing with the ICL. The rate of freezing inside the chamber is controlled by adjusting the flow rate of the cold ICL relative to the inlet brine. Both eutectic freeze and non-eutectic freezing are possible by controlling the temperature of the freezing chamber. The solid-liquid slurry created in the freezing chamber flows to a filter separation unit after passing through a separation column. The ice crystals are separated from the liquid by a rotary filter device. The collected ice is rinsed by a centrifugal unit. The rinsed ice is melted to produce fresh water. The developed technology resolves the issues common to current freeze-desalination systems while offering superior heat transfer performance ascribed to direct-contact between the cooling medium and the brine in addition to producing solid salt and pure water. The system operates under atmospheric pressure. The low operation temperature minimizes the corrosion issues common in systems exposed to brine. Also, due to the low operating temperature, there are no concerns about evaporation of volatile organics from the produced water during the treatment. The developed technology effectively eliminates the use of chemicals for pretreatment of brine. The system can be constructed from commercially available components. Extensive testing using a lab-scale prototype by using synthetic brine and actual produced water with TDS values from 50,000 ppm to 300,000 ppm showed that freshwater could be achieved with relatively small dependency on the quality of the input brine. The quality of the treated water was found to be mainly a function of the duration of the centrifugal process. Treated water with TDS as low as 200 ppm could be achieved in the tests. Thermo-economic modeling coupled with system level thermochemical process model were developed for prediction of the LCOW. The LCOW was predicted to be $\$$0.50 to $\$$0.90 per barrel of produced water. The potential markets for deployment of the technology are the oil and gas, in-land thermoelectric powerplants, RO desalination plants, and mining industry.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

Research in Theoretical Particle Physics: Neutrino Masses and The Origin of CP-Violation

The main goal of this project was to investigate the theories of neutrino masses at the low scale and the origin of CP-violation in physics beyond the Standard Model. The origin of neutrino masses is one of the most pressing issues in particle physics. In this proposal, we investigated the theories for neutrino masses based on local total lepton number. In these theories anomaly cancellation predicts the existence of extra fermions with lepton number, and one of them is a good dark matter (DM) candidate. The cosmological constraints on the DM relic density implies that the lepton number symmetry breaking scale must be below the multi-TeV scale. Therefore, one can hope to test the origin of neutrino masses in current or future experiments. We investigated in great detail the anomaly cancellation in these gauge theories, study the different mechanisms for neutrino masses, study the predictions for direct and indirect dark matter experiments. Since these theories predict new sources of CP violation, we investigated the predictions for the electric dipole moments (EDM). The Higgs decays and the different signatures at the Large Hadron Collider were investigated in great detail. Finally, we studied the possible baryogenesis mechanisms in these theories in agreement with the EDM, DM and collider constraints. The origin of CP violation in the Standard Model is unknown. In this proposal, we investigated the simplest mechanisms for spontaneous CP-violation to explain the CP-violation in the CKM matrix and the value of the QCD vacuum angle. We discussed the Nelson-Barr mechanism in gauge theories predicting vector-like quarks from anomaly cancellation such as theories for local baryon number. We discussed the Bento-Branco-Parada mechanism in gauge theories to explain the CP-Violation in the CKM matrix. We will investigate models with two Higgs doublets in the context of the minimal theory for quark-lepton unification. We will show the non-decoupling effects in the Higgs sector, the main constraints coming from flavour violating processes when the Yukawa couplings are related by the gauge symmetry in theories for quark-lepton unification. We will investigate the predictions for electric dipole moments in these theories. The possibility to have successful baryogenesis was investigated. Finally, we investigated the relation between CP-violation in the quark and leptonic sectors. The future results from these studies could help us to understand two main issues in physics beyond the Standard Model: The Origin of Neutrino Masses and CP-violation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Nitric Oxide Production from Nitrite plus Ascorbate during Ischemia upon Hippocampal Glutamate NMDA Receptor Stimulation

Nitric oxide (•NO), a diffusible free radical, is an intercellular messenger, playing a crucial role in several key brain physiological processes, including in neurovascular coupling (NVC). In the brain, glutamatergic activation of the neuronal nitric oxide synthase (nNOS) enzyme constitutes its main synthesis pathway. However, when oxygen (O2) supply is compromised, such as in stroke, ischemia, and aging, such •NO production pathway may be seriously impaired. In this context, evidence suggests that, as already observed in the gastric compartment, the reduction of nitrite by dietary compounds (such as ascorbate and polyphenols) or by specific enzymes may occur in the brain, constituting an important rescuing or complementary mechanism of •NO production. Here, using microsensors selective for •NO, we show that nitrite enhanced the •NO production in a concentration-dependent manner and in the presence of ascorbate evoked by N-methyl-D-aspartate (NMDA) and glutamate stimulation of rat hippocampal slices. Additionally, nitrite potentiated the •NO production induced by oxygen-glucose deprivation (OGD). Overall, these observations support the notion of a redox interaction of ascorbate with nitrite yielding •NO upon neuronal glutamatergic activation and given the critical role of NO as the direct mediator of neurovascular coupling may represents a key physiological mechanism by which •NO production for cerebral blood flow (CBF) responses to neuronal activation is sustained under hypoxic/acidic conditions in the brain.

Nunes, Carla (ORCID:0000000274521592)↗
DOE OSTI2021

Plasmonic Core–Shell–Satellites with Abundant Electromagnetic Hotspots for Highly Sensitive and Reproducible SERS Detection

In this work, we develop a Ag@Al2O3@Ag plasmonic core–shell–satellite (PCSS) to achieve highly sensitive and reproducible surface-enhanced Raman spectroscopy (SERS) detection of probe molecules. To fabricate PCSS nanostructures, we employ a simple hierarchical dewetting process of Ag films coupled with an atomic layer deposition (ALD) method for the Al2O3 shell. Compared to bare Ag nanoparticles, several advantages of fabricating PCSS nanostructures are discovered, including high surface roughness, high density of nanogaps between Ag core and Ag satellites, and nanogaps between adjacent Ag satellites. Finite-difference time-domain (FDTD) simulations of the PCSS nanostructure confirm an enhancement in the electromagnetic field intensity (hotspots) in the nanogap between the Ag core and the satellite generated by the Al2O3 shell, due to the strong core–satellite plasmonic coupling. The as-prepared PCSS-based SERS substrate demonstrates an enhancement factor (EF) of 1.7 × 107 and relative standard deviation (RSD) of ~7%, endowing our SERS platform with highly sensitive and reproducible detection of R6G molecules. We think that this method provides a simple approach for the fabrication of PCSS by a solid-state technique and a basis for developing a highly SERS-active substrate for practical applications.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Energetic Electrons Accelerated and Trapped in a Magnetic Bottle above a Solar Flare Arcade

Where and how flares efficiently accelerate charged particles remains an unresolved question. Recent studies revealed that a “magnetic bottle” structure, which forms near the bottom of a large-scale reconnection current sheet above the flare arcade, is an excellent candidate for confining and accelerating charged particles. However, further understanding its role requires linking the various observational signatures to the underlying coupled plasma and particle processes. Here we present the first study combining multiwavelength observations with data-informed macroscopic magnetohydrodynamics and particle modeling in a realistic eruptive flare geometry. The presence of an above-the-loop-top magnetic bottle structure is strongly supported by the observations, which feature not only a local minimum of magnetic field strength but also abruptly slowing plasma downflows. It also coincides with a compact above-the-loop-top hard X-ray source and an extended microwave source that bestrides the flare arcade. Spatially resolved spectral analysis suggests that nonthermal electrons are highly concentrated in this region. Our model returns synthetic emission signatures that are well matched to the observations. The results suggest that the energetic electrons are strongly trapped in the magnetic bottle region due to turbulence, with only a small fraction managing to escape. The electrons are primarily accelerated by plasma compression and facilitated by a fast-mode termination shock via the Fermi mechanism. Our results provide concrete support for the magnetic bottle as the primary electron acceleration site in eruptive solar flares. They also offer new insights into understanding the previously reported small population of flare-accelerated electrons entering interplanetary space.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Redox Poise during Rhodospirillum rubrum Phototrophic Growth Drives Large-scale Changes in Macromolecular Synthesis Pathways

During photoheterotrophic growth on organic substrates, purple nonsulfur photosynthetic bacteria like Rhodospirillum rubrum can acquire electrons by multiple means, including oxidation of organic substrates, oxidation of inorganic electron donors (e.g., H2), and by reverse electron flow from the photosynthetic electron transport chain. These electrons are stored as reduced electron-carrying cofactors (e.g., NAD(P)H and ferredoxin). The overall ratio of oxidized to reduced cofactors (e.g., NAD(P)+:NAD(P)H), or ’redox poise’, is difficult to understand or predict, as are the cellular processes for dissipating these reducing equivalents. Using physics-based models that capture mass action kinetics consistent with the thermodynamics of reactions and pathways, a range of redox conditions for heterophototrophic growth are evaluated, from conditions in which the NADP+/NADPH levels approach thermodynamic equilibrium to conditions in which the NADP+/NADPH ratio is far above the typical physiological values. Modeling predictions together with experimental measurements indicate that the redox poise of the cell results in large-scale changes in the activity of biosynthetic pathways and, thus, changes in cell macromolecule levels (DNA, RNA, proteins, and fatty acids). Furthermore, modeling predictions indicate that during phototrophic growth, reverse electron flow from the quinone pool is a minor contributor to the production of reduced cofactors (e.g., NAD(P)H) compared to other oxidative processes (H2 and carbon substrate oxidation). Instead, the quinone pool primarily operates to aid ATP production. The high level of ATP, in turn, drives reduction processes even when NADPH levels are relatively low compared to NADP+ by coupling ATP hydrolysis to the reductive processes. The model, in agreement with experimental measurements of macromolecule ratios of cells growing on different carbon substrates, indicates that the dynamics of nucleotide versus lipid and protein production is likely a significant mechanism of balancing oxidation and reduction in the cell.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Shallow- and deep-convection characteristics in the greater Houston, Texas, area using cell tracking methodology

The convective lifecycle, from initiation to maturity and dissipation, is driven by a combination of kinematic, thermodynamic, microphysical, and radiative processes that are strongly coupled and variable in time and space. Weather radars have been traditionally used to provide various convective-cloud characteristics. Here, we analyzed climatological convective-cell radar characteristics to obtain and assess the diurnal cycles of three convective-cell types – shallow, modest deep, and vigorous deep convective cells – that formed in the greater Houston area, using the National Weather Service radar from Houston, Texas, and a multi-cell identification and tracking algorithm. The examined dataset spans 4 years (2018–2021) and covers the warm-season months (June to September) in those years. The analysis showed clear diurnal cycles in cell initiation (CI) consistent with the sea breeze circulation and showed diurnal and normalized lifetime relationships in cell evolution parameters (e.g., maximum reflectivity, echo-top height, Geostationary Operational Environmental Satellite-16 (GOES-16) channel 13 brightness temperature, and the height of maximum reflectivity). The cell evolution is well represented by relationships between (1) the height and value of the maximum radar reflectivity, (2) the minimum GOES-16 channel 13 brightness temperature and the maximum vertically integrated liquid, (3) the maximum reflectivity and columnar-average reflectivity, and (4) the echo-top ascent rate and cell lifetime. The relationships presented herein help to identify the cell lifecycle stages such as early shallow convection, vigorous vertical development, anvil development, and convective core dissipation. GOES-16 Aerosol Optical Depth values are also used as a proxy for cell initiation aerosol concentrations to investigate any potential relationships between initiation location and aerosol concentration. Overall, no significant relationships between initiation location and aerosol concentration were found for the three cell types investigated, but there are some minor differences in the pre-CI aerosol optical depth for vigorous deep convective cells.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Feedback on Forty-year Long Clean-up Operations of a Contaminated Soil for Environmental Purpose - 20026

One of the CEA facilities in France is the place of miscellaneous mapping and clean-up operations since the 1980's. The final purpose of the conducted work deals with the environmental remediation of the building. In France, the absence of a regulatory framework for the management of sites polluted by radioactive substances, especially as the lack of release thresholds, has led to develop different approaches of soil radiological characterization. The geostatistical approach is part of them. Preliminary investigations are an unavoidable step prior to clearance and remediation. First, historical analysis enables to define site perimeters to investigate as well as to precisely locate expected contaminated parts. For radioactive pollutants, the media (water, soil, air) and transfer routes participate in the definition of the investigation perimeter and thus also need to be identified through a geological study. At the same time, functional analysis leads to choose the best measurement means which could provide useful information when combined with visual inspection. An equipment well adapted to easy-to-measure radionuclides as gamma emitters, such as gamma probe or on-site gamma spectrometry, is usually selected as a non-destructive assay. Coupled with spatial positions, processed on-site radiological data can be mapped. Data processing consists in mapping and kriging before a geostatistical analysis in order to identify the zones of interest to be targeted. Selected zones are not only linked to high level count rates but also to low level ones so that a pollution-free reference is known in the area. This data processing also enables to categorize waste types according to their origin and contamination levels. An expected volume of different waste then follows as well as the waste characterization equipments. In-depth investigations, e.g. core drillings and samplings, naturally ensue to aim at remediation optimization through the assessment of environmental impact. After clean-up operations, non-destructive assays as well as destructive samplings enable to check the effective clearance. Following this methodology in the present remediation work, results of the first non-destructive assay campaign highlighted a long-lived transuranic radionuclides contamination in the concrete flagstone. Successive non-destructive assays were then performed in the defined zones of interest. Count rates measured with a surface probe were first mapped. The best localization of destructive assays, i.e. core drillings at a 1 meter depth, came from this surface mapping followed by a geostatistical analysis. The concrete flagstone was totally removed. Then, the depth of soil to remove under the flagstone was optimized from the results of radiological activities measured in core drillings samples. A new surface mapping was then drawn after this partial soil excavation. The sand and demolition rubble samples, homogeneously constituted from the flagstone fragments, were analyzed by gamma spectrometry. Meanwhile, waste drums and other containers were also measured by gamma spectrometry and passive neutron measurement devices. All these results have brought miscellaneous pieces of information. In this context, major issues appear in terms of physical constraints such as the premises tininess as well as of soil sampling techniques and in terms of measurement performances assessment and the representativeness of homogeneous samples. The choice of radiological soil and waste characterization devices is described in the light of enhanced performances by keeping in mind radiation protection requirements and a willing to always optimize waste categories (very low level waste) and volume. The improvement of the characterization methods appears through this forty-year long work while sharing the feedback of encountered difficulties. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Solvated Nuclear–Electronic Orbital Structure and Dynamics

Nonadiabatic dynamical processes such as proton-coupled electron transfer and excited state intramolecular proton transfer have been the subject of much research. One of the promising theoretical methods to describe these processes is the nuclear–electronic orbital (NEO) approach. This approach inherently accounts for nuclear quantum effects within quantum chemistry calculations, and it has recently been extended to directly simulate nonadiabatic processes with the development of real-time NEO methods. These processes can also be significantly dependent on the surrounding chemical environment, however, and capturing the effects of the environment is often necessary for analyzing experimentally relevant systems. This work couples the NEO density functional theory and real-time time-dependent density functional theory approaches with solvation through the polarizable continuum model. Here, the effects of this coupling are investigated for ground state properties, solvent-dependent vibrational frequencies, and direct excited state intramolecular proton transfer dynamics.

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

Assessments of advanced reactor heat supply to high temperature industrial unit operations: Heat Engines and Heat Pumps

In this report, the feasibility of employing Advanced Nuclear Reactors to supply the necessary high-quality heat to Industrial Processes is investigated. As a first step, the most relevant industrial applications that could benefit from this coupling were identified. Four industries including petroleum refining (distillation, catalyst regeneration, hydro-steam cracking), chemical polymers, steel manufacturing (metals annealing/soaking) and cement manufacturing (calcination of limestone and dolomite industries) were considered. An overview of the heat duties, i.e., the temperatures and thermal powers required by representative plants of the different processes, was provided. Secondly, the thermodynamic conditions of the steam produced by representative designs of Advanced Reactor concepts (Liquid Metal Fast Breeder Reactors, Molten Salt Reactors, and High Temperature Gas-cooled Reactors) were summarized. The comparison between the requirements of the industrial processes and the capabilities of the candidate nuclear power plants showed that only a few processes could benefit from the coupling. Although Advanced Reactors operate at higher temperatures than Light Water Reactors, the thermodynamic conditions of the generated steam flow rates are generally not suitable for the selected applications. At the same time, data indicate that most of the processes can be suitably supplied if the temperature of steam is elevated to 900 °C. One possible solution to boost the quality of the steam from the nuclear island involves the use of a heat pump. Based on a technology similar to that found in refrigerators and air conditioners, heat pumps extract heat from a source, elevate its temperature and transfer it to where it is needed. A simple numerical example illustrating the viability of this method is described. Specifically, the energy conversion cycle of a conventional Pressurized Water Reactor was extended to incorporate a compressor to raise the heating value of the steam produced by the reactor. Two reference configurations, i.e., one where the industrial process coupled to the nuclear power plant exhausts saturated liquid and the other where it exhausts saturated steam, were considered. The thermal efficiencies achievable by increasing the heating value of the reactor steam in this way are significantly greater than by direct electric heating.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Real-Time Time-Dependent Nuclear–Electronic Orbital Approach: Dynamics beyond the Born–Oppenheimer Approximation

The quantum mechanical treatment of both electrons and nuclei is crucial in nona- diabatic dynamical processes such as proton-coupled electron transfer. The nuclear electronic orbital (NEO) method provides an elegant framework for including nuclear quantum effects beyond the Born-Oppenheimer approximation. To enable the study of nonequilibrium properties, we derive and implement a real-time NEO (RT-NEO) approach based on time-dependent Hatree-Fock or density functional theory, in which both the electronic and nuclear degrees of freedom are propagated in a time-dependent variational framework. Nuclear and electronic spectral features can be resolved from the time-dependent dipole moment computed using the RT-NEO method. The test cases demonstrate the dynamical interplay between the quantum nuclei and the electrons through vibronic coupling. Moreover, vibrational excitation in the RT-NEO approach is demonstrated by applying a resonant driving eld, and electronic excitation is demonstrated by simulating excited state intramolecular proton transfer. This work shows that the RT-NEO approach is a promising tool to study nonadiabatic quantum dynamical processes within a time-dependent variational description for the coupled electronic and nuclear degrees of freedom.

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

DECOVALEX-2023: Task C Final Report

The Full-scale Emplacement (FE) heater experiment at the Mont Terri Underground Rock Laboratory (URL) was designed and conducted by Nagra to replicate an emplacement tunnel of Nagra’s reference repository design at 1:1 scale. Alongside testing the technical feasibility of constructing disposal tunnels, emplacing waste containers in the tunnels and then backfilling them, the main goals of the FE experiment are (1) to obtain a better understanding of the coupled effects of induced thermo-hydro-mechanical (THM) processes that may occur and (2) to validate existing coupled THM models (Müller et al., 2017). A key aspect of ensuring safety for repositories located in low-permeability rock involves minimizing any damage to the rock itself, thereby preserving its integrity and promoting a stable environment Amongst a number of processes that could damage the rock is the increase in pore pressure due to thermal loading caused by heat emitted from the waste. To reduce the potential damage of the rock, it is important to analyse the evolution of heat over time due to the heat load of the containers and assess possible consequences by coupled THM models. The aim of Task C of DECOVALEX-2023 was to build 3D numerical models of the FE experiment, focussing in particular on the heating induced pore pressure change in the Opalinus Clay. Data from a large number of sensors were available from the FE experiment for model comparison. These sensors measured temperature and relative humidity in the bentonite around the heaters, and temperature, pressure and displacement/strain in the surrounding Opalinus clay. Data were available from the start of excavation (April 2012) up to August 2020 for most sensors (more than 5 years from the start of heating in December 2014). To fulfil the overall aim of the task, the work was broken down into a number of steps, starting with simpler models to build confidence in each team’s approach and then moving to more complex models that better represent the FE experiment. Step 0 consisted of 2D benchmark models, gradually increasing the number of processes that are represented from thermal (T) only models in Step 0a, to coupled thermal hydraulic (TH) models in Step 0b with a representation of changing porosity, to coupled thermo-hydro-mechanical (THM) models in Step 0c, where porosity changes are calculated by the mechanical model. A detailed specification of processes, parameters, initial and boundary conditions was provided for this step, with the ambition that all teams would work towards close agreement in their model results, thus building confidence in the model implementations. vi It was not straightforward to achieve agreement between the teams, so additional steps (Step 0b2, 0b3, 0c2, 0c3) were added along with derivation of some analytical solutions against which the models could be compared. The reasons for the differences between teams were investigated and found to be caused primarily by different conceptual model assumptions (including temperature dependence of the thermal expansion of water), different model formulations (including porosity evolution) and differences in modelled domain sizes, boundary conditions and grid discretisation. This demonstrates that comparisons between multiple modelling teams and/or comparison with analytical results and experimental data are highly beneficial in providing an indication of uncertainty in model predictions. At the conclusion of Step 0, almost all teams had achieved a close agreement in model results and those that had not achieved an agreement knew the reason for this. Step 1 moved from 2D models to 3D models of the FE experiment without adding technical features like shotcrete or EDZ, and only considering the heating phase. Initially the 3D model was tightly specified to continue to build confidence in the model implementations (Step 1a). The results of Step 1a were compared to the data from the FE-experiment without the teams seeing the data. The teams were then provided with a sub-set of the data from the FE-experiment and invited to consider how best to use the large dataset for model comparison (Step 1b). Teams were then asked to use the data provided to calibrate their models, only changing material property values rather than adding features or processes to their models (Step 1c). In Step 1, teams were asked to only model the heating phase of the experiment, so pressure in the Opalinus Clay was reported as change in pressure since the initial conditions were specified rather than modelled. The change from 2D to 3D models was accompanied by an increase in the dispersion of results between the teams. Some of this was resolved during the task, but some remained and is potentially due to model discretisation. Calibration of parameters was useful in improving the fit of the models to the data but the remaining differences indicated that the models were missing features or processes. In Step 2, the teams were asked to update their models with additional features and processes as well as calibrating parameters to try and improve the fit of the models to the data. Teams were encouraged to represent ventilation of the open FE tunnel prior to backfilling with heaters and bentonite and in Step 2, the absolute pressure in the Opalinus Clay was compared between the teams. Teams took different approaches, but there was consideration of adding shotcrete and an EDZ into the model, representing stress change during excavation and different approaches to modelling ventilation of the FE tunnel. Overall, the documented results showed a very good agreement for temperature. The results for porewater pressure evolution showed a significant improvement for most teams compared to Step 1c with a good agreement to the measurements for several teams whereas some teams overpredicted the pressure increase and others overpredicted the drainage effect especially for the sensors close to the heater. Step 3 was an opportunity for teams to use the models developed in Step 1 and Step 2 to make predictions about the temperature and pressure changes that will be expected at the FE experiment over the next few years in light of the planned changes in thermal output of the heaters.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗