CONSERVATION LAWS AND CONSERVED QUANTITIES FOR LAMINAR RADIAL JETS WITH SWIRL

Size: px
Start display at page:

Download "CONSERVATION LAWS AND CONSERVED QUANTITIES FOR LAMINAR RADIAL JETS WITH SWIRL"

Transcription

1 Mathematical and Computational Applications,Vol. 15, No. 4, pp , 21. c Association for Scientific Research CONSERVATION LAWS AND CONSERVED QUANTITIES FOR LAMINAR RADIAL JETS WITH SWIRL R. Naz 1, I. Naeem 2 and F.M. Mahomed 3 1 Centre for Mathematics and Statistical Sciences Lahore School of Economics, Lahore 532, Pakistan. 2 Department of Mathematics, School of Science and Engineering, LUMS, Lahore Cantt 54792, Pakistan. 3 School of Computational and Applied Mathematics, Centre for Differential Equations, Continuum Mechanics and Applications, Universit of the Witwatersrand, Wits 25, South Africa. rehananaz qau@ahoo.com, inqau@ahoo.com, Fazal.Mahomed@wits.ac.za Abstract- Conserved quantities pla a central role in the solution of jet flow problems. A sstematic wa of deriving conserved quantities for the radial jets with swirl is presented. The multiplier approach is used to derive the conservation laws for the sstem of three boundar laer equations for the velocit components and for the sstem of two partial differential equations for the stream function. When the swirl is zero or at a large distance from the orifice (at infinit), the boundar laer equations for the radial jets with swirl reduce to those of the purel radial jets. The conserved quantities for the radial liquid, free and wall jets with swirl are derived b integrating the conservation laws across the jets. Kewords- Radial jets with swirl, conservation laws, conserved quantit. 1. INTRODUCTION It is of great interest to stud the three dimensional boundar laer flows, for eample the radial jets with swirl due to their etensive range of applications. The concept of two-dimensional boundar laer flow of the purel radial liquid jet was introduced b Watson [1]. The two-dimensional boundar laer flow of the purel radial free jet was first considered b Schwarz [2]. Glauert in [3] studied the two-dimensional boundar laer flow of the purel radial wall jet. Later, Rile in [4] also discussed the two-dimensional boundar laer flow of purel radial liquid, free and wall jets. In all these jet flow problems, the boundar conditions are homogeneous and therefore the unknown eponent in the similarit solution cannot be obtained from the boundar conditions. B integrating Prandtl s momentum boundar laer equation across the jet and using the boundar conditions and the

2 R. Naz, I. Naeem and F. M. Mahomed 743 continuit equation, a condition, known as the conserved quantit, was derived. Rile [5] etended the problem of the purel radial jets to the radial jets with swirl. Later, O Nan and Schwarz [6] also studied the problem of the radial free jet with swirl. The governing equations consist of a continuit equation and two momentum equations in the radial and swirl directions respectivel. This problem consists of homogeneous boundar conditions and thus also requires conserved quantities. The conserved quantities for the radial jets with swirl were established using the same idea as was used in [1, 2, 3, 4] to construct the conserved quantities for the purel radial jets. B integrating the momentum equation in the swirl direction and using the boundar conditions and the continuit equation, one conserved quantit for the radial free jet with swirl and one for the radial wall jet with swirl were derived. The second conserved quantit for the radial free and wall jets with swirl were established b integrating the momentum equation in the radial direction and requiring that swirl is zero at =. The onl conserved quantit for the radial liquid jet with swirl [5] was obtained from a phsical argument based on conservation of mass flu in an incompressible fluid and is the same as for the purel radial liquid jet. The conserved quantit for each problem has in general phsical significance. Rile [5] established the conserved quantities for the radial liquid, free and wall jets with swirl. The onl conserved quantit for the radial liquid jet with swirl is total mass flu per radian which is constant. One of the conserved quantities for the radial free jet with swirl is the radial flu of angular momentum about =. The second conserved quantit for the radial free jet with swirl is the radial flu of momentum and is invariant onl when the swirl is zero or at large distances from the orifice. For a radial wall jet with swirl one conserved quantit is the flu of eterior momentum flu and is constant onl when the swirl is zero or at large distances from the orifice. There is no obvious phsical interpretation of the other conserved quantit for the radial wall jet with swirl. In this paper, we construct the conserved quantities for the radial liquid, free and wall jets with swirl b the method recentl introduced b Naz et al [7]. We use the multiplier approach introduced b Steudel [8] (see also [9, 1]) to derive the conservation laws for the sstem of three partial differential equations for the velocit components as well as for the sstem of two partial differential equations for the stream function. The radial jets with swirl behave as purel radial jets when the swirl is zero or at large distances from the orifice as becomes large the swirl component will tend to zero more rapidl than the radial velocit component. Then the boundar laer equations for radial jets with swirl reduce to that of purel radial jets. The conservation laws for the purel radial jets were derived in [7] with the

3 744 Conservation Laws for Laminar Radial Jets with Swirl help of the multiplier approach. The conserved quantities for the radial, liquid, free and wall jets with swirl are derived with the help of conservation laws. The radial jets (liquid, free and wall) with swirl satisf the same partial differential equations but the boundar conditions for each jet are different. The derivation of conservation laws depends onl on the partial differential equations but the derivation of the phsical conserved quantities depends also on the boundar conditions. The boundar conditions therefore determine which conserved vector is associated with which jet. 2. CONSERVATION LAWS FOR RADIAL LAMINAR JET We consider the radial liquid jet, free jet and wall jet with swirl. The radial liquid jet with swirl is formed when a circular jet of liquid, with an aiall smmetrical swirling component of velocit, strikes a plane boundar normall and spreads over it. The radial free jet with swirl is formed when two circular jets, each with an aiall-smmetrical swirling component of velocit, impinge on one another. When a circular jet, having aiall smmetric swirling component of velocit, strikes a plane surface normall which is surrounded b the same fluid as the jet and spreads out radiall over it then a radial wall jet with swirl is formed. All these jets were described b Rile [5]. The fluid in the jets is viscous and incompressible. In the free jet and wall jet the surrounding fluid consists of the same fluid as in the jet but for the liquid jet the surrounding fluid is a gas. The clindrical polar coordinates (, θ, ) are used with as radial coordinate and = is the ais of smmetr. All fluid variables are independent of θ. For the liquid jet the solid boundar is at = and the free surface is at = ϕ(). The wall jet has a solid boundar at = and the free jet is smmetrical about =. For all three jets an aiall smmetrical swirling component of velocit is considered. The boundar laer equations governing the flow in a radial laminar jet with swirl, in the absence of pressure gradient, are uu + vu w2 = νu, (1) uw + vw + uw = νw, (2) (u) + (v) =, (3) where u(, ) and v(, ) are the velocit components in and directions respectivel, w(, ) is an aiall-smmetrical swirling component of velocit and ν is the

4 R. Naz, I. Naeem and F. M. Mahomed 745 kinematic viscosit of the fluid. Equations (1) and (2) are the momentum equations in the radial and θ directions respectivel and (3) is the continuit equation. The boundar-laer equations (1)-(3) have been derived with the assumption that the swirl velocit w and radial velocit u are of same order of magnitude whereas the transverse velocit v is of lower order (see [5]). The stream function ψ(, ) satisfies [11] u = ψ, v = ψ. (4) The radial free jet with swirl is smmetrical about the plane =, therefore v(, ) =. For the liquid jet and the wall jet there is no suction or blowing of fluid at the solid boundar = and thus v(, ) =. Therefore ψ (, ) = for all three jets and thus ψ(, ) is a constant. We choose the constant as zero and hence ψ(, ) =. (5) Let O be an reference point on the line = in the aial plane θ = θ o, then the stream function at P (, θ o, ) is (see [11]) ψ(, ) = P O [u(, )d v(, )d]. (6) Integrating (6) with respect to from = to = along a straight line with kept fied ields ψ(, ) = u(, )d. (7) Now, ψ(, ) is finite onl if we assume u(, ) sufficientl rapidl as. and Equation (3) is identicall satisfied and equations (1) and (2) ield 1 ψ ψ 1 2 ψ2 1 ψ ψ νψ w 2 =, (8) 1 ψ w 1 ψ w ψ w νw =. (9) We will derive the conservation laws for the sstem (1)-(3) and also for the sstem (8)-(9) for the stream function, using the multiplier approach [8, 9].

5 746 Conservation Laws for Laminar Radial Jets with Swirl 2.1. Conservation laws for sstem of equations for velocit components Consider the multipliers of the form Λ 1 = Λ 1 (,, u, v, w), Λ 2 = Λ 2 (,, u, v, w)andλ 3 = Λ 3 (,, u, v, w) for the sstem of equations (1)-(3), such that ( ) ( Λ 1 uu + vu νu w2 + Λ 2 uw + vw + uw ) νw +Λ 3 (u + u + v ) = D T 1 + D T 2, (1) for all functions u(, ), v(, ) and w(, ) and not onl for the solutions of sstem. In (1), D and D are total derivative operators defined b D = + u u + v v + w w + u + v + w u v w +u + v + w +, (11) u v w D = + u u + v v + w w + u + v + w u v w +u + v + w +. (12) u v w The determining equations for the multipliers Λ 1, Λ 2 and Λ 3 are [ E u Λ 1 (uu + vu νu w2 ) + Λ 2(uw + vw + uw νw ) ] +Λ 3 (u + u + v ) =, (13) [ E v Λ 1 (uu + vu νu w2 ) + Λ 2(uw + vw + uw νw ) ] +Λ 3 (u + u + v ) =, (14) [ E w Λ 1 (uu + vu νu w2 ) + Λ 2(uw + vw + uw νw ) ] +Λ 3 (u + u + v ) =, (15) where E u, E v and E w are the standard Euler operators which annihilate divergence epressions: E u = u D D + D 2 + D D + D 2, (16) u u u u u

