Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled processes”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Rapid scalable fabrication of roll-to-roll slot-die coated flexible perovskite solar cells using intense pulse light annealing

Perovskite solar cells have achieved high efficiencies at the lab scale using batch processes that often may not be scalable. Processes that can be integrated into a continuous manufacturing platform, such as roll-to-roll, where several processes can be carried out sequentially are necessary to develop a new cost-competitive solar technology. In this work, intense pulsed light thermal post-processing is coupled with scalable slot-die coating and air knife evaporation in a high-speed roll-to-roll process of up to 2 m min –1 . The optimization of IPL energy parameters, web speeds, solution flow rates, and concentration of precursor inks yielded a champion power conversion efficiency of 11.24% for a 1 cm 2 roll-to-roll coated device. The overall processing time for a single layer in this process is reduced from minutes to seconds and the overall footprint of the manufacturing process is decreased. Furthermore, this demonstrates the promising capabilities for intensed pulsed light to reduce manufacturing costs and paves a roadmap for scalable and efficient PSCs on a roll-to-roll platform.

14 SOLAR ENERGY↗
DOE OSTI2021

BETO 2021 Peer Review - Process Monitoring and Predictions of Biorefinery Performance

Online process monitoring coupled with rapid predictive tools will be essential to refineries of the future to provide real-time feedback and process control on new, renewable feeds and their accompanying processes and products. This project will provide refinery operators with tools to predict product component concentrations in minutes from online, slip stream mass spectra, allowing for rapid detection of off-specification product. We will arrive at a template for predictive tool generation through the development of a specific tool as a starting point -- co-processing of pyrolysis oils and vacuum gas oil over fluid catalytic cracking (FCC) catalysts.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

Tree root nutrient uptake kinetics vary with nutrient availability, environmental conditions, and root traits: a global analysis

Summary: Root nutrient uptake by trees is a critical process that couples carbon and nutrient cycling in forest ecosystems. Yet, root nutrient uptake traits are poorly constrained, and the dynamics of this process are often not represented in models reflecting sparse measurements and understanding of root nutrient uptake physiology that lags those of aboveground physiology in forest ecosystems. Here, we present a global dataset of published nutrient uptake capacity and affinity values for tree species, with the goal of describing global patterns and evaluating responses to environmental drivers and associations with root traits. Furthermore, the dataset contains observations for ammonium, nitrate, and phosphate uptake spanning 77 tree species. Nutrient uptake capacity and affinity varied by more than an order of magnitude for each nutrient. Notably, tropical forests are underrepresented in these observations. Nutrient uptake capacity was generally diminished under nutrient enrichment but enhanced with soil warming and root–mycorrhizal colonization. The magnitude and direction of these effects can depend on the duration of exposure to a given treatment. Species with thinner roots had a tendency toward greater uptake capacity and affinity. Overall, root nutrient uptake traits are highly variable across tree species, yet they depend on environmental drivers and life-history strategies.

Km↗
DOE OSTI2022

Activation of Transition Metal (Fe, Co and Ni)-Oxide Nanoclusters by Nitrogen Defects in Carbon Nanotube for Selective CO 2 Reduction Reaction

The electrochemical carbon dioxide reduction reaction (CO2RR), which can produce value-added chemical feedstocks, is a proton-coupled-electron process with sluggish kinetics. Thus, highly efficient, cheap catalysts are urgently required. Transition metal oxides such as CoO x , FeO x , and NiO x are low-cost, low toxicity, and abundant materials for a wide range of electrochemical reactions, but are almost inert for CO 2 RR. Here, we report for the first time that nitrogen doped carbon nanotubes (N-CNT) have a surprising activation effect on the activity and selectivity of transition metal-oxide (MO x where M = Fe, Ni, and Co) nanoclusters for CO 2 RR. MO x supported on N-CNT, MOx/N-CNT, achieves a CO yield of 2.6–2.8 mmol cm -2 min -1 at an overpotential of -0.55V, which is two orders of magnitude higher than MO x supported on acid treated CNTs (MOx/O-CNT) and four times higher than pristine N-CNT. The faraday efficiency for electrochemical CO 2 -to-CO conversion is as high as 90.3% at overpotential of 0.44V. Both in-situ XAS measurements and DFT calculations disclose that MO x nanoclusters can be hydrated in CO 2 saturated KHCO 3 , and the N defects of N-CNT effectively stabilize these metal hydroxyl species under carbon dioxide reduction reaction conditions, which can split the water molecules and provide local protons to inhibit the poisoning of active sites under carbon dioxide reduction reaction conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Choosing the Best Modeling Platform for Radiological Risk Assessment Models - 20468

