The solution of the laminar-boundary-layer equation for the flat plate for velocity and temperature fields for variable ...


Material Information

The solution of the laminar-boundary-layer equation for the flat plate for velocity and temperature fields for variable physical properties and for the diffusion field at high concentration = Über die lösung der laminaren grenzschichtgleichung an der ebenen platte für geschwindigkeits- und temperaturfeld bei veränderlichen stoffwerten und für das diffusionsfeld bei höheren konzentrationen
Portion of title:
Über die lösung der laminaren grenzschichtgleichung an der ebenen platte für geschwindigkeits- und temperaturfeld bei veränderlichen stoffwerten und für das diffusionsfeld bei höheren konzentrationen
Physical Description:
19 p. : ; 28 cm.
Schuh, H ( Herbert )
United States -- National Advisory Committee for Aeronautics
National Advisory Committee for Aeronautics
Place of Publication:
Washington, D.C
Publication Date:


Subjects / Keywords:
Boundary layer   ( lcsh )
Laminar flow   ( lcsh )
federal government publication   ( marcgt )
bibliography   ( marcgt )
technical report   ( marcgt )
non-fiction   ( marcgt )


Includes bibliographic references (p. 15).
Statement of Responsibility:
by H. Schuh.
General Note:
"Technical memorandum 1275."
General Note:
"Report date May 1950."

Record Information

Source Institution:
University of Florida
Rights Management:
All applicable rights reserved by the source institution and holding location.
Resource Identifier:
aleph - 003689761
oclc - 76823498
sobekcm - AA00006222_00001
System ID:

Full Text
~kcA-mi~ci K







By H. Schuh


In connection with Pohlhausen's solution for the temperature field
on the flat plate, a series of formulas were indicated by means of
which the velocity and temperature field for variable physical charac-
teristics can be computed by an integral equation and an iteration
method based on it. With it, the following cases were solved: On the
assumption that the viscosity simply varies with the temperature while
the other fluid properties remain constant, the velocity and tempera-
ture field on the heated and cooled plate, respectively, was computed
at the Prandtl numbers 12.5 and 100 (viscous fluids). A closer study
of these two cases resulted In general relations: The calculations
for a gas of Pr number 0.7 (air) were conducted on the assumption that
all fluid properties vary with the temperature, and the velocities are
low enough for the heat of friction to be discounted. The result was
a thickening of the boundary layers, but no appreciable modification
in shearing stress or heat-transfer coefficient. The effects of
density and viscosity or density and heat conductivity have opposite
effect for velocity and temperature field and almost cancel one
another. Formulas allowing for the heat produced by the friction were
indicated, but no calculations were carried through in view of the
already existing report by Crocco. The methods of solution developed
here were finally applied also to the case of diffusion of admixtures,
where at higher concentration finite transverse velocities occur at the

*"Uber die Losung der laminaren Grenzschichtgleichung an der
ebenen Platte fur Geschwindigkeits- und Temperaturfeld bel
veranderlichen Stoffwerten und fur das Diffusionsfeld bel hoheren
Konzentrationen.4 Zentrale fir wissenschaftliches Berichtswesen der
Luftfahrtforschung des Generalluftzeugmei sters (ZWB) Berlin-Adlershof,
Forschungsbericht Nr. 1980, August 18, 1944.

NACA TM 1275


The laminar-boundary-layer equation for the flat plate lined up
with the flow and with constant fluid properties was solved by Prandtl
(reference 1) and Blasius (reference 2) for the flow field and by
Pohlhausen (reference 3) for the temperature field. Since temperature
and velocity field coincide when kinematic viscosity (V) and tempera-
ture conductivity (a) are identically equal Pr= y = ), Pohlhausen's
formula for the temperature field contains a solution for the velocity
field also in the form of an integral equation. Piercy and Preston
(reference 4), proceeding from a rough approximation, indicated that,
with the aid of this integral equation and an iteration method, the
well-known Blasius solution can be obtained in a few steps. This
method of solution has the advantage of being simple and requiring
relatively little time. It is as is shown in the following -
particularly suitable for boundary-layer calculations involving
variable fluid properties, because a first, and usually fairly close,
approximation, is already available in the solution for constant fluid

