Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sequestration”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Soil Organic Matter Dynamics in the Critical Zone

Soil is the thin, vital layer of Earth’s surface that forms the foundation of the Critical Zone and sustains life. This chapter explores the intricate dynamics of soil organic matter within the critical zone, focusing on three key thematic areas: deep soil organic matter, wildfire organic matter interactions, and organic matter erosion. Soil organic matter, although a small fraction of soil mass, plays a crucial role in soil function and ecosystem stability. The complexity of soil organic matter arises from its diverse chemical composition and interactions with minerals, which influence its persistence in the environment. This chapter begins by examining deep soil organic matter, which constitutes a significant portion of global carbon storage. We discuss how deep soil organic matter, typically isolated from surface processes, may become vulnerable to decomposition and carbon release due to disturbances. Next, we explore the impacts of fire on soil organic matter, particularly the formation and stability of pyrolyzed organic matter. The decomposition of pyrolyzer organic matter is influenced by its chemical composition and the surrounding environmental conditions, with implications for carbon cycling and soil fertility in post-fire ecosystems. Finally, the chapter addresses the role of erosion in soil organic matter dynamics. Erosion, accelerated by human activities, redistributes soil organic matter across landscapes, affecting its turnover and the broader biogeochemical cycles. We consider how micro-topographic features and erosion processes interact to influence soil OM stability and carbon sequestration. By synthesizing recent advances in these areas, this chapter provides a comprehensive overview of the complex and dynamic nature of soil organic matter in the critical zone, highlighting its importance in understanding the critical zone in the face of environmental change.

Moreland, Kimber↗
DOE OSTI2024

Transgenic Sugarcane–Oilcane: An Alternative Feedstock for the Production of Drop-in Fuel and Value-Added Bioproducts

The utilization of plants and other agricultural produce can partly offset petroleum dependency for energy requirements and can potentially provide sustainable solutions to global environmental problems. To this end, synthetic biology has shown great potential in developing transgenic bioenergy grasses such as sugarcane, sorghum, miscanthus, and energy cane that hyperaccumulate energy-rich lipid molecules in their vegetative tissues, such as leaves, stems, and roots. These perennial high-biomass transgenic C4 grasses are not targeted to grow on prime agricultural land and can be dedicatedly used to produce biofuels and other value-added bioproducts. Recently, sugarcane has been metabolically engineered to sequestrate carbon from juice towards biosynthesis of lipid molecules in the vegetative tissues. Transgenic sugarcane is referred to as “oilcane.” Transgenic sugarcane–oilcane has improved energy density due to an elevated lipid content in the vegetative tissues. Transgenic bioenergy crops cater to both cellulosic sugars and vegetative lipids, hence resulting in higher biofuel yield (biodiesel/renewable diesel and bioethanol) per unit area of cultivable land. In this book chapter, the genetic engineering of sugarcane and its bioprocessing are discussed to illustrate its development and use as an alternative feedstock for the production of biofuels (bio-jet fuel, biodiesel, and bioethanol) and value-added bioproducts.

Maitra, Shraddha↗
DOE OSTI2023

Assessing the Influence of Organic Carbon, Aging Time and Temperature on Bioaccessibility of Bifenthrin

Tenax extraction, a measure of chemical desorption rates from sediments, was used to evaluate the bioaccessibility of bifenthrin in two different sediments exposed to three temperatures aged over a 56-d holding period. A 24-h single-point Tenax extraction was used and parent 14 C-bifenthrin and polar metabolites were quantified in the sediment and Tenax. Bioaccessibility of bifenthrin was inversely related to the organic carbon (OC) content in the sediment, holding time, and temperature. Furthermore, sequestration of the bifenthrin into slowly desorbing fractions within the sediment appears to have decreased degradation of the parent compound into metabolites and decreased the amount of parent compound bioaccessible for uptake by the Tenax. These results suggest that the environmental risk of bifenthrin to aquatic species is greatest immediately after the pesticide enters a waterbody after runoff, for low-OC content sediments, and in areas or seasons where water temperatures are colder.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Effectiveness and outcomes of invasive species removal in Hawaiian streams

