Quick viewing(Text Mode)

Fluid Mechanics, Euler and Bernoulli Equations

Fluid Mechanics, Euler and Bernoulli Equations

FLUID , EULER AND BERNOULLI EQUATIONS

© M. Ragheb 3/30/2019

INTRODUCTION

th The early part of the 18 century saw the burgeoning of the field of theoretical mechanics pioneered by Leonhard Euler and the father and son Johann and . We introduce the equations of continuity and conservation of of fluid flow, from which we derive the Euler and Bernoulli equations. The Bernoulli equation is the most famous equation in . Its significance is that when the increases in a fluid stream, the decreases, and when the velocity decreases, the pressure increases. The Bernoulli equation is applied to the airfoil of a wind machine rotor, defining the , drag and thrust coefficients and the pitching .

Figure 1. Johannes Bernoulli, the father.

THE CONSERVATION OR

The continuity equation of fluid mechanics expresses the notion that mass cannot be created nor destroyed or that mass is conserved. It relates the flow field variables at a of the flow in terms of the fluid density and the fluid velocity vector, and is given by:

  .(V )  0 (1) t

We consider the vector identity resembling the of differentiation:

.(VVV )    .  .   (2) where the operator is noted to act on a vector quantity, and the operator acts on a scalar quantity. This allows us to rewrite the continuity equation as:

 VV.    .  0 (3) t

SUBSTANTIAL

We can use the substantial derivative:

D   (.)V  Convetive (4) Local Derivative Dt t Derivative where the partial derivative is called the local derivative and the term is called the convective derivative. In terms of the substantial derivative the continuity equation can be expressed as::

D  .0V  (5) Dt

MOMENTUM CONSERVATION OR EQUATION OF

Newton’s second law is frequently written in terms of an and a vectors as:

F ma (6)

A more general form describes the force vector as the rate of change of the momentum vector as:

d F () mV (7) dt

Its general form is written in term of volume and a surface over an arbitrary control volume v as:

()V dv ( VdSV . )    pd v  fd v  Fd v (8)      viscous vt S v v v

where the velocity vector is:

V uxˆˆ  vy  wzˆ (9)

The cartesian coordinates x, y and z components of the continuity equation are:

()up dv ( V . dS ) u   d v  f d v  ( F ) d v    x  x viscous vtxS v v v ()vp dv ( V . dS ) v   d v  f d v  ( F ) d v (10)    y  y viscous vtySv v v ()wp dv ( V . dS ) w   d v  f d v  ( F ) dv     z z viscous vvtzSvv

In this equation the product:

(.)V dS (11) is a scalar and has no components.

NAVIER STOKES EQUATIONS

By using the divergence or Gauss’s theorem the can be turned into a :

(V . dS ) u  (  uV ). dS    .(  uV ) d v SS v (VdSv . )  (  vVdS ).    .(  vVd ) v (12) SS v (V . dS ) w  (  wV ). dS    .(  wV ) d v SS v

The volume integrals over an arbitrary volume now yield:

()up  .(uV )    f  ( F ) txx x viscous ()vp  .(vV )    f  ( F ) (13) tyy y viscous ()wp  .(wV )    f  ( F ) tzz z viscous

These are known as the Navier-Stokes equations. They apply to the unsteady,

three dimensional flow of any fluid, compressible or incompressible, viscous or inviscid. In terms of the substantial derivative, the Navier-Stokes equations can be expressed as:

Du p   fF  () Dt x x x viscous Dv p   fF  () (14) Dt y y y viscous Dw p   fF  () Dt z z z viscous

EULER EQUATIONS

For a steady state flow the time partial vanish. For the viscous terms are equal to zero. In the absence of body the f , f , anf f terms x y z disappear. The Euler equations result as:

p .(uV )   x p .(vV )   (15) y p .(wV )   z

INVISCID

For an inviscid flow without body forces, the momentum conservation equations of fluid mechanics are:

Du p   Dt x Dv p   (16) Dt y Dw p   Dt z

These equations can also be written as:

()up Vu.()    tx ()vp Vv.()    (17) ty ()wp Vw.()    tz

For steady flow the partial time derivative vanishes, and we can write:

p Vu.()   x p Vv.()   (18) y p Vw.()   z

Expanding the gradient term, we get:

u  u  u  p u  v   w   x  y  z  x v  v  v  p u  v   w   (19) x  y  z  y w  w  w  p u  v   w   x  y  z  z

Rearranging, we get:

u  u  u1  p u v  w   x  y  z  x v  v  v1  p u v  w   (20) x  y  z  y w  w  w1  p u v  w   x  y  z  z

STREAMLINES DIFFERENTIAL EQUATIONS

The definition of a streamline in a flow is that it is parallel to the velocity vector. Hence the cross product of the directed of the streamline and the velocity vector is zero:

ds V 0 (21)

where: ds dx xˆˆ  dy y  dz zˆ

V u xˆˆ  v y  w zˆ

The cross product can be expanded in the form of a determinant as:

xˆˆ y zˆ ds V dx dy dz u v w xwdyˆˆ()()()  vdz  yudz  wdx  zvdxˆ  udy (22)  0

The vector being equal to zero, its components must be equal to zero yielding the differential equations for the streamline f(x,y,z) =0, as:

wdy v dz 0 u dz wdx 0 (23) v dx u dy 0

EULER’S EQUATION

Multiplying the flow equations respectively by dx, dy and dz, we get:

u  u  u1  p u dx v dx  w dx   dx x  y  z  x v  v  v1  p u dy v dy  w dy   dy (24) x  y  z  y w  w  w1  p u dz v dz  w dz   dz x  y  z  z

Using the streamline differential equations, we can write:

u  u  u1  p u dx u dy  w dz   dx x  y  z  x v  v  v1  p u dx v dy  w dz   dy (25) x  y  z  y w  w  w1  p u dx v dy  w dz   dz x  y  z  z

The differentials of functions u = u(x,y,z), v = v(x,y,z), w =w (x,y,z) are:

u  u  u du dx  dy  dz x  y  z v  v  v dv dx  dy  dz (26) x  y  z w  w  w dw dx  dy  dz x  y  z

This allows us to write:

1 p udu dx  x 1 p vdv dy (27)  y 1 p wdw dz  z

Through integration we can write:

11p d() u2  dx 2  x 11p d() v2  dy (28) 2  y 11p d() w2  dz 2  z

Adding the three last equations we get:

11p  p  p d()() u2 v 2  w 2   dx  dy  dz 2  x  y  z (29) 11 d() V2  dp 2 

From the last equation we can write a simple form of Euler’s equation as:

dp VdV (30)

Euler’s equation applies to an inviscid flow with no body forces. It relates the change in velocity along a streamline dV to the change in pressure dp along the same streamline.

BERNOULLI EQUATION, INCOMPRESSIBLE FLOW

Considering the case of incompressible flow, we can use integration to yield:

pV22 dp VdV pV11

 22 p2 p 1   V 2  V 1  (31) 2 11 p V22  p  V  constant 122 1 2 2

The relation between pressure and velocity in an inviscid incompressible flow was enunciated in the form of Bernoulli’s equation, first presented by Euler:

1 p V2 constant (32) 2

This equation is the most famous equation in fluid mechanics. Its significance is that when the velocity increases, the pressure decreases, and when the velocity decreases, the pressure increases. The of the terms in the equation are kinetic per unit volume. Even though it was derived from the momentum conservation equation, it is also a relation for the mechanical energy in an incompressible flow. It states that the done on a fluid by the pressure forces is equal to the change of of the flow. In fact it can be derived from the energy conservation equation of fluid flow. The fact that Bernoulli’s equation can be interpreted as Newton’s second law or an energy equation illustrates that the energy equation is redundant for the analysis of inviscid, incompressible flow.

ROTOR AIRFOIL

The sharp end of an airfoil shape at point B is designated as the trailing edge. The leading edge is the locos of the point A of the nose of the airfoil profile that is the farthest from the trailing edge. The distance L = AB is known as the chord of the profile, with AMB being the upper surface and ANB the lower surface. At any distance along the chord AB from the nose, points may be identified half way between the upper and lower surfaces whose locus is usually curved is called the camber or median line of the airfoil section. The incidence angle i is the angle between the chord and the air vector V at infinite upstream. The zero lift angle is the angle φ between the chord and the zero lift 0 line. The lift angle is the angle φ between the zero lift line and the air speed vector V at infinite upstream. It is conventional to take φ and i as positive and φ as negative. The 0

following relations hold:

ii ()      00 (33) i  ()  00    

Figure 2. Rotor airfoil geometry.

Another rotor blade parameter is the maximum thickness h, which is sometimes expressed as a fraction of the chord length and is called the thickness/chord ratio or relative thickness. The relative thickness can range from 3-20 percent, with the common wind machine rotors values covering the range of 10-15 percent. The spot along the chord where the maximum thickness occurs in airfoils covers a range of 20 - 60 percent of the chord from the leading edge, and in wind machines rotors this is around 30 percent.

FORCES ON A MOVING ROTOR IN A STILL

We can assume that the airfoil is at rest and the air moving at the same speed but in the opposite direction. In this case the aerodynamic force exerted on the rotor will not change in magnitude. The resultant force will depend only on the relative speed and the angle of attack. To simplify the analysis, the airfoil is taken at rest in a moving stream of air in an infinite upstream speed V. The pressure of the air on the external surface of the airfoil will not be uniform due to the effect of the Bernoulli force. On the upper surface there results a lower pressure, and on the lower surface an increase in the pressure. We can represent the pressure variation on the rotor surface by considering a line to the airfoil profile surface whose magnitude is Kp , given by:

pp K  0 (34) p 1 V 2 2 where: p is the static pressureat thebaseof the perpendicular tothe surface p is the pressureat theinfiniteupstream 0 V is the speed at theinfiniteupstream  is thedensity at theinfiniteupstream

The value of Kp is negative for the points at the upper surface of the airfoil, and positive for the lower surface.

Figure 3. Pressure profile on airfoil segment.

The resultant thrust force FT is inclined with respect to the relative speed direction and is given by:

1 F  C AV 2 (35) TT2 where: A isthearea equal tothe product of thechord by thelengthof therotor

CT isthetotal aeorodynamiccoefficient

The force F has two components. The first component is parallel to the velocity vector or the drag force:

1 F  C AV 2 (36) DD2

The second component is perpendicular to the velocity vector, or the lift force:

1 F  C AV 2 (37) LL2

These forces are perpendicular and we can apply the leading to the thrust force as:

2 2 2 FFFTDL (38)

Consequently:

22 FFFTDL (39)

In addition:

2 2 2 12   1 2   1 2  CTLD AV    C AV    C AV  2   2   2  (40) 2 2 2 CCCTDL where C is the thrust coefficient. T

Figure 4. Angle of attack and angle for an aircraft.

Figure 5. Vortices formation behind an airfoil at a steep angle of attack.

Figure 6. Flow separation and loss of lift at a steep angle of attack.

Figure *7. Generation and loss of lift on an airfoil.

Figure 8. Effect of angle of attack on lift and flow separation.

AERODYNAMIC

If M is the aerodynamic moment of the force F relative to the leading edge, we can define a pitching moment coefficient C from the expression: m

1 M  C ALV 2 (41) 2 m where: L isthechord length.

The aerodynamic forces on the rotor may be represented by a lift, a drag, and a pitching moment. At each value of the incidence angle there exists a particular point C about which the pitching moment of the aerodynamic force F is zero. This unique point is called the center of pressure. The aerodynamic effects on the airfoil section can be represented by the lift and the drag alone acting at that point. The position of the center of pressure relative to the leading edge is calculated from the ratio:

AC C CP m (42) AB CL

It is usually in the range of: 25 – 30 percent.

REFERENCES

1. Désiré Le Gouriérès, “Wind Plants, Theory and Design,” Pergamon Press, 1982. rd 2. John D. Anderson, Jr., “Fundamentals of ,” 3 edition, McGrawHill, 2001.

EXERCISE

1. A wind rotor airfoil is placed in the air flow at sea level conditions with a free stream 3 velocity of 10 m/s. The density at standard sea level conditions is 1.23 kg/m and the 5 2 pressure is 1.01 x 10 Newtons/m . At a point along the rotor airfoil the pressure is 0.90 5 2 x 10 Newtons/m . By applying Bernoulli’s equation estimate the velocity at this point. 1 2. The lift force on a rotor blade is given by: F  C AV 2 . The drag force is given by: LL2 1 F  C AV 2 . Derive expressions for: DD2 a) The thrust force FT, b) The thrust coefficient CT.