Gravitational effects on the deformation of a droplet adhering to a horizontal solid surface in shear flow
|
|
- Sylvia Mason
- 5 years ago
- Views:
Transcription
1 PHYSICS OF FLUIDS 19, Gravitational effects on the deformation of a droplet adhering to a horizontal solid surface in shear flow P. Dimitrakopoulos Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, Maryland 20742, USA Received 12 February 2007; accepted 30 October 2007; published online 19 December 2007 In this paper we investigate the gravitational effects on the deformation of a three-dimensional droplet adhering to a horizontal rough solid surface in steady shear Stokes flows. Our study considers both positive and negative Bond numbers for viscous and inviscid droplets. When the interfacial system is initially at hydrostatic equilibrium, our study shows that the Bond number affects the deformation of viscous droplets with moderate and large initial contact angles in a different way than those for small angles owing to the interplay between the viscous and surface tension forces. Inviscid droplets with different initial contact angles show similar behavior as the Bond number increases, i.e., their deformation is monotonically decreased owing to the monotonic decrease of the droplets height and thus the exerted pressure force. Our study identifies the gravitational effects of the onset of interfacial sliding, i.e., on the portions of the contact line which slide first due to violation of the hysteresis condition. When the interfacial system is not at hydrostatic equilibrium at the flow initiation, its dynamic evolution is more complicated owing to the combined action of the shear flow with the gravitational forcing due to the difference between the initial shape with the hydrostatic one. Our computational results are accompanied with an analysis of the forces on the droplet which provides physical insight and identifies the three-dimensional nature of the interfacial deformation American Institute of Physics. DOI: / I. INTRODUCTION The deformation, sliding, and dislodging of fluid droplets adherent to solid substrates show applications in numerous areas including distillation, spray coating, packed towers, and a variety of multiphase flow operations in the chemical process industry. Our interest in the problem focuses on viscous flows at low Reynolds number. This regime has relevance in coating operations, enhanced oil recovery, microfluidics, and biological systems. 1 4 The main characteristic of the interfacial attachment on rough and chemically inhomogeneous solid surfaces is the hysteresis effect where the contact line remains stationary as long as the contact angle satisfies the requirement R A 1 anywhere along the contact line. The advancing and receding angles, A and R, are properties of the specific interfacial system. 5,6 As a result, the entire process under external forcing i.e., a shearing flow may exhibit several stages including interfacial deformation with a stationary contact line, initial sliding with partial rearrangement of the contact line, and the final drop dislodging, depending on the strength of the external forcing. 7 Several studies have considered the drop deformation under fixed constant line conditions for two- and threedimensional interfaces at both low- and finite-reynolds numbers, e.g., Refs Recently, Dimitrakopoulos 7 also considered the deformation of three-dimensional attached droplets with fixed circular contact lines but focused his interest on the onset of interfacial sliding under contact angle hysteresis. The fundamental issues associated with the yield criteria for drop displacement from rigid boundaries i.e., the equilibrium conditions just before the final drop dislodging have been addressed in a series of papers by Dussan and co-workers, and by Dimitrakopoulos and Higdon. 5,6,16 18 The flow-induced droplet sliding has been considered by several computational studies mainly for systems with negligible contact angle hysteresis where even the slightest external forcing is able to cause interfacial sliding, e.g., Refs. 10, 12, 19, and 20. We emphasize that the aforementioned studies commonly ignored the gravitational effects on the drop deformation and sliding by considering fluids with the same density or very small droplets where the Bond number can be neglected due to its stronger length dependence compared to that for the capillary number as seen in Eq. 2 below. However, for many systems the gravitational effects cannot be neglected even for relatively small droplets. For example, a 10 l drop having the surface tension of water and a density difference equal to the density of water, shows a Bond number B d ±0.31. The corresponding Bond numbers for 50 and 100 l drops are B d ±0.91, We note that the aforementioned drop volumes exhibit a radius of O 1 mm when they acquire a spherical shape. The hydrostatic shapes of these droplets attached to horizontal solid surfaces are thus much deformed compared to the sphericalcap shape of zero-bond-number droplets. Therefore, the flow-induced deformation of such droplets may differ considerably from that for negligible Bond number. In the present study we consider the gravitational effects /2007/19 12 /122105/14/$ , American Institute of Physics
2 P. Dimitrakopoulos Phys. Fluids 19, on the dynamics of a three-dimensional droplet attached to a horizontal rough plane wall with a stationary contact line in the presence of a steady shear flow at low Reynolds number. The initial shape of the drop is axisymmetric and thus, at the flow initiation, the drop forms a constant contact angle 0 with R 0 A around a circular contact line. To study this problem, we employ our interfacial spectral boundary element algorithms described briefly in Sec. II. In particular, for the efficient determination of the droplet shape as a function of time, we employ our fully implicit time integration algorithm, while for the determination of the equilibrium deformation in a wide range of subcritical flow rates, we utilize our Newton method for equilibrium interfaces in Stokes flow. For the majority of our computations presented in Sec. III, the drops are at hydrostatic equilibrium for the desired Bond number at the instance of the flow initiation. This physical system may represent a droplet which has been attached to a solid surface under quiescent conditions, formed a stationary circular contact line, and acquired its hydrostatic shape for a given Bond number without violating the hysteresis condition, Eq. 1. Moreover, in Sec. III F we consider the case where a droplet may be attached to a solid surface without being at hydrostatic equilibrium at the flow initiation. Our work includes both positive and negative Bond numbers for viscous and inviscid droplets, and identifies the gravitational effects on the drop deformation and dynamics for different initial configurations i.e., different angles 0. Our computational results are accompanied with an analysis of the forces on the droplet interface which provide useful physical insight on the deformation behavior. We also identify the gravitational effects of the onset of interfacial sliding, i.e., on the portions of the contact line which slide first due to violation of the hysteresis condition, Eq. 1, in strong enough shear flows. As discussed in Sec. IV, our study identifies the three-dimensional nature of the deformation and sliding of adherent droplets owing to the ability of the drops to expand or contract in the cross-flow direction changing the interfacial dynamics. Our interest lies on the interfacial deformation under a stationary contact line and the onset of interfacial sliding that occur at low and moderate flow rates. Note that we do not consider interfacial breakup which may occur at high enough flow rates. By considering representative values for small, moderate, and large initial angles and via the associated scaling analysis, the physical understanding presented in this work covers a wide range of contact angles for low and moderate flow rates. As mentioned earlier, the physical system considered in this paper corresponds to the initial stage of the deformation/ displacement process where the droplet deforms but its contact line remains stationary. This system provides physical insight for the droplet deformation and dynamics which is useful in processes where we want to keep the droplets attached to a solid surface. In addition, our work provides information on the onset of interfacial sliding. We emphasize that, for the problem studied in this paper, the onset of interfacial sliding cannot be inferred from the knowledge of the equilibrium contact angle distribution at gradually increasing flow rates. In particular, the interfacial sliding may occur during the transient evolution, depending on the flow rate, the physical parameters of the material system i.e., R, A, the viscosity ratio and the Bond number as well as the initial angle 0. As we discuss in Sec. III, the onset of interfacial sliding is strongly affected by the drop s initial configuration and the Bond number; this can be exploited to initiate interfacial sliding at a desired portion of the contact line. II. MATHEMATICAL FORMULATION AND COMPUTATIONAL ALGORITHM We consider a three-dimensional droplet attached to a horizontal plane solid wall surrounded by a viscous fluid as shown later in Fig. 3 c. The droplet fluid 1 has density 1 and viscosity, while the surrounding fluid fluid 2 has density 2 and viscosity. The droplet size is specified by its volume V or equivalently by the radius a of a spherical droplet of volume 4 a 3 /3=V. The gravitational acceleration is g while the surface tension is assumed constant. At time t =0, a steady shear flow u = Gz,0,0, where G is the shear rate, is introduced into the system causing interfacial deformation. In our study the time is scaled with the flow time scale f =G 1. The dimensionless parameters which represent the strength of the viscous flow forces and the gravitational forces, with respect to the interfacial forces, are the capillary number Ca and Bond number B d, respectively, defined by Ca = Ga, B d = 1 2 ga 2. 2 Note that the length scale a used in the definition of these parameters is based on the droplet volume V since we care to identify the deformation of a specific fluid volume V which may be attached to a solid surface in many possible initial configurations and undergo simple shear flow at different flow rates or Ca. For the majority of our computations, before the initiation of the steady shear flow, the fluids are at hydrostatic equilibrium under a specific gravitational influence or Bond number B d while the drop forms a constant contact angle = 0 measured from within the drop phase around a circular contact line. In Sec. III F we also consider the case where initially the fluids may not be at hydrostatic equilibrium for a desired Bond number. Assuming low-reynolds-number conditions for both fluid phases, the flow over the drop interface S d may be described by the boundary integral formula, u x 0 u x 0 = S 2 S f 2 f T u 2 u n ds + S 1 S f 1 T u 1 n ds S d S f f T 1 u u n ds, 3
3 Gravitational effects on the deformation of a droplet Phys. Fluids 19, where S 1 is the portion of the solid surface wetted by the drop and S 2 is the remaining solid surface. 6 Note that S is the fundamental solution for the three-dimensional Stokes equations and T is the associated stress, while n is the unit normal which we choose to point into fluid 2. The coefficient takes values 4 1+, 4, and 4 for points x 0 on the surfaces S d, S 1, and S 2, respectively, while takes the value 4 for points x 0 on the surfaces S d and S 2, and zero for points on the surface S 1. It is of interest to note that the requirement for low- Reynolds number in the interior flow is not formally valid for very small viscosity ratios. For low-viscosity/inviscid droplets 1, our computations require that the inner shear stress and pressure change are very small compared to the exterior stresses. As discussed by Dimitrakopoulos and Higdon, 17 for 1 this condition remains valid for arbitrary interior Reynolds number, as long as the exterior Reynolds number is very small. See the Appendix in Ref. 17. At the interface, the boundary conditions on the velocity u and surface stress f are u 1 = u 2 = u, f = f 2 f 1 = n n g x n, where the subscripts designate quantities evaluated in fluids 1 and 2, respectively. The interfacial shape x evolves over time based on the kinematic condition at the interface, dx = u n n. dt 6 In order to determine the droplet shape as a function of time, we employ our fully implicit spectral boundary element algorithm for interfacial dynamics in Stokes flow. 21 Our method is based on a mathematically rigorous combination of implicit schemes with our Jacobian-free three-dimensional Newton method, 6 and thus it has strong stability properties which permit the utilization of large time steps, independent of the space discretization. In essence, our implicit interfacial algorithm removes the penalty of a large number of time steps required to monitor the deformation of the fluid interface associated with the commonly employed explicit timeintegration interfacial methodologies. For the efficient determination of the equilibrium deformation in a wide range of subcritical flow rates, we utilize our Newton method for equilibrium interfaces under flow conditions. 6 The numerical solution of our algorithms is achieved through an extension of the spectral boundary element method, 6,22 and thus our methodologies achieve exponential convergence in the numerical accuracy of determining the interfacial shape, without being affected by the disadvantage of the spectral methods for partial differential equations to create denser systems. The interested reader is referred to our relevant papers for more details on the fully implicit and the Newton interfacial methods. 6,21 The adherent droplets are described by employing a discretization of N E =24 quadrilateral spectral elements. The surface of the drop is projected onto a cube whose faces are 4 5 subdivided into a total of 11 elements as shown in Figs. 3 c and 4 c below. The wetted area S 1 on the solid surface is discretized into five elements while the outer solid surface S 2 is discretized into two rows of four elements each. In our computations, the outer solid surface which formally should extend to infinity covers an area of at least 10 times the radius of the contact line which is sufficient to produce negligible error in all cases. Our computational results presented in Secs. III A III E were derived by employing N B =7 9 basis points along each curvilinear direction of the spectral elements. In our fully implicit algorithm we utilize the third-order diagonally implicit Runge-Kutta scheme. For low and moderate viscosity ratio, 1, we mainly employed a time step t=0.1; for droplets with high viscosity ratio, =10, a time step t=1 was used. For our results presented in Sec. III F, which show large variations in short time periods and the contact line tends rapidly to be tangential to the solid surface, we employed smaller time steps i.e., t=0.02,0.01 and N B =8 10 basis points so that we have sufficient space discretization for accurate solution. In Sec. III F we also employed a multiple time-step approach as described in the caption of Figs. 17 and 19, shown below. The accuracy of our results was verified by employing smaller time steps and different basis points for several representative cases; these comparisons show that our results are accurate to at least three significant digits, and thus if we include them on the same plot, they show no difference. The three-dimensional droplet shapes presented here were derived from the actual spectral grid by spectrally interpolating to N B =20. The present problem depends on the Bond number B d, the capillary number Ca, the viscosity ratio, and the initial angle 0. Since our results for different low and moderate flow rates Ca show similar behavior, we restrict our presentation to Ca=0.10 so that we are able to identify the influence of the rest parameters. An exception is Sec. III E where we collect our results on the equilibrium deformation of viscous and inviscid droplets for a wide range of low and moderate subcritical flow rates. III. RESULTS A. Quiescent drop shape under gravity The first step in our computations was to determine the axisymmetric droplet shape for a given initial angle 0 with R 0 A and Bond number B d 0 under quiescent conditions, i.e., Ca= 0. This physical system may represent a droplet which has been attached on a solid surface under quiescent conditions, formed a stationary circular contact line, and acquired its hydrostatic shape for a given Bond number without violating the hysteresis condition, Eq. 1. To achieve this, we employed our Newton method for equilibrium interfaces under Stokes flow conditions. 6 The initial shape was that of a spherical cap with the desired initial angle 0, and by specifying Ca=0 we let our Newton method determine the quiescent shape under different positive and negative Bond numbers by seeking the optimal contact line shape under these conditions as described in Sec. 2.3 of our earlier study.
4 P. Dimitrakopoulos Phys. Fluids 19, FIG. 2. Variation of the drop surface S d scaled with its value for B d =0 with the Bond number, under quiescent conditions for 0 =50,90,120. FIG. 1. Profiles of an adherent droplet under quiescent conditions for a 0 =50 and B d = 0.5,0,0.5; b 0 =90 and B d = 0.5,0,0.5,1,1.5,2; c 0 =120 and B d = 0.3,0,0.5,1,1.5. Figure 1 shows the quiescent droplet profile i.e., theintersection of droplet surface with the plane y=0 for 0 =50, 90, 120 and several Bond numbers. As expected, a negative Bond number results in an upward-extended droplet with larger height and smaller length/width compared to those for B d =0. The opposite happens for positive Bond numbers where the flattened drops show a larger wetted area. Beyond the obvious changes on the drop s height and width, the Bond number also affects the surface area S d of the quiescent shape in a more complicated way as shown in Fig. 2. In particular, for 0 =90 the spherical cap corresponding to B d =0 shows the minimum surface area S d as expected. For smaller initial angles, S d increases monotonically with the Bond number while the opposite happens for larger initial angles. We note that the variation in the droplet s dimensions and surface area with the Bond number provides information on the forces acting on the interface as we discuss later. B. Droplets with 0 =90 and =1 We begin our investigation by considering the effects of Bond number on the deformation of adherent droplets with =1 and 0 =90. Figures 3 and 4 show the transient evolution of an adherent droplet for a negative B d = 0.5 and a positive B d =2 Bond number, respectively. In both cases, in the upstream portion of the interface, the contact angle is reduced monotonically with time until equilibrium is reached while the opposite happens in the downstream portion. These figures show that the time evolution of the interfacial shape and its properties is similar for both negative and positive Bond numbers, i.e., for both upward-extended and flatter droplets. Figures 3 b and 4 b show that the contact angle increases monotonically with the azimuthal angle from the upstream towards the downstream area of the contact line at all times. Thus, as commonly happens for circular contact lines, the smallest contact angle u at =180 and the largest one d at =0 can be used as a measurement of the interfacial deformation and to provide information on the onset of interfacial sliding e.g., see Ref. 7. Figure 5 a shows the time evolution of the upstream and downstream contact angles, u and d, for an adherent droplet with =1 and 0 =90 in simple shear flow with Ca=0.10, and for several Bond numbers i.e., B d = 0.5,0,0.5,1,1.5,2. This figure clearly reveals that Ca =0.10 is a subcritical flow rate for all these Bond numbers; after the initial transient evolution both angles reach equilibrium where the interfacial shape does not change with time. The contact angles show initially the same evolution for all Bond numbers; however later the droplet with the largest Bond number shows the smaller variation for both upstream and downstream contact angles, u and d, i.e., less deformation. This is also evident in Fig. 5 b where we plot the time evolution of the contact angle difference d u. The droplet profiles at equilibrium for the different Bond numbers are shown in Fig. 5 c. Therefore, as the Bond number increases, the flatter viscous droplet is more stable for a given flow rate. Equivalently, for specific advancing and receding contact angles, A and R, the flatter droplet requires a higher flow rate to start sliding on the solid surface. The explanation for the monotonic decrease in deformation as the Bond number increases becomes clear when one considers the forces acting on the droplet. For viscous drop-
5 Gravitational effects on the deformation of a droplet Phys. Fluids 19, FIG. 3. Deformation of an adherent droplet with =1 and 0 =90 in simple shear flow with Ca=0.10 and for B d = 0.5. a Droplet profile at times t =0,0.5,1,2,10; also included is the undisturbed droplet profile - - for B d =0. b The variation of the contact angle as a function of the azimuthal angle for times t=0,0.2,0.5,1,2,4,10. The azimuthal angle is measured with respect to the positive x-direction as usual. c Droplet shape at t=10, i.e., at equilibrium. FIG. 4. Deformation of an adherent droplet with =1 and 0 =90 in simple shear flow with Ca=0.10 and for B d =2. a Droplet profile at times t =0,0.5,1,2,10; also included is the undisturbed droplet profile - - for B d =0. b The variation of the contact angle as a function of the azimuthal angle for times t=0,0.2,0.5,1,2,4,10. c Droplet shape at t=10, i.e., at equilibrium. lets the dominant deforming force is the shear component F s of the hydrodynamic force on the interface; F s is proportional to the shear stress on the drop d and the drop s surface area S d, i.e., F s d S d. Trying to balance this deforming force is the component of surface tension force in the plane of the wall parallel to the flow direction; this restoring force is proportional to the width or radius w of the contact region and, for small angle difference d u, it scales as F cos u cos d w. As the Bond number increases, the droplet becomes flatter; 7 the increased width w of the contact region results in an increase for the surface tension force F and thus a smaller deformation. Our study also provides information on the portions of the contact line which will slide first. To show this, in Fig. 6 we plot the evolution of the contact angle upstream and downstream change, 0 u and d 0, for a negative and a positive Bond number. From this figure it is evident that for small enough hysteresis, during the transient evolution the drops show an early period where both upstream and downstream sliding are equally favorable. During this regime, if the droplet is placed on the solid forming an initial angle 0 R + A /2, then after some time the downstream por-
6 P. Dimitrakopoulos Phys. Fluids 19, FIG. 6. Time evolution of the contact angle upstream and downstream change, 0 u and d 0, for an adherent droplet with =1 and 0 =90 in simple shear flow with Ca=0.10. a B d = 0.5; b B d =0.5. FIG. 5. Deformation of an adherent droplet with =1 and 0 =90 in simple shear flow with Ca=0.10 for B d = 0.5,0,0.5,1,1.5,2. a Time evolution of the upstream and downstream contact angles, u and d. b Time evolution of the contact angle difference d u. c Droplet profile at time t=10 i.e., well past equilibrium ; also included is the undisturbed droplet profile - - for B d =0. tion of the drop will slide first. The opposite will happen if the initial angle is 0 R + A /2 where the upstream drop portion will begin moving on the solid substrate. For large enough hysteresis, drops with moderate initial angles 0 show a late downstream-favored regime, where the downstream sliding occurs for initial angles 0 0 * with 0 * R + A /2. Thus, with respect to the onset of interfacial sliding, droplets with negative or positive Bond numbers and for moderate initial angles e.g., 0 =90 show similar behavior with those for B d =0 presented in our earlier publication. 7 C. Influence of the initial angle We investigate now the gravitational effects on droplets with smaller and larger initial angles. Our computational results for adherent droplets with 0 =50 reveal that the deformation of droplets with small contact angles is affected differently by the Bond number compared to those with moderate angles. As shown in Fig. 7, instead of a clear monotonic decrease in deformation as the Bond number increases, a droplet with 0 =50 is not affected much by the Bond number. The force balance on these droplets may provide insight on their deformation behavior. Similar to our earlier discussion for 0 =90, as the Bond number increases the flatter droplet is accompanied with an increased width w of the contact region which results in an increase for the surface tension force F. However, for small contact angles, as the Bond number increases, the flatter droplet is accompanied with an increase in the surface area S d and thus an increased shear force F s. Thus these two opposite actions i.e., increase in w and S d counterbalance the effects of Bond number variation. We now turn our attention to viscous droplets with large initial angles. Figure 8 shows that for a droplet with 0 = 120, by increasing the Bond number, the flatter drop shows less deformation, in qualitative agreement with our
7 Gravitational effects on the deformation of a droplet Phys. Fluids 19, FIG. 7. Time evolution of the contact angle difference d u for an adherent droplet with =1 and 0 =50 in simple shear flow with Ca=0.10 and for B d = 0.5,0,0.5. results for 0 =90 presented earlier in Fig. 5. This figure also reveals that the negative Bond number B d = 0.3 is accompanied with a significant deformation especially at the upstream portion of the droplet. These conclusions are also evident from the deformed droplet profiles presented in Fig. 8 c, where it is shown that for droplets with large initial angles, the downstream portion of the interface tends to become tangential to the solid wall. The monotonic decrease in deformation as the Bond number increases can be explained based on the interfacial force balance. Similar to our earlier discussion for small and moderate contact angles, as the Bond number increases the flatter droplet is accompanied with an increased surface tension force F due to the increased width w of the contact region. For droplets with large contact angles, as the Bond number increases the flatter droplet is accompanied with a decrease in the surface area S d and thus a decreased shear force F s. Both actions increase in w and decrease in S d contribute to the monotonic decrease in deformation as the Bond number increases for droplets with large contact angles. They also justify the larger effects of the Bond number variation on the droplet deformation shown in Fig. 8 with respect to those for droplets with moderate angles presented in Fig. 5. With respect to the onset of interfacial sliding, our earlier study for B d =0 Ref. 7 revealed that droplets with large initial angles 0 120, after a rather small early equally favorable period, show a large period where downstream sliding is more favorable than the upstream sliding, i.e., downstream sliding occurs for initial angles 0 0 * with 0 * R + A /2. Our present results depicted in Fig. 9 reveal that the gravitational effects on the onset of interfacial sliding for droplets with large initial angles is more complicated. In particular, this figure demonstrates that a negative Bond number reinforces the behavior found for B d =0, i.e., an upward-extended droplet shows a much larger upstream change in the contact angle 0 u than a downstream change d 0 ; thus for such drops downstream sliding is much more favorable than upstream sliding. However, a positive Bond number for droplets with large initial angles 0 reverts to the FIG. 8. Deformation of an adherent droplet with =1 and 0 =120 in simple shear flow with Ca=0.10, for B d = 0.3,0,0.5,1,1.5. a Time evolution of the upstream and downstream contact angles, u and d. b Time evolution of the contact angle difference d u. c Droplet profile at time t=15 i.e., well past equilibrium for B d =0,0.5,1,1.5, and at time t=6.5 for B d = 0.3; also included is the undisturbed droplet profile - - for B d =0. behavior of B d =0. As our results for B d =1.5 presented in Fig. 9 indicate, the flatter droplet is more prone to slide at its upstream portion. D. Influence of the viscosity ratio In this section, we focus our attention on the effects of the viscosity ratio on the deformation of adherent droplets with positive or negative Bond numbers. As discussed earlier, the Bond number affects viscous droplets with small contact angles in a different way than those for moderate and large angles owing to the different variation of the drop surface S d and thus the shear force F s on the drop with the Bond number.
8 P. Dimitrakopoulos Phys. Fluids 19, FIG. 9. Deformation of an adherent droplet with =1 and 0 =120 in simple shear flow with Ca=0.10. Time evolution of the contact angle downstream and upstream change:, d 0 ; - -, 0 u for B d = 0.3,1.5. On the other hand, inviscid droplets/bubbles i.e., =0 exhibit similar behavior with the Bond number, independent of the contact angle, as illustrated in Fig. 10 where we present the time evolution of the contact angle difference d u for several B d and for 0 =50,90,120. In all cases, the droplet deformation is decreased as the Bond number increases, i.e., the flatter drops are more stable for all contact angles studied. In addition, droplets with large initial angles e.g., 0 =120 show significant deformation especially for small Bond numbers and require more time to reach equilibrium. The explanation for the monotonic decrease in deformation as the Bond number increases for all initial contact angles becomes clear when one considers the forces acting on the droplet. For inviscid droplets, the deforming force is the pressure force F p on the interface which is proportional to the drop s frontal area A d and the pressure change p over the drop F p A d p hw p, 8 where A d scales with the height h and the width w of the droplet. Comparing Eqs. 8 and 7, one readily realizes that as the Bond number increases, the increased width w of the droplet results in an increase of both the deforming force F p and the restoring surface tension force F. Therefore, for inviscid droplets, the monotonic decrease in deformation with the Bond number is mainly caused by the decrease of the droplet s height h which occurs for all contact angles. High-viscosity droplets show similar deformation with that for =1 droplets presented earlier. However, as the viscosity ratio increases from small values, the increased shear force on the droplet results in higher interfacial deformation, e.g., larger d u difference. This is clearly evident in Fig. 11 a where we present the evolution of the contact angle difference d u for a negative and a positive Bond number and for viscosity ratio =0,1,10. Observe that in this figure the time is scaled with 1+ which reflects the fact that the dynamics of droplets with different viscosity ratio evolve on FIG. 10. Deformation of an adherent droplet with =0 in simple shear flow with Ca=0.10. a Time evolution of the contact angle difference d u for 0 =50 and B d = 0.5,0,0.5. b As in a, but for 0 =90 and B d = 0.5,0,0.5,1. c As in a, but for 0 =120 and B d = 0.3,0,1.5. the surface tension time scale which may be described by = 1+ a/ = 1+ Ca f. The increased interfacial deformation for the high-viscosity droplets is also evident at the equilibrium profiles for droplets with negative and positive Bond numbers as seen in Figs. 11 b and 11 c. Itisof interest to note that, for subcritical flow rates, the droplet deformation reaches an asymptotic upper limit for very large viscosity ratios i.e., as discussed for the case B d =0 by Dimitrakopoulos 7 see Fig. 7 of the earlier study.
9 Gravitational effects on the deformation of a droplet Phys. Fluids 19, FIG. 12. Equilibrium deformation vs capillary number Ca for an adherent droplet with =1 and 0 =90. a Upstream and downstream contact angles, u and d, for B d = 0.5,0,0.5,1,1.5,2. b Contact angle downstream and upstream change:, d 0 ; - -, 0 u, for B d = 0.5,1. FIG. 11. Influence of the viscosity ratio on the deformation of an adherent droplet with 0 =90 in simple shear flow with Ca=0.10. Viscosity ratios: =0,1,10. a Time evolution of the contact angle difference d u for Bond number B d : 0.5 and b Droplet profile at equilibrium for B d = 0.5. c As in b, but for B d =0.5. E. Equilibrium deformation of adherent droplets In the previous three subsections we presented the dynamics of the interfacial deformation by considering both the transient evolution and the equilibrium conditions for a representative flow rate i.e., Ca=0.1. In this section, we consider the equilibrium deformation of attached droplets versus Ca for a wide range of low and moderate subcritical flow rates. To derive these equilibrium conditions, we employed our Newton method for equilibrium interfaces under Stokes flow conditions. 6 Our results here verify our earlier conclusions presented in Secs. III B III D while they provide numerical data over a range of low and moderate flow rates. Figure 12 a shows the equilibrium variation of the upstream and downstream contact angles, u and d, with the capillary number Ca for an adherent droplet with =1 and 0 =90 and for different Bond numbers. In agreement with our conclusions in Sec. III B, the equilibrium deformation decreases as the Bond number increases owing to the increased surface tension force which accompanies the flatter interfacial shape. In addition, at equilibrium the downstream portion of the drop shows a higher deformation than its upstream portion for both negative and positive Bond numbers as shown in Fig. 12 b. The influence of the Bond number on the equilibrium variation of the upstream and downstream contact angles for a viscous droplet with 0 =120 is similar to that for 0 =90 as shown in Fig. 13 a. As discussed in Sec. III C, for large initial angles 0, the monotonic decrease of the interfacial deformation with the Bond number results from both the increased surface tension force and the decreased shear force which accompanies the decreased surface area. In addition, Fig. 13 b shows that while negative Bond numbers cause a higher upstream deformation compared to the downstream deformation, this behavior is reverted for positive Bond numbers. The limited influence of the Bond number on the equilibrium variation of the upstream and downstream contact angles with the capillary number Ca for a viscous droplet with 0 =50 is shown in Fig. 14. As discussed in Sec. III C,
10 P. Dimitrakopoulos Phys. Fluids 19, FIG. 14. Variation of the equilibrium upstream and downstream contact angles, u and d, with the capillary number Ca for an adherent droplet with =1 and 0 =50. Bond number B d :, 0.5;, 0; - ---, 0.5. FIG. 13. Equilibrium deformation vs capillary number Ca for an adherent droplet with =1 and 0 =120. a Upstream and downstream contact angles, u and d, for B d = 0.3,0,0.5. c Contact angle downstream and upstream change:, d 0 ; - -, 0 u, for B d = 0.3,0.5. this is a consequence of the fact that, for small initial angles 0, both the shear force and the surface tension force increase with the Bond number. As discussed in Sec. III D, inviscid droplets or bubbles show a similar behavior with the Bond number independently of the initial contact angle 0 as clearly shown in Fig. 15 over a wide range of flow rates. In this case, the monotonic decrease of the drop s height with the Bond number results in a decrease of the deforming pressure force and thus of the interfacial deformation. Figure 16 shows the effects of the viscosity ratio on the droplet deformation for a negative and a positive Bond number. As the viscosity ratio increases from small values, the increased viscous force causes a higher equilibrium deformation, i.e., a higher variation of the upstream and downstream contact angles. F. Influence of the initial shape In the previous sections we have considered attached droplets which, before the initiation of the steady shear flow, are at hydrostatic equilibrium under a specific gravitational influence. However, a droplet may be attached to a solid surface without being at hydrostatic equilibrium under a desired Bond number at the instance of the flow initiation. In such a case, the drop may follow a different transient deformation compared to those presented so far. Figure 17 a shows the time evolution of the upstream and downstream contact angles, u and d, for an adherent droplet with =1 and B d =2 in simple shear flow with Ca =0.10. The initial drop shape is a spherical cap with 0 =90. For the same circular contact line and for this Bond number, the hydrostatic shape is flatter showing a contact angle The evolution of the downstream contact angles d is similar to that presented earlier, i.e., this angle shows a monotonic increase towards its equilibrium value. On the other hand, the upstream contact angle u shows an evolution which has not been observed so far. Instead of a monotonic decrease, this angle first increases towards a maximum of u near t=0.16, and then decreases towards its equilibrium value. The nontrivial contact angle evolution in the upstream portion of the contact line is demonstrated in Fig. 17 b, where we present the variation of the contact angle around the contact line for several times. This contact angle evolution can be understood if we consider that in this case the drop is under the influence of both the shear flow and the gravitational forcing due to the fact that the initial shape differs from the hydrostatic one. The former tends to decrease the upstream contact angle and increase the downstream angle while the latter tends to make the drop shape flatter. Since for B d =2 the hydrostatic shape shows a contact angle i.e., greater than that of the initial shape, the gravitational forcing tends to increase the contact angle all around the contact line. At the downstream portion of the contact line both actions tend to increase the contact angle. On the other hand, at the upstream portion of the contact line the two actions have opposite effects. The maximum in u at the early times reveals that initially the gravitational deformation due to the fact that the initial shape differs from the hydrostatic one dominates the flowinduced deformation, and the drop tends towards its hydrostatic shape; the latter is however never completely reached due to the flow-induced deformation. As shown in Fig. 17 a, after the influence of this gravitational forcing has ceased
11 Gravitational effects on the deformation of a droplet Phys. Fluids 19, FIG. 16. Upstream and downstream contact angles, u and d vs capillary number Ca for an adherent droplet with 0 =90 and for viscosity ratio: =0,1,10. a B d = 0.5; b B d =1. FIG. 15. Upstream and downstream contact angles, u and d vs capillary number Ca for an adherent droplet with =0. a 0 =50 and B d = 0.5,0,0.5. b 0 =90 and B d = 0.5,0,0.5,1,1.5,2. c 0 =120 and B d = 0.3,0,0.5,1,1.5. i.e., around t=1, the droplet follows the same evolution with that when the initial shape is the hydrostatic one. The effects of the gravitational forcing on the contact angle evolution are more pronounced when the drop s initial shape differs considerably from its hydrostatic shape. This is depicted in Fig. 18 where we present the deformation of a droplet whose initial shape is a spherical cap under the large Bond number B d =4. It is of interest to note that in this case the combined action of the shear flow and the increased gravitational forcing causes the downstream portion of the droplet to become rapidly tangential to the solid plane, i.e., the downstream contact angle shows a fast increase towards the limit of 180. We emphasize that, with respect to the onset of interfacial sliding, this type of experiment makes the downstream sliding much more favorable than the upstream sliding. To investigate the case where the initial drop shape is flatter than its hydrostatic configuration for a specific Bond number, in Fig. 19 we present the time evolution of the upstream and downstream contact angles, u and d, for an adherent droplet with =1 and B d = 2 in simple shear flow with Ca=0.10. The initial drop shape is a spherical cap with 0 =90 while for the same circular contact line and for this Bond number the hydrostatic shape is more upward-extended showing a contact angle In this case, the combined action of the shear flow and the gravitational forcing due to the fact that the initial shape differs from the hydrostatic one tend to decrease faster than the contact angle in the upstream portion of the contact line. In the downstream area, the two forces have opposite effects while at early times the gravitational deformation dominates the flow-induced deformation, and d shows a local minimum before its final increase towards its equilibrium value. We emphasize that this gravitational influence is restricted to early times only; afterwards the droplet tends to follow the same evolution with that when the initial shape is the hydro-
12 P. Dimitrakopoulos Phys. Fluids 19, FIG. 18. Deformation of an adherent droplet with =1 and B d =4 in simple shear flow with Ca=0.10; the initial drop shape is a spherical cap with 0 =90. a Time evolution of the upstream and downstream contact angles, u and d. b Droplet profile at times t=0,0.2,1,1.2. FIG. 17. Deformation of an adherent droplet with =1 and B d =2 in simple shear flow with Ca=0.10; the initial drop shape is a spherical cap with 0 =90. a Time evolution of the upstream and downstream contact angles, u and d. Also included are the results - - when the initial drop shape is the hydrostatic shape which shows 0 = b The variation of the contact angle as a function of the azimuthal angle for times t =0,0.08,0.16,0.5,1,2,10. A multiple time-step approach was used with t=0.02 for the time period 0,1 and t=0.1 for 1,10. static one. In addition, this type of experiments makes upstream interfacial sliding much more favorable than downstream sliding. For the same type of experiments, by reducing further the Bond number, the observed behavior is more pronounced as shown in Figs. 20 and 21 where we present the evolution of a droplet whose initial shape is a spherical cap for B d = 2 and B d = 4, respectively. As shown in Fig. 20 the Bond number B d = 2 causes a very fast decrease of the upstream contact angle and thus the upstream portion of the droplet becomes rapidly tangential to the solid surface. On the other hand, at the downstream portion the flow dominates eventually the gravitational forcing and thus after the early minimum, the downstream contact angle increases. By reducing further the Bond number to B d = 4 seen in Fig. 21, the increased negative gravitational forcing causes a very fast upward motion of the entire droplet and thus both the upstream and downstream portion of the droplet tend now to become tangential to the solid surface. IV. DISCUSSION In this paper we have considered the gravitational effects i.e., the influence of positive and negative Bond numbers on the deformation of viscous and inviscid droplets adhering to horizontal rough solid surfaces in steady shear Stokes flows. We have also identified the gravitational effects on the onset of interfacial sliding, i.e., on the portions of the contact line which slide first due to violation of the hysteresis condition. We have considered different initial configurations in- FIG. 19. Time evolution of the upstream and downstream contact angles, u and d, for an adherent droplet with =1 and B d = 1 in simple shear flow with Ca=0.10; the initial drop shape is a spherical cap with 0 =90. A multiple time-step approach was used with t=0.02 for the time period 0,2 and t=0.1 for 2,10.
13 Gravitational effects on the deformation of a droplet Phys. Fluids 19, FIG. 20. Deformation of an adherent droplet with =1 and B d = 2 in simple shear flow with Ca=0.10; the initial drop shape is a spherical cap with 0 =90. a Time evolution of the upstream and downstream contact angles, u and d. b Droplet profile at times t=0,0.4,0.8,1.2,1.6,2. FIG. 21. Deformation of an adherent droplet with =1 and B d = 4 in simple shear flow with Ca=0.10; the initial drop shape is a spherical cap with 0 =90. a Time evolution of the upstream and downstream contact angles, u and d. b Droplet profile at times t=0,0.2,0.4,0.5. cluding drops which are at hydrostatic equilibrium for the desired Bond number at the instance of the flow initiation as well as fluids which are not at hydrostatic equilibrium. When the interfacial system is at hydrostatic equilibrium at the initiation of a steady shear flow, our study shows that the Bond number affects the deformation of viscous droplets with moderate and large initial contact angles in a different way than those for small angles. In particular, the former shows a monotonic decrease in deformation as the Bond number increases as a result of the associated increase of the drop s width and thus the restoring surface tension force for the flatter, high Bond number, droplets. This behavior is more pronounced for droplets with large initial contact angles which are also accompanied with a decrease in their surface area and thus in the deforming shear force. On the other hand for small contact angles, as the Bond number increases, the flatter viscous droplet is accompanied with an increased shear force owing to the increase of the drop s surface area as well as with an increase in the surface tension force. These opposite actions counterbalance the effects of Bond number increase and thus the deformation of droplets with small contact angles is not affected much by the Bond number. Inviscid droplets with different initial contact angles show similar behavior as the Bond number increases, i.e., their deformation is monotonically decreased. In this case, the increase of the drop width for the flatter, high Bond number droplets results in an increase in the restoring surface tension force but also in the deforming force which for inviscid droplets/bubbles is the pressure force. Therefore, for inviscid droplets, the monotonic decrease in deformation with the Bond number is mainly caused by the decrease of the droplet s height which occurs for all contact angles. We emphasize that our forces analysis, which helps to explain our numerical results of the gravitational influence on the drop deformation, relies mostly on the effects of the Bond number on the drop geometric characteristics, i.e., its width, height, and surface area. However, the deforming shear and pressure forces are also affected by physical variables, i.e., the shear stress on the drop d and the pressure change p over the drop, respectively. These physical variables are implicit functions of the drop shape and are not expected to change our aforementioned analysis. This is supported by the scaling analysis for small angles which predicts that d G and p 0 G, i.e., they are not affected by the Bond number. With respect to the onset of interfacial sliding, our study shows that droplets under positive or negative Bond numbers show similar behavior to that for B d =0 identified in our earlier work. 7 In particular, droplets with small and moderate initial angles show an early period where both upstream and downstream sliding are equally favorable as well as a late downstream-favored period; by contrast, droplets with large
14 P. Dimitrakopoulos Phys. Fluids 19, initial angles, after a rather small early equally favorable period, show a large period where downstream sliding is more favorable than the upstream sliding. In addition, for droplets with large initial angles our present study shows that a negative Bond number reinforces the behavior found for B d =0, i.e., upward-extended droplets exhibit a more favorable downstream initial sliding. For the same case, positive Bond numbers weaken this behavior while large enough Bond numbers revert to the behavior of B d =0, i.e., they make upstream sliding more favorable. When the interfacial system is not at hydrostatic equilibrium at the initiation of a steady shear flow, its dynamic evolution is more complicated. In this case the droplet deformation is affected by both the shear flow and the gravitational forcing due to the fact that the initial shape differs from the hydrostatic one. When the drop is initially more upward-extended than its hydrostatic shape, both forces contribute to a fast increase of the downstream contact angle; however at the upstream portion of the contact line they have opposite effects. As a result at early times the gravitational deformation dominates the flow-induced deformation and initially the upstream contact angle increases towards a local maximum before its final flow-induced decrease. For large deviations from the hydrostatic shape, the increased gravitational forcing causes the downstream portion of the droplet to rapidly become tangential to the solid surface. Thus, this case makes the downstream sliding much more favorable than the upstream sliding. The opposite happens for droplets which are initially flatter than their hydrostatic configurations. Therefore this type of experiments can be employed to initiate interfacial sliding at a desired portion of the contact line. As discussed in Sec. I above, earlier studies have commonly ignored the gravitational effects on the drop deformation. One exception is the work by Dimitrakopoulos and Higdon 5 who considered the equilibrium conditions just before the final flow-induced drop displacement from rigid boundaries for zero and positive Bond numbers. Via computational work and asymptotic analysis based on small-angle lubrication theory, Dimitrakopoulos and Higdon 5 demonstrated that for a given flow rate, the deformation of viscous two-dimensional droplets is increased with the Bond number due to the corresponding increase of the drop s length. They also showed that the opposite happens for two-dimensional inviscid droplets due to the decrease of the drop height. Therefore, the present study has shown that the gravitational influence on a three-dimensional droplet differs considerably for that on a two-dimensional drop i.e., a cylindrical cap owing to the ability of the three-dimensional droplets to expand or contract in the lateral i.e., cross-flow direction changing the interfacial dynamics. In essence, our present work identifies the three-dimensional nature of the deformation and sliding of adherent droplets which has also been demonstrated in the earlier work of Dimitrakopoulos and Higdon on the yield criteria for displacement of threedimensional droplets, 6,16 18 and the recent experiments of Podgorski et al. 23 ACKNOWLEDGMENTS This work was supported in part by the National Science Foundation. Acknowledgment is made to the Donors of the American Chemical Society Petroleum Research Fund for partial support of this research. Some computations were performed on multiprocessor supercomputers provided by the National Center for Supercomputing Applications NCSA in Illinois. 1 J. Bear, Dynamics of Fluids in Porous Media Dover, New York, H. A. Stone, A. D. Stroock, and A. Ajdari, Engineering flows in small devices: Microfluidics toward a lab-on-a-chip, Annu. Rev. Fluid Mech. 36, V. Cristini and Y.-C. Tan, Theory and numerical simulation of droplet dynamics in complex flows: A review, Lab Chip 4, R. O. Hynes, Integrins, versatility, modulation and signaling in cell adhesion, Cell 69, P. Dimitrakopoulos and J. J. L. Higdon, Displacement of fluid droplets from solid surfaces in low-reynolds-number shear flows, J. Fluid Mech. 336, P. Dimitrakopoulos and J. J. L. Higdon, On the displacement of threedimensional fluid droplets from solid surfaces in low-reynolds-number shear flows, J. Fluid Mech. 377, P. Dimitrakopoulos, Deformation of a droplet adhering to a solid surface in shear flow: Onset of interfacial sliding, J. Fluid Mech. 580, J. Q. Feng and O. A. Basaran, Shear flow over a translationally symmetric cylindrical bubble pinned on a slot in a plane wall, J. Fluid Mech. 275, X. Li and C. Pozrikidis, Shear flow over a liquid drop adhering to a solid surface, J. Fluid Mech. 307, A. D. Schleizer and R. T. Bonnecaze, Displacement of a two-dimensional immiscible droplet adhering to a wall in shear- and pressure-driven flows, J. Fluid Mech. 383, S. Yon and C. Pozrikidis, Deformation of a liquid drop adhering to a plane wall: Significance of the drop viscosity and the effect of an insoluble surfactant, Phys. Fluids 11, P. D. M. Spelt, Shear flow past two-dimensional droplets pinned or moving on an adhering channel wall at moderate Reynolds numbers: A numerical study, J. Fluid Mech. 561, E. B. Dussan V. and R. T.-P. Chow, On the ability of drops or bubbles to stick to non-horizontal surfaces of solids, J. Fluid Mech. 137, E. B. Dussan V. On the ability of drops or bubbles to stick to nonhorizontal surfaces of solids. Part 2. Small drops or bubbles having contact angles of arbitrary size, J. Fluid Mech. 151, E. B. Dussan V. On the ability of drops to stick to surfaces of solids. Part 3. The influences of the motion of the surrounding fluid on dislodging drops, J. Fluid Mech. 174, P. Dimitrakopoulos and J. J. L. Higdon, On the gravitational displacement of three-dimensional fluid droplets from inclined solid surfaces, J. Fluid Mech. 395, P. Dimitrakopoulos and J. J. L. Higdon, On the displacement of threedimensional fluid droplets adhering to a plane wall in viscous pressuredriven flows, J. Fluid Mech. 435, P. Dimitrakopoulos and J. J. L. Higdon, On the displacement of threedimensional fluid bridges from solid surfaces in viscous pressure-driven flows, Phys. Fluids 15, Q. J. Kang, D. X. Zhang, and S. Y. Chen, Displacement of a threedimensional immiscible droplet in a duct, J. Fluid Mech. 545, J. Zhang, M. J. Miksis, and S. G. Bankoff, Nonlinear dynamics of a two-dimensional viscous drop under shear flow, Phys. Fluids 18, P. Dimitrakopoulos, Interfacial dynamics in Stokes flow via a threedimensional fully-implicit interfacial spectral boundary element algorithm, J. Comput. Phys. 225, G. P. Muldowney and J. J. L. Higdon, A spectral boundary element approach to three-dimensional Stokes flow, J. Fluid Mech. 298, T. Podgorski, J.-M. Flesselles, and L. Limat, Corners, cusps, and pearls in running drops, Phys. Rev. Lett. 87,
Deformation of a droplet adhering to a solid surface in shear flow: onset of interfacial sliding
J. Fluid Mech. (27), vol. 58, pp. 451 466. c 27 Cambridge University Press doi:1.117/s2211275721 Printed in the United Kingdom 451 Deformation of a droplet adhering to a solid surface in shear flow: onset
More informationOn the displacement of three-dimensional fluid droplets adhering to a plane wall in viscous pressure-driven flows
J. Fluid Mech. (2001), vol. 435, pp. 327 350. Printed in the United Kingdom c 2001 Cambridge University Press 327 On the displacement of three-dimensional fluid droplets adhering to a plane wall in viscous
More informationNature of the hemodynamic forces exerted on vascular endothelial cells or leukocytes adhering to the surface of blood vessels
PHYSICS OF FLUIDS 18, 087107 2006 Nature of the hemodynamic forces exerted on vascular endothelial cells or leukocytes adhering to the surface of blood vessels Yechun Wang and P. Dimitrakopoulos a Department
More informationFormation of Droplets and Bubbles in a Microfluidic T-junction. Scaling and Mechanism of Break-Up. Supplementary Information
Formation of Droplets and Bubbles in a Microfluidic T-junction Scaling and Mechanism of Break-Up Supplementary Information Piotr Garstecki 1,2 *, Michael J. Fuerstman 1, Howard A. Stone 3 and George M.
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 informationShell Balances in Fluid Mechanics
Shell Balances in Fluid Mechanics R. Shankar Subramanian Department of Chemical and Biomolecular Engineering Clarkson University When fluid flow occurs in a single direction everywhere in a system, shell
More informationFluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition
Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition C. Pozrikidis m Springer Contents Preface v 1 Introduction to Kinematics 1 1.1 Fluids and solids 1 1.2 Fluid parcels and flow
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 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 informationMicrofluidics 2 Surface tension, contact angle, capillary flow
MT-0.6081 Microfluidics and BioMEMS Microfluidics 2 Surface tension, contact angle, capillary flow 28.1.2017 Ville Jokinen Surface tension & Surface energy Work required to create new surface = surface
More informationFluid Mechanics Prof. S.K. Som Department of Mechanical Engineering Indian Institute of Technology, Kharagpur
Fluid Mechanics Prof. S.K. Som Department of Mechanical Engineering Indian Institute of Technology, Kharagpur Lecture - 42 Flows with a Free Surface Part II Good morning. I welcome you to this session
More information10.52 Mechanics of Fluids Spring 2006 Problem Set 3
10.52 Mechanics of Fluids Spring 2006 Problem Set 3 Problem 1 Mass transfer studies involving the transport of a solute from a gas to a liquid often involve the use of a laminar jet of liquid. The situation
More informationFluid Mechanics. du dy
FLUID MECHANICS Technical English - I 1 th week Fluid Mechanics FLUID STATICS FLUID DYNAMICS Fluid Statics or Hydrostatics is the study of fluids at rest. The main equation required for this is Newton's
More informationEnhancement of Heat Transfer by an Electric Field for a Drop Translating at Intermediate Reynolds Number
Rajkumar Subramanian M. A. Jog 1 e-mail: milind.jog@uc.edu Department of Mechanical, Industrial, and Nuclear Engineering, University of Cincinnati, Cincinnati, OH 45221-0072 Enhancement of Heat Transfer
More informationLinear stability of radial displacements in porous media: Influence of velocity-induced dispersion and concentration-dependent diffusion
PHYSICS OF FLUIDS VOLUME 16, NUMBER 10 OCTOBER 2004 Linear stability of radial displacements in porous media: Influence of velocity-induced dispersion and concentration-dependent diffusion A. Riaz Department
More informationThe Shape of a Rain Drop as determined from the Navier-Stokes equation John Caleb Speirs Classical Mechanics PHGN 505 December 12th, 2011
The Shape of a Rain Drop as determined from the Navier-Stokes equation John Caleb Speirs Classical Mechanics PHGN 505 December 12th, 2011 Derivation of Navier-Stokes Equation 1 The total stress tensor
More informationVibration of submillimeter-size supported droplets
PHYSICAL REVIEW E 73, 041602 2006 Vibration of submillimeter-size supported droplets Franck Celestini* and Richard Kofman Laboratoire de Physique de la Matière Condensée, UMR 6622, CNRS, Université de
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 informationREE Internal Fluid Flow Sheet 2 - Solution Fundamentals of Fluid Mechanics
REE 307 - Internal Fluid Flow Sheet 2 - Solution Fundamentals of Fluid Mechanics 1. Is the following flows physically possible, that is, satisfy the continuity equation? Substitute the expressions for
More informationUNIT IV BOUNDARY LAYER AND FLOW THROUGH PIPES Definition of boundary layer Thickness and classification Displacement and momentum thickness Development of laminar and turbulent flows in circular pipes
More informationColloidal Particles at Liquid Interfaces: An Introduction
1 Colloidal Particles at Liquid Interfaces: An Introduction Bernard P. Binks and Tommy S. Horozov Surfactant and Colloid Group, Department of Chemistry, University of Hull, Hull, HU6 7RX, UK 1.1 Some Basic
More informationDifferential relations for fluid flow
Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow
More informationPrinciples of Convection
Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid
More informationAP Physics Laboratory #6.1: Analyzing Terminal Velocity Using an Interesting Version of Atwood s Machine
AP Physics Laboratory #6.1: Analyzing Terminal Velocity Using an Interesting Version of Atwood s Machine Name: Date: Lab Partners: PURPOSE The purpose of this Laboratory is to study a system as it approaches
More informationExperimental and Theoretical Study of Motion of Drops on Horizontal Solid Surfaces with a Wettability Gradient Nadjoua Moumen
Experimental and Theoretical Study of Motion of Drops on Horizontal Solid Surfaces with a Wettability Gradient Nadjoua Moumen Department of Chemical and Biomolecular Engineering Clarkson University Outline
More informationch-01.qxd 8/4/04 2:33 PM Page 1 Part 1 Basic Principles of Open Channel Flows
ch-01.qxd 8/4/04 2:33 PM Page 1 Part 1 Basic Principles of Open Channel Flows ch-01.qxd 8/4/04 2:33 PM Page 3 Introduction 1 Summary The introduction chapter reviews briefly the basic fluid properties
More informationThermocapillary Migration of a Drop
Thermocapillary Migration of a Drop An Exact Solution with Newtonian Interfacial Rheology and Stretching/Shrinkage of Interfacial Area Elements for Small Marangoni Numbers R. BALASUBRAMANIAM a AND R. SHANKAR
More informationMECHANICAL PROPERTIES OF FLUIDS:
Important Definitions: MECHANICAL PROPERTIES OF FLUIDS: Fluid: A substance that can flow is called Fluid Both liquids and gases are fluids Pressure: The normal force acting per unit area of a surface is
More informationFigure 11.1: A fluid jet extruded where we define the dimensionless groups
11. Fluid Jets 11.1 The shape of a falling fluid jet Consider a circular orifice of a radius a ejecting a flux Q of fluid density ρ and kinematic viscosity ν (see Fig. 11.1). The resulting jet accelerates
More informationPhysics 3 Summer 1990 Lab 7 - Hydrodynamics
Physics 3 Summer 1990 Lab 7 - Hydrodynamics Theory Consider an ideal liquid, one which is incompressible and which has no internal friction, flowing through pipe of varying cross section as shown in figure
More informationFriction Factors and Drag Coefficients
Levicky 1 Friction Factors and Drag Coefficients Several equations that we have seen have included terms to represent dissipation of energy due to the viscous nature of fluid flow. For example, in the
More informations and FE X. A. Flow measurement B. properties C. statics D. impulse, and momentum equations E. Pipe and other internal flow 7% of FE Morning Session I
Fundamentals of Engineering (FE) Exam General Section Steven Burian Civil & Environmental Engineering October 26, 2010 s and FE X. A. Flow measurement B. properties C. statics D. impulse, and momentum
More informationExercise: concepts from chapter 10
Reading:, Ch 10 1) The flow of magma with a viscosity as great as 10 10 Pa s, let alone that of rock with a viscosity of 10 20 Pa s, is difficult to comprehend because our common eperience is with s like
More informationUpper and Lower Bounds for Interfacial Tension using Spinning Drop Devices. D. D. Joseph, M. Arney and G. Ma 1
Upper and Lower Bounds for Interfacial Tension using Spinning Drop Devices D. D. Joseph, M. Arney and G. Ma 1 In this note we show how to use spinning drop devices to determine lower and upper bounds for
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 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 informationInterfacial Flows of Contact Line Dynamics and Liquid Displacement in a Circular Microchannel
Proceedings of the 3 rd World Congress on Mechanical, Chemical, and Material Engineering (MCM'17) Rome, Italy June 8 10, 2017 Paper No. HTFF 159 ISSN: 2369-8136 DOI: 10.11159/htff17.159 Interfacial Flows
More informationIntroduction to Marine Hydrodynamics
1896 1920 1987 2006 Introduction to Marine Hydrodynamics (NA235) Department of Naval Architecture and Ocean Engineering School of Naval Architecture, Ocean & Civil Engineering First Assignment The first
More informationHYDRAULICS STAFF SELECTION COMMISSION CIVIL ENGINEERING STUDY MATERIAL HYDRAULICS
1 STAFF SELECTION COMMISSION CIVIL ENGINEERING STUDY MATERIAL Syllabus Hydraulics ( Fluid Mechanics ) Fluid properties, hydrostatics, measurements of flow, Bernoulli's theorem and its application, flow
More informationFE Fluids Review March 23, 2012 Steve Burian (Civil & Environmental Engineering)
Topic: Fluid Properties 1. If 6 m 3 of oil weighs 47 kn, calculate its specific weight, density, and specific gravity. 2. 10.0 L of an incompressible liquid exert a force of 20 N at the earth s surface.
More informationP = 1 3 (σ xx + σ yy + σ zz ) = F A. It is created by the bombardment of the surface by molecules of fluid.
CEE 3310 Thermodynamic Properties, Aug. 27, 2010 11 1.4 Review A fluid is a substance that can not support a shear stress. Liquids differ from gasses in that liquids that do not completely fill a container
More informationspreading of drops on soft surfaces
Supplementary Material on Electrically modulated dynamic spreading of drops on soft surfaces Ranabir Dey 1, Ashish Daga 1, Sunando DasGupta 2,3, Suman Chakraborty 1,3 1 Department of Mechanical Engineering,
More informationSlow viscous flow in a microchannel with similar and different superhydrophobic walls
Journal of Physics: Conference Series PAPER OPEN ACCESS Slow viscous flow in a microchannel with similar and different superhydrophobic walls To cite this article: A I Ageev and A N Osiptsov 2018 J. Phys.:
More informationPaper No. : 04 Paper Title: Unit Operations in Food Processing Module- 18: Circulation of fluids through porous bed
Paper No. : 04 Paper Title: Unit Operations in Food Processing Module- 18: Circulation of fluids through porous bed 18.1 Introduction A typical packed bed is a cylindrical column that is filled with a
More informationUNIT II CONVECTION HEAT TRANSFER
UNIT II CONVECTION HEAT TRANSFER Convection is the mode of heat transfer between a surface and a fluid moving over it. The energy transfer in convection is predominately due to the bulk motion of the fluid
More informationCapillary Phenomena in Assemblies of Parallel Cylinders
Capillary Phenomena in Assemblies of Parallel Cylinders I. Capillary Rise between Two Cylinders H. M. PRINCEN Lever Brothers Company, Research and Development Division, Edgewater, New Jersey 07020 Received
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 informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationDepartment of Mechanical Engineering
Department of Mechanical Engineering AMEE401 / AUTO400 Aerodynamics Instructor: Marios M. Fyrillas Email: eng.fm@fit.ac.cy HOMEWORK ASSIGNMENT #2 QUESTION 1 Clearly there are two mechanisms responsible
More informationINTRODUCTION TO FLUID MECHANICS June 27, 2013
INTRODUCTION TO FLUID MECHANICS June 27, 2013 PROBLEM 3 (1 hour) A perfect liquid of constant density ρ and constant viscosity µ fills the space between two infinite parallel walls separated by a distance
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 informationLiquids and solids are essentially incompressible substances and the variation of their density with pressure is usually negligible.
Properties of Fluids Intensive properties are those that are independent of the mass of a system i.e. temperature, pressure and density. Extensive properties are those whose values depend on the size of
More informationAdvanced Hydraulics Prof. Dr. Suresh A. Kartha Department of Civil Engineering Indian Institute of Technology, Guwahati
Advanced Hydraulics Prof. Dr. Suresh A. Kartha Department of Civil Engineering Indian Institute of Technology, Guwahati Module - 2 Uniform Flow Lecture - 1 Introduction to Uniform Flow Good morning everyone,
More informationLaboratory 9: The Viscosity of Liquids
Laboratory 9: The Viscosity of Liquids Introduction The essential difference between solids and fluids lies in the nature of their response to the socalled shearing stress. In solids, an elastic force
More informationR09. d water surface. Prove that the depth of pressure is equal to p +.
Code No:A109210105 R09 SET-1 B.Tech II Year - I Semester Examinations, December 2011 FLUID MECHANICS (CIVIL ENGINEERING) Time: 3 hours Max. Marks: 75 Answer any five questions All questions carry equal
More informationNUMERICAL SIMULATION OF MICRO-FILTRATION OF OIL-IN-WATER EMULSIONS. Tohid Darvishzadeh
NUMERICAL SIMULATION OF MICRO-FILTRATION OF OIL-IN-WATER EMULSIONS By Tohid Darvishzadeh A DISSERTATION Submitted to Michigan State University in partial fulfillment of the requirements for the degree
More informationModule 9: Packed beds Lecture 29: Drag, particles settling. Flow through a packed bed of solids. Drag. Criteria of settling.
Flow through a packed bed of solids Drag Criteria of settling Hindred settling file:///d /Web%20Course/Dr.%20Nishith%20Verma/local%20server/fluid_mechanics/lecture29/29_1.htm[5/9/2012 3:38:37 PM] Flow
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 informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationAn-Najah National University Civil Engineering Department. Fluid Mechanics. Chapter 1. General Introduction
1 An-Najah National University Civil Engineering Department Fluid Mechanics Chapter 1 General Introduction 2 What is Fluid Mechanics? Mechanics deals with the behavior of both stationary and moving bodies
More informationLECTURE 4: Marangoni flows
LECTURE 4: Marangoni flows Marangoni flows are those driven by surface tension gradients. In general, surface tension σ depends on both the temperature and chemical composition at the interface; consequently,
More informationUniversity of Hail Faculty of Engineering DEPARTMENT OF MECHANICAL ENGINEERING. ME Fluid Mechanics Lecture notes. Chapter 1
University of Hail Faculty of Engineering DEPARTMENT OF MECHANICAL ENGINEERING ME 311 - Fluid Mechanics Lecture notes Chapter 1 Introduction and fluid properties Prepared by : Dr. N. Ait Messaoudene Based
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 informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationStudy on Non-Uniqueness of Taylor Vortex Flow Changing Inner Cylinder Acceleration Time
World Journal of Mechanics, 2018, 8, 301-310 http://www.scirp.org/journal/wjm ISSN Online: 2160-0503 ISSN Print: 2160-049X Study on Non-Uniqueness of Taylor Vortex Flow Changing Inner Cylinder Acceleration
More informationTransport by convection. Coupling convection-diffusion
Transport by convection. Coupling convection-diffusion 24 mars 2017 1 When can we neglect diffusion? When the Peclet number is not very small we cannot ignore the convection term in the transport equation.
More informationLecture 2: Hydrodynamics at milli micrometer scale
1 at milli micrometer scale Introduction Flows at milli and micro meter scales are found in various fields, used for several processes and open up possibilities for new applications: Injection Engineering
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 informationFluid Dynamics for Ocean and Environmental Engineering Homework #2 Viscous Flow
OCEN 678-600 Fluid Dynamics for Ocean and Environmental Engineering Homework #2 Viscous Flow Date distributed : 9.18.2005 Date due : 9.29.2005 at 5:00 pm Return your solution either in class or in my mail
More informationAn Adsorption Desorption-Controlled Surfactant on a Deforming Droplet
JOURNAL OF COLLOID AND INTERFACE SCIENCE 208, 68 80 (1998) ARTICLE NO. CS985816 An Adsorption Desorption-Controlled Surfactant on a Deforming Droplet Charles D. Eggleton* and Kathleen J. Stebe,1 *Department
More informationInstability of a rotating thread in a second immiscible liquid
PHYSICS OF FLUIDS VOLUME 16, NUMBER 1 JANUARY 2004 Instability of a rotating thread in a second immiscible liquid J. Ashmore a) and H. A. Stone Division of Engineering and Applied Sciences, Harvard University,
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 informationIn the name of Allah the most beneficent the most merciful
In the name of Allah the most beneficent the most merciful Transient flows of Maxwell fluid with slip conditions Authors Dr. T. Hayat & Sahrish Zaib Introduction Non-Newtonian Newtonian fluid Newtonian
More informationTwo-Fluid Model 41. Simple isothermal two-fluid two-phase models for stratified flow:
Two-Fluid Model 41 If I have seen further it is by standing on the shoulders of giants. Isaac Newton, 1675 3 Two-Fluid Model Simple isothermal two-fluid two-phase models for stratified flow: Mass and momentum
More informationConvective Mass Transfer
Convective Mass Transfer Definition of convective mass transfer: The transport of material between a boundary surface and a moving fluid or between two immiscible moving fluids separated by a mobile interface
More information150A Review Session 2/13/2014 Fluid Statics. Pressure acts in all directions, normal to the surrounding surfaces
Fluid Statics Pressure acts in all directions, normal to the surrounding surfaces or Whenever a pressure difference is the driving force, use gauge pressure o Bernoulli equation o Momentum balance with
More informationUNIT II Real fluids. FMM / KRG / MECH / NPRCET Page 78. Laminar and turbulent flow
UNIT II Real fluids The flow of real fluids exhibits viscous effect that is they tend to "stick" to solid surfaces and have stresses within their body. You might remember from earlier in the course Newtons
More informationConvection. forced convection when the flow is caused by external means, such as by a fan, a pump, or atmospheric winds.
Convection The convection heat transfer mode is comprised of two mechanisms. In addition to energy transfer due to random molecular motion (diffusion), energy is also transferred by the bulk, or macroscopic,
More informationLecture 6 Friction. Friction Phenomena Types of Friction
Lecture 6 Friction Tangential forces generated between contacting surfaces are called friction forces and occur to some degree in the interaction between all real surfaces. whenever a tendency exists for
More informationCapillary-gravity waves: The effect of viscosity on the wave resistance
arxiv:cond-mat/9909148v1 [cond-mat.soft] 10 Sep 1999 Capillary-gravity waves: The effect of viscosity on the wave resistance D. Richard, E. Raphaël Collège de France Physique de la Matière Condensée URA
More informationChapter 6: Momentum Analysis
6-1 Introduction 6-2Newton s Law and Conservation of Momentum 6-3 Choosing a Control Volume 6-4 Forces Acting on a Control Volume 6-5Linear Momentum Equation 6-6 Angular Momentum 6-7 The Second Law of
More informationMECHANICAL PROPERTIES OF FLUIDS
CHAPTER-10 MECHANICAL PROPERTIES OF FLUIDS QUESTIONS 1 marks questions 1. What are fluids? 2. How are fluids different from solids? 3. Define thrust of a liquid. 4. Define liquid pressure. 5. Is pressure
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 informationANALYSIS OF FLOW IN A CONCENTRIC ANNULUS USING FINITE ELEMENT METHOD
Nigerian Journal of Technology (NIJOTECH) Vol 35, No 2, April 2016, pp 344 348 Copyright Faculty of Engineering, University of Nigeria, Nsukka, Print ISSN: 0331-8443, Electronic ISSN: 2467-8821 wwwnijotechcom
More informationNumerical Simulation of the Hagemann Entrainment Experiments
CCC Annual Report UIUC, August 14, 2013 Numerical Simulation of the Hagemann Entrainment Experiments Kenneth Swartz (BSME Student) Lance C. Hibbeler (Ph.D. Student) Department of Mechanical Science & Engineering
More information(Refer Slide Time 1:25)
Mechanical Measurements and Metrology Prof. S. P. Venkateshan Department of Mechanical Engineering Indian Institute of Technology, Madras Module - 2 Lecture - 24 Transient Response of Pressure Transducers
More informationAny correspondence concerning this service should be sent to The Strathprints Administrator:
Sullivan, J.M. and Wilson, S.K. and Duffy, B.R. 2008 Air-blown rivulet flow of a perfectly wetting fluid on an inclined substrate. In: Progress in Industrial Mathematics at ECMI 2006. Springer, pp. 774-778.
More informationInterfacial Instabilities in a Microfluidic Hele-Shaw Cell: Supplemental
Supplementary Material (ESI) for Soft Matter This journal is The Royal Society of Chemistry 2008 Interfacial Instabilities in a Microfluidic Hele-Shaw Cell: Supplemental Michinao Hashimoto 1, Piotr Garstecki
More informationPORE-SCALE PHASE FIELD MODEL OF TWO-PHASE FLOW IN POROUS MEDIUM
Excerpt from the Proceedings of the COMSOL Conference 2010 Paris PORE-SCALE PHASE FIELD MODEL OF TWO-PHASE FLOW IN POROUS MEDIUM Igor Bogdanov 1*, Sylvain Jardel 1, Anis Turki 1, Arjan Kamp 1 1 Open &
More informationFigure 3: Problem 7. (a) 0.9 m (b) 1.8 m (c) 2.7 m (d) 3.6 m
1. For the manometer shown in figure 1, if the absolute pressure at point A is 1.013 10 5 Pa, the absolute pressure at point B is (ρ water =10 3 kg/m 3, ρ Hg =13.56 10 3 kg/m 3, ρ oil = 800kg/m 3 ): (a)
More informationMicrofluidics 1 Basics, Laminar flow, shear and flow profiles
MT-0.6081 Microfluidics and BioMEMS Microfluidics 1 Basics, Laminar flow, shear and flow profiles 11.1.2017 Ville Jokinen Outline of the next 3 weeks: Today: Microfluidics 1: Laminar flow, flow profiles,
More informationMiddle East Technical University Department of Mechanical Engineering ME 305 Fluid Mechanics I Fall 2018 Section 4 (Dr.
Reading Assignments Middle East Technical University Department of Mechanical Engineering ME 305 Fluid Mechanics I Fall 2018 Section 4 (Dr. Sert) Study Set 1 You can find the answers of some of the following
More informationChapter 9 Fluids. Pressure
Chapter 9 Fluids States of Matter - Solid, liquid, gas. Fluids (liquids and gases) do not hold their shapes. In many cases we can think of liquids as being incompressible. Liquids do not change their volume
More informationThin Film Breakup and Rivulet Evolution Modeling
Thin Film Breakup and Rivulet Evolution Modeling Nicola Suzzi, Giulio Croce 2, Paola D Agaro 3 DIEG - Dipartimento di Ingegneria Elettrica, Gestionale e Meccanica - Universitá di Udine - Via delle Scienze
More informationINTRODUCTION DEFINITION OF FLUID. U p F FLUID IS A SUBSTANCE THAT CAN NOT SUPPORT SHEAR FORCES OF ANY MAGNITUDE WITHOUT CONTINUOUS DEFORMATION
INTRODUCTION DEFINITION OF FLUID plate solid F at t = 0 t > 0 = F/A plate U p F fluid t 0 t 1 t 2 t 3 FLUID IS A SUBSTANCE THAT CAN NOT SUPPORT SHEAR FORCES OF ANY MAGNITUDE WITHOUT CONTINUOUS DEFORMATION
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 informationCorrection of Lamb s dissipation calculation for the effects of viscosity on capillary-gravity waves
PHYSICS OF FLUIDS 19, 082105 2007 Correction of Lamb s dissipation calculation for the effects of viscosity on capillary-gravity waves J. C. Padrino and D. D. Joseph Aerospace Engineering and Mechanics
More informationTime dependent singularities in incompressible potential flow
Time dependent singularities in incompressible potential flow Daniel Weihs Department of Aerospace Engineering and Autonomous Systems Program Technion, Haifa 32000, Israel Abstract The flow-field resulting
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 informationME3560 Tentative Schedule Spring 2019
ME3560 Tentative Schedule Spring 2019 Week Number Date Lecture Topics Covered Prior to Lecture Read Section Assignment Prep Problems for Prep Probs. Must be Solved by 1 Monday 1/7/2019 1 Introduction to
More information