Search NASA⌕ Search

SEARCH · Search NASA

Results for “Parametric 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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Chemical looping air separation with Sr 0.8 Ca0.2Fe 0.9 Co 0.1 O 3-δ perovskite sorbent: Packed bed modeling, verification, and optimization

Chemical looping air separation (CLAS) represents a promising approach for efficient O 2 production from the air. This present study aims at optimizing the absorber/desorber operations and the separation process with extensive experimental validation. Specifically, a one-dimensional packed bed model was developed to investigate the CLAS operation with a Sr 0.8 Ca 0.2 Fe 0.9 Co 0.1 O 3-δ perovskite sorbent. The redox thermodynamics of perovskite sorbent was measured by TGA and then incorporated into a linear driving force model to describe the O 2 absorption and desorption rates. Both 4-step and 5-step air separation cycle configurations, with various cyclic structures, were performed in a subpilot-scale packed bed. The model predicted O2 purity and productivity were consistent with experimental results, supporting its accuracy and applicability. Parametric analysis and multi-objective optimization were further carried out to assess the performance of CLAS. Both O 2 purity and recovery increased monotonically with the cycle time, airflow rate, steam flow rate, and absorption pressure. Meanwhile, optimal O 2 productivity and power consumption can only be achieved by specific combinations of these parameters. The optimized results showed that CLAS can be highly competitive when compared to conventional pressure swing adsorption (PSA) or cryogenic distillation. The 5-step cycle configuration achieved a minimum power consumption of 118 kW·h for producing 1 ton O 2 with ≥ 95% purity. The maximum O 2 productivity reached 0.0932 g O2 /(g sorbent ·h) with 390 kW·h/ton O 2 of energy consumption (95% pure). The optimization results also indicate that CLAS can potentially be more efficient than cryogenic distillation even when the required O 2 purity is above 99%.

42 ENGINEERING↗
DOE OSTI2022

A microphysiological model of human trophoblast invasion during implantation

Abstract Successful establishment of pregnancy requires adhesion of an embryo to the endometrium and subsequent invasion into the maternal tissue. Abnormalities in this critical process of implantation and placentation lead to many pregnancy complications. Here we present a microenigneered system to model a complex sequence of orchestrated multicellular events that plays an essential role in early pregnancy. Our implantation-on-a-chip is capable of reconstructing the three-dimensional structural organization of the maternal-fetal interface to model the invasion of specialized fetal extravillous trophoblasts into the maternal uterus. Using primary human cells isolated from clinical specimens, we demonstrate in vivo-like directional migration of extravillous trophoblasts towards a microengineered maternal vessel and their interactions with the endothelium necessary for vascular remodeling. Through parametric variation of the cellular microenvironment and proteomic analysis of microengineered tissues, we show the important role of decidualized stromal cells as a regulator of extravillous trophoblast migration. Furthermore, our study reveals previously unknown effects of pre-implantation maternal immune cells on extravillous trophoblast invasion. This work represents a significant advance in our ability to model early human pregnancy, and may enable the development of advanced in vitro platforms for basic and clinical research of human reproduction.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

A multi‐variable framework for selecting WRF physics configurations at convection‐permitting scales: An Amazon wet‐season case study

Tropical convection over rainforests modulates atmospheric circulation and the energy and hydrological cycles across multiple scales. However, the scarcity of observations still limits our understanding of these processes. Although numerical models are utilized to investigate atmospheric physical processes, their performance depends on the choice of parameterizations and grid resolution. Here, in this work, we introduce and apply a multi‐variable, multi‐physics framework to evaluate and rank the Advanced Research Weather Research and Forecasting (WRF–ARW) configurations at 1‐km resolution over the central Amazon during the wet season. A 48‐member ensemble combines three land‐surface models (LSMs), four planetary boundary layer (PBL), and four microphysics (MP) schemes. Model performance for seven convection‐related near‐surface and boundary‐layer variables is assessed using Taylor diagrams and the Taylor skill score (TSS), analysis of variance (ANOVA)‐based sensitivity metrics, and non‐parametric rank tests. We then construct a combined, weighted TSS to identify configurations that are comparatively robust across variables. Results showed that most configurations reproduce near‐surface temperature, sensible and latent heat fluxes, and boundary‐layer height reasonably well, whereas humidity and rainfall remain challenging. LSM choice has the strongest impact on the surface fluxes and a secondary influence on near‐surface temperature and humidity, PBL schemes dominate boundary‐layer height, and MP schemes exert the largest control on rainfall. No single configuration is optimal for all variables, but the combination of Noah (LSM), Yonsei University (PBL), and Morrison (MP) emerges as the most robust configuration for this case, with the WRF single‐moment six‐class scheme (WSM6) providing a competitive, computationally cheaper MP alternative. The framework is general and can be applied to other regions, seasons, and convective regimes.

