Module-5: Hypersonic Boundary Layer theory. Lecture-20: Hypersonic boundary equation

Size: px
Start display at page:

Download "Module-5: Hypersonic Boundary Layer theory. Lecture-20: Hypersonic boundary equation"

Transcription

1 Module-5: Hypersonic Boundary Layer theory Lecture-0: Hypersonic boundary equation 0.1 Governing Equations for Viscous Flows The Navier-Stokes (NS) equaadtions are the governing equations for the viscous compressible flow and hence are the governing equations for hypersonic flows. This section deals with the basics of NS equations and its non-dimensionalization. Continuity Equation: Considering steady state conditions we have; ( ρu) ( ρv) + = 0 i.e. (0.1) X Momentum Equation: D( ρu) ρdu udρ = +, now since D ρ = 0; Dt Dt Dt Dt D( ρu) ρdu = ; therefore the L.H.S. simplifies to give X momentum equation in Dt Dt steady state conditions as; u u p v u u u ρu + ρ v = + µ + + µ + (0.) Y Momentum Equation: D( ρv) ρdv vdρ = +, now since D ρ = 0; Dt Dt Dt Dt D( ρv) ρdv = ; therefore the L.H.S. simplifies to give Y momentum equation in Dt Dt steady State Conditions as; Joint initiative of IITs and IISc Funded by MHRD Page 1 of 39

2 v v p v u v v ρu + ρ v = + µ + + µ + (0.3) Energy Equation: ( u + v ) ( u + v ) ( u + v ) D ρ e+ ρd e e Dρ + + = +, Dt Dt Dt Since D ρ = 0; Dt ( u + v ) ( u + v ) D ρ e+ ρd e+ =, Dt Dt Therefore the L.H.S. simplifies to give Energy Equation in steady State Conditions as; ( ) ( ) u e+ = ρq+ κ + κ + + τ + τ + τ + τ V T T pu pv u u v v xx xy xy yy (0.4) The above equations are written for steady, compressible, viscous, two dimensional flows in Cartesian coordinates. Where u and v are velocities in x and y directions respectively; e represents internal energy per unit mass and q represents the volumetric heating that might occur. All other notations carry their usual meaning. We can simplify the above set of equations via appropriate assumptions, and obtain approximate viscous flow results. Joint initiative of IITs and IISc Funded by MHRD Page of 39

3 0. Non Dimensional Form of Governing Equations The non dimensional form of Navier-Stokes equations can be obtained as follows. Here we have considered a two dimensional steady flow and ignored the normal stresses τ xx and τ yy. Reference variables of the flow can be used for nondimensionalization. Non Dimensional Variables: u = u V v = v V x x = c y y c p p = ρv = e = e CT v µ µ = µ κ κ = κ ρ = ρ ρ Where V, T, ρ, µ, κ are free stream parameters and c reference length. Therefore the Non Dimensional Equations are given as: Non Dimensional Continuity Equation: ( ρu) ( ρv) + = 0 (0.5) Non Dimensional X Momentum Equation: ρ u u 1 p 1 u v v u x ρ y γ M x Re y µ + = + + (0.6) Non Dimensional Y Momentum Equation: ρ v v 1 p 1 u v v u x ρ y γ M y Re x µ + = + + (0.7) Joint initiative of IITs and IISc Funded by MHRD Page 3 of 39

4 Non Dimensional Energy Equation: ( ) ( ) e e γ u v ( 1) M u u v v u v T ρ ρ γγ ρ ρ κ κ T + = Pr Re M v u v u + γ ( γ 1) µν + + µ u + (0.8) Re 0.3 Process of Non-dimensionlisation of Governing Equations:- Continuity Equation:- We Have, ( ρu) ( ρv) + = 0, With non-dimensional parameters ρ V ( ρu/ ρ V ) ( ρv/ ρ V ). 0 c + ( x/ c) ( y/ c) = ( ρu) ( ρv) + = 0 X & Y momentum Equations:- The process of non dimensionlisation is almost the same for both X & Y momentum equations, therefore we have for X momentum equation; u u p v u u u ρu + ρ v = + µ + + µ + Rewriting the above Equation with Non Dimensional parameters, we get; ρ V ρu ( u/ V ) ρu ( v/ V ) c + = ρ V ( x/ c) ρ V ( y/ c) p ( p/ p ) c ( x/ c) µ V µ ( v/ V ) ( u/ V ) µ ( v/ V ) ( u/ V ) c ( y/ c) µ ( x/ c) ( y/ c) ( x/ c) µ ( x/ c) ( y/ c) Joint initiative of IITs and IISc Funded by MHRD Page 4 of 39

5 Now dividing the above equation on both sides by the factor wise ρ V c ; we get term p c p First, the pressure term on Right Hand side (R.H.S); = c ρ V ρ V Since Mach Number can be written as M V = therefore; a p ρ V p = and furthermore, p ρ RT ρ M a = and a = γ RT ; the term reduces to p RT 1 = = ρ M a M γrt γm µ V c µ 1 Now the viscous term, = = c ρ V ρ V c Re Thus the X momentum Equation reduces to, ρ u u 1 p 1 u v v u x ρ y γ M x Re y µ + = + + & similarly Y momentum Equation reduces to, ρ v v 1 p 1 u v v u x ρ y γ M y Re x µ + = + + Normal stresses τ xx and τ yy are ignored for the sake of simplicity. Joint initiative of IITs and IISc Funded by MHRD Page 5 of 39

6 Lecture-1: Non-dimensionalisation of governing equations 1.1 Process of Non-dimensionlisation of Governing Equations Energy Equation: We know the energy equation as Eq. (0.4). Expanding the LHS of this equation, we get, ( u v ) ( u v ) ρu + e ρv e+ = T T ( pu) ( pv) ( uτ xx) ( uτ xy) ( vτ xy) ( vτ yy) κ + κ Now simplifying the above equation by neglecting normal stresses τ xx and τ yy, we get e e ( u v ) ( u v ) ρu ρv ρu + ρv + + = + T T ( pu) ( pv) v u v u κ + κ + + v µ + + u µ + Now non dimensionalising Left hand side (L.H.S.) & R.H.S. terms separately; L.H.S.= ρ V CT v ρu ( e/ CT v ) ρv ( e/ CT v ) ρ V CT v e e ρu ρv c + ρ V ( x/ c) ρ V ( y/ c) = + c Non-dimensional R.H.S., 1 st & nd term; 3 ρ V ρu ( u/ V ) + ( v/ V ) ρv ( u/ V ) + ( v/ V ) + c ρ V ( x/ c) ρ V ( y/ c) = ρ V ( u + v ) ( u + v ) ρu + ρv c 3 Joint initiative of IITs and IISc Funded by MHRD Page 6 of 39

7 Non-dimensional R.H.S., 3 rd & 4 th term; κ T κ ( T / T ) κ ( T / T ) κ T T T + = κ κ + c ( x/ c) κ ( x/ c) ( y/ c) κ ( y/ c) c Non-dimensional R.H.S., 6 th & 7 th term reduces to; µ V v µ ( v/ V ) ( u/ V ) µ V v u v µ + = + c ( x/ c) V µ ( x/ c) ( y/ c) c µ V u µ ( v/ V ) ( u/ V ) µ V v u u µ + = + c ( y/ c) V µ ( x/ c) ( y/ c) c Non-dimensional R.H.S., 5 th term reduces to; pv ( p/ pu / V ) pv ( pu) c ( x/ c) = c Now dividing the L.H.S. & R.H.S. by the factor as; c ρ V CT v, we obtain R.H.S. terms For term 1 st & nd ; V c V V ( 1) V ( 1) c ρ V CT CT RT a 3 ρ γ γγ = = = = M γγ v v ( 1) For term 3 rd & 4 th ; κ T c κ µ κγ γ = = = Re Pr c ρ V CvT cρ V Cv cρ V µ Cp For term 5 th ; p V c p RT ( γ 1) = = = ( γ 1) c ρ V CT v ρ CT v RT Joint initiative of IITs and IISc Funded by MHRD Page 7 of 39

8 For term 6 th & 7 th ; µ V c µ V 1 ( γ 1) V 1 γγ ( 1) V γγ ( 1) M = = = = c ρ V CT v cρ V CT v Re RT Re a Re Finally the Non dimensional Energy equation is given as; ( ) ( ) e e γ u v ( 1) M u u v v u v T ρ ρ γγ ρ ρ κ κ T + = Pr Re M v u v u + γ ( γ 1) µν + + µ u + Re The non-dimensional parameters of the Navier-Stokes equation are, Ratio of specific heats: C γ = C p v Mach number: M V = a Reynolds number: Re ρ V µ = c Prandtl number: Pr = µ Cp k These four dimensionless parameters are called similarity parameters, and are very important in determining the nature of a given viscous flow problem. Thermodynamic properties, as reflected by γ, are important for any high speed flow problem. A combination of thermodynamics and flow kinetic energy can be found in M, and it is known that, Flow Kinetic Energy M ~ Flow Internal Energy For Reynolds number, we have Re ~ Inertia Force Viscous Force Joint initiative of IITs and IISc Funded by MHRD Page 8 of 39

9 Prandtl number, appearing in the energy equation, has the physical interpretation as, Pr ~ Frictional dissipation Thermal conduction Here γ and Pr are the properties of gas while Re and M involve flow properties as well. Joint initiative of IITs and IISc Funded by MHRD Page 9 of 39