Crocco (reference 5) and von KArmAn and Tsien (reference 6) (see
also reference 9, 10) computed velocity and temperature field for
variable fluid properties. In both reports, the differential equations
are put in a different form from the elsewhere conventional boundary-
layer calculation by changing to new variables. Crocco obtains two
simultaneous differential equations of the second order which he solves
for a gas with the Prandtl number Pr = 0.725 (air). Von Karmnn and
Tsien treat the case of Pr = 1 and have to solve only one differential
equation, since then the temperature is related in a simple manner to
the velocity.

In the following, it is shown that a number of boundary-layer
problems for the flat plate can be solved in a comparatively simple
manner, involving merely quadrature, by means of the cited integral
equation and an iteration method.


The boundary-layer equations for velocity and temperature field at
the flat plate at variable density read (reference 3)
Cu bu = / ou
p aru- +p (1)

NACA TM 1275

c (ou) + Cp = (la)
ax 3Y

-pT v (r)r (2)

with u and v the speed in flow direction and at right angle to it,
T, the temperature, x the distance from the plate leading edge,
y the distance from the wall, o the density, P the viscosity, Cp the
specific heat, and X the heat conductivity. In the equation for the
temperature field, the heat produced by friction is, at first, not
taken into account; as long as the speeds are not excessive and the
temperature differences not too small, this is justified.

For constant density, equations (1) and (2) can be reduced to an
ordinary differential equation (reference 3) on the assumption that u
and T are a function merely of the one dimensionlesss)

coordinate g = Since the density depends only on the tempera-
ture, the idea suggests itself that the same simplification is possible
also for variable density. We put

