Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

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.

Abstract not provided.

Abstract not provided.

We are developing computational models to help understand manufacturing processes, final properties and aging of structural foam, polyurethane PMDI. Th e resulting model predictions of density and cure gradients from the manufacturing process will be used as input to foam heat transfer and mechanical models. BKC 44306 PMDI-10 and BKC 44307 PMDI-18 are the most prevalent foams used in structural parts. Experiments needed to parameterize models of the reaction kinetics and the equations of motion during the foam blowing stages were described for BKC 44306 PMDI-10 in the first of this report series (Mondy et al. 2014). BKC 44307 PMDI-18 is a new foam that will be used to make relatively dense structural supports via over packing. It uses a different catalyst than those in the BKC 44306 family of foams; hence, we expect that the reaction kineti cs models must be modified. Here we detail the experiments needed to characteriz e the reaction kinetics of BKC 44307 PMDI-18 and suggest parameters for the model based on these experiments. In additi on, the second part of this report describes data taken to provide input to the preliminary nonlinear visco elastic structural response model developed for BKC 44306 PMDI-10 foam. We show that the standard cu re schedule used by KCP does not fully cure the material, and, upon temperature elevation above 150°C, oxidation or decomposition reactions occur that alter the composition of the foam. These findings suggest that achieving a fully cured foam part with this formulation may be not be possible through therma l curing. As such, visco elastic characterization procedures developed for curing thermosets can provide only approximate material properties, since the state of the material continuously evolves during tests.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

A model system was developed for enabling a multiscale understanding of centrifugal-contactor liquid–liquid extraction.The system consisted of Nd(III) + xylenol orange in the aqueous phase buffered to pH =5.5 by KHP, and dodecane + thenoyltrifluroroacetone (HTTA) + tributyphosphate (TBP) in the organic phase. Diffusion constants were measured for neodymium in both the organic and aqueous phases, and the Nd(III) partition coefficients were measured at various HTTA and TBP concentrations. A microfluidic channel was used as a high-shear model environment to observe mass-transfer on a droplet scale with xylenol orange as the aqueous-phase metal indicator; mass-transfer rates were measured quantitatively in both diffusion and reaction limited regimes on the droplet scale. Lastly, the microfluidic results were comparable to observations made for the same system in a laboratory scale liquid–liquid centrifugal contactor, indicating that single drop microfluidic experiments can provide information on mass transfer in complicated flows and geometries.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Liquid droplets flowing through a rectangular microfluidic channel develop a vortical flow field due to the presence of shear forces from the surrounding fluid. In this paper, we present an experimental and computational study of droplet velocities and internal flow patterns in a rectangular pressure-driven flow for droplet diameters ranging from 0.1 to 2 times the channel height. Our study shows excellent agreement with asymptotic predictions of droplet and interfacial velocities for infinitesimally small droplets. As the droplet diameter nears the size of the channel height, the droplet velocity slows significantly, and the changing external flow field causes a qualitative change in the location of internal vortices. This behavior is relevant for future studies of mass transfer in microfluidic devices. © 2014 AIP Publishing LLC.

Abstract not provided.

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.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Goma 6.0 is a finite element program which excels in analyses of multiphysical processes, particularly those involving the major branches of mechanics (viz. fluid/solid mechanics, energy transport and chemical species transport). Goma is based on a full-Newton-coupled algorithm which allows for simultaneous solution of the governing principles, making the code ideally suited for problems involving closely coupled bulk mechanics and interfacial phenomena. Example applications include, but are not limited to, coating and polymer processing flows, super-alloy processing, welding/soldering, electrochemical processes, and solid-network or solution film drying. This document serves as a user's guide and reference.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

