Purely irrotational theories of the effects of viscosity and viscoelasticity on capillary instability of a liquid cylinder
|
|
- Myrtle Wood
- 5 years ago
- Views:
Transcription
1 . Non-Newtonian Fluid Mech. 19 ( Purely irrotational theories of the effects of viscosity and viscoelasticity on capillary instability of a liquid cylinder. Wang a, D.D. oseph a,, T. Funada b a Department of Aerospace Engineering and Mechanics, University of Minnesota, 110 Union St. SE, Minneapolis, MN 55455, USA b Department of Digital Engineering, Numazu College of Technology, 3600 Ooka, Numazu, Shizuoka , apan Received 3 September 004; received in revised form 19 March 005; accepted 0 April 005 Abstract Capillary instability of a liquid cylinder can arise when either the interior or exterior fluid is a gas of negligible density and viscosity. The shear stress must vanish at the gas-liquid interface but it does not vanish in irrotational flows. oseph and Wang [D.D. oseph,. Wang, The dissipation approximation and viscous potential flow,. Fluid Mech. 505 ( ] derived an additional viscous correction to the irrotational pressure. They argued that this pressure arises in a boundary layer induced by the unphysical discontinuity of the shear stress. Wang et al. [. Wang, D.D. oseph, T. Funada, Pressure correction for potential flow analysis of capillary instability of viscous fluids,. Fluid Mech. 5 ( ] showed that the dispersion relation for capillary instability in the Newtonian case is almost indistinguishable from the exact solution when the additional pressure contribution is included in the irrotational theory. Here we extend the formulation for the additional pressure to potential flows of viscoelastic fluids in flows governed by linearized equations, and apply this additional pressure to capillary instability of viscoelastic liquid filaments of effreys type. The shear stress at the gas liquid interface cannot be made to vanish in an irrotational theory, but the explicit effect of this uncompensated shear stress can be removed from the global equation for the evolution of the energy of disturbances. This line of thought allows us to present the additional pressure theory without appeal to boundary layers. The validity of this purely irrotational theory can be judged by comparison with the exact solutions of Navier Stokes equations. Here we show that our purely irrotational theory is in remarkably good agreement with the exact solution in linear analysis of the capillary instability of a viscoelastic liquid cylinder. 005 Elsevier B.. All rights reserved. Keywords: Capillary instability; iscoelastic potential flow; Additional pressure contribution; Dissipation method 1. Introduction Capillary instability of a liquid cylinder of mean radius R leading to capillary collapse can be described as a neckdown due to surface tension γ in which fluid is ejected from the throat of the neck, leading to a smaller neck and greater neckdown capillary force as seen in the diagram in Fig. 1. Capillary instability of Newtonian fluids was studied by Rayleigh [1] following earlier work by Plateau [] who showed that a long cylinder of liquid is unstable to disturbances with wave lengths greater than πr. The analysis of Rayleigh is based on potential flow of an inviscid liquid. Tomotika [3] studied the capillary instability and gave an exact Corresponding author. Tel.: ; fax: address: joseph@aem.umn.edu (D.D. oseph. normal mode solution of the linearized Navier Stokes equations. The linear analysis of capillary instability of viscoelastic fluids has been done by Middleman [4], Goldin et al. [5], Goren and Gottlieb [6]. They showed that the growth rates are larger for the viscoelastic fluids than for the equivalent Newtonian fluids. Funada and oseph [7,8] presented potential flow analyses of capillary instability of viscous and viscoelastic fluids. In their studies, the flow is assumed to be irrotational but the viscous and viscoelastic effects are retained (viscous or viscoelastic potential flow, PF. The viscous and viscoelastic stresses enter into the analyses through the normal stress balance at the interface. Funada and oseph compared their results based on potential flow to the unapproximated normal mode results (Tomotika [3]. They showed that the results with viscous and viscoelastic effects retained are in /$ see front matter 005 Elsevier B.. All rights reserved. doi: /j.jnnfm
2 . Wang et al. /. Non-Newtonian Fluid Mech. 19 ( where u = ue r + we z is the velocity, ρ is the density, p is the pressure, and τ is the extra stress. The extra stress may be modeled by effreys model τ + λ 1 τ = µ ( D + λ D, (3 Fig. 1. Capillary instability. The force γ/r drives fluid away from the throat, leading to collapse. better agreement with the unapproximated results than those assuming inviscid fluids. The capillary instability can be viewed as a free surface problem when either the interior or the exterior fluid is a gas of negligible density and viscosity. One difficulty in the potential flow analyses of free surface problems is that the non-zero irrotational shear stress violates the zero-shear-stress condition at the free surface. oseph and Wang [9] derived an additional viscous correction for the irrotational pressure, which arises in the boundary layer induced by the unphysical discontinuity of the shear stress. Wang et al. [10] applied this additional pressure contribution to the potential flow analysis of capillary instability of Newtonian fluids. They showed that the results computed with the additional pressure contribution are almost indistinguishable from the exact results. Here we extend the formulation for the additional pressure correction to potential flows of viscoelastic fluids in flows governed by linearized equations (viscoelastic correction of viscoelastic potential flow, CPF, and apply this additional pressure correction to capillary instability of viscoelastic liquid filaments of effreys type. The results are in remarkably good agreement with those obtained from the unapproximated normal mode analysis for viscoelastic fluids. The linear stability analysis given here and elsewhere indicates that the liquid jets are less stable with increasing elasticity, which contradicts the observation in experiments. A possible explanation of this contradiction is related to the linear stability analysis of a stressed filament at rest (Entov [11]. One difficulty is that a stressed filament at rest is not a permanent solution.. Linear stability equations and the exact solution (Tomotika [3] In an undisturbed rest state, the long cylinder of a viscoelastic liquid is surrounded by a gas of negligible density and viscosity. We use cylindrical coordinates (r, θ, z and consider small axisymmetric disturbances. The linearized governing equations of the interior liquid are u = 0, (1 ρ u = p + τ ( where D is the rate of strain tensor, µ is the viscosity, λ 1 and λ are the relaxation and retardation times, respectively. Suppose that we have normal mode solutions with the growth rate σ: τ = exp(σt τ and D = exp(σt D, (4 then (3 leads to τ = 1 + λ σ 1 + λ 1 σ µ D τ = 1 + λ σ µd. 1 + λ 1 σ (5 The momentum Eq. ( becomes ρ u ( 1 + λ σ = p λ 1 σ µd = p λ σ 1 + λ 1 σ µ u. (6 The shear and normal stress boundary conditions are ( 1 + λ σ u 1 + λ 1 σ µ z + w = 0; r (7 p λ σ 1 + λ 1 σ µ u r = γ ( η z + η R, (8 where η is the varicose displacement. The governing Eqs. (1 and (6 and boundary conditions (8 and (7 are the same as those for a Newtonian fluid except that (1 + λ σµ/(1 + λ 1 σ replaces µ. Following scales are used to construct dimensionless governing equations: the cylinder diameter D for length, U = γ/(ρd for velocity, T = D/U for time and p0 = ρu for pressure. The dimensionless momentum equation is (we use the same symbols for dimensionless variables u = p + ˆµ u, (9 where ˆµ = 1 + ˆλ σ ( ˆλ 1 σ with ˆλ 1 = λ 1 U D = λ 1 and γ U γ ρd 3 and ˆλ = λ D = λ ρd 3, (11 = ργd µ (1
3 108. Wang et al. /. Non-Newtonian Fluid Mech. 19 ( is the Reynolds number and 1/ is the Ohnesorge number. The dimensionless boundary conditions at the cylinder surface R = 0.5 are p + ˆµ u r = η z + η R ; (13 ( ˆµ u z + w = 0. (14 r A solution of (9 which satisfies both the boundary conditions (13 and (14 takes the following form: ψ = [A 1 ri 1 (kr + A ri 1 (k v r] exp(σt + ikz, w = 1 r u = 1 r ψ z, ψ r, (15 η = H exp(σt + ikz, (16 where k is the wave number and I 1 denotes the first kind modified Bessel function of the first order. Substitution of (15 and (16 into (13 and (14 leads to the solvability condition, which is given as the dispersion relation of σ: k I 1 (kr (k + kv I 1(k v R F 1 F = 0 (17 where ( F 1 = σi 0 (kr + ˆµk di1 (kr d(kr ( 1 k R k σ I 1(kR, (18 F = ˆµkk ( v di1 (k v R d(k v R ( 1 k R k σ I 1(k v R, (19 k + ( /ˆµσ. This solution satisfies the gov- with k v = erning equations and all the boundary conditions and is an exact solution iscoelastic potential flow (PF It is easy to show that the momentum Eq. (9 admits potential flow solutions. Take curl of Eq. (9 and use u = φ, we obtain φ = ( p + ˆµ φ. (0 Both sides of (0 are zero, therefore potential flow solutions are compatible in this problem. The pressure integral can also be easily obtained from (9, 1 In our former paper [10], the exact solution of the linearized equations was called the fully viscous flow (FF solution. ( φ = p p + ˆµ φ p p = φ, (1 where p p denotes the pressure from the potential flow solution and it is equal to the pressure from the inviscid potential flow. The potential flow solution is given by φ = AiI 0 (kr exp(σt + ikz, u = φ r, w = φ z, ( η = H exp(σt + ikz. (3 Substitution of the potential flow solution into the normal stress balance (13 leads to the dispersion relation I 0 (kr ( I 1 (kr σ 1 + ˆλ 1 σ + ( 1 + ˆλ σ [ σ k I0 (kr I 1 (kr 1 ] kr ( 1 (1 k R k + ˆλ 1 σ = 0, (4 which is a cubic equation of σ and has explicit solutions. When, Eq. (4 reduces to ( I 0 (kr 1 I 1 (kr σ = k R k, (5 which is the dispersion relation for inviscid potential flow (IPF solution. The IPF solution does not allow viscous or viscoelastic effects. 4. Dissipation and the formulation for the additional pressure contribution oseph and Wang [9] derived a viscous pressure contribution in addition to the irrotational pressure for the potential flow solutions of Newtonian fluids by considering the dissipation of energy. Here we extend the analysis to a viscoelastic fluid of effereys type in flows governed by linearized equations. We start from the momentum equation ρ du dt = T u ρ du dt where T is the total stress. It follows that ρ d ( 1 dt u u = ( T u, (6 = (T u u : T = (T u (D + :( p1 + ˆµµD = (T u D : ( p1 + ˆµµD = (T u ˆµµD : D. It follows that d dt = ( ρ u u d n (T udω ˆµµ D : D d, (7
4 . Wang et al. /. Non-Newtonian Fluid Mech. 19 ( where is the volume occupied by the viscoelastic fluid, Ω is the boundary of, and n is the outward normal of on Ω. We have shown that the potential flow is a solution of the momentum equation in this problem. Thus we can insert the velocity and stress tensor evaluated on the potential flow into (7 to obtain d dt = ( ρ u u d Ω [( p p + τ rr u + τ rz w]dω ˆµµ D : D d. (8 At the free surface, the potential flow leads to a non-zero irrotational shear stress and does not satisfy the zero-shearstress condition. We introduce a pressure contribution p c in addition to the irrotational pressure p p ; p c cancels out the power due to the unphysical irrotational shear stress in the energy equation and (7 becomes d dt = ( ρ u u d Ω [( p p p c + τ rr u]dω ˆµµ D : D d. (9 Comparing (8 and (9, we obtain τ rz wdω = ( p c u dω, (30 Ω Ω which is the same as the formulation for the additional pressure contribution as in the potential flow of a viscous Newtonian fluid (oseph and Wang [9]. However, the calculation of τ rz in viscoelastic fluids is different than in Newtonian fluids. The additional pressure contribution p c depends strongly on viscoelastic parameters and is determined solely by the irrotational flow. 5. The additional pressure contribution for capillary instability Now we consider the additional pressure contribution for the potential flow analysis of capillary instability. oseph and Wang [9] showed that in linearized problems, the governing equation for the additional pressure contribution is p c = 0. (31 It is easy to show that (31 holds for the viscoelastic fluid under consideration here. Solving (31, we obtain p c = C j ii 0 (jr exp(σt + ijz, (3 j=0 where C j are constants. With the additional pressure contribution, the normal stress balance becomes p p p c + ˆµ 1 u r = η z + η R, (33 which gives rise to [ {AσI 0 (kr + C k I 0 (kr + ˆµk A I 0 (kr I ]} 1(kR kr exp(σt + ikz + j k C j I 0 (jr exp(σt + ijz = A k ( 1 σ I 1(kR R k exp(σt + ikz. (34 By orthogonality of Fourier series, C j = 0ifj k. The coefficient C k can be determined using (30. The left hand side of (30 is τ rz w dω = ˆµ 4πlRAA k 3 I 0 (kri 1 (kr exp(σ + σ t, (35 where l is the length of one wave period and denotes conjugate variables. On the other hand, ( p c u dω = πlrc k A ki 0 (kri 1 (kr exp(σ + σ t. (36 It follows that C k = ( ˆµ/ Ak and p c = iak ˆµ I 0 (kr exp(σt + ikz. (37 Inserting C k into (34, we obtain [ σi 0 (kr + ˆµk I 0 (kr + ˆµk I 0 (kr I ] 1(kR kr = k ( 1 σ I 1(kR R k, which can be written as I 0 (kr ( I 1 (kr σ 1 + ˆλ 1 σ + ( 1 + ˆλ σ [ σ k I0 (kr I 1 (kr 1 ] kr ( 1 (1 k R k + ˆλ 1 σ = 0. (38 Eq. (38 is the dispersion relation from the viscoelastic correction of PF (CPF. If the pressure correction (37 is inserted back into the governing Eq. (9, we obtain u c = p c + ˆµ u c, (39
5 110. Wang et al. /. Non-Newtonian Fluid Mech. 19 ( where u c is the velocity correction induced by the pressure correction p c. We can find a potential flow solution u c = φ c, such that u c = φ c = 0 and φ c = p c. (40 It can be readily shown that φ c = i ˆµ Ak I 0 (kr exp(σt + ikz. (41 σ Thus the pressure correction p c which is proportional to 1/ induces a velocity correction proportional to 1/. This velocity correction gives rise to uncompensated shear stress proportional to 1 which may induce a new pressure correction now proportional to 1. In this way we may generate, successively, irrotational solutions proportional to increasing powers of 1/. We believe that only the first pressure correction proportional to 1/ is of physical significance; the higher order corrections are not considered in the normal stress balance ( Comparison of the growth rate We compare the dispersion relation (38 from CPF with (4 from PF, (5 from IPF and (17 from the exact solution. Eqs. (17, (4, (5 and (38 are solved by numerical methods for the growth rate σ and the values of σ are compared. First we examine two practical cases: % PAA in air and % PEO in air (following Funada and oseph [8]. We choose the diameter of the fluid cylinder to be 1 cm. The σ versus k plots for % PAA and % PEO are shown in Figs. and 3, respectively. These figures show that the results from CPF are almost indistinguishable from the exact solution, whereas IPF and PF overestimates σ significantly. Fig. 3. The growth rate σ vs. k from inviscid potential flow (IPF, viscoelastic potential flow (PF, viscoelastic correction of PF (CPF and the exact solution. The results of the exact solution and CPF are almost the same. The fluid is % PEO, ρ = 0.99gcm 1, µ = 350 P, γ = 63.0 dyn cm 1, λ 1 = 0.1 s, λ = 0s, = , ˆλ 1 = Capillary instability is controlled by three dimensionless numbers:, ˆλ 1, and ˆλ. We vary these parameters and present the computed growth rate in Figs The Reynolds number ranges from 10 4 to 10 4, ˆλ 1 ranges from 0.1 to 1000, and ˆλ ranges from 0 to 100. In all the cases, the growth rates from CPF are in excellent agreement with the exact solution, indicating that our additional pressure contribution is valid for a wide range of controlling parameters. Figs. 4 and 5 show that the growth rates increase with ˆλ 1 when and ˆλ are fixed. Comparing Figs. 5 and 6, it can be seen that the effect of ˆλ is opposite to that of ˆλ 1 ; the growth rates decreases with ˆλ. When ˆλ 1 = ˆλ, the fluid becomes Newtonian. When the Reynolds number is as high as 10 4 (Fig. 7, IPF and PF slightly over-estimate the maximum growth rate whereas the CPF results are almost the same as the exact solution. Fig.. The growth rate σ vs. k from inviscid potential flow (IPF, viscoelastic potential flow (PF, viscoelastic correction of PF (CPF and the exact solution. The growth rates for the exact solution and CPF are almost the same. The fluid is % PAA, ρ = 0.99gcm 1, µ = 96 P, γ = 45.0 dyn cm 1, λ 1 = s, λ = 0s, = , ˆλ 1 = Fig. 4. The growth rate σ vs. k from inviscid potential flow (IPF, viscoelastic potential flow (PF, viscoelastic correction of PF (CPF the exact solution. = 10 4, ˆλ 1 = 0.1, ˆλ = 0.
6 . Wang et al. /. Non-Newtonian Fluid Mech. 19 ( Fig. 5. The growth rate σ vs. k from inviscid potential flow (IPF, viscoelastic potential flow (PF, viscoelastic correction of PF (CPF and the exact solution. = 10 4, ˆλ 1 = 1000, ˆλ = 0. Fig. 7. The growth rate σ vs. k from inviscid potential flow (IPF, viscoelastic potential flow (PF, viscoelastic correction of PF (CPF and the exact solution. = 10 4, ˆλ 1 = 0.1, ˆλ = 0. When the Reynolds number is large, viscoelastic effects are relatively small, and the four curves are close; but differences among them can be seen near the peak growth rate. The inset is the amplified plot for the region near the peak growth rate. CPF is the best approximation to the exact solution. 7. Comparison of the stream function Fig. 6. The growth rate σ vs. k from inviscid potential flow (IPF, viscoelastic potential flow (PF, viscoelastic correction of PF (CPF and the exact solution. = 10 4, ˆλ 1 = 1000, ˆλ = 100. In Table 1 we present the maximum growth rate σ m and the associated wave number k m computed from PF, CPF and the exact solution. The value of σ m given by PF is several times larger than the exact result when is small. CPF gives excellent approximation to the values of σ m and k m in all the cases. Next we compare the stream functions from PF, CPF and the exact solution at the same wave number. The wave number chosen for the comparison is k m at which the maximum growth rate σ m occurs in the exact solution. The relation between the constants A 1 and A in the exact stream function (15 and A in the potential flow solution ( must be established before one can compare the stream functions. Here this relation is obtained by assuming that the magnitude of the disturbance H is the same in the exact solution and in the potential flow solution. We use a superscript E for quantities appearing in the exact solution and (15 and (16 are rewritten as ψ E = [A E 1 ri 1(kr + A E ri 1(k v r] exp(σ E t + ikz, (4 η E = H E exp(σ E t + ikz. (43 Table 1 Maximum growth rate σ m and the associated wave number k m for viscoelastic potential flow (PF, viscoelastic correction of PF (CPF and the exact solution in Figs. 7 Figure PF CPF Exact solution k m σ m k m σ m k m σ m 3.439e e 0.05e 01.74e 0.135e 01.78e e 01.83e e e e e e e e 0 3.3e e e e e e e e e e e 0.101e e 0.186e e e e e e e e 01 For inviscid potential flow (IPF solution, k m = and σ m = in all the six cases.
7 11. Wang et al. /. Non-Newtonian Fluid Mech. 19 ( Table The growth rate σ computed from viscoelastic potential flow (PF, viscoelastic correction of PF (CPF and the exact solution at the same wave number k m ˆλ 1 ˆλ k m σ PF σ CPF σ E In the exact solution, k m is the wave number for the maximum growth rate. The relation between A E 1 and AE is determined by the zeroshear-stress condition at r R: A E = k I 1 (kr (k + kv I 1(k v R AE 1. (44 Therefore, we can write the stream function as ψ E = A E 1 [I r k ] I 1 (kr 1 (kr (k + kv I 1(k v R I 1(k v r exp(σ E t + ikz. (45 The amplitude of the disturbance H E is related to A E 1 through the kinematic condition: H E = ik (1 σ E k k + kv I 1 (kra E 1. (46 Now we consider the potential flow solution which is indicated by a superscript P. The stream function and the disturbance are given by ψ P = A P ri 1 (kr exp(σ P t + ikz, (47 η P = H P exp(σ P t + ikz, (48 respectively. The amplitude of the disturbance H P is related to A P through the kinematic condition: H P = ik σ P AP I 1 (kr. (49 We assume that the amplitude of the disturbance is the same in the exact solution and the potential flow solution. Thus H E = H P and it follows that A P = A E σ P 1 (1 σ E k k + kv. (50 Then the stream function of the potential flow can be written as ψ P = A E σ P 1 (1 σ E k k + kv ri 1 (kr exp(σ P t + ikz. (51 Now we can compare (45 and (51. The stream function is decomposed into two parts, the exponential function depending on t and z and the rest part depending on r. Since we are comparing the stream functions at the same wave number k m, the comparison of the exponential function is equivalent to the comparison of the growth rate. In Table, we list the values of the growth rate σ computed from PF, CPF and the exact solution. In all the cases, the growth rate from PF is larger than the exact result, whereas the growth rate from CPF is very close to the exact result. The rest part of the stream function depends on r and we define SF(r = σpf σ E SF(r = σcpf σ E (1 k k + k v (1 k k + k v ri 1 (kr for PF; (5 ri 1 (kr [ k ] I 1 (kr SF(r = r I 1 (kr (k + kv I 1(k v R I 1(k v r for CPF; (53 for the exact solution. (54 Three examples for the comparison of the function SF(r are shown in Figs The curves for SF(r are very close to straight lines, indicating power functions. This can also be Fig. 8. The part of the stream function depending on r defined in (5 (54 for viscoelastic potential flow (PF, viscoelastic correction of PF (CPF and the exact solution respectively. The fluid is % PAA, = , ˆλ 1 = 0.63, ˆλ = 0. The wave number for the maximum growth rate k m = is chosen for the comparison.
8 . Wang et al. /. Non-Newtonian Fluid Mech. 19 ( SF(r = (1 σcpf σ E k kr k + kv + O(r4 for CPF. (57 For (54, we also expand I 1 (kr and I 1 (k v R and keep only the leading term, which gives rise to SF(r = kr k (kr/ k v r (k + kv (k vr/ + O(r4 = (1 k kr k + kv + O(r4 for the exact solution. (58 Fig. 9. The part of the stream function depending on r defined in (5 (54 for viscoelastic potential flow (PF, viscoelastic correction of PF (CPF and the exact solution respectively. In this case, = 10 4, ˆλ 1 = 1000, ˆλ = 0. The wave number for the maximum growth rate k m = is chosen for the comparison. seen from (5 (54. The expansion of the modified Bessel function gives I 1 (kr = kr + k3 r k5 r O(r7. (55 Higher order terms of r may be neglected because 0 r 0.5 inside the cylinder. If we only keep the first term in the expansion, the stream functions (5 and (53 become, respectively SF(r = σpf σ E (1 k k + k v kr + O(r4 for PF; (56 Eqs. (56 (58 show that the functions SF(r are approximately quadratic functions for small r and this is confirmed in Figs The comparison of the leading terms of SF(r depends directly on the growth rate σ PF, σ CPF and σ E. Since σ PF >σ E, the curves for SF(r of PF are higher than those for the exact solution. On the other hand, σ CPF is very close to σ E and the curves for CPF and the exact solution almost overlap. Combining the comparison of the growth rate in Table and the comparison of the function SF(r infigs. 8 10, we show that the stream function given by CPF is in remarkably good agreement with the exact solution. This result indicates that the vorticity plays a small role in the exact solution and our CPF solution, which is based solely on potential flow, can give an excellent approximation to the flow. 8. Discussion Chang et al. [1] did a long wave study of the stretching dynamics of bead-string filaments for FENE and Oldroyd- B fluids. They also did a long wave study of linear stability and their results can be compared to ours. To this end, we first convert the parameters used by Chang et al. to the parameters used by us. In the notation of Chang et al. Ca is the capillary number, We is the Weissenberg number and S is the retardation number. We linearize the stress equation of Chang et al. and reduce it to a form comparable to our effreys model (3, then the relation between We and S used by Chang et al. and ˆλ 1 and ˆλ used by us is revealed. After taking the different length and time scales into account, we can express the parameters in Chang et al. in terms of our parameters Ca =, We = 4ˆλ 1, and S = ˆλ ˆλ 1. (59 Fig. 10. The part of the stream function depending on r defined in (5 (54 for viscoelastic potential flow (PF, viscoelastic correction of PF (CPF and the exact solution respectively. In this case, = 10 4, ˆλ 1 = 1000, ˆλ = 100. The wave number for the maximum growth rate k m = is chosen for the comparison. Then the dispersion relation given by the linear stability analysis of Chang et al. (their Eq. (16 can be written as ( ˆλ 1 σ k ˆλ [ 3k σ ] + k ˆλ 1 4 (4 k σ k 4 (4 k = 0. (60
9 114. Wang et al. /. Non-Newtonian Fluid Mech. 19 ( Now we consider the dispersion relation (38 from the CPF method. The dimensionless radius R = 1/ and the Bessel functions can be expanded for small k I 0 (kr I 1 (kr = 4 k + k 8 k O(k5, I 0 (kr I 1 (kr 1 kr = 6 k + k 4 k O(k5. (61 Inserting (61 into (38, we can obtain (1 + k + ˆλ 1 σ 3 + [1 + k 3 [ 3k ( 1 + k ˆλ ( 3k 1 + k 4 ] k ˆλ 1 4 (4 k σ ] σ k 4 (4 k + O(k 4 = 0. (6 The expansion of the Bessel functions can also be applied to the exact solution and the result will be compared to (60 and (6. After some arrangement, the dispersion relation (17 of the exact solution can be written as 4k 3 [ k v I0 (k v R ˆµ I 1 (k v R 1 k v R [ I0 (kr I 1 (kr 1 ] k3 kr σ + (k + kv k ( 1 σ R k ] k (k + k v ˆµ ( 1 R k (k + k v σ I 0(kR I 1 (kr = 0. (63 After expanding the Bessel functions as power series of k,we obtain ( ˆλ k σ 3 + (1 + k ˆλ 3k σ [ 3k ] + k ˆλ 1 4 (4 k σ k 4 (4 k + O(k 4 = 0. (64 The dispersion relation (64 given by the exact solution is different from both (60 given by Chang et al. and (6 given by our CPF method; the first order differences are O(k in both cases. The differences between (64 and (60 are two k /3 terms in the coefficients of σ 3 and σ ; the differences between (64 and (6 are two k /4 terms in the coefficients of σ and σ. The limit of a Newtonian fluid can be obtained by letting ˆλ 1 = ˆλ = 0. Then the dispersion relations (60, (6 and (64 reduce to, respectively σ + 3k σ k 4 (4 k = 0 for Chang et al.; (65 Fig. 11. The growth rate σ as a function of k computed using (60 given by Chang et al., (38 given by our CPF method, and (17 given by the exact solution. The fluid is water with ρ = 1000 kg/m 3, µ = kg/(m s, γ = N/m. The diameter of the liquid cylinder is assumed to be 0.01 m and the Reynolds number is = ( (1 + k σ + 3k 1 + k σ k (4 k + O(k 4 = 0 for CPF; (66 (1 + k σ + 3k σ k 3 4 (4 k + O(k 4 = 0 for the exact solution. (67 The first order differences among the dispersion relations (65 (67 are O(k. The difference between (65 and (67 is a k /3 term in the coefficient of σ ; the difference between (66 and (67 is a k /4 term in the coefficient of σ. In Figs , we plot the growth rate σ as a function of k computed using (60 given by Chang et al., (38 given by our Fig. 1. The growth rate σ as a function of k computed using (60 given by Chang et al., (38 given by our CPF method, and (17 given by the exact solution. The fluid is a Newtonian fluid SO10000 oil with ρ = 969 kg/m 3, µ = 10 kg/(m s, γ = 0.01 N/m. The Reynolds number is =
10 . Wang et al. /. Non-Newtonian Fluid Mech. 19 ( Fig. 13. The growth rate σ as a function of k computed using (60 given by Chang et al., (38 given by our CPF method, and (17 given by the exact solution. The fluid is % PAA with = , ˆλ 1 = 0.63, ˆλ = 0. CPF method, and (17 given by the exact solution. Both Newtonian fluids and viscoelastic fluids are compared. The limit of Newtonian fluids is achieved by setting ˆλ 1 and ˆλ to be zero in (60, (38 and (17. There is almost no difference between the three curves when k is close to zero, and small differences can be seen when k is close to. The dispersion relation of Chang et al. is in better agreement with the exact solution when the Reynolds number is small (Figs. 1 15, whereas our CPF is in better agreement with the exact solution when is large (Fig. 11. In this work, linear stability analysis of the capillary instability of a viscoelastic thread is carried out under the assumption that the flow is irrotational. The non-zero irrotational shear stress at the surface of the liquid thread does not agree with the zero-shear-stress condition. We derive a pressure contribution in addition to the irrotational pressure. This additional pressure contribution depends on the viscoelastic pa- Fig. 15. The growth rate σ as a function of k computed using (60 given by Chang et al., (38 given by our CPF method, and (17 given by the exact solution. The fluid is a viscoelastic fluid with = 10 4, ˆλ 1 = 1000, ˆλ = 100. rameters and cancels out the power due to the uncompensated irrotational shear stress in the energy equation. We include the additional pressure contribution, the irrotational pressure and the extra stress evaluated using the irrotational flow in the normal stress balance at the surface, then a dispersion relation is obtained. We call this approach as the viscoelastic correction of the viscoelastic potential flow (CPF. The comparison of the growth rate and the stream function show that the CPF solution is an excellent approximation to the exact solution. The dispersion relation given by CPF is also compared to that obtained by Chang et al. [1] using a long wave approximation. The differences between the two dispersion relations are negligible when the wave number k is small and both dispersion relations are in remarkably good agreement with the exact solution. Acknowledgement This work was supported in part by the NSF under grants from Chemical Transport Systems. References Fig. 14. The growth rate σ as a function of k computed using (60 given by Chang et al., (38 given by our CPF method, and (17 given by the exact solution. The fluid is a viscoelastic fluid with = 10 4, ˆλ 1 = 1000, ˆλ = 0. [1] L. Rayleigh, On the capillary phenomena of jets, Proc. R. Soc. London A 9 ( [].A.F. Plateau, Statique experimentale et theorique des liquide soumis aux seules forces moleculaire, vol. ii, 1873, p. 31. [3] S. Tomotika, On the instability of a cylindrical thread of a viscous liquid surrounded by another viscous fluid, Proc. R. Soc. London A 150 ( [4] S. Middleman, Stability of a viscoelastic jet, Chem. Eng. Sci. 0 ( [5] M. Goldin,. Yerushalmi, R. Pfeffer, R. Shinnar, Breakup of a laminar capillary jet of a viscoelastic fluid,. Fluid Mech. 38 (
11 116. Wang et al. /. Non-Newtonian Fluid Mech. 19 ( [6] S.L. Goren, M. Gottlieb, Surface-tension-driven breakup of viscoelstic liquid threads,. Fluid Mech. 10 ( [7] T. Funada, D.D. oseph, iscous potential flow analysis of capillary instability, Int.. Multiphase Flow 8 ( [8] T. Funada, D.D. oseph, iscoelastic potential flow analysis of capillary instability,. Non-Newtonian Fluid Mech. 111 ( [9] D.D. oseph,. Wang, The dissipation approximation and viscous potential flow,. Fluid Mech. 505 ( [10]. Wang, D.D. oseph, T. Funada, Pressure correction for potential flow analysis of capillary instability of viscous fluids,. Fluid Mech. 5 ( [11].M. Entov, On the stability of capillary jets of elastoviscous fluids,. Eng. Phys. Thermophys. 34 ( [1] H.-C. Chang, E.A. Demekhin, E. Kalaidin, Iterated stretching of viscoelastic jets, Phys. Fluids 11 (
Pressure corrections for viscoelastic potential flow analysis of capillary instability
ve-july29-4.tex 1 Pressure corrections for viscoelastic potential flow analysis of capillary instability J. Wang, D. D. Joseph and T. Funada Department of Aerospace Engineering and Mechanics, University
More informationPurely irrotational theories of the effects of viscosity and viscoelasticity on capillary instability of a liquid cylinder
ve-march17.tex 1 Purely irrotational theories of the effects of viscosity and viscoelasticity on capillary instability of a liquid cylinder J. Wang, D. D. Joseph and T. Funada Department of Aerospace Engineering
More informationViscous contributions to the pressure for potential flow analysis of capillary instability of two viscous fluids
PHYSICS OF FLUIDS 17, 052105 2005 Viscous contributions to the pressure for potential flow analysis of capillary instability of two viscous fluids. Wang and D. D. oseph a Department of Aerospace Engineering
More informationLinear analysis of three-dimensional instability of non-newtonian liquid jets
J. Fluid Mech. (2006), vol. 559, pp. 451 459. c 2006 Cambridge University Press doi:10.1017/s0022112006000413 Printed in the United Kingdom 451 Linear analysis of three-dimensional instability of non-newtonian
More informationViscoelastic Potential Flow Analysis of Capillary Instability
VE-CapInstab/viscoelastic-Rev2-NNFM.tex Viscoelastic Potential Flow Analysis of Capillary Instability T. Funada a and D.D. oseph b a Department of Digital Engineering, Numazu College of Technology, Ooka
More informationHigher Orders Instability of a Hollow Jet Endowed with Surface Tension
Mechanics and Mechanical Engineering Vol. 2, No. (2008) 69 78 c Technical University of Lodz Higher Orders Instability of a Hollow Jet Endowed with Surface Tension Ahmed E. Radwan Mathematics Department,
More information2 D.D. Joseph To make things simple, consider flow in two dimensions over a body obeying the equations ρ ρ v = 0;
Accepted for publication in J. Fluid Mech. 1 Viscous Potential Flow By D.D. Joseph Department of Aerospace Engineering and Mechanics, University of Minnesota, MN 55455 USA Email: joseph@aem.umn.edu (Received
More informationPlateau-Rayleigh Instability of a Cylinder of Viscous Liquid (Rayleigh vs. Chandrasekhar) L. Pekker FujiFilm Dimatix Inc., Lebanon NH USA
Plateau-Rayleigh Instability of a Cylinder of Viscous Liquid (Rayleigh vs. Chandrasekhar) L. Pekker FujiFilm Dimatix Inc., Lebanon NH 03766 USA Abstract In 1892, in his classical work, L. Rayleigh considered
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 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 informationLiquid Jet Breakup at Low Weber Number: A Survey
International Journal of Engineering Research and Technology. ISSN 0974-3154 Volume 6, Number 6 (2013), pp. 727-732 International Research Publication House http://www.irphouse.com Liquid Jet Breakup at
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 informationViscous potential flow analysis of stress induced cavitation in an aperture flow
cavitation-july9.tex Viscous potential flow analysis of stress induced cavitation in an aperture flow T. Funada, J. Wang and D. D. Joseph Department of Digital Engineering, Numazu College of Technology,
More informationCapillary Collapse and Rupture. Abstract
Capillary Collapse and Rupture T.S. Lundgren University of Minnesota, Minneapolis, MN 55455 D.D. Joseph University of Minnesota, Minneapolis, MN 55455 (June 9, 1998) Abstract The breakup of a liquid capillary
More informationStability of a Liquid Jet into Incompressible Gases and Liquids
LiquidJet-06-25bw.tex 1 Stability of a Liquid Jet into Incompressible Gases and Liquids T. Funada, 1 D.D. Joseph 23 and S. Yamashita 1 1 Department of Digital Engineering, Numazu College of Technology,
More informationPredicting Breakup Characteristics of Liquid Jets Disturbed by Practical Piezoelectric Devices
ILASS Americas 2th Annual Conference on Liquid Atomization and Spray Systems, Chicago, IL, May 27 Predicting Breakup Characteristics of Liquid Jets Disturbed by Practical Piezoelectric Devices M. Rohani,
More informationOn the displacement of two immiscible Stokes fluids in a 3D Hele-Shaw cell
On the displacement of two immiscible Stokes fluids in a 3D Hele-Shaw cell Gelu Paşa Abstract. In this paper we study the linear stability of the displacement of two incompressible Stokes fluids in a 3D
More informationMODELING ON THE BREAKUP OF VISCO-ELASTIC LIQUID FOR EFFERVESCENT ATOMIZATION
1446 THERMAL SCIENCE, Year 2012, Vol. 16, No. 5, pp. 1446-1450 MODELING ON THE BREAKUP OF VISCO-ELASTIC LIQUID FOR EFFERVESCENT ATOMIZATION by Li-Juan QIAN * China Jiliang University, Hangzhou, China Short
More informationViscosity of magmas containing highly deformable bubbles
Journal of Volcanology and Geothermal Research 105 (2001) 19±24 www.elsevier.nl/locate/jvolgeores Viscosity of magmas containing highly deformable bubbles M. Manga a, *, M. Loewenberg b a Department of
More informationPotential flow of a second-order fluid over a sphere or an ellipse
J. Fluid Mech. (4), vol. 5, pp. 5. c 4 Cambridge niversity Press DOI:.7/S4954 Printed in the nited Kingdom Potential flow of a second-order fluid over a sphere or an ellipse By J. W A N G AND D. D. J O
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 informationThe effects of gravity on the capillary instability in tubes
J. Fluid Mech. (2006), vol. 556, pp. 217 226. c 2006 Cambridge University Press doi:10.1017/s0022112006009505 Printed in the United Kingdom 217 The effects of gravity on the capillary instability in tubes
More information3 Generation and diffusion of vorticity
Version date: March 22, 21 1 3 Generation and diffusion of vorticity 3.1 The vorticity equation We start from Navier Stokes: u t + u u = 1 ρ p + ν 2 u 1) where we have not included a term describing a
More informationAcceleration of a liquid drop suddenly exposed to a high-speed airstream
International Journal of Multiphase Flow 30 (2004) 217 224 www.elsevier.com/locate/ijmulflow Brief communication Acceleration of a liquid drop suddenly exposed to a high-speed airstream C. Ortiz, D.D.
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 informationInterfacial hoop stress and viscoelastic free surface flow instability. Michael D. Graham University of Wisconsin-Madison
Interfacial hoop stress and viscoelastic free surface flow instability Michael D. Graham University of Wisconsin-Madison Free surface instabilities of viscoelastic flows Eccentric cylinders (Varela-Lopez
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 informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationIterated stretching of viscoelastic jets
PHYSICS OF FLUIDS VOLUME 11, NUMBER 7 JULY 1999 Iterated stretching of viscoelastic jets Hsueh-Chia Chang, a) Evgeny A. Demekhin, and Evgeny Kalaidin Department of Chemical Engineering, University of Notre
More informationLecture notes Breakup of cylindrical jets Singularities and self-similar solutions
Lecture notes Breakup of cylindrical jets Singularities and self-similar solutions by Stephane Popinet and Arnaud Antkowiak June 8, 2011 Table of contents 1 Equations of motion for axisymmetric jets...........................
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 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 informationReview of fluid dynamics
Chapter 2 Review of fluid dynamics 2.1 Preliminaries ome basic concepts: A fluid is a substance that deforms continuously under stress. A Material olume is a tagged region that moves with the fluid. Hence
More informationRayleigh Taylor instability of viscoelastic drops at high Weber numbers
J. Fluid Mech. (2002), vol. 453, pp. 109 132. c 2002 Cambridge University Press DOI: 10.1017/S0022112001006802 Printed in the United Kingdom 109 Rayleigh Taylor instability of viscoelastic drops at high
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 informationUnsteady Flow of a Newtonian Fluid in a Contracting and Expanding Pipe
Unsteady Flow of a Newtonian Fluid in a Contracting and Expanding Pipe T S L Radhika**, M B Srinivas, T Raja Rani*, A. Karthik BITS Pilani- Hyderabad campus, Hyderabad, Telangana, India. *MTC, Muscat,
More informationUNIVERSITY OF EAST ANGLIA
UNIVERSITY OF EAST ANGLIA School of Mathematics May/June UG Examination 2007 2008 FLUIDS DYNAMICS WITH ADVANCED TOPICS Time allowed: 3 hours Attempt question ONE and FOUR other questions. Candidates must
More informationIrrotational Faraday Waves on a Viscous Fluid
faraday-dec27.tex 1 Irrotational Faraday Waves on a Viscous Fluid T.Funada, J.Wang, D.D.Joseph, & N.Tashiro Department of Digital Engineering, Numazu National College of Technology, 36 Ooka, Numazu, Shizuoka,
More informationLecture 1. Hydrodynamic Stability F. H. Busse NotesbyA.Alexakis&E.Evstatiev
Lecture Hydrodynamic Stability F. H. Busse NotesbyA.Alexakis&E.Evstatiev Introduction In many cases in nature, like in the Earth s atmosphere, in the interior of stars and planets, one sees the appearance
More informationChapter 6: Incompressible Inviscid Flow
Chapter 6: Incompressible Inviscid Flow 6-1 Introduction 6-2 Nondimensionalization of the NSE 6-3 Creeping Flow 6-4 Inviscid Regions of Flow 6-5 Irrotational Flow Approximation 6-6 Elementary Planar Irrotational
More informationLecture 3: The Navier-Stokes Equations: Topological aspects
Lecture 3: The Navier-Stokes Equations: Topological aspects September 9, 2015 1 Goal Topology is the branch of math wich studies shape-changing objects; objects which can transform one into another without
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 informationEllipsoidal model of the rise of a Taylor bubble in a round tube
Ellipsoidal model of the rise of a Taylor bubble in a round tube T.Funada, 1 D.D.oseph, T.Maehara 1 and S.Yamashita 1 1 Department of Digital Engineering, Numazu College of Technology, 3600 Ooka, Numazu,
More informationApplication of He s homotopy perturbation method to boundary layer flow and convection heat transfer over a flat plate
Physics Letters A 37 007) 33 38 www.elsevier.com/locate/pla Application of He s homotopy perturbation method to boundary layer flow and convection heat transfer over a flat plate M. Esmaeilpour, D.D. Ganji
More informationVortex motion. Wasilij Barsukow, July 1, 2016
The concept of vorticity We call Vortex motion Wasilij Barsukow, mail@sturzhang.de July, 206 ω = v vorticity. It is a measure of the swirlyness of the flow, but is also present in shear flows where the
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
More information1. Introduction, tensors, kinematics
1. Introduction, tensors, kinematics Content: Introduction to fluids, Cartesian tensors, vector algebra using tensor notation, operators in tensor form, Eulerian and Lagrangian description of scalar and
More informationStability of Shear Flow
Stability of Shear Flow notes by Zhan Wang and Sam Potter Revised by FW WHOI GFD Lecture 3 June, 011 A look at energy stability, valid for all amplitudes, and linear stability for shear flows. 1 Nonlinear
More informationHydromagnetic Self-gravitating Stability of Streaming Fluid Cylinder. With Longitudinal Magnetic Field
Applied Mathematical Sciences, Vol. 8, 2014, no. 40, 1969-1978 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.312721 Hydromagnetic Self-gravitating Stability of Streaming Fluid Cylinder
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 informationChristel Hohenegger A simple model for ketchup-like liquid, its numerical challenges and limitations April 7, 2011
Notes by: Andy Thaler Christel Hohenegger A simple model for ketchup-like liquid, its numerical challenges and limitations April 7, 2011 Many complex fluids are shear-thinning. Such a fluid has a shear
More informationNumerical study of flow of Oldroyd-3-Constant fluids in a straight duct with square cross-section
Korea-Australia Rheology Journal Vol. 19, No. 2, August 2007 pp. 67-73 Numerical study of flow of Oldroyd-3-Constant fluids in a straight duct with square cross-section Mingkan Zhang, Xinrong Shen, Jianfeng
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 informationDESIGN OF ATOMIZERS AND BURNERS FOR COAL-WATER SLURRY COMBUSTION. Grant Number: DE-FG22-95PC95105 Progress Report for Period 1/1/97-3/31/97
DESIGN OF ATOMIZERS AND BURNERS FOR COAL-WATER SLURRY COMBUSTION Grant Number: DE-FG22-95PC95105 Progress Report for Period 1/1/97-3/31/97 A. Mansour, and N. Chigier Spray Systems Technology Center Dept.
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 informationRobust Monolithic - Multigrid FEM Solver for Three Fields Formulation Rising from non-newtonian Flow Problems
Robust Monolithic - Multigrid FEM Solver for Three Fields Formulation Rising from non-newtonian Flow Problems M. Aaqib Afaq Institute for Applied Mathematics and Numerics (LSIII) TU Dortmund 13 July 2017
More informationKelvin Helmholtz Instability
Kelvin Helmholtz Instability A. Salih Department of Aerospace Engineering Indian Institute of Space Science and Technology, Thiruvananthapuram November 00 One of the most well known instabilities in fluid
More informationMOMENTUM TRANSPORT Velocity Distributions in Turbulent Flow
TRANSPORT PHENOMENA MOMENTUM TRANSPORT Velocity Distributions in Turbulent Flow Introduction to Turbulent Flow 1. Comparisons of laminar and turbulent flows 2. Time-smoothed equations of change for incompressible
More informationDrop formation and breakup of low viscosity elastic fluids: Effects of molecular weight and concentration
PHYSICS OF FLUIDS 18, 043101 2006 Drop formation and breakup of low viscosity elastic fluids: Effects of molecular weight and concentration V. Tirtaatmadja Particulate Fluids Processing Centre, Department
More informationExam in Fluid Mechanics 5C1214
Eam in Fluid Mechanics 5C1214 Final eam in course 5C1214 13/01 2004 09-13 in Q24 Eaminer: Prof. Dan Henningson The point value of each question is given in parenthesis and you need more than 20 points
More informationNavier-Stokes Equation: Principle of Conservation of Momentum
Navier-tokes Equation: Principle of Conservation of Momentum R. hankar ubramanian Department of Chemical and Biomolecular Engineering Clarkson University Newton formulated the principle of conservation
More informationIntroduction to Fluid Mechanics
Introduction to Fluid Mechanics Tien-Tsan Shieh April 16, 2009 What is a Fluid? The key distinction between a fluid and a solid lies in the mode of resistance to change of shape. The fluid, unlike the
More informationStability of a viscoelastic thread immersed in a confined region
Stability of a viscoelastic thread immersed in a confined region A.Y. Gunawan Departemen Matematika, Institut Teknologi Bandung, Indonesia J. Molenaar, A.A.F. van de Ven Department of Mathematics and Computer
More information1. Comparison of stability analysis to previous work
. Comparison of stability analysis to previous work The stability problem (6.4) can be understood in the context of previous work. Benjamin (957) and Yih (963) have studied the stability of fluid flowing
More informationC. Project Description
C. Project Description The purpose of this proposal is to carry out the extension of mathematical studies of inviscid potential flow to viscous and viscoelastic potential flow. My claim is that inviscid
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 informationVisualization of polymer relaxation in viscoelastic turbulent micro-channel flow
Supplementary Information for Visualization of polymer relaxation in viscoelastic turbulent micro-channel flow Authors: J. Tai, C. P. Lim, Y. C. Lam Correspondence to: MYClam@ntu.edu.sg This document includes:
More informationUnit operations of chemical engineering
1 Unit operations of chemical engineering Fourth year Chemical Engineering Department College of Engineering AL-Qadesyia University Lecturer: 2 3 Syllabus 1) Boundary layer theory 2) Transfer of heat,
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 informationLiquid jet in a high Mach number air stream
jet-compressible-gas-july25.tex Liquid jet in a high Mach number air stream T.Funada, D.D.Joseph M.Saitoh & S.Yamashita Department of Digital Engineering, Numazu College of Technology, Ooka 36, Numazu,
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 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 informationOn the exact solution of incompressible viscous flows with variable viscosity
Advances in Fluid Mechanics IX 481 On the exact solution of incompressible viscous flows with variable viscosity A. Fatsis 1, J. Statharas, A. Panoutsopoulou 3 & N. Vlachakis 1 1 Technological University
More informationExcerpt from the Proceedings of the COMSOL Users Conference 2006 Boston
Using Comsol Multiphysics to Model Viscoelastic Fluid Flow Bruce A. Finlayson, Professor Emeritus Department of Chemical Engineering University of Washington, Seattle, WA 98195-1750 finlayson@cheme.washington.edu
More informationON THE EFFECTIVENESS OF HEAT GENERATION/ABSORPTION ON HEAT TRANSFER IN A STAGNATION POINT FLOW OF A MICROPOLAR FLUID OVER A STRETCHING SURFACE
5 Kragujevac J. Sci. 3 (29) 5-9. UDC 532.5:536.24 ON THE EFFECTIVENESS OF HEAT GENERATION/ABSORPTION ON HEAT TRANSFER IN A STAGNATION POINT FLOW OF A MICROPOLAR FLUID OVER A STRETCHING SURFACE Hazem A.
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 informationFlow Transition in Plane Couette Flow
Flow Transition in Plane Couette Flow Hua-Shu Dou 1,, Boo Cheong Khoo, and Khoon Seng Yeo 1 Temasek Laboratories, National University of Singapore, Singapore 11960 Fluid Mechanics Division, Department
More informationLecture 9 Laminar Diffusion Flame Configurations
Lecture 9 Laminar Diffusion Flame Configurations 9.-1 Different Flame Geometries and Single Droplet Burning Solutions for the velocities and the mixture fraction fields for some typical laminar flame configurations.
More informationQuick Recapitulation of Fluid Mechanics
Quick Recapitulation of Fluid Mechanics Amey Joshi 07-Feb-018 1 Equations of ideal fluids onsider a volume element of a fluid of density ρ. If there are no sources or sinks in, the mass in it will change
More informationFundamentals of Fluid Dynamics: Waves in Fluids
Fundamentals of Fluid Dynamics: Waves in Fluids Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/ tzielins/ Institute
More informationPaper ID ICLASS EXPERIMENTS ON BREAKUP OF WATER-IN-DIESEL COMPOUND JETS
ICLASS-2006 Aug.27-Sept.1, 2006, Kyoto, Japan Paper ID ICLASS06-047 EXPERIMENTS ON BREAKUP OF WATER-IN-DIESEL COMPOUND JETS Sheng-Lin Chiu 1, Rong-Horng Chen 2, Jen-Yung Pu 1 and Ta-Hui Lin 1,* 1 Department
More informationRheometry. II.1 Introduction
II Rheometry II.1 Introduction Structured materials are generally composed of microstructures dispersed in a homogeneous phase [30]. These materials usually have a yield stress, i.e. a threshold stress
More informationAdvantages of a Finite Extensible Nonlinear Elastic Potential in Lattice Boltzmann Simulations
The Hilltop Review Volume 7 Issue 1 Winter 2014 Article 10 December 2014 Advantages of a Finite Extensible Nonlinear Elastic Potential in Lattice Boltzmann Simulations Tai-Hsien Wu Western Michigan University
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
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 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 informationFLUID MECHANICS. Chapter 9 Flow over Immersed Bodies
FLUID MECHANICS Chapter 9 Flow over Immersed Bodies CHAP 9. FLOW OVER IMMERSED BODIES CONTENTS 9.1 General External Flow Characteristics 9.3 Drag 9.4 Lift 9.1 General External Flow Characteristics 9.1.1
More informationLiquid Sloshing in a Rotating, Laterally Oscillating Cylindrical Container
Universal Journal of Mechanical Engineering 5(3): 97-101, 2017 DOI: 10.13189/ujme.2017.050304 http://www.hrpub.org Liquid Sloshing in a Rotating, Laterally Oscillating Cylindrical Container Yusuke Saito,
More informationDetailed Outline, M E 521: Foundations of Fluid Mechanics I
Detailed Outline, M E 521: Foundations of Fluid Mechanics I I. Introduction and Review A. Notation 1. Vectors 2. Second-order tensors 3. Volume vs. velocity 4. Del operator B. Chapter 1: Review of Basic
More informationTutorial 10. Boundary layer theory
Tutorial 10 Boundary layer theory 1. If the velocity distribution law in a laminar boundary layer over a flat plate is assumes to be of the form, determine the velocity distribution law. At y = 0, u= 0
More informationExperiments at the University of Minnesota (draft 2)
Experiments at the University of Minnesota (draft 2) September 17, 2001 Studies of migration and lift and of the orientation of particles in shear flows Experiments to determine positions of spherical
More informationTurbulence - Theory and Modelling GROUP-STUDIES:
Lund Institute of Technology Department of Energy Sciences Division of Fluid Mechanics Robert Szasz, tel 046-0480 Johan Revstedt, tel 046-43 0 Turbulence - Theory and Modelling GROUP-STUDIES: Turbulence
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 informationViscoelasticity. Basic Notions & Examples. Formalism for Linear Viscoelasticity. Simple Models & Mechanical Analogies. Non-linear behavior
Viscoelasticity Basic Notions & Examples Formalism for Linear Viscoelasticity Simple Models & Mechanical Analogies Non-linear behavior Viscoelastic Behavior Generic Viscoelasticity: exhibition of both
More informationBOUNDARY LAYER ANALYSIS WITH NAVIER-STOKES EQUATION IN 2D CHANNEL FLOW
Proceedings of,, BOUNDARY LAYER ANALYSIS WITH NAVIER-STOKES EQUATION IN 2D CHANNEL FLOW Yunho Jang Department of Mechanical and Industrial Engineering University of Massachusetts Amherst, MA 01002 Email:
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 informationFluid Mechanics II Viscosity and shear stresses
Fluid Mechanics II Viscosity and shear stresses Shear stresses in a Newtonian fluid A fluid at rest can not resist shearing forces. Under the action of such forces it deforms continuously, however small
More information5 The Oldroyd-B fluid
5 The Oldroyd-B fluid Last time we started from a microscopic dumbbell with a linear entropic spring, and derived the Oldroyd-B equations: A u = u ρ + u u = σ 2 pi + η u + u 3 + u A A u u A = τ Note that
More informationInstability of a high-speed submerged elastic jet
J. Fluid Mech. (1995), od. 288, pp. 311-324 Copyright @ 1995 Cambridge University Press 311 Instability of a high-speed submerged elastic jet By J. M. RALLISON AND E. J. HINCH Department of Applied Mathematics
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 information