Publications

The SNL Sierra Mechanics code suite is designed to enable simulation of complex multiphysicsscenarios. The code suite is composed of several specialized applications which can operate either instandalone mode or coupled with each other. Arpeggio is a supported utility that enables loose couplingof the various Sierra Mechanics applications by providing access to Framework services that facilitatethe coupling.

Presented in this document is a portion of the tests that exist in the Sierra Thermal/Fluids verificationtest suite. Each of these tests is run nightly with the Sierra/TF code suite and the results of the testchecked under mesh refinement against the correct analytic result. For each of the tests presented in thisdocument the test setup, derivation of the analytic solution, and comparison of the code results to theanalytic solution is provided.

The SIERRA Low Mach Module: Fuego, henceforth referred to as Fuego, is the key element of theASC fire environment simulation project. The fire environment simulation project is directed atcharacterizing both open large-scale pool fires and building enclosure fires. Fuego represents theturbulent, buoyantly-driven incompressible flow, heat transfer, mass transfer, combustion, soot, andabsorption coefficient model portion of the simulation software.

Abstract not provided.

Prediction of flow, transport, and deformation in fractured and porous media is critical to improving our scientific understanding of coupled thermal-hydrological-mechanical processes related to subsurface energy storage and recovery, nonproliferation, and nuclear waste storage. Especially, earth rock response to changes in pressure and stress has remained a critically challenging task. In this work, we advance computational capabilities for coupled processes in fractured and porous media using Sandia Sierra Multiphysics software through verification and validation problems such as poro-elasticity, elasto-plasticity and thermo-poroelasticity. We apply Sierra software for geologic carbon storage, fluid injection/extraction, and enhanced geothermal systems. We also significantly improve machine learning approaches through latent space and self-supervised learning. Additionally, we develop new experimental technique for evaluating dynamics of compacted soils at an intermediate scale. Overall, this project will enable us to systematically measure and control the earth system response to changes in stress and pressure due to subsurface energy activities.

Accurate modeling of subsurface flow and transport processes is vital as the prevalence of subsurface activities such as carbon sequestration, geothermal recovery, and nuclear waste disposal increases. Computational modeling of these problems leverages poroelasticity theory, which describes coupled fluid flow and mechanical deformation. Although fully coupled monolithic schemes are accurate for coupled problems, they can demand significant computational resources for large problems. In this work, a fixed stress scheme is implemented into the Sandia Sierra Multiphysics toolkit. Two implementation methods, along with the fully coupled method, are verified with one-dimensional (1D) Terzaghi, 2D Mandel, and 3D Cryer sphere benchmark problems. The impact of a range of material parameters and convergence tolerances on numerical accuracy and efficiency was evaluated. Overall the fixed stress schemes achieved acceptable numerical accuracy and efficiency compared to the fully coupled scheme. However, the accuracy of the fixed stress scheme tends to decrease with low permeable cases, requiring the finer tolerance to achieve a desired numerical accuracy. For the fully coupled scheme, high numerical accuracy was observed in most of cases except a low permeability case where an order of magnitude finer tolerance was required for accurate results. Finally, a two-layer Terzaghi problem and an injection–production well system were used to demonstrate the applicability of findings from the benchmark problems for more realistic conditions over a range of permeability. Simulation results suggest that the fixed stress scheme provides accurate solutions for all cases considered with the proper adjustment of the tolerance. This work clearly demonstrates the robustness of the fixed stress scheme for coupled poroelastic problems, while a cautious selection of numerical tolerance may be required under certain conditions with low permeable materials.

To mitigate adverse effects from molten corium following a reactor pressure vessel failure (RPVF), some new reactor designs employ a core catcher and a sacrificial material (SM), such as ceramic or concrete, to stabilize the molten corium and avoid containment breach. Existing reactors cannot easily be modified to include these SMs but could be modified to allow injectable cooling materials. Current reactor designs are limited to using water to stabilize the corium, but this can create other issues such as reaction of water with the concrete forming hydrogen gas. The novel SM proposed here is a granular carbonate mineral that can be used in existing light water reactor plants. The granular carbonate will decompose when exposed to heat, inducing an endothermic reaction to quickly solidify the corium in place and producing a mineral oxide and carbon dioxide. Corium spreading is a complex process strongly influenced by coupled chemical reactions, including decay heat from the corium, phase change, and reactions between the concrete containment and available water. A recently completed Sandia National Laboratories laboratory directed research and development (LDRD) project focused on two research areas: experiments to demonstrate the feasibility of the novel SM concept, and modeling activities to determine the potential applications of the concept to actual nuclear plants. Small-scale experiments using lead oxide (PbO) as a surrogate for molten corium demonstrate that the reaction of the SM with molten PbO results in a fast solidification of the melt due to the endothermic carbonate decomposition reaction and the formation of open pore structures in the solidified PbO from CO2 released during the decomposition. A simplified carbonate decomposition model was developed to predict thermal decomposition of carbonate mineral in contact with corium. This model was incorporated into MELCOR, a severe accident nuclear reactor code. A full-plant MELCOR simulation suggests that by the introduction of SM to the reactor cavity prior to RPVF ex-vessel accident progression, e.g., core-concrete interaction and core spreading on the containment floor, could be delayed by at least 15 h; this may be enough for additional accident management to be implemented to alleviate the situation.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Melting and flowing of aluminum alloys is a challenging problem for computational codes. Unlike most common substances, the surface of an aluminum melt exhibits rapid oxidation and elemental migration, and like a bag filled with water can remain 2-dimensionally unruptured while the metal inside is flowing. Much of the historical work in this area focuses on friction welding and neglects the surface behavior due to the high stress of the application. We are concerned with low-stress melting applications, in which the bag behavior is more relevant. Adapting models and measurements from the literature, we have developed a formulation for the viscous behavior of the melt based on an abstraction of historical measurement, and a construct for the bag behavior. These models are implemented and demonstrated in a 3D level-set multi-phase solver package, SIERRA/Aria. A series of increasingly complex simulation scenarios are illustrated that help verify implementation of the models in conjunction with other required model components like convection, radiation, gravity, and surface interactions.