Radiological risk assessments, in the form of performance or safety assessments, are often required under regulations or guidance for remediation of contaminated land, decommissioning of contaminated buildings or structures, and radioactive waste disposal. These risk assessments are usually supported by fate and transport models that address decay and ingrowth of radionuclides, as well as their movement through engineered systems and the natural environment. These models are often projected thousands, or more, years into the future, largely because the radioactive species change through decay and ingrowth, and hence the magnitude of the radioactive effect changes with time. There are many computer codes that are available to address this type of modeling. They range from addressing specific pathways or processes such as infiltration of water, groundwater, surface water, air, biota, diffusion and advection of water and gases, to those that try to couple all processes together to evaluate the impact of fate and transport through the entire system to places in space and time to which access is assumed. These different types of codes are sometimes separated with the monikers process-level and systems-level codes, although it is often not clear that this separation does justice to the capabilities of the many codes that are available to evaluate fate and transport of radionuclides. The focus of this paper is the latter group of modeling codes. Several systems level modeling codes exist and are used. There are differences between these codes in terms of utility, flexibility, complexity and cost. The purpose of this paper is to compare a few of these codes in the context of work currently being performed by the International Atomic Energy Agency (IAEA) Modeling and Data for Radiological Impact Assessments (MODARIA) II Working Group 1 (WG1). The MODARIA II WG1's main focus is how stakeholder engaged decision analysis can, or should, be applied to radiological contamination problems so that better, longstanding, sustainable, solutions are reached. However, the WG1 also recognizes the potential impact of the modeling tools that are chosen to address radiological risk, which is often a primary objective of decision making for radiological problems. Other objectives might also be important, such as constraining costs, obtaining financing, minimizing impact on ecosystems, saving cultural resources, saving jobs, farmland, environmental justice, etc., in a full decision analysis for a given radiological contamination problem, but none of these other objectives have the same types of complex modeling needs as minimize radiological dose. Consequently, a further focus of WG1 is to evaluate the potential impacts on decision making of the choice of fate and transport, and risk assessment, modeling codes that are used to support decision making. The WG1 will produce a report at the end of 2020 that will focus on an approach to effective decision making and stakeholder engagement. The report will also consider the role that performance assessment modeling should play in the decision-making process, including the impact of the choice of modeling tools or computer codes on risk-informed decision making. Several sites around the World have been made available by Member States for these model comparisons, and several modeling tools have been considered. However, the focus of this paper is on two of the sites, one in Belgium and one in Ukraine, and on three of the tools: NORMALYSA (NORM And Legacy Site Assessment); GoldSim{sup C}, and AMBER{sup C}. The final report from this working group will also cover other modeling tools, including RESRAD, and PC-Cream{sup R}. (authors)

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

Chemo‐Mechanical Coupling in Hydrogels: Dynamics in the Diffusion‐Limited Regime

