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.
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.
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.
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.
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.
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.
Solidification and blood flow seemingly have little in common, but each involves a fluid in contact with a deformable solid. In these systems, the solid-fluid interface moves as the solid advects and deforms, often traversing the entire domain of interest. Currently, these problems cannot be simulated without innumerable expensive remeshing steps, mesh manipulations or decoupling the solid and fluid motion. Despite the wealth of progress recently made in mechanics modeling, this glaring inadequacy persists. We propose a new technique that tracks the interface implicitly and circumvents the need for remeshing and remapping the solution onto the new mesh. The solid-fluid boundary is tracked with a level set algorithm that changes the equation type dynamically depending on the phases present. This novel approach to coupled mechanics problems promises to give accurate stresses, displacements and velocities in both phases, simultaneously.