Recent advances in drilling technology, especially horizontal drilling, have prompted a renewed interest in the use of closed loop geothermal energy extraction systems. Deeply placed closed loops in hot wet or dry rock reservoirs offer the potential to exploit the vast thermal energy in the subsurface. To better understand the potential and limitations for recovering thermal and mechanical energy from closed-loop geothermal systems (CLGS), a collaborative study is underway to investigate an array of system configurations, working fluids, geothermal reservoir characteristics, operational periods, and heat transfer enhancements (Parisi et al., 2021; White et al., 2021). This paper presents numerical results for the heat exchange between a closed loop system (single U-tube) circulating water as the working fluid in a hot rock reservoir. The characteristics of the reservoir are based on the Frontier Observatory for Research in Geothermal Energy (FORGE) site, near Milford Utah. To determine optimal system configurations, a mechanical (electrical) objective function is defined for a bounded optimization study over a specified design space. The objective function includes a surface plant thermal to mechanical energy conversion factor, pump work, and an energy drilling capital cost. To complement the optimization results, detailed parametric studies are also performed. The numerical model is built using the Sandia National Laboratories (SNL) massively parallel Sierra computational framework, while the optimization and parametric studies are driven using the SNL Dakota software package. Together, the optimization and parametric studies presented in this paper will help assess the impact of CLGS parameters (e.g., flow rate, tubing length and diameter, insulation length, etc.) on CLGS performance and optimal energy recovery.

Abstract not provided.

This project combines several new concepts to create a boundary layer transition prediction capability that is suitable for analyzing modern hypersonic flight vehicles. The first new concept is the use of ''optimization'' methods to detect the hydrodynamic instabilities that cause boundary layer transition; the use of this method removes the need for many limiting assumptions of other methods and enables quantification of the interactions between boundary layer instabilities and the flow field imperfections that generate them. The second new concept is the execution of transition analysis within a conventional hypersonics CFD code, using the same mesh and numerical schemes for the transition analysis and the laminar flow simulation. This feature enables rapid execution of transition analysis with less user oversight required and no interpolation steps needed.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The catastrophic nuclear reactor accident at Fukushima damaged public confidence in nuclear energy and a demand for new engineered safety features that could mitigate or prevent radiation releases to the environment in the future. We have developed a novel use of sacrificial material (SM) to prevent the molten corium from breaching containment during accidents as well as a validated, novel, high-fidelity modeling capability to design and optimize the proposed concept. Some new reactor designs employ a core catcher and a SM, such as ceramic or concrete, to slow the molten corium and avoid the breach of the containment. However, existing reactors cannot easily be modified to include these SMs but could be modified to allow injectable cooling materials (current designs are limited to water). The SM proposed in this Laboratory Development Research and Development (LDRD) project is based on granular carbonate minerals that can be used in existing light water reactor plants. This new SM will induce an endothermic reaction to quickly freeze the corium in place, with minimal hydrogen explosion and maximum radionuclide retention. Because corium spreading is a complex process strongly influenced by coupled chemical reactions (with underlying containment material and especially with the proposed SM), decay heat and phase change. No existing tool is available for modeling such a complex process. This LDRD project focused on two research areas: experiments to demonstrate the feasibility of the novel SM concept, and modeling activities to determine the potential applications of the concept to actual nuclear plants. We have demonstrated small-scale to large-scaled experiments using lead oxide (Pb0) as surrogate for molten corium, which showed that the reaction of the SM with molten Pb0 results in a fast solidification of the melt and the formation of open pore structures in the solidified Pb0 because of CO₂ released from the carbonate decomposition.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