Hydrogels are characterized by substantial volume changes in response to external stimuli, making them promising candidates for developing smart materials with enhanced adaptability and responsiveness. By integrating chemical reactions, hydrogels acquire dynamic and tunable responsiveness to external stimuli through chemo‐mechanical coupling, expanding their potential in emerging applications. However, capturing their transient behavior remains challenging due to the complex interplay of chemical reactions, solvent transport, and polymer network deformation. Classical theories capture equilibrium swelling but fail to describe time‐dependent phenomena. To address this, a time‐dependent continuum model is developed that explicitly couples these processes. Volume phase transition in hydrogels with homogeneous chemical reactions is investigated, then the effects of reaction kinetics on these transitions are analyzed. The coupling of these mechanisms is further explored through a study of transient mechanical instabilities. To illustrate the impact of distinct reaction and diffusion timescales, a photo‐active gripper is studied for robotic applications. Finally, a photo‐active microswimmer that exhibits non‐reciprocal motion is proposed to highlight the solvent diffusion for locomotion at the micro‐ and nano‐ scales. The work establishes that transient dynamics of chemo‐mechanical hydrogels generate functions not accessible by steady state models and provides a predictive platform for designing adaptive materials in emerging applications.

chemo-mechanical coupling↗
NASA NTRS1981

The energetic ion composition spectrometer /EICS/ for the Dynamics Explorer-A

The energetic ion composition spectrometer (EICS) on the DE-A spacecraft is a high sensitivity, high resolution (M/Delta M greater than 10 FWHM at focus) instrument. It covers the entire mass range from less than 1 amu/e to greater than 150 amu/e in 64 mass channels, at each of 32 energy per charge steps covering the range from 0 to approximately 17 keV/e. The measurements obtained from this instrument will be used primarily to investigate the coupling between the magnetosphere and the ionosphere via the auroral acceleration region. Additional uses of the data will be to investigate the ring current loss processes, the coupling of the dayside cusp particles to the magnetosheath, and the tracing of artificially injected ions.

Shelley, E. G.↗
DOE OSTI2020

Minimizing Timing Jitter’s Impact on Ground-Penetrating Radar Array Coupling Signals

This article presents a novel investigation into the effect of receiver timing jitter on the quality of imaging of ground-penetrating radar systems. We explicitly show the impact of timing jitter throughout an entire example processing pipeline. The process of coupling removal in systems whose receivers have a random Gaussian-distributed jitter with a high standard deviation leaves randomly distributed residue in the radar image. This residue, called jitter-induced coupling noise, extends even into the expected target domain which may cause false alarms. This residue approximately follows a sum of Gaussian and Gamma distributions based on the higher order derivatives of the coupling signal. Once migrated to create an image, the coupling residue is best described by a Weibull distribution. We show that simple filtering does not remove enough residue and propose two preventative design specifications called the linear slope and migrated probability specifications. Antenna and system designers should use these specifications to assess the true severity of coupling signals in a given jitter environment.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Energy Efficient Thermoplastic Composite Manufacturing