Despite aquatic invasive species (AIS) being a widely recognized threat to native biodiversity in streams on subtropical and tropical oceanic islands across the Pacific, AIS management has been constrained by concerns that methods of removal could result in more harm than good. In this study, we assessed the effectiveness and outcomes of combining hand-netting and streamside sequestration of native species with AIS removal via three-pass electrofishing tailored for Hawaiian Island watersheds. Removals were conducted in 13 watersheds on the island of O‘ahu following a before-after, control-impact design to assess how AIS removal and post-removal community reassembly varied by surface flow and target species of interest. We found that removals resulted in a sharp reduction in AIS densities without altering native species densities. Removal efficacy was greatest in streams with intermediate mean discharge irrespective of target group or pre-removal AIS densities. Long-term monitoring demonstrated that removals led to persistent shifts in community composition reflecting sustained reductions in AIS densities. These findings indicate that electrofishing can be a valuable tool for AIS management in Hawaiian and other subtropical Pacific island streams, so long as steps are taken to minimize potential harm to at-risk and culturally important native species.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Assessing microbial residues in soil as a potential carbon sink and moderator of carbon use efficiency

Abstract A longstanding assumption of glucose tracing experiments is that all glucose is microbially utilized during short incubations of ≤2 days to become microbial biomass or carbon dioxide. Carbon use efficiency (CUE) estimates have consequently ignored the formation of residues (non-living microbial products) although such materials could represent an important sink of glucose that is prone to stabilization as soil organic matter. We examined the dynamics of microbial residue formation from a short tracer experiment with frequent samplings over 72 h, and conducted a meta-analysis of previously published glucose tracing studies to assess the generality of these experimental results. Both our experiment and meta-analysis indicated 30–34% of amended glucose-C ( 13 C or 14 C) was in the form of residues within the first 6 h of substrate addition. We expand the conventional efficiency calculation to include residues in both the numerator and denominator of efficiency, thereby deriving a novel metric of the potential persistence of glucose-C in soil as living microbial biomass plus residues (‘carbon stabilization efficiency’). This new metric indicates nearly 40% of amended glucose-C persists in soil 180 days after amendment, the majority as non-biomass residues. Starting microbial biomass and clay content emerge as critical factors that positively promote such long term stabilization of labile C. Rapid residue production supports the conclusion that non-growth maintenance activity can illicit high demands for C in soil, perhaps equaling that directed towards growth, and that residues may have an underestimated role in the cycling and sequestration potential of C in soil.

Geyer, Kevin (ORCID:0000000252182745)↗
DOE OSTI2022

TOUGH3-FLAC3D: a modeling approach for parallel computing of fluid flow and geomechanics

The recent development of the TOUGH3 code allows for a faster and more reliable fluid flow simulator. At the same time, new versions of FLAC3D are released periodically, allowing for new features and faster execution. In this paper, we present the first implementation of the coupling between TOUGH3 and FLAC3Dv6/7, maintaining parallel computing capabilities for the coupled fluid flow and geomechanical codes. We compare the newly developed version with analytical solutions and with the previous approach, and provide some performance analysis on different meshes and varying the number of running processors. Finally, we present two case studies related to fault reactivation during CO 2 sequestration and nuclear waste disposal. The use of parallel computing allows for meshes with a larger number of elements, and hence more detailed understanding of thermo-hydro-mechanical processes occurring at depth.

58 GEOSCIENCES↗
DOE OSTI2024

An integrated approach to derive relative permeability from capillary pressure

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability curve can be determined. Relative permeability is a numerical description of the interaction between two or more fluids and the porous media. It is a critical parameter for various tools that characterize subsurface multiphase flow systems, such as numerical simulation for carbon sequestration, oil and gas development, and groundwater contamination remediation. Therefore, it is critical to get a good statistical distribution of relative permeability in the porous media under study. Empirical formula for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility that is required to match laboratory-derived relative permeability curves. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability curves. Mercury intrusion capillary pressure (MICP) data from the Morrow 'B' Sandstone coupled with interfacial tension and contact angle measurements were used to create a suite of relative permeability curves. Furthermore, these curves were then calibrated to a small sample of existing laboratory curves to elucidate common fitting parameters for the formation that were then used to create relative permeability curves from MICP data that does not have an associated laboratory-measured relative permeability curve.