6 R. Naz, I. Naeem and F. M. Mahomed 747 and E v = v D D + D 2 + D D + D 2, (17) v v v v v E w = w D D + D 2 + D D + D 2. (18) w w w w w The epansion of equations (13)-(15) ields Λ 1u (uu + vu νu w2 ) + Λ 2u(uw + vw + uw νw ) + Λ 3u (u + u + v ) +Λ 1 u + Λ 2 [w + w ] + Λ 3 D [uλ 1 + Λ 3 ] D [vλ 1 ] νd 2 (Λ 1 ) =, (19) Λ 1v (uu + vu νu w2 ) + Λ 2v(uw + vw + uw νw ) + Λ 3v (u + u + v ) +Λ 1 u + Λ 2 w D (Λ 3 ) =, (2) Λ 1w (uu + vu νu w2 ) + Λ 2w(uw + vw + uw νw ) + Λ 3w (u + u + v ) 2w Λ 1 + u Λ 2 D [uλ 2 ] D [vλ 2 ] νd 2 (Λ 2 ) =. (21) Equations (19)-(21) must be satisfied for all functions u(, ), v(, ) and w(, ). After epansion the highest order derivative terms in sstem (19)-(21) are secondorder derivative terms. Therefore, equating the coefficients of u, v and w to zero, we obtain Λ 1u =, Λ 1v =, Λ 1w =, Λ 2u =, Λ 2v =, Λ 2w =, (22) which in turn results in Λ 1 = A(, ), Λ 2 = B(, ), Λ 3 = C(,, u, v, w). (23) Equations (19)-(21), with the help of (23), reduce to v C v +v (A C u )+w (C w B) uc u Bw +ua +C +va +νa =, (24) u C v + u (A C u ) w (C w B) + uc v C =, (25) (u + v )(C w B) + uc w 2w A + u B ub vb νb =. (26) Equating the first-order derivative terms u, u and w in (25) to zero, we obtain C v =, A C u =, C w B =, (27)

7 748 Conservation Laws for Laminar Radial Jets with Swirl which ields C = u A(, ) + w B(, ) + D(, ). (28) Equations (24)-(26) become ( 2u A A ) + va + w ( B 2B ) + D + νa =, (29) ua + wb + D =, (3) ( u B 2B ) + vb + 2w A + νb =. (31) The separation of equations (29)-(31) with respect to u, v and w finall gives Λ 1 =, Λ 2 = c 2 2, Λ 3 = c 1 + c 2 w, (32) where c 1 and c 2 are constants. It is not difficult to construct the conserved vectors b elementar manipulations with the help of the multipliers (32). Equations (1) and (32) give [uu +vu νu w2 ]+c 2 2 [uw +vw + uw νw ]+(c 1 +c 2 w)[u+u +v ] = D [c 1 u + c 2 2 uw] + D [c 1 v + c 2 2 (vw νw )], (33) for arbitrar functions u(, ), v(, ) and w(, ) and therefore when u(, ), v(, ) and w(, ) are the solutions of the sstem of differential equations (1)-(3), then Thus D [c 1 u + c 2 2 uw] + D [c 1 v + c 2 2 (vw νw )] =. (34) T 1 = u, T 2 = v, (35) T 1 = 2 uw, T 2 = 2 (vw νw ), (36) are the conserved vectors for the sstem of differential equations (1)-(3) corresponding to multipliers of the form Λ 1 = Λ 1 (,, u, v, w) and Λ 2 = Λ 2 (,, u, v, w) Conservation Laws for the stream function formulation Consider multipliers of the form Λ 1 = Λ 1 (,, ψ, w) and Λ 2 = Λ 2 (,, ψ, w) in order to derive the conservation laws for the sstem (8)-(9). The determining equations for the multipliers Λ 1 and Λ 2 are E ψ [Λ 1 ( 1 ψ ψ 1 2 ψ2 1 ψ ψ νψ w 2 )

8 R. Naz, I. Naeem and F. M. Mahomed 749 +Λ 2 ( 1 ψ w 1 ψ w + 1 ] ψ w νw 2 ) =, (37) and where E w [Λ 1 ( 1 ψ ψ 1 2 ψ2 1 ψ ψ νψ w 2 ) +Λ 2 ( 1 ψ w 1 ψ w + 1 ] ψ w νw 2 ) =, (38) E ψ = ψ D D + D 2 + D D + D 2, (39) ψ ψ ψ ψ ψ E w = w D D + D 2 + D D + D 2, (4) w w w w w are the standard Euler operators. In (39) and (4), D and D are the total derivative operators and are defined b D = + ψ ψ + w w + ψ + w + ψ + w +, (41) ψ w ψ w D = + ψ ψ + w w + ψ + w + ψ + w +. (42) ψ w ψ w Epansion of equations (37) and (38) ields Λ 1ψ [ 1 ψ ψ 1 2 ψ2 1 ψ ψ νψ w 2 ] + Λ 2ψ [ 1 ψ w 1 ψ w ψ w νw ] + D [ 1 Λ 1ψ + 1 Λ 2w ] D [Λ 1 ( 1 ψ 2 2 ψ ) +Λ 2 ( w + w 2 )] + D D [ 1 Λ 1ψ ] D 2 [ 1 Λ 1ψ ] + νd 3 [Λ 1 ] =, (43) Λ 1w [ 1 ψ ψ 1 2 ψ2 1 ψ ψ νψ w 2 ] + Λ 2w [ 1 ψ w 1 ψ w ψ w νw ] 2wΛ 1 + Λ 2 ψ 2 D [ Λ 2 ψ ] + D [ Λ 2 ψ ] νd[λ 2 2 ] =. (44) Equating the coefficients of ψ and w in (44) to zero, we obtain Λ 1w =, Λ 2w =. (45) Now, since Λ 1 and Λ 2 are not functions of w, therefore separating the rest of the terms in (44) (in epanded form) with respect to w ields Λ 1 =. (46)

9 75 Conservation Laws for Laminar Radial Jets with Swirl Equations (43) and (44) reduce to Λ 2ψ [ 1 ψ w 1 ψ w ψ w νw ] 2w Λ w [Λ 2 + ψ Λ 2ψ ] and ( w + w 2 )[Λ 2 + ψ Λ 2ψ ] =, (47) 2ψ Λ 2 2 ψ [Λ 2 + ψ Λ 2ψ ] + ψ [Λ 2 + ψ Λ 2ψ ] ν[λ 2 + ψ Λ 2ψ + ψ Λ 2ψ + ψ Λ 2ψ + ψ 2 Λ 2ψψ ] =. (48) From equation (47), the coefficient of w (or in (48) the coefficient of ψ ) ields Λ 2ψ = and thus Λ 2 = A(, ). (49) Equations (47) and (48) become A w + w [ 2A A 2 ] + w A =, (5) 2 ψ A + ψ [ 2A 2 A ] νa =, (51) which finall ield A = c 3 2 and Λ 1 =, Λ 2 = c 3 2, (52) where c 3 is an arbitrar constant. The multipliers Λ 1 and Λ 2 for the sstem (8)-(9) satisf Λ 1 [ 1 ψ ψ 1 2 ψ2 1 ψ ψ νψ w 2 ] +Λ 2 [ 1 ψ w 1 ψ w ψ w νw ] (53) = D T 1 + D T 2, for all functions ψ(, ) and w(, ). From equations (52) and (53), we have [ 1 ψ ψ 1 2 ψ2 1 ψ ψ νψ w 2 ]+c 3 2 [ 1 ψ w 1 ψ w ψ w νw ] = D [c 3 wψ ] + D [c 3 ( wψ ν 2 w )], (54)

