Nanoparticles are now more than ever being used to tailor materials function and performance in differentiating technologies because of their profound effect on thermo-physical, mechanical and optical properties. The most feasible way to disperse particles in a bulk material or control their packing at a substrate is through fluidization in a carrier, followed by solidification through solvent evaporation/drying/curing/sintering. Unfortunately processing particles as concentrated, fluidized suspensions into useful products remains an art largely because the effect of particle shape and volume fraction on fluidic properties and suspension stability remains unexplored in a regime where particle-particle interaction mechanics is prevalent. To achieve a stronger scientific understanding of the factors that control nanoparticle dispersion and rheology we have developed a multiscale modeling approach to bridge scales between atomistic and molecular-level forces active in dense nanoparticle suspensions. At the largest length scale, two 'coarse-grained' numerical techniques have been developed and implemented to provide for high-fidelity numerical simulations of the rheological response and dispersion characteristics typical in a processing flow. The first is a coupled Navier-Stokes/discrete element method in which the background solvent is treated by finite element methods. The second is a particle based method known as stochastic rotational dynamics. These two methods provide a new capability representing a 'bridge' between the molecular scale and the engineering scale, allowing the study of fluid-nanoparticle systems over a wide range of length and timescales as well as particle concentrations. To validate these new methodologies, multi-million atoms simulations explicitly including the solvent have been carried out. These simulations have been vital in establishing the necessary 'subgrid' models for accurate prediction at a larger scale and refining the two coarse-grained methodologies.
This paper presents continuum simulations of viscous polymer flow during nanoimprint lithography (NIL) for embossing tools having irregular spacings and sizes. Simulations varied non-uniform embossing tool geometry to distinguish geometric quantities governing cavity filling order, polymer peak deformation, and global mold filling times. A characteristic NIL velocity predicts cavity filling order. In general, small cavities fill more quickly than large cavities, while cavity spacing modulates polymer deformation mode. Individual cavity size, not total filling volume, dominates replication time, with large differences in individual cavity size resulting in non-uniform, squeeze flow filling. High density features can be modeled as a solid indenter in squeeze flow to accurately predict polymer flow and allow for optimization of wafer-scale replication. The present simulations make it possible to design imprint templates capable of distributing pressure evenly across the mold surface and facilitating symmetric polymer flow over large areas to prevent mold deformation and non-uniform residual layer thickness.
This paper presents continuum simulations of polymer flow during nanoimprint lithography (NIL). The simulations capture the underlying physics of polymer flow from the nanometer to millimeter length scale and examine geometry and thermophysical process quantities affecting cavity filling. Variations in embossing tool geometry and polymer film thickness during viscous flow distinguish different flow driving mechanisms. Three parameters can predict polymer deformation mode: cavity width to polymer thickness ratio, polymer supply ratio, and Capillary number. The ratio of cavity width to initial polymer film thickness determines vertically or laterally dominant deformation. The ratio of indenter width to residual film thickness measures polymer supply beneath the indenter which determines Stokes or squeeze flow. The local geometry ratios can predict a fill time based on laminar flow between plates, Stokes flow, or squeeze flow. Characteristic NIL capillary number based on geometry-dependent fill time distinguishes between capillary or viscous driven flows. The three parameters predict filling modes observed in published studies of NIL deformation over nanometer to millimeter length scales. The work seeks to establish process design rules for NIL and to provide tools for the rational design of NIL master templates, resist polymers, and process parameters.
This article presents continuum simulations of viscous polymer flow during nanoimprint lithography (NIL) for embossing tools having irregular spacings and sizes. Simulations vary nonuniform embossing tool geometry to distinguish geometric quantities governing cavity filling order, polymer peak deformation, and global mold filling times. A characteristic NIL velocity predicts cavity filling order. In general, small cavities fill more quickly than large cavities, while cavity spacing modulates polymer deformation mode. Individual cavity size, not total filling volume, dominates replication time, with large differences in individual cavity size resulting in nonuniform, squeeze flow filling. High density features can be modeled as a solid indenter in squeeze flow to accurately predict polymer flow and allow for optimization of wafer-scale replication. The present simulations make it possible to design imprint templates capable of distributing pressure evenly across the mold surface and facilitating symmetric polymer flow over large areas to prevent mold deformation and nonuniform residual layer thickness.
This paper presents continuum simulations of polymer flow during nanoimprint lithography (NIL). The simulations capture the underlying physics of polymer flow from the nanometer to millimeter length scale and examine geometry and thermophysical process quantities affecting cavity filling. Variations in embossing tool geometry and polymer film thickness during viscous flow distinguish different flow driving mechanisms. Three parameters can predict polymer deformation mode: cavity width to polymer thickness ratio, polymer supply ratio and capillary number. The ratio of cavity width to initial polymer film thickness determines vertically or laterally dominant deformation. The ratio of indenter width to residual film thickness measures polymer supply beneath the indenter which determines Stokes or squeeze flow. The local geometry ratios can predict a fill time based on laminar flow between plates, Stokes flow, or squeeze flow. A characteristic NIL capillary number based on geometry-dependent fill time distinguishes between capillary- or viscous-driven flows. The three parameters predict filling modes observed in published studies of NIL deformation over nanometer to millimeter length scales. The work seeks to establish process design rules for NIL and to provide tools for the rational design of NIL master templates, resist polymers and process parameters.
This report summarizes research advances pursued with award funding issued by the DOE to Drexel University through the Presidential Early Career Award (PECASE) program. Professor Rich Cairncross was the recipient of this award in 1997. With it he pursued two related research topics under Sandia's guidance that address the outstanding issue of fluid-structural interactions of liquids with deformable solid materials, focusing mainly on the ubiquitous dynamic wetting problem. The project focus in the first four years was aimed at deriving a predictive numerical modeling approach for the motion of the dynamic contact line on a deformable substrate. A formulation of physical model equations was derived in the context of the Galerkin finite element method in an arbitrary Lagrangian/Eulerian (ALE) frame of reference. The formulation was successfully integrated in Sandia's Goma finite element code and tested on several technologically important thin-film coating problems. The model equations, the finite-element implementation, and results from several applications are given in this report. In the last year of the five-year project the same physical concepts were extended towards the problem of capillary imbibition in deformable porous media. A synopsis of this preliminary modeling and experimental effort is also discussed.
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.
Encapsulation is a common process used in manufacturing most non-nuclear components including: firing sets, neutron generators, trajectory sensing signal generators (TSSGs), arming, fusing and firing devices (AF and Fs), radars, programmers, connectors, and batteries. Encapsulation is used to contain high voltage, to mitigate stress and vibration and to protect against moisture. The purpose of the ASCI Encapsulation project is to develop a simulation capability that will allow us to aid in the encapsulation design process, especially for neutron generators. The introduction of an encapsulant poses many problems because of the need to balance ease of processing and properties necessary to achieve the design benefits such as tailored encapsulant properties, optimized cure schedule and reduced failure rates. Encapsulants can fail through fracture or delamination as a result of cure shrinkage, thermally induced residual stresses, voids or incomplete component embedding and particle gradients. Manufacturing design requirements include (1) maintaining uniform composition of particles in order to maintain the desired thermal coefficient of expansion (CTE) and density, (2) mitigating void formation during mold fill, (3) mitigating cure and thermally induced stresses during cure and cool down, and (4) eliminating delamination and fracture due to cure shrinkage/thermal strains. The first two require modeling of the fluid phase, and it is proposed to use the finite element code GOMA to accomplish this. The latter two require modeling of the solid state; however, ideally the effects of particle distribution would be included in the calculations, and thus initial conditions would be set from GOMA predictions. These models, once they are verified and validated, will be transitioned into the SIERRA framework and the ARIA code. This will facilitate exchange of data with the solid mechanics calculations in SIERRA/ADAGIO.
Living systems exhibit form and function on multiple length scales and at multiple locations. In order to mimic such natural structures, it is necessary to develop efficient strategies for assembling hierarchical materials. Conventional photolithography, although ubiquitous in the fabrication of microelectronics and microelectromechanical systems, is impractical for defining feature sizes below 0.1 micrometres and poorly suited to pattern chemical functionality. Recently, so-called 'soft' lithographic approaches have been combined with surfactant and particulate templating procedures to create materials with multiple levels of structural order. But the materials thus formed have been limited primarily to oxides with no specific functionality, and the associated processing times have ranged from hours to days. Here, using a self-assembling 'ink', we combine silica-surfactant self-assembly with three rapid printing procedures-pen lithography, ink-jet printing, and dip-coating of patterned self-assembled monolayers-to form functional, hierarchically organized structures in seconds. The rapid-prototyping procedures we describe are simple, employ readily available equipment, and provide a link between computer-aided design and self-assembled nanostructures. We expect that the ability to form arbitrary functional designs on arbitrary surfaces will be of practical importance for directly writing sensor arrays and fluidic or photonic systems.
Electron-Beam (EBeam) melting furnaces are routine used to minimize the occurrence of second-phase particles in the processing of segregation-sensitive alloys. As one part of the process, a circulating electron beam impinges the surface of a crucible melt pool to help control the shape of the solidification front below. By modeling melt pool hydrodynamics, heat transfer, and the shape of solidification boundaries, we plan to optimize the dwell pattern of the beam so that the material solidifies with a composition as spatially homogeneous as possible. Both two- and three-dimensional models are being pursued with FIDAP 5.02, the former serving as a test bed for various degrees of model sophistication. A heat flux distribution is specified on the top of the domain to simulate the EBeam dwell pattern. In two dimensions it is found that an inertially-driven recirculation in the melt pool interacts with a counter-rotating buoyancy-driven recirculation, and that both recirculations are influenced heavily by surface tension gradients on the melt-pool surface. In three dimensions the inertial cell decays quickly with distance from the position of the inlet stream, causing the fluid to precess the crucible. Ingot macrosegregation patterns for a U-6 wt. % Nb alloy are calculated with the Flemings-Mehrabian equation of solute redistribution; the sensitivity of these patterns to EBeam dwell pattern is explored.
Dip coating is the primary means of depositing sol-gel films for precision optical coatings. Sols are typically multicomponent systems consisting of an inorganic phase dispersed in a solvent mixture, with each component differing in volatility and surface tension. This, together with slow coating speeds (<1cm/s), makes analysis of the coating process complicated; unlike most high-speed coating methods, solvent evaporation, evolving rheology, and surface tension gradients alter significantly the fluid mechanics of the deposition stage. We set out to understand these phenomena with computer-aided predictions of the flow and species transport fields. The underlying theory involves mass, momentum, and species transport on a domain of unknown shape, with models and constitutive equations for vapor-liquid equilibria and surface tension. Due accounting is made for the unknown position of the free surface, which locates according to the capillary hydrodynamic forces and solvent loss by evaporation. Predictions of the effects of mass transfer, hydrodynamics, and surface tension gradients on final film thickness are compared with ellipsometry measurements of film thickness on a laboratory pilot coater. Although quantitative agreement is still lacking, both experiment and theory reveal that the film profile near the drying line takes on a parabolic shape. 2 refs., 2 figs.