Publications

162 Results
Skip to search filters

The kinetics of polyurethane structural foam formation: Foaming and polymerization

AIChE Journal

Rao, Rekha R.; Mondy, L.A.; Long, Kevin N.; Celina, Mathias C.; Wyatt, Nicholas B.; Roberts, Christine C.; Soehnel, Melissa M.; Brunini, Victor B.

Kinetic models have been developed to understand the manufacturing of polymeric foams, which evolve from low viscosity Newtonian liquids, to bubbly liquids, finally producing solid foam. Closed-form kinetics are formulated and parameterized for PMDI-10, a fast curing polyurethane, including polymerization and foaming. PMDI-10 is chemically blown, where water and isocyanate react to form carbon dioxide. The isocyanate reacts with polyol in a competing reaction, producing polymer. Our approach is unique, although it builds on our previous work and the polymerization literature. This kinetic model follows a simplified mathematical formalism that decouples foaming and curing, including an evolving glass transition temperature to represent vitrification. This approach is based on IR, DSC, and volume evolution data, where we observed that the isocyanate is always in excess and does not affect the kinetics. The kinetics are suitable for implementation into a computational fluid dynamics framework, which will be explored in subsequent articles. © 2017 American Institute of Chemical Engineers AIChE J, 63: 2945–2957, 2017.

More Details

Modeling Manufacturing Impacts on Aging and Reliability of Polyurethane Foams

Rao, Rekha R.; Roberts, Christine C.; Mondy, L.A.; Soehnel, Melissa M.; Johnson, Kyle J.; Lorenzo, Henry T.

Polyurethane is a complex multiphase material that evolves from a viscous liquid to a system of percolating bubbles, which are created via a CO2 generating reaction. The continuous phase polymerizes to a solid during the foaming process generating heat. Foams introduced into a mold increase their volume up to tenfold, and the dynamics of the expansion process may lead to voids and will produce gradients in density and degree of polymerization. These inhomogeneities can lead to structural stability issues upon aging. For instance, structural components in weapon systems have been shown to change shape as they age depending on their molding history, which can threaten critical tolerances. The purpose of this project is to develop a Cradle-to-Grave multiphysics model, which allows us to predict the material properties of foam from its birth through aging in the stockpile, where its dimensional stability is important.

More Details

A kinetic approach to modeling the manufacture of high density strucutral foam: Foaming and polymerization

Rao, Rekha R.; Mondy, L.A.; Noble, David R.; Brunini, Victor B.; Roberts, Christine C.; Long, Kevin N.; Soehnel, Melissa M.; Celina, Mathias C.; Wyatt, Nicholas B.; Thompson, Kyle R.

We are studying PMDI polyurethane with a fast catalyst, such that filling and polymerization occur simultaneously. The foam is over-packed to tw ice or more of its free rise density to reach the density of interest. Our approach is to co mbine model development closely with experiments to discover new physics, to parameterize models and to validate the models once they have been developed. The model must be able to repres ent the expansion, filling, curing, and final foam properties. PMDI is chemically blown foam, wh ere carbon dioxide is pr oduced via the reaction of water and isocyanate. The isocyanate also re acts with polyol in a competing reaction, which produces the polymer. A new kinetic model is developed and implemented, which follows a simplified mathematical formalism that decouple s these two reactions. The model predicts the polymerization reaction via condensation chemis try, where vitrification and glass transition temperature evolution must be included to correctly predict this quantity. The foam gas generation kinetics are determined by tracking the molar concentration of both water and carbon dioxide. Understanding the therma l history and loads on the foam due to exothermicity and oven heating is very important to the results, since the kinetics and ma terial properties are all very sensitive to temperature. The conservation eq uations, including the e quations of motion, an energy balance, and thr ee rate equations are solved via a stabilized finite element method. We assume generalized-Newtonian rheology that is dependent on the cure, gas fraction, and temperature. The conservation equations are comb ined with a level set method to determine the location of the free surface over time. Results from the model are compared to experimental flow visualization data and post-te st CT data for the density. Seve ral geometries are investigated including a mock encapsulation part, two configur ations of a mock stru ctural part, and a bar geometry to specifically test the density model. We have found that the model predicts both average density and filling profiles well. However, it under predicts density gradients, especially in the gravity direction. Thoughts on m odel improvements are also discussed.

More Details

Multilayer Coextrusion of Polymer Composites to Develop Organic Capacitors

International Polymer Processing

Mondy, L.A.; Bieg, Lothar F.; Spangler, Scott W.; Stavig, Mark E.; Schroeder, John L.; Rao, Rekha R.; DiAntonio, Christopher D.

Multilayer coextrusion is applied to produce a tape containing layers of alternating electrical properties to demonstrate the potential for using coextrusion to manufacture capacitors. To obtain the desired properties, we develop two filled polymer systems, one for conductive layers and one for dielectric layers. We describe numerical models used to help determine the material and processing parameters that impact processing and layer stability. These models help quantify the critical ratios of densities and viscosities of the two layers to maintain stable layers, as well as the effect of increasing the flow rate of one of the two materials. The conducting polymer is based on polystyrene filled with a blend of low-melting-point eutectic metal and nickel particulate filler, as described by Mrozek et al. (2010). The appropriate concentrations of fillers are determined by balancing measured conductivity with processability in a twin screw extruder. Based on results of the numerical models and estimates of the viscosity of emulsions and suspensions, a dielectric layer composed of polystyrene filled with barium titanate is formulated. Despite the fact that the density of the dielectric filler is less than the metallic filler of the conductive phase, as well as rheological measurements that later showed that the dielectric formulation is not an ideal match to the viscosity of the conductive material, the two materials can be successfully coextruded if the flow rates of the two materials are not identical. A measurable capacitance of the layered structure is obtained.

