Misplaced Pages

Pressure coefficient

Article snapshot taken from Wikipedia with creative commons attribution-sharealike license. Give it a read and then ask your questions in the chat. We can research this topic together.
(Redirected from Pressure distribution) Dimensionless number describing relative pressures in a fluid flow field

In fluid dynamics, the pressure coefficient is a dimensionless number which describes the relative pressures throughout a flow field. The pressure coefficient is used in aerodynamics and hydrodynamics. Every point in a fluid flow field has its own unique pressure coefficient, Cp.

In many situations in aerodynamics and hydrodynamics, the pressure coefficient at a point near a body is independent of body size. Consequently, an engineering model can be tested in a wind tunnel or water tunnel, pressure coefficients can be determined at critical locations around the model, and these pressure coefficients can be used with confidence to predict the fluid pressure at those critical locations around a full-size aircraft or boat.

Definition

The pressure coefficient is a parameter for studying both incompressible/compressible fluids such as water and air. The relationship between the dimensionless coefficient and the dimensional numbers is

C p = p p 1 2 ρ V 2 {\displaystyle C_{p}={p-p_{\infty } \over {\frac {1}{2}}\rho _{\infty }V_{\infty }^{2}}}

where:

p {\displaystyle p} is the static pressure at the point at which pressure coefficient is being evaluated
p {\displaystyle p_{\infty }} is the static pressure in the freestream (i.e. remote from any disturbance)
ρ {\displaystyle \rho _{\infty }} is the freestream fluid density (Air at sea level and 15 °C is 1.225 k g / m 3 {\displaystyle {\rm {kg/m^{3}}}} )
V {\displaystyle V_{\infty }} is the freestream velocity of the fluid, or the velocity of the body through the fluid

Incompressible flow

Main article: Incompressible flow

Using Bernoulli's equation, the pressure coefficient can be further simplified for potential flows (inviscid, and steady):

C p | M 0 = p p p 0 p = 1 ( u u ) 2 {\displaystyle C_{p}|_{M\,\approx \,0}={p-p_{\infty } \over p_{0}-p_{\infty }}={1-{\bigg (}{\frac {u}{u_{\infty }}}{\bigg )}^{2}}}

where:

u {\displaystyle u} is the flow speed at the point at which pressure coefficient is being evaluated
M {\displaystyle M} is the Mach number, which is taken in the limit of zero
p 0 {\displaystyle p_{0}} is the flow's stagnation pressure

This relationship is valid for the flow of incompressible fluids where variations in speed and pressure are sufficiently small that variations in fluid density can be neglected. This assumption is commonly made in engineering practice when the Mach number is less than about 0.3.

  • C p {\displaystyle C_{p}} of zero indicates the pressure is the same as the freestream pressure.
  • C p {\displaystyle C_{p}} of one corresponds to the stagnation pressure and indicates a stagnation point.
  • the most negative values of C p {\displaystyle C_{p}} in a liquid flow can be summed to the cavitation number to give the cavitation margin. If this margin is positive, the flow is locally fully liquid, while if it is zero or negative the flow is cavitating or gas.

Locations where C p = 1 {\displaystyle C_{p}=-1} are significant in the design of gliders because this indicates a suitable location for a "Total energy" port for supply of signal pressure to the Variometer, a special Vertical Speed Indicator which reacts to vertical movements of the atmosphere but does not react to vertical maneuvering of the glider.

In an incompressible fluid flow field around a body, there will be points having positive pressure coefficients up to one, and negative pressure coefficients including coefficients less than minus one.

Compressible flow

Main article: Compressible flow

In the flow of compressible fluids such as air, and particularly the high-speed flow of compressible fluids, 1 2 ρ v 2 {\displaystyle {{\frac {1}{2}}\rho v^{2}}} (the dynamic pressure) is no longer an accurate measure of the difference between stagnation pressure and static pressure. Also, the familiar relationship that stagnation pressure is equal to total pressure does not always hold true. (It is always true in isentropic flow, but the presence of shock waves can cause the flow to depart from isentropic.) As a result, pressure coefficients can be greater than one in compressible flow.