U T T) 1 (o (3)

where U is the velocity at the edge of the boundary layer,
To and T1 the wall temperature and the temperature at the edge of the
boundary layer, respectively. The quantity Vk in the dimensionless 5
denotes the kinematic viscosity for the fixed temperature Tk, for
which in suitable manner the wall temperature (k = 0), or the tempera-
ture at the edge of the boundary layer (k = 1), is chosen. The
boundary conditions for flow and temperature field read

y=0 =0 m=0 0 =
y -> = 1 0 = 1


p( x (5)

NACA TM 1275

where the subscript k denotes the density at temperature Tk, gives
by (la)

pv out ou dt (6)

hence, by (1), after introduction of (5) and (6), the more suitable
_Ld (/W f = 2j dt (7)

From (7), regarded as differential equation for the quantity cp and
f temporarily as a known function of e, the following expression for
a) is derived

J(f) 1 0 d (8)
J(m) = O, T

This disposes of the integration constant from consideration of the
boundary condition (4). Likewise, there Is afforded for the dimen-
sionless temperature 0 the expression

K(E) Prk I
S= K) -e di (9)
K(0) X

where Prk is the Prandtl number with the density at temperature Tk.

For constant density (rp = = X = 1), velocity and temperature
field are independent of each other and (9) gives the Pohlhausen
expression (reference 3) for the temperature field, which represents
the solution for the velocity field at Pr = l(V = a). When the
velocity field is known, the solution for the temperature by (9) is
obtainable by simple quadrature. But the calculation of the velocity
field runs into difficulties, at first, because in (8) the still
unknown velocity appears on the right-hand side In the expression
for f. The methods of solution by Piercy and Preston proceed from a
random approximation for a with which f and J(t) in (8) are

NACA TM 1275

computed. The improved value m obtained forms the starting point for
the next step, etc. Figure 1 represents the several steps of this
approximation method. The intentionally rough approximation c = 1
over the entire boundary layer was chosen as original solution; the
corresponding first approximation c is given by the error integral.
After the third approximation, the shearing stress shows a mere
difference of h.5 percent from the exact value. Instead of continuing
the process mechanically, the final solution to be expected was
estimated from the variation of the previously computed approximations
and utilized as basis for the subsequent step; the solution wu
contained but a -percent error in shearing stress.

With this method of solution, the improvement effected by each
step can be estimated according to order of magnitude. The
equations (8) and (9) are identical for constant density and Pr = 1.
Assuming that the approximate solution for w was such that for each
individual value w the corresponding t coordinate differed by a
constant factor e from the E coordinate of the exact solution, the
effect of factor X is then obviously just as great as that of
quantity Pr for the temperature field. Pohlhausen found, on the
basis of his numerical calculations, that the heat-transfer coefficient
is proportional to \P-r, thus the shearing stress at the wall is
afflicted at each new step by an error of only about one-third of the
error of the preceding step.

For variable density, the discussed solution steps of "mathe-
matical nature" can be combined with the steps of physical nature :

Step 1: as starting point the known solutions for constant
density are assumed:

(a) The Blasius solution (reference 2) for the velocity

(b) Pohlhausen's method for the temperature field

Step 2:

(a) Calculation of velocity profile by (8), the temperature
variation being based on the density of the tempera-
ture profile according to step l(b)

(b) Calculation of temperature field by (9) with the velocity
profile according to step 2(a); relation of density to
temperature as in step 2(a)

lit took a subsidiary worker 10 hours to reach the final solution of
the velocity field in figure 1.

NACA TM 1275

The process is repeated till the final solution is sufficiently
exact, usually requiring three to four steps.

A few general remarks about the influence of the temperature
variability of the physical properties.- The flow with constant
physical properties can be regarded as first approximation, and the
problem is then to ascertain the differences which are produced by
variable physical properties. The quality of this approximation
depends, of course, on the temperature assumed for the physical
properties at the isothermal" flow. Choosing the wall temperature or
the temperature at the edge of the boundary layer as reference tempera-
ture for the isothermal flow so results on the basis of physical
point of view as well as on the basis of the equations that an increase
of the viscosity or density inside the boundary layer is accompanied by
an increase in the resistance; similarly, an increase in heat con-
ductivity and density effects a greater heat transfer. But the
magnitude of the effect of variability of the separate physical prop-
erties is contingent upon the ratio of the boundary-layer thickness of
the temperature and velocity field. (The ratio of both is proportional
according to Pohlhausen.)

This is illustrated by the following case, which is, at the same
time, of practical importance. The temperature boundary layer is
assumed very small compared to the flow boundary layer; consequently,
the variation of the physical properties within the thermal boundary
layer can be disregarded for the shearing stress and the latter
computed as if the temperature at the edge of the boundary layer
reaches to the wall. The same holds for the velocity profile, with the
exception of a small area within the thermal boundary layer, where the
velocity profile by the viscosity variation is deformed correspondingly.
But for the temperature profile this area is exactly decisive.

From the equality for the shearing stresses the velocity gradients
at the walls are:

[l o1 for Pr '-- (10)

the subscript 11 denotes the isothermall" flow with the physical prop-
erties at temperature Tl. The variability of density is noneffective for
the field of flow, in this instance. It can be mathematically derived from
the formulas (8) and (9). The ratios for the temperature field are
discussed in the next chapter by means of the two examples.

NACA TM 1275


In accordance with the physical properties of viscous fluids, the
velocity and temperature field were computed on the assumption that
only the .viscosity should change with the temperature by the following

(Tkk Tcb

where b and Tc are constants, chosen so as to reproduce the tem-
perature variation as closely as possible. The subscript k is to
be 0 or 1, depending upon the choice of the physical properties
in the dimensionless t The choice was b = 3 (viscous lubricating
oil) and the two cases of a heated and cooled plate computed with
and 8 and Pr = 12.5 and 100; it thus concerned identically

great temperature differences of the same fluid, since Pr is
for the present formed with the physical properties at wall temperature.
Choosing To as reference temperature gives by (11)

= .M.= 3 (12)
o 1 +i

The result of the calculation by the iteration method of the preceding
section is seen in figures 2 and 3. In both graphs, the dimensionless
wall distances o and l1, formed with o and ~1, are plotted to

the scale 1:V and V:1, respectively, so that the actual wall
distance y is the same for both abscissas. Besides the solution m,
which took three steps to compute, the isothermal velocity profiles
(0)) and (Um)1 at constant density at temperature' T and T1 are
shown plotted against the dimensionless coordinates o and 1".2
In the subsequent compilation To and a denote the shearing stress
at the wall and the heat-transfer coefficient ( = )o );

For the isothermall" temperature profiles (e) and (0) the
Prandtl numbers at temperatures T and T1 must be inserted. For
example, in figure 2: Pr = 12.5 and Pr = 100.
o 1

NACA TM 1275

(To and (a)) are the corresponding values in isothermal flow with
the viscosity at temperature T1 and T .

Table I

oo oC
(T0 a P( / d0 \

Heated wall 0.841 1.20 12.5 1.58 1.84

Cooled wall 1.08 0.98 100 0.255 3.01

Although in both cases the thermal boundary-layer thickness is far
from small compared to the flow boundary layer, the shearing stress
can still be computed satisfactorily by the isothermal formula with the
viscosity of the wall temperature. The assumptions to equation (10)
are thus shown for Pr > 10.

The conditions are more complicated for the heat-transfer
coefficient; from (9), it follows that the heat-transfer coefficient

a is proportional to a(Pr) ~\ wherein, according to Pohlhausen,

a is, with high accuracy, assumed as 0.664 V ;. Bearing in mind that
Pr = V it follows that the heat-transfer coefficient is inversely
proportional to the sixth root of the viscosity. Since all physical
properties except the viscosity have been assumed constant, there
results, when it is referred once to the wall temperature, the other
time to the temperature at the edge of the boundary layer
()o 1/6 (13)

A comparison with the foregoing tabulation indicates that (a)1
supplies a poorer approximation for the heat-transfer coefficient
than (a)o; this is readily explained by the variation of the velocity
profile (figures 2 and 3). It is to be expected that the conditions
are similar at higher Prandtl numbers.

NACA TM 1275

Another reference point for the heat-transfer coefficient is found
in the velocity gradient at the wall; by (10) and allowing for (3),
there follows

o I1- (14)
pL 11 Po oL 1 o

the subscripts io and il denoting isothermal flow at tempera-
ture To and T1. These relations are confirmed in figures 2 and 3.

With these formulas, limits can be indicated for the heat-transfer
coefficients (figures 2 and 3). One is given according to (13)
by (a)l, because the velocity profile (wu) yields at all points
higher velocities at cooled and lower velocities at heated wall. The
other limit is given by a velocity profile of isothermal form, where
the abscissa scale is so modified that its gradient at the wall agrees
with the actual velocity distribution. From the remark about the
convergence of the method of solution in II, it follows then that the
heat-transfer coefficient is proportional to the third root of the
velocity gradient at the wall; for this extreme value, the second

equation of (14) gives: ( (a)o. Summed up, the limits of the

heat-transfer coefficients, by a change in viscosity, are

(a.) a (15)

the upper signs applying to heated, the lower to cooled wall.

Hence, the following approximate rule for viscous fluids (Pr>10):
For computing the resistance, the physical properties are referred to
the temperature at the edge of the boundary layer; for heat transfer,
to the wall temperature.



In the -500 to 1400 temperature range the physical properties
of the air can be represented by the following formulas

NACA TM 1275

= ITO .780

p = K T-1

X = K T0"821

T = temperature in absolute degrees. With

T T1


there results

= = [ +(i 078

and similar expressions for r and X.
1 1
The calculations for a heated plate and = and 2 showed
only moderate differences in the velocity and temperature field from
the form for isothermal flow (Table 2). For the investigation of the
conditions at higher temperature differences, the case T1 = 20
and To = 6200 C was computed. The velocity and temperature fields
already exhibit, according to figure 4, appreciable differences from
the form for constant physical properties; to and l1 are formed
with the physical properties at temperatures To and T1, respectively.
This results in a substantial thickening of the boundary layer for both
fields; nevertheless, wall shearing stress and heat-transfer coefficient
indicate only minor departures from the values for constant physical

Table 2

_/cku \0 \e 'o '''o a a,
Heating ( d) To T a a

1 0.575 0.90 1.02 1.00 1.01 .00o

S=1 .514 .420 1.05 1.00 1.02 1.00

To = 620 C
.286 .235 1.11 .93 1.03 .96
T1 = 200 C

NACA TM 1275

The explanation for it is that in air the growth of the viscosity with
the temperature acts in the sense of a resistance increase, the drop
in density in the sense of a resistance decrease, and both effects
practically cancel one another at Pr = 0.7, where thermal and flow
boundary-layer thickness are about equally great. The conditions for
the temperature field are almost identical, because the heat con-
ductivity and the viscosity are similarly affected by the temperature.

The frictional heat can be allowed for in similar manner;

equation (2) contains then an additive term 4 C
side, and the solution reads

S1+B(-)Ag() B(S )

A(t) = f le-R() d

Pr, AT,, tT
B(t) = 2 Prk e '
(Tk To) T

R() = Prk dt
k 0 X

on the right-hand

~. e'2 ep( d e- R() d

e 2cp


The iteration method can be applied again, although a little more
paper work is involved. For constant physical properties, equation (16)
reduces to Eckert's solution (reference 8). The thermometer problem
(vanishing temperature gradient at the wall) can also be solved by
suitable variation of the integration constant. In view of Crocco's
calculations for a gas with Pr = 0.725, it was decided not to cal-
culate any model problems by the new method.

NACA TM 1275


The concentration field for the problem of diffusion at the flat
plate can be calculated in the same manner as the temperature field.3,4
The differential equation reads

3c + v = k 2c (17)
u vx ,y 6y2

where k is the diffusion factor and c the concentration which is
defined as quantity of gas or vapor per unit volume. The physical
properties are regarded as constant, but it is also taken into account
that for greater concentrations the velocity v at the wall no longer
disappears, as already pointed out by Nusselt (reference 7). When
fluid from a wall is vaporized, say by a gas such as gas flowing along
a wetted wall, substance passes continuously into the flow.
Hence v(O)>0 at the wall. When, on the other hand, vapor condenses
at the wall or when air containing ammonia, for example, passes over
blotting paper impregnated with hydrochloric acid, it results
in v(0)<0., The boundary conditions for v are according to the
equations (100) and (101) of reference (7);

k- Jc 1 v(0) (18)
co\ io 2- 1
where c is the concentration of the gas or vapor, for which the wall
is permeable, co the concentration at the wall, p the corresponding
partial pressure, and p the total pressure.

Introduction of the flow velocity U and the dimensionless 5
results in

(vj 1 k cl c dC c co (19)
U 2 co x c -co

3Eckert reported a solution of this problem at the 1943 meeting
of the VDT committee for heat research in Bayreuth, where an approxima-
tion method similar to Pohlhausen's method for the flow boundary layer
was used.
4 "
Damkohler's estimate for the present problem was published in
Z. fur Elektrochemie, 1942, p. 178.

NACA TM 1275

where co and cl are the concentration at the wall and at the edge
of the boundary layer. Similarly to (9), it gives

vU v 2]
E w d [W 2

and similarly to (9), the solution for the concentration field

(f( )-M)dt

L(6) = 0


f(E) = 2 u-d
J0 U

where 1 is a quantity analogous to the Prandtl number.
velocity w, simply put 1 = 1 in equation (21).

The concentration gradient at the wall is contained
can, in the first instance, be solved for any M values

N =- -l o computed with the aid of the value

o The velocity and concentration
the solution for %.) The velocity and concentration
*t 0

calculated M



To obtain the

in M; but (21)
and the quantity

obtained from

fields for the

and N values are represented in figures 5 and 6, the

concentration gradient at the wall, in figure 7. M >0 denotes evapora-
tion at the plate; M<0, condensation and absorption at the plate; the

value 0.6 chosen for the quantity 1 is applicable in good approxima-
tion for the diffusion of water and ammonia in air. Strictly speaking,
for the specified higher concentrations, the density and viscosity of
the two fluids are dependent on the concentration; and the diffusion
factor, on the temperature. Cases of that kind can be calculated
with the aid of the described method. If the diffusion is

5According to Ten Bosch: Die Warmeubertragung, Berlin 1936,
pp. 189 and 257.

M= k cl co LCo
Vc o o d
S- -
0 I



NACA TM 1275

accompanied by a heat transfer, the solution for the concentration field
can equally be applied to the temperature field with good approximation.
In the same way, the heat transfer can be derived from the solution for
the concentration field, when air is exhausted or blown at the plate
with transverse velocities at the wall corresponding to equation (20).


After completion of the calculations the writer received knowledge
of a report by Schlichting and Bussmann (reference 11) about the
velocity profile at the flat plate for exhaustion where the transverse
velocity at the wall was expressed by

(x) = -

Between the present value M and C the following relation exists
(see also (19) and (20)).

C =-M

The present velocity distributions agree to about 1 percent with those
calculated by Schlichting (by a different method), with exception of
M = 1, where the writer plainly chose too few approximation steps and
the differences are therefore a little greater. For the present
calculation three steps were usually sufficient.

Translated by J. Vanier
National Advisory Committee
for Aeronautics

NACA TM 1275


1. Prandtl, L.: Ueber Flussigkeitsbewegung bel sehr kleiner Relbung.
Verhandl. d. III. Intern. Math. Kongr. Heidelberg 1904, p. 484.
(Available as NACA TM 452.)
2. Blasius, H.: Grenzschichten in Flussigkeiten mit kleiner Reibung.
Z. Math. Phys. Bd. 56, 1908, p. 1,. (Available as NACA TM 1256.)

3. Pohlhausen, E.: Der Warmeaustausch zwischen festep Korpern und
Flussigkeiten mit kleiner Reibung und kleiner Warmeleitung.
ZAMM Heft 2, 1921, p. 115.

4. Piercy, N. A. V., and Preston, J. H:: A Simple Solution of the
Flat Plate Problem of Skin Friction and Heat Transfer. Phil.
Mag. May 1936, p. 995.

5. Crocco, L.: Sullo strato limited laminare nel gas lungo una parete
piana. Rend. Circ. Mathem. Palermo Bd. LXITI, 1940/41.

6. Von Karman, Th., and Tsien, H. S.: Boundary Layer in Compressible
Fluids. Jour. Aero. Sciences, Vol. 5, No. 6, 1938, p. 227.

7. Nueselt, W.: Warmeubergang, Diffusion und Verdunstung. ZAMM
Band 10, Heft 2, 1930, p. 105.

8. Eckert, E., and Drewvtz, 0.: Der Warmeubergang an eine mit grosser
Geschwindigkeit langs angestromte Platte. Forsch. Ing.-Wes.
Bd. 11, 1940, p. 116. (Available as NACA TM 1045.)

9. Hantzsche, W., and Wendt, H.: Zum Kompressibilitatseinfluss bei
der laminaren Grenzschicht der ebenen Platte. Jahrb. 1940 der
Dtsch. Luftfahrtforschung, p. I 517.

10. Hantzsche, W., and Wendt, H.: Die laminare Grenzschicht der ebenen
Platte mit und ohne Warmeubergang hunter Berucksichtigung der
Kompressibilitat. Jahirb. 1942 der Dtsch. Luftfahrtforschung,
p. I 40.

11. Schlichting, H., and Bussmann, K.: Exakte Losungen fi die lamlnare
Grenzschicht mit Absaugung und Ausblasen. Schriften der Dtsch.
Akademie der Luftfahrtforschung, Heft 2, 1943, P. 25.




/.0 2.0 3.0

Figure 1.- The several approximations for computing the velocity
distribution at the flat plate by the method of Preston and Piercy
(constant physical quantities).

/.0 2.0 3.0 = YF vo
2 0 vx

Figure 2.- Velocity and temperature distribution at a heated plate
for variable viscosity. Viscosity exponent b = 3. (u)o, (e)o'
and (u)1, ()1 isother,.al velocity and temperature distributions,
v and v1 kinematic viscosity at wall temperature To and
temperature T1 at the edge of the boundary layer.



NACA TM 1275

02 0.4 0.6 = o.
0 A onA

Figure 3.- Velocity and temperature distribution at the
Viscosity exponent b = 3. (u)o, (e)o, and (w)1, (e)l
velocity and temperature distributions, vo and v,
viscosity at wall temperature To and temperature
the edge of the boundary layer.

cooled plate.
T1 at

I.O 2.0o _?o v,
I 2 xv,

Figure 4.- Velocity and temperature distribution at a heated plate
for Pr = 0.7 (air); all physical quantities constant with

NACA TM 1275

b. 2.0

0.6 ______
2 0 51,____ __A -

Figure 5.- Velocity field at diffusion with higher concentrations, where
finite transverse velocities occur at the wall (see text for equations (18)
to (20)).

to (20)).

Figure 6.- Concentration distribution to figure 5.

NACA TM 1275

2.0 /.

AM dC k v p

\ k c,-co (C

dC --. ...

d ^?'0
1.0 0.- 5


-I/ / I3 4 5 6 7

-/.o --0.5

Figure 7.- Concentration gradient at the wall and quantity M plotted
against N (see text for equations (19) and (20)).

NACA-Langley 5-19-50 -B25


1 .8
de 8

o ,c

k0 O


r- I
1 4 g
0rd t <4l

0 i 4 >

0 O 0
4 P,- 0

r-1H 1

B ^ ^

1X1 41



a a 4


o 4

+3 a)

q-1 -I H

0 w rm

0 T< 14


E- (4
E-4 0















4) W)P 0 a -P"
0 i-l -H P r F4,--4

H O (D $4

-0 -P $3 H 1) P bo .t-i
S- o 0 d d a>
S-p 0 rl 0 c k M

P m -P r 0) O
0 a 0 0o Hp. ,
0 Op 4 p >o
O(4OO o i-l M
04 H0 H U

O'4).ar ao a pa
0 u i flH 0 0 0 00
0 r-I w -- 4) 9) >3 0 H +

0 0p 0. 4- m -
0 O m .T o e
c r_ rd ( 4 0 l
4- -p 00a SO 4) l
SR0 .0 0, +3o
Cd u) r-4 ? i 44-4 4-4 m

S 4.- r-i p a )4 0
I Hl r-4 0 0- V> u 4m O
0- i r-4 U r-4 ad M 0 -r,
$dO aO i 4 > 4 o
-tl4e -c 0 3 H v (W M
H N 0 4 rS P+S P h 0 -
0 0 P- F4 > o a>
r OMa p o| mI --
P 00 0 o Cd + 0

C0 r- d : r o p Cd P 40
P cc o Vc o 0 04 P fl-
P- 0 O r 4) > I -r (O
4-'J 0-P 4 0 4-) +- 0 -iW

) $3 (U ;3 4 +3 rO C d
0 4 Q- F W
) HN uUm r O+ r-a
O-, 14 P a -+.- 0 0