Amazon rainforest↗
DOE OSTI2020

Linear two-pool models are insufficient to infer soil organic matter decomposition temperature sensitivity from incubations

Abstract Terrestrial carbon (C)-climate feedbacks depend strongly on how soil organic matter (SOM) decomposition responds to temperature. This dependency is often represented in land models by the parameter Q 10 , which quantifies the relative increase of microbial soil respiration per 10 °C temperature increase. Many studies have conducted paired laboratory soil incubations and inferred “active” and “slow” pool Q 10 values by fitting linear two-pool models to measured respiration time series. Using a recently published incubation study (Qin et al. in Sci Adv 5(7):eaau1218, 2019) as an example, here we first show that the very high parametric equifinality of the linear two-pool models may render such incubation-based Q 10 estimates unreliable. In particular, we show that, accompanied by the uncertain initial active pool size, the slow pool Q 10 can span a very wide range, including values as high as 100, although all parameter combinations are producing almost equally good model fit with respect to the observations. This result is robust whether or not interactions between the active and slow pools are considered (typically these interactions are not considered when interpreting incubation data, but are part of the predictive soil carbon models). This very large parametric equifinality in the context of interpreting incubation data is consistent with the poor temporal extrapolation capability of linear multi-pool models identified in recent studies. Next, using a microbe-explicit SOM model (RESOM), we show that the inferred two pools and their associated parameters (e.g., Q 10 ) could be artificial constructs and are therefore unreliable concepts for integration into predictive models. We finally discuss uncertainties in applying linear two-pool (or more generally multiple-pool) models to estimate SOM decomposition parameters such as temperature sensitivities from laboratory incubations. We also propose new observations and model structures that could enable better process understanding and more robust predictive capabilities of soil carbon dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Thermodynamic modeling of countercurrent chemical looping reverse water gas shift process for redox material screening

The reverse water gas shift (RWGS) reaction is a key pathway for CO 2 utilization, particularly within Power-to-X process chains aimed at sustainable fuel and chemical production. Countercurrent chemical looping (CL-RWGS) using non-stoichiometric oxides can overcome equilibrium limitations of conventional RWGS reactors, enabling significantly higher CO 2 conversions. However, modeling the limiting performance of such systems is challenging due to their multiphase nature and coupled spatial and temporal variation in chemical composition. In this work, we present a discretized batch equilibrium model that simulates CL-RWGS reactors as a series of localized equilibrium exchanges between gas and solid elements. The model is numerically stable, computationally efficient, and free of kinetic source terms, making it well-suited for parametric studies and system-level integration. It is validated against established convection–diffusion models and shown to predict reasonable upper bounds on experimental results. Application of the model to a range of oxygen carrier materials identifies cerium–zirconium solid solutions, particularly Ce 0.80 Zr 0.20 O 2 , as a promising class offering superior oxygen storage characteristics compared to state-of-the-art La 0.6 Sr 0.4 FeO 3 . This framework provides a robust platform for materials screening, reactor sizing, and performance optimization in chemical looping systems. The model implementation is available as open-source software to support further research and development.

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

Evidential Deep Learning: Enhancing Predictive Uncertainty Estimation for Earth System Science Applications