Perturbation theory

The pressure coefficient C p {\displaystyle C_{p}} can be estimated for irrotational and isentropic flow by introducing the potential Φ {\displaystyle \Phi } and the perturbation potential ϕ {\displaystyle \phi } , normalized by the free-stream velocity u {\displaystyle u_{\infty }}

Φ = u x + ϕ ( x , y , z ) {\displaystyle \Phi =u_{\infty }x+\phi (x,y,z)}

Using Bernoulli's equation,

Φ t + Φ Φ 2 + γ γ 1 p ρ = constant {\displaystyle {\frac {\partial \Phi }{\partial t}}+{\frac {\nabla \Phi \cdot \nabla \Phi }{2}}+{\frac {\gamma }{\gamma -1}}{\frac {p}{\rho }}={\text{constant}}}

which can be rewritten as

Φ t + Φ Φ 2 + a 2 γ 1 = constant {\displaystyle {\frac {\partial \Phi }{\partial t}}+{\frac {\nabla \Phi \cdot \nabla \Phi }{2}}+{\frac {a^{2}}{\gamma -1}}={\text{constant}}}

where a {\displaystyle a} is the sound speed.

The pressure coefficient becomes

C p = p p γ 2 p M 2 = 2 γ M 2 [ ( a a ) 2 γ γ 1 1 ] = 2 γ M 2 [ ( γ 1 a 2 ( u 2 2 Φ t Φ Φ 2 ) + 1 ) γ γ 1 1 ] 2 γ M 2 [ ( 1 γ 1 a 2 ( ϕ t + u ϕ x ) ) γ γ 1 1 ] 2 ϕ t u 2 2 ϕ x u {\displaystyle {\begin{aligned}C_{p}&={\frac {p-p_{\infty }}{{\frac {\gamma }{2}}p_{\infty }M^{2}}}={\frac {2}{\gamma M^{2}}}\left\\&={\frac {2}{\gamma M^{2}}}\left\\&\approx {\frac {2}{\gamma M^{2}}}\left\\&\approx -{\frac {2\phi _{t}}{u_{\infty }^{2}}}-{\frac {2\phi _{x}}{u_{\infty }}}\end{aligned}}}

where a {\displaystyle a_{\infty }} is the far-field sound speed.

Local piston theory

The classical piston theory is a powerful aerodynamic tool. From the use of the momentum equation and the assumption of isentropic perturbations, one obtains the following basic piston theory formula for the surface pressure:

p = p ( 1 + γ 1 2 w a ) 2 γ γ 1 {\displaystyle p=p_{\infty }\left(1+{\frac {\gamma -1}{2}}{\frac {w}{a}}\right)^{\frac {2\gamma }{\gamma -1}}}

where w {\displaystyle w} is the downwash speed and a {\displaystyle a} is the sound speed.

C p = p p γ 2 p M 2 = 2 γ M 2 [ ( 1 + γ 1 2 w a ) 2 γ γ 1 1 ] {\displaystyle C_{p}={\frac {p-p_{\infty }}{{\frac {\gamma }{2}}p_{\infty }M^{2}}}={\frac {2}{\gamma M^{2}}}\left}

The surface is defined as

F ( x , y , z , t ) = z f ( x , y , t ) = 0 {\displaystyle F(x,y,z,t)=z-f(x,y,t)=0}

The slip velocity boundary condition leads to

F | F | ( u + ϕ x , ϕ y , ϕ z ) = V wall F | F | = F t 1 | F | {\displaystyle {\frac {\nabla F}{|\nabla F|}}(u_{\infty }+\phi _{x},\phi _{y},\phi _{z})=V_{\text{wall}}\cdot {\frac {\nabla F}{|\nabla F|}}=-{\frac {\partial F}{\partial t}}{\frac {1}{|\nabla F|}}}

The downwash speed w {\displaystyle w} is approximated as

w = f t + u f x {\displaystyle w={\frac {\partial f}{\partial t}}+u_{\infty }{\frac {\partial f}{\partial x}}}

Hypersonic flow

Main article: Hypersonic speed