.Ht O M ) >i Sp r-

OW O Od Cd
d H 0d 0 P
to 0 co r-

$i.1-4 0) a) 4 0) P 4-)
H a CC -4 4-+

O )P M 4) -P 4 1- H
Oa 0 4A M 0 c P u 4

aH rdo$ d- (D
P r. o P P +1 u m S,

m 0i- 0 kr R u
A a .-p4 m a r- 14

0d -P 04 ) q0
4) ao-H o 40 S OI
.0 C(0 0 N 0D

kcT) M 0 > 5-i 04
4-4 H3 c i a w- (a

s S-1) mo o >

q O- M>4 H -. H 4
Bl OW meM aI
*r -r4 L) U 4i 4 > S H

ra -iP O .4 3 aj V Do aP P
4' 0 CS tb 13 m 'A l p o oM

Up0 -H CB u P Cd

k U O r. X q : o a
4J C O (D -P d d 3 D r- l dS
m 0 d ( 0 l c)

0hACR m-H O
V4 h I3 O- 4 4 cu M 0 (D
4O-P4 Pi^ n H L r-i r mo

)M a 0e m 171> r-l



S- o m (U E
a1) 0 C 0d (a
O- H' 0 -P 01
r 001H0 4
a4 Q, M 0 ) -4 11 0 C -i
0a C. S V) p ;:5 C 0 "
*H 0 ; E-4 () 0 -4 H 0 13
-t P. o i-H ) 4- ; 0 0i
d 0 do ( L
p I, U 3 p, a (11 4_1 I S.
0 4-) Z ,, -P Oi tIO

F1 + d 0 .,-i o :0 o ?-4
aCc4 410 l cV 01
O 0 s co a M d
o cdO 4 c (1 cd

r_ c0 a 0 -4 O 0 $-q
0 C 4 r i -ii 0- 0 0
- 01 4c cd 0 >
om 0o a) rU H o
-,O a. P x H 003

o' w. a ac a 0 :S
_0014 0 0 0T0
-4 o o o --
0 cc 1a1 > 0 a:
0iW4-.' V0. HO iin

o r O O P
0 '- r u o u+ o o 4)

P- 0 w> 0> m w i0
o m: rc 0 Or d )