10 Lecture-: Order of magnitude estimate.1 Boundary conditions An important difference between inviscid and viscous flows can be seen explicitly in the boundary conditions at the wall. The usual boundary condition for an inviscid flow is no mass transfer through the wall which mathematically gets expressed as the normal component of velocity to be zero at the wall. This boundary condition is termed as free slip along the wall. This boundary condition gets added with the cancellation of tangential velocity at the wall due to the existence of friction. This boundary condition is termed as no slip along the wall. Therefore both the components of velocity become zero for viscous wall boundary condition, that is, Wall boundary condition: u=v=0 If there is mass transfer at the wall, then we have to express the normal velocity at the wall as per the known mass flow rate to satisfy the mass conservation equation. However tangential component of velocity will still remain zero at the wall if wall has zero velocity while for conservation of momentum. There are two types of boundary conditions related with the energy equation. In one of them, wall is treated with isothermal wall temperature where the known temperature is assigned at the wall as, Constant wall temperature boundary condition: T=T w Here T w is the specified wall temperature. For non uniform temperature distribution along the surface we have, Variable wall temperature boundary condition: T=T w (s) here T w (s) is the specified wall temperature variation as a function of distance along the surface (s). This boundary condition is very much suitable for high conductivity wall materials so as to keep the wall at known constant temperature variation. Joint initiative of IITs and IISc Funded by MHRD Page 10 of 39

11 However in case of insulators, where thermal conductivity is very low, the wall temperature usually remains unknown. In such cases, wall heat flux is treated as zero or wall is treated as adiabatic wall. The mathematical representation of this boundary condition is, Adiabatic wall boundary condition: q = k w T n = 0 Here q w is the wall heat flux. Moreover in some cased wall heat flux distribution can be apriorily known. Therefore q w or the wall heat transfer rate should be specified as the boundary condition. This wall heat flux is dependent on temperature gradient normal to the wall in the gas immediately above the wall.. Application to boundary layer flow Consider the boundary layer along a flat plate of length c as shown in Fig..1. V δ c Fig..1 Hypersonic flow over flat plate. A thin layer of fluid is assumed to be decelerated in the presence of the wall. This assumption leads to the mathematical expressionδ < c. Here δ is the local boundary layer thickness. Apart from this, for hypersonic flow, we can also assume that, that are; v<< u and << y Now consider the Continuity Equation in Non Dimensional form; ( ρu) ( ρv) + = 0 Joint initiative of IITs and IISc Funded by MHRD Page 11 of 39

12 Here u varies from 0 at the wall to 1 at the edge of the boundary layer. Therefore we can consider that u has order of magnitude equal to 1. It is mathematically represented as O(1). On similar lines, we can as well mention for density as, ρ = O(1). Actually the x-coordinate of all the points in the fluid domain vary from 0 to c which is length of the plate. Therefore the non-dimensional x length scale can be represented as x =O(1). However the y co-ordinate of all the points at a particular x- location varies from 0 toδ where δ is the local boundary layer thickness. Hence the non-dimensional length scale y is smaller magnitude in comparison with other length scales. This can be represented as y =O( δ / c ). For unit flat plate length, we have y =O(δ ). Therefore from the continuity equation in terms of order of magnitude is, [ 0(1) ][ 0(1) ] [ 0(1) ] v + = 0 0(1) 0( δ ) From the above equation it is clear that v must be of an order of magnitude equal to the local boundary layer thickness, δ, i.e. v =O(δ ). We know the non-dimensional form of X-momentum Eq. (0.6). Consider the order of magnitude form of each term as, u ρu = 0(1) x u ρv = 0(1) y p = 0(1) v µ = 0(1) u 1 µ = 0 δ Thus the order of magnitude equation for X momentum can be written as, (1) + 0(1) = 0(1) + 0(1) γm Re + δ Joint initiative of IITs and IISc Funded by MHRD Page 1 of 39

13 Lets assume that the Reynolds number is large. Therefore the term with Reynolds number in the denominator is of small magnitude which can be mathematically mentioned as, 1 Re = 0( ) δ Therefore the above equation now becomes; 1 1 0(1) + 0(1) = 0(1) + 0( δ ) 0(1) γm + δ It is clear from this figure that product of 0( δ )[0(1)] 0( δ ) = has very low order of magnitude in comparison with the rest of the terms in the same equation. This term actually is 1 Re v µ in X momentum equation (Eq. 0.6). Since this term is very small in magnitude we can neglect it. Therefore the non-dimensional X- momentum equation can now be written as, u u 1 p 1 ρu + ρv = + M Re x y γ x y u µ y (.1) The same in the dimensional form is, u u p u ρu + ρ v = + µ This is X-momentum is valid for high Reynolds number flows having thin boundary layer at the wall. Consider the Y momentum equation in non dimensional form given by Eq. (0.7); ρ v v 1 p 1 u v v u x ρ y γ M y Re x µ + = + + Lets do the order of magnitude analysis of this equation. 1 p 1 + = + + γm y δ 0( δ ) 0( δ ) 0( δ ) 0(1) Joint initiative of IITs and IISc Funded by MHRD Page 13 of 39

14 Here all the terms are of very small magnitude 0( δ ) except the term of pressure gradient if γ M is of unity order of magnitude. Hence that term should also have p very small order of magnitude. It implies that = 0( δ ). Therefore, from Y- y momentum equation given by Eq. 0.7 gets transformed for the boundary layer theory as; p = 0 y (.) This equation clearly states that pressure at a particular X location does not change with Y- coordinate such that the gradient of pressure remains zero in the boundary layer. Therefore the pressure is only function of X co-ordinate, p= px ( ) = pe( x), where pe( x ) is the pressure distribution outside the boundary layer. However if freestream Mach number, M, is very large so as to have 1/ M 0( ) γ = δ. In such cases, the pressure gradient in normal direction cab be large and still satisfy the Y- momentum equation. Hence for large Mach numbers, p/ y might be large enough to be expressed as 0(1). Hence pressure is not constant in the direction normal to the wall for hypersonic flows. Joint initiative of IITs and IISc Funded by MHRD Page 14 of 39

15 Lecture-3: Boundary layer equations 3.1 Boundary layer equations Let s derive the energy equation under the assumption of very thin boundary layer for very high Reynolds number hypersonic flows. We know that the non-dimensional energy equation is given by Eq This energy equation is for total energy which means the summation of kinetic energy and internal energy. We have neglected the potential energy of the fluid particle. Therefore let us derive the energy equation for kinetic energy alone. Consider the X-momentum and Y-momentum equations given by Eq. (0.) and (0.3). in the dimensional form. Multiplying the X momentum equation by u velocity, we get; Du D u p τ τ ρu = ρ = u + u + u Dt Dt ( / ) xx yx Similarly multiply Y momentum equation by v velocity, we get; Dv D v p ρv = ρ = v + v + v Dt Dt ( / ) τ xy τ yy Adding the above both the equations we have; Du ( / + v /) DV ( /) p p xx yx xy yy u v u τ τ ρ ρ v τ τ = = Dt Dt (3.1) But the energy equation given by Eq. (0.4), ( u v ) ( u v ) T T ρu + e ρv + e κ = + κ y ( pu) ( pv) ( uτ xx) ( uτ xy) ( vτ xy) ( vτ yy) Joint initiative of IITs and IISc Funded by MHRD Page 15 of 39

16 Replacing the L.H.S. in terms of the substantial derivative form we get, V V De ( + V / ) ρu e+ + ρv e+ = ρ = RHS... y Dt (3.) Subtracting Eq. (3.1) from (3.), we get, the heat energy equation is; D() e T T ( pu) ( pv) ρ = κ + κ + Dt ( u τ xx) ( u τ xy) ( v τ xy) ( v τ yy) p p xx yx xy yy u v u τ τ v τ τ Let s simplify the L.H.S., De () e e LHS... = ρ = ρu + ρv Dt From the definition of enthalpy we have, e= h p/ ρ Therefore the heat energy equation is, ( h p/ ρ) ( h p/ ρ) h h p/ ρ p/ ρ ρu + ρv = ρu + ρv ρu ρv, Further simplification of this equation leads to, h h p / ρ p / ρ h h ( pu) ( pv) ( ρu) ( ρv) ρu + ρv ρu ρv = ρu + ρv + pρ + Cancellation of similar terms from both the sides we have, L.H.S.= h h ρu + ρv Now, R.H.S. = T T ( uτ xx) ( uτ xy) ( vτ xy) ( vτ yy) κ + κ Joint initiative of IITs and IISc Funded by MHRD Page 16 of 39

17 p p τ xx τ yx τ xy τ yy + u + v u + v + Simplifying the viscous terms, ( u τ xx) ( u τ xy) ( v τ xy) ( v τ yy) xx yx xy yy u v u v u τ τ v τ τ = τ xx + τ yy + τ yx + u v u v u τ v xx τ yy τ yx λ µ u v u v = Therefore we have heat energy equation as, h h ρu + ρv T T = κ + κ p p + u + v u v λ µ u v u v p Now from Y momentum Boundary Layer Equation we have, = 0. Let us carry out the order of magnitude analysis for the above Energy Equation; 1 1 0(1) + 0(1) = 0(1) + + 0(1) + λ ( 0(1) + 0(1) ) + µ * 0(1) + * 0(1) + + δ 0( δ ) δ From the above equation it is very much clear that temperature gradient along Y axis has larger magnitude as compared to the temperature gradient along X axis. Moreover, we know from the momentum equation, that, v<< u & << for boundary layer assumption. Therefore gradients of lower magnitude, Energy Equation can be expressed as h h T dpe u ρu ρv κ u µ + = + + dx (3.3) Joint initiative of IITs and IISc Funded by MHRD Page 17 of 39