10 R. Naz, I. Naeem and F. M. Mahomed 751 for arbitrar ψ(, ) and w(, ). sstem (8)-(9), we have When ψ(, ) and w(, ) are solutions of the D [c 3 wψ ] + D [c 3 ( wψ ν 2 w )] =. (55) Thus T 1 = wψ, T 2 = wψ ν 2 w, (56) are the components of a conserved vector for the sstem (8)-(9). Notice that the conserved vector (56) is equivalent to the conserved vector (36). The multiplier approach gives multipliers onl for the local conserved vectors. The multipliers for the non-local conserved vectors cannot be obtained with the help of the multiplier approach as presented. Notice that if we multipl equation (8) with zero and equations (9) with 2 ψ, we have [ 1 ψ ψ 1 2 ψ2 1 ψ ψ νψ w 2 ]+ 2 ψ[ 1 ψ w 1 ψ w ψ w νw ] ] = D [ψwψ ] + D [ wψψ ν 2 ψw + ν 2 w ψ d, (57) for arbitrar ψ(, ) and w(, ). When ψ(, ) and w(, ) are solutions of the sstem (8)-(9), we have ] D [ψwψ ] + D [ wψψ ν 2 ψw + ν 2 w ψ d =. (58) Thus T 1 = ψwψ, T 2 = wψψ ν 2 ψw + ν 2 w ψ d, (59) are the components of a non-local conserved vector for the sstem (8)-(9). It ma be possible that the non-local conservation law (59) which does not arise from local multipliers can be derived through a potential sstem associated with the sstem (1), (2) and (3) or the sstem (8) and (9) Conservation laws for radial jets at large distance from the orifice The radial jets with swirl at large distance from the orifice, correspond to a radial jet with zero swirl. When, w, then the sstem (1)-(3) reduces to uu + vu = νu, (6) (u) + (v) =, (61)

11 752 Conservation Laws for Laminar Radial Jets with Swirl which is same as the boundar laer equations for purel radial jet flows. The sstem (8)-(9) for the stream function, when, reduces to a single third-order partial differential equation 1 ψ ψ 1 2 ψ2 1 ψ ψ νψ =. (62) The local conserved vectors for sstem (6)-(61) derived in [7] are T 1 = u, T 2 = v, (63) T 1 = u 2, T 2 = uv νu. (64) The conserved vectors for the third-order partial differential equation (62) (see [7]) are and T 1 = 1 ψ2, T 2 = 1 ψ ψ νψ, (65) T 1 = 1 ψψ2, T 2 = 1 ψψ ψ + ν 2 ψ2 νψψ. (66) The conserved vectors given in equations (64) and (65) are equivalent. 3. CONSERVED QUANTITIES FOR RADIAL LIQUID JETS In this section, we derive the conserved quantities for the radial liquid, free and wall jets with swirl b the technique used in [7]. We first consider the conserved vectors (35), (36) and (64) for the sstem of equations for the velocit components and derive the conserved quantities for the radial liquid and the free jets with swirl. Then we consider the stream function formulation to give an alternative derivation of the conserved quantities for the radial free jet with swirl and also derive the conserved quantities for the radial wall jet with swirl. All the conserved vectors (T 1, T 2 ) derived here depend on u(, ), v(, ), w(, ) or ψ(, ) and satisf D T 1 + D T 2 = T 1 + T 2. (67) But for a conserved vector D T 1 + D T 2 = and therefore equation (67) ields T 1 + T 2 =. (68)

12 R. Naz, I. Naeem and F. M. Mahomed 753 The conserved quantities for the radial liquid, free and wall jets with swirl are obtained b integrating equation (68) for the corresponding conservation law across the jet and imposing the boundar conditions Conserved quantit for the radial liquid jet with swirl The conserved quantit for the radial liquid jet with swirl is derived with the help of conserved vector (35). The derivation for conserved quantit for radial liquid jet with swirl is same as in [7] for the purel radial jet flow. The boundar conditions for the radial liquid jet with swirl are = : u =, v =, w =, (69) = ϕ() : u =, w =. (7) The velocit component v(, ϕ()) is where v(, ϕ()) = D [ϕ()] = u(, ϕ())dϕ() Dt d, (71) D Dt = + u(, ) t + v(, ) (72) is the material time derivative. The conserved quantit for the radial liquid jet with swirl is obtained b integrating (68) with respect to from = to = ϕ() keeping as fied. For the conserved vector (35), we obtain ϕ() [ (u) + (v)]d =. (73) Differentiating (73) under integration sign [12], we have d d ϕ() [u]d u(, ϕ()) dϕ() d + [v(, )]ϕ() =. (74) The boundar condition (69) for v(, ) and epression (71) for v(, ϕ()), reduces (74) to ϕ() ud = constant, independent of, (75) which gives that the total mass flu is constant along the jet. Therefore F = ϕ() ud (76)

13 754 Conservation Laws for Laminar Radial Jets with Swirl is the conserved quantit for the radial liquid jet with swirl [5] Conserved quantities for the radial free jet with swirl The boundar conditions for the radial free jet with swirl are = : v =, u =, w =, (77) = ± : u =, u =, w =, w =. (78) For the conserved vector (36), equation (68) on integration with respect to from = to = with kept as fied during the integration, ields d d [2 (uw)d] + [ ] 2 vw ν 2 w =. (79) We assume that v(, ± ) is finite. The boundar conditions (77) and (78) on (79) gives 2 2 uwd = constant, independent of. (8) The conserved quantit for the radial free jet with swirl is F = 2 2 uwd, (81) where F is a constant independent of. If ρ is the densit of fluid, then the constant ρf is the total radial flu of angular momentum about = (Rile [5]). The second conserved quantit which is valid onl when there is no swirl or at large distance from the orifice can be obtained from the conserved vector (64). Integrate (68) with respect to for the conserved vector (64), we obtain d d [ u 2 d] + [uv νu ] =. (82) The boundar condition (78) and v(, ± ) is finite, ields Thus 2 G = 2 u 2 d = constant, independent of. (83) u 2 d (84) and ρg is the second conserved quantit for the radial free jet with swirl and is valid onl at large distance from orifice.

14 R. Naz, I. Naeem and F. M. Mahomed 755 The conserved quantities for the radial free jet with swirl can also be derived for the stream function formulation. We now outline an alternative derivation of the conserved quantities for the radial free jet with swirl. Also we derive the conserved quantities for the radial wall jet with swirl. In terms of the stream function the boundar conditions (77) and (78) take the following form: = : ψ =, ψ =, w =, (85) = ± : ψ =, ψ =, w =, w =. (86) Integrate equation (68) with respect to form = to = keeping as fied. For the conserved vector (56), we have d d [ wψ d] + [ ] wψ ν 2 w =. (87) Equation (87) with the help of boundar condition (86) and requirement ψ (, ± ) = v(, ± ) is finite, gives 2 wψ d = constant, independent of. (88) Equation (88) is equivalent to (8) and thus we finall obtain the conserved quantit F given in (81). The second conserved quantit which is valid onl when there is no swirl or at large distance from the orifice can be obtained from conserved vector (65). Integration of (68) with respect to with kept as fied for the conserved vector (65) gives d d [ 1 ψ d] 2 + [ 1 ψ ψ νψ ] =. (89) Equation (89) with the help of boundar condition (86) and the requirement ψ (, ± ) = v(, ± ) is finite, ields 2 ψ 2 d = constant, independent of. (9) Equation (9) is equivalent to (83) and hence we obtain the conserved quantit (84) Conserved quantities for the radial wall jet with swirl For the radial wall jet with swirl the boundar conditions are = : u =, v =, w =, (91)

15 756 Conservation Laws for Laminar Radial Jets with Swirl = : u =, u =, w =, w =. (92) In terms of the stream function, we obtain = : ψ =, ψ =, w =, (93) = : ψ =, ψ =, w =, w =. (94) B integrating equation (68) with respect to from = to = keeping as fied and considering the conserved vector (59), we have d d [ ψwψ d] + [ ] wψψ ν 2 ψw + ν 2 w ψ d =. (95) But ψ(, ) = whereas ψ(, ), ψ (, ) and w (, ) are assumed to be finite. Thus boundar conditions (93) and (94) reduce equation (95) to [ d ] ψwψ d + ν 2 w ψ d =. (96) d Integrating equation (96) b parts and using boundar conditions, we deduce ( d ( ) ψ wψ d )d = ν 2 wψ d, (97) d and in terms of the velocit components, we have d ( ) ( 3 u wud )d = ν 3 wu d. (98) d To make the right side of equation (98) zero, Rile [5] assumed that u/w is a function of onl and this condition is satisfied b the similarit solution. Thus in the similarit solution regime, we have ( ) 3 u wud d = constant, independent of. (99) Therefore M = 3 ( ) u wud d, (1) is a conserved quantit. Rile [5] used ρ 2 M as a conserved quantit for the radial wall jet with swirl. There is no obvious phsical significance of this conserved quantit. The second conserved quantit is derived with no swirl or at large distance from the orifice and the conserved vector (66) is used in its derivation. Integrating (68)

16 R. Naz, I. Naeem and F. M. Mahomed 757 with respect to from = to = keeping as fied during the integration, we have d d [ 1 ψψ 2 d] + [ 1 ψψ ψ + ν 2 ψ2 νψψ ] The boundar conditions (93) and (94) ield 1 =. (11) ψψ 2 d = constant, independent of, (12) we have used the requirements ψ(, ) =, ψ(, ) and ψ (, ) are finite. Integrating (12) b parts, we obtain N = 1 or equivalentl N = 2 ( ) ψ ψ 2 d d, (13) ( u u 2 d )d, (14) which is valid onl for no swirl or at large distance from orifice. The constant ρ 2 N is interpreted as the flu of eterior momentum flu Comparison between conservation laws for radial jets with swirl and purel radial jet flows The results for radial jets with swirl and for purel radial jets derived in [7] are displaed in Table 1. The boundar laer equations for the velocit components as well as for the stream function formulation are given in Table 1 for both the radial jets with swirl and for the purel radial jets. The conservation laws for the purel radial jet flows hold for the radial jets with swirl but onl at large distances from the orifice (as ). The multiplier approach on the sstem of equations for the velocit components for purel radial jets ields two local conserved vectors and for the stream function formulation also two local conserved vectors are obtained. The multiplier approach on the sstem of equations for the velocit components for radial jets with swirl ields three conserved vectors, one of these conserved vectors is valid onl when. For the stream function formulation for the radial jets with swirl four conserved vectors are obtained and two of them are valid onl when. The conserved vectors for the radial jets with swirl as and for the purel radial jets are the same. The conserved quantit for the radial liquid jet with swirl and the purel radial liquid jet are the same. The onl conserved quantit for the purel radial free jet

17 758 Conservation Laws for Laminar Radial Jets with Swirl is that the total radial flu of momentum is constant in the jet. The radial free jet with swirl has two conserved quantities. One conserved quantit for the radial free jet with swirl is the total radial flu of angular momentum about = which is constant. The second conserved quantit is the same as for the purel radial free jet but it holds onl as. The conserved quantit for a purel radial wall jet is the flu of eterior momentum flu in the wall jet which is constant and holds for the radial wall jet with swirl at. Another conserved quantit for the radial wall jet with swirl is established with the help of a non-local conserved vector for the sstem of two partial differential equations for the stream function formulation. There is no obvious phsical interpretation of that conserved quantit.

18 R. Naz, I. Naeem and F. M. Mahomed 759 Boundar laer equations Radial jets with swirl Purel radial jets uu + vu w2 = νu uu + vu = νu uw + vw + uw = νw (u) + (v) = (u) + (v) = Conserved vectors T 1 = 2 uw for T 2 = 2 (vw νw ) velocit components T 1 = u T 1 = u T 2 = v T 2 = v As T 1 = u 2 T 1 = u 2 T 2 = (uv νu ) T 2 = (uv νu ) Stream function 1 ψ 1 ψ 2 2 1ψ ψ 1 ψ 1 ψ 2 2 formulation νψ w 2 = 1ψ ψ νψ = 1 ψ w 1ψ w + 1 ψ 2 w νw = Conserved vectors for stream function T 1 = wψ T 2 = wψ ν 2 w T 1 = ψwψ T 2 = wψψ ν 2 ψw + ν2 w ψ As T 1 = 1 ψ2 T 2 = 1ψ ψ νψ T 1 = 1 ψψ2 T 2 = 1ψψ ψ + ν 2 ψ2 νψψ T 1 = 1 ψ2 T 2 = 1ψ ψ νψ T 1 = 1 ψψ2 T 2 = 1ψψ ψ + ν 2 ψ2 νψψ

19 76 Conservation Laws for Laminar Radial Jets with Swirl Conserved quantities For liquid jet For free jet ϕ() ud ϕ() ud 2 2 uwd 2 u 2 d as 2 u 2 d For wall jet ( 3 u wud )d ( 2 ) u u 2 d d as ( 2 u u 2 d )d Table 1: Comparison between radial jets with swirl and purel radial jets 4. CONCLUDING REMARKS The conserved quantities for the radial liquid, free and wall jets with swirl derived b Rile [5] can be constructed with the help of conservation laws. The radial jets with swirl behave as purel radial jets at large distances from the orifice. At large distances from the orifice as, the swirling component w and the boundar laer equations for radial jets with swirl reduce to that of purel radial jets. Therefore, as the conservation laws for radial jets with swirl coincide with those of the purel radial jets. The multiplier approach gave two local conservation laws for the sstem of equations for the velocit components. One of the conserved vectors gave the conserved quantit for the radial liquid jet with swirl and the second conserved vector ielded one conserved quantit for the radial free jet with swirl. Also, another conserved vector for this sstem was obtained as. The second conserved quantit for the radial free jet with swirl was obtained from that conserved vector. The conserved quantities for the wall jet with swirl cannot be obtained from the conserved vectors for the sstem of equations for the velocit components. For the sstem of two partial differential equations for the stream function formulation one local and one non-local conserved vector were obtained. The local conserved vector gave one conserved quantit for the radial free jet with swirl and the non-local conserved vector was used to establish one conserved quantit for the radial wall jet with swirl. Also, two conserved vectors for the stream function for-

20 R. Naz, I. Naeem and F. M. Mahomed 761 mulation were obtained as which were used to derive the second conserved quantities for the radial free and wall jets with swirl. It is of interest to notice that the non-local conserved vector for the sstem of two partial differential equations for the stream function formulation cannot be obtained b the multiplier approach. The reason is that the multiplier approach onl gives multipliers for the local conserved vectors. It ma be possible that the non-local conservation law which does not arise from local multipliers can be derived through a potential sstem associated with the sstem (1)-(3) or the sstem (8)-(9). 5. REFERENCES [1] E. J. Watson, The radial spread of a liquid jet over a horizontal plane, J. Fluid Mech, 2, , [2] W. H. Schwarz, The radial free jet, Chem. Eng, Sci, 18, , [3] M. B. Glauert, The wall jet, J. Fluid Mech,1, , [4] N. Rile, Asmptotic epansions in radial jets, J. Math. Phs, 41, , [5] N. Rile, Radial jets with swirl, Part I. Incompressible Flow, Quart. J. Mech. and Appl. Math, 15, , [6] M. O Nan, W. H. Schwarz, The swirling radial free jet, Appl. Sci. Res. 15, , [7] R. Naz, D. P. Mason and F. M. Mahomed, Conservation laws and conserved quantities for laminar two-dimensional and radial jets, Nonlinear Analsis: Real World Applications, 1, , 29. [8] H. Steudel, Uber die Zuordnung Zwischen Invarianzeigenschaften und Erhaltungssatzen, Zeit Naturforsch, 17, , [9] P. J. Olver, Applications of Lie Groups to Differential Equations, Springer, New York, [1] R. Naz, F. M. Mahomed, D. P. Mason, Comparison of different approaches to conservation laws for some partial differential equations in fluid mechanics, Appl. Math.Comp, 25, , 28. [11] G. K. Batchelor, An Introduction to Fluid Dnamics, Cambridge Universit Press, 75-79, [12] R. P. Gillespie, Integration, Oliver and Bod, Edinburgh, , 1967.

THE HEATED LAMINAR VERTICAL JET IN A LIQUID WITH POWER-LAW TEMPERATURE DEPENDENCE OF DENSITY. V. A. Sharifulin.

THE HEATED LAMINAR VERTICAL JET IN A LIQUID WITH POWER-LAW TEMPERATURE DEPENDENCE OF DENSITY. V. A. Sharifulin. THE HEATED LAMINAR VERTICAL JET IN A LIQUID WITH POWER-LAW TEMPERATURE DEPENDENCE OF DENSITY 1. Introduction V. A. Sharifulin Perm State Technical Universit, Perm, Russia e-mail: sharifulin@perm.ru Water

More information

Fluid Mechanics II. Newton s second law applied to a control volume

Fluid Mechanics II. Newton s second law applied to a control volume Fluid Mechanics II Stead flow momentum equation Newton s second law applied to a control volume Fluids, either in a static or dnamic motion state, impose forces on immersed bodies and confining boundaries.

More information

Joule Heating Effects on MHD Natural Convection Flows in Presence of Pressure Stress Work and Viscous Dissipation from a Horizontal Circular Cylinder

Joule Heating Effects on MHD Natural Convection Flows in Presence of Pressure Stress Work and Viscous Dissipation from a Horizontal Circular Cylinder Journal of Applied Fluid Mechanics, Vol. 7, No., pp. 7-3, 04. Available online at www.jafmonline.net, ISSN 735-357, EISSN 735-3645. Joule Heating Effects on MHD Natural Convection Flows in Presence of

More information

Available online at ScienceDirect. Procedia Engineering 90 (2014 )

Available online at   ScienceDirect. Procedia Engineering 90 (2014 ) Available online at.sciencedirect.com ScienceDirect Procedia Engineering 9 (14 383 388 1th International Conference on Mechanical Engineering, ICME 13 Effects of volumetric heat source and temperature

More information

Incompressible Viscous Flows

Incompressible Viscous Flows Incompressible Viscous Flows Q. Choose the crect answer (i) The maximum velocit of a one-dimensional incompressible full developed viscous flow between two fixed parallel plates is 6m/s. The mean velocit

More information

Symmetry Solutions of a Third-Order Ordinary Differential Equation which Arises from Prandtl Boundary Layer Equations

Symmetry Solutions of a Third-Order Ordinary Differential Equation which Arises from Prandtl Boundary Layer Equations Journal of Nonlinear Mathematical Physics Volume 15, Supplement 1 (2008), 179 191 Article Symmetry Solutions of a Third-Order Ordinary Differential Equation which Arises from Prandtl Boundary Layer Equations

More information

1.1 The Equations of Motion

1.1 The Equations of Motion 1.1 The Equations of Motion In Book I, balance of forces and moments acting on an component was enforced in order to ensure that the component was in equilibrium. Here, allowance is made for stresses which

More information

CONSERVATION OF ANGULAR MOMENTUM FOR A CONTINUUM

CONSERVATION OF ANGULAR MOMENTUM FOR A CONTINUUM Chapter 4 CONSERVATION OF ANGULAR MOMENTUM FOR A CONTINUUM Figure 4.1: 4.1 Conservation of Angular Momentum Angular momentum is defined as the moment of the linear momentum about some spatial reference

More information

International Journal of Mathematical Engineering and Science ISSN : Volume 1 Issue 2

International Journal of Mathematical Engineering and Science ISSN : Volume 1 Issue 2 ISSN : 77-698 Volume Issue The Mathematical Structure and Analsis of an MHD Flow in Porous Media () () Naji Qatanani and Mai Musmar () Department of Mathematics, An-Najah National Universit Nablus-Palestine

More information

UNIT 4 HEAT TRANSFER BY CONVECTION

UNIT 4 HEAT TRANSFER BY CONVECTION UNIT 4 HEAT TRANSFER BY CONVECTION 4.1 Introduction to convection 4. Convection boundar laers 4..1 Hdrodnamic boundar laer over flat plate 4.. Thermal boundar laer 4..3 Concentration boundar laer 4.3 Dimensional

More information

Conservation of Linear Momentum

Conservation of Linear Momentum Conservation of Linear Momentum Once we have determined the continuit equation in di erential form we proceed to derive the momentum equation in di erential form. We start b writing the integral form of

More information

Apply mass and momentum conservation to a differential control volume. Simple classical solutions of NS equations

Apply mass and momentum conservation to a differential control volume. Simple classical solutions of NS equations Module 5: Navier-Stokes Equations: Appl mass and momentum conservation to a differential control volume Derive the Navier-Stokes equations Simple classical solutions of NS equations Use dimensional analsis

More information

A first contact with STAR-CCM+ Comparing analytical and finite volume solutions with STAR-CCM+ simulations

A first contact with STAR-CCM+ Comparing analytical and finite volume solutions with STAR-CCM+ simulations A first contact with STAR-CCM+ Comparing analtical and finite volume solutions with STAR-CCM+ simulations Michael Heer Analtical Finite volumes STAR-CCM+ What is ParisTech? ParisTech is a consortium of

More information

UNCORRECTED. To recognise the rules of a number of common algebraic relations: y = x 1 y 2 = x

UNCORRECTED. To recognise the rules of a number of common algebraic relations: y = x 1 y 2 = x 5A galler of graphs Objectives To recognise the rules of a number of common algebraic relations: = = = (rectangular hperbola) + = (circle). To be able to sketch the graphs of these relations. To be able

More information

CHAPTER-III CONVECTION IN A POROUS MEDIUM WITH EFFECT OF MAGNETIC FIELD, VARIABLE FLUID PROPERTIES AND VARYING WALL TEMPERATURE

CHAPTER-III CONVECTION IN A POROUS MEDIUM WITH EFFECT OF MAGNETIC FIELD, VARIABLE FLUID PROPERTIES AND VARYING WALL TEMPERATURE CHAPER-III CONVECION IN A POROUS MEDIUM WIH EFFEC OF MAGNEIC FIELD, VARIABLE FLUID PROPERIES AND VARYING WALL EMPERAURE 3.1. INRODUCION Heat transer studies in porous media ind applications in several

More information

Chapter 2 Basic Conservation Equations for Laminar Convection

Chapter 2 Basic Conservation Equations for Laminar Convection Chapter Basic Conservation Equations for Laminar Convection Abstract In this chapter, the basic conservation equations related to laminar fluid flow conservation equations are introduced. On this basis,

More information

MMJ1153 COMPUTATIONAL METHOD IN SOLID MECHANICS PRELIMINARIES TO FEM

MMJ1153 COMPUTATIONAL METHOD IN SOLID MECHANICS PRELIMINARIES TO FEM B Course Content: A INTRODUCTION AND OVERVIEW Numerical method and Computer-Aided Engineering; Phsical problems; Mathematical models; Finite element method;. B Elements and nodes, natural coordinates,

More information

7. TURBULENCE SPRING 2019

7. TURBULENCE SPRING 2019 7. TRBLENCE SPRING 2019 7.1 What is turbulence? 7.2 Momentum transfer in laminar and turbulent flow 7.3 Turbulence notation 7.4 Effect of turbulence on the mean flow 7.5 Turbulence generation and transport

More information

International Journal of Applied Mathematics and Physics, 3(2), July-December 2011, pp Global Research Publications, India

International Journal of Applied Mathematics and Physics, 3(2), July-December 2011, pp Global Research Publications, India International Journal of Applied Mathematics and Phsics, 3(), Jul-December 0, pp. 55-67 Global Research Publications, India Effects of Chemical Reaction with Heat and Mass Transfer on Peristaltic Flow

More information

Unsteady free MHD convection flow past a vertical porous plate in slip-flow regime under fluctuating thermal and mass diffusion *

Unsteady free MHD convection flow past a vertical porous plate in slip-flow regime under fluctuating thermal and mass diffusion * CHAPTER 4 Unstead free MHD convection flow past a vertical porous plate in slip-flow regime under fluctuating thermal and mass diffusion * 4. Nomenclature A : Suction parameter C : Dimensionless species

More information

Double reduction of partial differential equations with applications to laminar jets and wakes

Double reduction of partial differential equations with applications to laminar jets and wakes Double reduction of partial differential equations with applications to laminar jets and wakes LADY NOMVULA KOKELA A Dissertation submitted to the Faculty of Science, University of the Witwatersrand, Johannesburg,

More information

Hurricane Modeling E XPANDING C ALCULUS H ORIZON

Hurricane Modeling E XPANDING C ALCULUS H ORIZON Februar 5, 2009 :4 Hurricane Modeling E XPANDING THE Sheet number Page number can magenta ellow black C ALCULUS H ORIZON Hurricane Modeling... Each ear population centers throughout the world are ravaged

More information

Numerical Analysis of Heat Transfer in the Unsteady Flow of a non- Newtonian Fluid over a Rotating cylinder

Numerical Analysis of Heat Transfer in the Unsteady Flow of a non- Newtonian Fluid over a Rotating cylinder Bulletin of Environment, Pharmacolog and Life Sciences Bull. Env.Pharmacol. Life Sci., Vol 4 [Spl issue 1] 215: 318-323 214 Academ for Environment and Life Sciences, India Online ISSN 2277-188 Journal

More information

Magnetohydrodynamic Fluid Flow Between Two Parallel Infinite Plates Subjected To An Inclined Magnetic Field Under Pressure Gradient

Magnetohydrodynamic Fluid Flow Between Two Parallel Infinite Plates Subjected To An Inclined Magnetic Field Under Pressure Gradient Journal of Multidisciplinar Engineering Science and Technolog (JMEST) Magnetohdrodnamic Fluid Flow Between Two Parallel Infinite Plates Subjected To An Inclined Magnetic Field Under Pressure Gradient Agnes

More information

Numerical Solution of Non-Darcian Effects on Natural Convection in a Rectangular Porous Enclosure with Heated Walls

Numerical Solution of Non-Darcian Effects on Natural Convection in a Rectangular Porous Enclosure with Heated Walls International Journal of Advanced Mechanical Engineering. ISSN 2250-3234 Volume 8, Number 1 (18), pp. 71-86 Research India Publications http://www.ripublication.com Numerical Solution of Non-Darcian Effects

More information

Module 1 : The equation of continuity. Lecture 4: Fourier s Law of Heat Conduction

Module 1 : The equation of continuity. Lecture 4: Fourier s Law of Heat Conduction 1 Module 1 : The equation of continuit Lecture 4: Fourier s Law of Heat Conduction NPTEL, IIT Kharagpur, Prof. Saikat Chakrabort, Department of Chemical Engineering Fourier s Law of Heat Conduction According

More information

4 Strain true strain engineering strain plane strain strain transformation formulae

4 Strain true strain engineering strain plane strain strain transformation formulae 4 Strain The concept of strain is introduced in this Chapter. The approimation to the true strain of the engineering strain is discussed. The practical case of two dimensional plane strain is discussed,

More information

LECTURE NOTES - VIII. Prof. Dr. Atıl BULU

LECTURE NOTES - VIII. Prof. Dr. Atıl BULU LECTURE NOTES - VIII «LUID MECHNICS» Istanbul Technical Universit College of Civil Engineering Civil Engineering Department Hdraulics Division CHPTER 8 DIMENSIONL NLYSIS 8. INTRODUCTION Dimensional analsis

More information

LATERAL BUCKLING ANALYSIS OF ANGLED FRAMES WITH THIN-WALLED I-BEAMS

LATERAL BUCKLING ANALYSIS OF ANGLED FRAMES WITH THIN-WALLED I-BEAMS Journal of arine Science and J.-D. Technolog, Yau: ateral Vol. Buckling 17, No. Analsis 1, pp. 9-33 of Angled (009) Frames with Thin-Walled I-Beams 9 ATERA BUCKING ANAYSIS OF ANGED FRAES WITH THIN-WAED

More information

KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS

KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS Chapter 8 KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS Figure 8.1: 195 196 CHAPTER 8. KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS 8.1 Motivation In Chapter 3, the conservation of linear momentum for a

More information

UNSTEADY LOW REYNOLDS NUMBER FLOW PAST TWO ROTATING CIRCULAR CYLINDERS BY A VORTEX METHOD

UNSTEADY LOW REYNOLDS NUMBER FLOW PAST TWO ROTATING CIRCULAR CYLINDERS BY A VORTEX METHOD Proceedings of the 3rd ASME/JSME Joint Fluids Engineering Conference Jul 8-23, 999, San Francisco, California FEDSM99-8 UNSTEADY LOW REYNOLDS NUMBER FLOW PAST TWO ROTATING CIRCULAR CYLINDERS BY A VORTEX

More information

Conservation of Linear Momentum for a Differential Control Volume

Conservation of Linear Momentum for a Differential Control Volume Conservation of Linear Momentum for a Differential Control Volume When we applied the rate-form of the conservation of mass equation to a differential control volume (open sstem in Cartesian coordinates,

More information

Boundary-Fitted Coordinates!

Boundary-Fitted Coordinates! Computational Fluid Dnamics http:wwwndedu~gtrggvacfdcourse Computational Fluid Dnamics Computational Methods or Domains with Comple BoundariesI Grétar Trggvason Spring For most engineering problems it

More information

DIGITAL CORRELATION OF FIRST ORDER SPACE TIME IN A FLUCTUATING MEDIUM

DIGITAL CORRELATION OF FIRST ORDER SPACE TIME IN A FLUCTUATING MEDIUM DIGITAL CORRELATION OF FIRST ORDER SPACE TIME IN A FLUCTUATING MEDIUM Budi Santoso Center For Partnership in Nuclear Technolog, National Nuclear Energ Agenc (BATAN) Puspiptek, Serpong ABSTRACT DIGITAL

More information

Magnetic Field and Chemical Reaction Effects on Convective Flow of

Magnetic Field and Chemical Reaction Effects on Convective Flow of Communications in Applied Sciences ISSN 221-7372 Volume 1, Number 1, 213, 161-187 Magnetic Field and Chemical Reaction Effects on Convective Flow of Dust Viscous Fluid P. Mohan Krishna 1, Dr.V.Sugunamma

More information

5. Nonholonomic constraint Mechanics of Manipulation

5. Nonholonomic constraint Mechanics of Manipulation 5. Nonholonomic constraint Mechanics of Manipulation Matt Mason matt.mason@cs.cmu.edu http://www.cs.cmu.edu/~mason Carnegie Mellon Lecture 5. Mechanics of Manipulation p.1 Lecture 5. Nonholonomic constraint.

More information

Chemical Reaction, Radiation and Dufour Effects on Casson Magneto Hydro Dynamics Fluid Flow over A Vertical Plate with Heat Source/Sink

Chemical Reaction, Radiation and Dufour Effects on Casson Magneto Hydro Dynamics Fluid Flow over A Vertical Plate with Heat Source/Sink Global Journal of Pure and Applied Mathematics. ISSN 973-768 Volume, Number (6), pp. 9- Research India Publications http://www.ripublication.com Chemical Reaction, Radiation and Dufour Effects on Casson

More information

Elliptic Equations. Chapter Definitions. Contents. 4.2 Properties of Laplace s and Poisson s Equations

Elliptic Equations. Chapter Definitions. Contents. 4.2 Properties of Laplace s and Poisson s Equations 5 4. Properties of Laplace s and Poisson s Equations Chapter 4 Elliptic Equations Contents. Neumann conditions the normal derivative, / = n u is prescribed on the boundar second BP. In this case we have

More information

Chapter 6 Laminar External Flow

Chapter 6 Laminar External Flow Chapter 6 aminar Eternal Flow Contents 1 Thermal Boundary ayer 1 2 Scale analysis 2 2.1 Case 1: δ t > δ (Thermal B.. is larger than the velocity B..) 3 2.2 Case 2: δ t < δ (Thermal B.. is smaller than

More information

ME 321: FLUID MECHANICS-I

ME 321: FLUID MECHANICS-I 6/07/08 ME 3: LUID MECHANI-I Dr. A.B.M. Toufique Hasan Professor Department of Mechanical Engineering Bangladesh Universit of Engineering & Technolog (BUET), Dhaka Lecture- 4/07/08 Momentum Principle teacher.buet.ac.bd/toufiquehasan/

More information

Nonlinear Systems Examples Sheet: Solutions

Nonlinear Systems Examples Sheet: Solutions Nonlinear Sstems Eamples Sheet: Solutions Mark Cannon, Michaelmas Term 7 Equilibrium points. (a). Solving ẋ =sin 4 3 =for gives =as an equilibrium point. This is the onl equilibrium because there is onl

More information

Exercise solutions: concepts from chapter 7

Exercise solutions: concepts from chapter 7 f () = -N F = +N f (1) = +4N Fundamentals of Structural Geolog 1) In the following eercise we consider some of the phsical quantities used in the stud of particle dnamics and review their relationships

More information

CHAPTER 2: Partial Derivatives. 2.2 Increments and Differential

CHAPTER 2: Partial Derivatives. 2.2 Increments and Differential CHAPTER : Partial Derivatives.1 Definition of a Partial Derivative. Increments and Differential.3 Chain Rules.4 Local Etrema.5 Absolute Etrema 1 Chapter : Partial Derivatives.1 Definition of a Partial

More information

Mathematical Modeling on MHD Free Convective Couette Flow between Two Vertical Porous Plates with Chemical Reaction and Soret Effect

Mathematical Modeling on MHD Free Convective Couette Flow between Two Vertical Porous Plates with Chemical Reaction and Soret Effect International Journal of Advanced search in Phsical Science (IJARPS) Volume, Issue 5, September 4, PP 43-63 ISSN 349-7874 (Print) & ISSN 349-788 (Online) www.arcjournals.org Mathematical Modeling on MHD

More information

Symmetry Arguments and the Role They Play in Using Gauss Law

Symmetry Arguments and the Role They Play in Using Gauss Law Smmetr Arguments and the Role The la in Using Gauss Law K. M. Westerberg (9/2005) Smmetr plas a ver important role in science in general, and phsics in particular. Arguments based on smmetr can often simplif

More information

Effects of wall properties and heat transfer on the peristaltic transport of a jeffrey fluid in a channel

Effects of wall properties and heat transfer on the peristaltic transport of a jeffrey fluid in a channel Available online at www.pelagiaresearchlibrar.com Advances in Applied Science Research,, 4(6):59-7 ISSN: 976-86 CODEN (USA): AASRFC Effects of wall properties and heat transfer on the peristaltic transport

More information

AN EFFICIENT MOVING MESH METHOD FOR A MODEL OF TURBULENT FLOW IN CIRCULAR TUBES *

AN EFFICIENT MOVING MESH METHOD FOR A MODEL OF TURBULENT FLOW IN CIRCULAR TUBES * Journal of Computational Mathematics, Vol.27, No.2-3, 29, 388 399. AN EFFICIENT MOVING MESH METHOD FOR A MODEL OF TURBULENT FLOW IN CIRCULAR TUBES * Yin Yang Hunan Ke Laborator for Computation and Simulation

More information

Bubble curtain system for prevention of seawater intrusion in coastal aquifers

Bubble curtain system for prevention of seawater intrusion in coastal aquifers Available online at www.pelagiaresearchlibrar.com Advances in Applied Science Research, 013, 4(4):81-90 ISSN: 0976-8610 CODEN (USA): AASRFC Bubble curtain sstem for prevention of seawater intrusion in

More information

Chapter 4 Transport of Pollutants

Chapter 4 Transport of Pollutants 4- Introduction Phs. 645: Environmental Phsics Phsics Department Yarmouk Universit hapter 4 Transport of Pollutants - e cannot avoid the production of pollutants. hat can we do? - Transform pollutants

More information

Candidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request.

Candidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request. UNIVERSITY OF EAST ANGLIA School of Mathematics Spring Semester Examination 2004 FLUID DYNAMICS Time allowed: 3 hours Attempt Question 1 and FOUR other questions. Candidates must show on each answer book

More information

MAE 323: Chapter 4. Plane Stress and Plane Strain. The Stress Equilibrium Equation

MAE 323: Chapter 4. Plane Stress and Plane Strain. The Stress Equilibrium Equation The Stress Equilibrium Equation As we mentioned in Chapter 2, using the Galerkin formulation and a choice of shape functions, we can derive a discretized form of most differential equations. In Structural

More information

Velocity Limit in DPD Simulations of Wall-Bounded Flows

Velocity Limit in DPD Simulations of Wall-Bounded Flows Velocit Limit in DPD Simulations of Wall-Bounded Flows Dmitr A. Fedosov, Igor V. Pivkin and George Em Karniadakis Division of Applied Mathematics, Brown Universit, Providence, RI 2912 USA Abstract Dissipative

More information

THE GENERAL ELASTICITY PROBLEM IN SOLIDS

THE GENERAL ELASTICITY PROBLEM IN SOLIDS Chapter 10 TH GNRAL LASTICITY PROBLM IN SOLIDS In Chapters 3-5 and 8-9, we have developed equilibrium, kinematic and constitutive equations for a general three-dimensional elastic deformable solid bod.

More information

Selected Papers from the WSEAS Conferences in Spain, September 2008 Santander, Cantabria, Spain, September 23-25, 2008

Selected Papers from the WSEAS Conferences in Spain, September 2008 Santander, Cantabria, Spain, September 23-25, 2008 Comparisons between Velocit Profiles According to the Modified and the Navier-Stokes Equations of Motion and the Eperimental Measurements for Laminar Boundar Laer Over a Flat Plate Matthew J. Inkman, and

More information

Aircraft Structures Structural & Loading Discontinuities

Aircraft Structures Structural & Loading Discontinuities Universit of Liège Aerospace & Mechanical Engineering Aircraft Structures Structural & Loading Discontinuities Ludovic Noels Computational & Multiscale Mechanics of Materials CM3 http://www.ltas-cm3.ulg.ac.be/

More information

A Simple Turbulence Closure Model

A Simple Turbulence Closure Model A Simple Turbulence Closure Model Atmospheric Sciences 6150 1 Cartesian Tensor Notation Reynolds decomposition of velocity: Mean velocity: Turbulent velocity: Gradient operator: Advection operator: V =

More information

NONLINEAR DYNAMICS AND CHAOS. Numerical integration. Stability analysis

NONLINEAR DYNAMICS AND CHAOS. Numerical integration. Stability analysis LECTURE 3: FLOWS NONLINEAR DYNAMICS AND CHAOS Patrick E McSharr Sstems Analsis, Modelling & Prediction Group www.eng.o.ac.uk/samp patrick@mcsharr.net Tel: +44 83 74 Numerical integration Stabilit analsis

More information

Infinitesimal Rotations

Infinitesimal Rotations Universit of Connecticut DigitalCommons@UConn Chemistr Education Materials Department of Chemistr Januar 007 Infinitesimal Rotations Carl W. David Universit of Connecticut, Carl.David@uconn.edu Follow

More information

On the Effect of the Boundary Conditions and the Collision Model on Rarefied Gas Flows

On the Effect of the Boundary Conditions and the Collision Model on Rarefied Gas Flows On the Effect of the Boundar Conditions and the Collision Model on Rarefied Gas Flows M. Vargas a, S. Stefanov a and D. Valougeorgis b a Institute of Mechanics, Bulgarian Academ of Sciences, Acad. G. Bonchev

More information

The Basic (Physical) Problem: Experiments. Internal Gravity Waves and Hyperbolic Boundary-Value Problems. Internal Gravity Waves. Governing Equations

The Basic (Physical) Problem: Experiments. Internal Gravity Waves and Hyperbolic Boundary-Value Problems. Internal Gravity Waves. Governing Equations The Basic (Phsical) Problem: Eperiments Internal Gravit Waves and Hperbolic Boundar-Value Problems P. A. Martin Colorado chool of Mines Collaborator: tefan Llewelln mith, UC an Diego D. E. Mowbra & B..

More information

The Gas-assisted Expelled Fluid Flow in the Front of a Long Bubble in a Channel

The Gas-assisted Expelled Fluid Flow in the Front of a Long Bubble in a Channel C. H. Hsu, P. C. Chen,. Y. ung, G. C. uo The Gas-assisted Epelled Fluid Flow in the Front of a Long Bubble in a Channel C.H. Hsu 1, P.C. Chen 3,.Y. ung, and G.C. uo 1 1 Department of Mechanical Engineering

More information

Analysis of condenser shell side pressure drop based on the mechanical energy loss

Analysis of condenser shell side pressure drop based on the mechanical energy loss Article Calorifics December 01 Vol.57 No.36: 4718475 doi: 10.1007/s11434-01-5558-4 Analsis of condenser shell side pressure drop based on the mechanical energ loss ZENG Hui, MENG JiAn & LI ZhiXin * Ke

More information

MAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design

MAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design MAE4700/5700 Finite Element Analsis for Mechanical and Aerospace Design Cornell Universit, Fall 2009 Nicholas Zabaras Materials Process Design and Control Laborator Sible School of Mechanical and Aerospace

More information

Three-dimensional numerical simulation of a vortex ring impacting a bump

Three-dimensional numerical simulation of a vortex ring impacting a bump HEOREICAL & APPLIED MECHANICS LEERS 4, 032004 (2014) hree-dimensional numerical simulation of a vorte ring impacting a bump Heng Ren, 1, a) 2, Changue Xu b) 1) China Electronics echnolog Group Corporation

