Numerical manifold method modeling of coupled processes in fractured geological media at multiple scales
Not Available
58 GEOSCIENCES↗
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.
Coupling processes related to the Sun-weather problem
Physical mechanisms for coupling the energetics of solar activity to meteorological responses are reviewed. Although several hypotheses have been advanced, none can be said to be sufficiently complete to be applied to weather or climate prediction. Solar activity indicators potentially useful for forecasting are identified, including sunspots, solar flares, and magnetic sector boundary crossings. Additional experiments, studies, and analyses are required before Sun-weather concepts can be utilized for predicting meteorological responses.
Goldberg, R. A.↗
International collaboration in disposal research: comparative modeling of coupled processes in the DECOVALEX project
Abstract. This presentation gives an overview of an international research collaboration for advancing the understanding and modeling of coupled thermo-hydro-mechanical-chemical (THMC) processes in geological systems. The creation of the international DECOVALEX project, now running for more than 25 years, was motivated by the recognition that prediction of these coupled effects is an essential part of the performance and safety assessment of geologic disposal systems for radioactive waste and spent nuclear fuel. DECOVALEX emphasizes joint analysis and comparative modeling of state-of-the-art field and laboratory experiments, across a range of host rock options and repository designs. Participating research teams are from radioactive waste management organizations, national research institutes, regulatory agencies, universities, and consulting groups, providing a wide range of perspectives and solutions to these complex problems. The presentation provides examples of the research contributions made collectively in past DECOVALEX tasks and also touches on the unique modeling challenges tackled in the ongoing project phase, referred to as DECOVALEX-2023. The current phase comprises 17 partner organizations, about 50 modeling teams, and 7 modeling tasks, which cover a broad portfolio from fundamental studies on gas migration to full-scale in situ heater experiments in different host rocks to performance assessment studies. Together, these examples illustrate that the insight and scientific knowledge gained within the DECOVALEX project would not have been possible if one group had studied these complex THMC modeling challenges alone rather than within a truly collaborative setting.
Birkholzer, Jens T.↗
Bentonite buffer under high temperature: laboratory experiments and coupled process modeling
Abstract. Bentonite buffer in a geological repository will be simultaneously heated from decaying radioactive waste and hydrated from the surrounding host rock, triggering complex and coupled THMC (thermal–hydrological–mechanical–chemical) processes. Understanding the THMC behavior of bentonite-based engineered barrier system (EBS) is key to the evaluation and prediction of its long-term performance. Studies on the THMC process have been focused on conditions under 100 ∘C, as most design concepts impose a thermal limit of 100 ∘C in bentonite. Recently, studies under high-temperature conditions have been conducted to evaluate the possibility of raising the thermal limit and expanding the data/knowledge base to increase the confidence level. In this abstract, we present a series of bench-scale laboratory experiments at high temperatures (up to 200 ∘C) and the corresponding modeling work. Two sets of column tests were conducted, and each set consisted of two test columns: a control column undergoing only hydration (non-heated) and an experiment column experiencing both heating and hydration (heated). During the experiment, frequent X-ray computed tomography (CT) images were collected to provide a 3D visualization of the density distribution and present the spatiotemporal evolution of (1) hydration/dehydration, (2) clay swelling/shrinkage, (3) displacement, and (4) mineral precipitation. The two sets of tests differ with respect to several experimental conditions, such as bentonite type, compacted density and water content, water chemistry, and hydration pressure, but the important difference is that the first set used bentonite powder with a dry density of 1.28 g cm−3, whereas the second set used granulated bentonite (mixture of pellets and powder) with a dry density of 1.45–1.5 g cm−3. In both sets of experiments, a comprehensive post-dismantling characterization of bentonite samples was carried out after the column tests had been running for 1.5 years. Comparing non-heated and heated columns, the temperature gradient led to lower degree of homogenization of bentonite after bentonite became fully saturated; comparing the first and second sets, granulated and powdered bentonite exhibited drastically different hydration behavior. A THM model with a 2D axisymmetric grid system was used to interpret the data from the first set of tests. The model considers the combined impact of saturation, fluid pressure, and porosity change due to swelling/compression on the spatiotemporal distribution of bulk density and movement of the thermocouple modules. Observations from the tests help us understand the early perturbation of bentonite buffer under high temperature, and data from these tests improve the calibration of key constitutive hydrological and mechanical models and, therefore, enhance the modeling capability with respect to calculating the long-term evolution of bentonite buffer.
Zheng, Liange↗
The role of nonlocalities in magnetosphere-ionosphere coupling processes
Microinstabilities are believed to play a crucial role in the physics of magnetosphere-ionosphere coupling. The current driven ion cyclotron instability is a very important microinstability in this respect. A nonlocal formalism is given for studying the ion cyclotron instability in a more realistic magnetospheric environment than is available in the widely used local theory. This formalism includes the magnetic shear produced self-consistently by the field aligned currents and the finite extent of such currents. Significant departures from the local theory are noted.
Ganguli, G.↗
Investigation of Coupled Processes in Fractures and the Bordering Matrix via a Micro‐Continuum Reactive Transport Model
Abstract In multi‐mineral fractured rocks, the altered porous layer on the fracture surface resulting from preferential dissolution of the fast‐reacting minerals can have profound impacts on subsequent chemical‐physical alteration of the fractures. This study adopts the micro‐continuum approach to provide further understanding of reactive transport processes in the altered layer (AL), and mass exchanges with the bordering matrix and fracture. The modeling framework couples the Darcy‐Brinkman‐Stokes (DBS) solver in COMSOL Multiphysics and the geochemical modeling capability of CrunchFlow. Three‐dimensional steady state simulations with systematically varied chemical‐physical parameters of the AL were performed to examine the impacts of individual factors and processes. Our simulation results confirm previous observations that dissolution of the fast‐reacting mineral (i.e., calcite) is largely controlled by diffusion across the AL. We also show that dissolution of the slow‐reacting mineral (i.e., dolomite), which controls AL development and fracture enlargement, increases with surface area and has a complex dependence on different local rate‐limiting processes. In particular, advection can result in evident spatial variations in the local dissolution rates of dolomite, although it does not affect the bulk chemistry significantly. The difference in the spatial patterns between simulations with and without advection in the AL is more noticeable in the locations with smaller apertures, with up to 20% difference in local reaction rates. Therefore, it is important to include a full depiction of advection, diffusion, and reactions for accurately capturing local dynamics that control long‐term fracture evolution.
58 GEOSCIENCES↗
Balancing sugar recovery and inhibitor generation during energycane processing: Coupling cryogenic grinding with hydrothermal pretreatment at low temperatures
Pretreatment of lignocellulosic biomass at high temperatures or with oxidizing chemicals generate various inhibitors that restrict the efficient bioconversion of sugars in subsequent steps. The present study systematically investigates individual and combinatorial effects of pretreatment parameters on the generation of inhibitors. Herein, a plot between pretreatment temperature and inhibitor revealed optimum pretreatment temperature for energycane bagasse i.e., 170 °C beyond which total inhibitor production increased exponentially. No inhibitor production was observed on mechanical processing i.e., disk milling/cryogenic grinding of biomass. Evaluation of response surface regression exhibited that biomass solids loading has a significant effect on inhibitor generation at higher temperatures. The concentrations of certain inhibitors such as acetic acid, furfurals, and HMF increased more than 3-folds on doubling the solids loading. Furthermore, a novel low-severity approach of low-temperature hydrothermal pretreatment coupled with cryogenic grinding for lignocellulosic biomasses has been introduced which improved sugar yields while maintaining a low inhibitor concentration.
09 BIOMASS FUELS↗
Multigrid reduction preconditioning framework for coupled processes in porous and fractured media
Not Available
02 PETROLEUM↗
Approaches for the Simulation of Coupled Processes in Evolving Fractured Porous Media Enabled by Exascale Computing
Models have historically represented fractured porous media with continuum descriptions that characterize the media using bulk parameters. The impact of small-scale features is not captured in these models, although they may be controlling the performance of subsurface applications. Pore-scale models can simulate processes in small-scale features by representing the pore space geometry explicitly but are computationally expensive for large domains. The alternative multiscale approach entails the combination of pore-scale and continuum-scale descriptions in a single framework. We use Chombo-Crunch, a computational capability that discretizes complex geometries with an adaptive, embedded boundary method to contrast these two approaches. Chombo-Crunch takes advantage of recent computational performance and memory bandwidth improvements resulting from the emergence of exascale computing resources. These combined improvements enable the efficient simulation of reactive transport in fractured media with a high degree of fidelity and the ability to capture the control small-scale processes exert on the overall medium evolution.
42 ENGINEERING↗
A high-order language for a system of closely coupled processing elements
The research reported in this paper was occasioned by the requirements on part of the Real-Time Digital Simulator (RTDS) project under way at NASA Lewis Research Center. The RTDS simulation scheme employs a network of CPUs running lock-step cycles in the parallel computations of jet airplane simulations. Their need for a high order language (HOL) that would allow non-experts to write simulation applications and that could be implemented on a possibly varying network can best be fulfilled by using the programming language Ada. We describe how the simulation problems can be modeled in Ada, how to map a single, multi-processing Ada program into code for individual processors, regardless of network reconfiguration, and why some Ada language features are particulary well-suited to network simulations.
Feyock, S.↗
DE-1 observations of solar wind-magnetosphere coupling processes in the polar cusp
Measurements of hot plasma distribution functions in the midaltitude polar cusp by the Dynamics Explorer-1 spacecraft have provided new information on the processes of magnetosheath plasma injection and transport and the dependence of these phenomena on solar-wind parameters. Results on magnetosheath plasma injection have included: (1) localization of the region of plasma injection to the high-latitude cusp-magnetosheath interface, (2) movement of this region toward lower latitudes as B(z) becomes more southward, (3) maintenance of quasi-neutrality of suprathermal plasmas within the cusp, and (4) support for merging-type injection at the high-latitude cusp boundary for northward IMF. Plasma flow measurements from Dynamics Explorer-1 have led to a global IMF-dependent convection model which contains both merging and viscous driving forces.
Burch, J. L.↗
REVEALING THE COUPLING PROCESS BETWEEN AEROSOL, PBL, AND CLOUD: IDENTIFICATION AND MECHANISMS
Explore the source record for details and available documents.
Geosciences↗
Understanding Uncertainty in Microstructure Evolution and Constitutive Properties in Additive Process Modeling
Coupled process–microstructure–property modeling, and understanding the sources of uncertainty and their propagation toward error in part property prediction, are key steps toward full utilization of additive manufacturing (AM) for predictable quality part development. The OpenFOAM model for process conditions, the ExaCA model for as-solidified grain structure, and the ExaConstit model for constitutive mechanical properties are used as part of the ExaAM modeling framework to examine a few of the various sources of uncertainty in the modeling workflow. In addition to “random” uncertainty (due to random number generation in the orientations and locations of grains present), the heterogeneous nucleation density N 0 and the mean substrate grain spacing S 0 are varied to examine their impact of grain area development as a function of build height in the simulated microstructure. While mean grain area after 1 mm of build is found to be sensitive to N 0 and S 0 , particularly at small N 0 and large S 0 (despite some convergence toward similar values), the resulting grain shapes and overall textures develop in a reasonably similar manner. As a result of these similar textures, ExaConstit simulation using ExaCA representative volume elements (RVEs) from various permutations of N 0 , S 0 , and location within the build resulted in similar yield stress, stress–strain curve shape, and stress triaxiality distributions. It is concluded that for this particular material and scan pattern, 15 layers is sufficient for ExaCA texture and ExaConstit predicted properties to become relatively independent of additional layer simulation, provided that reasonable estimates for N 0 and S 0 are used. However, additional layers of ExaCA will need to be run to obtain mean grain areas independent of build height and baseplate structure.
36 MATERIALS SCIENCE↗
A Physics-Constrained Deep Learning Model for Simulating Multiphase Flow in 3D Heterogeneous Porous Media
Physics-based simulators for multiphase flow in porous media emulate nonlinear processes with coupled physics, and usually require extensive computational resources for software development, maintenance and simulation execution. As a result, a huge demand exists for fast modeling of coupled processes in a wide range of subsurface applications including geological sequestration, hydrocarbon recovery and geothermal energy extraction. In this work, an efficient physics-constrained deep learning model is developed for solving multiphase flow in 3-Dimensional (3D) heterogeneous porous media. The model fully leverages the spatial topology predictive capability of convolutional neural networks, specifically U-Net with successive contracting and expansive steps, and is coupled with an efficient continuity-based smoother to predict flow responses that need spatial continuity. Furthermore, the transient regions are penalized to steer the training process such that the model can accurately capture flow in these regions. The model takes inputs including properties of porous media, fluid properties and well controls, and predicts the temporal-spatial evolution of the state variables (pressure and saturation). While maintaining the continuity of fluid flow, the 3D spatial domain is decomposed into 2D images for reducing training cost, and the decomposition results in an increased number of training data samples and better training efficiency. Additionally, a surrogate model is separately constructed as a postprocessor to calculate well flow rate based on the predictions of state variables from the deep learning model. We use the example of CO 2 injection into saline aquifers, and apply the physics-constrained deep learning model that is trained from physics-based simulation data and emulates the physics process. The model performs prediction with a speedup of ~ 1400 times compared to physics-based simulations, and the average temporal errors of predicted pressure and saturation plumes are 0.27% and 0.099% respectively. Furthermore, water production rate is efficiently predicted by a surrogate model for well flow rate, with a mean error less than 5%. Therefore, with its unique scheme to cope with the fidelity in fluid flow in porous media, the physics-constrained deep learning model can become an efficient predictive model for computationally demanding inverse problems or other coupled processes.
58 GEOSCIENCES↗
Integrated Process Testing of MSR Salt Spill Accidents
Part of the licensing process for new nuclear reactors requires vendors to assess the potential consequences of identified accident scenarios using accident progression modeling. The accident scenario that will likely be evaluated by all molten salt reactor (MSR) developers is a spill of radionuclide-bearing fuel salt onto the reactor containment floor (i.e., a salt spill accident). The development of accident progression models requires experimental data to inform which processes to incorporate, to enable the calculation of parameters to model these processes, and to validate the model predictions. The data should quantify the sensitivities of key processes (e.g., molten salt spreading, heat transfer, containment structure corrosion, radionuclide vaporization, and aerosol generation) towards the initial conditions of the spill, the ambient environment, and the features of the containment. In addition, results from integrated process tests that quantify coupled processes are required to validate systems-level models. This report documents results from integrated process tests conducted on simulated molten salt spill accidents. The generated experimental data simultaneously quantify the heat transfer behavior of the spilled salt, compositional changes to the bulk salt, and the release of surrogate fission products from the spilled salt as aerosol particles. All tests that were conducted used FLiNaK doped with surrogate fission products, and the variables that were evaluated included the initial salt temperature and the concentration of surrogate fission products present in the salt. The major accomplishments of this work include identifying surrogate fuel salt compositions that provide insight into the dispersal behavior of radionuclides of potential significance to the source term, employing previously developed methods and measurement techniques to simultaneously measure key processes, generating data on the coupled processes of molten salt heat transfer and surrogate fission product release as aerosol particles, demonstrating new test methods for real-time monitoring of the flow rate of the spill and aerosol size quantification in an argon atmosphere, and developing a mass transfer model for cesium and iodine release from molten FLiNaK to provide insight into aerosol formation by vapor condensation. The same methodology applied herein can be employed to study different salt compositions of interest to MSR developers, different environmental conditions, and other variables that are relevant to postulated accident scenarios. The insights gained from these integrated process tests conducted at a laboratory scale will be incorporated into future integral effects tests conducted at an engineering scale.
20 FOSSIL-FUELED POWER PLANTS↗
A process-based evaluation of biases in extratropical stratosphere–troposphere coupling in subseasonal forecast systems
Abstract. Two-way coupling between the stratosphere and troposphere is recognized as an important source of subseasonal-to-seasonal (S2S) predictability and can open windows of opportunity for improved forecasts. Model biases can, however, lead to a poor representation of such coupling processes; drifts in a model's circulation related to model biases, resolution, and parameterizations have the potential to feed back on the circulation and affect stratosphere–troposphere coupling. We introduce a set of diagnostics using readily available data that can be used to reveal these biases and then apply these diagnostics to 22 S2S forecast systems. In the Northern Hemisphere, nearly all S2S forecast systems underestimate the strength of the observed upward coupling from the troposphere to the stratosphere, downward coupling within the stratosphere, and the persistence of lower-stratospheric temperature anomalies. While downward coupling from the lower stratosphere to the near surface is well represented in the multi-model ensemble mean, there is substantial intermodel spread likely related to how well each model represents tropospheric stationary waves. In the Southern Hemisphere, the stratospheric vortex is oversensitive to upward-propagating wave flux in the forecast systems. Forecast systems generally overestimate the strength of downward coupling from the lower stratosphere to the troposphere, even as most underestimate the radiative persistence in the lower stratosphere. In both hemispheres, models with higher lids and a better representation of tropospheric quasi-stationary waves generally perform better at simulating these coupling processes.
Garfinkel, Chaim I. (ORCID:000000017258666X)↗
Measurement Acceptance Region (MAR) Assessment Results Based On Sludge Batch 10 Projections from December 2021, February 2022 and August 2022
The Defense Waste Processing Facility (DWPF) is currently preparing to initiate processing of Sludge Batch 10 (SB10), which is comprised of material from Tanks 11H, 13H, 15H, and 26F, Alternate Feed Stock-2 and Sodium Reactor Experiment material from H-Canyon. Frit 473 (8B 2 O 3 -8Li 2 O-5Na 2 O-79SiO 2 , in weight percent) was recommended for sludge-only (SO) and coupled processing with the Salt Waste Processing Facility (SWPF) based on previous assessments of SB10 projections with the DWPF Product Composition Control System (PCCS) glass property models and their associated Measurement Acceptance Region (MAR) constraints. Due to the lower processing rate of Sludge Batch 9 (SB9), the heel in Tank 40 is anticipated to be greater than 40 inches at the projected start of SB10 processing to meet the upcoming Accelerated Basin Deinventory addition dates in Tank 51 for Sludge Batch 11. In December 2021 and February 2022, Savannah River Remediation System Planning provided updated SB10 Tank 40 blend projections based on heels of 74 inches, 84 inches, 94 inches and 103.4 inches. Savannah River Mission Completion (SRMC) subsequently pursued Wash Cycle Y to further reduce the total sulfur in the sludge batch and increase processing flexibility at DWPF. In August 2022, SRMC System Planning provided an updated SB10 Tank 40 blend projection with a 76 inch-heel representing Decant Y2 and 60 kilogallons of bearing water inleakage that is anticipated during the Tank 51 to Tank 40 transfer. The objectives of this task were to: • Determine the impact on the operating windows for SO and coupled processing • Determine whether any composition gaps exist between the already completed SB10 variability study and the reprojected SB10 glass composition region • Compare the SB10 reprojected glass composition region to the DWPF PCCS model development and validation ranges to ensure that compositional gaps do not exist between the data sets. This report documents the results of these evaluations. Calculation-based frit assessments were performed using the DWPF PCCS glass property models and their associated MAR constraints. Evaluated parameters for coupled processing included the following transfer volumes per Sludge Receipt and Adjustment Tank batch: 5700 gallons of Tank 40 sludge, 2400-4500 gallons of the SWPF monosodium titanate (MST) and sludge solids stream, and 15,000 gallons of strip effluent. Based on these MAR assessment results, Frit 473 remains viable for SB10 processing. A target waste loading (WL) of 36% is possible for SO operation and single strike (0.4 g/L MST) coupled processing up to 600 mg/L of SB9 insoluble sludge solids at a nominal transfer volume of ~2800 gallons. Increasing the single strike transfer volume to 4500 gallons may reduce the maximum WL below 40%. Operating windows are 12 percentage points for SO processing and 14-16 percentage points for coupled processing. Frit 625 allows for a target WL of 36% and is acceptable for use during the SB9 to SB10 transition to deplete remaining inventory as needed. Operating windows are 8 percentage points for SO processing and 12-14 percentage points for coupled processing. Predictive PCCS evaluations performed at DWPF will provide insight into batch-specific acceptability at desired WLs for compositions having expected oxide ratios during processing versus the extreme vertices (corner points) evaluated in this study. The reprojected SB10 glass composition region generally overlaps the previously evaluated SB10 variability study composition region. Thus, the minor composition shift of these updated SB10 projections indicates that no additional glasses are necessary to demonstrate acceptability relative to the chemical durability of the Environmental Assessment benchmark glass and predictability using the current PCCS models for durability. Based on a comparison of the PCCS model development and validation data to the reprojected SB10 glass composition region, the viscosity and liquidus temperature models will reliably predict SB10 compositions. No additional glasses are necessary to demonstrate predictability of these models.
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