Integral Boundary Layer Solution

4.11.1 Thwaites method Assumptions and formulation

Integral Boundary Layer Solution

The Thwaites method [21], [11], [12] solves the classical low-speed laminar boundary layer problem by as­suming specific functional forms for Cf, H, and the entire righthand side of the von Karman equation (4.28). This is first manipulated by multiplying 2Req x [equation (4.28)], where

Подпись: T Integral Boundary Layer Solution Подпись: d U d{n/9) Подпись: 62ue s due vs ue ds Integral Boundary Layer Solution Подпись: d 2U d{n/9)2 Подпись: (4.83)

where T is a normalized wall shear, and A is a normalized edge velocity gradient or equivalently a normal­ized pressure gradient. For the power-law flow case ue = Csa these are also products of Falkner-Skan solution parameters, and also of the normalized velocity profile derivatives at the wall.

Their numerical values are listed in Table 4.1. Also listed is the entire righthand side parameter combination Fq in the manipulated von Karman equation (4.81).

For power-law flows, all the Falkner-Skan parameters (columns of Table 4.1) are functions of the expo­nent a, or more generally functions of each other. The key assumption of Thwaites’s method is that the A), Ffl(A) relations arc valid for any boundary layer flow, not just a power-law flow. In effect this assumes that the boundary layer profile U{iy, s) at each s location has the same shape as one of the Falkner- Skan profiles. That is, if A is known at a location s in a general flow, then H, T, Fq at that location are also immediately known from the corresponding row in Table 4.1.

Explicit integration

Rather than use the Falkner-Skan solutions in Table 4.1, Thwaites examined a number of other theoretical flows to quantify these functions. For Fe he chose

Fe(A) = 0.45 – 6A (4.84)

which doesn’t quite match the values in Table 4.1, but it is close. His H and T were provided in tabulated form. They are shown in Figure 4.23, and are closely approximated by the following convenient curve fits.

Подпись:Я, л, = 2.61-4ЛА + 14А* + р24^_

Подпись: Figure 4.23: Thwaites’s values for parameters H and T versus A, plotted as symbols. Curve-fit functions H(A) and T(A) given by (4.85) and (4.86) are plotted as solid lines. Falkner-Skan U (n) profile shapes corresponding to the H values are drawn on top.

T, a, = Um + l-m-SA’-jOTP

Strictly speaking, equation (4.81) could be numerically integrated using only the T(A) and H(A) functions inserted in its righthand side. However, the fact that their combination Fe(A) is very nearly linear allows a simpler explicit solution. When the approximate Fe(A) given by (4.84) is inserted for the entire righthand side of equation (4.81), and the entire equation is then multiplied by vu5, its lefthand side becomes a perfect differential which can then be explicitly integrated.

Подпись: ue d(6|2) v ds = 0.45 - - 6A| vu5 d(«e6^2) ds = 0.45 v u5 92(s) 1 ul(s) u^(so) e2(so) + 0.45 v ( u5e(s') ds' so (4.87)

Integral Boundary Layer Solution Integral Boundary Layer Solution Подпись: (4.88) (4.89) (4.90) (4.91)

If ue(s) and v are provided, the corresponding 9(s) can be immediately obtained by evaluation of the integral in (4.87), numerically if necessary. This 9(s) is then used in the closure relations to obtain the remaining quantities of interest.

In general, the integral in (4.87) is evaluated from some initial location s0 where 9(s0) must be specified. But if this s0 is the sharp leading edge of a Blasius flow, then 9(so) = 0. Alternatively, if s0 is at a blunt leading edge stagnation point, then ue(s0) = 0, in which case 9(s0) in (4.87) is immaterial. Hence, in both of these typical situations no initial data is required.