More Details

The importance of experimental design on measurement of dynamic interfacial tension and interfacial rheology in diffusion-limited surfactant systems

Colloids and Surfaces A: Physicochemical and Engineering Aspects

Reichert, Matthew D.; Alvarez, Nicolas J.; Brooks, Carlton F.; Grillet, Anne M.; Mondy, L.A.; Anna, Shelley L.; Walker, Lynn M.

Pendant bubble and drop devices are invaluable tools in understanding surfactant behavior at fluid-fluid interfaces. The simple instrumentation and analysis are used widely to determine adsorption isotherms, transport parameters, and interfacial rheology. However, much of the analysis performed is developed for planar interfaces. The application of a planar analysis to drops and bubbles (curved interfaces) can lead to erroneous and unphysical results. We revisit this analysis for a well-studied surfactant system at air-water interfaces over a wide range of curvatures as applied to both expansion/contraction experiments and interfacial elasticity measurements. The impact of curvature and transport on measured properties is quantified and compared to other scaling relationships in the literature. The results provide tools to design interfacial experiments for accurate determination of isotherm, transport and elastic properties.

More Details

Experiments to populate and validate a processing model for polyurethane foam. BKC 44306 PMDI-10

Mondy, L.A.; Bauer, Stephen J.; Hileman, Michael B.; Thompson, Kyle R.; Rao, Rekha R.; Shelden, Bion S.; Soehnel, Melissa M.; O'Hern, Timothy J.; Grillet, Anne M.; Celina, Mathias C.; Wyatt, Nicholas B.; Russick, Edward M.

We are developing computational models to elucidate the expansion and dynamic filling process of a polyurethane foam, PMDI. The polyurethane of interest is chemically blown, where carbon dioxide is produced via the reaction of water, the blowing agent, and isocyanate. The isocyanate also reacts with polyol in a competing reaction, which produces the polymer. Here we detail the experiments needed to populate a processing model and provide parameters for the model based on these experiments. The model entails solving the conservation equations, including the equations of motion, an energy balance, and two rate equations for the polymerization and foaming reactions, following a simplified mathematical formalism that decouples these two reactions. Parameters for the polymerization kinetics model are reported based on infrared spectrophotometry. Parameters describing the gas generating reaction are reported based on measurements of volume, temperature and pressure evolution with time. A foam rheology model is proposed and parameters determined through steady-shear and oscillatory tests. Heat of reaction and heat capacity are determined through differential scanning calorimetry. Thermal conductivity of the foam as a function of density is measured using a transient method based on the theory of the transient plane source technique. Finally, density variations of the resulting solid foam in several simple geometries are directly measured by sectioning and sampling mass, as well as through x-ray computed tomography. These density measurements will be useful for model validation once the complete model is implemented in an engineering code.

More Details

Thermophysical properties of BKC 44306 and BKC 44307 PMDI urethane solid and foams

Bauer, Stephen J.; Flint, Gregory M.; Mondy, L.A.

Accurate knowledge of thermophysical properties of urethane foam is considered extremely important for meaningful models and analyses to be developed of scenarios wherein the foam is heated. Its performance at temperature requires a solid understanding of the foam material properties at temperature. Also, foam properties vary with density/porosity. An experimental program to determine the thermal properties of the two foams and their parent solid urethane was developed in order to support development of a predictive model relating density and thermal properties from first principles. Thermal properties (thermal conductivity, diffusivity, and specific heat) of the foam were found to vary with temperatures from 26°C to 90°C. Thermal conductivity generally increases with increasing temperature for a given initial density and ranges from .0433 W/mK at 26°C to .0811 W/mK at 90°C; thermal diffusivity generally decreases with increasing temperature for a given initial density and ranges from .4101 mm2/s at 26°C to .1263 mm2/s at 90°C; and specific heat generally increases with increasing temperature for a given initial density and ranges from .1078 MJ/m3K at 26°C to .6323 MJ/m3K at 90°C. Thermal properties of the solid urethane were also found to vary with temperatures from 26°C to 90°C. Average thermal conductivity generally increases with increasing temperature for a given initial density and ranges from 0.126 to 0.131 W/mK at 26°C to 0.153 to 0.157 W/mK at 90°C; average thermal diffusivity generally decreases with increasing temperature for a given initial density and ranges from 0.142 to 0.147 mm2/s at 26°C to 0.124 to 0.125 mm2/s at 90°C; and average specific heat generally increases with increasing temperature for a given initial density and ranges from 0.889 to 0.899 MJ/m3K to 1.229 to 1.274 MJ/m3K at 90°C. The density of both foam and solid urethane decreased with increasing temperature.

More Details

New composite separator pellet to increase power density and reduce size of thermal batteries