Abstract Robust quantification of predictive uncertainty is a critical addition needed for machine learning applied to weather and climate problems to improve the understanding of what is driving prediction sensitivity. Ensembles of machine learning models provide predictive uncertainty estimates in a conceptually simple way but require multiple models for training and prediction, increasing computational cost and latency. Parametric deep learning can estimate uncertainty with one model by predicting the parameters of a probability distribution but does not account for epistemic uncertainty. Evidential deep learning, a technique that extends parametric deep learning to higher-order distributions, can account for both aleatoric and epistemic uncertainties with one model. This study compares the uncertainty derived from evidential neural networks to that obtained from ensembles. Through applications of the classification of winter precipitation type and regression of surface-layer fluxes, we show evidential deep learning models attaining predictive accuracy rivaling standard methods while robustly quantifying both sources of uncertainty. We evaluate the uncertainty in terms of how well the predictions are calibrated and how well the uncertainty correlates with prediction error. Analyses of uncertainty in the context of the inputs reveal sensitivities to underlying meteorological processes, facilitating interpretation of the models. The conceptual simplicity, interpretability, and computational efficiency of evidential neural networks make them highly extensible, offering a promising approach for reliable and practical uncertainty quantification in Earth system science modeling. To encourage broader adoption of evidential deep learning, we have developed a new Python package, Machine Integration and Learning for Earth Systems (MILES) group Generalized Uncertainty for Earth System Science (GUESS) (MILES-GUESS) ( https://github.com/ai2es/miles-guess ), that enables users to train and evaluate both evidential and ensemble deep learning. Significance Statement This study demonstrates a new technique, evidential deep learning, for robust and computationally efficient uncertainty quantification in modeling the Earth system. The method integrates probabilistic principles into deep neural networks, enabling the estimation of both aleatoric uncertainty from noisy data and epistemic uncertainty from model limitations using a single model. Our analyses reveal how decomposing these uncertainties provides valuable insights into reliability, accuracy, and model shortcomings. We show that the approach can rival standard methods in classification and regression tasks within atmospheric science while offering practical advantages such as computational efficiency. With further advances, evidential networks have the potential to enhance risk assessment and decision-making across meteorology by improving uncertainty quantification, a longstanding challenge. This work establishes a strong foundation and motivation for the broader adoption of evidential learning, where properly quantifying uncertainties is critical yet lacking.

Schreck, John S.↗
DOE OSTI2022

Predicted Impacts of Pt and Ionomer Distributions on Low-Pt-Loaded PEMFC Performance

Low-cost, high performance proton exchange membrane fuel cells (PEMFCs) have been difficult to develop due to limited understanding of coupled processes in the cathode catalyst layer (CCL). Low-Pt-loaded PEMFCs suffer losses beyond those predicted solely due to reduced catalyst area. Although consensus links these losses to thin ionomer films in the CCL, a precise mechanistic explanation remains elusive. In this publication, we present a physically based PEMFC model with novel structure-property relationships for thin-film Nafion, validated against PEMFC data with low Pt loading. Results suggest that flooding exacerbates kinetic limitations in low-loaded PEMFCs, shifting the Faradaic current distribution. As current density increases, protons travel further into the CCL, resulting in higher Ohmic overpotentials. We also present a parametric study of CCL design parameters. We find that graded Pt and ionomer loadings reduce Ohmic losses and flooding, but individually do not provide significant improvements. However, a dual-graded CCL (i.e., graded Pt and ionomer) is predicted to significantly improve the maximum power density and limiting current compared to uniformly loaded CCLs. This work highlights the importance of accurate transport parameters for thin-film Nafion and provides a pathway to low-cost PEMFCs via precise control of CCL microstructures.

08 HYDROGEN↗
DOE OSTI2025

Coupled geomechanical investigation of depletion-induced fault reactivation

Fault reactivation during subsurface fluid production pose significant challenges to safe and sustainable resource extraction. Here, this study presents a three-dimensional coupled geomechanical framework to investigate the processes driving fault reactivation, capturing the interactions between reservoir dynamics and geomechanical responses. Verification against theoretical estimations based on linear poroelasticity confirms the model's capacity in representing reservoir background stress responses. However, the study reveals that relying solely on background stress states can underestimate or overestimate fault reactivation potential, emphasizing the importance of including localized stress perturbations such as differential compaction and stress redistribution. Applied to a fault (M1) inspired by the geological characteristics of the Groningen field, the model shows slip initiation at 2965 m depth with 16.0 MPa depletion, aligning with field observations where seismicity occurred at approximately 3 km depth after 15.8 MPa depletion. Parametric studies reveal: (1) inelastic reservoir compaction delays fault reactivation and mitigates fault slip by reducing stress concentration, (2) higher intermediate in-situ stress magnitudes decrease the Coulomb Failure Stress (CFS) increase rate and reduce fault slip, (3) larger fault offsets amplify shear stress near the offset zone, promoting earlier reactivation and longer rupture propagation, and (4) fault permeability significantly influences pressure diffusion, with low-permeability faults leading to sharper stress changes and earlier fault destabilization. These insights highlight the critical role of geological and mechanical parameters in fault reactivation and provide a predictive framework for mitigating induced seismicity risks.

Coupled geomechanical modeling↗
DOE OSTI2023

Laser additive manufacturing for infrastructure repair: A case study of a deteriorated steel bridge beam

As our transportation infrastructure ages, its deterioration is becoming a growing concern. Rapid and effective methods for repair can reduce the need for time-and cost-intensive replacement of structures such as roadway bridges. In this work, laser-based additive manufacturing (AM) was developed for the re-pair of damaged beams in steel bridges. A corroded structural beam of low-carbon A36 steel was selected for this case study and 316L stainless steel was used for repair due to its excellent mechanical properties and corrosion resistance. Simulated repair experiments were performed on specimens extracted from the beam, and prepared with various groove shapes (Rectangular, Trapezoidal, and U-shaped). More in-depth parametric and microstructural studies were carried out with a series of U-shaped specimens repaired with the linear laser input energy varied from 19.5 to 30.0 J mm -1 . The tensile properties of the repaired specimens and interface structures were also investigated in detail. It was found that the repaired specimens had a slightly higher tensile strength despite a reduction in tensile elongation compared to the A36 steel base metal. The influence of interfacial microstructure and heterogeneity on the mechanical proper-ties of the repaired specimens was discussed. Furthermore, our work suggests the promising potential to employ AM for structural repair and provides fundamental insights into processing-structure-property relationships in laser AM-repaired materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Generative adversarial networks for scintillation signal simulation in EXO-200

Generative Adversarial Networks trained on samples of simulated or actual events have been proposed as a way of generating large simulated datasets at a reduced computational cost. In this work, a novel approach to perform the simulation of photodetector signals from the time projection chamber of the EXO-200 experiment is demonstrated. The method is based on a Wasserstein Generative Adversarial Network — a deep learning technique allowing for implicit non-parametric estimation of the population distribution for a given set of objects. Our network is trained on real calibration data using raw scintillation waveforms as input. We find that it is able to produce high-quality simulated waveforms an order of magnitude faster than the traditional simulation approach and, importantly, generalize from the training sample and discern salient high-level features of the data. In particular, the network correctly deduces position dependency of scintillation light response in the detector and correctly recognizes dead photodetector channels. Furthermore, the network output is then integrated into the EXO-200 analysis framework to show that the standard EXO-200 reconstruction routine processes the simulated waveforms to produce energy distributions comparable to that of real waveforms. Finally, the remaining discrepancies and potential ways to improve the approach further are highlighted.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Systems for Transducing Entanglement Between Photons and Electron Spins

The transfer of quantum properties and entanglement between photons and matter (transduction) is a major challenge in QIS, especially with regard to networking quantum information processing. In order to overcome this challenge, it is important to know how non-classical light (single or entangled photons) interacts with matter. This project specifically addresses fundamental questions about the interaction of single and entangled photons with molecules and nanostructures to produce electrons with prescribed spin polarizations and pairs of electrons with entangled spins, thus resulting in transduction by addressing two major objectives: (2) In Objective 1, we will investigate the formation of polarization-encoded radical (spin qubit) pairs and molecular triplet states photoexcited by single quantum emitter light sources. (2) In Objective 2, we will demonstrate the use of two entangled photons, generated by spontaneous parametric down conversion (SPDC) or by the coherent biexcitonic states of perovskite NPLs, to produce two entangled radical (spin qubit) pairs or molecular triplet states.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

A Fast and Accurate Reduced-Order Model for High-Intensity Transferred Arc Discharges

Arc discharges are widely used in welding, plasma smelting, and other industrial processes, where variations in operating conditions strongly affect arc stability, temperature distribution, and energy transfer. Accurate modeling of these phenomena typically requires computationally expensive high-fidelity simulations. This study presents a hierarchy of three arc discharge models with progressively reduced physical fidelity. The high-fidelity model provides a fully physics-resolved reference, the reduced-order model reproduces the dominant thermal and flow characteristics of the arc with reasonable accuracy, and the Elenbaas-Heller model captures key trends in a simplified, rapid formulation. This hierarchy demonstrates that reduced-order model can effectively balance predictive fidelity and computational efficiency, providing practical tools for arc simulation and parametric studies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Toward Robust and Routine Determination of Mw for Small Earthquakes: Application to the 2020 Mw 5.7 Magna, Utah, Seismic Sequence

To better characterize seismic hazard, particularly, for induced seismicity, there is an increasing interest in methods to estimate moment magnitude (M w ) for small earthquakes. M w is generally preferred over other magnitude types, but, it is difficult to estimate M w for earthquakes with local magnitude (M L ) <3-3.5, using conventional moment tensor (MT) inversion. The 2020 M ww 5.7 Magna, Utah, seismic sequence provides an opportunity to illustrate and evaluate the value of spectral methods for this purpose. Starting with a high-quality seismic catalog of 2103 earthquakes (M L <5.6), we estimate M w using two independent spectral methods—one based on direct waves, yielding M w,direct , and the other based on coda waves, yielding M w,coda . For the direct-wave method, we present a non-parametric (NP) inversion scheme that solves for apparent geometrical spreading, G(R), and site effects (S), similar to other NP procedures that have been used to calibrate regional M L scales. The NP inversion is constrained using M ws derived from MTs for nine events in the Magna sequence. We recover statistically robust and physically reasonable G(R) and S and compute M w,direct for 635 Magna earthquakes down to M L 0.7. For the coda-wave method, we consider two separate calibration schemes involving previous MT solutions and compute M w,coda for 311 earthquakes down to M L 1.0. For 280 of the events that were processed with both methods—M w,direct and M w,coda —are strongly correlated (r = 0.98), with a mean difference of only 0.05. We compare M w,direct and M w,coda with M L and find reasonably good agreement for M L <3.6 with the theoretically predicted relationship of M w =(2/3)M L +C, in which C is a regional constant. Our results imply that seismic network operators can use spectral-based M w estimates to replace M L estimates for events with M L ≥1.0, and possibly smaller. The main requirement is the existence of a small number of MT solutions for calibration purposes.

58 GEOSCIENCES↗
DOE OSTI2020

Modeling occupancy-driven building loads for large and diversified building stocks through the use of parametric schedules

Building energy modeling provides a fundamental tool to assess the potential for energy efficiency to contribute to reducing world energy consumption and global emissions. Occupancy-related operations are a key source of uncertainty for building energy analysis, particularly for aggregated building stocks. At a district or city level, it is critical to estimate aggregated power load profiles for sizing power grid infrastructure, power plant capacity allocation, and energy efficiency measures. The stochastic nature of behavior-related operations complicates the creation of models that accurately capture building load profiles for entire building stocks. This research introduces a new methodology called parametric schedules to model occupancy-driven schedules for large and diverse building stocks. In contrast to computationally expensive methodologies proposed in the literature, our work does not use a recursive time-consuming step. Occupancy is estimated by the extrapolation of operation times directly from metered electric consumption data; occupancy-related schedules are stochastically assigned to each building model, guaranteeing diversity of operation times in the stock. Our procedure has been tested on a large, diversified data-set of 25,000 commercial buildings in Los Angeles, California. It proved to be able to adequately represent the stochastic schedules diversity of the stock and to refine the stock calibration process by 1%. This innovative approach represents a useful asset for utility companies, grid operators, urban planners, and balancing authorities, seeking to improve building stock modeling and better estimate the impact of energy conservation measures. – This work is part of a larger stock modeling tool called ComStock, which is under development by NREL.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Reduction of the neutron imaginary potential off the stability line and its possible impact on neutron capture rates

The effect of the isovector imaginary optical potential has been studied experimentally by using the particle-evaporation technique for the 11 B + 48 Ca reaction with a 21.8 MeV 11 B beam. Spectra of neutron, proton, and α particles emitted from the neutron-rich compound nucleus 59 Mn have been measured and analyzed with traditional optical model potentials with their original parametrizations as well as with adjusted isovector imaginary components. It is shown that the isovector component of the imaginary potential is indispensable in the reproduction of proton and α-particle yields from this reaction and even needs to be enhanced compared with the suggestions of the original model parametrizations. Furthermore, this can lead to important consequences for astrophysical reaction-rate calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

A Conceptual Design of the Reactor Cavity Cooling System for the Horizontal Compact High Temperature Gas Reactor (HC-HTGR)

Horizontal Compact High Temperature Gas Reactor (HC-HTGR) is being designed by a multi-disciplinary team of nuclear, mechanical, and structural engineers under the support of a DOE-NE Advanced Reactor Demonstration Program’s Advanced Reactor Concepts-20 (ARC-20) award. The objective of this ARC-20 project is to deliver a conceptual design for the proposed HC-HTGR in 3 years and support its commercialization as a safe, low-cost HTGR. Argonne National Laboratory (Argonne) is responsible for the design and analysis of the reactor cavity cooling system (RCCS) as a safety system for passive decay heat removal of the reactor concept. This report documents the design study to derive a conceptual design study of the RCCS for the HC-HTGR. It includes the identification of the functions and requirements of the HC-HTGR RCCS, design analyses including high-level design consideration and the calculations for optimizing design space of the system with supporting component-level analysis to inform the material selection and performance of the water panel, the description of the conceptual design of the HC-HTGR RCCS derived based on the analyses results, and performance evaluation of the conceptual RCCS for the HC-HTGR. A detailed concept of the RCCS has been identified and high-level system requirements has been developed for the HC-HTGR. Design space focusing on the natural circulation loop portion of the RCCS has been investigated to optimize the system performance. The initial baseline dimensions were firstly derived based on the scoping calculations. A component level design analysis was conducted for the water panel to inform the material selection and to assess its conduction performance. A preliminary system-level performance analysis was performed for the 1/8th of the compartment of the initial baseline design of the RCCS using RELAP5-3D. To improve the system thermal performance, the RCCS design has been updated by exploring various design options by design parametric analyses. Based on the results, the conceptual design of the RCCS for the HC-HTGR has been derived, which satisfies the target performance of ~1 MWt at the elevated vessel wall temperature conditions. Transient simulations were conducted for the conceptual RCCS design for the HC-HTGR under various operation modes and heat load conditions using RELAP5-3D. The system dynamics in different operating states was investigated and the system performance under transients of interest was evaluated. The results demonstrated the overall system feasibility that the RCCS design maintains structures temperatures lower than maximum allowable temperature with sufficient system inventory without any active heat removal in the design process with certain transients addressed. The HC-HTGR RCCS will have additional design updates of subsystems or optimization of the system components during the preliminary and final design phases. Since the entire plant has not been integrated yet, this delivered conceptual design is subject to changes for integration, that require additional conceptual design activities and Quality and Assurance implementation (Q&A). The performance assessment of the RCCS for the HC-HTGR will be then revisited and optimized to finalize the system design, and the RCCS integrated primary system analysis will be utilized to simulate selective accident scenarios of interest where efforts are currently undergoing in the project.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Modeling of fluid injection-induced fault reactivation using coupled fluid flow and mechanical interface model

The present study is aimed at developing a numerical model to reproduce coupled hydro-mechanical processes associated with fault reactivation by fluid injection in low permeability rock, as part of the DECOVALEX-2019 project Task B. We proposed a modeling approach for simulating the processes using the TOUGH-FLAC simulator, and modeled a fault reactivation experiment conducted at Mont Terri Rock Laboratory in Switzerland. The first step of the study involved benchmark calculations considering a simplified fault plane and geometry. Fluid flow along a fault was modeled using elements of aperture-sized thickness on the basis of Darcy's law and the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by zero-thickness interface elements in FLAC3D upon which a slip and/or separation is allowed. A methodology to connect a TOUGH2 volume element to a FLAC3D interface element was developed for handling the hydro-mechanical interactions on the fault during fluid injection. Two different fault models for describing the evolutions of hydraulic aperture by elastic fracture opening and failure-induced aperture increase were considered in the benchmark calculations. In the coupling process, the changes in geometrical features and hydrological properties induced by mechanical deformation were continuously updated. The transient responses of the fault and host rock to stepwise pressurization were examined during the simulation. The hydro-mechanical behavior, including the injection flow rate, pressure distribution around the borehole, stress conditions, and displacements in normal and shear directions were monitored in the surrounding rock and along the fault. The results of benchmark calculations suggest that the developed model reasonably represents the hydro-mechanical behavior of a fault and the surrounding rock. This modeling approach was applied to the fault reactivation experiment of the Mont Terri Rock Laboratory. In this interpretive modeling, a parametric study was conducted to examine the effects of input parameters regarding in situ stress and fault properties on the hydro-mechanical responses of the fault to water injection. Then, an optimal parameter set to reproduce the field experiment results was chosen by trial-and-error. The injection flow rate and pressure response during fault reactivation closely matched those obtained at the site, which indicates the capability of the model to appropriately capture the progressive pathway evolution during fault reactivation tests at the site. Finally, the anchor displacements were overestimated by the model, but a fair agreement was obtained in terms of the order of magnitude and the variation tendency.

58 GEOSCIENCES↗
DOE OSTI2020

Quantifying Uncertainties In Neutron-Alpha Scattering With Chiral Nucleon-Nucleon And Three-Nucleon Forces

We report that modern ab initio theory combined with high-quality nucleon-nucleon (NN) and three-nucleon (3N) interactions from chiral effective field theory (EFT) can provide a predictive description of low-energy light-nuclei reactions relevant for astrophysics and fusion-energy applications. However, the high cost of computations has so far impeded a complete analysis of the uncertainty budget of such calculations. Starting from NN potentials up to fifth order (N 4 LO) combined with leading-order 3N forces, we study how the order-by-order convergence of the chiral expansion and confidence intervals for the 3N contact and contact-plus-one-pion-exchange low-energy constants (c E and c D ) contribute to the overall uncertainty budget of many-body calculations of neutron- 4 He (n-α) elastic scattering. We compute structure and reaction observables for three-, four-, and five-nucleon systems within the ab initio frameworks of the no-core shell model and no-core shell model with continuum. Using a small set of design runs, we construct a Gaussian process model (GPM) that acts as a statistical emulator for the theory. With this, we gain insight into how uncertainties in the 3N low-energy constants propagate throughout the calculation and determine the Bayesian posterior distribution of these parameters with Markov-Chain Monte Carlo. We find rapidly converging n-α phase shifts with respect to the chiral order. With the adopted leading-order 3N force, calculations based on the NN interaction at N 4 LO of Entem, Machleidt, and Nosyk are unable to reproduce the experimental phase shifts in the 3/2 - channel within the estimated chiral truncation errors. Closer agreement with empirical data is found when using an older parametrization of the NN interaction at order N 3 LO, and the position and width of the P-wave resonances can be used to reduce the uncertainty of the 3N low-energy constants. The present results point to a lack of spin-orbit strength when the newer parametrization of the chiral NN force up to fifth order is combined with the leading-order 3N force. The inclusion of higher-order 3N-force terms may be required to recover the missing strength. GPMs can act as fast and accurate emulators of ab initio many-body calculations of low-energy scattering and reactions of light nuclei, opening the way to a robust quantification of theoretical uncertainties grounded in the description of the underlying chiral Hamiltonian.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