58 GEOSCIENCES↗
DOE OSTI2022

Improving sustainable tropical forest management with voluntary carbon markets

Due to a rapidly changing climate, voluntary carbon markets are gaining momentum and should be leveraged to improve and expand tropical sustainable forest management plans, limiting carbon emissions and enhancing critical carbon sinks. By sequestering more carbon than any other terrestrial ecosystem — ~1 Pg C yr –1 — tropical forests provide crucial natural climate solutions and opportunities in the evolving voluntary carbon market. Here, we argue that some issues with the current sustainable management of tropical forests can be addressed using carbon-focused sustainable forest management (SFM + C) to leverage financial resources for tropical forest carbon storage and sequestration. We suggest an extended harvest cycle in SFM + C and calculate an associated potential increase in aboveground carbon stocks of commercial timber of 1.26 Mg C ha –1 after each cycle in the Brazilian Amazon. The additional carbon storage due to a longer harvest cycle can generate carbon credits worth 152.6 (SD 9.2) US dollars per hectare in 40 years. Considering an average cost of 180 BRL per m3 of commercial timber delivered to the sawmill, an SFM + C plan with a 40-year cycle could generate 28.7% (SD 2.5) more profit than 35-year cycles by combining timber and carbon revenues. A robust carbon price could incentivize the further extension of harvest cycles, providing a monetary return that offsets the opportunity cost intrinsic to harvesting under longer cycles. Lastly, we highlight research needs to support tropical SFM + C, which can be part of a global collective effort to limit global warming to below 2 °C above pre-industrial levels.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Rapid in situ nutrient element distribution in plants and soils using laser-induced breakdown spectroscopy ($\mathrm{LIBS}$)

The aim of this study is to develop and test the applicability of a rapid in situ plant chemistry profiling technique to determine elemental composition of small-volume plant and soil samples obtained from a woody bioenergy crop species, Populus trichocarpa. Expanding the research tools available to characterize the nutrient element correlations among plant tissue types and soil depths is a critical need in the path of understanding productivity and adaptation of plants to variations in external abiotic and biotic factors and developing sustainable perennial bioenergy crops that are co-optimized for biomass valorization aboveground and carbon sequestration belowground. Several plant root, stem, and soil samples were tested using laser-induced breakdown spectroscopy (LIBS) to evaluate the presence and distribution of nutrient elements. Samples were tested as collected and after being dried and cross sectioned to evaluate the effectiveness of using LIBS for in situ analysis on plant samples. The collected LIBS spectra show the elemental peaks were the same in both the as collected and prepared samples for roots and stems. Qualitative amounts of elements such as H, C, N, O, Li, Na, Mg, K, Ca, Fe, Al, and Si were able to be identified rapidly in raw samples. Here in this report we demonstrate suitability of LIBS in obtaining rapid, in situ, elemental distribution in plant and soil samples, utilizing only small sample volumes and minimal sample preparation. This demonstration opens up a new rapid phenotyping avenue necessary to fill the asymmetrical knowledge gaps in belowground performance of plant systems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Breakthrough innovations in carbon dioxide mineralization for a sustainable future

Greenhouse gas emissions and climate change concerns have prompted worldwide initiatives to lower carbon dioxide (CO 2 ) levels and prevent them from rising in the atmosphere, thereby controlling global warming. Effective CO 2 management through carbon capture and storage is essential for safe and permanent storage, as well as synchronically meeting carbon reduction targets. Lowering CO 2 emissions through carbon utilization can develop a wide range of new businesses for energy security, material production, and sustainability. CO 2 mineralization is one of the most promising strategies for producing thermodynamically stable solid calcium or magnesium carbonates for long-term sequestration using simple chemical reactions. Current advancements in CO 2 mineralization technologies, focusing on pathways and mechanisms using different industrial solid wastes, including natural minerals as feedstocks, are briefly discussed. However, the operating costs, energy consumption, reaction rates, and material management are major barriers to the application of these technologies in CO 2 mineralization. Further, the optimization of operating parameters, tailor-made equipment, and smooth supply of waste feedstocks require more attention to make the carbon mineralization process economically and commercially viable. Here, carbonation mechanisms, technological options to expedite mineral carbonation, environmental impacts, and prospects of CO 2 mineralization technologies are critically evaluated to suggest a pathway for mitigating climate change in the future. The integration of industrial wastes and brine with the CO 2 mineralization process can unlock its potential for the development of novel chemical pathways for the synthesis of calcium or magnesium carbonates, valuable metal recovery, and contribution to sustainability goals while reducing the impact of global warming.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A Fresh Take: Seasonal Changes in Terrestrial Freshwater Inputs Impact Salt Marsh Hydrology and Vegetation Dynamics