When the core is breached during a severe nuclear accident, a molten mixture of nuclear fuel, cladding, and structural supports is discharged from the reactor vessel. This molten mixture of ceramic and metal is often referred to as “corium”. Predicting the flow and solidification of corium poses challenges for numerical models due to the presence of large Peclet numbers when convective transport dominates the physics. Here, we utilize a control volume finite-element method (CVEM) discretization to stabilize the advection dominated flow and heat transport. This CVFEM approach is coupled with the conformal decomposition finite-element method (CDFEM), which tracks the corium/air interface on an existing background mesh. CDFEM is a sharp-interface method, allowing the direct discretization of the corium front. This CVFEM-CDFEM approach is used to model the spreading of molten corium in both two- and three-dimensions. The CVFEM approach is briefly motivated in a comparison with a streamwise upwind/Petrov-Galerkin (SUPG) stabilized finite-element method, which was not able to suppress spurious temperature oscillations in the simulations. Our model is compared directly with the FARO L26 corium spreading experiments and with previous numerical simulations, showing both quantitative and qualitative agreement with those studies.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Imaging techniques for the analysis of porous structures have revolutionized our ability to quantitatively characterize geomaterials. Digital representations of rock from CT images and physics modeling based on these pore structures provide the opportunity to further advance our quantitative understanding of fluid flow, geomechanics, and geochemistry, and the emergence of coupled behaviors. Additive manufacturing, commonly known as 3D printing, has revolutionized production of custom parts with complex internal geometries. For the geosciences, recent advances in 3D printing technology may be co-opted to print reproducible porous structures derived from CT-imaging of actual rocks for experimental testing. The use of 3D printed microstructure allows us to surmount typical problems associated with sample-to-sample heterogeneity that plague rock physics testing and to test material response independent from pore-structure variability. Together, imaging, digital rocks and 3D printing potentially enables a new workflow for understanding coupled geophysical processes in a real, but well-defined setting circumventing typical issues associated with reproducibility, enabling full characterization and thus connection of physical phenomena to structure. Here we report on our research exploring the possibilities that these technologies can bring to geosciences for coupled multiscale experimental and numerical analysis using 3D printed fractured rock specimens.

Abstract not provided.

The migration and trapping of supercritical CO₂ (scCO₂) in geologic carbon storage is strongly dependent on the geometry and wettability of the pore network in the reservoir rock. During displacement, resident fluids may become trapped in the pits of a rough pore surface forming an immiscible two-phase fluid interface with the invading fluid, allowing apparent slip flow at this interface. We present a two-phase fluid dynamics model, including interfacial tension, to characterize the impact of mineral surface roughness on this slip flow. We show that the slip flow can be cast in more familiar terms as a contact-angle (wettability)-dependent effective permeability to the invading fluid, a nondimensional measurement which relates the interfacial slip to the pore geometry. The analysis shows the surface roughness-induced slip flow can effectively increase or decrease this effective permeability, depending on the wettability and roughness of the mineral surfaces. Configurations of the pore geometry where interfacial slip has a tangible influence on permeability have been identified. The results suggest that for large roughness features, permeability to CO₂ may be enhanced by approximately 30% during drainage, while the permeability to brine during reimbibition may be enhanced or diminished by 60%, depending on the contact angle with the mineral surfaces and degrees of roughness. For smaller roughness features, the changes in permeability through interfacial slip are small. As a result, a much larger range of effective permeabilities are suggested for general fluid pairs and contact angles, including occlusion of the pore by the trapped phase.

The migration and trapping of supercritical CO2 (scCO2) in geologic carbon storage is strongly dependent on the geometry and wettability of the pore network in the reservoir rock. During displacement, resident fluids may become trapped in the pits of a rough pore surface forming an immiscible two-phase fluid interface with the invading fluid, allowing apparent slip flow at this interface. We present a two-phase fluid dynamics model, including interfacial tension, to characterize the impact of mineral surface roughness on this slip flow. We show that the slip flow can be cast in more familiar terms as a contact-angle (wettability)-dependent effective permeability to the invading fluid, a nondimensional measurement which relates the interfacial slip to the pore geometry. The analysis shows the surface roughness-induced slip flow can effectively increase or decrease this effective permeability, depending on the wettability and roughness of the mineral surfaces. Configurations of the pore geometry where interfacial slip has a tangible influence on permeability have been identified. The results suggest that for large roughness features, permeability to CO2 may be enhanced by approximately 30% during drainage, while the permeability to brine during reimbibition may be enhanced or diminished by 60%, depending on the contact angle with the mineral surfaces and degrees of roughness. For smaller roughness features, the changes in permeability through interfacial slip are small. A much larger range of effective permeabilities are suggested for general fluid pairs and contact angles, including occlusion of the pore by the trapped phase.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

An enriched finite element method is described for capillary hydrodynamics including dynamic wetting. The method is enriched via the Conformal Decomposition Finite Element Method (CDFEM). Two formulations are described, one with first-order accuracy and one with second-order accuracy in time. Both formulations utilize a semi-implicit form for the surface tension that is shown to effectively circumvent the explicit capillary time step limit. Sharp interface boundary conditions are developed for capturing the dynamic contact angle as the fluid interface moves along the wall. By virtue of the CDFEM, the contact line is free to move without risk of mesh tangling, but is sharply captured. Multiple problems are used to demonstrate the effectiveness of the methods.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.