Brazing and soldering are metallurgical joining techniques that use a wetting molten metal to create a joint between two faying surfaces. The quality of the brazing process depends strongly on the wetting properties of the molten filler metal, namely the surface tension and contact angle, and the resulting joint can be susceptible to various defects, such as run-out and underfill, if the material properties or joining conditions are not suitable. In this work, we implement a finite element simulation to predict the formation of such defects in braze processes. This model incorporates both fluid–structure interaction through an arbitrary Eulerian–Lagrangian technique and free surface wetting through conformal decomposition finite element modeling. Upon validating our numerical simulations against experimental run-out studies on a silver-Kovar system, we then use the model to predict run-out and underfill in systems with variable surface tension, contact angles, and applied pressure. Finally, we consider variable joint/surface geometries and show how different geometrical configurations can help to mitigate run-out. This work aims to understand how brazing defects arise and validate a coupled wetting and fluid–structure interaction simulation that can be used for other industrial problems.
A new particle-based reweighting method is developed and demonstrated in the Aleph Particle-in-Cell with Direct Simulation Monte Carlo (PIC-DSMC) program. Novel splitting and merging algorithms ensure that modified particles maintain physically consistent positions and velocities. This method allows a single reweighting simulation to efficiently model plasma evolution over orders of magnitude variation in density, while accurately preserving energy distribution functions (EDFs). Demonstrations on electrostatic sheath and collisional rate dynamics show that reweighting simulations achieve accuracy comparable to fixed weight simulations with substantial computational time savings. This highly performant reweighting method is recommended for modeling plasma applications that require accurate resolution of EDFs or exhibit significant density variations in time or space.
We describe a data-driven, multiscale technique to model reactive wetting of a silver–aluminum alloy on a Kovar™ (Fe-Ni-Co alloy) surface. We employ molecular dynamics simulations to elucidate the dependence of surface tension and wetting angle on the drop's composition and temperature. A design of computational experiments is used to efficiently generate training data of surface tension and wetting angle from a limited number of molecular dynamics simulations. The simulation results are used to parameterize models of the material's wetting properties and compute the uncertainty in the models due to limited data. The data-driven models are incorporated into an engineering-scale (continuum) model of a silver–aluminum sessile drop on a Kovar™ substrate. Model predictions of the wetting angle are compared with experiments of pure silver spreading on Kovar™ to quantify the model-form errors introduced by the limited training data versus the simplifications inherent in the molecular dynamics simulations. The paper presents innovations in the determination of “convergence” of noisy MD simulations before they are used to extract the wetting angle and surface tension, and the construction of their models which approximate physio-chemical processes that are left unresolved by the engineering-scale model. Together, these constitute a multiscale approach that integrates molecular-scale information into continuum scale models.
As additive manufacturing (AM) has become a reliable method for creating complex and unique hardware rapidly, the quality assurance of printed parts remains a priority. In situ process monitoring offers an approach for performing quality control while simultaneously minimizing post-production inspection. For extrusion printing processes, direct linkages between extrusion pressure fluctuations and print defects can be established by integrating pressure sensors onto the print head. In this work, the sensitivity of process monitoring is tested using engineered spherical defects. Pressure and force sensors located near an ink reservoir and just before the nozzle are shown to assist in identification of air bubbles, changes in height between the print head and build surface, clogs, and particle aggregates with a detection threshold of 60–70% of the nozzle diameter. Visual evidence of printed bead distortion is quantified using optical image analysis and correlated to pressure measurements. Importantly, this methodology provides an ability to monitor the quality of AM parts produced by extrusion printing methods and can be accomplished using commonly available pressure-sensing equipment.
Carbon capture is essential to meeting climate change mitigation goals. One approach currently being commercialized utilizes liquid-based solvents to capture CO2 directly from the atmosphere but is limited by slow absorption of CO2 into the liquid. Improved air/solvent liquid mixing increases CO2 absorption rate, and this increased CO2 absorption efficiency allows for smaller carbon capture systems with lower capital costs and better economic viability. In this project, we study the use of passive micromixers fabricated by metal additive manufacturing. The micromixer’s small-scale surface geometric features perturb and mix the liquid film to enhance mass transfer and CO2 absorption. In this project, we evaluated this hypothesis through computational and experimental studies. Computational investigations focused on developing capabilities to simulate thin film (~ 100μm) fluid flow on rough surfaces. Such thin films are in a surface-tension dominated regime and simulations in this regime are prone to instabilities. Improvements to the Nalu code completed in this project resulted in a 10x timestep stability improvement for these problems.
Many important engineering and scientific applications such as cement slurries, foams, crude oil, and granular avalanches involve the concept of yield stress. Therefore, modeling yield stress fluids in different flow configurations, including the accurate prediction of the yield surface, is important. In this paper, we present a computational model based on the finite element method to study the flow of yield stress fluids in a thin mold and compare the results with data from flow visualization experiments. We use the level set method to describe the interface between the filling fluid and air. We use polypropylene glycol as a model Newtonian fluid and Carbopol for the model yield stress fluid, as the Carbopol solution demonstrates yielding without thixotropy. To describe the yielding and shear-thinning behavior, we use a generalized Newtonian constitutive equation with a Bingham–Carreau–Yasuda form. We compare the results obtained from the mold filling experiments with the results from the three-dimensional (3D) model and from a reduced-order Hele-Shaw (HS) model that is two-dimensional, including the effect of shear-thinning along the thin direction only approximately. We show that both the 3D and the HS model can capture the experimental meniscus shape reasonably well for all the fluids considered at three different flow rates. This indicates that the shape evolution is insensitive to the dimensionality of the model. However, the viscosity and yield surfaces predicted by the 3D and HS models are different. The HS model underestimates the high viscosity and unyielded regions compared to the estimation by the 3D model.
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.
We propose developing numerical methods informed by novel experimental diagnostics that transition from solid-to-fluid, while accurately predicting the stress and deformation regardless of phase.
Filtration, pressure drop and quantitative fit of N95 respirators were robust to several decontamination methods including vaporous hydrogen peroxide, wet heat, bleach, and ultraviolet light. Bleach may not have penetrated the hydrophobic outer layers of the N95 respirator. Isopropyl alcohol and detergent both severely degraded the electrostatic charge of the electret filtration layer. First data in N95 respirators that the loss of filtration efficiency was directly correlated with loss of surface potential on the filtration layer. The pressure drop was unchanged, so loss of filtration efficacy would not be apparent during a user seal check. Mechanical straps degrade with repeated mechanical cycling during extended use. Decontamination did not appear to degrade the elastic straps. Significant loss of strap elasticity would be apparent during a user negative pressure seal check.
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.
Sandia Materials Science Investment Area contributed to the SARS-CoV-2 virus and COVID-19 disease which represent the most significant pandemic threat in over 100 years. We completed a series of 7, short duration projects to provide innovative materials science research and development in analytical techniques to aid the neutralization of COVID-19 on multiple surfaces, approaches to rapidly decontaminate personal protective equipment, and pareto assessment of construction materials for manufacturing personal protective equipment. The developed capabilities and processes through this research can help US medical personnel, government installations and assets, first responders, state and local governments, and multiple federal agencies address the COVID-19 Pandemic.
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.
This study evaluated gamma irradiation for sterilization and reuse of two models of N95 respirators after gamma radiation sterilization as a method to increase availability of N95 respirators during a shortage. The Sandia National Laboratories Gamma Irradiation Facility was used to irradiate two different models of N95 filtering facepiece respirators at doses ranging from 0 kGy(tissue) to 50 kGy(tissue). The following tests were used to determine the efficacy of the respirator after irradiation sterilization: Ambient Aerosol Condensation Nuclei Counter Quantitative Fit Test, tensile test, strain cycling, oscillatory dynamic mechanical analysis, microscopic image analysis of fiber layers, and electrostatic field measurements. Both of the respirator models exhibited statistically significant changes after gamma irradiation as shown by the Quantitative Fit Test, electrostatic testing and the aerosol testing. The change in electrostatic capability of the filter reduced the efficiency of challenging particles near the 200 nm size by approximately 40-50%. Both tested respirators showed statistically significant changes associated with gamma sterilization. However, our results indicate that choices in materials and manufacturing methods to achieve N95 filtration lead to different magnitudes of damage when exposed to gamma radiation at sterilization relevant doses. This damage results in lower filtration performance. While our sample size (2 different types of respirators) was small, we did observe a change in electrostatic properties on a filter layer that coincided with the failure on the Quantitative Fit Test.
This study evaluated gamma irradiation for sterilization and reuse of two models of N95 respirators after gamma radiation sterilization as a method to increase availability of N95 respirators during a shortage. The Sandia National Laboratories Gamma Irradiation Facility was used to irradiate two different models of N95 filtering facepiece respirators at doses ranging from 0 kGy(tissue) to 50 kGy(tissue). The following tests were used to determine the efficacy of the respirator after irradiation sterilization: Ambient Aerosol Condensation Nuclei Counter Quantitative Fit Test, tensile test, strain cycling, oscillatory dynamic mechanical analysis, microscopic image analysis of fiber layers, and electrostatic field measurements. Both of the respirator models exhibited statistically significant changes after gamma irradiation as shown by the Quantitative Fit Test, electrostatic testing and the aerosol testing. The change in electrostatic charge of the filter was correlated with a reduction in capturing particles near the 200 nm size by approximately 40-50%. Both tested respirators showed statistically significant changes associated with gamma sterilization. However, our results indicate that choices in materials and manufacturing methods to achieve N95 filtration lead to different magnitudes of damage when exposed to gamma radiation at sterilization relevant doses. This damage results in lower filtration performance. While our sample size (2 different types of respirators) was small, we did observe a change in electrostatic properties on a filter layer that coincided with the failure on the Quantitative Fit Test and reduction in aerosol filtering efficiency. Key Words: N95 respirators, respirators, airborne transmission, pandemic prevention, COVID-19, gamma sterilization
Sandia National Laboratories currently has 27 COVID-related Laboratory Directed Research & Development (LDRD) projects focused on helping the nation during the pandemic. These LDRD projects cross many disciplines including bioscience, computing & information sciences, engineering science, materials science, nanodevices & microsystems, and radiation effects & high energy density science.
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.
Understanding the stress development in fluids as they transition to solids is not well-understood. Computational models are needed to represent "birthing stress" for multiphysics applications such as polymer encapsulation around sensitive electronics and additive manufacturing where these stresses can lead to defects such as cracking and voids. The local stress state is also critical to understand and predict the net shape of parts formed in the liquid phase. In this one-year exploratory LDRD, we have worked towards a novel experimental diagnostic to measure the fluid rheology, degree of solidification, and the solid stress development simultaneously. We debugged and made viable a "first-generation" Rheo-Raman system and used it to characterize two types of solidifying systems: paraffin wax, which crystalizes as it solidifies, and thermoset polymers, which form a network of covalent bonds. We used the paraffin wax as a model system to perform flow visualization studies and did some preliminary modeling of the experiment, to demonstrate the inadequacy of the current modeling approaches. This work will inform an advanced fluid constitutive equation that includes a yield stress, temperature dependence, and an evolving viscosity when we pursue the full proposal, which was funded for FY20.
Lithium batteries provide high energy density storage with applications ranging from consumer electronics to electric vehicles. However, they have a limited lifespan and experience capacity loss with aging. Multiple mechanisms contribute to battery aging. The battery binder plays two important roles in the electrodes, and the damage it sustains during cycling may play a role in the degradation of the overall battery performance. Mechanical stress during battery operations occurs as a result of the swelling and shrinking of the electrodes because of the movement of lithium with cycling. The yield stress of the swollen polyvinylidene fluoride carbon black (PVDFCB) binder was measured at approximately 4MPa for PVDF with carbon black CB weight fractions between 10-30% swollen in propylene carbonate. This is far less stress than is typically experienced in an electrode during cycling. The effects of this permanent damage to the binder were explored by measuring the conductivity loss with strains in excess of the binder yield.