18 The mass, X-momentum, Y-momentum and energy equations for the boundary layer are respectively as follows, ρu ρv + = 0 u u p u ρu + ρ v = + µ p = 0 h h T dpe u ρu ρv κ u µ + = + + dx These equations are non-linear. However assumptions of boundary layer theory make solution procedure simpler. Apart from this the pressure is only function of X- coordinate hence can be represented using an ordinary differential equation rather than a partial differential equation. The variables u, v, p, ρ, T and h are the unknowns in these equations. However, p can be known from p=p e (x). Rest of the variables likeµ and κ are properties of fluid and are temperature dependant. The following perfect gas relations should also be used to complete the set of equations. p = ρrt h= CT p Boundary conditions to be considered to solve above equations are T At Wall: y=0, u=0, v=0, T=T w or = 0 n wall) w (adiabatic At boundary layer edge: y, u u e, T T e Here subscript e stands for the values measured or known at the edge of the boundary layer. The boundary layer equations are valid for compressible subsonic or supersonic flow. In case of application for hypersonic flows, it should again be noted that the Y- momentum equation p/ y = 0 should be changed for high Mach number cases. Joint initiative of IITs and IISc Funded by MHRD Page 18 of 39

19 Lecture-4: Similarity solution for boundary layer equation 4.1 Similarity solution of compressible boundary layer equation. Boundary layer theory makes it convenient to reduce the complexity of basic governing equation which has been derived using the order of magnitude estimate for the non-dimensional form of the governing equations. Hence following are the mass, momentum and energy equation which are to be solved for compressible boundary layer. ρu ρv + = 0 (4.1) u u p u ρu + ρ v = + µ (4.) p = 0 y (4.3) h h T dpe u ρu ρv κ u µ + = + + dx (4.4) These equations are derived for the X-Y coordinate system. Hence the variation of all the properties is assumed to be dependant on X and Y co-ordinates. However, lets transform the dependence of all the variables from X and Y to new dependant variables ( ξη, ). This transformation ensures the self similar solution for the velocity profile where u = u( η) and independent of ξ. This transformation has the following dependant variables as, x = eue edx (4.1) ξ ρ µ 0 ue v η = ρdy ξ (4.) 0 Joint initiative of IITs and IISc Funded by MHRD Page 19 of 39

20 Here ρ e, ue and µ e are the density, velocity and viscosity coefficients at the edge of the boundary layer and are functions of x only. Therefore ξ = ξ( x). These special variables chosen for transformation should be implemented for the governing equations derived especially for the boundary layer. This would lead to new form of the same governing equations. Following are the basic steps involved in this transformation. Step 1. Replacement of derivatives of independent variables. The new independent variables are expressed in terms of the old independent variables using Eq. (4.1) and (4.). Now we have to express the derivatives of them in terms of new independent variables. ξ η = + ξ η ξ η = + ξ η (4.3) (4.4) From the definition of ξ = ξ( x) given by Eq. (4.1) we can write, ξ = ρ eue µ e (4.5) ξ = 0 (4.6) From the definition of η given by Eq. (4.) we can write, η = ρ ξ η ue (4.7) Substituting Eqs. (4.5)-(4.7) into Eqs. (4.3) and (4.4), we can the derivatives as, η = + ρ eue µ e ξ x η (4.8) = ρ ξ η ue (4.9) Joint initiative of IITs and IISc Funded by MHRD Page 0 of 39

21 We know the definition of the stream function ψ defined as ψ = ρ u ψ = ρv (4.10) (4.11) The X momentum boundary layer equation given by Eq. (4.) in terms of ψ is, ψ u ψ u dpe u = + µ dx (4.1) Let s introduce the derivatives from eq. (4.8) and (4.9) in the eq. (4.1) we get, ueρ ψ u η u ψ η ψ ueρ u ρ eueµ e + ρ eueµ e + ξ η ξ x η ξ x η ξ η dpe ueρ ueρµ u = ρ eueµ e + d ξ ξ η ξ η (4.13) Multiplying Eq. (4.13) by ξ / ueρ, we get, ψ u η u ψ η ψ u ρ eueµ e ρ eueµ + e + η ξ x η ξ x η η ρ e dpe ueρµ u = ξ µ e + ρ dξ η ξ η (4.14) Step. Replacement of dependent variables. Let us define a function f of ξ and η, f ( ξη, ), such that u f = ue η f (4.15) Joint initiative of IITs and IISc Funded by MHRD Page 1 of 39

22 Here prime denotes the partial derivative of f with respect to η. We know that the velocity at the boundary layer edge is function of X alone. Hence it is function of ξ only ( ue= ue( ξ )). Therefore we can get the following u due f = f + ue ξ dξ ξ (4.16) u = η uf e (4.17) These two steps make it easy to express the governing equations of the thin boundary layer in terms of new variables so as to make it easier to solve. Joint initiative of IITs and IISc Funded by MHRD Page of 39

23 Lecture-5: Transformation of boundary layer equations 5.1 Transformation of boundary layer momentum equation. We have already seen the two steps involved in getting simpler boundary layer equation. The steps followed to those are seen here. Step 3. Expression of f in terms of ψ. The stream function ψ has been introduced to reduce the number of dependant variables of the governing equations. Therefore let us express the new dependent variable f ( ξη, ) in terms of stream function. From eq. (4.9), (4.10) and (4.15) we can write the following, ue ρ ψ = ρ u= ρ fu ξ η e ψ or = ξ f η (5.1) Integrating this equation with respect to η, we get, ψ = ξ f + F( ξ) (5.) Here F( ξ ) is any arbitrary function of ξ. However, from the definition of stream function, we know that difference in stream function casts mass flow rate. Therefore for the stream function should be anchored to zero at the non-blowing wall, ψξ (,0) = 0. Hence f = 0 and F( ξ ) = 0 ensure the zero value of ψ at the wall. This makes it clear that the any arbitrary function, represented by F( ξ ) must be zero, which leads to, ψ = ξ f (5.3) We will also have, ψ f 1 = ξ + f ξ ξ ξ (5.4) Joint initiative of IITs and IISc Funded by MHRD Page 3 of 39

24 Step 4. Derivation for final equation. We can derive the final expression using Eqs. (4.16)- (4.17) and (5.1) and substituting into Eq. (4.14) which is X-momentum boundary layer equation, we get, due f η ξ f ρ eueµ e f + ue + uf e dξ ξ f 1 η ρ eueµ e ξ + f + ξ f uf e ξ ξ ρ e dpe ue ρµ = ξ µ e + f (5.5) ρ dξ η ξ The Euler equation for the outer flow is, (5.6) Using Eq. (5.6) we can get the Eq. (5.5) as ue f ' η ξρ eueµ e( f ') ξ f ' ρ eueµ e ξµ eff + + ξ ξ = f ρ eueµ e η ξρ eue µ ef ff ξµ ef f ξ ξ Therefore, ue f ' e e e e e e ( ) ξρ u µ f ' + ξ f ' ρ u µ ξ ξ ( e) due ue ρ ρµ = ξ ue µ e + f ρ dξ η ξ (5.7) Joint initiative of IITs and IISc Funded by MHRD Page 4 of 39

25 It can be seen that the term involving η / does not appear in the Eq. (5.7). This the main reason for non evaluation of the same equation is, η / x earlier explicitly. Simplified form of 1 ( ) due f f ρ e 1 due f f f 1 ρµ + = + f ue dξ ξ ξ ρue dξ η ξρµ e e (5.8) Let s introduce a variable, C = ρµ / ρ eµ e and obtain the final form of the X- momentum equation in the transformed state. ξ ρ du f f Cf + ff = f + ξ f f u ρ dξ ξ ξ e e ( ) ( ) e (5.9) This transformed equation is valid for steady compressible flow in the thin boundary layer. The Y-momentum boundary layer y-momentum equation, gets transformed as, p = 0 η (5.10) 5. Transformation of boundary layer energy equation We can obtain the energy equation in the transformed variables using same strategy incorporated earlier. Let s substitute Eq. (4.8)- (4.11) in the L.H.S. of Eq. (4.4), we obtain the transformed R.H.S. term of energy equation as ψ h ψ h T dpe u κ u µ = + + dx (5.11) Using Eq.(4.8), we can express the first term on LHS as, ψ h ue ρ ψ h h ρ eu η = eµ e + y x ξ η ξ η (5.1) Joint initiative of IITs and IISc Funded by MHRD Page 5 of 39

26 Let s define the static enthalpy as the non dimensional variable, h g = g( ξη, ) = (5.13) he Hence we can write as, h g he = h e + g ξ ξ ξ (5.14) h = η hg e (5.15) Here g = g/ η, since he= he( x) hence its derivation with respect to η is zero. Therefore using Eq. (5.14), (5.15), we can re-express the first term of LHS of Eq. (5.1) as, ψ h ue ρ ψ h h ρ eu η = eµ e + y x ξ η ξ η ueρ g he η = ξ f ρ eueµ e he + g + hg e ξ ξ ξ g he η = ( ueρ f ) ρ eueµ ehe + ρ eueµ eg + heg ξ ξ ψ h g he η = ρρ e ue µ ehf e + ρρ e ue µ efg + ueρf hg e ξ ξ (5.16) Consider the second term of L.HS. of Eq.(5.1) using Eqs. (4.8) and (4.9) in the same way, we get, ψ h ue h eu ψ η ρ ψ ρ = eµ e + ξ η ξ η (5.17) Joint initiative of IITs and IISc Funded by MHRD Page 6 of 39

