Publications

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.

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.

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.