Prandtl–Glauert transformation explained

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

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

\phixx+\phiyy+\phizz=

2
M
infty

\phixx(inflowfield)

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

Vinftynx+\phiyny+\phiznz=0(onbodysurface)

Minfty

is the freestream Mach number, and

nx,ny,nz

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

Vinfty

which is assumed here to be along

x

.

\vec{V}=\nabla\phi+Vinfty\hat{x}=(Vinfty+\phix)\hat{x}+\phiy\hat{y}+\phiz\hat{z}

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

|\nabla\phi|\llVinfty

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)

\phix
Vinfty

\right]

2
M
infty

<1

The Prandtl–Glauert (PG) transformation uses the Prandtl–Glauert factor

\beta\equiv\sqrt{1-

2}
M
infty
. It consists of scaling down all y and z dimensions and angle of attack by the factor of

\beta,

the potential by

\beta2,

and the x component of the normal vectors by

\beta

:

\begin{align} \bar{x}&=x\\ \bar{y}&=\betay\\ \bar{z}&=\betaz\\ \bar{\alpha}&=\beta\alpha\\ \bar{\phi}&=\beta2\phi \end{align}

This

\bar{x}\bar{y}\bar{z}

geometry will then have normal vectors whose x components are reduced by

\beta

from the original ones:

\begin{align} \bar{n}\bar{x

} &= \beta n_x \\\bar_ &= n_y \\\bar_ &= n_z\end

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(inflowfield)

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

Vinfty\bar{n}\bar{x

} + \bar_ \bar_ + \bar_ \bar_ = 0 \quad \mbox

This is the incompressible potential-flow problem about the transformed

\bar{x}\bar{y}\bar{z}

geometry. 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_, \bar_ in the transformed space. The physical linearized pressure coefficient is then obtained by the inverse transformation

Cp=-2

\phix
Vinfty

=-

2
\beta2
\bar{\phi
\bar{x
}} = \frac \bar_p

which is known as Göthert's rule

Results

For two-dimensional flow, the net result is that

Cp

and also the lift and moment coefficients

cl,cm

are increased by the factor

1/\beta

:

\begin{align} Cp&=

Cp0
\beta

\\ cl&=

cl0
\beta

\\ cm&=

cm0
\beta

\end{align}

where

Cp0,cl0,cm0

are the incompressible-flow values for the original (unscaled)

xyz

geometry. This 2D-only result is known as the Prandtl Rule.

For three-dimensional flows, these simple

1/\beta

scalings do NOT apply. Instead, it is necessary to work with the scaled

\bar{x}\bar{y}\bar{z}

geometry as given above, and use the Göthert's Rule to compute the

Cp

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

CL=

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

cl0=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.

History

The interest in compressibility research emerged after the WWI, when the aircraft propeller tips started to reach M=0.8. Ludwig Prandtl had taught the transformation in his lectures about 1922, however the first rigorous proof was published in 1928 by Hermann Glauert. The introduction of this relation allowed the design of aircraft which were able to operate in higher subsonic speed areas. Originally all these results were developed for 2D flow. Göthert eventually realized in 1946 that the geometric distortion induced by 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 in 1925. 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

See main article: Prandtl–Glauert singularity.

Near the sonic speed

Minfty\simeq1

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

Minfty\simeq1.

The more correct nonlinear equation does not exhibit the singularity.

See also

References

Sources