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.
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.
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.
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.