Mondy, L.A.; Evans, Lindsey E.; Roberts, Christine C.; Grillet, Anne M.; Soehnel, Melissa M.; Barringer, David A.; DiAntonio, Christopher D.; Chavez, Tom C.; Ingersoll, David I.; Hughes, Lindsey G.

We show that it is possible to manufacture strong macroporous ceramic films that can be backfilled with electrolyte to form rigid separator pellets suitable for use in thermal batteries. Several new ceramic manufacturing processes are developed to produce sintered magnesium oxide foams with connected porosities of over 80% by volume and with sufficient strength to withstand the battery manufacturing steps. The effects of processing parameters are quantified, and methods to imbibe electrolyte into the ceramic scaffold demonstrated. Preliminary single cell battery testing show that some of our first generation pellets exhibit longer voltage life with comparable resistance at the critical early times to that exhibited by a traditional pressed pellets. Although more development work is needed to optimize the processes to create these rigid separator pellets, the results indicate the potential of such ceramic separator pellets to be equal, if not superior to, current pressed pellets. Furthermore, they could be a replacement for critical material that is no longer available, as well as improving battery separator strength, decreasing production costs, and leading to shorter battery stacks for long-life batteries.

More Details

Computational Mechanics for Heterogeneous Materials

Baczewski, Andrew D.; Yarrington, Cole Y.; Bond, Stephen D.; Erikson, William W.; Lehoucq, Richard B.; Mondy, L.A.; Noble, David R.; Pierce, Flint P.; Roberts, Christine C.; Van Swol, Frank

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.

More Details

Size and structure of Chlorella zofingiensis /FeCl 3 flocs in a shear flow: Algae Floc Structure

Biotechnology and Bioengineering

Wyatt, Nicholas B.; O'Hern, Timothy J.; Shelden, Bion S.; Hughes, Lindsey G.; Mondy, L.A.

Flocculation is a promising method to overcome the economic hurdle to separation of algae from its growth medium in large scale operations. But, understanding of the floc structure and the effects of shear on the floc structure are crucial to the large scale implementation of this technique. The floc structure is important because it determines, in large part, the density and settling behavior of the algae. Freshwater algae floc size distributions and fractal dimensions are presented as a function of applied shear rate in a Couette cell using ferric chloride as a flocculant. Comparisons are made with measurements made for a polystyrene microparticle model system taken here as well as reported literature results. The algae floc size distributions are found to be self-preserving with respect to shear rate, consistent with literature data for polystyrene. Moreover, three fractal dimensions are calculated which quantitatively characterize the complexity of the floc structure. Low shear rates result in large, relatively dense packed flocs which elongate and fracture as the shear rate is increased. Our results presented here provide crucial information for economically implementing flocculation as a large scale algae harvesting strategy.

More Details

First-principles flocculation as the key to low energy algal biofuels processing

Hewson, John C.; Mondy, L.A.; Murton, Jaclyn K.; O'Hern, Timothy J.; Parchert, Kylea J.; Pohl, Phillip I.; Williams, Cecelia V.; Wyatt, Nicholas B.; Barringer, David A.; Pierce, Flint P.; Brady, Patrick V.; Dwyer, Brian P.; Grillet, Anne M.; Hankins, M.G.; Hughes, Lindsey G.; Lechman, Jeremy B.

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.

More Details

Finite element analysis of multilayer coextrusion

Rao, Rekha R.; Mondy, L.A.; Schunk, Randy; Hopkins, Matthew M.

Multilayer coextrusion has become a popular commercial process for producing complex polymeric products from soda bottles to reflective coatings. A numerical model of a multilayer coextrusion process is developed based on a finite element discretization and two different free-surface methods, an arbitrary-Lagrangian-Eulerian (ALE) moving mesh implementation and an Eulerian level set method, to understand the moving boundary problem associated with the polymer-polymer interface. The goal of this work is to have a numerical capability suitable for optimizing and troubleshooting the coextrusion process, circumventing flow instabilities such as ribbing and barring, and reducing variability in layer thickness. Though these instabilities can be both viscous and elastic in nature, for this work a generalized Newtonian description of the fluid is used. Models of varying degrees of complexity are investigated including stability analysis and direct three-dimensional finite element free surface approaches. The results of this work show how critical modeling can be to reduce build test cycles, improve material choices, and guide mold design.

More Details

Creation and characterization of magnesium oxide macroporous ceramics

Materials Engineering and Sciences Division - Core Programming Topic at the 2011 AIChE Annual Meeting

Mondy, L.A.; DiAntonio, Christopher D.; Chavez, Tom C.; Hughes, Lindsey G.; Grillet, Anne M.; Roberts, Christine C.; Ingersoll, David I.

We examine several methods to create a sheet of magnesium oxide (MgO) macroporous ceramic material via tape casting. These methods include the approach pioneered by Akartuna et al.1 in which an oil/water emulsion is stabilized by surface-modified metal oxide particles at the droplet interfaces. Upon drying, a scaffold of the self-assembled particles is strong enough to be removed from the substrate material and sintered. We find that this method can be used with MgO particles surface modified by short amphiphilic molecules. This approach is compared with two more traditional methods to induce structure into a green ceramic: 1) creation of an MgO ceramic slip with added pore formers, and 2) sponge impregnation of a reticulated foam with the MgO slip. Green and sintered samples made using each method are hardness tested and results compared for several densities of the final ceramics. Optical and SEM images of the materials are shown.