More information

Chemical Reaction and Viscous Dissipation Effects on an Unsteady MHD Flow of Heat and Mass Transfer along a Porous Flat Plate

Chemical Reaction and Viscous Dissipation Effects on an Unsteady MHD Flow of Heat and Mass Transfer along a Porous Flat Plate International Journal of Scientific and Innovative Mathematical Research (IJSIMR) Volume 2 Issue 2 December 24 PP 977-988 ISSN 2347-37X (Print) & ISSN 2347-342 (Online) www.arcjournals.org Chemical Reaction

More information

Conservative fields and potential functions. (Sect. 16.3) The line integral of a vector field along a curve.

Conservative fields and potential functions. (Sect. 16.3) The line integral of a vector field along a curve. onservative fields and potential functions. (Sect. 16.3) eview: Line integral of a vector field. onservative fields. The line integral of conservative fields. Finding the potential of a conservative field.

More information

Chapter 18 KINETICS OF RIGID BODIES IN THREE DIMENSIONS. The two fundamental equations for the motion of a system of particles .

Chapter 18 KINETICS OF RIGID BODIES IN THREE DIMENSIONS. The two fundamental equations for the motion of a system of particles . hapter 18 KINETIS F RIID DIES IN THREE DIMENSINS The to fundamental equations for the motion of a sstem of particles ΣF = ma ΣM = H H provide the foundation for three dimensional analsis, just as the do

