Thin Shear Layer Approximation Thin Shear Layer Approximation

1  Viscous flow equations

Steady incompressible laminar viscous flow is described by the Navier-Stokes Equations together with appropriate boundary conditions.
Ñu
=
0
u ·Ñu
=
-Ñp/r + n Ñ2 u
in freestream:
u   =  u¥
on body:
u   =   ®
0
 
(1)
Near a solid wall, the tangential velocity goes to its zero no-slip value through a thin thear layer of thickness d, which varies depending on the distance L to the leading edge as shown in Figure . At large Reynolds numbers ue L/n >> 1, the shear layer becomes very slender such that d/L << 1.

thin_shear_layer.gif

Figure 1: Thin shear layer geometry and velocities.

2  Scales and ordering

For a slender shear layer, scales for the velocities and velocity gradients can be obtained from the edge velocity values ue, ve.
u
=
O(ue)
v
=
O(ve)

x
=
O æ
ç
è
1
L
ö
÷
ø

y
=
O æ
ç
è
1
d
ö
÷
ø
1
r
p
x
=
O æ
ç
è
ue2
L
ö
÷
ø
The last pressure gradient scale follows from the Bernoulli relation which is valid outside the shear layer.

Ordering the terms in the continuity equation

u
x
+
v
y
=
0
ue
L
ve
d
and the self-evident requirement that the only two terms in the equation must be of the same order, produces the scaling rule for the vertical velocity v.
ve
=
O æ
ç
è
ue d
L
ö
÷
ø

The terms in the x-momentum equation have the following ordering.

u u
x
+
v u
y
=
- 1
r
p
x
+
n 2 u
x2
+
n 2 u
y2
ue2
L
ue2
L
ue2
L
n ue
L2
n ue
d2
Since d << L, the first viscous term can be clearly neglected over the second. The remaining terms then produce an order estimate for the shear layer fineness ratio d/L.
ue2
L
=
O æ
ç
è
ue d
L
ö
÷
ø
d
L
=
O é
ê
ë
æ
ç
è
n
ue L
ö
÷
ø
1/2

 
ù
ú
û
   =   O æ
ç
è
1
ReL1/2
ö
÷
ø

Ordering the terms in the y-momentum equation

u v
x
+
v v
y
=
- 1
r
p
y
+
n 2 v
x2
+
n 2 v
y2
d
L
ue2
L
d
L
ue2
L
1
r
p
y
n d
L
ue
L2
n d
L
ue
d2
produces the conclusion that the normal pressure gradient is or higher order than the streamwise pressure gradient.
1
r
p
y
=
O æ
ç
è
d
L
ue2
L
ö
÷
ø
   =   O æ
ç
è
d
L
1
r
p
x
ö
÷
ø
The difference in pressure Dpy across the shear layer is clearly negligible to the differences in pressure Dpx along the shear layer.
Dpy    ~   d p
y
   =   O æ
ç
è
ue2 d2
L2
ö
÷
ø
=
O æ
ç
è
d2
L2
Dpx ö
÷
ø

For a curved wall, the leading-order approximation to the normal pressure gradient becomes

- 1
r
p
y
=
O ( ue2 k)
where k is the wall curvature. In this case,
Dpy    ~   d p
y
   =   O ( ue2 kd)
=
O ( kdDpx )
and the pressure differences across the shear layer are still small, since kd << 1 is typical of most shear layer situations.

This analysis indicates that the pressure can be assumed constant across the shear layer. The pressure gradient term can then be expressed only in terms of the tangential edge velocity.

p(x,y)   @  pe(x)
- 1
r
p
x
   @ - 1
r
d pe
d x
   =   ue   d ue
d x

3  Thin Shear Layer Equations

The above analysis shows that for large Reynolds numbers, the steady incompressible Navier-Stokes equations can be reduced to the simpler Thin Shear Layer Equations.
u
x
 +  v
y
=
0
u u
x
 + v u
y
=
ue d ue
d x
 + n  2 u
y2
(2)
The derivation was based on considering of the flat-plate geometry in Figure 1, although these equations are also valid for boundary layers on curved walls and free shear layers and wakes along curved streamlines, provided the

In such cases the coordinate x is the arc length along the layer, y is the transverse distance, and u,v are the respective streamwise and transverse velocity components as shown in Figure .

curved_tsl.gif

Figure 2: Curved shear layers and coordinates.

4  Thin Shear Layer Categories

Equations (2) are valid for all types of shear layers, which are distinguished by different initial and boundary conditions. The general initial condition at xo is
u(xo,y)   =  uo(y)
v(xo,y)   =  vo(y)
The various possible downstream boundary conditions for x > xo are given below.

4.1  Wall boundary layer

bc_wallbl.gif

edge:
u(x,ye) = ue(x)
wall:
u(x,0) = 0 ,       v(x,0) = 0

4.2  Wake

bc_wake.gif

For a symmetric wake, suitable boundary conditions are:

edge:
u(x,ye) = ue(x)
centerline:
u/y(x,0) = 0,       v(x,0) = 0
while for a general assymetric wake they are:
top edge:
u(x,ye+) = ue(x)
``centerline¢¢:
v(x,yc) = vc(x)
bottom edge:
u(x,ye-) = ue(x)
The ``centerline'' position yc and specified vertical velocity vc are both arbitrary, provided vc << ue. This condition merely serves to position the x,y coordinate system within the shear layer flowfield.

4.3  Mixing layer

bc_mixlayer.gif

top edge:
u(x,ye+) = ue+(x)
bottom edge:
v(x,yc) = vc(x)
bottom edge:
u(x,ye-) = ue-(x)
This differs from the wake in that the top and bottom potential flows have different velocities and hence different total pressure.

4.4  Jet

bc_jet.gif

top edge:
u(x,ye+) = 0
bottom edge:
v(x,yc) = vc(x)
bottom edge:
u(x,ye-) = 0