More Details

Surface rheology and interface stability

Mondy, L.A.; Brooks, Carlton F.; Grillet, Anne M.; Moffat, Harry K.; Koehler, Timothy P.; Yaklin, Melissa A.; Cote, Raymond O.; Castaneda, Jaime N.

We have developed a mature laboratory at Sandia to measure interfacial rheology, using a combination of home-built, commercially available, and customized commercial tools. An Interfacial Shear Rheometer (KSV ISR-400) was modified and the software improved to increase sensitivity and reliability. Another shear rheometer, a TA Instruments AR-G2, was equipped with a du Nouey ring, bicone geometry, and a double wall ring. These interfacial attachments were compared to each other and to the ISR. The best results with the AR-G2 were obtained with the du Nouey ring. A Micro-Interfacial Rheometer (MIR) was developed in house to obtain the much higher sensitivity given by a smaller probe. However, it was found to be difficult to apply this technique for highly elastic surfaces. Interfaces also exhibit dilatational rheology when the interface changes area, such as occurs when bubbles grow or shrink. To measure this rheological response we developed a Surface Dilatational Rheometer (SDR), in which changes in surface tension with surface area are measured during the oscillation of the volume of a pendant drop or bubble. All instruments were tested with various surfactant solutions to determine the limitations of each. In addition, foaming capability and foam stability were tested and compared with the rheology data. It was found that there was no clear correlation of surface rheology with foaming/defoaming with different types of surfactants, but, within a family of surfactants, rheology could predict the foam stability. Diffusion of surfactants to the interface and the behavior of polyelectrolytes were two subjects studied with the new equipment. Finally, surface rheological terms were added to a finite element Navier-Stokes solver and preliminary testing of the code completed. Recommendations for improved implementation were given. When completed we plan to use the computations to better interpret the experimental data and account for the effects of the underlying bulk fluid.

More Details

Investigation of interfacial rheology & foam stability

Grillet, Anne M.; Koehler, Timothy P.; Reichert, Matthew D.; Yaklin, Melissa A.; Brooks, Carlton F.; Mondy, L.A.; Castaneda, Jaime N.; Cote, Raymond O.

The rheology at gas-liquid interfaces strongly influences the stability and dynamics of foams and emulsions. Several experimental techniques are employed to characterize the rheology at liquid-gas interfaces with an emphasis on the non-Newtonian behavior of surfactant-laden interfaces. The focus is to relate the interfacial rheology to the foamability and foam stability of various aqueous systems. An interfacial stress rheometer (ISR) is used to measure the steady and dynamic rheology by applying an external magnetic field to actuate a magnetic needle suspended at the interface. Results are compared with those from a double wall ring attachment to a rotational rheometer (TA Instruments AR-G2). Micro-interfacial rheology (MIR) is also performed using optical tweezers to manipulate suspended microparticle probes at the interface to investigate the steady and dynamic rheology. Additionally, a surface dilatational rheometer (SDR) is used to periodically oscillate the volume of a pendant drop or buoyant bubble. Applying the Young-Laplace equation to the drop shape, a time-dependent surface tension can be calculated and used to determine the effective dilatational viscosity of an interface. Using the ISR, double wall ring, SDR, and MIR, a wide range of sensitivity in surface forces (fN to nN) can be explored as each experimental method has different sensitivities. Measurements will be compared to foam stability.

More Details

Liquid organic foams for formulation optimization : an assessment of foam linear viscoelasticity and its temporal dependence

Kropka, Jamie M.; Celina, Mathias C.; Mondy, L.A.

Liquid foams are viscoelastic liquids, exhibiting a fast relaxation attributed to local bubble motions and a slow response due to structural evolution of the intrinsically unstable system. In this work, these processes are examined in unique organic foams that differ from the typically investigated aqueous systems in two major ways: the organic foams (1) posses a much higher continuous phase viscosity and (2) exhibit a coarsening response that involves coalescence of cells. The transient and dynamic relaxation responses of the organic foams are evaluated and discussed in relation to the response of aqueous foams. The change in the foam response with increasing gas fraction, from that of a Newtonian liquid to one that is strongly viscoelastic, is also presented. In addition, the temporal dependencies of the linear viscoelastic response are assessed in the context of the foam structural evolution. These foams and characterization techniques provide a basis for testing stabilization mechanisms in epoxy-based foams for encapsulation applications.

More Details

Multilayer co-extrusion technique for developing high energy density organic devices

Mondy, L.A.; Rao, Rekha R.; Bieg, Lothar F.; Schneider, Duane A.; Stavig, Mark E.; Schroeder, John L.; Winter, Michael R.

The purpose of this project is to develop multi-layered co-extrusion (MLCE) capabilities at Sandia National Laboratories to produce multifunctional polymeric structures. Multi-layered structures containing layers of alternating electrical, mechanical, optical, or structural properties can be applied to a variety of potential applications including energy storage, optics, sensors, mechanical, and barrier applications relevant to the internal and external community. To obtain the desired properties, fillers must be added to the polymer materials that are much smaller than the end layer thickness. We developed two filled polymer systems, one for conductive layers and one for dielectric layers and demonstrated the potential for using MLCE to manufacture capacitors. We also developed numerical models to help determine the material and processing parameters that impact processing and layer stability.

More Details

Measurements of wall slip during rise of a physically blown foam