More information

Closed form expressions for the gravitational inner multipole moments of homogeneous elementary solids

Closed form expressions for the gravitational inner multipole moments of homogeneous elementary solids Closed form epressions for the gravitational inner multipole moments of homogeneous elementar solids Julian Stirling 1,2, and Stephan Schlamminger 1 1 National Institute of Standards and Technolog, 1 Bureau

More information

An elementary model for the validation of flamelet approximations in non-premixed turbulent combustion

An elementary model for the validation of flamelet approximations in non-premixed turbulent combustion Combust. Theor Modelling 4 () 89. Printed in the UK PII: S364-783()975- An elementar model for the validation of flamelet approimations in non-premied turbulent combustion A Bourliou and A J Majda Département

More information

Large Eddy Simulation of the flow and heat transfer in a half-corrugated channel with various wave amplitudes

Large Eddy Simulation of the flow and heat transfer in a half-corrugated channel with various wave amplitudes Int. J Heat Mass Transfer, Vol. 76 pp. 432-446, 214 Large Edd Simulation of the flow and heat transfer in a half-corrugated channel with various wave amplitudes Maram Mirzaei, Ahmad Sohankar, Lars Davidson,

More information

First Order Equations

First Order Equations 10 1 Linear and Semilinear Equations Chapter First Order Equations Contents 1 Linear and Semilinear Equations 9 Quasilinear Equations 19 3 Wave Equation 6 4 Sstems of Equations 31 1 Linear and Semilinear

