The breakup of liquid drops is an important phenomenology for many applications. We approach this problem with the objective to improve methods for handling the modeling of the impulse and impact dispersal of liquids in transportation accident scenarios. These scenarios can be distinguished from many other simpler problems, due to the quantity of liquid and the complexity of the intermediate liquid morphology. These differences necessitate alternative approaches to the problem. We have recently implemented a model for inter-particle forces between particles in a Lagrangian/Eulerian CFD code. The inter-particle force model is inspired by molecular dynamics methods, and employs a Lennard-Jones (LJ) attractive force and a spring-based repulsive force. The LJ parameters are related to the bulk fluid properties through a theoretical relationship model. Methods are necessary for modifying the single particle aerodynamic drag term, depending on the new notion of particle connectivity. We want to evaluate these methods for potential utilization in practical simulations. Classical breakup tests for drops in flows suggest a critical Weber number relating to the onset of breakup for a drop. We seek to replicate these data with our model methods as a preliminary step before deploying the method in larger scale practical environments.
Safety basis analysts throughout the U.S. Department of Energy (DOE) complex rely heavily on the information provided in the DOE Handbook, DOE - HDBK - 3010, Airborne Release Fractions/Rates and Respirable Fractions for Nonreactor Nuclear Facilities, to determine radionuclide source terms. In calculating source terms, analysts tend to use the DOE Handbook's bounding values on airborne release fractions (ARFs) and respirable fractions (RFs) for various categories of insults (representing potential accident release categories). This is typically due to both time constraints and the avoidance of regulatory critique. Unfortunately, these bounding ARFs/RFs represent extremely conservative values. Moreover, they were derived from very limited small-scale bench/laboratory experiments and/or from engineered judgment. Thus, the basis for the data may not be representative of the actual unique accident conditions and configurations being evaluated. The goal of this research is to develop a more accurate and defensible method to determine bounding values for the DOE Handbook using state-of-art multi-physics-based computer codes. This enables us to better understand the fundamental physics and phenomena associated with the types of accidents in the handbook. In this year, this research included improvements of the high-fidelity codes to model particle resuspension and multi-component evaporation for fire scenarios. We also began to model ceramic fragmentation experiments, and to reanalyze the liquid fire and powder release experiments that were done last year. The results show that the added physics better describes the fragmentation phenomena. Thus, this work provides a low-cost method to establish physics-justified safety bounds by taking into account specific geometries and conditions that may not have been previously measured and/or are too costly to perform.
In this work we have presented a particle resuspension model implemented in the SNL code SIERRA/Fuego, which can be used to model particle dispersal and resuspension from surfaces. The method demonstrated is applicable to a class of particles, but would require additional parametric fits or physics models for extension to other applications, such as wetted particles or walls. We have demonstrated the importance of turbulent variations in the wall shear stress when considering resuspension, and implemented both shear stress variation models and stochastic resuspension models (not shown in this work). These models can be used in simulations with of physically realistic scenarios to augment lab-scale DOE Handbook data for airborne release fractions and respirable fractions in order to provide confidences for safety analysts and facility designers to apply in their analyses at DOE sites. Future work on this topic will involve validation of the presented model against experimental data and extension of the empirical models to be applicable to different classes of particles and surfaces.
In this work we have presented a particle resuspension model implemented in the SNL code SIERRA/Fuego, which can be used to model particle dispersal and resuspension from surfaces. The method demonstrated is applicable to a class of particles, but would require additional parametric fits or physics models for extension to other applications, such as wetted particles or walls. We have demonstrated the importance of turbulent variations in the wall shear stress when considering resuspension, and implemented both shear stress variation models and stochastic resuspension models (not shown in this work). These models can be used in simulations with of physically realistic scenarios to augment lab-scale DOE Handbook data for airborne release fractions and respirable fractions in order to provide confidences for safety analysts and facility designers to apply in their analyses at DOE sites. Future work on this topic will involve validation of the presented model against experimental data and extension of the empirical models to be applicable to different classes of particles and surfaces.
Thermal analysts address a wide variety of applications requiring the simulation of radiation heat transfer phenomena. There are gaps in the currently available modeling capabilities. Addressing these gaps would allow for the consideration of additional physics and increase confidence in simulation predictions. This document outlines a five year plan to address the current and future needs of the analyst community with regards to modeling radiation heat transfer processes. This plan represents a significant multi-year effort that must be supported on an ongoing basis.
We present a review and critique of several methods for the simulation of the dynamics of colloidal suspensions at the mesoscale. We focus particularly on simulation techniques for hydrodynamic interactions, including implicit solvents (Fast Lubrication Dynamics, an approximation to Stokesian Dynamics) and explicit/particle-based solvents (Multi-Particle Collision Dynamics and Dissipative Particle Dynamics). Several variants of each method are compared quantitatively for the canonical system of monodisperse hard spheres, with a particular focus on diffusion characteristics, as well as shear rheology and microstructure. In all cases, we attempt to match the relevant properties of a well-characterized solvent, which turns out to be challenging for the explicit solvent models. Reasonable quantitative agreement is observed among all methods, but overall the Fast Lubrication Dynamics technique shows the best accuracy and performance. We also devote significant discussion to the extension of these methods to more complex situations of interest in industrial applications, including models for non-Newtonian solvent rheology, non-spherical particles, drying and curing of solvent and flows in complex geometries. This work identifies research challenges and motivates future efforts to develop techniques for quantitative, predictive simulations of industrially relevant colloidal suspension processes.
Predicting the behavior of solid fuels in response to a fire is a complex endeavor. Heterogeneity, charring, and intumescence are a few examples of the many challenges presented by some common materials. If one desires to employ a 3-dimensional computational fluid dynamics (CFD) model for fire, an accurate solid combustion model for materials at the domain boundary is often desirable. Methods for such modeling are not currently mature, and this is a current topic of research. For some practical problems, it may be acceptable to abstract the surface combustible material as a 1-dimensional reacting boundary condition. This approach has the advantage of being a relatively simple model, and may provide acceptably accurate predictions for problems of interest. Such a model has recently been implemented in Sandia's low-Mach number CFD code for reacting flows, the SIERRA/FUEGO code. Theory for the implemented model is presented. The thermal transport component of the model is verified by approximating a 1-D conduction problem with a closed form solution. The code is further demonstrated by predicting the fire behavior of a block of burning plexiglas (PMMA). The predictions are compared to the reported data from a corresponding experimental program. The predictions are also used to evaluate the sensitivity of model parameters through a sensitivity study using the same test configuration.
The subject of this work is the development of models for the numerical simulation of matter, momentum, and energy balance in heterogeneous materials. These are materials that consist of multiple phases or species or that are structured on some (perhaps many) scale(s). By computational mechanics we mean to refer generally to the standard type of modeling that is done at the level of macroscopic balance laws (mass, momentum, energy). We will refer to the flow or flux of these quantities in a generalized sense as transport. At issue here are the forms of the governing equations in these complex materials which are potentially strongly inhomogeneous below some correlation length scale and are yet homogeneous on larger length scales. The question then becomes one of how to model this behavior and what are the proper multi-scale equations to capture the transport mechanisms across scales. To address this we look to the area of generalized stochastic process that underlie the transport processes in homogeneous materials. The archetypal example being the relationship between a random walk or Brownian motion stochastic processes and the associated Fokker-Planck or diffusion equation. Here we are interested in how this classical setting changes when inhomogeneities or correlations in structure are introduced into the problem. Aspects of non-classical behavior need to be addressed, such as non-Fickian behavior of the mean-squared-displacement (MSD) and non-Gaussian behavior of the underlying probability distribution of jumps. We present an experimental technique and apparatus built to investigate some of these issues. We also discuss diffusive processes in inhomogeneous systems, and the role of the chemical potential in diffusion of hard spheres is considered. Also, the relevance to liquid metal solutions is considered. Finally we present an example of how inhomogeneities in material microstructure introduce fluctuations at the meso-scale for a thermal conduction problem. These fluctuations due to random microstructures also provide a means of characterizing the aleatory uncertainty in material properties at the mesoscale.
This document summarizes a three year Laboratory Directed Research and Development (LDRD) program effort to improve our understanding of algal flocculation with a key to overcoming harvesting as a techno-economic barrier to algal biofuels. Flocculation is limited by the concentrations of deprotonated functional groups on the algal cell surface. Favorable charged groups on the surfaces of precipitates that form in solution and the interaction of both with ions in the water can favor flocculation. Measurements of algae cell-surface functional groups are reported and related to the quantity of flocculant required. Deprotonation of surface groups and complexation of surface groups with ions from the growth media are predicted in the context of PHREEQC. The understanding of surface chemistry is linked to boundaries of effective flocculation. We show that the phase-space of effective flocculation can be expanded by more frequent alga-alga or floc-floc collisions. The collision frequency is dependent on the floc structure, described in the fractal sense. The fractal floc structure is shown to depend on the rate of shear mixing. We present both experimental measurements of the floc structure variation and simulations using LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator). Both show a densification of the flocs with increasing shear. The LAMMPS results show a combined change in the fractal dimension and a change in the coordination number leading to stronger flocs.