Many materials of interest to Sandia transition from fluid to solid or have regions of both phases coexisting simultaneously. Currently there are, unfortunately, no material models that can accurately predict this material response. This is relevant to applications that "birth stress" related to geoscience, nuclear safety, manufacturing, energy production and bioscience. Accurately capturing solidification and residual stress enables fully predictive simulations of the evolving front shape or final product. Accurately resolving flow of proppants or blood could reduce environmental impact or lead to better treatments for heart attacks, thrombosis, or aneurism. We will address a science question in this proposal: When does residual stress develop during the critical transition from liquid to solid and how does it affect material deformation? Our hypothesis is that these early phases of stress development are critical to predictive simulation of material performance, net shape, and aging. In this project, we use advanced constitutive models with yield stress to represent both fluid and solid behavior simultaneously. The report provides an abbreviated description of the results from our LDRD "Stress Birth and Death: Disruptive Computational Mechanics and Novel Diagnostics for Fluid-to-Solid Transitions," since we have written four papers that document the work in detail and which we reference. We give highlights of the work and describe the gravitationally driven flow visualization experiment on a model yield stress fluid, Carbopol, at various concentrations and flow rates. We were able to collapse the data on a single master curve by showing it was self-similar. We also describe the Carbopol rheology and the constitutive equations of interest including the Bingham-Carreau-Yasuda model, the Saramito model, and the HB-Saramito model including parameter estimation for the shear and oscillatory rheology. We present several computational models including the 3D moving mesh simulations of both the Saramito models and Bingham-Carreau-Yasuda (BCY) model. We also show results from the BCY model using a 3D level set method and two different ways of handling reduced order Hele-Shaw modeling for generalized Newtonian fluids. We present some first ever two-dimensional results for the modified Jeffries Kamani-Donley-Rogers constitutive equation developed during this project. We include some recent results with a successful Saramito-level set coupling that allows us to tackle problems with complex geometries like mold filling in a thin gap with an obstacle, without the need for remeshing or remapping. We report on some experiments for curing systems where fluorescent particles are used to track material flow. These experiments were carried out in an oven on Sylgard 184 as a model polymerizing system. We conclude the report with a summary of accomplishments and some thoughts on follow-on work.
In response to personal protective equipment (PPE) shortages in the United States due to the Coronavirus Disease 2019, two models of N95 respirators were evaluated for reuse after gamma radiation sterilization. Gamma sterilization is attractive for PPE reuse because it can sterilize large quantities of material through hermetically sealed packaging, providing safety and logistic benefits. The Gamma Irradiation Facility at Sandia National Laboratories was used to irradiate N95 filtering facepiece respirators to a sterilization dose of 25 kGy(tissue). Aerosol particle filtration performance testing and electrostatic field measurements were used to determine the efficacy of the respirators after irradiation. Both respirator models exhibited statistically significant decreases in particle filtering efficiencies and electrostatic potential after irradiation. The largest decrease in capture efficiency was 40–50% and peaked near the 200 nm particle size. The key contribution of this effort is correlating the electrostatic potential change of individual filtration layer of the respirator with the decrease filtration efficiency after irradiation. This observation occurred in both variations of N95 respirator that we tested. Electrostatic potential measurement of the filtration layer is a key indicator for predicting filtration efficiency loss.
Shortages of N95 respirators for use by medical personnel have driven consideration of novel conservation strategies, including decontamination for reuse and extended use. Decontamination methods listed as promising by the Centers for Disease Control and Prevention (CDC) (vaporous hydrogen peroxide (VHP), wet heat, ultraviolet irradiation (UVI)) and several methods considered for low resource environments (bleach, isopropyl alcohol and detergent/soap) were studied for two commonly used surgical N95 respirators (3M™ 1860 and 1870+ Aura™). Although N95 filtration performance depends on the electrostatically charged electret filtration layer, the impact of decontamination on this layer is largely unexplored. As such, respirator performance following decontamination was assessed based on the fit, filtration efficiency, and pressure drop, along with the relationship between (1) surface charge of the electret layer, and (2) elastic properties of the straps. Decontamination with VHP, wet heat, UVI, and bleach did not degrade fit and filtration performance or electret charge. Isopropyl alcohol and soap significantly degraded fit, filtration performance, and electret charge. Pressure drop across the respirators was unchanged. Modest degradation of N95 strap elasticity was observed in mechanical fatigue testing, a model for repeated donnings and doffings. CDC recommended decontamination methods including VHP, wet heat, and UV light did not degrade N95 respirator fit or filtration performance in these tests. Extended use of N95 respirators may degrade strap elasticity, but a loss of face seal integrity should be apparent during user seal checks. NIOSH recommends performing user seal checks after every donning to detect loss of appropriate fit. Decontamination methods which degrade electret charge such as alcohols or detergents should not be used on N95 respirators. The loss of N95 performance due to electret degradation would not be apparent to a respirator user or evident during a negative pressure user seal check.
N95 respirators became scarce to the general public in mid-to-late March of 2020 due to the SARS-CoV-2 epidemic. By mid-April of 2020, most states in the United States were requiring face coverings to be worn while in public enclosed places and in busy outdoor areas where groups of people were in close proximity. Many resorted to cloth masks, homemade masks, procedure masks obtained through online purchases, and other ad-hoc means. Thus, there was and still is a need to determine the aerosol filtration efficacy of commonly available materials that can be used for homemade mask construction. This study focused on non- woven polymeric fabrics that are readily available for homemade mask construction. The conclusion of this study is that non-woven materials that carry a high electric charge or those that can easily acquire charge had the highest aerosol filtration efficiency per unit of pressure drop. Future work should examine a wider variety of these materials and determine the maximum pressure drop that a nominal homemade mask can withstand before a significant portion of airflow is diverted around the mask. More broadly, a better understanding of the charge state on non-woven materials and impact of that charge state on filtration efficiency is needed.
Kinetic simulations of plasma phenomena during and after formation of the conductive plasma channel of a nanosecond pulse discharge are analyzed and compared to existing experimental measurements. Particle-in-cell with direct simulation Monte Carlo collisions (PIC-DSMC) modeling is used to analyze a discharge in helium at 200 Torr and 300 K over a 1 cm gap. The analysis focuses on physics that would not be reproduced by fluid models commonly used at this high number density and collisionality, specifically non-local and stochastic phenomena. Similar analysis could be used to improve the predictive capability of lower fidelity or reduced order models. First, the modeling results compare favorably with experimental measurements of electron number density, temperature, and 1D electron energy distribution function at the same conditions. Second, it is shown that the ionization wave propagates in a stochastic, stepwise manner, dependent on rare, random ionization events ahead of the ionization wave when the ionization fraction in front of the ionization wave is very low, analagous to the stochastic branching of streamers in 3D. Third, analysis shows high-energy runaway electrons accelerated in the cathode layer produce electron densities in the negative glow region over an order of magnitude above those in the positive column. Future work to develop reduced order models of these two phenomena would improve the accuracy of fluid plasma models without the cost of PIC-DSMC simulations.
Glass-ceramic seals may be the future of hermetic connectors at Sandia National Laboratories. They have been shown capable of surviving higher temperatures and pressures than amorphous glass seals. More advanced finite-element material models are required to enable model-based design and provide evidence that the hermetic connectors can meet design requirements. Glass-ceramics are composite materials with both crystalline and amorphous phases. The latter gives rise to (non-linearly) viscoelastic behavior. Given their complex microstructures, glass-ceramics may be thermorheologically complex, a behavior outside the scope of currently implemented constitutive models at Sandia. However, it was desired to assess if the Simplified Potential Energy Clock (SPEC) model is capable of capturing the material response. Available data for SL 16.8 glass-ceramic was used to calibrate the SPEC model. Model accuracy was assessed by comparing model predictions with shear moduli temperature dependence and high temperature 3-point bend creep data. It is shown that the model can predict the temperature dependence of the shear moduli and 3- point bend creep data. Analysis of the results is presented. Suggestions for future experiments and model development are presented. Though further calibration is likely necessary, SPEC has been shown capable of modeling glass-ceramic behavior in the glass transition region but requires further analysis below the transition region.
We present selected results from a series of Open Stack thermal battery tests performed in FY14 and FY15 and discuss our findings. These tests were meant to provide validation data for the comprehensive thermal battery simulation tools currently under development in Sierra/Aria under known conditions compared with as-manufactured batteries. We are able to satisfy this original objective in the present study for some test conditions. Measurements from each test include: nominal stack pressure (axial stress) vs. time in the cold state and during battery ignition, battery voltage vs. time against a prescribed current draw with periodic pulses, and images transverse to the battery axis from which cell displacements are computed. Six battery configurations were evaluated: 3, 5, and 10 cell stacks sandwiched between 4 layers of the materials used for axial thermal insulation, either Fiberfrax Board or MinK. In addition to the results from 3, 5, and 10 cell stacks with either in-line Fiberfrax Board or MinK insulation, a series of cell-free “control” tests were performed that show the inherent settling and stress relaxation based on the interaction between the insulation and heat pellets alone.
Battery performance, while observed at the macroscale, is primarily governed by the bicontinuous mesoscale network of the active particles and a polymeric conductive binder in its electrodes. Manufacturing processes affect this mesostructure, and therefore battery performance, in ways that are not always clear outside of empirical relationships. Directly studying the role of the mesostructure is difficult due to the small particle sizes (a few microns) and large mesoscale structures. Mesoscale simulation, however, is an emerging technique that allows the investigation into how particle-scale phenomena affect electrode behavior. In this manuscript, we discuss our computational approach for modeling electrochemical, mechanical, and thermal phenomena of lithium-ion batteries at the mesoscale. We review our recent and ongoing simulation investigations and discuss a path forward for additional simulation insights.
The transient transport of electrolytes in thermally-activated batteries is studied using electron probe micro-analysis (EPMA), demonstrating the robust capability of EPMA as a useful tool for studying and quantifying mass transport within porous materials, particularly in difficult environments where classical flow measurements are challenging. By tracking the mobility of bromine and potassium ions from the electrolyte stored within the separator into the lithium silicon anode and iron disulfide cathode, we are able to quantify the transport mechanisms and physical properties of the electrodes including permeability and tortuosity. Due to the micron to submicron scale porous structure of the initially dry anode, a fast capillary pressure driven flow is observed into the anode from which we are able to set a lower bound on the permeability of 10-1 mDarcy. The transport into the cathode is diffusion-limited because the cathode originally contained some electrolyte before activation. Using a transient one-dimensional diffusion model, we estimate the tortuosity of the cathode electrode to be 2.8 ± 0.8.
As LiCoO2 cathodes are charged, delithiation of the LiCoO2 active material leads to an increase in the lattice spacing, causing swelling of the particles. When these particles are packed into a bicontinuous, percolated network, as is the case in a battery electrode, this swelling leads to the generation of significant mechanical stress. In this study we performed coupled electrochemical-mechanical simulations of the charging of a LiCoO2 cathode in order to elucidate the mechanisms of stress generation and the effect of charge rate and microstructure on these stresses. Energy dispersive spectroscopy combined with scanning electron microscopy imaging was used to create 3D reconstructions of a LiCoO2 cathode, and the Conformal Decomposition Finite Element Method is used to automatically generate computational meshes on this reconstructed microstructure. Replacement of the ideal solution Fickian diffusion model, typically used in battery simulations, with a more general non-ideal solution model shows substantially smaller gradients of lithium within particles than is typically observed in the literature. Using this more general model, lithium gradients only appear at states of charge where the open-circuit voltage is relatively constant. While lithium gradients do affect the mechanical stress state in the particles, the maximum stresses are always found in the fully-charged state and are strongly affected by the local details of the microstructure and particle-to-particle contacts. These coupled electrochemical-mechanical simulations begin to yield insight into the partitioning of volume change between reducing pore space and macroscopically swelling the electrode. Finally, preliminary studies that include the presence of the polymeric binder suggest that it can greatly impact stress generation and that it is an important area for future research.
The contribution from loss of Li+ inventory to capacity fade is described for slow rates (C/10) and long-term cycling (up to 80 cycles). It was found through electrochemical testing and ex-situ Raman analysis that at these slow rates, the entirety of capacity loss up to 80 cycles can be explained by loss of Li+ inventory in the cell. The Raman spectrum of LiCoO2 is sensitive to the state of lithiation and can therefore be leveraged to quantify the state of lithiation for individual particles. With these Raman derived estimates, the lithiation state of the cathode in the discharged state is compared to electrochemical data as a function of cycle number. High correlation is found between Raman quantifications of cycleable lithium and the capacity fade. Additionally, the linear relationship between discharge capacity and cell overpotential suggests that the loss of capacity stems from an impedance rise of the electrodes, which based on Li inventory losses, is caused by SEI formation and repair.
The polymer-composite binder used in lithium-ion battery electrodes must both hold the electrodes together and augment their electrical conductivity while subjected to mechanical stresses caused by active material volume changes due to lithiation and delithiation. We have discovered that cyclic mechanical stresses cause significant degradation in the binder electrical conductivity. After just 160 mechanical cycles, the conductivity of polyvinylidene fluoride (PVDF):carbon black binder dropped between 45-75%. This degradation in binder conductivity has been shown to be quite general, occurring over a range of carbon black concentrations, with and without absorbed electrolyte solvent and for different polymer manufacturers. Mechanical cycling of lithium cobalt oxide (LiCoO2 ) cathodes caused a similar degradation, reducing the effective electrical conductivity by 30-40%. Mesoscale simulations on a reconstructed experimental cathode geometry predicted the binder conductivity degradation will have a proportional impact on cathode electrical conductivity, in qualitative agreement with the experimental measurements. Finally, ohmic resistance measurements were made on complete batteries. Direct comparisons between electrochemical cycling and mechanical cycling show consistent trends in the conductivity decline. This evidence supports a new mechanism for performance decline of rechargeable lithium-ion batteries during operation - electrochemically-induced mechanical stresses that degrade binder conductivity, increasing the internal resistance of the battery with cycling.
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.
We present a phenomenological constitutive model that describes the macroscopic behavior of pressed-pellet materials used in molten salt batteries. Such materials include separators, cathodes, and anodes. The purpose of this model is to describe the inelastic deformation associated with the melting of a key constituent, the electrolyte. At room temperature, all constituents of these materials are solid and do not transport cations so that the battery is inert. As the battery is heated, the electrolyte, a constituent typically present in the separator and cathode, melts and conducts charge by flowing through the solid skeletons of the anode, cathode, and separator. The electrochemical circuit is closed in this hot state of the battery. The focus of this report is on the thermal-mechanical behavior of the separator, which typically exhibits the most deformation of the three pellets during the process of activating a molten salt battery. Separator materials are composed of a compressed mixture of a powdered electrolyte, an inert binder phase, and void space. When the electrolyte melts, macroscopically one observes both a change in volume and shape of the separator that depends on the applied boundary conditions during the melt transition. Although porous flow plays a critical role in the battery mechanics and electrochemistry, the focus of this report is on separator behavior under flow-free conditions in which the total mass of electrolyte is static within the pellet. Specific poromechanics effects such as capillary pressure, pressure-saturation, and electrolyte transport between layers are not considered. Instead, a phenomenological model is presented to describe all such behaviors including the melting transition of the electrolyte, loss of void space, and isochoric plasticity associated with the binder phase rearrangement. The model is appropriate for use finite element analysis under finite deformation and finite temperature change conditions. The model reasonably describes the stress dependent volume and shape change associated with dead load compression and spring-type boundary conditions; the latter is relevant in molten salt batteries. Future work will transition the model towards describing the solid skeleton of the separator in the traditional poromechanics context.
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.
Lithium-ion battery electrodes rely on a percolated network of solid particles and binder that must maintain a high electronic conductivity in order to function. Coupled mechanical and electrochemical simulations may be able to elucidate the mechanisms for capacity fade. We present a framework for coupled simulations of electrode mechanics that includes swelling, deformation, and stress generation driven by lithium intercalation. These simulations are performed at the mesoscale, which requires 3D reconstruction of the electrode microstructure from experimental imaging or particle size distributions. We present a novel approach for utilizing these complex reconstructions within a finite element code. A mechanical model that involves anisotropic swelling in response to lithium intercalation drives the deformation. Stresses arise from small-scale particle features and lithium concentration gradients. However, we demonstrate, for the first time, that the largest stresses arise from particle-to-particle contacts, making it important to accurately represent the electrode microstructure on the multi-particle scale. Including anisotropy in the swelling mechanics adds considerably more complexity to the stresses and can significantly enhance peak particle stresses. Shear forces arise at contacts due to the misorientation of the lattice structure. These simulations will be used to study mechanical degradation of the electrode structure through charge/discharge cycles.
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.
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.
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.
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.
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.
In this presentation we examine the accuracy and performance of a suite of discrete-element-modeling approaches to predicting equilibrium and dynamic rheological properties of polystyrene suspensions. What distinguishes each approach presented is the methodology of handling the solvent hydrodynamics. Specifically, we compare stochastic rotation dynamics (SRD), fast lubrication dynamics (FLD) and dissipative particle dynamics (DPD). Method-to-method comparisons are made as well as comparisons with experimental data. Quantities examined are equilibrium structure properties (e.g. pair-distribution function), equilibrium dynamic properties (e.g. short- and long-time diffusivities), and dynamic response (e.g. steady shear viscosity). In all approaches we deploy the DLVO potential for colloid-colloid interactions. Comparisons are made over a range of volume fractions and salt concentrations. Our results reveal the utility of such methods for long-time diffusivity prediction can be dubious in certain ranges of volume fraction, and other discoveries regarding the best formulation to use in predicting rheological response.
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.
Lack of robust manufacturing capabilities have limited our ability to make tailored materials with useful optical and thermal properties. For example, traditional methods such as spontaneous self-assembly of spheres cannot generate the complex structures required to produce a full bandgap photonic crystals. The goal of this work was to develop and demonstrate novel methods of directed self-assembly of nanomaterials using optical and electric fields. To achieve this aim, our work employed laser tweezers, a technology that enables non-invasive optical manipulation of particles, from glass microspheres to gold nanoparticles. Laser tweezers were used to create ordered materials with either complex crystal structures or using aspherical building blocks.
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.
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.
Laser tweezers optical trapping provides a unique noninvasive capability to trap and manipulate particles in solution at the focal point of a laser beam passed through a microscope objective. Additionally, combined with image analysis, interaction forces between colloidal particles can be quantitatively measured. By looking at the displacement of particles within the laser trap due to the presence of a neighboring particle or looking at the relative diffusion of two particles held near each other by optical traps, interparticle interaction forces ranging from pico- to femto-Newtons can be measured. Understanding interaction forces is critical for predicting the behavior of particle dispersions including dispersion stability and flow rheology. Using a new analysis method proposed by Sainis, Germain, and Dufresne, we can simultaneously calculate the interparticle velocity and particle diffusivity which allows direct calculation of the interparticle potential for the particles. By applying this versatile tool, we measure difference in interactions between various phospholipid bilayers that have been coated onto silica spheres as a new type of solid supported liposome. We measure bilayer interactions of several cell membrane lipids under various environmental conditions such as pH and ionic strength and compare the results with those obtained for empty liposomes. These results provide insight into the role of bilayer fluctuations in liposome fusion, which is of fundamental interest to liposome based drug delivery schemes.
Optical tweezers has become a powerful and common tool for sensitive determination of electrostatic interactions between colloidal particles. Recently, two techniques, 'blinking' tweezers and direct force measurements, have become increasingly prevalent in investigations of inter-particle potentials. The 'blinking' tweezers method acquires physical statistics of particle trajectories to determine drift velocities, diffusion coefficients, and ultimately colloidal forces as a function of the center-center separation of two particles. Direct force measurements monitor the position of a particle relative to the center of an optical trap as the separation distance between two continuously trapped particles is gradually decreased. As the particles near each other, the displacement from the trap center for each particle increases proportional to the inter-particle force. Although commonly employed in the investigation of interactions of colloidal particles, there exists no direct comparison of these experimental methods in the literature. In this study, an experimental apparatus was developed capable of performing both methods and is used to quantify electrostatic potentials between particles in several particle/solvent systems. Comparisons are drawn between the experiments conducted using the two measurement techniques, theory, and existing literature. Forces are quantified on the femto-Newton scale and results agree well with literature values.
A method to measure interfacial mechanical properties at high temperatures and in a controlled atmosphere has been developed to study anodized aluminum surface coatings at temperatures where the interior aluminum alloy is molten. This is the first time that the coating strength has been studied under these conditions. We have investigated the effects of ambient atmosphere, temperature, and surface finish on coating strength for samples of aluminum alloy 7075. Surprisingly, the effective Young's modulus or strength of the coating when tested in air was twice as high as when samples were tested in an inert nitrogen or argon atmosphere. Additionally, the effective Young's modulus of the anodized coating increased with temperature in an air atmosphere but was independent of temperature in an inert atmosphere. The effect of surface finish was also examined. Sandblasting the surface prior to anodization was found to increase the strength of the anodized coating with the greatest enhancement noted for a nitrogen atmosphere. Machining marks were not found to significantly affect the strength.
Here, we have developed a novel flow chamber which imposes a controlled axisymmetric stagnation flow to enable the study of external flow effects on coalescence dynamics. This system allows for the first time the precise positioning of a drop in a three dimensional flow and additionally enforces a highly symmetric flow around the drop. We focus on the study of a single drop approaching a stationary flat plane as this is analogous to two drops approaching each other. A single drop is created and then guided along the unsteady center line of a stagnation flow. The real time computer control algorithm analyzes video images of the drop in two orthogonal planes and manipulates flow restricting valves along the four outlets of the flow. We demonstrate using particle image velocimetry that the computer control not only controls the drop position but also ensures a symmetric flow inside the flow chamber. Finally, this chamber will enable a detailed investigation of the drainage of the thin film between the drop and the lower surface in order to probe the effect of external flow on coalescence.
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.
We have investigated a novel emulsion interfacial filter that is applicable for a wide range of materials, from nano-particles to cells and bacteria. This technology uses the interface between the two immiscible phases as the active surface area for adsorption of targeted materials. We showed that emulsion interfaces can effectively collect and trap materials from aqueous solution. We tested two aqueous systems, a bovine serum albumin (BSA) solution and coal bed methane produced water (CBMPW). Using a pendant drop technique to monitor the interfacial tension, we demonstrated that materials in both samples were adsorbed to the liquid-liquid interface, and did not readily desorb. A prototype system was built to test the emulsion interfacial filter concept. For the BSA system, a protein assay showed a progressive decrease in the residual BSA concentration as the sample was processed. Based on the initial prototype operation, we propose an improved system design.
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.
The goal of this project was to design, build and test novel diagnostics to probe the effect of hydrodynamic forces on coalescence dynamics. Our investigation focused on how a drop coalesces onto a flat surface which is analogous to two drops coalescing, but more amenable to precise experimental measurements. We designed and built a flow cell to create an axisymmetric compression flow which brings a drop onto a flat surface. A computer-controlled system manipulates the flow to steer the drop and maintain a symmetric flow. Particle image velocimetry was performed to confirm that the control system was delivering a well conditioned flow. To examine the dynamics of the coalescence, we implemented an interferometry capability to measure the drainage of the thin film between the drop and the surface during the coalescence process. A semi-automated analysis routine was developed which converts the dynamic interferogram series into drop shape evolution data.
Nearly every manufacturing and many technologies central to Sandia's business involve physical processes controlled by interfacial wetting. Interfacial forces, e.g. conjoining/disjoining pressure, electrostatics, and capillary condensation, are ubiquitous and can surpass and even dominate bulk inertial or viscous effects on a continuum level. Moreover, the statics and dynamics of three-phase contact lines exhibit a wide range of complex behavior, such as contact angle hysteresis due to surface roughness, surface reaction, or compositional heterogeneities. These thermodynamically and kinetically driven interactions are essential to the development of new materials and processes. A detailed understanding was developed for the factors controlling wettability in multicomponent systems from computational modeling tools, and experimental diagnostics for systems, and processes dominated by interfacial effects. Wettability probed by dynamic advancing and receding contact angle measurements, ellipsometry, and direct determination of the capillary and disjoining forces. Molecular scale experiments determined the relationships between the fundamental interactions between molecular species and with the substrate. Atomistic simulations studied the equilibrium concentration profiles near the solid and vapor interfaces and tested the basic assumptions used in the continuum approaches. These simulations provide guidance in developing constitutive equations, which more accurately take into account the effects of surface induced phase separation and concentration gradients near the three-phase contact line. The development of these accurate models for dynamic multicomponent wetting allows improvement in science based engineering of manufacturing processes previously developed through costly trial and error by varying material formulation and geometry modification.
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.
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.
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.
Filling operations, in which a viscous fluid displaces a gas in a complex geometry, occur with surprising frequency in many manufacturing processes. Difficulties in generating accurate models of these processes involve accurately capturing the interfacial boundary as it undergoes large motions and deformations, preventing dispersion and mass-loss during the computation, and robustly accounting for the effects of surface tension and wetting phenomena. This paper presents a numerical capturing algorithm using level set theory and finite element approximation. Important aspects of this work are addressing issues of mass-conservation and the presence of wetting effects. We have applied our methodology to a three-dimension model of a complicated filling problem. The simulated results are compared to experimental flow visualization data taken for filling of UCON oil in the identical geometry. Comparison of simulation and experiment indicates that the simulation conserved mass adequately and the simulated interface shape was in approximate agreement with experiment. Differences seen were largely attributed to inaccuracies in the wetting line model.