This work summarizes the findings of a reduced order model (ROM) study performed using Sierra ROM module Pressio_Aria on Sandia National Laboratories' (SNL) Crash-Burn L2 milestone thermal model with pristine geometry. Comparisons are made to full order model (FOM) results for this same Crash-Burn model using Sierra multiphysics module Aria.
This report details model development, theory, and a literature review focusing on the emission of contaminants on solid substrates in fires. This is the final report from a 2-year Nuclear Safety Research and Development (NSRD) project. The work represents progress towards a goal of having modeling and simulation capabilities that are sufficiently mature and accurate that they can be utilized in place of physical tests for determining safe handling practices. At present, the guidelines for safety are largely empirically based, derived from a survey of existing datasets. This particular report details the development, verification and calibration of a number of code improvements that have been implemented in the SIERRA suite of codes, and the application of those codes to three different experimental scenarios that have been subject of prior tests. The first scenario involves a contaminated PMMA slab, which is exposed to heat. The modeling involved a novel method for simulating the viscous diffusion of the particles in the slab. The second scenario involved a small pool fire of contaminated combustible liquid mimicking historical tests and finds that the release of contaminants has a high functionality with the height of the liquid in the container. The third scenario involves the burning of a contaminated tray of shredded cellulose. A novel release mechanism was formulated based on predicted progress of the decomposition of the cellulose, and while the model was found to result in release that can be tuned to match the experiments, some modifications to the model are desirable to achieve quantitative accuracy.
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.
In response to the global SARS-CoV-2 transmission pandemic, Sandia National Laboratories Rapid Lab-Directed Research and Development COVID-19 initiative has deployed a multi-physics, droplet-laden, turbulent low-Mach simulation tool to model pathogen-containing water droplets that emanate from synthetic human coughing and breathing. The low-Mach turbulent large-eddy simulation-based Eulerian/point-particle Lagrangian methodology directly couples mass, momentum, energy, and species to capture droplet evaporation physics that supports the ability to distinguish between droplets that deposit and those that persist in the environment. The cough mechanism is modeled as a pulsed spray with a prescribed log-normal droplet size distribution. Simulations demonstrate direct droplet deposition lengths in excess of three meters while the persistence of droplet nuclei entrained within a buoyant plume is noted. Including the effect of protective barriers demonstrates effective mitigation of large-droplet transport. For coughs into a protective barrier, jet impingement and large-scale recirculation can drive droplets vertically and back toward the subject while supporting persistence of droplet nuclei. Simulations in quiescent conditions demonstrate droplet preferential concentrations due to the coupling between vortex ring shedding and the subsequent advection of a series of three-dimensional rings that tilt and rise vertically due to a misalignment between the initial principle vortex trajectory and gravity. These resolved coughing simulations note vortex ring formation, roll-up and breakdown, while entraining droplet nuclei for large distances and time scales.
This interim report details model development, theory, and a literature review focusing on the evaporation induced entrainment (sub-boiling) of contaminated liquids. Entrainment from a variety of sources is the topic of DOE Handbook 3010, and this report deals more broadly with fire related airborne sources of contaminants in hazardous operations. Relatively few studies have examined sub-boiling behavior in the past, however, it can be a phenomenon that presents a fire related risk under hazardous operations. Molecular dynamics simulations are used to infer the gaseous evolution of coordinated complexes, and a model for a water/plutonium/nitrate system is deduced from the simulation results by evaluating the statistical trends of the results. Questions remain as to the chemical reactivity and longevity of entrained species. A generalized computer model capability and simple analytical model assumptions are developed for predicting the results of these and other (boiling and solid entrainment) scenarios. Verification related predictions using these models are illustrated.
This milestone campaign was focused on coupling Sandia physics codes SIERRA low Mach module Fuego and RAMSES Boltzmann transport code Sceptre(Scefire). Fuego enables simulation of low Mach, turbulent, reacting, particle laden flows on unstructured meshes using CVFEM for abnormal thermal environments throughout SNL and the larger national security community. Sceptre provides simulation for photon, neutron, and charged particle transport on unstructured meshes using Discontinuous Galerkin for radiation effects calculations at SNL and elsewhere. Coupling these ”best of breed” codes enables efficient modeling of thermal/fluid environments with radiation transport, including fires (pool, propellant, composite) as well as those with directed radiant fluxes. We seek to improve the experience of Fuego users who require radiation transport capabilities in two ways. The first is performance. We achieve this through leveraging additional computational resources for Scefire, reducing calculation times while leaving unaffected resources for fluid physics. This approach is new to Fuego, which previously utilized the same resources for both fluid and radiation solutions. The second improvement enables new radiation capabilities, including spectral (banded) radiation, beam boundary sources, and alternate radiation solvers (i.e. Pn). This summary provides an overview of these achievements.
This report presents multi-phase modeling approaches that are developed for simulating rubble fire scenarios similar to a large-scale rubble pool fire test at Sandia National Laboratories using composite materials and jet fuel. The rubble pool fire test burnt oddly shaped combustible solid objects submerged in liquid fuel. As an intermediate step toward a full scale rubble fire simulation, various model improvement tasks were performed. For modeling solid decomposition, a multi-step degradation model was used for canonical verification problems and the Chemical Percolation for Devolatilization (CPD) approach was implemented. Capabilities of Lagrangian particle approach has been extended such that a group of particles may represent a solid bulk. For gas-liquid interface, the volume of fluid (VOF) technique was implemented and relevant physics were added. The developed tools offer a potential for simulating three-phase (gas, liquid, and solid) combustion applications.
The goal of this milestone is to demonstrate effective coupling between the Sierra low-Mach module Fuego and the RAMSES Boltzmann transport (particle and radiation) code Sceptre.
With growing use of carbon fiber-epoxy in transportation systems, it is important to understand fire reaction properties of the composite to ensure passenger safety. Recently, a micro-scale pyrolysis study and macro-scale fire tests were performed using carbon fiber-epoxy at Sandia National Laboratories. Current work focuses on numerical modeling of the material conversion, pyrolysis, and gas-phase combustion that replicate the experiments. Large-eddy simulations (LES) and eddy-dissipation concept (EDC) approach are incorporated in the gas phase along with multiple relevant reaction model methods in the solid phase. The numerical methods that use multi-step pyrolysis rate expressions are validated by thermogravimetric analysis (TGA) results. The pyrolyzed fuel components participate in gas-phase combustion using a turbulent combustion model. The multi-phase combustion capability was further assessed using two cases: a single particle reaction and a solid panel exposed to strong radiant heat. The panel fire test indicates that the model accurately reproduces panel temperature profile while a weaker oxidation is predicted.