SlideShare a Scribd company logo

G04444146

•
•
I
IOSR-JEN

This document summarizes research on the steady flow of thin liquid films on an inclined solid surface. It presents the governing equations for this system using the Navier-Stokes equations and continuity equation. It describes the boundary conditions and introduces dimensionless variables. The solution of the governing equations for the film thickness is obtained numerically.

Related slideshows

CFD
CFDCFD
CFD
 
Exact Solutions for MHD Flow of a Viscoelastic Fluid with the Fractional Bur...
Exact Solutions for MHD Flow of a Viscoelastic Fluid with the Fractional Bur...Exact Solutions for MHD Flow of a Viscoelastic Fluid with the Fractional Bur...
Exact Solutions for MHD Flow of a Viscoelastic Fluid with the Fractional Bur...
 
www.ijerd.com
www.ijerd.comwww.ijerd.com
www.ijerd.com
 
1 of 6
Download to read offline
IOSR Journal of Engineering (IOSRJEN) www.iosrjen.org ISSN (e): 2250-3021, ISSN (p): 2278-8719 Vol. 04, Issue 04 (April. 2014), ||V4|| PP 41-46 International organization of Scientific Research 41 | P a g e Steady Flow of Thin Liquid Films on an Inclined Solid surface Joseph G. Abdulahad 1 , Salih A. Derbaz2 1,2 Department of Mathematics, Faculty of Science, University of Zakho,Duhok, Kurdistan Region, Iraq Abstract: - In this paper, we consider the thinning process of an inclined thin liquid film over a solid boundary with an inclination angle to the horizontal in gravity driven flow. Throughout this work, we assumed that the fluid thickness is constant far behind the front and we neglect the thickness of the film at the beginning of the motion. The differential equation of the film thickness is obtained analytically and the solution of equation that represents the film thickness is obtained numerically by using Rung-Kutta method. Keywords: - Thin Liquid Films , Navier-Stokes equations , continuity equation . I. INTRODUCTION We present here some of the theoretical aspects of the instability development in an inclined thin liquid films on a solid surface in two dimensional coordinate system with an inclination angle  to the horizontal . There are different types of phenomena that can occur, such as drainage, details of rupture, non–Newtonian surface properties in moving contract lines in thin liquid films [1]. These phenomena can help to describe the physical processes that occur in our real world. [2] have studied the case of contact line instabilities of thin liquid films but with constant flux configuration and also they considered some global models of a moving contact lines. [3] studied the thin liquid films flowing down the inverted substrate in three dimensional flow. [4] investigated the dynamics of an inclined thin liquid films of variable thickness in steady and unsteady cases and when the film is stationary and uniform. [5] considered the stability of thin liquid films and sessile droplets. The stability of the contact line of thin fluid film flows with constant flux configuration is considered by [6]. [7] considered the spreading of thin liquid films with small surface tension in the case when the flow is unsteady. [8] have studied the drainage of thin liquid films on an inclined solid surfacefor unsteadynflow by using similarity solution. In this paper we investigate the drainage of the inclined thin liquid films. The solution of the governing equations of the liquid film thickness is obtained numerically. II. GOVERNING EQUATIONS: Let  wuq , denotes the fluid velocity, where u and w are the velocity components in x and z directions respectively. Let  txhz , be the equation of the inclined thin liquid films as shown in Figure (1) and the flow is in x direction. The continuity equation is given by: 0      z w x u (1) and from the incompressibility condition, we have x u zw    (2) and this insures that x u   is a function of x only. The Navier-Stokes equations in x and z directions respectively for an inclined thin liquid film are given by:  sin2 2 2 2 g z u x u x P z u w x u u                          (3) And  cos2 2 2 2 g z w x w z P z w w x w u                          (4) where  ,  and P are the density, viscosity of fluid and P the pressure. In lubrication theory the inertia terms can be neglected and the Navier-Stokes equations (3) and (4) become  sin2 2 2 2 g z u x u x P               (5)
Steady Flow of Thin Liquid Films on an Inclined Solid surface International organization of Scientific Research 42 | P a g e and  cos2 2 2 2 g z w x w z P               (6) III. BOUNDARY CONDITIONS The boundary conditions to be imposed on the bounding surfaces are as follows: III.1 The no-slip condition gives at 0z , 0u (7) III. The shear stress condition on the surface vanishes, that is: at hz  , 0   z u (8) Furthermore the normal stress condition on the surface of the film is given by kw  (9) Where  is the surface tension and k is the curvature of the surface which is given by 2/32 2 2 ))(1((/ x h x h k       (10) The pressure is related to the normal stress by the formula s z w wP )(2     (11) Now from the continuity equation (1), we have    s z w )( s x u )(    , (12) and thus equation (11) is then become s x u wP )(2     (13) since 2 )( x h   is so small for thin liquid films, thus equation (10) is then reduces to give 2 2 x h k    , and equation (13) is then become s x u x h P )(22 2        (14) Differentiate equation (14), we get s x u x h x P )(2 2 2 3 3          . (15) It is to be noted that x P   is a function of x only on the surface of the liquid film z = h ( x ). Now from the conservation of mass and since the free surface is a stream line, the derivative following the motion ( the material or the substantial derivative) Dt DF must be vanished on z = h ( x ) and thus, for steady flow we have 0      z F w x F u (16) Where F(x,z) = z–h(x) (17) from equations (16) and (17), we have 0   w x h u (18) now substitute equation (2) into equation (18), we obtain 0      x u z x h u (19) on the surface of the liquid film z = h( x ), equation (19) then gives 0      x u h x h u (20)
Steady Flow of Thin Liquid Films on an Inclined Solid surface International organization of Scientific Research 43 | P a g e IV. Dimensional analysis: We now introduce the following non-dimensional variables as follow: xLx  , zHz  , uUu  , wUw  , P L MU P  , (21) hHh  and 1 L H  where L, H and U are the characteristics. By using the dimensionless variables the continuity equation (1) then gives 0      z w x u (22) Also equations (14) in non-dimensional variables has the form s x u x h Ca P )(22 2        (23) Where  U Ca  is the capillary number. Furthermore the Navier-Stokes equations (5) and (6) in non-dimensional variables are respectively given by   sin 2 2 2 22 2 Ca B z u x u x P          (24) and  cos2 2 2 2 2 Ca B z w x w z P          (25) Where   2 gL B  is the bond number. Also equation (20) has the following non-dimensional form 0      x u h x h u or 0)(    hu x or Re)( hu (26) Where  UL Re is the Reynolds number and     is the kinematic viscosity. From (23) and on the surface of the film, we have s x u x h Cax P )(2 2 2 3 3          (27) From equations (24) and (27), we get s x u x h CaCa B z u x u )(2sin 2 2 2 3 3 2 2 22 2               (28) Or 0sin 2 )(3 2 2 22 2 3 3             Ca B z u x u x h Ca s (29) Integrating equation (29) with respect to z , we get )( 2 )sin)(3( 122 2 3 3 xA z u z Ca B x u x h Ca s             (30)
Steady Flow of Thin Liquid Films on an Inclined Solid surface International organization of Scientific Research 44 | P a g e The boundary conditions (7) and (8) in non- dimensional variables have the form: At 0,0  uz (31) At 0,)(     z u xhz (32) By using the boundary condition (32), equation (30) reduces to give )()sin)(3( 12 2 3 3 xAh Ca B x u x h Ca s         (33) And thus equation (30), gives 0 2 )()sin)(3( 22 2 3 3          z u hz Ca B x u x h Ca s    (34) Integrating equation (34), we get )( 2 ) 2 ()sin)(3( 22 2 2 2 3 3 xAuzh z Ca B x u x h Ca s          (35) From the boundary condition (31) and equation (35), we have   02 xA , And thus equation (35), then gives 0 2 ) 2 ()sin)(3( 2 2 2 2 3 3       uzh z Ca B x u x h Ca s    (36) Since the first term in equation (36) is very small so we can neglected and thus equation (36) reduces to give 0 2 ) 2 (]sin)(3[ 2 2 2 2    uzh z Ca B x u s   (37) From equation (26), we have ) 2 (Re)( 3 '''2 2 2 2 h hhhh x u s      (38) Substituting equation (38) into equation (37), we get 0 2 ) 2 (]sin)2(Re3[ 3 2 2 3'''2 2  huzh z h Ca B hhhh   (39) From (26) and equation (39), we obtain 0 Re2 ) 2 (]sin)2(Re3[ 2 2 2 3'''2 2  hzh z h Ca B hhhh    (40) On the surface of the film z = h(x), equation (40) then gives 0 Re2 ) 2 (]sin)2(Re3[ 2 2 2 3'''2 2  h h h Ca B hhhh    Or 0 Re4 sin)2(Re3 2 3'''2 2     h Ca B hhhh (41) Equation (41) represents the non-linear ordinary differential equation for the thickness of the inclined thin liquid film on a solid boundary in ),( hx plane for steady flow. Also from (39), we have ) 2 ()sin 2 Re3 2 Re3 ( 22 2 '2''2 2 zh z Ca B h h h h u    (42) And the average velocity over the film thickness is given by  h av dzu h u 0 1 (43) And thus from equations (42) and (43) we have 2 2 '2'' 2 sin 6 Re 2 Re 2 h Ca B hhhuav      (44) Some solution curves of equation (41) are presented in figure (2)and shows that the thickness of liquid film for increases as the angle of inclination to the horizontal increases, moreover figure (3) show the thickness of some liquid films namely linseed oil, olive oil and glycerinfor =1.57.
Steady Flow of Thin Liquid Films on an Inclined Solid surface International organization of Scientific Research 45 | P a g e V. FIGURES AND TABLES Fig. (1): Sketch of the flow in two–dimensional geometry, where the capillary ridge is just behind the flow front. (water) (Glyserin) (Carbon Tetrachloride) Fig.(2): Thickness of thin liquid films for different inclination angles = 0.174 = 0.52 = 1.04 = 1.57 Fig.(3): Thickness of thin liquid films for linseed oil,olive oil and glycerin for = 1.57. linseed oil olive oil glycerin VI. CONCLUSION The dynamics of a free surface liquid films flowing over an inclined solid boundary with an inclination angle to the horizontal is very useful in industrial coating and spinning processes. The non- linear differential equation that governs such flow for steady case in non- dimensional form is obtained and the solution curves for different inclination angles shows that the thickness of the thin liquid films increases as the inclination angle increases..
Steady Flow of Thin Liquid Films on an Inclined Solid surface International organization of Scientific Research 46 | P a g e REFERENCES [1] Bertozzi, A.L. (1998). The mathematics of moving contact lines in thin liquid film, Amer. Math. Soc. Pp. 689–697, 45. [2] Diez, J., Kondic L. and Bertozzi, A.L. (2001). Global models for moving contact lines, Phys. Rev., pp. 011208, 63. [3] Lin. T. S., Kondic, L. and Filippove, A.(2012). Thin films flowing down inverted Substrated three dimensional flow, Physics of fluids, Vol.24, pp. 1-18. [4] Karmina, K. A., (2013). Fluid flow and stability analysis in certain thin liquid films. M.Sc. thesis, University of Zakho. [5] Drofler, F., Rauscher M. and Dietrich S.(2013). Stability of thin liquid films and sessile droplets under confinement, Eur. Phys. J. E.Vol.20, pp. 1-14 [6] Kondic L (2003). Instabilities in gravity driven flow of thin fluid films, Siam Review, vol. 45, No. 1, pp. 95–115. [7] Moriarty J.A., and Schwartz L.W.(1991). Unsteady spreading of thin liquid films with small surface tension, American institute of Physics, vol3, No.5, pp. 733-742. [8] Joseph G. A., and Salih A. D., (2014) .The drainage of thin liquid films on an inclined solid surface .IOSR- JM, vol 10, pp 73-80.

More Related Content

G04444146

  • 1. IOSR Journal of Engineering (IOSRJEN) www.iosrjen.org ISSN (e): 2250-3021, ISSN (p): 2278-8719 Vol. 04, Issue 04 (April. 2014), ||V4|| PP 41-46 International organization of Scientific Research 41 | P a g e Steady Flow of Thin Liquid Films on an Inclined Solid surface Joseph G. Abdulahad 1 , Salih A. Derbaz2 1,2 Department of Mathematics, Faculty of Science, University of Zakho,Duhok, Kurdistan Region, Iraq Abstract: - In this paper, we consider the thinning process of an inclined thin liquid film over a solid boundary with an inclination angle to the horizontal in gravity driven flow. Throughout this work, we assumed that the fluid thickness is constant far behind the front and we neglect the thickness of the film at the beginning of the motion. The differential equation of the film thickness is obtained analytically and the solution of equation that represents the film thickness is obtained numerically by using Rung-Kutta method. Keywords: - Thin Liquid Films , Navier-Stokes equations , continuity equation . I. INTRODUCTION We present here some of the theoretical aspects of the instability development in an inclined thin liquid films on a solid surface in two dimensional coordinate system with an inclination angle  to the horizontal . There are different types of phenomena that can occur, such as drainage, details of rupture, non–Newtonian surface properties in moving contract lines in thin liquid films [1]. These phenomena can help to describe the physical processes that occur in our real world. [2] have studied the case of contact line instabilities of thin liquid films but with constant flux configuration and also they considered some global models of a moving contact lines. [3] studied the thin liquid films flowing down the inverted substrate in three dimensional flow. [4] investigated the dynamics of an inclined thin liquid films of variable thickness in steady and unsteady cases and when the film is stationary and uniform. [5] considered the stability of thin liquid films and sessile droplets. The stability of the contact line of thin fluid film flows with constant flux configuration is considered by [6]. [7] considered the spreading of thin liquid films with small surface tension in the case when the flow is unsteady. [8] have studied the drainage of thin liquid films on an inclined solid surfacefor unsteadynflow by using similarity solution. In this paper we investigate the drainage of the inclined thin liquid films. The solution of the governing equations of the liquid film thickness is obtained numerically. II. GOVERNING EQUATIONS: Let  wuq , denotes the fluid velocity, where u and w are the velocity components in x and z directions respectively. Let  txhz , be the equation of the inclined thin liquid films as shown in Figure (1) and the flow is in x direction. The continuity equation is given by: 0      z w x u (1) and from the incompressibility condition, we have x u zw    (2) and this insures that x u   is a function of x only. The Navier-Stokes equations in x and z directions respectively for an inclined thin liquid film are given by:  sin2 2 2 2 g z u x u x P z u w x u u                          (3) And  cos2 2 2 2 g z w x w z P z w w x w u                          (4) where  ,  and P are the density, viscosity of fluid and P the pressure. In lubrication theory the inertia terms can be neglected and the Navier-Stokes equations (3) and (4) become  sin2 2 2 2 g z u x u x P               (5)
  • 2. Steady Flow of Thin Liquid Films on an Inclined Solid surface International organization of Scientific Research 42 | P a g e and  cos2 2 2 2 g z w x w z P               (6) III. BOUNDARY CONDITIONS The boundary conditions to be imposed on the bounding surfaces are as follows: III.1 The no-slip condition gives at 0z , 0u (7) III. The shear stress condition on the surface vanishes, that is: at hz  , 0   z u (8) Furthermore the normal stress condition on the surface of the film is given by kw  (9) Where  is the surface tension and k is the curvature of the surface which is given by 2/32 2 2 ))(1((/ x h x h k       (10) The pressure is related to the normal stress by the formula s z w wP )(2     (11) Now from the continuity equation (1), we have    s z w )( s x u )(    , (12) and thus equation (11) is then become s x u wP )(2     (13) since 2 )( x h   is so small for thin liquid films, thus equation (10) is then reduces to give 2 2 x h k    , and equation (13) is then become s x u x h P )(22 2        (14) Differentiate equation (14), we get s x u x h x P )(2 2 2 3 3          . (15) It is to be noted that x P   is a function of x only on the surface of the liquid film z = h ( x ). Now from the conservation of mass and since the free surface is a stream line, the derivative following the motion ( the material or the substantial derivative) Dt DF must be vanished on z = h ( x ) and thus, for steady flow we have 0      z F w x F u (16) Where F(x,z) = z–h(x) (17) from equations (16) and (17), we have 0   w x h u (18) now substitute equation (2) into equation (18), we obtain 0      x u z x h u (19) on the surface of the liquid film z = h( x ), equation (19) then gives 0      x u h x h u (20)
  • 3. Steady Flow of Thin Liquid Films on an Inclined Solid surface International organization of Scientific Research 43 | P a g e IV. Dimensional analysis: We now introduce the following non-dimensional variables as follow: xLx  , zHz  , uUu  , wUw  , P L MU P  , (21) hHh  and 1 L H  where L, H and U are the characteristics. By using the dimensionless variables the continuity equation (1) then gives 0      z w x u (22) Also equations (14) in non-dimensional variables has the form s x u x h Ca P )(22 2        (23) Where  U Ca  is the capillary number. Furthermore the Navier-Stokes equations (5) and (6) in non-dimensional variables are respectively given by   sin 2 2 2 22 2 Ca B z u x u x P          (24) and  cos2 2 2 2 2 Ca B z w x w z P          (25) Where   2 gL B  is the bond number. Also equation (20) has the following non-dimensional form 0      x u h x h u or 0)(    hu x or Re)( hu (26) Where  UL Re is the Reynolds number and     is the kinematic viscosity. From (23) and on the surface of the film, we have s x u x h Cax P )(2 2 2 3 3          (27) From equations (24) and (27), we get s x u x h CaCa B z u x u )(2sin 2 2 2 3 3 2 2 22 2               (28) Or 0sin 2 )(3 2 2 22 2 3 3             Ca B z u x u x h Ca s (29) Integrating equation (29) with respect to z , we get )( 2 )sin)(3( 122 2 3 3 xA z u z Ca B x u x h Ca s             (30)
  • 4. Steady Flow of Thin Liquid Films on an Inclined Solid surface International organization of Scientific Research 44 | P a g e The boundary conditions (7) and (8) in non- dimensional variables have the form: At 0,0  uz (31) At 0,)(     z u xhz (32) By using the boundary condition (32), equation (30) reduces to give )()sin)(3( 12 2 3 3 xAh Ca B x u x h Ca s         (33) And thus equation (30), gives 0 2 )()sin)(3( 22 2 3 3          z u hz Ca B x u x h Ca s    (34) Integrating equation (34), we get )( 2 ) 2 ()sin)(3( 22 2 2 2 3 3 xAuzh z Ca B x u x h Ca s          (35) From the boundary condition (31) and equation (35), we have   02 xA , And thus equation (35), then gives 0 2 ) 2 ()sin)(3( 2 2 2 2 3 3       uzh z Ca B x u x h Ca s    (36) Since the first term in equation (36) is very small so we can neglected and thus equation (36) reduces to give 0 2 ) 2 (]sin)(3[ 2 2 2 2    uzh z Ca B x u s   (37) From equation (26), we have ) 2 (Re)( 3 '''2 2 2 2 h hhhh x u s      (38) Substituting equation (38) into equation (37), we get 0 2 ) 2 (]sin)2(Re3[ 3 2 2 3'''2 2  huzh z h Ca B hhhh   (39) From (26) and equation (39), we obtain 0 Re2 ) 2 (]sin)2(Re3[ 2 2 2 3'''2 2  hzh z h Ca B hhhh    (40) On the surface of the film z = h(x), equation (40) then gives 0 Re2 ) 2 (]sin)2(Re3[ 2 2 2 3'''2 2  h h h Ca B hhhh    Or 0 Re4 sin)2(Re3 2 3'''2 2     h Ca B hhhh (41) Equation (41) represents the non-linear ordinary differential equation for the thickness of the inclined thin liquid film on a solid boundary in ),( hx plane for steady flow. Also from (39), we have ) 2 ()sin 2 Re3 2 Re3 ( 22 2 '2''2 2 zh z Ca B h h h h u    (42) And the average velocity over the film thickness is given by  h av dzu h u 0 1 (43) And thus from equations (42) and (43) we have 2 2 '2'' 2 sin 6 Re 2 Re 2 h Ca B hhhuav      (44) Some solution curves of equation (41) are presented in figure (2)and shows that the thickness of liquid film for increases as the angle of inclination to the horizontal increases, moreover figure (3) show the thickness of some liquid films namely linseed oil, olive oil and glycerinfor =1.57.
  • 5. Steady Flow of Thin Liquid Films on an Inclined Solid surface International organization of Scientific Research 45 | P a g e V. FIGURES AND TABLES Fig. (1): Sketch of the flow in two–dimensional geometry, where the capillary ridge is just behind the flow front. (water) (Glyserin) (Carbon Tetrachloride) Fig.(2): Thickness of thin liquid films for different inclination angles = 0.174 = 0.52 = 1.04 = 1.57 Fig.(3): Thickness of thin liquid films for linseed oil,olive oil and glycerin for = 1.57. linseed oil olive oil glycerin VI. CONCLUSION The dynamics of a free surface liquid films flowing over an inclined solid boundary with an inclination angle to the horizontal is very useful in industrial coating and spinning processes. The non- linear differential equation that governs such flow for steady case in non- dimensional form is obtained and the solution curves for different inclination angles shows that the thickness of the thin liquid films increases as the inclination angle increases..
  • 6. Steady Flow of Thin Liquid Films on an Inclined Solid surface International organization of Scientific Research 46 | P a g e REFERENCES [1] Bertozzi, A.L. (1998). The mathematics of moving contact lines in thin liquid film, Amer. Math. Soc. Pp. 689–697, 45. [2] Diez, J., Kondic L. and Bertozzi, A.L. (2001). Global models for moving contact lines, Phys. Rev., pp. 011208, 63. [3] Lin. T. S., Kondic, L. and Filippove, A.(2012). Thin films flowing down inverted Substrated three dimensional flow, Physics of fluids, Vol.24, pp. 1-18. [4] Karmina, K. A., (2013). Fluid flow and stability analysis in certain thin liquid films. M.Sc. thesis, University of Zakho. [5] Drofler, F., Rauscher M. and Dietrich S.(2013). Stability of thin liquid films and sessile droplets under confinement, Eur. Phys. J. E.Vol.20, pp. 1-14 [6] Kondic L (2003). Instabilities in gravity driven flow of thin fluid films, Siam Review, vol. 45, No. 1, pp. 95–115. [7] Moriarty J.A., and Schwartz L.W.(1991). Unsteady spreading of thin liquid films with small surface tension, American institute of Physics, vol3, No.5, pp. 733-742. [8] Joseph G. A., and Salih A. D., (2014) .The drainage of thin liquid films on an inclined solid surface .IOSR- JM, vol 10, pp 73-80.