AIP Conference Proceedings

Bourdon, Christopher B.; Bourdon, Christopher B.; Grillet, Anne M.; Mondy, L.A.; Rao, Rekha R.

Polymeric foam systems are widely used in industrial applications due to their low weight and abilities to thermally insulate and isolate vibration. However, processing of these foams is still not well understood at a fundamental level. The precursor foam of interest starts off as a liquid phase emulsion of blowing agent in a thermosetting polymer. As the material is heated either by an external oven or by the exothermic reaction from internal polymerization of the suspending fluid, the blowing agent boils to produce gas bubbles and a foamy material. A series of experiments have been performed to allow observation of the foaming process and the collection of temperature, rise rate, and microstructural data. Microfocus video is used in conjunction with particle image velocimetry (PIV) to elucidate the boundary condition at the wall. These data provide input to a continuum level finite element model of the blowing process. PIV is used to measure the slip velocity of foams with a volume fraction range of 0.50 to 0.71. These results are in agreement with theoretical predictions which suggest that at high volume fractions the bubbles would exhibit jamming behavior and slip at the wall. At these volume fractions, the slip velocity profile has a shear profile shape near the side walls and a plug flow shape at the center. The shape of the velocity profile is in agreement with previous experimental work investigating different foam systems. As time increases, the available blowing agent decreases, the volume fraction increases, the viscosity increases, and the average slip velocity decreases, but the slip velocity profile maintains the plug-shear shape. © 2008 American Institute of Physics.

More Details

Foam process models

Rao, Rekha R.; Mondy, L.A.; Moffat, Harry K.; Noble, David R.; Notz, Patrick N.; Adolf, Douglas B.

In this report, we summarize our work on developing a production level foam processing computational model suitable for predicting the self-expansion of foam in complex geometries. The model is based on a finite element representation of the equations of motion, with the movement of the free surface represented using the level set method, and has been implemented in SIERRA/ARIA. An empirically based time- and temperature-dependent density model is used to encapsulate the complex physics of foam nucleation and growth in a numerically tractable model. The change in density with time is at the heart of the foam self-expansion as it creates the motion of the foam. This continuum-level model uses an homogenized description of foam, which does not include the gas explicitly. Results from the model are compared to temperature-instrumented flow visualization experiments giving the location of the foam front as a function of time for our EFAR model system.

More Details

Experiments for foam model development and validation

Mondy, L.A.; Gorby, Allen D.; Cote, Raymond O.; Castaneda, Jaime N.; Thompson, Kyle R.; Rao, Rekha R.; Moffat, Harry K.; Kraynik, Andrew M.; Russick, Edward M.; Adolf, Douglas B.; Grillet, Anne M.; Brotherton, Christopher M.; Bourdon, Christopher B.

A series of experiments has been performed to allow observation of the foaming process and the collection of temperature, rise rate, and microstructural data. Microfocus video is used in conjunction with particle image velocimetry (PIV) to elucidate the boundary condition at the wall. Rheology, reaction kinetics and density measurements complement the flow visualization. X-ray computed tomography (CT) is used to examine the cured foams to determine density gradients. These data provide input to a continuum level finite element model of the blowing process.

More Details

Pressure-driven and free-rise foam flow

Mondy, L.A.; Kropka, Jamie M.; Celina, Mathias C.; Rao, Rekha R.; Brotherton, Christopher M.; Bourdon, Christopher B.; Noble, David R.; Moffat, Harry K.; Grillet, Anne M.; Kraynik, Andrew M.; Leming, Sarah L.

Many weapons components (e.g. firing sets) are encapsulated with blown foams. Foam is a strong lightweight material--good compromise between conflicting needs of structural stability and electronic function. Current foaming processes can lead to unacceptable voids, property variations, cracking, and slipped schedules which is a long-standing issue. Predicting the process is not currently possible because the material is polymerizing and multiphase with changing microstructure. The goals of this project is: (1) Produce uniform encapsulant consistently and improve processability; (2) Eliminate metering issues/voids; (3) Lower residual stresses, exotherm to protect electronics; and (4) Maintain desired properties--lightweight, strong, no delamination/cracking, and ease of removal. The summary of achievements in the first year are: (1) Developed patentable chemical foaming chemistry - TA; (2) Developed persistent non-curing foam for systematic evaluation of fundamental physics of foams--Initial testing of non-curing foam shows that surfactants very important; (3) Identified foam stability strategy using a stacked reaction scheme; (4) Developed foam rheology methodologies and shear apparatuses--Began testing candidates for shear stability; (5) Began development of computational model; and (6) Development of methodology and collection of property measurements/boundary conditions for input to computational model.

More Details

Wetting and free surface flow modeling for potting and encapsulation

Adolf, Douglas B.; Castaneda, Jaime N.; Kraynik, Andrew M.; Noble, David R.; Sun, Amy C.; Cote, Raymond O.; Grillet, Anne M.; Notz, Patrick N.; Brooks, Carlton F.; Givler, R.C.; Hopkins, Matthew M.; Mondy, L.A.; Rao, Rekha R.

As part of an effort to reduce costs and improve quality control in encapsulation and potting processes the Technology Initiative Project ''Defect Free Manufacturing and Assembly'' has completed a computational modeling study of flows representative of those seen in these processes. Flow solutions are obtained using a coupled, finite-element-based, numerical method based on the GOMA/ARIA suite of Sandia flow solvers. The evolution of the free surface is solved with an advanced level set algorithm. This approach incorporates novel methods for representing surface tension and wetting forces that affect the evolution of the free surface. In addition, two commercially available codes, ProCAST and MOLDFLOW, are also used on geometries representing encapsulation processes at the Kansas City Plant. Visual observations of the flow in several geometries are recorded in the laboratory and compared to the models. Wetting properties for the materials in these experiments are measured using a unique flowthrough goniometer.

More Details

Modeling injection molding of net-shape active ceramic components

Rao, Rekha R.; Brooks, Carlton F.; Cote, Raymond O.; Castaneda, Jaime N.; Mondy, L.A.; Noble, David R.; Hopkins, Matthew M.; Notz, Patrick N.; Halbleib, Laura L.; Yang, Pin Y.; Burns, George B.; Grillet, Anne M.

To reduce costs and hazardous wastes associated with the production of lead-based active ceramic components, an injection molding process is being investigated to replace the current machining process. Here, lead zirconate titanate (PZT) ceramic particles are suspended in a thermoplastic resin and are injected into a mold and allowed to cool. The part is then bisque fired and sintered to complete the densification process. To help design this new process we use a finite element model to describe the injection molding of the ceramic paste. Flow solutions are obtained using a coupled, finite-element based, Newton-Raphson numerical method based on the GOMA/ARIA suite of Sandia flow solvers. The evolution of the free surface is solved with an advanced level set algorithm. This approach incorporates novel methods for representing surface tension and wetting forces that affect the evolution of the free surface. Thermal, rheological, and wetting properties of the PZT paste are measured for use as input to the model. The viscosity of the PZT is highly dependent both on temperature and shear rate. One challenge in modeling the injection process is coming up with appropriate constitutive equations that capture relevant phenomenology without being too computationally complex. For this reason we model the material as a Carreau fluid and a WLF temperature dependence. Two-dimensional (2D) modeling is performed to explore the effects of the shear in isothermal conditions. Results indicate that very low viscosity regions exist near walls and that these results look similar in terms of meniscus shape and fill times to a simple Newtonian constitutive equation at the shear-thinned viscosity for the paste. These results allow us to pick a representative viscosity to use in fully three-dimensional (3D) simulation, which because of numerical complexities are restricted to using a Newtonian constitutive equation. Further 2D modeling at nonisothermal conditions shows that the choice of representative Newtonian viscosity is dependent on the amount of heating of the initially room temperature mold. An early 3D transient model shows that the initial design of the distributor is sub-optimal. However, these simulations take several months to run on 4 processors of an HP workstation using a preconditioner/solver combination of ILUT/GMRES with fill factors of 3 and PSPG stabilization. Therefore, several modifications to the distributor geometry and orientations of the vents and molds have been investigated using much faster 3D steady-state simulations. The pressure distribution for these steady-state calculations is examined for three different distributor designs to see if this can indicate which geometry has the superior design. The second modification, with a longer distributor, is shown to have flatter, more monotonic isobars perpendicular to the flow direction indicating a better filling process. The effects of the distributor modifications, as well as effects of the mold orientation, have also been examined with laboratory experiments in which the flow of a viscous Newtonian oil entering transparent molds is recorded visually. Here, the flow front is flatter and voids are reduced for the second geometry compared to the original geometry. A horizontal orientation, as opposed to the planned vertical orientation, results in fewer voids. Recently, the Navier-Stokes equations have been stabilized with the Dohrman-Bochev PSPP stabilization method, allowing us to calculate transient 3D simulations with computational times on the order of days instead of months. Validation simulations are performed and compared to the experiments. Many of the trends of the experiments are captured by the level set modeling, though quantitative agreement is lacking mainly due to the high value of the gas phase viscosity necessary for numerical stability, though physically unrealistic. More correct trends are predicted for the vertical model than the horizontal model, which is serendipitous as the actual mold is held in a vertical geometry. The full, transient mold filling calculations indicate that the flow front is flatter and voids may be reduced for the second geometry compared to the original geometry. The validated model is used to predict mold filling for the actual process with the material properties for the PZT paste, the original distributor geometry, and the mold in a vertical orientation. This calculation shows that voids may be trapped at the four corners of the mold opposite the distributor.

More Details

Maximizing storage rates and capacity of carbon dioxide sequestration in saline reservoirs

WIT Transactions on Ecology and the Environment

Abou-Sayed, A.; Guo, Q.; Graham, A.L.; Mondy, L.A.; Ingber, M.S.; Mammoli, A.A.

The Kyoto Accords have been signed by 140 nations in order to significantly reduce carbon dioxide emissions into the atmosphere in the medium to long term. In order to achieve this goal without drastic reductions in fossil fuel usage, carbon dioxide must be removed from the atmosphere and stored in acceptable reservoirs. Research has been undertaken to develop economical new technologies for the transfer and storage of carbon dioxide in saline aquifers. In order to maximize the storage rate, the aquifer is first hydraulically fractured in a conventional well stimulation treatment with a slurry containing solid proppant. Well fracturing would increase the injection volume flowrate greatly. In addition, there are several ancillary benefits including extension of the reservoir early storage volume by moving the carbon dioxide further from the well. This extended reach would mitigate the problems with the buoyant plume and increase the surface area between the carbon dioxide and the formation facilitating absorption. A life-cycle cost estimate has been performed showing the benefits of this approach compared to injection without fracturing.

More Details

Mold-filling experiments for validation of modeling encapsulation. Part 1, "wine glass" mold

Mondy, L.A.; Grillet, Anne M.; Cote, Raymond O.; Castaneda, Jaime N.

The C6 project 'Encapsulation Processes' has been designed to obtain experimental measurements for discovery of phenomena critical to improving these processes, as well as data required in the verification and validation plan (Rao et al. 2001) for model validation of flow in progressively complex geometries. We have observed and recorded the flow of clear, Newtonian liquids and opaque, rheologically complex suspensions in two mold geometries. The first geometry is a simple wineglass geometry in a cylinder and is reported here in Part 1. The results in a more realistic encapsulation geometry are reported in Part 2.

More Details

The effects of process parameters on injection-molded PZT ceramics part fabrication- compounding process rheology

Yang, Pin Y.; Mondy, L.A.; Burns, George B.

Solid solutions of lead-based perovskites are the backbone materials of the piezoelectric components for transducer, actuator, and resonator applications. These components, typically small in size, are fabricated from large sintered ceramic slugs using grinding and lapping processes. These operations increase manufacturing costs and produce a large hazardous waste stream, especially when component size decreases. To reduce costs and hazardous wastes associated with the production of these components, an injection molding technique is being investigated to replace the machining processes. The first step in the new technique is to compound an organic carrier with a ceramic powder. The organic carrier is a thermoplastic based system composed of a main carrier, a binder, and a surfactant. Understanding the rheology of the compounded material is necessary to minimize the creation of defects such as voids or cavities during the injection-molding process. An experiment was performed to model the effects of changes in the composition and processing of the material on the rheological behavior. Factors studied included: the surfactant of the organic carrier system, the solid loading of the compounded material, and compounding time. The effects of these factors on the viscosity of the material were investigated.

More Details

Particle migration rates in a Couette apparatus

Proposed for publication in Journal of Mathematics.

Mondy, L.A.

Bulk migration of particles towards regions of lower shear occurs in suspensions of neutrally buoyant spheres in Newtonian fluids undergoing creeping flow in the annular region between two rotating, coaxial cylinders (a wide-gap Couette). For a monomodal suspension of spheres in a viscous fluid, dimensional analysis indicates that the rate of migration at a given concentration should scale with the square of the sphere radius. However, a previous experimental study showed that the rate of migration of spherical particles at 50% volume concentration actually scaled with the sphere radius to approximately the 2.9 power.

More Details

Particle migration rates in a Couette apparatus

Advances in Fluid Mechanics

Hsiao, S.C.; Christensen, D.; Ingber, Marc S.; Mondy, L.A.; Altobelli, S.A.

Bulk migration of particles towards regions of lower shear occurs in suspensions of neutrally buoyant spheres in Newtonian fluids undergoing creeping flow in the annular region between two rotating, coaxial cylinders (a wide-gap Couette). For a monomodal suspension of spheres in a viscous fluid, dimensional analysis indicates that the rate of migration at a given concentration should scale with the square of the sphere radius. However, a previous experimental study [12] showed that the rate of migration of spherical particles at 50% volume concentration actually scaled with the sphere radius to approximately the 2.9 power. In the current study, a series of experiments is performed to extend the previous study to two new concentrations, namely, 35% and 42.5%. The new study indicates that the power scaling decreases with concentration.

More Details

Three-dimensional boundary element simulations of concentrated suspensions in a spinning ball rheometer

Proposed for publication in Communications in Numerical Methods in Engineering.

Grillet, Anne M.; Mondy, L.A.

The spinning ball rheometer has been proposed as a method to measure rheological properties of concentrated suspensions. Recent experiments have shown that the measured extra torque on the spinning ball decreases as the radius of the spinning ball becomes comparable to the size of the suspended particle. We have performed a series of three dimensional boundary element calculations of the rheometer geometry to probe the microstructure effects that contribute to the apparent 'slip.' We will present a series of snap-shot results as well as several transient calculations which are compared to the available experimental data. The computational limitations of these large-scale simulations shall also be discussed.

More Details

Apparent slip at the surface of a ball spinning in a concentrated suspension

Proposed for publication in the Journal of Fluid Mechanics.

Mondy, L.A.; Grillet, Anne M.; Henfling, John F.

The couple on a ball rotating relative to an otherwise quiescent suspension of comparably-sized, neutrally buoyant spheres is studied both experimentally and numerically. Apparent 'slip' relative to the analytical solution for a sphere spinning in a Newtonian fluid (based upon the viscosity of the suspension) is determined in suspensions with volume fractions c ranging from 0.03 to 0.50. This apparent slip results in a decrease of the measured torque on the spinning ball when the radius of the ball becomes comparable with that of the suspended spheres. Over the range of our data, the slip becomes more pronounced as the concentration c increases. At c = 0.25, three-dimensional boundary-element simulations agree well with the experimental data. Moreover, at c = 0.03, good agreement exists between such calculations and theoretical predictions of rotary slip in dilute suspensions.

More Details

Chemically selective NMR imaging of a 3-component (solid-solid-liquid) sedimenting system

Journal of Magnetic Resonance

Beyea, Steven D.; Altobelli, Stephen A.; Mondy, L.A.

A novel magnetic resonance imaging (MRI) technique which resolves the separate components of the evolving vertical concentration profiles of 3-component non-colloidal suspensions is described. This method exploits the sensitivity of MRI to chemical differences between the three phases to directly image the fluid phase and one of the solid phases, with the third phase obtained by subtraction. 19F spin-echo imaging of a polytetrafluoroethylene (PTFE) oil was interlaced with 1H SPRITE imaging of low-density polyethylene (LDPE) particles. The third phase was comprised of borosilicate glass spheres, which were not visible while imaging the PTFE or LDPE phases. The method is demonstrated by performing measurements on 2-phase materials containing only the floating (LDPE) particles, with the results contrasted to the experimental behaviour of the individual phases in the full 3-phase system. All experiments were performed using nearly monodisperse particles, with initial suspension volume fractions, φi, of 0.1. © 2003 Elsevier Science (USA). All rights reserved.

More Details

NMR measurements and simulations of particle migration in non-Newtonian fluids

Chemical Engineering Communications

Rao, Rekha R.; Mondy, L.A.; Baer, Thomas A.; Altobelli, Stephen A.; Stephens, Thomas S.

Shear-induced migration of particles is studied during the slow flow of suspensions of neutrally buoyant spheres, at 50% particle volume fraction, in an inelastic but shear-thinning, suspending fluid. The suspension is flowing in between a rotating inner cylinder and a stationary outer cylinder. The conditions are such that nonhydrodynamic effects are negligible. Nuclear magnetic resonance (NMR) imaging demonstrates that the movement of particles away from the high shear rate region is more pronounced than for a Newtonian suspending liquid. We test a continuum constitutive model for the evolution of particle concentration in a flowing suspension proposed by Phillips et al., but extended to shear-thinning, suspending fluids. The fluid constitutive equation is Carreau-like in its shear-thinning behavior but also varies with the local particle concentration. The model captures many of the trends found in the experimental data, but does not yet agree quantitatively. In fact, quantitative agreement with a diffusive flux constitutive equation would be impossible without the addition of another fitting parameter that may depend on the shear-thinning nature of the suspending fluid. Because of this, we feel that the Phillips model may be fundamentally inadequate for simulating flows of particles in non-Newtonian suspending fluids without the introduction of a normal stress correction or other augmenting terms.

More Details

Verification and Validation of Encapsulation Flow Models in GOMA, Version 1.1

Mondy, L.A.; Rao, Rekha R.; Schunk, Randy; Sackinger, Philip A.; Adolf, Douglas B.; Rao, Rekha R.

Encapsulation is a common process used in manufacturing most non-nuclear components including: firing sets, neutron generators, trajectory sensing signal generators (TSSGs), arming, fusing and firing devices (AF and Fs), radars, programmers, connectors, and batteries. Encapsulation is used to contain high voltage, to mitigate stress and vibration and to protect against moisture. The purpose of the ASCI Encapsulation project is to develop a simulation capability that will allow us to aid in the encapsulation design process, especially for neutron generators. The introduction of an encapsulant poses many problems because of the need to balance ease of processing and properties necessary to achieve the design benefits such as tailored encapsulant properties, optimized cure schedule and reduced failure rates. Encapsulants can fail through fracture or delamination as a result of cure shrinkage, thermally induced residual stresses, voids or incomplete component embedding and particle gradients. Manufacturing design requirements include (1) maintaining uniform composition of particles in order to maintain the desired thermal coefficient of expansion (CTE) and density, (2) mitigating void formation during mold fill, (3) mitigating cure and thermally induced stresses during cure and cool down, and (4) eliminating delamination and fracture due to cure shrinkage/thermal strains. The first two require modeling of the fluid phase, and it is proposed to use the finite element code GOMA to accomplish this. The latter two require modeling of the solid state; however, ideally the effects of particle distribution would be included in the calculations, and thus initial conditions would be set from GOMA predictions. These models, once they are verified and validated, will be transitioned into the SIERRA framework and the ARIA code. This will facilitate exchange of data with the solid mechanics calculations in SIERRA/ADAGIO.

More Details

Massively parallel boundary integral element method modeling of particles in low-Reynolds number Newtonian fluid flows

Advances in High Performance Computing

Ingber, M.S.; Subia, S.R.; Mondy, L.A.

The analysis of many complex multiphase fluid flow systems is based on a scale decoupling procedure. At the macroscale, continuum models are used to perform large scale simulations. At the mesoscale, statistical homogenization theory is used to derive continuum models based on representative volume elements (RVEs). At the microscale, small scale features such as interfacial properties are analyzed to be incorporated into mesoscale simulations. In this research, mesoscopic simulations of hard particles suspended in a Newtonian fluid undergoing nonlinear shear flow are performed using a boundary element method. To obtain an RVE at higher concentrations, several hundred particles are included in the simulations putting considerable demands on the computational resources both in terms of CPU and memory. Parallel computing provides a viable platform to study these large multiphase systems. The implementation of a portable, parallel computer code based on the boundary element method using a block-block data distribution is discussed in this paper. The code employs updated direct-solver technologies that make use of dual-processor compute nodes.

More Details
162 Results
162 Results