On the Dewetting of Liquefied Metal Nanostructures
|
|
- Donna Webster
- 5 years ago
- Views:
Transcription
1 On the Dewetting of Liquefied Metal Nanostructures Shahriar Afkhami and Lou Kondic Department of Mathematical Sciences, New Jersey Institute of Technology, Newark, NJ USA June 10, 2013 Abstract. Direct numerical simulations of liquefied metal nanostructures dewetting a substrate are carried out. Full three-dimensional Navier-Stokes equations are solved and a volume-of-fluid method is used for tracking and locating the interface. Substrate wettability is varied to study the influence of the solid/liquid interaction. The effects of initial geometry on the retraction dynamics is numerically investigated. It is shown that the dewetting velocity increases with increasing the contact angle, and that the retraction dynamics is governed by an elaborate interplay of initial geometry, inertial and capillary forces, and the dewetting phenomena. Numerical results are presented of the dewetting of nanoscale Cu and Au liquefied structures on a substrate. Keywords: Liquefied Metal Nanostructures, Dewetting, Navier-Stokes, Surface Tension, Volume-of-Fluid 1. Introduction Modeling of fluid problems characterized by length scale measured in 10s to 100s of nanometers involves considering a number of issues, in particular when free boundaries are involved. This modeling is further complicated in the setups involving the presence of fluid/solid interfaces. One member of a family of problems where these effects are relevant involves evolution of liquefied metals. Very recently, these materials have been extensively considered for the synthesis and assembly on nanoscale (see e.g. (Habenicht et al., 2005; Habenicht et al., 2008; Roberts et al., 2013)). The dewetting of liquefied gold (Au) nanoparticles are experimentally considered by Habenicht et al. (2005; 2008), for example. In these physical experiments, the long time dynamics of the nanodroplet is governed by the competition of the inertial, capillary, and surface forces, giving rise to droplet contraction toward a global energy minimizer; if the inertial forces overcome the adhesion energy and dissipative losses, the nanodroplet is ejected from the substrate on a time scale of the order of 10 ns. The small length and short time scales involved in this process make the experimental study of the dewetting dynamics a daunting task. Discrete molecular dynamics (MD) approaches can be utilized to simulate the systems on nanoscale (Fuentes-Cabrera et al., 2011a; Fuentesc 2013 Kluwer Academic Publishers. Printed in the Netherlands. JEnggMaths_KA-v2.tex; 10/06/2013; 18:32; p.1
2 2 Cabrera et al., 2011b). Fuentes-Cabrera et al., for example, carried out MD simulation using a Lennard-Jones (LJ) interface potential to study various wetting regimes of nanoscale liquefied copper (Cu) nanodisks on graphite (Fuentes-Cabrera et al., 2011a). They showed that the dewetting velocities are consistent with inertial flow estimates, and that thicker and smaller radius nanodisks have proportionately lower velocities. They have also found that for a particular choice of the LJ potential, the dewetting nanoparticles jump off from the graphitic substrate with velocities that are similar to those observed experimentally in (Habenicht et al., 2005) for Au liquefied nanostructures on glass. While the MD approach is very powerful, it is computationally expensive. An alternative approach is to use continuum modeling although it is unclear that the assumptions of continuum fluid mechanics are satisfied on these very short spatial and temporal scales. In addition to concerns related to validity of continuum models on nanoscale, modeling of fluid structures on substrates involves additional degree of complexity since one needs to carefully include the liquid/solid interaction forces, which may be dominant on such short length scales. Time-dependent evolution necessarily involves consideration of contact lines, where liquid, solid, and gas phase meet and modeling of the physics in the vicinity of contact lines is demanding and challenging, due to well known issues related to the so-called contact line singularity (de Gennes, 1985; Haley and Miksis, 1991). Due to all these complexities, typically the problem is considered within the long-wave theory, although in the context of liquid metal films this approach is questionable due to large contact angles involved. Despite the concerns about its validity, this approach has turned out to be very useful, providing significant insight to the dynamics and stability properties of fluids on nanoscale, including liquefied metals (Trice et al., 2007; Ajaev and Willis, 2003; Kondic et al., 2009). An important property of the numerical schemes that needs to be developed for NS equations on nanoscale is the ability to accurately and robustly deal with the fluid-fluid as well as fluid-solid interactions. Furthermore, capillary forces are often very important on nanoscale, and therefore the efficacy of a numerical scheme is closely dependent on the ability of both the surface tension disctretization and the flow solver in minimizing the introduction and amplification of errors. Accurate computation of the geometrical properties of the interface such as interfacial normals and curvature, sharp surface force representation, and accurate and stable simulation are key ingredients in feasible computation on nanoscale. Here, we use the volume-of-fluid (VoF) interface capturing methodology, which is a widespread approach for representing fluid interfaces and, when coupled to a flow solver, is also used to compute JEnggMaths_KA-v2.tex; 10/06/2013; 18:32; p.2
3 quantities related to surface tension that enter the flow calculation. Until recently, however, the VoF method was not generally deemed appropriate for the study of phenomena in which surface tension is the driving force (e.g. the wetting and dewetting of solid surfaces, as occurs in droplet and spray technologies as well as in various natural phenomena). However, recent improvements to calculating curvature and applying the surface tension force appear to resolve this issue (Afkhami and Bussmann, 2008; Afkhami and Bussmann, 2009; Afkhami et al., 2009). In this paper, we use numerical simulations based on full threedimensional (3D) NS equations based on a VoF method to study the dewetting of liquefied nano-sized structures on substrates with varying wettability. The motion of the fluid both inside the drop and in the surrounding flow are numerically modeled. Unlike most of other theoretical and computational studies, the Reynolds number and the contact angle need not be small. In this work, we explore whether our continuum formulation will result in simulations consistent with MD simulations using LJ potentials. We initially focus our simulations on the dewetting behavior of liquid Cu nanodisks and ring-shaped nanostructures. We will then consider Au triangular nanostructures, that do not possess radial symmetry Mathematical Model The equations of conservation of mass, u = 0, and momentum, D(ρu)/Dt = ( pi+ τ)+f st, (NS equations), govern the fluid motion in the drop and the surrounding gas. Here u is the velocity vector, p is the pressure, ρ is the fluid density, τ is the shear stress tensor, and F st is the surface tension force (per unit volume). The shear stress tensor is defined as τ = µ( u + ( u) T ), where µ represents the fluid dynamic viscosity. The flow equations are written in an Eulerian frame of reference, and thus a solution of these equations must be coupled with a methodology for following the deforming fluid-fluid interface. Here, the VoF volume tracking algorithm is implemented (Hirt and Nichols, 1981; Gueyffier et al., 1999; Afkhami and Bussmann, 2009). Volume tracking requires the introduction of a scalar function f, defined as f = 1 inside the drop and f = 0 outside the drop, for a two fluid system. Since f is passively advected with the flow, it satisfies the advection equation, Df/Dt = 0. Density and viscosity are then evaluated via volume-weighted formulae as ρ = ρ g + (ρ d ρ g )f and µ = µ g + (µ d µ g )f, respectively, where subscripts d and g refer to the drop and gas phases, respectively. Surface tension is reformulated as an equivalent JEnggMaths_KA-v2.tex; 10/06/2013; 18:32; p.3
4 4 Figure 1. Adaptive mesh: the mesh is adaptively refined around the contact line (solid line) and areas of high capillary pressure. Color contours depict the pressure distribution colored according to the reference pressure in the gas phase, with a maximum value colored in dark red a minimum value colored in dark blue. volume force as originally proposed in the continuum surface force (CSF) model (Brackbill et al., 1992) as, F st = σκδ S n, where σ is a constant interfacial tension and δ S denotes the Dirac delta function for the surface separating the fluids. The curvature κ and unit normal n are geometric characteristics of the surface and are described in terms of f and computed with a second-order height-function (HF) method (Sussman, 2003; Cummins et al., 2005; Afkhami and Bussmann, 2008; Francois et al., 2006). If partial slip is allowed at the contact line, we will impose the Navier-slip boundary condition at the (static) solid surface y = 0, (u, w) y=0 = λ (u, w)/ y y=0, where λ is the slip length. To investigate the effects of the relevant dimensionless parameter, we nondimensionalize the governing equations using a velocity scale V = σ/(ρl), where L is the appropriate length scale. The resulting time scale, T = ρl 3 /σ, is proportional to the time it takes for the capillary wave to propagate across the drop. With these choice of scales, the NS equations now can be written in nondimensional form D(ρu)/Dt = ( pi + Oh 1 τ) + F st. Oh = η/ ρσl is the Ohnesorge number, representing the ratio of the time it takes for momentum to diffuse across the droplet to the period of the capillary wave. 3. Numerical Model A complete description of the numerical methodologies can be found in (Popinet, 2003; Afkhami and Bussmann, 2008; Afkhami and Bussmann, 2009). Our in-house code is built based on an early version of the opensource flow solver Gerris (see (Popinet, 2003)). Briefly, the domain is JEnggMaths_KA-v2.tex; 10/06/2013; 18:32; p.4
5 5 (a) (b) (c) Figure 2. Initial structures: (a) Disk, (b) ring, and (c) triangular liquid nanostructures deposited on a substrate. Initial contact angle θ 0 = 90 (a) and (c), and 60 (b). discretized adaptively using the oct-tree finite volumes arranged hierarchically (see Fig. 1). Primitive variables are collocated at cell centers; normal velocities are also specified at faces, as in the MAC grid. The flow equations are discretized using a projection method based on a fractional-step scheme. In the fractional-step projection method, an interim velocity is computed first; this velocity field is projected onto a divergence-free velocity field, with the pressure field obtained as the solution of a Poisson equation. Advection terms are discretized using a second-order upwind scheme (Bell et al., 1989). As the face-centered velocities are exactly divergence-free, the volume fractions are advected using these velocities. The numerical algorithm for the imposition of the contact angle is described in details in (Afkhami and Bussmann, 2009). In this algorithm, the value of the contact angle affects the overall flow calculation in two ways: it defines the orientation of the VoF reconstruction in cells that contain the contact line and it influences the calculation of F st by affecting the curvature calculated in cells at and near the contact line to balance the surface tension forces on the contact line as in Young s relation (Young, 1805): σ cos θ eq = σ sg σ sl where σ sg and σ sl are interfacial tensions for the solid-gas and solidliquid surfaces, respectively, and θ eq is the equilibrium contact angle. If the actual contact angle is out of equilibrium, the resulting mismatch in curvature generates a force thus driving the contact line to a configuration that satisfies Young s relation Computational Setup We consider three geometrical configurations. Figure 2 shows the initial configurations for the case of a Cu disk (a), a Cu ring (b), and a Au triangular (c) nanostructures. We denote the disk geometry with an initial height h 0, and radius R 0 ; the ring geometry with an initial outer radius R out and inner radius R in ; the triangle geometry with an initial long side a, short side b, and thickness h 0 ; the initial contact angle is denoted as θ 0. The computational domain is L x L y L z and the JEnggMaths_KA-v2.tex; 10/06/2013; 18:32; p.5
6 6 maximum grid resolution of the adaptive mesh is denoted by. For all the simulations, an open boundary condition (pressure and velocity gradient equal zero) is imposed at the top and a symmetry boundary condition is imposed on the lateral sides. For nanodisk and nanoring liquid structures, a free-slip boundary condition is specified at the substrate in accordance with the MD simulations in (Fuentes-Cabrera et al., 2011a). For triangular nanostructures, the Navier-slip boundary condition with λ = 3 nm is imposed, so that the computed detachment velocities are consistent with the experimental results (Afkhami and Kondic, 2013; Habenicht et al., 2005). Table I. Overview of the sets of parameters used in the numerical simulations (Fuentes-Cabrera et al., 2011b; Habenicht et al., 2005). ρ(kg m 3 ) η(kg m 1 s 1 ) σ(kg s 2 ) l v(nm) Air Cu (@ 1500 K) Au (@ 1500 K) Results 4.1. Dewetting of Cu Nanostructures Nanodisks We first consider a Cu disk of h 0 = 15 Å, R 0 = 150 Å, and θ 0 = 90, when the equilibrium contact angle θeq = 80 is imposed. The nanodisk dewetts the substrate and collapses into a spherical cap with an equilibrium base radius Req = 8.65 nm corresponding to a contact angle of 80. We set L x = L y = L z = 375 Å. Numerical simulations are carried out on an adaptive mesh (see Fig. 1) with a maximum grid resolution of = 2.93 Å. We compute the retraction of the disk and consequent coalescence into a nanodroplet. We track the position of the edge of the retracting rim R(t) as a function of time. Figure 3 shows the results of direct simulations ( ). We also plot the simulation results with = 1.46 Å ( ) to show convergence with mesh refinement. After 140 ps (not shown in Fig. 3), the Cu nanodrop reaches an equilibrium state with R eq = 8.63 nm. The final stage is followed by the contact line oscillations about the equilibrium state; a signature of inertially-dominated convection. To gain additional insight regarding the effects that drive JEnggMaths_KA-v2.tex; 10/06/2013; 18:32; p.6
7 R(t) (nm) t (ps) Figure 3. Retraction of a Cu nanodisk on the substrate (h 0 = 15 Å, R 0 = 150 Å, θ 0 = 90 ). The front position, R(t), as a function of time for θeq = 80 ; direct numerical simulation = 2.93 Å ( ), and 1.46 Å ( ) for Oh The effect of increasing the Oh number on the retraction of the nanodisk is also plotted ( ) for Oh 1.1. oscillations, we have carried out additional set of simulations with modified inertial effects. This was implemented by varying incompressible fluid density, but more generally the dynamics can be characterized by Ohnesorge number, here defined as Oh = η/ ρσr 0. For the Cu nanodisk, Oh 0.35, suggesting that inertial effects are important; the same conclusion can be reached by considering an intrinsic length scale, l v = η 2 /(ρσ), above which inertial effects become significant; for liquid Cu, l v 1.78 nm, therefore smaller than typical length scales considered here. As an example, we also present the computational results for Oh 1.1 ( ) in Fig. 3. As shown, reduced inertial effects eliminate the final stage oscillations. Figure 4(a-c) shows sequences of snapshots of the simulations at various times for θeq = 80, 115, and 140. For θ eq = 80, viscous forces overcome the capillary force during dewetting until the drop assumes its equilibrium shape. For θ eq = 115, the nanodisk quickly dewets the surface and at about 40 ps the retracting rim has already collapsed into a droplet that begins to elongate in the y(vertical)-direction. The collapsing nanodrop then reaches to a maximum elongation at around 100 ps. After that, the droplet rebounds and settles at an equilibrium state with θ eq = 115 at around 250 ps. We also notice that the contraction in the x(horizontal)-direction and the elongation in the y-direction are now more pronounced compared to what was found for θ eq = 80. A different type of evolution occurs when θ eq is increased beyond 130. JEnggMaths_KA-v2.tex; 10/06/2013; 18:32; p.7
8 8 30 ps 100 ps 190 ps (a) (b) (c) (d) Figure 4. Snapshots of the Cu nanodisk with h 0 = 15 Å, R 0 = 150 Å, and θ 0 = 90. Direct numerical simulations of NS equations for θ eq = 80 (a), 115 (b), and 140 (c). MD simulations (d) employing a Lennard-Jones potential corresponding to θ eq = 138.5, from (Fuentes-Cabrera et al., 2011a). For θ eq = 140, the contraction is so fast that the nanodrop jumps off the surface, following elongation in the y-direction. This process begins with the accumulation of the liquid material towards the center of nanodroplet, leading to a vertical movement of the center of mass. The nanodroplet detaches from the substrate at about 100 ps. Figure 4(d) shows the results of MD simulations in (Fuentes-Cabrera et al., 2011a), when employing a Lennard-Jones potential that corresponds to θ eq = We find quantitative agreement with our continuum simulations. JEnggMaths_KA-v2.tex; 10/06/2013; 18:32; p.8
9 (a) (b) x (nm) x (nm) Figure 5. x-y cross sections of a Cu disk (h 0 = 15 Å, R 0 = 150 Å, θ eq = 140 ) at t = 0, 5, 10, 20, 30, 40, 110 ps with (a) Oh 0.35 and (b) Oh 1.1. The case (a) results in a drop that eventually detaches from the substrate as in Fig. 4(c), while for (b) the final outcome is a sessile droplet. Figure 5(a-b) further illustrates the relevance of inertial effects. Here, we compare the x-y cross sections for the case θ eq = 140 when Oh 0.35 and 1.1. As shown, reduced inertial effects eliminate the ejection of the nanodrop from the surface. For Oh 0.35 (Fig. 5(a)), a capillary ridge is formed at the initial stage and propagates from the nanodisk rim along the drop-gas interface progressing to the top of the droplet, causing a rapid rise in the droplet height and eventual ejection, while for Oh 1.1 (Fig. 5(b)), the capillary ridge is less pronounced, leading to a weaker collapse and no ejection. We also study the effect of the equilibrium contact angle, θ eq, on the retraction time and the upward velocity. Table II presents the variation of the retraction time and the upward velocity with θeq. As shown, the retraction dynamics is slower for smaller θ eq, and the upward velocity increases with increasing θ eq. These values are fully consistent with the values obtained in (Fuentes-Cabrera et al., 2011a) using MD simulations. It is also instructive to compare the velocity field for various θ eq. Figure 6 shows the velocity field at 20 ps for θ eq =140 and at 40 ps for θ eq =80 where the position of the edge of the rim is the same in both cases. As shown, the magnitude of the velocities are larger for θ eq =140 than of those for θ eq =80. In this figure, the background depicts the JEnggMaths_KA-v2.tex; 10/06/2013; 18:32; p.9
10 10 Table II. Effect of contact angle on the retraction time and the upward velocity. θ eq collapse time (ps) upward velocity (m s 1 ) (a) (b) Figure 6. Velocity vectors of a Cu nanodisk dewetting (h 0 = 15 Å, R 0 = 150 Å, θ 0 = 90 ) (a) at 20 ps for θ eq=80 and (b) at 40 ps for θ eq=140 ; the front has moved about the same distance in both cases. Color contours show the pressure field distribution colored according to the reference pressure in the gas phase, with a maximum value colored in dark red a minimum value colored in dark blue. adaptive mesh and the color contours are the pressure field distribution. Figure 7 shows the front velocity, Ṙ(t), at the early time as θ eq is varied. The results show a fast initial retraction followed by a slow contraction, indicating a transition from an inertially dominant regime to a viscous regime. The early-time inertial regime is regulated by the generation and propagation of the capillary wave initiated at the contact line (see e.g. Fig. 5(a)). Figure 7 also shows that the contact angle condition can influence the rapid inertial flow Nanorings We next consider the collapsing of Cu liquefied nanorings. When a ring of liquid is placed on a substrate, it retracts and collapses into a static droplet due to capillary effects (González et al., 2013). We show here the direct computational results for the dewetting of Cu liquefied nanorings as θ eq is varied. We consider the ring geometry with an initial outer radius R out = 150 nm, inner radius R in = 75 nm, and θ 0 = 60 (see Fig. 2). The radii and the θ 0 are chosen so that the drop volume remains almost unchanged from the previous nanodisk JEnggMaths_KA-v2.tex; 10/06/2013; 18:32; p.10
11 R (t) (m s -1 ) t (ps) Figure 7. Early-time retraction velocities of a Cu nanodisk dewetting (h 0 = 15 Å, R 0 = 150 Å, θ 0 = 90 ). θ eq = 80 ( ), 115 ( ), 130 ( ), and 140 ( ). case. The nanoring structure is out of equilibrium and the capillary pressures resulted from the inner and outer curvature imbalance drives an inward flow which eventually leads to a collapse of the ring into a static drop. We note that here Oh 0.4 calculated based on the average radius, i.e. (R out + R in )/2 = 75 nm. We set L x = L y = L z = 375 Å. Numerical simulations are carried out on an adaptive mesh (see Fig. 1) with a maximum grid resolution of = 2.93 Å. Figure 8 shows the time evolution of the nanorings for θ eq = 45, 60, 75, and 90. In particular, we are interested in the time, τ, it takes for a ring to collapse, defined as the time when the inner radius goes to zero. Figure 9 shows τ as a function of θ eq. We find that a larger θ eq leads to a smaller collapse time. We also observe that τ scales approximately linearly with θ eq Dewetting of Au Nanostructures Nanotriangles Recently, Habenicht et al. (Habenicht et al., 2005; Habenicht et al., 2008) experimentally examined whether it would be possible to begin with initially deformed droplets, in the form of solid flat triangular gold structures, on the surface and observe dewetting induced detachment. They considered Au triangular structures on the SiO 2 substrate. When melted by laser irradiation, the energetically unfavorable initial structure coalesces into a droplet. The transformation of surface energy to JEnggMaths_KA-v2.tex; 10/06/2013; 18:32; p.11
12 12 20 ps 45 ps 76 ps (a) 20 ps 45 ps 70 ps (b) 20 ps 45 ps 66 ps (c) 20 ps 45 ps 63 ps (d) Figure 8. Snapshots of the Cu nanoring with R out = 150 Å, R in = 75 Å, and θ 0 = 60, for θ eq = 45 (a), 60 (b), 75 (c), and 90 (d). The frames on very right show the drop at the time at which the collapse takes place, with inner radius going to zero. Oh 0.4. kinetic energy drives the motion and under appropriate conditions, detachment of the resulting droplets was observed. Figure 10 (Fig. 1(c) in (Habenicht et al., 2005)) shows sequences of images of the development of the dewetting process from flat triangles toward spherical structures. Here, we simulate the experiments of Habenicht et al. (2005). The initial triangle geometry has equal sides a = 405 nm and the initial contact angle θ 0 = 90 with θ eq = 140 (equilibrium contact angle of liquid Au on SiO 2 ). The computational domain is L x = L y = L z = 494 nm and the maximum grid = 1.5 nm. At the substrate, the Navierslip with slip length of 3 nm is imposed; this particular value is obtained JEnggMaths_KA-v2.tex; 10/06/2013; 18:32; p.12
13 τ (ps) θ eq Figure 9. The effect of θ eq on the collapse time, τ, for a Cu nanoring on the substrate (R out = 150 Å, R in = 75 Å, and θ 0 = 60 ). Oh 0.4. Figure 10. Fig. 1(c) in (Habenicht et al., 2005). Initially solid triangular gold structure with a thickness of 50 nm and a side length of 800 nm dewets on the substrate and coalesces into a droplet. θ eq = 131. by direct comparison with the experiments of Habenicht et al. (2005). Figure 11 represents the droplet velocity as a function of the initial thickness h 0. Our numerical results reproduce very well the experimental findings in (Habenicht et al., 2005). In particular, we note that both experiments and simulations show decreasing ejection velocity as the initial thickness is increased. This trend can be intuitively understood by realizing that the upward momentum is gained by the transformation of surface to kinetic energy. Thus, increasing the initial thickness results in a lower upward velocity by decreasing the surface to volume ratio. We observe that both the simulations and the experimental results show approximately linear dependence of the ejection velocity on h 0, with the simulation predicting zero ejection velocity for h nm. Figure 12 shows sequences of simulations of the dewetting of initially triangular nanostructures when initial thickness is varied. The consequent detachment of the droplet at various times is also demonstrated. As shown, the dewetting process begins at the corners, where JEnggMaths_KA-v2.tex; 10/06/2013; 18:32; p.13
14 14 Velocity (m s -1 ) h 0 (nm) Figure 11. Droplet detachment velocity for liquefied triangular Au nanostructures with a side length of 405 nm and various thicknesses, numerical results ( ) and experimental measurements ( ). Numerical results are obtained using the slip length of 3 nm. The numerical simulation predicts zero ejection velocity for h nm. the curvature is high. Due to a high surface tension force and large equilibrium contact angle, the liquid starts to accumulate at the corners. Vertices accumulate more mass than the straights sides. The humps at the corners then coalesce into a droplet. Finally, owing to a low viscosity of liquid gold (l v 0.9 nm), inertia dominates over viscous dissipation giving rise to an upward movement that leads to droplets detaching from the surface with velocities of tens of meters per second (see Fig. 11). As shown, increasing the initial thickness will result in a slower dynamics and an increase in the collapsing time. The time scale of the retraction process is consistent with the experiments of (Habenicht et al., 2005) measuring a detachment time scale of the order of 10 ns. Figure 13 shows the pressure distribution at the (substrate) x-z plane for (a) h 0 = 24 nm and (b) 47 nm. Solid line represents the line of contact and the background demonstrates the adaptive grid used. The high surface tension force due to high curvature increases the capillary pressure at the corners. This pressure gradient results in a transversal velocity driving the traveling rims. The time of which all the rims coalesce depends on the pressure gradient, explaining the difference of collapsing times shown in Fig. 12. As also shown, the sides of the triangle reach the center first, followed by vertices. This asynchronicity is the reason for excitation of contact line oscillations during the dewetting process, as illustrated in Fig. 13. Finally, we explore the effect of the initial shape. We consider an isosceles nanotriangle geometry with an initial long side a = 488 nm, short side b = 247 nm, h 0 = 24 nm, θ 0 = 90, and θ eq = 145. Figure JEnggMaths_KA-v2.tex; 10/06/2013; 18:32; p.14
15 ns 2 ns 5 ns (a) 1.2 ns 2 ns 7 ns (b) 1.2 ns 2 ns 11 ns (c) 1.2 ns 2 ns 13 ns (d) Figure 12. Snapshots of the Au equilateral nanotriangles with a = 405 nm, θ 0 = 90, and θ eq = 140, for h 0 = 24 nm (a), 47 nm (b), 95 nm (c), and 142 nm (d). The nanotriangular structure collapses into a droplet and detaches from the substrate (see Fig. 11 for corresponding ejection velocities). 14(a) shows the sequences of the numerical simulation. As shown, the dewetting behavior is different compared to dewetting dynamics of an equilateral triangle, although both processes lead to the ejection of the formed droplet from the substrate. Figure 14(b) shows the pressure distribution at x-z plane. It shows a higher capillary pressure at the vertex of the long edges due to a higher curvature. However, the triangle small edge vertices arrive faster at the center than the vertex of the longer sides. This mismatch excites the oscillation of the droplet translating into a tumbling movement after the droplet leaves the substrate. JEnggMaths_KA-v2.tex; 10/06/2013; 18:32; p.15
16 16 5. Conclusion We show that the direct numerical simulation based on Navier-Stokes equations can be utilized to study dewetting of liquefied metal nanostructures. In particular, we find that early-time dynamics of the dewetting of a Cu nanodisk is inertially dominated, and that viscous dissipation can not be expected to be the dominant resistance to the short-time dewetting. Furthermore, we demonstrate that the equilibrium contact angle, θ eq, plays a key role in the dewetting dynamics. For sufficiently large θ eq, the retraction and the consequent collapse into a nanodroplet is followed by the detachment of a nanodrop from the substrate with an upward velocity that is an increasing function of θ eq. Considering Cu nanorings, we show that θ eq influences the retraction velocity and the collapsing time. Our numerical results based on a VoF formulation is consistent with the results of the MD simulations offering a potential venue for modeling dewetting phenomena on nanoscale. We also directly simulate the dewetting process of initially triangular Au nanostrucutres that coalesce into nanodroplets which in turn detach from the surface. We demonstrate that the dynamics of the dewetting behavior is consistent with the experimental observations and that the computed detachment velocity and time are also in agreement with experimental measurements. We expect our findings to be of direct relevance to nanomaterial processes involving rapid dewetting. Acknowledgements The results presented here are a part of larger research project involving Kyle Mahady (NJIT), Javier Diez and Alejandro Gonzalez from UNCPBA (Argentina), and Jason Fowlkes, Miguel Fuentes-Cabrera and Philip Rack from ORNL and U. Tennessee. We acknowledge extensive discussion with all members of the research team. This work was supported in part by the NSF grant CBET References A. Habenicht, M. Olapinski, F. Burmeister, P. Leiderer, and J. Boneberg. Jumping Nanodroplets. Science, 309:2043, J. Boneberg, A. Habenicht, D. Benner, P. Leiderer, M. Trautvetter, C. Pfahler, A. Plettl, and P. Ziemann. Jumping nanodroplets: a new route towards metallic nano-particles. Appl. Phys. A, 93:415, N. A. Roberts, J. D. Fowlkes, K. Mahady, S. Afkhami, L. Kondic, and P. D. Rack. Directed Assembly of One- and Two-Dimensional Nanoparticle Arrays from Pulsed JEnggMaths_KA-v2.tex; 10/06/2013; 18:32; p.16
17 Laser Induced Dewetting of Square Waveforms. ACS Appl. Mater. Interfaces, 5:4450, M. Fuentes-Cabrera, B.H. Rhodes, J.D. Fowlkes, A. López-Benzanilla, H. Terrones, M.L. Simpson, and P.D. Rack. Molecular dynamics study of the dewetting of copper on graphite and graphene: Implications for nanoscale self-assembly. Phys. Rev. E, 83:041603, M. Fuentes-Cabrera, B.H. Rhodes, M.I. Baskes, H. Terrones, J.D. Fowlkes, M.L. Simpson, and P.D. Rack. Controlling the Velocity of Jumping Nanodroplets Via Their Initial Shape and Temperature. ACS Nano, 5:7130, P.G. de Gennes. Wetting: Statics and Dynamics. Rev. Mod. Phys., 57:827, P. J. Haley and M. J. Miksis. The effect of the contact line on droplet spreading. J. Fluid Mech., 223:57, J. Trice, D. Thomas, C. Favazza, R. Sureshkumar, and R. Kalyanaraman. Pulsedlaser-induced dewetting in nanoscopic metal films: Theory and experiments. Phys. Rev. B, 75:235439, V.S. Ajaev and D.A. Willis. Thermocapillary flow and rupture in films of molten metal on a substrate. Phys. Fluids, 15:3144, L. Kondic, J. Diez, P. Rack, Y. Guan, and J. Fowlkes. Nanoparticle assembly via the dewetting of patterned thin metal lines: Understanding the instability mechanism. Phys. Rev. E, 79:026302, S. Afkhami and M. Bussmann. Height functions for applying contact angles to 2D VOF simulations. Int. J. Numer. Meth. Fluids, 57:453, S. Afkhami and M. Bussmann. Height functions for applying contact angles to 3D VOF simulations. Int. J. Numer. Meth. Fluids, 61:827, S. Afkhami, S. Zaleski, and M. Bussmann. A Mesh-Dependent Model for Applying Dynamic Contact Angles to VOF Simulations. J. Comput. Phys., 228:5370, C. W. Hirt and B. D. Nichols. Volume of Fluid VOF Method for the Dynamics of Free Boundaries. J. Comput. Phys., 39:201, D. Gueyffier, J. Li, A. Nadim, R. Scardovelli, and S. Zaleski. Volume-of-fluid interface tracking and smoothed surface stress methods for three-dimensional flows. J. Comput. Phys., 152:423, J. U. Brackbill, D. B. Kothe, and C. Zemach. A continuum method for modeling surface tension. J. Comput. Phys., 100:335, M. Sussman. A second order coupled level set and volume-of-fluid method for computing growth and collapse of vapor bubbles. J. Comput. Phys, 187:110, S. J. Cummins, M. M. Francois, and D. B. Kothe. Estimating curvature from volume fractions. Comp. Struct., 83:425, M. M. Francois, S. J. Cummins, E. D. Dendy, D. B. Kothe, J. M. Sicilian, and M. W. Williams. A balanced-force algorithm for continuous and sharp interfacial surface tension models within a volume tracking framework. J. Comput. Phys., 213:141, S. Popinet. Gerris: a tree-based adaptive solver for the incompressible Euler equations in complex geometries. J. Comput. Phys., 190:572, J. B. Bell, P. Colella, and H. M. Glaz. A second-order projection method for the incompressible Navier-Stokes equations. J. Comput. Phys., 85(2):257, T. Young. An essay on the cohesion of fluids. Philos. Trans. R. Soc. London, 95, S. Afkhami and L. Kondic. Numerical simulation of ejected molten metalnanoparticles liquefied by laser irradiation: Interplay of geometry and dewetting. under revision: Phys. Rev. Lett., JEnggMaths_KA-v2.tex; 10/06/2013; 18:32; p.17
18 18 A. G. González, J. A. Diez, and L. Kondic. Stability of a liquid ring on a substrate. J. Fluid Mech., 718:246, JEnggMaths_KA-v2.tex; 10/06/2013; 18:32; p.18
19 ns 1.2 ns 2 ns 3 ns 4 ns 4.5 ns (a) 0.6 ns 1.2 ns 2 ns 3 ns 4 ns 5 ns (b) Figure 13. Contours of the pressure field distribution and the adaptive mesh at x-z plane for the Au equilateral nanotriangles with a = 405 nm, θ 0 = 90, and θ eq = 140. (a) h 0 = 24 nm and (b) 47 nm. The high capillary pressure at the corners are shown. The pressure distribution is colored according to the reference pressure in the gas phase, with a maximum value colored in dark red a minimum value colored in dark blue. JEnggMaths_KA-v2.tex; 10/06/2013; 18:32; p.19
20 ns 1.2 ns 2 ns 4.1 ns (a) Figure 14. (a) Snapshots of the Au isosceles nanotriangle structure with a = 488 nm, b = 247 nm, h 0 = 24 nm, θ 0 = 90, and θ eq = 145. The nanotriangular structure collapses into a droplet and detaches from the substrate at around 5 ns. (b) Contours of the pressure field distribution and the adaptive mesh at x-z plane. The pressure distribution is colored according to the reference pressure in the gas phase, with a maximum value colored in dark red a minimum value colored in dark blue. (b) JEnggMaths_KA-v2.tex; 10/06/2013; 18:32; p.20
On the dewetting of liquefied metal nanostructures
J Eng Math DOI 10.1007/s10665-013-9685-y On the dewetting of liquefied metal nanostructures Shahriar Afkhami Lou Kondic Received: 11 June 2013 / Accepted: 23 December 2013 Springer Science+Business Media
More informationReduction of parasitic currents in the DNS VOF code FS3D
M. Boger a J. Schlottke b C.-D. Munz a B. Weigand b Reduction of parasitic currents in the DNS VOF code FS3D Stuttgart, March 2010 a Institut für Aerodynamik und Gasdynamik, Universität Stuttgart, Pfaffenwaldring
More informationA semi-implicit finite volume implementation of the CSF method for treating surface tension in interfacial flows
A semi-implicit finite volume implementation of the CSF method for treating surface tension in interfacial flows M. Raessi, M. Bussmann*, and J. Mostaghimi Department of Mechanical and Industrial Engineering,
More informationHeight functions for applying contact angles to 2D VOF simulations
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids (2007) Published online in Wiley InterScience (www.interscience.wiley.com). DOI: 0.002/fld.65 Height functions for applying
More informationarxiv: v1 [physics.flu-dyn] 3 May 2018
Breakup of finite-size liquid filaments: Transition from no-breakup to breakup including substrate effects arxiv:1805.01558v1 [physics.flu-dyn] 3 May 2018 A. Dziedzic, 1 M. Nakrani, 1 B. Ezra, 1 M. Syed,
More informationOn the Breakup of Patterned Nanoscale Copper. Rings into Nanoparticles: Competing Instability and. Transport Mechanisms
On the Breakup of Patterned Nanoscale Copper Rings into Nanoparticles: Competing Instability and Transport Mechanisms Yueying Wu, Jason D. Fowlkes, Philip D. Rack,, Javier A. Diez, and Lou Kondic Department
More information, where the -function is equal to:
Paper ID ILASS08-000 ILASS08-9-4 ILASS 2008 Sep. 8-10, 2008, Como Lake, Italy BINARY COLLISION BETWEEN UNEQUAL SIZED DROPLETS. A NUMERICAL INVESTIGATION. N. Nikolopoulos 1, A. Theodorakakos 2 and G. Bergeles
More informationInstabilities in the Flow of Thin Liquid Films
Instabilities in the Flow of Thin Liquid Films Lou Kondic Department of Mathematical Sciences Center for Applied Mathematics and Statistics New Jersey Institute of Technology Presented at Annual Meeting
More informationIMPINGEMENT OF A DROPLET ONTO A DRY WALL: A NUMERICAL INVESTIGATION
IMPINGEMENT OF A DROPLET ONTO A DRY WALL:. Introduction A NUMERICAL INVESTIGATION N. Nikolopoulos, and G. Bergeles Department Mechanical Engineering Nat. Technical University of Athens 57 Zografos, Greece
More informationInvestigation of an implicit solver for the simulation of bubble oscillations using Basilisk
Investigation of an implicit solver for the simulation of bubble oscillations using Basilisk D. Fuster, and S. Popinet Sorbonne Universités, UPMC Univ Paris 6, CNRS, UMR 79 Institut Jean Le Rond d Alembert,
More informationOblique Drop Impact on Deep and Shallow Liquid
1 2 3 4 5 6 7 8 9 11 10 12 13 14 15 16 17 18 19 20 21 22 23 24 25 Commun. Comput. Phys. doi: 10.4208/cicp.XXX.XXX Oblique Drop Impact on Deep and Shallow Liquid B. Ray 1, G. Biswas 1,2, and A. Sharma 3
More informationA numerical analysis of drop impact on liquid film by using a level set method
Journal of Mechanical Science and Technology 25 (10) (2011) 2567~2572 wwwspringerlinkcom/content/1738-494x DOI 101007/s12206-011-0613-7 A numerical analysis of drop impact on liquid film by using a level
More informationA surfactant-conserving volume-of-fluid method for interfacial flows with insoluble surfactant
A surfactant-conserving volume-of-fluid method for interfacial flows with insoluble surfactant Ashley J. James Department of Aerospace Engineering and Mechanics, University of Minnesota John Lowengrub
More informationarxiv: v1 [physics.flu-dyn] 1 Mar 2017
Direct numerical simulation of variable surface tension flows using a Volume-of-Fluid method arxiv:1703.00327v1 [physics.flu-dyn] 1 Mar 2017 Ivana Seric, Shahriar Afkhami, Lou Kondic Department of Mathematical
More informationHeight function interface reconstruction algorithm for the simulation of boiling flows
Computational Methods in Multiphase Flow VI 69 Height function interface reconstruction algorithm for the simulation of boiling flows M. Magnini & B. Pulvirenti Dipartimento di Ingegneria Energetica, Nucleare
More informationApplication of the immersed boundary method to simulate flows inside and outside the nozzles
Application of the immersed boundary method to simulate flows inside and outside the nozzles E. Noël, A. Berlemont, J. Cousin 1, T. Ménard UMR 6614 - CORIA, Université et INSA de Rouen, France emeline.noel@coria.fr,
More informationUnderstanding Magnetic Field Gradient Effect From a Liquid Metal Droplet Movement
Donghong Gao* Neil B. Morley Vijay Dhir Department of Mechanical and Aerospace Engineering, University of California, Los Angeles, CA 90095 Understanding Magnetic Field Gradient Effect From a Liquid Metal
More informationDrop Impact on a Wet Surface: Computational Investigation of Gravity and Drop Shape
Drop Impact on a Wet Surface: Computational Investigation of Gravity and Drop Shape MURAT DINC and DONALD D. GRAY Department of Civil and Environmental Engineering West Virginia University P.O. Box 6103,
More informationCAPILLARY PHENOMENA ON A LIQUID SURFACE
45 CAPILLARY PHENOMENA ON A LIQUID SURFACE Mohammad Ali and Akira Umemura Department of Aerospace Engineering, Graduate School of Engineering, Nagoya University Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan
More informationNumerical simulations of drop impacts
Numerical simulations of drop impacts Christophe Josserand Institut D Alembert, CNRS-UPMC L. Duchemin, Z. Jian, P. Ray and S. Zaleski Numerical simulations of drop impacts Christophe Josserand Institut
More informationDirect Numerical Simulation of Single Bubble Rising in Viscous Stagnant Liquid
Direct Numerical Simulation of Single Bubble Rising in Viscous Stagnant Liquid Nima. Samkhaniani, Azar. Ajami, Mohammad Hassan. Kayhani, Ali. Sarreshteh Dari Abstract In this paper, direct numerical simulation
More informationApplying Contact Angle to a 2D Multiphase Smoothed Particle Hydrodynamics Model
Applying Contact Angle to a 2D Multiphase Smoothed Particle Hydrodynamics Model Amirsaman Farrokhpanah Graduate Research Assistant, Student Member of ASME Centre for Advanced Coating Technologies Mechanical
More informationThe effects of confinement and inertia on the production of droplets
Rheologica Acta manuscript No. (will be inserted by the editor) Y. Renardy The effects of confinement and inertia on the production of droplets Received: date / Accepted: date Abstract Recent experiments
More informationAbsorption of gas by a falling liquid film
Absorption of gas by a falling liquid film Christoph Albert Dieter Bothe Mathematical Modeling and Analysis Center of Smart Interfaces/ IRTG 1529 Darmstadt University of Technology 4th Japanese-German
More informationISCST shall not be responsible for statements or opinions contained in papers or printed in its publications.
Modeling of Drop Motion on Solid Surfaces with Wettability Gradients J. B. McLaughlin, Sp. S. Saravanan, N. Moumen, and R. S. Subramanian Department of Chemical Engineering Clarkson University Potsdam,
More informationFUNDAMENTAL STUDY OF BINGHAM FLUID BY MEANS OF DAM-BREAK FLOW MODEL
Annual Journal of Hydraulic Engineering, JSCE, Vol.54, 2010, February FUNDAMENTAL STUDY OF BINGHAM FLUID BY MEANS OF DAM-BREAK FLOW MODEL How Tion PUAY1 and Takashi HOSODA2 1 Member of JSCE, Phd Student,
More informationModeling the combined effect of surface roughness and shear rate on slip flow of simple fluids
Modeling the combined effect of surface roughness and shear rate on slip flow of simple fluids Anoosheh Niavarani and Nikolai Priezjev www.egr.msu.edu/~niavaran November 2009 A. Niavarani and N.V. Priezjev,
More informationNumerical Studies of Droplet Deformation and Break-up
ILASS Americas 14th Annual Conference on Liquid Atomization and Spray Systems, Dearborn, MI, May 2001 Numerical Studies of Droplet Deformation and Break-up B. T. Helenbrook Department of Mechanical and
More informationNumerical Solution of Integral Equations in Solidification and Melting with Spherical Symmetry
Numerical Solution of Integral Equations in Solidification and Melting with Spherical Symmetry V. S. Ajaev and J. Tausch 2 Southern Methodist University ajaev@smu.edu 2 Southern Methodist University tausch@smu.edu
More informationA unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation
A unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation S. Bordère a and J.-P. Caltagirone b a. CNRS, Univ. Bordeaux, ICMCB,
More informationAn OpenFOAM-based electro-hydrodynamical model
An OpenFOAM-based electro-hydrodynamical model Ivo Roghair, Dirk van den Ende, Frieder Mugele Department of Science and Technology, University of Twente, Enschede, The Netherlands Keywords: modelling,
More informationBoundary Conditions in Fluid Mechanics
Boundary Conditions in Fluid Mechanics R. Shankar Subramanian Department of Chemical and Biomolecular Engineering Clarkson University The governing equations for the velocity and pressure fields are partial
More informationDifferential criterion of a bubble collapse in viscous liquids
PHYSICAL REVIEW E VOLUME 60, NUMBER 1 JULY 1999 Differential criterion of a bubble collapse in viscous liquids Vladislav A. Bogoyavlenskiy* Low Temperature Physics Department, Moscow State University,
More informationSimulating Interfacial Tension of a Falling. Drop in a Moving Mesh Framework
Simulating Interfacial Tension of a Falling Drop in a Moving Mesh Framework Anja R. Paschedag a,, Blair Perot b a TU Berlin, Institute of Chemical Engineering, 10623 Berlin, Germany b University of Massachusetts,
More informationx j r i V i,j+1/2 r Ci,j Ui+1/2,j U i-1/2,j Vi,j-1/2
Merging of drops to form bamboo waves Yuriko Y. Renardy and Jie Li Department of Mathematics and ICAM Virginia Polytechnic Institute and State University Blacksburg, VA -, U.S.A. May, Abstract Topological
More informationChapter 10. Solids and Fluids
Chapter 10 Solids and Fluids Surface Tension Net force on molecule A is zero Pulled equally in all directions Net force on B is not zero No molecules above to act on it Pulled toward the center of the
More informationNumerical simulation of wave breaking in turbulent two-phase Couette flow
Center for Turbulence Research Annual Research Briefs 2012 171 Numerical simulation of wave breaking in turbulent two-phase Couette flow By D. Kim, A. Mani AND P. Moin 1. Motivation and objectives When
More informationLevel Set and Phase Field Methods: Application to Moving Interfaces and Two-Phase Fluid Flows
Level Set and Phase Field Methods: Application to Moving Interfaces and Two-Phase Fluid Flows Abstract Maged Ismail Claremont Graduate University Level Set and Phase Field methods are well-known interface-capturing
More informationMolecular Dynamics Study on the Binary Collision of Nanometer-Sized Droplets of Liquid Argon
Molecular Dynamics Study on the Binary Collision Bull. Korean Chem. Soc. 2011, Vol. 32, No. 6 2027 DOI 10.5012/bkcs.2011.32.6.2027 Molecular Dynamics Study on the Binary Collision of Nanometer-Sized Droplets
More informationCENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer
CENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer You are assigned to design a fallingcylinder viscometer to measure the viscosity of Newtonian liquids. A schematic
More informationDYNAMICS OF LIQUEFIED SEDIMENT FLOW. Advances in Natural and Technological Hazards Research Vol. 19
DYNAMICS OF LIQUEFIED SEDIMENT FLOW Advances in Natural and Technological Hazards Research Vol. 9 THE DYNAMICS OF LIQUEFIED SEDIMENT FLOW UNDERGOING PROGRESSIVE SOLIDIFICATION S. SASSA Disaster Prevention
More informationFluid Mechanics Theory I
Fluid Mechanics Theory I Last Class: 1. Introduction 2. MicroTAS or Lab on a Chip 3. Microfluidics Length Scale 4. Fundamentals 5. Different Aspects of Microfluidcs Today s Contents: 1. Introduction to
More informationDerivation of continuum models for the moving contact line problem based on thermodynamic principles. Abstract
Derivation of continuum models for the moving contact line problem based on thermodynamic principles Weiqing Ren Courant Institute of Mathematical Sciences, New York University, New York, NY 002, USA Weinan
More informationcontact line dynamics
contact line dynamics part 2: hydrodynamics dynamic contact angle? lubrication: Cox-Voinov theory maximum speed for instability corner shape? dimensional analysis: speed U position r viscosity η pressure
More informationFlow Field and Oscillation Frequency of a Rotating Liquid Droplet
Flow Field and Oscillation Frequency of a Rotating Liquid Droplet TADASHI WATANABE Center for Computational Science and e-systems Japan Atomic Energy Agency (JAEA) Tokai-mura, Naka-gun, Ibaraki-ken, 319-1195
More informationW ettability is a fundamental property of a solid surface, and plays important roles in daily life, industry
OPEN SUBJECT AREAS: COMPUTATIONAL METHODS BIOSURFACES Received 14 October 2013 Accepted 15 January 2014 Published 3 February 2014 Correspondence and requests for materials should be addressed to H.L. (lihuilmy@hotmail.
More informationSimulation of Droplet Ejection for a Piezoelectric Inkjet Printing Device
Materials Transactions, Vol. 45, No. 3 (2004) pp. 893 to 899 #2004 The Japan Institute of Metals Simulation of Droplet Ejection for a Piezoelectric Inkjet Printing Device Hsuan-Chung Wu 1, Huey-Jiuan Lin
More informationViscous non-linear theory of Richtmyer-Meshkov Instability. Abstract
Viscous non-linear theory of Richtmyer-Meshkov Instability Pierre Carles and Stéphane Popinet Laboratoire de Modélisation en Mécanique, Université Pierre et Marie Curie, Case 162, 4 place Jussieu, 75252
More informationarxiv: v4 [physics.flu-dyn] 27 Oct 2018
Transition in a numerical model of contact line dynamics and forced dewetting arxiv:1703.07038v4 [physics.flu-dyn] 27 Oct 2018 S. Afkhami 1a, J. Buongiorno b, A. Guion b, S. Popinet c, Y. Saade c, R. Scardovelli
More informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
More informationA numerical study on the effects of cavitation on orifice flow
PHSICS OF FLUIDS, A numerical study on the effects of cavitation on orifice flow S. Dabiri, W. A. Sirignano, and D. D. Joseph, University of California, Irvine, California 9697, USA University of Minnesota,
More informationTitle. Author(s)Nonomura, Makiko; Kobayashi, Ryo; Nishiura, Yasumasa. CitationJournal of the Physics Society Japan, 72(10): Issue Date
Title Periodic Precipitation during Droplet Evaporation on AuthorsNonomura, Makiko; Kobayashi, Ryo; Nishiura, Yasumasa CitationJournal of the Physics Society Japan, 721: 268-2 Issue Date 23-1 Doc URL http://hdl.handle.net/2115/6
More informationAn Introduction to Theories of Turbulence. James Glimm Stony Brook University
An Introduction to Theories of Turbulence James Glimm Stony Brook University Topics not included (recent papers/theses, open for discussion during this visit) 1. Turbulent combustion 2. Turbulent mixing
More informationGRAVITY-DRIVEN MOTION OF A SWARM OF BUBBLES IN A VERTICAL PIPE
27th International Conference on Parallel Computational Fluid Dynamics Parallel CFD2015 GRAVITY-DRIVEN MOTION OF A SWARM OF BUBBLES IN A VERTICAL PIPE Néstor Balcázar, Oriol Lehmkuhl, Jesús Castro, Joaquim
More informationDirect Numerical Simulations of Gas-Liquid Flows
Direct Numerical Simulations of Gas-Liquid Flows 1 Gretar Tryggvason*; 1 Jiacai Lu; 2 Ming Ma 1 Johns Hopkins University, Baltimore, MD, USA; 2 University of Notre Dame, Notre Dame, IN, USA Introduction
More informationA FEM STUDY ON THE INFLUENCE OF THE GEOMETRIC CHARACTERISTICS OF METALLIC FILMS IRRADIATED BY NANOSECOND LASER PULSES
8 th GRACM International Congress on Computational Mechanics Volos, 12 July 15 July 2015 A FEM STUDY ON THE INFLUENCE OF THE GEOMETRIC CHARACTERISTICS OF METALLIC FILMS IRRADIATED BY NANOSECOND LASER PULSES
More informationNumerical study of turbulent two-phase Couette flow
Center for Turbulence Research Annual Research Briefs 2014 41 Numerical study of turbulent two-phase Couette flow By A.Y. Ovsyannikov, D. Kim, A. Mani AND P. Moin 1. Motivation and objectives The motivation
More informationNumerical simulation of breakup of a viscous drop in simple shear flow through a volume-of-fluid method
Numerical simulation of breakup of a viscous drop in simple shear flow through a volume-of-fluid method Jie Li, Yuriko Y. Renardy, and Michael Renardy Citation: Physics of Fluids (1994-present) 12, 269
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
More informationA Discontinuous Galerkin Conservative Level Set Scheme for Simulating Turbulent Primary Atomization
ILASS-Americas 3rd Annual Conference on Liquid Atomization and Spray Systems, Ventura, CA, May 011 A Discontinuous Galerkin Conservative Level Set Scheme for Simulating Turbulent Primary Atomization M.
More informationModel Studies on Slag-Metal Entrainment in Gas Stirred Ladles
Model Studies on Slag-Metal Entrainment in Gas Stirred Ladles Anand Senguttuvan Supervisor Gordon A Irons 1 Approach to Simulate Slag Metal Entrainment using Computational Fluid Dynamics Introduction &
More informationNonlinear oscillations and rotations of a liquid droplet
Nonlinear oscillations and rotations of a liquid droplet Tadashi Watanabe watanabe.tadashi66@jaea.go.jp Abstract Nonlinear oscillations and rotations of a liquid droplet are simulated numerically by solving
More informationDynamics of Transient Liquid Injection:
Dynamics of Transient Liquid Injection: K-H instability, vorticity dynamics, R-T instability, capillary action, and cavitation William A. Sirignano University of California, Irvine -- Round liquid columns
More informationSIMULATION OF A 2D GRANULAR COLUMN COLLAPSE ON A RIGID BED
1 SIMULATION OF A 2D GRANULAR COLUMN COLLAPSE ON A RIGID BED WITH LATERAL FRICTIONAL EFFECTS High slope results and comparison with experimental data Nathan Martin1, Ioan Ionescu2, Anne Mangeney1,3 François
More informationSupplementary Information for Engineering and Analysis of Surface Interactions in a Microfluidic Herringbone Micromixer
Supplementary Information for Engineering and Analysis of Surface Interactions in a Microfluidic Herringbone Micromixer Thomas P. Forbes and Jason G. Kralj National Institute of Standards and Technology,
More informationLIQUID FILM THICKNESS OF OSCILLATING FLOW IN A MICRO TUBE
Proceedings of the ASME/JSME 2011 8th Thermal Engineering Joint Conference AJTEC2011 March 13-17, 2011, Honolulu, Hawaii, USA AJTEC2011-44190 LIQUID FILM THICKNESS OF OSCILLATING FLOW IN A MICRO TUBE Youngbae
More information12.1 Viscous potential flow (VPF)
1 Energy equation for irrotational theories of gas-liquid flow:: viscous potential flow (VPF), viscous potential flow with pressure correction (VCVPF), dissipation method (DM) 1.1 Viscous potential flow
More informationCharacterization of Low Weber Number Post-Impact Drop-Spread. Dynamics by a Damped Harmonic System Model
Characterization of Low Weber Number Post-Impact Drop-Spread Dynamics by a Damped Harmonic System Model A Thesis Submitted to the Division of Research and Advanced Studies of the University of Cincinnati
More informationNumerical Analysis of Metal Transfer in Gas Metal Arc Welding
Numerical Analysis of Metal Transfer in Gas Metal Arc Welding G. WANG, P.G. HUANG, and Y.M. ZHANG The present article describes a numerical procedure to simulate metal transfer and the model will be used
More informationOn the breakup of fluid films of finite and infinite extent
PHYSICS OF FLUIDS 19, 072107 2007 On the breakup of fluid films of finite and infinite extent Javier A. Diez a Instituto de Física Arroyo Seco, Universidad Nacional del Centro de la Provincia de Buenos
More informationDriven large contact angle droplets on chemically heterogeneous substrates
October 2012 EPL, 100 (2012) 16002 doi: 10.1209/0295-5075/100/16002 www.epljournal.org Driven large contact angle droplets on chemically heterogeneous substrates D. Herde 1,U.Thiele 2,S.Herminghaus 1 and
More informationModule 3: "Thin Film Hydrodynamics" Lecture 11: "" The Lecture Contains: Micro and Nano Scale Hydrodynamics with and without Free Surfaces
The Lecture Contains: Micro and Nano Scale Hydrodynamics with and without Free Surfaces Order of Magnitude Analysis file:///e /courses/colloid_interface_science/lecture11/11_1.htm[6/16/2012 1:39:56 PM]
More information1.3 Molecular Level Presentation
1.3.1 Introduction A molecule is the smallest chemical unit of a substance that is capable of stable, independent existence. Not all substances are composed of molecules. Some substances are composed of
More informationFinite Element Model of a complex Glass Forming Process as a Tool for Control Optimization
Excerpt from the Proceedings of the COMSOL Conference 29 Milan Finite Element Model of a complex Glass Forming Process as a Tool for Control Optimization Felix Sawo and Thomas Bernard Fraunhofer Institute
More informationNonlinear oscillations and collapse of elongated bubbles subject to weak viscous effects
PHYSICS OF FLUIDS 17, 102101 2005 Nonlinear oscillations and collapse of elongated bubbles subject to weak viscous effects Kostas Tsiglifis a and Nikos A. Pelekasis b Department of Mechanical and Industrial
More informationInstabilities in the flow of thin films on heterogeneous surfaces
PHYSICS OF FLUIDS VOLUME 16, NUMBER 9 SEPTEMBER 2004 Instabilities in the flow of thin films on heterogeneous surfaces Lou Kondic a) Department of Mathematical Sciences & Center for Applied Mathematics
More informationImplementation of a symmetry-preserving discretization in Gerris
Implementation of a symmetry-preserving discretization in Gerris Daniel Fuster Cols: Pierre Sagaut, Stephane Popinet Université Pierre et Marie Curie, Institut Jean Le Rond D Alembert Introduction 10/11:
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationSimulation numérique de l atomisation. Institut Jean Le Rond d Alembert, CNRS & Université Pierre et Marie Curie -- UPMC Paris 6
Simulation numérique de l atomisation Stéphane Zaleski & Yue Ling Institut Jean Le Rond d Alembert, CNRS & Université Pierre et Marie Curie -- UPMC Paris 6 web site http://www.ida.upmc.fr/~zaleski 1/79
More informationA charge-conservative approach for simulating electrohydrodynamic two-phase flows using Volume-Of-Fluid
A charge-conservative approach for simulating electrohydrodynamic two-phase flows using Volume-Of-Fluid J. M. López-Herrera,, S. Popinet b, M. A. Herrada a Dept. Ingeniería Aerospacial y Mecánica de Fluidos,
More informationMULTI-SCALE VOLUME OF FLUID MODELLING OF DROPLET COALESCENCE
Ninth International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 10-12 December 2012 MULTI-SCALE VOLUME OF FLUID MODELLING OF DROPLET COALESCENCE Lachlan R. MASON,
More informationA comparison of viscoelastic stress wakes for 2D and 3D Newtonian drop deformation in a viscoelastic matrix under shear
A comparison of viscoelastic stress wakes for 2D and 3D Newtonian drop deformation in a viscoelastic matrix under shear S. Afkhami, P. Yue, and Y. Renardy Department of Mathematics, 460 McBryde Hall, Virginia
More informationHydrodynamics of wetting phenomena. Jacco Snoeijer PHYSICS OF FLUIDS
Hydrodynamics of wetting phenomena Jacco Snoeijer PHYSICS OF FLUIDS Outline 1. Creeping flow: hydrodynamics at low Reynolds numbers (2 hrs) 2. Thin films and lubrication flows (3 hrs + problem session
More informationSupplementary table I. Table of contact angles of the different solutions on the surfaces used here. Supplementary Notes
1 Supplementary Figure 1. Sketch of the experimental setup (not to scale) : it consists of a thin mylar sheet (0, 02 4 3cm 3 ) held fixed vertically. The spacing y 0 between the glass plate and the upper
More informationEffect of density ratio on the secondary breakup: A numerical study
ICLASS 2018, 14 th Triennial International Conference on Liquid Atomization and Spray Systems, Chicago, IL, USA, July 22-26, 2018 Effect of density ratio on the secondary breakup: A numerical study Suhas.
More information7 The Navier-Stokes Equations
18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and
More informationFEM-Level Set Techniques for Multiphase Flow --- Some recent results
FEM-Level Set Techniques for Multiphase Flow --- Some recent results ENUMATH09, Uppsala Stefan Turek, Otto Mierka, Dmitri Kuzmin, Shuren Hysing Institut für Angewandte Mathematik, TU Dortmund http://www.mathematik.tu-dortmund.de/ls3
More informationFLOW ASSURANCE: DROP COALESCENCE IN THE PRESENCE OF SURFACTANTS
FLOW ASSURANCE: DROP COALESCENCE IN THE PRESENCE OF SURFACTANTS Vishrut Garg and Osman A. Basaran Davidson School of Chemical Engineering Purdue University With special thanks to: Krish Sambath (now at
More informationMaximal deformation of an impacting drop
J. Fluid Mech. (24), vol. 57, pp. 99 28. c 24 Cambridge University Press DOI:.7/S222494 Printed in the United Kingdom 99 Maximal deformation of an impacting drop By CHRISTOPHE CLANET, CÉDRIC BÉGUIN, DENIS
More informationAn Overview of Fluid Animation. Christopher Batty March 11, 2014
An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.
More informationNumerical Simulation of Film Flow over an Inclined Plate: Effects of Solvent Properties and Contact Angle
Numerical Simulation of Film Flow over an Inclined Plate: Effects of Solvent Properties and Contact Angle Janine Carney and Rajesh Singh Multiphase Flow Science Workshop August 5-6, 214 Lakeview Golf Resort
More informationDetailed 3D modelling of mass transfer processes in two phase flows with dynamic interfaces
Detailed 3D modelling of mass transfer processes in two phase flows with dynamic interfaces D. Darmana, N.G. Deen, J.A.M. Kuipers Fundamentals of Chemical Reaction Engineering, Faculty of Science and Technology,
More information- Marine Hydrodynamics. Lecture 4. Knowns Equations # Unknowns # (conservation of mass) (conservation of momentum)
2.20 - Marine Hydrodynamics, Spring 2005 Lecture 4 2.20 - Marine Hydrodynamics Lecture 4 Introduction Governing Equations so far: Knowns Equations # Unknowns # density ρ( x, t) Continuity 1 velocities
More informationInterpreting Differential Equations of Transport Phenomena
Interpreting Differential Equations of Transport Phenomena There are a number of techniques generally useful in interpreting and simplifying the mathematical description of physical problems. Here we introduce
More informationThe Role of Splatting Effect in High Schmidt Number Turbulent Mass Transfer Across an Air-Water Interface
Turbulence, Heat and Mass Transfer 4 K. Hanjalic, Y. Nagano and M. Tummers (Editors) 3 Begell House, Inc. The Role of Splatting Effect in High Schmidt Number Turbulent Mass Transfer Across an Air-Water
More informationPHYSICS OF FLUID SPREADING ON ROUGH SURFACES
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 5, Supp, Pages 85 92 c 2008 Institute for Scientific Computing and Information PHYSICS OF FLUID SPREADING ON ROUGH SURFACES K. M. HAY AND
More informationFluid Animation. Christopher Batty November 17, 2011
Fluid Animation Christopher Batty November 17, 2011 What distinguishes fluids? What distinguishes fluids? No preferred shape Always flows when force is applied Deforms to fit its container Internal forces
More informationNUMERICAL INVESTIGATION OF THERMOCAPILLARY INDUCED MOTION OF A LIQUID SLUG IN A CAPILLARY TUBE
Proceedings of the Asian Conference on Thermal Sciences 2017, 1st ACTS March 26-30, 2017, Jeju Island, Korea ACTS-P00786 NUMERICAL INVESTIGATION OF THERMOCAPILLARY INDUCED MOTION OF A LIQUID SLUG IN A
More informationCHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE
CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE In this chapter, the governing equations for the proposed numerical model with discretisation methods are presented. Spiral
More information