More information

ENGI 4430 Advanced Calculus for Engineering Faculty of Engineering and Applied Science Problem Set 3 Solutions [Multiple Integration; Lines of Force]

ENGI 4430 Advanced Calculus for Engineering Faculty of Engineering and Applied Science Problem Set 3 Solutions [Multiple Integration; Lines of Force] ENGI 44 Advanced Calculus for Engineering Facult of Engineering and Applied Science Problem Set Solutions [Multiple Integration; Lines of Force]. Evaluate D da over the triangular region D that is bounded

More information

Directional derivatives and gradient vectors (Sect. 14.5). Directional derivative of functions of two variables.

Directional derivatives and gradient vectors (Sect. 14.5). Directional derivative of functions of two variables. Directional derivatives and gradient vectors (Sect. 14.5). Directional derivative of functions of two variables. Partial derivatives and directional derivatives. Directional derivative of functions of

More information

SAMPLE. A Gallery of Graphs. To recognise the rules of a number of common algebraic relationships: y = x 1,

SAMPLE. A Gallery of Graphs. To recognise the rules of a number of common algebraic relationships: y = x 1, Objectives C H A P T E R 5 A Galler of Graphs To recognise the rules of a number of common algebraic relationships: =, =, = / and + =. To be able to sketch the graphs and simple transformations of these