The objective of the project was to establish an effective and affordable method to lay-up and consolidate/join large thermoplastic composite aerospace structure with cycle times measured in minutes rather than hours. Composite airplane designs have proven efficient and effective however future potential product production rates are challenging what the current systems can efficiently produce due to material lay-down constraints and extended thermal cycle times. The ability to lay-up then rapidly heat, consolidate, and cool large complex composite structures plus very accurately tool them (i.e. matching CTE of composite materials) along with very precise thermal control is a difficult challenge. Current systems such as autoclave processing of thermoset materials require long cycle times due to method of heating and the large associated thermal masses. These extended cycle times inhibit the ability to meet higher rate production scenarios due to the need for multiple sets of equipment and tools. Thermoplastic composite materials were used to facilitate more rapid cycle times via the elimination of a need for a cure dwell at temperature. Also, utilized induction heating along with smart susceptors to enable the quick cycle times needed while providing precise intrinsic thermal control. The intrinsic control of the induction heating process via the smart susceptors along with the use of laser assisted fiber placement were the key innovations areas developed. The ability to rapidly heat and then precisely control the temperature of the consolidation component was paramount to the success of this project. It is this unique processing attribute coupled with the thermoplastic material characteristics along with rapid lay-up methods that are novel and advantageous. While many requirements and influencing factors decide the materials and processes utilized for future airplane construction, the forecast of accelerated production rates and the recent performance successes of composites in airplane construction provide an opportunity for this processing technology to have significant influence. During the course of this project a number of key accomplishments were completed. They are as follows: The needed laser assisted fiber placement capabilities have been developed and validated on a component of meaningful scale and complexity. Subsequently, a thermoplastic preform of the part scale-up part design was laser assisted fiber placed and made availalbe for consolidation. An induction consolidation system was sized, designed, and fabricated resulting in the full capability to consolidate the scale-up component. The tool for induction consolidating of the scale-up component was designed and fabricated. An aluminum pressure bladder was welded and formed for use as the pressurization membrane to perform the needed application of pressure on the part at temperature to achieve consolidation. The induction consolidation system consisting of the restraint, induction power supply, induction tool with integrated induction coils and smart susceptor liners, along with the aluminum bladder were successfully used to consolidate a large thermoplastic skin. The capability of accomplishing rapid heat-up rates and cool down rates for consolidation of large thermoplastic skins with precise thermal control and even pressure application at the consolidation temperature was validated. The general conclusion is that the use of induction consolidation with smart susceptors is an effective method to rapidly and efficiently consolidate large thermoplastic composite components and enable affordable high rate manufacturing of composites. More specific conclusions generated by this project are as follows: A restraint type induction consolidation system provides an affordable option over the standard press when using the induction consolidation processing method. Aluminum bladders can be formed to the complex shapes needed and supply uniform pressure for consolidation and co-consolidation of thermoplastic composite materials. Straightforward modifications can be made to existing fiber placement robotic systems to enable their use for thermoplastic lay-up applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Charge transfer driven by ultrafast spin transition in a CoFe Prussian blue analogue

Photoinduced charge-transfer is an important process in nature and technology and is responsible for the emergence of exotic functionalities, such as magnetic order for cyanide-bridged bimetallic coordination networks. Despite its broad interest and intensive developments in chemistry and material sciences, the atomic-scale description of the initial photoinduced process, which couples intermetallic charge-transfer and spin transition, has been debated for decades; it has been beyond reach due to its extreme speed. In this work we study this process in a prototype cyanide-bridged CoFe system by femtosecond X-ray and optical absorption spectroscopies, enabling the disentanglement of ultrafast electronic and structural dynamics. Our results demonstrate that it is the spin transition that occurs first on the Co site within ~50 fs, and it is this that drives the subsequent Fe-to-Co charge-transfer within ~200 fs. This study represents a step towards understanding and controlling charge-transfer-based functions using light.

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

Using the optimal combined index weight ratio to improve the probability of anomaly detection in big area additive manufacturing

Big Area Additive Manufacturing (BAAM) of composites requires significant time, energy, and material, so it is critical to reduce production inefficiencies to make functional parts without multiple iterations. Statistical process control coupled with Principal Component Analysis (PCA) is a powerful technique that provides a quick, computationally inexpensive, and intuitive way for operators to detect defects that form in a manufacturing process without massive datasets. Recently, a combined index that is a weighted sum of the Hotelling's T 2 and squared residual error statistics has been proposed that can be monitored in one chart, improving interpretation accuracy and simplicity. However, the literature does not offer a formal method to optimise the weights. Here, we introduce two new approaches to the traditional weight selection approach using simulated and BAAM image data. Approach 1 uses a theoretically motivated optimum inspired by probabilistic principal component analysis. Approach 2 systematically varies the ratio of the weights to find the optimum. We show that approach 1 delivers optimal anomaly detection performance in select cases while approach 2 fares better in practice. Surprisingly, we also show that choosing a more complex PCA model has a minimal negative impact on anomaly detection performance compared to a more simplistic model.

3-dimensional printing↗
NASA NTRS1985

The upper atmosphere research satellite