Abstract Salt marshes exist at the terrestrial-marine interface, providing important ecosystem services such as nutrient cycling and carbon sequestration. Tidal inputs play a dominant role in salt marsh porewater mixing, and terrestrially derived freshwater inputs are increasingly recognized as important sources of water and solutes to intertidal wetlands. However, there remains a critical gap in understanding the role of freshwater inputs on salt marsh hydrology, and how this may impact marsh subsurface salinity and plant productivity. Here, we address this knowledge gap by examining the hydrologic behavior, porewater salinity, and pickleweed ( Sarcocornia pacifica also known as Salicornia pacifica) plant productivity along a salt marsh transect in an estuary along the central coast of California. Through the installation of a suite of hydrometric sensors and routine porewater sampling and vegetation surveys, we sought to understand how seasonal changes in terrestrial freshwater inputs impact salt marsh ecohydrologic processes. We found that salt marsh porewater salinity, shallow subsurface saturation, and pickleweed productivity are closely coupled with elevated upland water level during the winter and spring, and more influenced by tidal inputs during the summer and fall. This seasonal response indicates a switch in salt marsh hydrologic connectivity with the terrestrial upland that impacts ecosystem functioning. Through elucidating the interannual impacts of drought on salt marsh hydrology, we found that the severity of drought and historical precipitation can impact contemporary hydrologic behavior and the duration and timing of the upland-marsh hydrologic connectivity. This implies that the sensitivity of salt marshes to climate change involves a complex interaction between sea level rise and freshwater inputs that vary at seasonal to interannual timescales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Practical Guide to Measuring Wetland Carbon Pools and Fluxes

Abstract Wetlands cover a small portion of the world, but have disproportionate influence on global carbon (C) sequestration, carbon dioxide and methane emissions, and aquatic C fluxes. However, the underlying biogeochemical processes that affect wetland C pools and fluxes are complex and dynamic, making measurements of wetland C challenging. Over decades of research, many observational, experimental, and analytical approaches have been developed to understand and quantify pools and fluxes of wetland C. Sampling approaches range in their representation of wetland C from short to long timeframes and local to landscape spatial scales. This review summarizes common and cutting-edge methodological approaches for quantifying wetland C pools and fluxes. We first define each of the major C pools and fluxes and provide rationale for their importance to wetland C dynamics. For each approach, we clarify what component of wetland C is measured and its spatial and temporal representativeness and constraints. We describe practical considerations for each approach, such as where and when an approach is typically used, who can conduct the measurements (expertise, training requirements), and how approaches are conducted, including considerations on equipment complexity and costs. Finally, we review key covariates and ancillary measurements that enhance the interpretation of findings and facilitate model development. The protocols that we describe to measure soil, water, vegetation, and gases are also relevant for related disciplines such as ecology. Improved quality and consistency of data collection and reporting across studies will help reduce global uncertainties and develop management strategies to use wetlands as nature-based climate solutions.

Bansal, Sheel (ORCID:0000000312331707)↗
DOE OSTI2024

Advancing quantitative models of soil microbiology, ecology, and biochemistry