OO 0+ a U 0 P, 0 a)
4 d 0 4 W > )0 (l W '
-4 001 bo-4r r4

6 0 P. ( 0 P.
;I 0 s 0r0 C 4 > 0
0 >0-.- 0 0 >0 4
H0 Q 4 O 0 +'

o r-" ?04 a1 0 d00 0a V 4
14o 10 O' 00- i a O ()
P..-1ir--4 a) > ) rn ri-l a)i U

E-i1 a U W B- Ca
a W0 01 0. i 0 a)C
(,-j1-i 4i +J .r< i1 O
rl OM a0. U-r-Md 4 r> l
Vk-i o Q-v r0V m0a piP.

0 0O I)
a)U 0 +o 0

0 4 H c ).o 4

S0 .-4 P010W 2
0- D 1 (U + l *l r :3 0
Oir4 rd -a(3 4 C

o P ) 'd 0 W
P) r -4 0 ) r_ Id al c

Ao o 0 -Ha 0 ri 0o o
0 4a )O -4 -'O P M0

(LI ci o C ci cd v}
S0 ti 0 w P O

l P. 4 r $l 0 ^ $1
4 co (U h Cd _Hi
0O r d diP

o Hna) co o : cad
So 0 a *VHH cIW
;14 a 0C 4 0S 0 a)

0 or 0 c0 0 -0
*0 .1kr4>1 > ,D ID < i"
) +3 ^ ^i +) ID 0r-
S l R +H- 0 0 0
a -I o 00 o. o
"C C 0 14 rH .- m D0 0

,. 001d > 0> +-) r
MIc O O q > w (U

0 ro 4> M a 8 0 >z F1
; 42 01 0 .- C C
W 0 > 0 1 +4 > 0 r-
s a I +) ) r4 r
-,0 00 0 0 ) P a
0 r- c : Q 0 w + a
i (0 0 crrl O C O 010 O'-H 00r4 C)
P. r4H Wa t ) > a) i-H H )D ;4 U
4-)H 0 00 .0 -P001
SO *n r- hk +O 3 Hl Hd A r
S(u 'R d a a) BR P
a P: k M H4-) (Ua)C
WL 4-) -H W.0 0 ICd -.0 -4 d)
Ci14 -4 -i m 00 1 n
H- 0 -1 r 0 -H C rl k H-
d C-i 0 -H cd 0a m 00 pI p

3 1262 08105 053 5