The Upper Atmosphere Research Satellite (UARS) will provide, for the first time, data on a global basis for the study of the physical processes acting witin and upon the stratosphere, mesosphere, and lower thermosphere. Specifically, the areas of scientifc study to be addressed are energy input and loss, photochemistry, dynamics, and the coupling among processes and between atmospheric regions. The UARS is a single observatory consisting of a multimission modular spacecraft (MMS) and an instrument module containing 10 scientific instruments. The satellite will be Shuttle launched and placed in a 57 deg inclined orbit at 600 km altitude. A Central Data Handling Facility (CDHF) will receive data from the satellite and process these data into atmospheric quantities for use by the science team. The 'processed' data will be stored at the CDHF and will be available via communication lines for analysis by the investigators at their home laboratories using remote computers. Together with other satellite programs, balloons, sounding rockets, and laboratory efforts, UARS will make available the opportunity for extensive coordination of data devoted to solar terrestrial study.

Burr, P. T.↗
DOE OSTI2025

Development of an ERT‐Based Framework for Bentonite Buffers Monitoring From Laboratory Tests: 2. Quantitative Moisture Dynamics Estimation Model

Abstract The long‐term containment of high‐level radioactive waste in geological disposal repositories relies on Engineered Barrier Systems (EBS), with bentonite clay emerging as a candidate material due to its unique properties. Understanding moisture dynamics within bentonite buffers is crucial for EBS performance, as it directly influences the material's swelling capacity, thermal and hydraulic conductivity, mechanical properties, and long‐term evolution under complex thermal‐hydrological‐mechanical (THM) processes. This study develops an advanced Electrical Resistivity Tomography (ERT)‐based framework to quantitatively monitor moisture dynamics under THM conditions. Our framework extends the Waxman‐Smits model to incorporate the coupled effects of temperature, water content, fluid chemistry, and mechanical changes on bentonite's electrical properties. Utilizing HotBENT‐Lab data from our companion paper, which includes electrical conductivity, CT density, and thermocouple measurements, this study offers a novel methodological framework bridging different scales of the model. Our results show that the extended model can estimate water content from ERT data, capturing spatial and temporal variations in moisture distribution within bentonite columns. However, the model tends to overestimate water content compared to CT density‐derived measurements. We address this discrepancy by incorporating a simplified swelling effect model, which improves agreement between ERT and CT density‐based water content estimates. We also discuss model limitations, including simplified treatment of swelling and micropore effects, and propose a conceptual framework for transitioning from laboratory to field applications, addressing challenges such as parameter scalability, field validation methods, and integration of diverse data sources. This ERT‐based framework can potentially advance real‐world moisture monitoring of bentonite‐based EBS in nuclear waste repositories. Plain Language Summary Safely containing high‐level radioactive waste depends on barriers made from materials like bentonite clay, which is effective because it swells and seals in the waste. To ensure these barriers work well over time, it's important to understand how moisture moves through the clay. Our study developed a new method using ERT to monitor moisture levels in bentonite under conditions that mimic those in actual storage sites, including changes in temperature, water content, and mechanical stress. This study improved an existing model to better account for how these factors affect the clay, allowing us to create more accurate moisture maps. Initially, the proposed model overestimated the amount of water in the clay, but its accuracy was improved by factoring in how the clay swells when wet. This study also identified some limitations of the model and suggested ways to adapt it for use in real‐world waste storage sites. This new approach could lead to better monitoring and safety checks for nuclear waste storage systems, helping to ensure long‐term containment. Key Points This work develops an ERT‐based framework extending the Waxman‐Smits model to monitor bentonite moisture dynamics during coupled THM processes The extended model accurately estimates water content from Electrical Resistivity Tomography data, incorporating swelling effects to improve precision This work proposes a conceptual framework for transitioning from laboratory to field applications, advancing EBS monitoring in nuclear waste repositories

Chen, Hang↗
NASA NTRS1974

Study of ceramic products and processing techniques in space