In hypersonic flow, the pressure coefficient can be accurately calculated for a vehicle using Newton's corpuscular theory of fluid motion, which is inaccurate for low-speed flow and relies on three assumptions:

  1. The flow can be modeled as a stream of particles in rectilinear motion
  2. Upon impact with a surface, all normal momentum is lost
  3. All tangential momentum is conserved, and flow follows the body

For a freestream velocity V {\displaystyle V_{\infty }} impacting a surface of area A {\displaystyle A} , which is inclined at an angle θ {\displaystyle \theta } relative to the freestream, the change in normal momentum is V sin θ {\displaystyle V_{\infty }\sin \theta } and the mass flux incident on the surface is ρ V A sin θ {\displaystyle \rho _{\infty }V_{\infty }A\sin \theta } , with ρ {\displaystyle \rho _{\infty }} being the freestream air density. Then the momentum flux, equal to the force exerted on the surface F {\displaystyle F} , from Newton's second law is equal to:

F = ( ρ V A sin θ ) ( V sin θ ) = ρ V 2 A sin 2 θ {\displaystyle F=(\rho _{\infty }V_{\infty }A\sin \theta )(V_{\infty }\sin \theta )=\rho _{\infty }V_{\infty }^{2}A\sin ^{2}\theta }

Dividing by the surface area, it is clear that the force per unit area is equal to the pressure difference between the surface pressure p {\displaystyle p} and the freestream pressure p {\displaystyle p_{\infty }} , leading to the relation:

F A = p p = ρ V 2 sin 2 θ p p 1 2 ρ V 2 = 2 sin 2 θ {\displaystyle {\frac {F}{A}}=p-p_{\infty }=\rho _{\infty }V_{\infty }^{2}\sin ^{2}\theta \implies {\frac {p-p_{\infty }}{{\frac {1}{2}}\rho _{\infty }V_{\infty }^{2}}}=2\sin ^{2}\theta }

The last equation may be identified as the pressure coefficient, meaning that Newtonian theory predicts that the pressure coefficient in hypersonic flow is:

C p = 2 sin 2 θ {\displaystyle C_{p}=2\sin ^{2}\theta }

For very high speed flows, and vehicles with sharp surfaces, the Newtonian theory works very well.

Modified Newtonian law

A modification to the Newtonian theory, specifically for blunt bodies, was proposed by Lester Lees:

C p = C p , max sin 2 θ {\displaystyle C_{p}=C_{p,\max }\sin ^{2}\theta }

where C p , max {\displaystyle C_{p,\max }} is the maximum value of the pressure coefficient at the stagnation point behind a normal shock wave:

C p , max = p o p 1 2 ρ V 2 = p 1 2 ρ V 2 ( p o p 1 ) = 2 γ M 2 ( p o p 1 ) {\displaystyle C_{p,\max }={\frac {p_{o}-p_{\infty }}{{\frac {1}{2}}\rho _{\infty }V_{\infty }^{2}}}={\frac {p_{\infty }}{{\frac {1}{2}}\rho _{\infty }V_{\infty }^{2}}}\left({\frac {p_{o}}{p_{\infty }}}-1\right)={\frac {2}{\gamma M_{\infty }^{2}}}\left({\frac {p_{o}}{p_{\infty }}}-1\right)}

where p o {\displaystyle p_{o}} is the stagnation pressure and γ {\displaystyle \gamma } is the ratio of specific heats. The last relation is obtained from the ideal gas law p = ρ R T {\displaystyle p=\rho RT} , Mach number M = V / a {\displaystyle M=V/a} , and speed of sound a = γ R T {\displaystyle a={\sqrt {\gamma RT}}} . The Rayleigh pitot tube formula for a calorically perfect normal shock says that the ratio of the stagnation and freestream pressure is:

p o p = [ ( γ + 1 ) 2 M 2 4 γ M 2 2 ( γ 1 ) ] γ / ( γ 1 ) [ γ ( 2 M 2 1 ) + 1 γ + 1 ] {\displaystyle {\frac {p_{o}}{p_{\infty }}}=\left^{\gamma /(\gamma -1)}\left}

Therefore, it follows that the maximum pressure coefficient for the Modified Newtonian law is:

C p , max = 2 γ M 2 { [ ( γ + 1 ) 2 M 2 4 γ M 2 2 ( γ 1 ) ] γ / ( γ 1 ) [ γ ( 2 M 2 1 ) + 1 γ + 1 ] 1 } {\displaystyle C_{p,\max }={\frac {2}{\gamma M_{\infty }^{2}}}\left\{\left^{\gamma /(\gamma -1)}\left-1\right\}}

In the limit when M {\displaystyle M_{\infty }\rightarrow \infty } , the maximum pressure coefficient becomes:

C p , max = [ ( γ + 1 ) 2 4 γ ] γ / ( γ 1 ) ( 4 γ + 1 ) {\displaystyle C_{p,\max }=\left^{\gamma /(\gamma -1)}\left({\frac {4}{\gamma +1}}\right)}

And as γ 1 {\displaystyle \gamma \rightarrow 1} , C p , max = 2 {\displaystyle C_{p,\max }=2} , recovering the pressure coefficient from Newtonian theory at very high speeds. The modified Newtonian theory is substantially more accurate than the Newtonian model for calculating the pressure distribution over blunt bodies.

Pressure distribution

An airfoil at a given angle of attack will have what is called a pressure distribution. This pressure distribution is simply the pressure at all points around an airfoil. Typically, graphs of these distributions are drawn so that negative numbers are higher on the graph, as the C p {\displaystyle C_{p}} for the upper surface of the airfoil will usually be farther below zero and will hence be the top line on the graph.

Relationship with aerodynamic coefficients

All the three aerodynamic coefficients are integrals of the pressure coefficient curve along the chord. The coefficient of lift for a two-dimensional airfoil section with strictly horizontal surfaces can be calculated from the coefficient of pressure distribution by integration, or calculating the area between the lines on the distribution. This expression is not suitable for direct numeric integration using the panel method of lift approximation, as it does not take into account the direction of pressure-induced lift. This equation is true only for zero angle of attack.

C l = 1 x T E x L E x L E x T E ( C p l ( x ) C p u ( x ) ) d x {\displaystyle C_{l}={\frac {1}{x_{TE}-x_{LE}}}\int \limits _{x_{LE}}^{x_{TE}}\left(C_{p_{l}}(x)-C_{p_{u}}(x)\right)dx}

where:

C p l {\displaystyle C_{p_{l}}} is pressure coefficient on the lower surface
C p u {\displaystyle C_{p_{u}}} is pressure coefficient on the upper surface
x L E {\displaystyle x_{LE}} is the leading edge location
x T E {\displaystyle x_{TE}} is the trailing edge location

When the lower surface C p {\displaystyle C_{p}} is higher (more negative) on the distribution it counts as a negative area as this will be producing down force rather than lift.

See also

References

  1. L. J. Clancy (1975) Aerodynamics, § 3.6, Pitman Publishing Limited, London. ISBN 0-273-01120-0
  2. Abbott and Von Doenhoff, Theory of Wing Sections, equation 2.24
  3. Anderson, John D. Fundamentals of Aerodynamics. 4th ed. New York: McGraw Hill, 2007. 219.
  4. https://thesis.library.caltech.edu/608/1/Scherer_lr_1950.pdf
  5. ^ Anderson, Jr., John D. (2019). Hypersonic and High-Temperature Gas Dynamics. AIAA Education Series (3rd ed.). American Institute of Aeronautics and Astronautics. pp. 58–67. ISBN 978-1-62410-514-2.
  6. Lees, Lester (1955). "Hypersonic Flow". Journal of Spacecraft and Rockets. 40 (5). Los Angeles: Institute of Aeronautical Sciences: 241–276. doi:10.2514/2.6897. ISSN 0022-4650.

Further reading

  • Abbott, I.H. and Von Doenhoff, A.E. (1959) Theory of Wing Sections, Dover Publications, Inc. New York, Standard Book No. 486-60586-8
  • Anderson, John D (2001) Fundamentals of Aerodynamic 3rd Edition, McGraw-Hill. ISBN 0-07-237335-0
Categories: