Prandtl–Glauert transformation

From Infogalactic: the planetary knowledge core
Jump to: navigation, search

Lua error in package.lua at line 80: module 'strict' not found.

The Prandtl–Glauert transformation is a mathematical technique which allows solving certain compressible flow problems by incompressible-flow calculation methods. It also allows applying incompressible-flow data to compressible-flow cases.

Mathematical formulation

File:Prandtl-glauert-transform.png
Plot of the inverse Prandtl-Glauert factor 1/\beta as a function of freestream Mach number. Notice the infinite limit at Mach 1.

Inviscid compressible flow over slender bodies is governed by linearized compressible small-disturbance potential equation:[1]


\phi_{xx} \,+\, \phi_{yy} \,+\, \phi_{zz} \;=\; M_\infty^2 \phi_{xx} \;\;\;\; \mbox{(in flowfield)}

together with the small-disturbance flow-tangency boundary condition.


V_\infty \, n_x \,+\, \phi_y \, n_y \,+\, \phi_z \, n_z \;=\; 0 \;\;\;\; \mbox{(on body surface)}

M_\infty is the freestream Mach number, and n_x, n_y, n_z are the surface-normal vector components. The unknown variable is the perturbation potential \phi(x,y,z), and the total velocity is given by its gradient plus the freestream velocity V_\infty which is assumed here to be along x.


\vec{V} \;=\; \nabla \phi  + V_\infty \hat{x} \;=\; (V_\infty + \phi_x) \,\hat{x} \,+\, \phi_y \,\hat{y} \,+\, \phi_z \,\hat{z}

The above formulation is valid only if the small-disturbance approximation applies ,[2]


| \nabla \phi | \ll V_\infty

and in addition that there is no transonic flow, approximately stated by the requirement that the local Mach number not exceed unity.


\left[1 + (\gamma+1) \frac{\phi_x}{V_\infty} \right] M_\infty^2 \;<\; 1

The Prandtl-Glauert (PG) transformation uses the Prandtl-Glauert Factor \beta \equiv \sqrt{1 - M_\infty^2}. It consists of scaling down all y and z dimensions and angle of attack by the factor of \beta, and the potential by \beta^2.


\begin{array}{rcl}
\bar{x} &=& x \\
\bar{y} &=& \beta y \\
\bar{z} &=& \beta z \\
\bar{\alpha} &=& \beta \alpha \\
\bar{\phi} &=& \beta^2 \phi
\end{array}

The small-disturbance potential equation then transforms to the Laplace equation,


\bar{\phi}_{\bar{x}\bar{x}} \,+\, \bar{\phi}_{\bar{y}\bar{y}} \,+\, \bar{\phi}_{\bar{z}\bar{z}} \;=\; 0  \;\;\;\; \mbox{(in flowfield)}

and the flow-tangency boundary condition retains the same form.


V_\infty \, \bar{n}_{\bar{x}} \,+\, \bar{\phi}_{\bar{y}} \, \bar{n}_{\bar{y}} \,+\, \bar{\phi}_{\bar{z}} \, \bar{n}_{\bar{z}} \;=\; 0  \;\;\;\; \mbox{(on body surface)}

This is the incompressible potential-flow problem about the transformed geometry with surface normal vector components \bar{n}_{\bar{x}},\bar{n}_{\bar{y}},\bar{n}_{\bar{z}}. It can be solved by incompressible methods, such as thin airfoil theory, vortex lattice methods, panel methods, etc. The result is the transformed perturbation potential \bar{\phi} or its gradient components \bar{\phi}_{\bar{x}},\bar{\phi}_{\bar{y}},\bar{\phi}_{\bar{z}} in the transformed space. The physical linearized pressure coefficient is then obtained by the inverse transformation


C_p \;=\; - 2 \frac{\phi_x}{V_\infty} \;=\;
-\frac{2}{\beta^2} \frac{\bar{\phi}_{\bar{x}}}{V_\infty}

which is known as Göthert's Rule [3]

Results