Rapid advances in soil biology are increasingly reflected in mathematical models. From molecular to global scales, models contribute to fundamental understanding and prediction of critical soil processes, such as carbon sequestration. Recent models explicitly incorporate microbial control over soil carbon and nutrient cycling, an approach that offers greater biological realism but raises new challenges in model parameterization and stability. Although they have proliferated at community and ecosystem scales, microbial-explicit models have yet to be incorporated into most Earth system models. Applying soil models across scales requires integration with data to support parameterization and validation using approaches such as Bayesian data assimilation. In conclusion, to catalyze future progress, we recommend breaking down barriers between modeling and empirical disciplines while broadening access to the computational infrastructure that supports soil biological modeling.

Xie, Wally↗
DOE OSTI2020

Onset dynamics of air-water menisci on rock fracture surfaces

Information on menisci dynamics and equilibrium contact angles, θ e , is needed for modeling multiphase flow of fluids in geologic formations. The wettability of fracture faces is relevant to applications such as waterflooding in enhanced oil recovery and evaluation of caprock integrity for geologic carbon sequestration. We investigated the vertical capillary rise of air-water menisci on exposed fracture faces for a sandstone, a shale, and three granites. The experimental design approximated that employed in the Wilhelmy plate method. Replicate fracture faces were prepared by fracturing cylindrical rock cores using the Brazilian method and splitting the fractured cores apart along their fracture planes. Meniscus onset was imaged from a cross-sectional perspective using dynamic neutron radiography, and quantified with change point analysis. An existing model for meniscus onset fitted the experimental data very well (median R 2 = 0.96). Capillary rise was typified by t 2/3 scaling (where t is time), followed by a constant equilibrium height, z e . The average time taken for the menisci to reach z e was 0.39 s, with no differences between rock types. The menisci achieved a greater vertical extent on the fracture faces of the sedimentary rocks as compared to those of the granites. Apparent θ e values were calculated from the best fit estimates of z e using the Wilhelmy plate equation for a smooth and planar vertical surface. Despite the fractal roughness of some of the fracture faces, there was a significant linear relationship (R 2 = 0.80) between the apparent θ e values and intrinsic θ e values measured on polished surfaces of the same rocks using the sessile drop method. The regression equation indicated the apparent θ e values were consistently greater than the intrinsic values. This overestimation may be due to averaging of pixels at the air-water-solid interface, a required step in the image analysis. Further exploration of the range of applicability and measurement accuracy of this proposed new technique would be valuable.

58 GEOSCIENCES↗
DOE OSTI2020

The role of the spatial heterogeneity and correlation length of surface wettability on two-phase flow in a CO 2 -water-rock system

This study characterized and modeled heterogeneous surface wettability in sandstone and investigated the role of spatial heterogeneity and correlation length of surface wettability on relative permeability in a supercritical CO 2 (scCO 2 )-brine-rock system. Understanding the role of wettability heterogeneity on relative permeability is essential to geological CO 2 sequestration, oil and gas recovery, and contaminated groundwater remediation. Although numerous studies have attempted to understand the influences of surface wettability, capillary number (Ca), and viscosity ratio, the role of the spatial variation and correlation length of surface wettability on two-phase flow in three-dimensional (3D) porous media has not been unraveled due to the challenges in the measurement and representation of realistic rock surface wettability. In this work, we conducted in-situ measurements of surface contact angle (CA) in a Bentheimer sandstone after CO 2 flooding using micro-computed tomography (micro-CT), and found that the pore-scale CA distribution on rock surfaces followed a log-normal distribution associated with a spatial correlation length. Based on the statistical information from CT scanning, a Gaussian random field was used to model CA distributions that had desired standard deviations and spatial correlation lengths, which were then adjusted within a certain range of values for sensitivity analyses to study their combined effects on the two-phase flow in the porous medium using the lattice Boltzmann (LB) method. The LB two-phase flow simulation was accelerated using hybrid, multicore parallel computing to overcome the challenges in simulating multiphase flow in a large 3D domain having 800 × 800 × 600 nodes. The simulation results showed that the surface wettability heterogeneity (i.e., standard deviation of CA) had a lesser effect on the relative permeability of the wetting fluid (water) but a more significant impact on the relative permeability of the non-wetting fluid (scCO 2 ). The Corey model was used to fit the LB-simulated relative permeability curves of water and scCO 2 and showed that the variations in the relative permeability curves for both water and scCO 2 increased as the standard deviation and spatial correlation length of CA increased. This study illustrated that the assumption of homogeneous surface wettability may cause errors in multiphase flow simulations. Furthermore, the impacts of both the standard deviation and spatial correlation length of CAs should be accounted for. This is the first study that explored the spatial correlation lengths associated with CA distributions on sandstone surfaces and comprehensively investigated the roles of both spatial variation and correlation length of CA on two-phase flow properties in 3D porous media. The optimized LB multiphase flow model was proved a powerful tool to study the interplays and combined effects of these statistical parameters, which had critical applications in numerous natural and engineering processes that involved multiphase flow in porous media.

58 GEOSCIENCES↗
DOE OSTI2021

Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.

Preconditioner, Nonlinear, Porous media, Multiphas↗
DOE OSTI2022

Non-intrusive reduced order modeling of natural convection in porous media using convolutional autoencoders: Comparison with linear subspace techniques

Natural convection in porous media is a highly nonlinear multiphysical problem relevant to many engineering applications (e.g., the process of CO 2 sequestration). Here, we extend and present a non-intrusive reduced order model of natural convection in porous media employing deep convolutional autoencoders for the compression and reconstruction and either radial basis function (RBF) interpolation or artificial neural networks (ANNs) for mapping parameters of partial differential equations (PDEs) on the corresponding nonlinear manifolds. To benchmark our approach, we also describe linear compression and reconstruction processes relying on proper orthogonal decomposition (POD) and ANNs. Further, we present comprehensive comparisons among different models through three benchmark problems. The reduced order models, linear and nonlinear approaches, are much faster than the finite element model, obtaining a maximum speed-up of 7 × 10 6 because our framework is not bound by the Courant–Friedrichs–Lewy condition; hence, it could deliver quantities of interest at any given time contrary to the finite element model. Our model’s accuracy still lies within a relative error of 7% in the worst-case scenario. We illustrate that, in specific settings, the nonlinear approach outperforms its linear counterpart and vice versa. We hypothesize that a visual comparison between principal component analysis (PCA) and t-Distributed Stochastic Neighbor Embedding (t-SNE) could indicate which method will perform better prior to employing any specific compression strategy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Newton trust-region methods with primary variable switching for simulating high temperature multiphase porous media flow

Coupling multiphase flow with energy transport due to high temperature heat sources introduces significant new challenges since boiling and condensation processes can lead to dry-out conditions with subsequent re-wetting. The transition between two-phase and single-phase behavior can require changes to the primary dependent variables adding discontinuities as well as extending constitutive nonlinear relations to extreme physical conditions. Practical simulations of large-scale engineered domains lead to Jacobian systems with a very large number of unknowns that must be solved efficiently using iterative methods in parallel on high-performance computers. Performance assessment of potential nuclear repositories, carbon sequestration sites and geothermal reservoirs can require numerous Monte-Carlo simulations to explore uncertainty in material properties, boundary conditions, and failure scenarios. Due to the numerical challenges, standard NR iteration may not converge over the range of required simulations and require more sophisticated optimization method like trust-region. In this study, we use the open-source simulator PFLOTRAN for the important practical problem of the safety assessment of future nuclear waste repositories in the U.S. DOE geologic disposal safety assessment Framework. The simulator applies the PETSc parallel framework and a backward Euler, finite volume discretization. We demonstrate failure of the conventional NR method and the success of trust-region modifications to Newton’s method for a series of test problems of increasing complexity. Trust-region methods essentially modify the Newton step size and direction under some circumstances where the standard NR iteration can cause the solution to diverge or oscillate. Furthermore, we show how the Newton Trust-Region method can be adapted for Primary Variable Switching (PVS) when the multiphase state changes due to boiling or condensation. The simulations with high-temperature heat sources which led to extreme nonlinear processes with many state changes in the domain did not converge with NR, but they do complete successfully with the trust-region methods modified for PVS. This implementation effectively decreased weeks of simulation time needing manual adjustments to complete a simulation down to a day. Finally, we show the strong scalability of the methods on a single node and multiple nodes in an HPC cluster.

54 ENVIRONMENTAL SCIENCES↗