27 Now substituting Eq.(5.15),(5.1) and (5.4) in Eq. (5.17), we have ψ h ψ h f 1 η ueρ = ρ eueµ e ξ + f + ξ f hg e ξ ξ ξ f ρρ e ue µ ehefg η = ρρ e ueµ ehg e + + ρuhgf e e ξ ξ (5.18) Therefore complete L.H.S. of Eq.(5.11) using Eq. (5.16) and (5.18) is ψ h ψ h LHS... = g he η = ρρ e ue µ ehf e + ρρ e ue µ efg + ueρf hg e ξ ξ f ρρ e ue µ ehefg η ρρ e ueµ ehg e ρuhgf e e ξ ξ g he f hefg LHS... = ρρ e ue µ e hf e fg hg e + ξ ξ ξ ξ (5.19) Joint initiative of IITs and IISc Funded by MHRD Page 7 of 39

28 Lecture-6: Similarity solution of boundary layer equation 6.1 Transformation of boundary layer energy equation and similarity solution The simpler forms of X-momentum and Y-momentum have been obtained. The energy equation has also been obtained in the transformed form for LHS. Let s consider the RHS of Eq. (5.11) for the transformation. RHS of the Eq. (5.11) is as. T dpe u RHS... κ u µ = + + dx (6.1) Lets consider the first term of this RHS and express the temperature in terms of enthalpy as T = h/ Cp we get, T ueρ κueρ h κ = ξ η ξ η Cp (6.) If we consider the gas to be calorically perfect and variation of C p to be negligible then, we get T ueρ ρκue h ueρ ρκue κ = = hg e ξ η ξ η Cp ξ η ξcp T hu e eρ ρµ κ = g ξ η Pr (6.3) For the nd term of Eq.(6.1), we know, dpe dpe u = uf e ρ eueµ e dx d ξ dpe dpe due 3 due u = uf e ρ eueµ e ρ eueµ ef ρeue ρ eueµ ef dx dξ = = dξ dξ (6.4) Joint initiative of IITs and IISc Funded by MHRD Page 8 of 39

29 For the 3 rd term, from Eq.(6.1), we get 4 e e e u u ρ u u ρµ u ρµ µ = µ = = ξ ξ ξ ξ ( uf e ) ( f ) (6.5) Therefore the complete form of RHS of Eq. (6.1), 4 e e eρ ρµ 3 e eρ µ e e e T dp u h u du u RHS... = κ + u + µ = g ρ u µ f + f y dx y ξ η Pr dξ ξ (6.6) Therefore the complete form of Eq. (6.1) using Eq. (5.19) and (6.3)-(6.5) is ( ) Dividing both sides by ρρ u e e e µ, we get g he f hfg e he ρµ ρ euf e due ueρµ hf e + fg hg e = g + f ξ ξ ξ ξ ξρ eµ e η Pr ρ dξ ξρ eµ e ( ) Multiplying both sides by ξ / he, introducing the term C = ρµ / ρ eµ ethe above equation simplifies to, g ξ he f C ξρ euef due ue C ξ f + f g ξg fg = g + f ξ he ξ ξ Pr ρhe dξ he Finally rearranging the terms we get, ( ) C g f g he f ρ euef due ue C g + fg = ξ f g f Pr + + ξ he ξ ξ ρhe dξ he ( ) (6.7) Hence, the governing Eqs. (5.9), (5.10) and (6.7) form the transformed set of equations for compressible hypersonic boundary layer. The boundary conditions to be considered for solution are, η = 0, f = f ' = 0 and g = gw for isothermal wall boundary condition While, η = 0, f = f ' = 0 and g ' = 0 for adiabatic wall boundary condition Joint initiative of IITs and IISc Funded by MHRD Page 9 of 39

30 The boundary condition at the edge of the boundary layer is, η, f ' = 1 and g = 1 We can solve the transformed boundary layer equations which are partial differential equations. The necessary outcome for this is the velocity and enthalpy variation in the boundary layer. Since pressure is assumed to have same variation in the boundary layer as that at the outer flow, we can always evaluate the other thermodynamic properties in the boundary layer. Moreover the main parameters which we can evaluate are skin friction coefficient and Stanton Number. We know that the skin friction coefficient is defined as, C f = τw 1 ρeu e Where, the wall shear can be calculated as, u τ W = µ w Hence, C f u = µ ρ w eu e w Using the expression given by Eq. (4.9) we can write the equation for skin friction coefficient as, C f u ρ u = ρ e w µ w eu e ξ η w C f = ρ u e e u eρw µ w f ( ξ,0) ξ C f µρ w w = f ( ξ,0) (6.8) ρ ξ e Joint initiative of IITs and IISc Funded by MHRD Page 30 of 39

31 We can also calculate the Stanton number from the solution of boundary layer equations. Here Stanton number is defined as, S t qw = ρ u ( h h ) e e aw w Here h aw is the adiabatic wall enthalpy, which is the enthalpy at the wall when there is no heat transfer from the fluid to the wall. As well, h w is the wall static enthalpy corresponding to the temperature of the wall at isothermal wall temperature condition. Here, q W is the heat flux to the wall from the fluid which can be evaluated as, T = qw k y w or qw k h = cp w Assumption made here is the constancy of specific heat or calorically perfect gas. Thus using this expression for heatflux and from Eq. (4.9) and (5.13) we can calculate the Stanton number as, S t 1 k h = ρeue( haw hw ) cp w S S t t 1 k ueρhe h = ρeue( haw hw ) cp ξ η 1 k ueρhe = g'( ξ, 0) ρ u ( h h ) c ξ e e aw w p w (6.9) Joint initiative of IITs and IISc Funded by MHRD Page 31 of 39

32 Lecture-7: Hypersonic flow over flat plate 7.1 Hypersonic flow over flat plate. The equations (5.9)(5.10) and (6.7) form the boundary layer equations. Let s consider these equations for solution of hypersonic flow over flat plate. As we had seen earlier, there are two prominent boundary conditions considered for this flow over flat plate as, Isothermal wall, T w =const T Or an adiabatic wall, = 0 w All the freestream variables of the variables at the edge of the boundary layer are now assumed be constant. Hence The variables ue, he, ρ e are of constant values and independent of ξ and η. Therefore the governing equations reduce to, f f ( Cf ) + ff = ξ f f ξ ξ (7.1) C g f ue C g fg ξ f g + = f Pr ξ ξ he ( ) (7.) These equations are the partial differential equations. Let s assume that the functions f and g are functions of Eqs.(4.40) and (4.41) are still partial differential equations. Let us assume that f and g are the functions of η only. Hence f and g are independent of ξ. For these assumptions the same equations reduce to, ( Cf ) + ff = 0 (7.3) C ue C g + fg + ( f ) = 0 Pr he (7.4) These equations are single independent variable equations, hence are the nonlinear ordinary differential equations. These equations are valid for a compressible boundary layer over a flat plate with constant wall conditions. Here both the constants Joint initiative of IITs and IISc Funded by MHRD Page 3 of 39

33 C = ρµ / ρ eµ e and Pr = µ cp / κ are meant for the local values in the boundary layer. We can you shooting technique for solving these equations. Thus obtained velocity and thermal boundary layers can be used to obtain the wall heat flux and shear stress in turn the Stanton number and skin friction coefficient. Here the ratio of Stanton number and skin friction coefficient can be approximated as the function of Prandtl s number by Reynolds analogy. 7. Hypersonic flow around a stagnation point Hypersonic flow around a blunt body forms a stagnation point which can also be evaluated using boundary layer equations. Consider hypersonic flow around a blunt body which marks a stagnation region, as sketched in Fig Consider the flow to be D flow for simplicity; hence the span of the cylinder is infinity. Let X be the direction of freestream flow and R be the radius of curvature at the surface. Fig. 7.1 Hypersonic flow around the stagnation point region [1]. Joint initiative of IITs and IISc Funded by MHRD Page 33 of 39

34 Let us consider that f and g from Eq. (5.9) and (6.7) are functions of η alone, Therefore, = = = 0 This leads to the following changed in Eq. (5.9) and (6.7) as, (CCff ) + ffff = ξξ [(ff ) ρρ ee ] dduu ee uu ee ρρ dddd (7.5) And CC PPPP gg + ffff = ξξ[ ρρ eeuu ee ρρh ee ff dduu ee dddd ] CC uu ee h ee (ff ) (7.6) These equations are still ξξ dependent. Moreover, we can assume that the velocity at the edge of the boundary layer, uu ee is very small and static enthalpy at the edge of the boundary layer is h ee = h 0 (stagnation enthalpy). These facts lead to the assumption that, uu ee h ee 0 (7.7) We can as well assume that, the flow velocity is the boundary layer at edge of the boundary layer at the stagnation point behind the normal shock is low as to be considered in the incompressible flow regime. Hence, we can use the result of incompressible and inviscid flow at the stagnation point, which expresses the boundary layer velocity as, uu ee = dduu ee xx (7.8) dddd ss Here dduu ee is the velocity gradient at the stagnation point external to the boundary dddd ss layer. Using these assumptions we can re-express the ξξ as, xx ξξ = ρρ ee uu ee μμ ee dddd = ρρ ee μμ ee dduu ee xxxxxx 0 0 dddd ss Or xx ξξ = ρρ ee μμ ee dduu ee xx dddd ss (7.9) Joint initiative of IITs and IISc Funded by MHRD Page 34 of 39

