Q J Mechanics Appl Math 1958 WATSON 302 25
Q J Mechanics Appl Math 1958 WATSON 302 25
Q J Mechanics Appl Math 1958 WATSON 302 25
SUMMARY
Exact solutions of the Navier-Stokes equations are derived by a Laplace-transform
technique for two-dimensional, incompressible flow past an infinite plane porous wall.
It is assumed that the flow is independent of the distance parallel to the wall and
that the velocity component normal to the wall is constant. A general formula is
derived for the velocity distribution as a function of the given free-stream velocity
and, from this, general formulae for the skin friction and displacement thickness are
obtained.
Among the particular cases considered are (i)flowin which the free-stream velocity
is changed impulsively from one value to another, (ii) uniform acceleration of the
main stream from a given value of the velocity, and (iii)flowin which the free-stream
velocity consists of a mean value with a superimposed decaying oscillation. In the
first case a secondary boundary layer of a Rayleigh impulsive flow type is created
immediately after the impulsive increase in free-stream velocity and this secondary
layer increases in thickness with time until the final steady state is reached. In the
second case also, a secondary boundary layer is formed at the instant when the
free-stream velocity begins to accelerate and ultimately the skin friction consists
of a quasi-steady part and a part proportional to the product of the free-stream
acceleration and a ' virtual mass' equal to the mass offluidin the displacement area.
In the final case the fluctuating skin friction is found to have a phase advance over
the fluctuating free-stream velocity. For sufficiently large frequencies this phase
advance is JIT and the amplitude of the skin friction fluctuation is proportional to
the square root of the frequency and to the amplitude of the free-stream velocity
fluctuation.
1. Introduction
LIGHTHTLL (1) has investigated the effect of fluctuations of the external
stream on the boundary layer of a two-dimensional body. In addition,
Stuart (2) has obtained an exact solution of the Navier-Stokes equations
for the case of fluctuating flow past an infinite plane wall with constant
suction, and has compared some of his results with LighthilTs. This paper
is an extension of Stuart's work to the case of general unsteady flows past
the wall, but with the suction still held constant.
The generalization of the velocity field is effected by means of a Laplacetransform technique; many forms of the time-dependent external stream
[Quart. Journ. Mech. and Applied Math., Vol. XI, Pt. 3, 1958]
303
304
J. WATSON
2. List of symbols
x
y
t
tc
p
u
v
U(t)
Uo
vw
p
/x
v
Q(y)
f(t), g(y, t)
T
Tc
T)
Q.
8 10
81
r,n
(2 f
305
1erf a;.
fcosa;da;
(2TT)* J
S(x)
1
a;*
J
T 1 fsinaw
Fresnel integral -7-
8u
du
- +
l(8p\ . -8H'
8y=-p[-8x)+V--^
8t
p\8yj
"
8v
8y=
(3>1)
du
dU .
82u
and
5092.43
(3.4)
3. General theory
Let x denote the distance along a two-dimensional, infinite, plane porous
wall, y the distance normal to it, u and v the corresponding components
of velocity, t the time, p the pressure, p the density, and v the kinematic
viscosity. For two-dimensional incompressible flow which is independent
of x, the Navier-Stokes equations and the equation of continuity are
306
J. WATSON
T = tvl/tv = vZ/48*0)
(3.6)
v = y\vJlv,
where 810 is the steady displacement thickness {v/\vw\). Substituting (3.5) in
(3.2) and utilizing the fact that 0 is given by (3.4) we obtain
)
8g(r,,T) _
d*g(r,,T)
4
/ U ) + 4
(37)
e+ioo
The inverse is
x(T) =
p
2-m JJ
c too
A formal solution of (3.7) by the Laplace transform technique will now be
derived. To do this we note that the transform of the first term on the
left-hand side of (3.7) is
00
PJ
dT =
where (?/) is given by (3.4) and where f(t) is a given arbitrary function
and g(y, t) is a function to be determined. The present paper deals with
an investigation of this generalized case. J
Let us first make the transformation
provided that e~pTg(-q, T) -> 0 as T -> 00, where g(r),p) is the transform
of g. Similarly, the transform off'{T) is pf(p), provided that e-Tf{T) -+ 0
as I 1 ^- 0- The Laplace transform of (3.7) is therefore
d2g
dg
d? + dj-Z9--4f
p ,
(3 8)
f After this work was completed, the author found that the solution when U(t) = 0
for t < 0 has been given by Hasimoto (4).
J There is also a solution of (3.1) in which the steady part of the free stream has a
uniform shear. The unsteady part g of the velocity profile is independent of the shear
(Curie (5)) and so is given by the solution of thi3 paper for zero shear.
307
-x^p^^p)
is
-r
f x1(X)x2(TA) dX.
_oo
(y)_e-iij 1 .
f
2ni J
epTe-i^a+p)
c ico
+
c'+ioo
~T
[ epTpe-Wp*?!. (c* =
2T7&
4 .-i-n.-r
21
o
dp
c'+ioo
[I ^vliTln^^^-^v^P
^^ JJ
c'-ioo
= S(T)
where S(T) is the delta function and H(T) is the unit function. The
composition-product rule then gives
00
, T) = J f(T-X)J8(A)= f(T)-\-r
A-/(T-A)exp - A + -^)\dX.
(3.11)
(3-9)
g - > / as T) -> oo J
If & = K^+V(l+?')} anc ^ that branch of the square root is chosen which
makes *J(l+p) = + 1 when p = 0, then the functionale~ki>) satisfies
(3.8) and the boundary conditions (3.9); the real part of k is positive for
all values of p with this choice of the square root. I t follows that
g=f(l-e-kv)
(3.10)
is the unique solution of (3.8) subject to (3.9).
It is now necessary to determine g. We first note that (3.10) is the
product of two functions of p one of which, f(p), has a known inverse,/(T).
Thus the function g is best obtained by using the composition-product rule,
namely that the inverse of
308
J. WATSON
^J
v = vw
for finite T. In such a case, it can be shown that g(rj, T) as given by (3.11)
satisfies (3.7) and, at points where/(T) is continuous, that it satisfies the
boundary conditions. Consequently (3.5), where g(y, t) is given by (3.11) and
(3.6), satisfies (3.1) and the corresponding boundary conditions at points
of continuity of/( T); on physical grounds the solution is unique. Functions,
f(T), which do not satisfy the above conditions may, however, lead to a
solution; in such cases it must be checked that (3.11) satisfies (3.7) and the
boundary conditions. From equations (3.4) and (3.11) the velocity distribution, (3.5), is given by
U(T)/U0 =
.
(3-12)
where T, -q are related to t, y by (3.6). It can be shown that if the freestream velocity ultimately becomes steady, then the velocity distribution
reaches a form similar to the undisturbed velocity distribution; that is if
lim U(T) = U(oo), where U(co) is a finite constant, then
6
which will now be evaluated from (3.12) under the assumption that the
integral
o
309
CO
CO
1 + ^ J A-*/(T-A)e-^A - 2 J/(T-A)erfcA* dX ,
oo
oo
-*
(3.14)
It now remains to derive a formula for the skin friction, TW, which is
^
,3.15,
0
II
J,=o
V-*O r]U0
since u -> 0 as -q -> 0. From (3.12) and by making use of equation (B,2),
this becomes
,
rJPU0\vw\ = 1+/(T) + _L J A-i{/(T)-/(r-A)}e-A dA, (3.16)
o
provided that the integral in (3.16) is absolutely convergent. This is so if
at the point T
\f(T)-f(T-X)\ < C\X\k for |A| < 8,
(3.17)
where C is a positive constant and k > . We note that a function which
satisfies (3.17) is necessarily continuous at the point T; while a function
which has a finite derivative at T satisfies (3.17) and the skin friction is
given by (3.16).
4. Applications of general theory
In the following examples all the conditions on f(T) = JJ(T)jU0\ are
satisfied in order that the velocity distribution, the displacement thickness
and the skin friction be given by (3.12), (3.14), and (3.16) respectively.
4.1. Periodic velocity field
Application is made to the case f(T) = eeinT, where fi = ^vw/vl,, in
order to demonstrate that certain results given by Stuart (2) can be derived
from the results given in section 3. The second equation in (3.12) becomes
pl-(\
+ --)\
dX
[[ \\ 16A/J
16A
477* JJ
o
Uo
CO
e e 4TT*
?
o
where a = (1-f-^)* with re(a) > 0 and b = 77. Using the evaluation of
the integral given in Appendix B, equation (B, 4), we obtain
u/U0 =
le-v+eeinT-eeinT-!Te-W1+in>
477*77
= l-e-v+eeinT(l-e-hv),
(4.1)
310
J. WATSON
i
^
o
oo
= l+eeinTh,
(4.2)
1
where h = ^{1-f (1+iQ) }, by use of equation (B, 6) of Appendix B. This
agrees with Stuart's result.
4.2. Impulsive velocity fields
(a) Single impulse
. The case of a single impulse at T = 0, namely
U(T) = Uo for T < 0;
U(T) = *7O(1+A) for T > 0
where A is a constant, is now considered. This corresponds to .
f(T) = AH(T),
where A is constant and H{T) is the unit function. The velocity is given
by (3.12) and, because of the factor H{T\), the integrand is zero for
A > T. Consequently it is only necessary to integrate from 0 to T when
T is positive; for negative values of T the integrand is identically zero.
Thus, provided T =fc 0 the velocity has the form
(*")
o
L
\
/
\ -*
/J
forT ^ 0. It is easily seen from this that, as T -* oo, u/U0-> (1e~7')(l+A)
as we would expect from Appendix A with f7(oo) = /0(l-j-A).
c(tiClT
311
(4.6)
o
so that (4.4) becomes
(4.7)
for T # 0.
Finally, the skin friction, T^, given by (3.16) is, for
A*
fe- A dA = f e-A dA
C -AdA
JJ IT J "1? J
rJPU0\vw\
for T # 0.
A graph of the skin friction, !,, for A = 0-2 is given in Fig. 1. Eventually
TW reaches the steady value r8 which, from (4.9), is given by
rJpUoK\ = 1+A.
(4.10)
It can be seen from Fig. 1 that at T = 0 the skin friction initially goes to
plus infinity and then decreases to its limiting value as T increases. The
reason for these properties of the skin friction can be seen as follows. Just
i - ^ _ = 77*erfT*.
(4.5)
o
Together with this result, a double integration by parts of the second
integral leads to
312
J. WATSON
-0-1
I
0
0-1
0-2
0-3
0-4
0-5
0-6
0-7
0-8
0-0
1-0
I-1
1-2
1-3
1-4
t-S
1-2
313
/0Aerf(?7/4Ti),
TI = pU0\vw\.
(4.11)
Tc= |7lkZjL_lW,
(4.12)
(4.13)
tc = 8?0/v,
(4.14)
U{T) = luo{l+b)
[Uo
(T<0)
(0<T<To)
(TQ < T)
forces'. Henco outside this layer, though still in the boundary layer, the disturbances to
the viscous forces are negligible, and the interaction is quasi-inviscid.
and is thus like the Rayleigh boundary layer for impulsive flow from rest.
It is of some interest to define a non-dimensional time, Tc, which is
characteristic of the time taken for the skin friction to change from one
steady value to another. In the case under consideration this can be denned
as follows. Let rx be the undisturbed skin friction, namely,
314
so that
J. WATSON
(0
(T < 0)
f(T)=f*(T)-f*(T-T0),
(4.15)
f*(T) = MI(T).
, T) = A^(T)[l-ie^erfc^-^j-ierfc^+yi)j (4.16)
From (3.12) it is seen that the part of the velocity distribution which depends
on T is linear i n / so that by (4.15) it is composed of the part arising from
/ *{T) together with the part arising from / *{TTQ). But, from (4.3), the
part arising from/ *(T) is
i-J
t-erfr } -2TerfcTH.
(4.20)
(4.21)
A graph of the skin friction, rw, for A = 0-2 and To = 0-04 is given in
Fig. 2. For T < To, <p{TT0) = 0 from (4.19), and (4.18) and (4.9) are
identical so that, for T < 0-04, Fig. 2 represents the same skin friction as
Fig. 1. By the same reasoning as in section 4.2 (a), a decrease in free-stream
velocity from f70(l+A) to Uo at T = TQ causes,,thefskin friction to go to
minus infinity and then to increase to its limiting value. From (4.18) and
(4.19) this steady limiting value, rs, is
rJpV0\vu\ = 1.
(4.22)
Again, as in section 4.2 (a), a new secondary boundary layer is created at
T = To the thickness of which is of order v^TT^ for small positive
(TTo) and as time progresses it increases in thickness and modifies the
flow to the ultimate steady velocity U0(le"^). By the same reasoning as
315
in section 4.2 (a), theflowoutside this new boundary layer may be regarded
as inviscid as far as the impulse at T = To is concerned and for T > To the
velocity there is decreased by Uo A so that at T = To the velocity for TJ > 0
is impulsively decreased by this amount.
ym 12
Uo
0-02
0-C4
006
007
008
FIG. 2. Behaviour of skin friction in flow with two impulsive changes of velocity
/lU.loul
316
J. WATSON
I A *exp< (A-|!- aA = o
and equation (B, 1) gives
i2.
+ r)}] (4.25)
J
0
r
,
-2KH(T) j (y-A)erfcA^ dA
1
2KTH(T)
{l+KTH(T)}
TWT}[1+
T
_2KH(T)
r xh_x
dX_2KTH{T)
erfcAjiA
0
T
+ 2KH(T) f\ertc\id\
0
. (4.26)
J
Of the integrals on the right-hand side of (4.26), the first is given by (4.5), the
second can be derived from (4.5) and is
T
(4.27)
the third is given by (4.6), and the fourth can be obtained by a double
integration by parts which, after using (4.27), gives
^ i ( 4 . 2 8 )
\
T
2TT*
4TT'
o
so that (4.24) may be written in the form
317
for all T.
Finally, from (3.16) the skin friction, rw, is, for finite T,
.T
= l+KTH(T)
JJ
oo
^ I
X~h~x dX . (4.30)
'
But
CO
(4.32)
PU0\vJ
TJPU0\VJ
KH(T)\ f
^
x
+ -V-^
X~h~
dX+T
477
(4.29)
318
J. WATSON
but for large T with free-stream velocity U0{\-\-KT) the skin friction
from (4.32) is
rw = PU0\vw\(l+KT)+ipU0\vJK
f.
(4.33)
01
0-2
0-3
0-4
OS
-j-
and Tg is the skin friction if the main stream were constant at U(t). It should
be noted from (4.29) that as T -* oo, SJi^l/v tends to the steady form,
namely unity; thus Sj -> 810. Consequently the displacement thickness
approaches a quasi-steady state in which, although thefluidis accelerating,
the steady displacement thickness applies.
If the outer stream is at rest and the wall is moving then the force required
to move the wall is given by the drag. The force required for a length I and
unit breadth is thus
where A10 is the 'displacement area' (cf. Lighthill (1)). The term Fq repre-
319
sents the drag which would arise if theflowwere steady at the instantaneous
velocity U(t). We can thus call it the quasi-steady drag. The second term,
however, is the product of a 'virtual mass', pA10, with the acceleration,
dU/dt. This 'virtual mass' arises because, as the wall is accelerating, fluid
is accelerated as well through the action of viscosity (cf. Lighthill (1)).
(Uo
7(7) = U0(l+KT)
[U0(l+KT0)
or
f{T) = KT
(To < T)
H(T)K(TT0)H(TT0)
f*{T) = KTH{T).
Precisely the same argument can be applied as in section 4.2 (6) to determine the velocity distribution, skin friction, and displacement thickness;
these are, respectively,
ulU0 = l-e-i+<x(v, T)-ac(r], T-To),
(4.35)
where
^ - ^ - f f l & ^ + r * ) } ] , - (4.36)
T
)
,
0
T
0(T) = \K H{T)[T+\eTiT*+TT-We- +T'erf
0
where
and
SiKI/v =
T*]
(4.37)
(4.38)
[l+y(T)-y(T-T0)]l[l+KTH(T)-K(T-T0)H(T-T0)],
where
.
<4-39>
I
I
ai
I
2I
I8
l
I
y(T) = ZJI(2 )[7r-!Z e- +j7r-*3 *e- -2 erfo2 *+3 erf ^ - i e r f J 1 * ]
(4-40)
for all finite T.
A graph of rw for K = 4 and To = is given in Fig. 4. For T . < To,
P(TT0) = 0 from (4.38), and (4.37) and (4.31) are identical so that, for
T < \, Fig. 4 represents the same skin friction as Fig. 3. As in section 4.3 (a)
it can be shown that the skin friction decreases parabolically for small
positive values of (TTo). For T > TQ there is no pressure gradient so
320
J . WATSON
that the flow will ultimately become steady; it can readily be shown that
the ultimate velocity distribution and skin friction are given by
= U0(l+KT0)(l-e-i) = U(oo)(l-e-n)
and rJpU0\vw\ = \-\-KT0 respectively. For K = 4, To = these are
= %U0{l-e-v)
and rJPU0\vu\ = 2.
of
(4.41)
f
=
i _ e-v
-L
*
1,-
^
I
-1
where a = (1K2-\-iQ.)* with the sqnare root chosen so that re(a) > 0,
321
X (e-yi-
)*}]. (4.43)
Consider firstly the particular case K = 1. Then the free-stream velocity
is
U(T)= U0[l+AH(T)e-TsinQT],
(4.44)
and the complex skin friction is deduced from (4.43) as
. rJPU0\vw\ =
l-&AH(T)e-^(nT)^+e+ePn(in)*ecf{(iClT)i}].
But Jahnke and Emde (8) give the result, with a change of notation,
%
322
J. WATSON
rJPU0\vw\ = l+^Ae-^TH(T)
where
flT)},
(4.47)
, CIT) = [
(4-48)
- t at na n -
'.
'.,
'
B{Cl) = [Q+(2Q)i
2*+Q*
B = Q*,
OL=\TT.
(4.50)
Thus (4.47) represents the skin friction exactly when K2 = 1 and for
Q. ^> |1 K2\ when K2 =fc 1. Since a is positive, the skin friction has a
phase lead over "the free-stream velocity fluctuation. Note that, since
G(x), S(x) - as x--*- oo, for large D.T, (4.48) may be replaced by
U(T) =
(4.51)
But Stuart's (2) result for high frequency may be written in the form
V(T)=
U0[l+esinQT]
(4.52)
= l+eQ*sin(QT+|77)J
which formally corresponds to (4.51) if e = Ae~K*T.. Thus the skin friction
(4.51) behaves as though the main stream were a pure oscillation about a
mean, with amplitude factor Ae~KtT. By comparison with quasi-steady
analogies obtained for other cases, we may refer to this case as quasiharmonic. In other words this quasi-harmonic property follows from the
fact that, for large Cl, the logarithmic decrement of the second term of the
skin friction (4.51),
TJP^OKI
TT
lo gg ll 00 e = 0-4342
04342 ,
323
10
l'&
70
values more and more closely as T approaches 1, so that the limiting values
of B and a are attained for quite small values of T. This is due to the large
value of Q considered. The values of B and a which result from using the
more approximate form (4.50) are B = 3-16 and a. = 0-785, and these do
not differ greatly from the values B = 3-93, a = 0-605.
5. Acknowledgements
The author is greatly indebted to his colleague, Dr. J. T. Stuart, at whose
suggestion the work was undertaken, for much valuable advice and criticism. The work described above was carried out in the Aerodynamics
Division of the National Physical Laboratory, and this paper is published
on the recommendation of the Aeronautical Research Council and by
permission of the Director of the Laboratory.
APPENDIX A
Suppose that the free-stream velocity ultimately becomes steady, that is,
lim U(T) = C7(oo) or lim f(T) =/(oo),
where C7(oo),/(oo) are finite constants. Then, from (3.12),
lan
lu\ =
\TT)
-4T7
^T
lim
r-*oo
0-5
324
J. WATSON
Let
x (A^)=A-*exp(-(A
+^ ) j .
(A, 2)
ff{T-\)x
A,
(A, 3)
The suppositions that /(oo) is finite and f(T) is bounded in any finite interval
imply that/(T) is bounded in oo < T < oo so that
Ao
Ao
{
dA < e
(A, 4)
{/(T-A)-/(o)}x dX
A.
IA;
for XQ sufficiently large (A is a positive constant). Having chosen some Ao to satisfy
(A, 4), then, since |/(TA)/(oo)| < e for T sufficiently large,
(A, 5)
loo
for T sufficiently large, where K is a finite positive constant for 0 < t] < oo. Thus
from (A, 1), (A, 2), (A, 3), (A, 4), and (A, 5) for 0 < rj < oo
lim
lim u = (l-e-")C7(oo).
(A, 6)
T> oo
lim
U(T) = (l-e-")U(-oo).
(A, 7)
APPENDIX B
>Some integrals used in the text
For re(a) > 0, b > 0, x > 0,
x
7(*,o,6) = J e x p ( - (aA+j)JA-i dA
o
( ^ ) ]
(B.Jl
J{x,a,b) =
o
(B,2
325
J(oo,a,6) = je-2ab.
(B,4)
(B, 5)
77*
1.
2.
3.
4.
5.
6.
7.
8.
/(oo,o,0) = .
'
(B.6)
REFERENCES
M. J. LIGHTHTLL, Proc. Roy. Soc. A, 224 (1954) 1.
J. T. STUART, ibid. 231 (1955) 116.
E. J. WATSON, ibid. 104.
H. HASIMOTO, J. Phys. Soc. Japan, 12 (1957) 68.
N. CUKLE, Private communication.
B. VAN DEE POL and H. BKEMMER, Operational Calculus (Cambridge, 1955).
M. J. LIGHTHILL, Proc. Roy. Soc. A, 217 (1953) 478.
E. JAHNKE and F. EMDE, Tables of Functions (Dover, 1945).
and