An analysis of the solidification kinetics of beta alumina in a zero-gravity environment was carried out, using computer-simulation techniques, in order to assess the feasibility of producing high-quality single crystals of this material in space. The two coupled transport processes included were movement of the solid-liquid interface and diffusion of sodium atoms in the melt. Results of the simulation indicate that appreciable crystal-growth rates can be attained in space. Considerations were also made of the advantages offered by high-quality single crystals of beta alumina for use as a solid electrolyte; these clearly indicate that space-grown materials are superior in many respects to analogous terrestrially-grown crystals. Likewise, economic considerations, based on the rapidly expanding technological applications for beta alumina and related fast ionic conductors, reveal that the many superior qualities of space-grown material justify the added expense and experimental detail associated with space processing.

Markworth, A. J.↗
DOE OSTI2026

Data for Zheng et al. (2025), "AquaMEND: Reconciling multiple impacts of salinization on soil carbon biogeochemistry"

Soil salinization, exacerbated by climate change, poses a global threat to coastal ecosystems and soil function. Salinity affects soil carbon cycling by directly impacting microbial activity and indirectly altering soil physicochemical properties, but current models inadequately represent these complexities. This dataset contains the observational and modeling data from Zheng et al. (2025), which described a process-based modeling framework that couples soil solution chemistry with microbial carbon cycling reactions to study the impacts of soil salinization. This conceptual model is implemented numerically into the open-source geochemical program PHREEQC 3.0 (Parkhurst and Appelo, 2013). This dataset consists of: - Figure2_AquaMEND_salinity_buffer: Contains model simulation outputs to assess the impact of three different cation exchange and surface complexation processes on salinity buffering (Fig. 2 from Zheng et al. 2025). - Figure3_Salinity_function: Contains salinity function fitting for literature data (Fig. 3 from Zheng et al. 2025). - Figure4_AquaMEND_microbial_mechanisms: Contains model simulation outputs for testing various microbial process-based hypotheses related to soil salinization, including microbial mortality, carbon use efficiency (CUE), extracellular enzyme activity, and other microbial mechanisms (Fig. 4 from Zheng et al. 2025). - Figure5_AquaMEND_Redox: Contains on model simulation outputs to evaluate shifts among key redox processes, such as aerobic respiration, sulfate reduction, and methanogenesis (Fig.5 from Zheng et al. 2025). - Figure6_AquaMEND_sorption: Contains on model simulation outputs for investigating the effects of salinity on dissolved organic matter (DOM) sorption and desorption processes (Fig. 6 from Zheng et al. 2025). - Figure7_AquaMEND_process_couple: Contains on model simulation outputs for exploring coupled biotic-abiotic processes and their interactions (Fig. 7 from Zheng et al. 2025). - data: Includes datasets used to develop salinity response functions and evaluate salinity buffering capacity. Datasets for MEND model calibration. - database: Contains the `.dat` file required by PHREEQC for model execution. - README.md: A Markdown plain text file describing the computational tools and directories. Files are a mixture of plain text CSV (comma-separated value) and plain text *.dat files written by the model; no special software is required to read them.

EARTH SCIENCE > AGRICULTURE > SOILS > SOIL SALINIT↗
DOE OSTI2026

Data, model inputs, and analysis scripts associated with a manuscript on stream intermittency controls across spatial scales in Pacific Northwest watersheds

NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the manuscript "Hydroclimatic Memory and Watershed Template Shape Stream Intermittency: Multi-scale Attribution Using Process-based Simulation and Explainable ML" by Niroula et al. (2026), submitted to Water Resources Research (WRR). The study investigates the dominant controls on stream intermittency across local, reach, and watershed scales using a coupled process-based simulation and explainable machine-learning framework. Long-term daily simulations from the Advanced Terrestrial Simulator (ATS) were used to generate wetness states and ponded-depth responses over river-corridor cells. These ATS outputs were then aggregated across scales and used to train XGBoost (eXtreme Gradient Boosting) models. SHAP (SHapley Additive exPlanations) was applied to quantify the relative importance of hydroclimatic forcings, watershed template attributes, and antecedent-memory effects in shaping intermittency behavior. The analysis is carried out for three contrasting Pacific Northwest watersheds: Oak Creek (OCW), American River Watershed (ARW), and H.J. Andrews (HJA). Across these testbeds, the package contains ATS-ready watershed inputs, ATS run configuration and selected output files, model-evaluation data products, intermittency-analysis datasets, machine-learning target-feature tables, SHAP outputs, and notebooks used to organize, analyze, and visualize results. At a high level, the package documents a workflow in which ATS provides the physically based simulation backbone and explainable machine learning is used as a post-processing attribution tool. The contents are intended to support interpretation of the manuscript figures and results, provide context for how intermittency metrics were generated at multiple scales, and preserve the key artifacts needed to understand and reuse the analysis workflow. The package contains a high-level directory summary file (`summary.txt`) and four main content folders (1) `evaluation_plots` contains evaluation figures and supporting evaluation datasets; (2) `intermittency_plots` contains intermittency-focused analysis notebook and prepared datasets; (3) `ml-training-and-shap_values_plots` contains ML training inputs, SHAP outputs, and figure-generation notebooks; and (4) `watershed_mesh_and_ats_input` contains ATS model setup materials, forcing inputs, geometry, and selected run files. More specifically, the `evaluation_plots` folder contains the notebook used for ATS evaluation plotting and site-specific evaluation datasets. These include evapotranspiration and water-balance products for three watersheds, as well as an Oak Creek field-measurement discharge file. The `intermittency_plots` folder contains the notebook used for intermittency analysis and the prepared datasets used to analyze intermittent and non-intermittent wetness behavior across the study watersheds. The `ml-training-and-shap_values_plots` folder contains notebooks and outputs for the machine-learning and explainability workflow. This includes the main XGBoost and SHAP notebook(s), a beeswarm plotting notebook, target-feature tables for machine-learning training, SHAP summary tables, and per-sample SHAP value archives. The `watershed_mesh_and_ats_input` folder contains ATS-related watershed inputs and supporting materials. This includes mesh and shape products, ATS-readable LAI and meteorological forcing inputs, selected ATS spinup and transient-run files, and a watershed workflow example notebook. Subdirectories are organized by watershed where applicable.All files are .cpg (codepage files), .csv (comma-separated values), .dbf (database files), .exo (Exodus mesh format), .h5 (HDF5 format), .ipynb (Jupyter notebooks), .pkl (Python pickle), .prj (projection files), .sh (shell scripts), .shp (shapefile geometry), .shx (shapefile index), .txt (text files), or .xml (markup data).

Advanced Terrestrial Simulator↗
DOE OSTI2023

Implementation and verification of a user-defined element (UEL) for coupled thermal-hydraulic-mechanical-chemical (THMC) processes in saturated geological media

Efficient and accurate modeling of the coupled thermal-hydraulic-mechanical-chemical (THMC) processes in various rock formations is indispensable for designing energy geo-structures such as underground repositories for high-level nuclear wastes. Here, this work focuses on developing and verifying an implicit finite element solver for generic coupled THMC problems in geological settings. Starting from the mass, momentum, and energy balance laws, a specialized set of governing equations and a thermoporoelastic constitutive model is derived. This system is then solved by an implicit finite element (FE) scheme. Specifically, the residuals and the Jacobians are scripted in a user-defined element (UEL) subroutine which is then combined with the general-purpose FE software Abaqus Standard to solve initial-boundary value problems. Considering the complexity of the system, the UEL development follows a stepwise manner by first solving the coupled hydraulic-mechanical (HM) and thermal-hydraulic-mechanical (THM) equations before moving on to the full THMC problem. Each implementation step consists of at least one verification test by comparing computed results with closed-form analytical solutions to ensure that the various coupling effects are correctly realized. To demonstrate the robustness of the algorithm and to validate the UEL, a three-dimensional case study is performed with reference to the in-situ heating test of ATLAS at Belgium in 1980s. A hypothetical radionuclide leakage event is then simulated by activating the chemical-concentration degree of freedom and prescribing a constant high concentration at the heater's surface. The model predicts a limited contaminated regime after six years considering both diffusion and advection effects on species transport.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