35 The velocity gradient required for calculation here can be evaluated as, dduu ee dddd = dduu ee dddd dddd dddd = (dduu ee dddd ) (dddd dddd ) (7.10) But from the definition of ξξ, Eq. (4.1) we know that, dddd dddd = ρρ eeuu ee μμ ee (7.11) Substituting this Eq. (7.11) into (7.10), we get dduu ee = 1 dduu ee dddd ρρ ee uu ee μμ ee dddd (7.1) From Eq. (7.11) into (7.1), we write, dduu ee dddd ss = 1 ρρ ee μμ ee xx Or, dduu ee dddd = ss 1 μμ ee ρρ ee (dduu ee dddd) ss xx dduu ee dddd ss Let s consider the term (ξξ uu ee ) dduu ee dddd which appears in the Eq We can rewrite this term and can derive for the same as, ξξ uu ee dduu ee ) ss xx ] = [ρρ eeμμ ee (dduu ee dddd dddd (dduu ee dddd ) ss xx 1 ρρ ee μμ ee xx = 1 Similarly consider the term (ξξ uu ee )(ρρ ee uu ee ρρh ee ) dduu ee dddd appearing in Eq. (7.). We can re-arrange this term as well, ξξ ρρ eeuu ee dduu ee = ρρ ee ρρ ρρh ee dddd ρρh ee uu ee dduu ee xx ee dddd ss dduu ee xx 1 = ρρ ee dduu ee dddd ss ρρ ee μμ ee xx ρρh ee dddd ss xx However at the stagnation point, x=0. This fact leads to, ξξ ρρ eeuu ee ρρh ee dduu ee dddd = 0 Moreover, for a calorically perfect gas, we have, ρρ ee = pp ee TT ee = pp ee h = h gg ρρ pp TT pp h ee h ee Joint initiative of IITs and IISc Funded by MHRD Page 35 of 39

36 Using all these short expression for various terms involved in Eq.(7.1) and (7.) we can get these equations as, (CCCC ) + ffff = (ff ) gg (7.13) CC PPPP gg + ffgg = 0 (7.14) These equations are special equation for stagnation point flow. These equations are independent of ξξ. We can use numerical techniques shooting technique to solve these equations. Joint initiative of IITs and IISc Funded by MHRD Page 36 of 39

37 Lecture-8: Stagnation point flow field 8.1 Hypersonic flow around a stagnation point. The known correlation for stagnation point heat flux for cylinder is Cylinder: (8.1) This equation is valid for D configurations. Moreover for, sphere or axi-symmetric configuration the stagnation point heat flux can be predicted using, Sphere: 0.763PPPP 0.6 (ρρ ee μμ ee ) 1 dduu ee (h dddd aaaa h ww ) (8.) These equations are called as Fay and Riddle equations for stagnation point heat flux prediction. The stagnation point heat flux is more for sphere in comparison for that of the cylinder of same diameter. The main reason for this discrimination is the dimensionality of the flow. The hypersonic flow is two dimensional for flow over cylinder hence it has two possible direction for passing over the cylinder however the flow over the sphere is three dimensional. Due to the extra available dimension, flow passes around the object easily hence the shock stand off distance and the boundary layer thickness decrease for the sphere in comparison with the cylinder of same diameter. The decreased boundary layer thickness increases the gradient and hence the shear stress and heat flux at the wall for sphere in comparison with the cylinder. The closer observation to the Eq. (8.1) and (8.) suggests that the wall heat flux is propotional to the square root of the stream wise velocity gradient,( dduu ee ) along the dddd stagnation streamline. dp = ρ u du e e e e Hence, du 1 dp e e = (8.3) dx ρ u dx e e Joint initiative of IITs and IISc Funded by MHRD Page 37 of 39

38 We can evaluate the pressure gradient of this equation from the known pressure variation given by Newtonian method. C p = sin θ Here θθ is flow deflection angle, the angle between a tangent at any point on the surface and the freestream direction. Let s define as the angle between freestream velocity and the local surface normal. Hence the pressure distribution gets transformed to, C p = cos ϕ From definition of pressure coefficient we have, Or, pp ee pp qq = cccccc cos p = q ϕ + p e Differentiating this equation w.r.t X we get, dpe 4q = dϕ cosϕsin ϕ dx ρ u dx e e (8.4) Substituting the Eq. (8.4) in to the Eq. (8.3) we get, dduu ee dddd = 4qq cccccc ssssss dd uu ee ρρ ee dddd (8.5) All the terms involved in this equation can be evaluated using various approximations. uu ee = dduu ee xx (8.6) dddd ss Joint initiative of IITs and IISc Funded by MHRD Page 38 of 39

39 cccccc 1 ssssss xx RR dd dddd = 1 RR (8.7) qq = 1 (pp ee pp ) (8.8) Now, substituting Eqs. (8.6) (8.8) in (8.5), we get, Or dduu ee dddd = (pp ee pp ) xx ρρ ee xx RR 1 RR dduu ee dddd = 1 RR (pp ee pp ) ρρ (8.9) This is the approximate expression for velocity gradient encountered in the Eq. (8.1) or (8.). The expression clearly suggests that the wall heat flux is inversely proportional to the nose radius or radius at the stagnation point. This is the main reason for having hypersonic vehicles being blunt nosed to reduce the heat load at the compromise of the drag. Joint initiative of IITs and IISc Funded by MHRD Page 39 of 39

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 15-Convective Heat Transfer Fausto Arpino f.arpino@unicas.it Introduction In conduction problems the convection entered the analysis

More information

Introduction to Aerospace Engineering

Introduction to Aerospace Engineering Introduction to Aerospace Engineering Lecture slides Challenge the future 3-0-0 Introduction to Aerospace Engineering Aerodynamics 5 & 6 Prof. H. Bijl ir. N. Timmer Delft University of Technology 5. Compressibility

More information

Governing Equations of Fluid Dynamics

Governing Equations of Fluid Dynamics Chapter 3 Governing Equations of Fluid Dynamics The starting point of any numerical simulation are the governing equations of the physics of the problem to be solved. In this chapter, we first present

More information

Multiscale Hydrodynamic Phenomena

Multiscale Hydrodynamic Phenomena M2, Fluid mechanics 2014/2015 Friday, December 5th, 2014 Multiscale Hydrodynamic Phenomena Part I. : 90 minutes, NO documents 1. Quick Questions In few words : 1.1 What is dominant balance? 1.2 What is

More information

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t) IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common

More information

Module 2: Governing Equations and Hypersonic Relations

Module 2: Governing Equations and Hypersonic Relations Module 2: Governing Equations and Hypersonic Relations Lecture -2: Mass Conservation Equation 2.1 The Differential Equation for mass conservation: Let consider an infinitely small elemental control volume

More information

CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION

CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,

More information

UNIT IV BOUNDARY LAYER AND FLOW THROUGH PIPES Definition of boundary layer Thickness and classification Displacement and momentum thickness Development of laminar and turbulent flows in circular pipes

More information

Viscous flow along a wall

Viscous flow along a wall Chapter 8 Viscous flow along a wall 8. The no-slip condition All liquids and gases are viscous and, as a consequence, a fluid near a solid boundary sticks to the boundary. The tendency for a liquid or

More information

To study the motion of a perfect gas, the conservation equations of continuity

To study the motion of a perfect gas, the conservation equations of continuity Chapter 1 Ideal Gas Flow The Navier-Stokes equations To study the motion of a perfect gas, the conservation equations of continuity ρ + (ρ v = 0, (1.1 momentum ρ D v Dt = p+ τ +ρ f m, (1.2 and energy ρ

More information

CHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW

CHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW CHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW 4.1 Introduction Boundary layer concept (Prandtl 1904): Eliminate selected terms in the governing equations Two key questions (1) What are the

More information

2. Getting Ready for Computational Aerodynamics: Fluid Mechanics Foundations

2. Getting Ready for Computational Aerodynamics: Fluid Mechanics Foundations . Getting Ready for Computational Aerodynamics: Fluid Mechanics Foundations We need to review the governing equations of fluid mechanics before examining the methods of computational aerodynamics in detail.

More information

6.2 Governing Equations for Natural Convection

6.2 Governing Equations for Natural Convection 6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed

More information

Chapter 9 Flow over Immersed Bodies

Chapter 9 Flow over Immersed Bodies 57:00 Mechanics of Fluids and Transport Processes Chapter 9 Professor Fred Stern Fall 009 1 Chapter 9 Flow over Immersed Bodies Fluid flows are broadly categorized: 1. Internal flows such as ducts/pipes,

More information

PEMP ACD2505. M.S. Ramaiah School of Advanced Studies, Bengaluru

PEMP ACD2505. M.S. Ramaiah School of Advanced Studies, Bengaluru Governing Equations of Fluid Flow Session delivered by: M. Sivapragasam 1 Session Objectives -- At the end of this session the delegate would have understood The principle of conservation laws Different

More information

Chapter 6: Incompressible Inviscid Flow

Chapter 6: Incompressible Inviscid Flow Chapter 6: Incompressible Inviscid Flow 6-1 Introduction 6-2 Nondimensionalization of the NSE 6-3 Creeping Flow 6-4 Inviscid Regions of Flow 6-5 Irrotational Flow Approximation 6-6 Elementary Planar Irrotational

More information

Conservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion

Conservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation

More information

Compressible Duct Flow with Friction

Compressible Duct Flow with Friction Compressible Duct Flow with Friction We treat only the effect of friction, neglecting area change and heat transfer. The basic assumptions are 1. Steady one-dimensional adiabatic flow 2. Perfect gas with

More information

FLUID MECHANICS. Chapter 9 Flow over Immersed Bodies

FLUID MECHANICS. Chapter 9 Flow over Immersed Bodies FLUID MECHANICS Chapter 9 Flow over Immersed Bodies CHAP 9. FLOW OVER IMMERSED BODIES CONTENTS 9.1 General External Flow Characteristics 9.3 Drag 9.4 Lift 9.1 General External Flow Characteristics 9.1.1

More information

ENGR Heat Transfer II

ENGR Heat Transfer II ENGR 7901 - Heat Transfer II Convective Heat Transfer 1 Introduction In this portion of the course we will examine convection heat transfer principles. We are now interested in how to predict the value

More information

Several forms of the equations of motion

Several forms of the equations of motion Chapter 6 Several forms of the equations of motion 6.1 The Navier-Stokes equations Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,

More information

Introduction to Turbomachinery

Introduction to Turbomachinery 1. Coordinate System Introduction to Turbomachinery Since there are stationary and rotating blades in turbomachines, they tend to form a cylindrical form, represented in three directions; 1. Axial 2. Radial

More information

Lecture1: Characteristics of Hypersonic Atmosphere

Lecture1: Characteristics of Hypersonic Atmosphere Module 1: Hypersonic Atmosphere Lecture1: Characteristics of Hypersonic Atmosphere 1.1 Introduction Hypersonic flight has special traits, some of which are seen in every hypersonic flight. Presence of

More information

Applied Gas Dynamics Flow With Friction and Heat Transfer

Applied Gas Dynamics Flow With Friction and Heat Transfer Applied Gas Dynamics Flow With Friction and Heat Transfer Ethirajan Rathakrishnan Applied Gas Dynamics, John Wiley & Sons (Asia) Pte Ltd c 2010 Ethirajan Rathakrishnan 1 / 121 Introduction So far, we have

More information

Review for Exam Hyunse Yoon, Ph.D. Adjunct Assistant Professor Department of Mechanical Engineering, University of Iowa

Review for Exam Hyunse Yoon, Ph.D. Adjunct Assistant Professor Department of Mechanical Engineering, University of Iowa Review for Exam2 11. 13. 2015 Hyunse Yoon, Ph.D. Adjunct Assistant Professor Department of Mechanical Engineering, University of Iowa Assistant Research Scientist IIHR-Hydroscience & Engineering, University

More information

3. FORMS OF GOVERNING EQUATIONS IN CFD

3. FORMS OF GOVERNING EQUATIONS IN CFD 3. FORMS OF GOVERNING EQUATIONS IN CFD 3.1. Governing and model equations in CFD Fluid flows are governed by the Navier-Stokes equations (N-S), which simpler, inviscid, form is the Euler equations. For

More information

Chapter 9 Flow over Immersed Bodies

Chapter 9 Flow over Immersed Bodies 57:00 Mechanics of Fluids and Transport Processes Chapter 9 Professor Fred Stern Fall 014 1 Chapter 9 Flow over Immersed Bodies Fluid flows are broadly categorized: 1. Internal flows such as ducts/pipes,

More information

Fundamental Concepts of Convection : Flow and Thermal Considerations. Chapter Six and Appendix D Sections 6.1 through 6.8 and D.1 through D.

Fundamental Concepts of Convection : Flow and Thermal Considerations. Chapter Six and Appendix D Sections 6.1 through 6.8 and D.1 through D. Fundamental Concepts of Convection : Flow and Thermal Considerations Chapter Six and Appendix D Sections 6.1 through 6.8 and D.1 through D.3 6.1 Boundary Layers: Physical Features Velocity Boundary Layer

More information

Introduction to Aerodynamics. Dr. Guven Aerospace Engineer (P.hD)

Introduction to Aerodynamics. Dr. Guven Aerospace Engineer (P.hD) Introduction to Aerodynamics Dr. Guven Aerospace Engineer (P.hD) Aerodynamic Forces All aerodynamic forces are generated wither through pressure distribution or a shear stress distribution on a body. The

More information

Getting started: CFD notation

Getting started: CFD notation PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =

More information

AE/ME 339. K. M. Isaac Professor of Aerospace Engineering. 12/21/01 topic7_ns_equations 1

AE/ME 339. K. M. Isaac Professor of Aerospace Engineering. 12/21/01 topic7_ns_equations 1 AE/ME 339 Professor of Aerospace Engineering 12/21/01 topic7_ns_equations 1 Continuity equation Governing equation summary Non-conservation form D Dt. V 0.(2.29) Conservation form ( V ) 0...(2.33) t 12/21/01

More information

Shock and Expansion Waves

Shock and Expansion Waves Chapter For the solution of the Euler equations to represent adequately a given large-reynolds-number flow, we need to consider in general the existence of discontinuity surfaces, across which the fluid

More information

Fundamentals of Fluid Dynamics: Elementary Viscous Flow

Fundamentals of Fluid Dynamics: Elementary Viscous Flow Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research

More information

Introduction to Fluid Mechanics

Introduction to Fluid Mechanics Introduction to Fluid Mechanics Tien-Tsan Shieh April 16, 2009 What is a Fluid? The key distinction between a fluid and a solid lies in the mode of resistance to change of shape. The fluid, unlike the

More information

Steady waves in compressible flow

Steady waves in compressible flow Chapter Steady waves in compressible flow. Oblique shock waves Figure. shows an oblique shock wave produced when a supersonic flow is deflected by an angle. Figure.: Flow geometry near a plane oblique

More information

Fundamentals of compressible and viscous flow analysis - Part II

Fundamentals of compressible and viscous flow analysis - Part II Fundamentals of compressible and viscous flow analysis - Part II Lectures 3, 4, 5 Instantaneous and averaged temperature contours in a shock-boundary layer interaction. Taken from (Pasquariello et al.,

More information

ρ Du i Dt = p x i together with the continuity equation = 0, x i

ρ Du i Dt = p x i together with the continuity equation = 0, x i 1 DIMENSIONAL ANALYSIS AND SCALING Observation 1: Consider the flow past a sphere: U a y x ρ, µ Figure 1: Flow past a sphere. Far away from the sphere of radius a, the fluid has a uniform velocity, u =

More information

Chapter 5. The Differential Forms of the Fundamental Laws

Chapter 5. The Differential Forms of the Fundamental Laws Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations

More information

Chapter 9: Differential Analysis

Chapter 9: Differential Analysis 9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control

More information

AA210A Fundamentals of Compressible Flow. Chapter 5 -The conservation equations

AA210A Fundamentals of Compressible Flow. Chapter 5 -The conservation equations AA210A Fundamentals of Compressible Flow Chapter 5 -The conservation equations 1 5.1 Leibniz rule for differentiation of integrals Differentiation under the integral sign. According to the fundamental

More information

Masters in Mechanical Engineering. Problems of incompressible viscous flow. 2µ dx y(y h)+ U h y 0 < y < h,

Masters in Mechanical Engineering. Problems of incompressible viscous flow. 2µ dx y(y h)+ U h y 0 < y < h, Masters in Mechanical Engineering Problems of incompressible viscous flow 1. Consider the laminar Couette flow between two infinite flat plates (lower plate (y = 0) with no velocity and top plate (y =

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

58:160 Intermediate Fluid Mechanics Bluff Body Professor Fred Stern Fall 2014

58:160 Intermediate Fluid Mechanics Bluff Body Professor Fred Stern Fall 2014 Professor Fred Stern Fall 04 Chapter 7 Bluff Body Fluid flows are broadly categorized:. Internal flows such as ducts/pipes, turbomachinery, open channel/river, which are bounded by walls or fluid interfaces:

More information

AE/ME 339. Computational Fluid Dynamics (CFD) K. M. Isaac. Momentum equation. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.

AE/ME 339. Computational Fluid Dynamics (CFD) K. M. Isaac. Momentum equation. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept. AE/ME 339 Computational Fluid Dynamics (CFD) 9//005 Topic7_NS_ F0 1 Momentum equation 9//005 Topic7_NS_ F0 1 Consider the moving fluid element model shown in Figure.b Basis is Newton s nd Law which says

More information

Chapter 9: Differential Analysis of Fluid Flow

Chapter 9: Differential Analysis of Fluid Flow of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known

More information

6. Laminar and turbulent boundary layers

6. Laminar and turbulent boundary layers 6. Laminar and turbulent boundary layers John Richard Thome 8 avril 2008 John Richard Thome (LTCM - SGM - EPFL) Heat transfer - Convection 8 avril 2008 1 / 34 6.1 Some introductory ideas Figure 6.1 A boundary

More information

Differential relations for fluid flow

Differential relations for fluid flow Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow

More information

Principles of Convection

Principles of Convection Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid

More information

Variable Normalization (nondimensionalization and scaling) for Navier-Stokes equations: a practical guide

Variable Normalization (nondimensionalization and scaling) for Navier-Stokes equations: a practical guide TR-CFD-13-77 Variable Normalization (nondimensionalization and scaling) for Navier-Stokes equations: a practical guide Marc Montagnac Technical Report V. 1.0 October 013 Copyrighted by the author(s) Centre

More information

UNIVERSITY OF EAST ANGLIA

UNIVERSITY OF EAST ANGLIA UNIVERSITY OF EAST ANGLIA School of Mathematics May/June UG Examination 2007 2008 FLUIDS DYNAMICS WITH ADVANCED TOPICS Time allowed: 3 hours Attempt question ONE and FOUR other questions. Candidates must

More information

1. Fluid properties and concept of continuum

1. Fluid properties and concept of continuum 1. Fluid properties and concept of continuum Fluid: a substance that deforms continuously under an applied shear stress (e.g., air, water, upper mantle...) Not practical/possible to treat fluid mechanics

More information

AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS

AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29

More information

Advanced Heat and Mass Transfer by Amir Faghri, Yuwen Zhang, and John R. Howell

Advanced Heat and Mass Transfer by Amir Faghri, Yuwen Zhang, and John R. Howell Laminar external natural convection on vertical and horizontal flat plates, over horizontal and vertical cylinders and sphere, as well as plumes, wakes and other types of free flow will be discussed in

More information

Introduction to Gas Dynamics All Lecture Slides

Introduction to Gas Dynamics All Lecture Slides Introduction to Gas Dynamics All Lecture Slides Teknillinen Korkeakoulu / Helsinki University of Technology Autumn 009 1 Compressible flow Zeroth law of thermodynamics 3 First law of thermodynamics 4 Equation

More information

n v molecules will pass per unit time through the area from left to

n v molecules will pass per unit time through the area from left to 3 iscosity and Heat Conduction in Gas Dynamics Equations of One-Dimensional Gas Flow The dissipative processes - viscosity (internal friction) and heat conduction - are connected with existence of molecular

More information

Section 2: Lecture 1 Integral Form of the Conservation Equations for Compressible Flow

Section 2: Lecture 1 Integral Form of the Conservation Equations for Compressible Flow Section 2: Lecture 1 Integral Form of the Conservation Equations for Compressible Flow Anderson: Chapter 2 pp. 41-54 1 Equation of State: Section 1 Review p = R g T " > R g = R u M w - R u = 8314.4126

More information

Review for Exam Hyunse Yoon, Ph.D. Assistant Research Scientist IIHR-Hydroscience & Engineering University of Iowa

Review for Exam Hyunse Yoon, Ph.D. Assistant Research Scientist IIHR-Hydroscience & Engineering University of Iowa 57:020 Fluids Mechanics Fall2013 1 Review for Exam3 12. 11. 2013 Hyunse Yoon, Ph.D. Assistant Research Scientist IIHR-Hydroscience & Engineering University of Iowa 57:020 Fluids Mechanics Fall2013 2 Chapter

More information

Rocket Thermodynamics

Rocket Thermodynamics Rocket Thermodynamics PROFESSOR CHRIS CHATWIN LECTURE FOR SATELLITE AND SPACE SYSTEMS MSC UNIVERSITY OF SUSSEX SCHOOL OF ENGINEERING & INFORMATICS 25 TH APRIL 2017 Thermodynamics of Chemical Rockets ΣForce

More information

Effective Boundary Conditions for Continuum Method of Investigation of Rarefied Gas Flow over Blunt Body

Effective Boundary Conditions for Continuum Method of Investigation of Rarefied Gas Flow over Blunt Body Effective Boundary Conditions for Continuum Method of Investigation of Rarefied Gas Flow over Blunt Body I.G. Brykina a, B.V. Rogov b, I.L. Semenov c, and G.A. Tirskiy a a Institute of Mechanics, Moscow

More information

2. FLUID-FLOW EQUATIONS SPRING 2019

2. FLUID-FLOW EQUATIONS SPRING 2019 2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid

More information

AA210A Fundamentals of Compressible Flow. Chapter 1 - Introduction to fluid flow

AA210A Fundamentals of Compressible Flow. Chapter 1 - Introduction to fluid flow AA210A Fundamentals of Compressible Flow Chapter 1 - Introduction to fluid flow 1 1.2 Conservation of mass Mass flux in the x-direction [ ρu ] = M L 3 L T = M L 2 T Momentum per unit volume Mass per unit

More information

Detailed Derivation of Fanno Flow Relationships

Detailed Derivation of Fanno Flow Relationships Detailed Derivation of Fanno Flow Relationships Matthew MacLean, Ph.D. Version. October 03 v. fixed a sign error on Eq. 5 ! My motivation for writing this note results from preparing course notes for the

More information

6.1 Momentum Equation for Frictionless Flow: Euler s Equation The equations of motion for frictionless flow, called Euler s

6.1 Momentum Equation for Frictionless Flow: Euler s Equation The equations of motion for frictionless flow, called Euler s Chapter 6 INCOMPRESSIBLE INVISCID FLOW All real fluids possess viscosity. However in many flow cases it is reasonable to neglect the effects of viscosity. It is useful to investigate the dynamics of an

More information

Summary of Dimensionless Numbers of Fluid Mechanics and Heat Transfer

Summary of Dimensionless Numbers of Fluid Mechanics and Heat Transfer 1. Nusselt number Summary of Dimensionless Numbers of Fluid Mechanics and Heat Transfer Average Nusselt number: convective heat transfer Nu L = conductive heat transfer = hl where L is the characteristic

More information

Lesson 6 Review of fundamentals: Fluid flow

Lesson 6 Review of fundamentals: Fluid flow Lesson 6 Review of fundamentals: Fluid flow The specific objective of this lesson is to conduct a brief review of the fundamentals of fluid flow and present: A general equation for conservation of mass

More information

Gasdynamics 1-D compressible, inviscid, stationary, adiabatic flows

Gasdynamics 1-D compressible, inviscid, stationary, adiabatic flows Gasdynamics 1-D compressible, inviscid, stationary, adiabatic flows 1st law of thermodynamics ρ const Kontrollfläche 1 2 m u 2 u 1 z Q 12 +P 12 = ṁ } h 2 h {{} 1 Enthalpy Q 12 + 1 2 (u2 2 u2 1 }{{} ) +

More information

Problem 4.3. Problem 4.4

Problem 4.3. Problem 4.4 Problem 4.3 Problem 4.4 Problem 4.5 Problem 4.6 Problem 4.7 This is forced convection flow over a streamlined body. Viscous (velocity) boundary layer approximations can be made if the Reynolds number Re

More information

Notes 4: Differential Form of the Conservation Equations

Notes 4: Differential Form of the Conservation Equations Low Speed Aerodynamics Notes 4: Differential Form of the Conservation Equations Deriving Conservation Equations From the Laws of Physics Physical Laws Fluids, being matter, must obey the laws of Physics.

More information

CONVECTIVE HEAT TRANSFER

CONVECTIVE HEAT TRANSFER CONVECTIVE HEAT TRANSFER Mohammad Goharkhah Department of Mechanical Engineering, Sahand Unversity of Technology, Tabriz, Iran CHAPTER 4 HEAT TRANSFER IN CHANNEL FLOW BASIC CONCEPTS BASIC CONCEPTS Laminar

More information

Entry Aerodynamics MARYLAND U N I V E R S I T Y O F. Entry Aerodynamics. ENAE Launch and Entry Vehicle Design

Entry Aerodynamics MARYLAND U N I V E R S I T Y O F. Entry Aerodynamics. ENAE Launch and Entry Vehicle Design Atmospheric Regimes on Entry Basic fluid parameters Definition of Mean Free Path Rarified gas Newtonian flow Continuum Newtonian flow (hypersonics) 2014 David L. Akin - All rights reserved http://spacecraft.ssl.umd.edu

More information

Fluid Mechanics II Viscosity and shear stresses

Fluid Mechanics II Viscosity and shear stresses Fluid Mechanics II Viscosity and shear stresses Shear stresses in a Newtonian fluid A fluid at rest can not resist shearing forces. Under the action of such forces it deforms continuously, however small

More information

INDIAN INSTITUTE OF TECHNOOGY, KHARAGPUR Date: -- AN No. of Students: 5 Sub. No.: ME64/ME64 Time: Hours Full Marks: 6 Mid Autumn Semester Examination Sub. Name: Convective Heat and Mass Transfer Instructions:

More information

UNIVERSITY OF OSLO. Faculty of Mathematics and Natural Sciences. MEK4300/9300 Viscous flow og turbulence

UNIVERSITY OF OSLO. Faculty of Mathematics and Natural Sciences. MEK4300/9300 Viscous flow og turbulence UNIVERSITY OF OSLO Faculty of Mathematics and Natural Sciences Examination in: Day of examination: Friday 15. June 212 Examination hours: 9. 13. This problem set consists of 5 pages. Appendices: Permitted

More information

CONVECTIVE HEAT TRANSFER

CONVECTIVE HEAT TRANSFER CONVECTIVE HEAT TRANSFER Mohammad Goharkhah Department of Mechanical Engineering, Sahand Unversity of Technology, Tabriz, Iran CHAPTER 3 LAMINAR BOUNDARY LAYER FLOW LAMINAR BOUNDARY LAYER FLOW Boundary

More information

Quick Recapitulation of Fluid Mechanics

Quick Recapitulation of Fluid Mechanics Quick Recapitulation of Fluid Mechanics Amey Joshi 07-Feb-018 1 Equations of ideal fluids onsider a volume element of a fluid of density ρ. If there are no sources or sinks in, the mass in it will change

More information

Flight Vehicle Terminology

Flight Vehicle Terminology Flight Vehicle Terminology 1.0 Axes Systems There are 3 axes systems which can be used in Aeronautics, Aerodynamics & Flight Mechanics: Ground Axes G(x 0, y 0, z 0 ) Body Axes G(x, y, z) Aerodynamic Axes

More information

Interpreting Differential Equations of Transport Phenomena

Interpreting Differential Equations of Transport Phenomena Interpreting Differential Equations of Transport Phenomena There are a number of techniques generally useful in interpreting and simplifying the mathematical description of physical problems. Here we introduce

More information

Thermodynamics ENGR360-MEP112 LECTURE 7

Thermodynamics ENGR360-MEP112 LECTURE 7 Thermodynamics ENGR360-MEP11 LECTURE 7 Thermodynamics ENGR360/MEP11 Objectives: 1. Conservation of mass principle.. Conservation of energy principle applied to control volumes (first law of thermodynamics).

More information

In this section, mathematical description of the motion of fluid elements moving in a flow field is

In this section, mathematical description of the motion of fluid elements moving in a flow field is Jun. 05, 015 Chapter 6. Differential Analysis of Fluid Flow 6.1 Fluid Element Kinematics In this section, mathematical description of the motion of fluid elements moving in a flow field is given. A small

More information

7.2 Sublimation. The following assumptions are made in order to solve the problem: Sublimation Over a Flat Plate in a Parallel Flow

7.2 Sublimation. The following assumptions are made in order to solve the problem: Sublimation Over a Flat Plate in a Parallel Flow 7..1 Sublimation Over a Flat Plate in a Parallel Flow The following assumptions are made in order to solve the problem: 1.. 3. The flat plate is very thin and so the thermal resistance along the flat plate

More information

Boundary-Layer Theory

Boundary-Layer Theory Hermann Schlichting Klaus Gersten Boundary-Layer Theory With contributions from Egon Krause and Herbert Oertel Jr. Translated by Katherine Mayes 8th Revised and Enlarged Edition With 287 Figures and 22

More information

Objectives. Conservation of mass principle: Mass Equation The Bernoulli equation Conservation of energy principle: Energy equation

Objectives. Conservation of mass principle: Mass Equation The Bernoulli equation Conservation of energy principle: Energy equation Objectives Conservation of mass principle: Mass Equation The Bernoulli equation Conservation of energy principle: Energy equation Conservation of Mass Conservation of Mass Mass, like energy, is a conserved

More information

MIXED CONVECTION SLIP FLOW WITH TEMPERATURE JUMP ALONG A MOVING PLATE IN PRESENCE OF FREE STREAM

MIXED CONVECTION SLIP FLOW WITH TEMPERATURE JUMP ALONG A MOVING PLATE IN PRESENCE OF FREE STREAM THERMAL SCIENCE, Year 015, Vol. 19, No. 1, pp. 119-18 119 MIXED CONVECTION SLIP FLOW WITH TEMPERATURE JUMP ALONG A MOVING PLATE IN PRESENCE OF FREE STREAM by Gurminder SINGH *a and Oluwole Daniel MAKINDE

More information

The Bernoulli Equation

The Bernoulli Equation The Bernoulli Equation The most used and the most abused equation in fluid mechanics. Newton s Second Law: F = ma In general, most real flows are 3-D, unsteady (x, y, z, t; r,θ, z, t; etc) Let consider

More information

Fundamentals of Aerodynamics

Fundamentals of Aerodynamics Fundamentals of Aerodynamics Fourth Edition John D. Anderson, Jr. Curator of Aerodynamics National Air and Space Museum Smithsonian Institution and Professor Emeritus University of Maryland Me Graw Hill

More information

High Speed Aerodynamics. Copyright 2009 Narayanan Komerath

High Speed Aerodynamics. Copyright 2009 Narayanan Komerath Welcome to High Speed Aerodynamics 1 Lift, drag and pitching moment? Linearized Potential Flow Transformations Compressible Boundary Layer WHAT IS HIGH SPEED AERODYNAMICS? Airfoil section? Thin airfoil

More information

Thin airfoil theory. Chapter Compressible potential flow The full potential equation

Thin airfoil theory. Chapter Compressible potential flow The full potential equation hapter 4 Thin airfoil theory 4. ompressible potential flow 4.. The full potential equation In compressible flow, both the lift and drag of a thin airfoil can be determined to a reasonable level of accuracy

More information

Notes #4a MAE 533, Fluid Mechanics

Notes #4a MAE 533, Fluid Mechanics Notes #4a MAE 533, Fluid Mechanics S. H. Lam lam@princeton.edu http://www.princeton.edu/ lam October 23, 1998 1 The One-dimensional Continuity Equation The one-dimensional steady flow continuity equation

More information

Review for Exam Hyunse Yoon, Ph.D. Adjunct Assistant Professor Department of Mechanical Engineering, University of Iowa

Review for Exam Hyunse Yoon, Ph.D. Adjunct Assistant Professor Department of Mechanical Engineering, University of Iowa Review for Exam3 12. 9. 2015 Hyunse Yoon, Ph.D. Adjunct Assistant Professor Department of Mechanical Engineering, University of Iowa Assistant Research Scientist IIHR-Hydroscience & Engineering, University

More information

ENGR Heat Transfer II

ENGR Heat Transfer II ENGR 7901 - Heat Transfer II External Flows 1 Introduction In this chapter we will consider several fundamental flows, namely: the flat plate, the cylinder, the sphere, several other body shapes, and banks

More information

Aerothermodynamics of High Speed Flows

Aerothermodynamics of High Speed Flows Aerothermodynamics of High Speed Flows Lecture 1: Introduction G. Dimitriadis 1 The sound barrier Supersonic aerodynamics and aircraft design go hand in hand Aspects of supersonic flow theory were developed

More information

APPLICATION OF THE DEFECT FORMULATION TO THE INCOMPRESSIBLE TURBULENT BOUNDARY LAYER

APPLICATION OF THE DEFECT FORMULATION TO THE INCOMPRESSIBLE TURBULENT BOUNDARY LAYER APPLICATION OF THE DEFECT FORMULATION TO THE INCOMPRESSIBLE TURBULENT BOUNDARY LAYER O. ROUZAUD ONERA OAt1 29 avenue de la Division Leclerc - B.P. 72 92322 CHATILLON Cedex - France AND B. AUPOIX AND J.-PH.

More information

Chapter 6 Fundamental Concepts of Convection

Chapter 6 Fundamental Concepts of Convection Chapter 6 Fundamental Concepts of Convection 6.1 The Convection Boundary Layers Velocity boundary layer: τ surface shear stress: s = μ u local friction coeff.: C f y y=0 τ s ρu / (6.) (6.1) Thermal boundary

More information

UNIVERSITY of LIMERICK

UNIVERSITY of LIMERICK UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH Faculty of Science and Engineering END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MA4607 SEMESTER: Autumn 2012-13 MODULE TITLE: Introduction to Fluids DURATION OF

More information

FORCED CONVECTION FILM CONDENSATION OF DOWNWARD-FLOWING VAPOR ON HORIZONTAL TUBE WITH WALL SUCTION EFFECT

FORCED CONVECTION FILM CONDENSATION OF DOWNWARD-FLOWING VAPOR ON HORIZONTAL TUBE WITH WALL SUCTION EFFECT Journal of Marine Science and Technology, Vol., No. 5, pp. 5-57 () 5 DOI:.69/JMST--5- FORCED CONVECTION FILM CONDENSATION OF DOWNWARD-FLOWING VAPOR ON HORIZONTAL TUBE WITH WALL SUCTION EFFECT Tong-Bou

More information

Lecture 3: 1. Lecture 3.

Lecture 3: 1. Lecture 3. Lecture 3: 1 Lecture 3. Lecture 3: 2 Plan for today Summary of the key points of the last lecture. Review of vector and tensor products : the dot product (or inner product ) and the cross product (or vector

More information

Transport processes. 7. Semester Chemical Engineering Civil Engineering

Transport processes. 7. Semester Chemical Engineering Civil Engineering Transport processes 7. Semester Chemical Engineering Civil Engineering 1. Elementary Fluid Dynamics 2. Fluid Kinematics 3. Finite Control Volume Analysis 4. Differential Analysis of Fluid Flow 5. Viscous

More information

Aerodynamics. Lecture 1: Introduction - Equations of Motion G. Dimitriadis

Aerodynamics. Lecture 1: Introduction - Equations of Motion G. Dimitriadis Aerodynamics Lecture 1: Introduction - Equations of Motion G. Dimitriadis Definition Aerodynamics is the science that analyses the flow of air around solid bodies The basis of aerodynamics is fluid dynamics

More information

2 Law of conservation of energy

2 Law of conservation of energy 1 Newtonian viscous Fluid 1 Newtonian fluid For a Newtonian we already have shown that σ ij = pδ ij + λd k,k δ ij + 2µD ij where λ and µ are called viscosity coefficient. For a fluid under rigid body motion

More information

Notes #6 MAE 533, Fluid Mechanics

Notes #6 MAE 533, Fluid Mechanics Notes #6 MAE 533, Fluid Mechanics S. H. Lam lam@princeton.edu http://www.princeton.edu/ lam October 1, 1998 1 Different Ways of Representing T The speed of sound, a, is formally defined as ( p/ ρ) s. It

More information