More information

Simulation of Acoustic and Vibro-Acoustic Problems in LS-DYNA using Boundary Element Method

Simulation of Acoustic and Vibro-Acoustic Problems in LS-DYNA using Boundary Element Method 10 th International LS-DYNA Users Conference Simulation Technolog (2) Simulation of Acoustic and Vibro-Acoustic Problems in LS-DYNA using Boundar Element Method Yun Huang Livermore Software Technolog Corporation

More information

Effects of wavelength and amplitude of a wavy cylinder in cross-flow at low Reynolds numbers

Effects of wavelength and amplitude of a wavy cylinder in cross-flow at low Reynolds numbers J. Fluid Mech. (2009), vol. 620, pp. 195 220. c 2009 Cambridge Universit Press doi:10.1017/s0022112008004217 Printed in the United Kingdom 195 Effects of wavelength and amplitude of a wav clinder in cross-flow

More information

Submarine sand ripples formation in a viscous fluid: 2D and 3D linear stability analysis

Submarine sand ripples formation in a viscous fluid: 2D and 3D linear stability analysis Marine Sandwave and River Dune Dnamics 1 & April 4 - Enschede, the Netherlands Submarine sand ripples formation in a viscous fluid: D and 3D linear stabilit analsis V. Langlois (1) and A. Valance (1) Groupe

More information

SOLUTIONS FOR THEORETICAL COMPETITION Theoretical Question 1 (10 points) 1A (3.5 points)

SOLUTIONS FOR THEORETICAL COMPETITION Theoretical Question 1 (10 points) 1A (3.5 points) II International Zhautkov Olmpiad/Theoretical Competition/Solutions Page 1/10 SOLUTIONS FOR THEORETICAL COMPETITION Theoretical Question 1 (10 points) 1A (.5 points) m v Mu m = + w + ( u w ) + mgr It is

