Boundary Layer-Inviscid Coupling: Appendix
Boundary Layer-Inviscid Coupling: Appendix
Boundary Layer-Inviscid Coupling: Appendix
Fig. El. Flow past an aerofoil. Typical problem for boundary layer-inviscid flow coupling.
Appendix E 303
can be used for different regions. The semi-inverse coupling is introduced here mainly to stabilize the solution in the turbulent region close to separation. f~ipul-e I .2shows the flow diagrams for the present boundary layer-inviscid coupling. Further details on the Thwaites compressible method and semi-inverse coupling can be found in the references discussed in Sec. 6.12, Chapter 6 (Le Balleur and coworkers).
304 Appendix E
In Fig. E.2, Cp is the coefficient of pressure; s the coordinate along the surface; S the boundary layer thickness; 0 the momentum thickness; C the skin friction coefficient; H the velocity profile shape parameter; p the density; V , the transpiration velocity; K* is a factor developed from stability analysis; the subscript Z I marks the viscous boundary layer region; S* the displacement thickness; the superscript i indicates inviscid region and the superscript m indicates the current iteration. Following are useful relations for some of the above quantities:
H = -0 ,
S*
S * = / ~ ( l - ~ ) d n ,
K * = 2x0 E
/?=
(E.l)
where n is the normal direction from the wing surface. We have the following equations to be solved in the integral boundary layer lagentrainment model.
Continuity
u,
ds
Momentum
d0 Cf -=--(H+2-M2)-ds 2
Lag-entrainment
0 du,
u, ds
OS=F[L((CT)&-Aq5)+
ds
H
u,
+HI
(1
(--) 0 du,
u~ ds
EQ
19 du, ds
+ 0.1M2)
0.o2ce F=
0.8Cf + c,+ 3
(0.01
+ C,)
In the above equations, H and HI are the velocity profile shape parameters defined as H=j/r(l-%)dn,
HI =s - s* 0
C , is the entrainment coefficient; ut)the mean component of the streamwise velocity at the edge of the boundary layer; M the Mach number; C, the shear stress coefficient; X the scaling factor on the dissipation length; the subscripts EQ and EQ, denote respectively the equilibrium conditions and equilibrium conditions in the absence of secondary influences on the turbulence structure.
Appendix E 305
Once the above equations are solved, the transpiration velocity V , is calculated as shown in Fig. E.2 and is added to the standard Euler boundary conditions on the wall and plays the role of a surface source. The coupling continues until convergence. In practice, in one coupling cycle, several Euler iterations are carried out for each boundary layer solution.