For two-dimensional flow, the net result is that the C_p and also the lift and moment coefficients c_l, c_m are increased by the factor 1/\beta over the incompressible-flow values:


\begin{array}{rcl}
C_p &= \displaystyle \frac{C_{p0}} {\beta} \\
c_l &= \displaystyle \frac{c_{l0}} {\beta} \\
c_m &= \displaystyle \frac{c_{m0}} {\beta}
\end{array}

where C_{p0}, c_{l0}, c_{m0} are the incompressible-flow values. This 2D-only result is known as the Prandtl Rule.[4]

For three-dimensional flows, these simple 1/\beta scalings do NOT apply. Instead, it is necessary to scale the xyz geometry as given above, and use the Göthert's Rule to compute the C_p and subsequently the forces and moments. No simple results are possible, except in special cases. For example, using Lifting-Line Theory for a flat elliptical wing, the lift coefficient is

 C_L = \frac{2 \pi \alpha}{\beta + 2/AR}

where AR is the wing's aspect ratio. Note that in the 2D case where AR → ∞ this reduces to the 2D case, since in incompressible 2D flow for a flat airfoil we have c_{l0} = 2 \pi \alpha, as given by Thin airfoil theory.

Limitations

The PG transformation works well for all freestream Mach numbers up to 0.7 or so, or once transonic flow starts to appear.[5]

History

Ludwig Prandtl had been teaching this transformation in his lectures for a while, however the first publication was in 1928 by Hermann Glauert.[6] The introduction of this relation allowed the design of aircraft which were able to operate in higher subsonic speed areas.[7] Originally all these results were developed for 2D flow. B.H. Göthert[8] then pointed out that the geometry distortion of the PG transformation renders the simple 2D Prandtl Rule invalid for 3D, and properly stated the full 3D problem as described above.

The PG transformation was extended by Jakob Ackeret to supersonic-freestream flows. Like for the subsonic case, the supersonic case is valid only if there are no transonic effect, which requires that the body be slender and the freestream Mach is sufficiently far above unity.

Singularity

<templatestyles src="Module:Hatnote/styles.css"></templatestyles>

Near the sonic speed M_\infty \simeq 1 the PG transformation features a singularity. The singularity is also called the Prandtl–Glauert singularity, and the flow resistance is calculated to approach infinity. In reality, aerodynamic and thermodynamic perturbations get amplified strongly near the sonic speed, but a singularity does not occur. An explanation for this is that the linearized small-disturbance potential equation above is not valid, since it assumes that there are only small variations in Mach number within the flow and absence of compression shocks and thus is missing certain nonlinear terms. However, these become relevant as soon as any part of the flow field accelerates above the speed of sound, and become essential near M_\infty \simeq 1. The more correct nonlinear equation does not exhibit the singularity.

See also

Notes

  1. Kuethe, A.M. and Chow, C.Y., Foundations of Aerodynamics, Wiley, 1976
  2. Shapiro, A.H., Compressible Fluid Flow I, Wiley, 1953
  3. Göthert, B.H. Plane and Three-Dimensional Flow at High Subsonic Speeds (Extension of the Prandtl Rule). NACA TM 1105, 1946.
  4. Erich Truckenbrodt: Fluidmechanik Band 2, 4. Auflage, Springer Verlag, 1996, p. 178-179
  5. Shapiro, A.H., Compressible Fluid Flow I, Wiley, 1953
  6. H. Glauert, The Effect of Compressibility on the Lift of an Aerofoil. Proc. Roy. Soc. London. VOL. CXVIII, 1928, p. 113–119.
  7. Meier, H.-U.: Die Entwicklung des Pfeilflügels, eine technische Herausforderung, Ludwig Prandtl memorial lecture, GAMM 2005, March 28th - April 1st 2005, Universität Luxemburg, Kapitel 1
  8. Göthert, B.H. Plane and Three-Dimensional Flow at High Subsonic Speeds (Extension of the Prandtl Rule). NACA TM 1105, 1946.