More information

VISCO-ELASTIC FLUID FLOW WITH HEAT AND MASS TRASNFER IN A VERTICAL CHANNEL THROUGH A POROUS MEDIUM

VISCO-ELASTIC FLUID FLOW WITH HEAT AND MASS TRASNFER IN A VERTICAL CHANNEL THROUGH A POROUS MEDIUM Volume 2, No. 1, Januar 214 Journal of Global Research in Mathematical Archives RESEARCH PAPER Available online at http://www.jgrma.info VISCO-ELASTIC FLUID FLOW WITH HEAT AND MASS TRASNFER IN A VERTICAL

More information

Fluid Mechanics for International Engineers HW #4: Conservation of Linear Momentum and Conservation of Energy

Fluid Mechanics for International Engineers HW #4: Conservation of Linear Momentum and Conservation of Energy 2141-365 Fluid Mechanics for International Engineers 1 Problem 1 RTT and Time Rate of Change of Linear Momentum and The Corresponding Eternal Force Notation: Here a material volume (MV) is referred to

More information

A Simple Turbulence Closure Model. Atmospheric Sciences 6150

A Simple Turbulence Closure Model. Atmospheric Sciences 6150 A Simple Turbulence Closure Model Atmospheric Sciences 6150 1 Cartesian Tensor Notation Reynolds decomposition of velocity: V = V + v V = U i + u i Mean velocity: V = Ui + V j + W k =(U, V, W ) U i =(U

More information

F11AE1 1. C = ρν r r. r u z r

F11AE1 1. C = ρν r r. r u z r F11AE1 1 Question 1 20 Marks) Consider an infinite horizontal pipe with circular cross-section of radius a, whose centre line is aligned along the z-axis; see Figure 1. Assume no-slip boundary conditions

More information

Heat and Mass Transfer Flow of MHD Viscous Dissipative Fluid in a Channel with a Stretching and Porous Plate

Heat and Mass Transfer Flow of MHD Viscous Dissipative Fluid in a Channel with a Stretching and Porous Plate Appl. Math. Inf. Sci. Lett. 5, No. 3, 8-87 (27) 8 Applied Mathematics & Information Sciences Letters An International Journal http://dx.doi.org/.8576/amisl/53 Heat and Mass Transfer Flow of MHD Viscous

More information

Mechanics Departmental Exam Last updated November 2013

Mechanics Departmental Exam Last updated November 2013 Mechanics Departmental Eam Last updated November 213 1. Two satellites are moving about each other in circular orbits under the influence of their mutual gravitational attractions. The satellites have

More information

AE/ME 339. K. M. Isaac Professor of Aerospace Engineering. December 21, 2001 topic13_grid_generation 1

AE/ME 339. K. M. Isaac Professor of Aerospace Engineering. December 21, 2001 topic13_grid_generation 1 AE/ME 339 Professor of Aerospace Engineering December 21, 2001 topic13_grid_generation 1 The basic idea behind grid generation is the creation of the transformation laws between the phsical space and the

More information

Film thickness Hydrodynamic pressure Liquid saturation pressure or dissolved gases saturation pressure. dy. Mass flow rates per unit length

Film thickness Hydrodynamic pressure Liquid saturation pressure or dissolved gases saturation pressure. dy. Mass flow rates per unit length NOTES DERITION OF THE CLSSICL REYNOLDS EQTION FOR THIN FIL FLOWS Te lecture presents te derivation of te Renolds equation of classical lubrication teor. Consider a liquid flowing troug a tin film region

More information

MATHEMATICS 200 December 2013 Final Exam Solutions

MATHEMATICS 200 December 2013 Final Exam Solutions MATHEMATICS 2 December 21 Final Eam Solutions 1. Short Answer Problems. Show our work. Not all questions are of equal difficult. Simplif our answers as much as possible in this question. (a) The line L

More information

VISCOUS FLOW DUE TO A SHRINKING SHEET

VISCOUS FLOW DUE TO A SHRINKING SHEET QUARTERLY OF APPLIED MATHEMATICS VOLUME, NUMBER 0 XXXX XXXX, PAGES 000 000 S 0000-0000(XX)0000-0 VISCOUS FLOW DUE TO A SHRINKING SHEET By M. MIKLAVČIČ (Department of Mathematics, Michigan State University,

More information

Potential flow of a second-order fluid over a tri-axial ellipsoid

Potential flow of a second-order fluid over a tri-axial ellipsoid Potential flow of a second-order fluid over a tri-axial ellipsoid F. Viana 1, T. Funada, D. D. Joseph 1, N. Tashiro & Y. Sonoda 1 Department of Aerospace Engineering and Mechanics, Universit of Minnesota,

More information

BASE VECTORS FOR SOLVING OF PARTIAL DIFFERENTIAL EQUATIONS

BASE VECTORS FOR SOLVING OF PARTIAL DIFFERENTIAL EQUATIONS BASE VECTORS FOR SOLVING OF PARTIAL DIFFERENTIAL EQUATIONS J. Roubal, V. Havlena Department of Control Engineering, Facult of Electrical Engineering, Czech Technical Universit in Prague Abstract The distributed

More information

HYDROMAGNETIC DIVERGENT CHANNEL FLOW OF A VISCO- ELASTIC ELECTRICALLY CONDUCTING FLUID

HYDROMAGNETIC DIVERGENT CHANNEL FLOW OF A VISCO- ELASTIC ELECTRICALLY CONDUCTING FLUID Rita Choudhury et al. / International Journal o Engineering Science and Technology (IJEST) HYDROAGNETIC DIVERGENT CHANNEL FLOW OF A VISCO- ELASTIC ELECTRICALLY CONDUCTING FLUID RITA CHOUDHURY Department

More information

Transformation of kinematical quantities from rotating into static coordinate system

Transformation of kinematical quantities from rotating into static coordinate system Transformation of kinematical quantities from rotating into static coordinate sstem Dimitar G Stoanov Facult of Engineering and Pedagog in Sliven, Technical Universit of Sofia 59, Bourgasko Shaussee Blvd,

More information

Ordinary Differential Equations of First Order

Ordinary Differential Equations of First Order CHAPTER 1 Ordinar Differential Equations of First Order 1.1 INTRODUCTION Differential equations pla an indispensable role in science technolog because man phsical laws relations can be described mathematicall

More information

Simultaneous Orthogonal Rotations Angle

Simultaneous Orthogonal Rotations Angle ELEKTROTEHNIŠKI VESTNIK 8(1-2): -11, 2011 ENGLISH EDITION Simultaneous Orthogonal Rotations Angle Sašo Tomažič 1, Sara Stančin 2 Facult of Electrical Engineering, Universit of Ljubljana 1 E-mail: saso.tomaic@fe.uni-lj.si

More information

ESCI 485 Air/sea Interaction Lesson 5 Oceanic Boundary Layer

ESCI 485 Air/sea Interaction Lesson 5 Oceanic Boundary Layer ESCI 485 Air/sea Interaction Lesson 5 Oceanic Boundar Laer References: Descriptive Phsical Oceanograph, Pickard and Emer Introductor Dnamical Oceanograph, Pond and Pickard Principles of Ocean Phsics, Apel

More information

arxiv:physics/ v1 [physics.ed-ph] 25 Apr 2000

arxiv:physics/ v1 [physics.ed-ph] 25 Apr 2000 Matching Conditions on Capillar Ripples: Polariation Estudio de las condiciones de empalme para las oscilaciones de intercara entre dos fluidos. Polariación. arxiv:phsics/0004064v1 [phsics.ed-ph] 25 Apr

More information

Liquid Impact on a Permeable Solid Body

Liquid Impact on a Permeable Solid Body Liquid Impact on a Permeable Solid Bod Yuri A. Semenov a, Guo Xiong Wu a a Department of Mechanical Engineering, Universit College London, Torrington Place, London WCE 7JE, UK Abstract The free surface

More information

L O N G I T U D I N A L V O R T I C E S I N G R A N U L A R F L O W S

L O N G I T U D I N A L V O R T I C E S I N G R A N U L A R F L O W S L O N G I T U D I N A L V O R T I C E S I N G R A N U L A R F L O W S Yoël FORTERRE, Olivier POULIQUEN Institut Universitaire des Sstèmes Thermiques et Industriels (UMR 6595 CNRS) 5, rue Enrico Fermi 13453

More information

Heat and Mass Transfer on Unsteady M.H.D Oscillatory Flow of Non-Newtonian Fluid through Porous Medium in Parallel Plate Channel

Heat and Mass Transfer on Unsteady M.H.D Oscillatory Flow of Non-Newtonian Fluid through Porous Medium in Parallel Plate Channel Heat and Mass Transfer on Unstead M.H.D Oscillator Flow of Non-Newtonian Fluid through Porous Medium in Parallel Plate Channel Dr. G. Kathaani, Dr. P. Sambasivudu, D. M. Praveen Babu Assistant Professor,

More information

Solutions to Two Interesting Problems

Solutions to Two Interesting Problems Solutions to Two Interesting Problems Save the Lemming On each square of an n n chessboard is an arrow pointing to one of its eight neighbors (or off the board, if it s an edge square). However, arrows

More information