A thin flow-focusing microfluidic channel is evaluated for generating monodisperse liquid droplets. The microfluidic device is used in its native state, which is hydrophilic, or treated with OTS to make it hydrophobic. Having both hydrophilic and hydrophobic surfaces allows for creation of both oil-in-water and water-in-oil emulsions, facilitating a large parameter study of viscosity ratios (droplet fluid/continuous fluid) ranging from 0.05 to 96 and flow rate ratios (droplet fluid/continuous fluid) ranging from 0.01 to 2 in one geometry. The hydrophilic chip provides a partially-wetting surface (contact angle less than 90°) for the inner fluid. This surface, combined with the unusually thin channel height, promotes a flow regime where the inner fluid wets the top and bottom of the channel in the orifice and a stable jet is formed. Through confocal microscopy, this fluid stabilization is shown to be highly influenced by the contact angle of the liquids in the channel. Non-wetting jets undergo breakup and produce drops when the jet is comparable to or smaller than the channel thickness. In contrast, partially-wetting jets undergo breakup only when they are much smaller than the channel thickness. Drop sizes are found to scale with a modified capillary number based on the total flow rate regardless of wetting behavior. © The Royal Society of Chemistry.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Continuous jets of non-Newtonian fluids impinging on a fluid surface exhibit instabilities from jet buckling and coiling at low Reynolds numbers to delayed die swell, mounding, and air entrainment at higher Reynolds numbers. Filling containers with complex fluids is an important process for many industries, where the need for high throughput requires operating at high Reynolds numbers. In this regime, air entrainment can produce a visually unappealing product, causing a major quality control issue. Just prior to the onset of air entrainment, however, there exists an ideal filling regime which we term " planar filling," as it is characterized by a relatively flat free surface that maintains its shape over time. In this paper, we create a steady-state, 2-D axisymmetric finite element model to study the transition from planar filling to the onset of air entrainment in a container filling process with generalized-Newtonian fluids. We use this model to explore the operating window for Newtonian and shear-thinning (or, more generally, deformation-rate-thinning) fluids, demonstrating that the flow behavior is characterized by a balance between inertial, viscous, and gravitational forces, as characterized by the Reynolds and Froude numbers. A scaling analysis suggests that the relevant parameters for calculating these dimensionless numbers are located where the jet impacts the liquid surface, and simulations show that the transition from planar filling to air entrainment often occurs when Re~O(10). We found that the bottom and side surfaces of the container drastically influence this transition to entrainment, stabilizing the flow. © 2011 Elsevier B.V.

Abstract not provided.

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.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Imported oil exacerabates our trade deficit and funds anti-American regimes. Nuclear Energy (NE) is a demonstrated technology with high efficiency. NE's two biggest political detriments are possible accidents and nuclear waste disposal. For NE policy, proliferation is the biggest obstacle. Nuclear waste can be reduced through reprocessing, where fuel rods are separated into various streams, some of which can be reused in reactors. Current process developed in the 1950s is dirty and expensive, U/Pu separation is the most critical. Fuel rods are sheared and dissolved in acid to extract fissile material in a centrifugal contactor. Plants have many contacts in series with other separations. We have taken a science and simulation-based approach to develop a modern reprocessing plant. Models of reprocessing plants are needed to support nuclear materials accountancy, nonproliferation, plant design, and plant scale-up.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

When a fluid jet impinges on a solid substrate, a variety of behaviors may occur around the impact region. One example is mounding, where the fluid enters the impact region faster than it can flow away, forming a mound of fluid above the main surface. For some operating conditions, this mound can destabilize and buckle, entraining air in the mound. Other behaviors include submerging flow, where the jet impinges into an otherwise steady pool of liquid, entraining a thin air layer as it enters the pool. This impact region is one of very high shear rates and as such, complex fluids behave very differently than do Newtonian fluids. In this work, we attempt to characterize this range of behavior for Newtonian and non-Newtonian fluids using dimensionless parameters. We model the fluid as a modified Bingham-Carreau-Yasuda fluid, which exhibits the full range of pseudoplastic flow properties throughout the impact region. Additionally, we study viscoelastic effects through the use of the Giesekus model. Both 2-D and 3-D numerical simulations are performed using a variety of finite element method techniques for tracking the jet interface, including Arbitrary Lagrangian Eulerian (ALE), diffuse level sets, and a conformal decomposition finite element method (CDFEM). The presence of shear-thinning characteristics drastically reduces unstable mounding behavior, yet can lead to air entrainment through the submerging flow regime. We construct an operating map to understand for what flow parameters mounding and submerging flows will occur, and how the fluid rheology affects these behaviors. This study has many implications in high-speed industrial bottle filling applications.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

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.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

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.

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.

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.

Abstract not provided.

Abstract not provided.

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.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

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.

Abstract not provided.