Int.3.

NeatMassTmnsfer.Vol.14.~P.1047-1056.
PergamonPress 1971. PrintedinGreat Britain

ANALYTICAL METHODS OF SOLUTION OF

CONJUGATED PROBLEMS IN CONVECTIVE
HEAT TRANSFER
A. V. LUIKOV, V. A. ALEKSASHENKO and A. A. ALEKSASHENKO
Heat and Mass Transfer Institute, Byelorussian Academy of Sciences, Minsk, U.S.S.R.

(Received 5 May 1970)

Atfwraqn+-R ]gaHHOa pa6ioTe lSpi?ACTaBJieHbl ME?TOj&l peUieHSiR 3aAaY KOHBeKTHBHOfO

TWUIOO6MeHa C YWTOM paCIlpOCTpaHeHMR TeJljla B TBepEOM TeJI@, conpk%Kacammemca C
&BWKKy~etiCR CpeAOti. 3TBT MeTORb H3:~B~~ibI.M~TO~3M~lp~~eIlIIRCOtIpfUKeHHbIX saaaq.
R YaCTHOCTM,B pa6oTe paCCM3TpllB3eTCHTeIbJIOOCiMeKhpMllaMllHapHOM TeYeHBki HtMp;KOCTH
B Kpyl'JrOtiH l"mOCKOti Tpy6ax C YYeTOM fiMCCMnaqr?kl MeXaHnYeCKOfi 3HepmH. KpOMe TOPO,
paCCMaTpHBaIOTCf3 K8K CTaL&EOHapHbIe, TaK ct KWT3~HOHapHble 3aRaYH TenzOO6MeHa lip&i
06TeKaK~~ n~3CT~K~C~~lMae~O~ ~K~KOCTb~.~p~3TOM Ten~OnepeHOC B~K~KOCT~ BO BCeX
CJIyYaKX paCCM3TpKBaeTCU BO B3a~MOCBX3~ C RepeHOCOE TeRna B CTeHKe TBepAOrO Ti?JIa.
Ha OCHOBe IIOJIyYeHHOrO aHa~~33pe~eHK~BBO~~TC~HOB~~ Kp&iTeptll&KOTOpbIZtXapaKTe-
pe3yeT BJIlfRHlileTeIIJIO&i3RYWKWX CBOtiCTB CTeHKH Ha IIpOl&?CC TeIIJIOO6MeHa. flnH
KJIJIlOCTp3~KHpaCCM3TpABaeTCfi HeCKOJIbKO KOHKpeTHldX IIpklMepOB.

NOMENCLATURE dimensionless spatial co-

temperature ; ordinates (internal prob-
dimensionless temperature lem) ;
(internal problem) ; dimensionless spatial co-
d~ensionless tem~rature ordinates (external prob-
(external problem); lem) ;
temperature at the solid- enthalpy ;
fluid (gas) interface ; tube radius, R, - R, thick-
dimensionless temperature ness of tube wall ;
at the conjugation boundary plate thickness ;
(internal problem); thermal conductivity ;
dimensionless temperature heat capacity ;
at the solid-gas interface Prandtl number ;
(external problem); PC&t number ;
temperature of external tube Mach number ;
surface ; parameter of generalized
dimensionless temperature Fourier sine tr~sfo~ation ;
of external tube surface ; confluent hypergeometric
spatial coordinates ; function ;
1047
IO48 A. V. LUIKOV, V. A. ALEKSASHENKO and A. A. ALEKSASHENKO

h 00. retardation temperature ; All these works are characterized by a general

ho recovery temperature : way of treatment: first, the problem is solved for
t’ l> mean velocity ; the temperature distribution in the boundary
210 = :; A, = f$F,(q) = 1 + $1 - r)
!ayer of the main stream under the prescribed
conditions at the surface, and the heat transfer
coefficient is calculated. Next comes considera-
for a plane tube ; tion of the heat transfer process in the solid.
At the solid-fluid interface there are formulated
v. = 2v, ; A, = F2; F&) = 1 - /?TIn q the so-called third-kind boundary conditions
involving the heat transfer coefficient calculated
for a circular tube ; beforehand.
Thus, there is made an attempt to describe a
complexity of heat transfer processes between
the solid and the main stream by a single heat-
transfer coefficient (many empirical and semi-
4v; Pf empirical formulae have been obtained for its
Hz------. Y, = )I@, ,);
cx,(=Y determination).
In such a statement the mutual thermal effect
&U.b; Re,=!!?; k=!E; of the solid and the fluid is not allowed for, i.e.
1’ V C
”
heat transfer process appears to be independent
ofthe solid properties (thermophysical properties,
J.3 _ 2Nu(q. hw- 4
h = dimensions, source distribution, etc.).
J(Rq)’ w h,; This statement does not seem to be physically
strict, It is worth noting that, as it has become
4‘1

;= 0.55-“.s $ pdy. known nowadays [5,6], the third-kind boundary

conditions are not valid for many cases, since
Subscripts aid Superscripts they lead to contradictory or even physically
senseless results [ 5,,61.
K, refers to main flow parameters ; Consideration of convective heat problems
$ characteristic of a solid ; as conjugated problems [7] appears to be
IV. refers to the surface in a gas flow ; physically more strict, i.e. the energy equations
m= 1, 21 m = 1, a plane tube; m = 2, a for fluid and solid are attacked together with
circular tube. those of hydrodynamics; the temperatures and
heat fluxes at the solid-fluid interface being con-
sidered equal (the fourth-kind boundary con-
CONVECTIVE heat transfer phenomena are of ditions are formulated).
great practical importance in many fields of This is the case when the mutual thermal
modem technology (missile and aircraft engineer- effect of the solid body and the fluid is already
ing, atomic power engineering). taken into account contrary to the alternative
Heat transfer problems in present-day high- statement.
speed aeroplanes involving considerable heating At present a number of publications devoted
of the aerodynamic construction have become of to the solution of conjugated problems [7-91
particular significance. are available.
A great deal of research efforts, both experi- The common points in them are as follows:
mental and theoretical, has been devoted to only particular problems have been solved by
these problems [l-6]. various methods that required physically not
CONVECTIVE HEAT TRANSFER 1049

warranted simplifications in their statement ;

neither new physical effects that should be
expected have been obtained nor general methods
of the conjugated problems have been developed.
Recently [lo, 111 general methods for both
internal and external conjugated heat transfer with the boundary conditions
problems have been developed. These allowed
analytical solutions to be obtained in a rather &l1, = 0 = &I0 &I, * m = A,(1 - @) + 1 ;
(2)
general statement for the cases of real physical
phenomena as well as new physical effects. ae, 0 %Il,=1 = &fm ; (3)
The present paper reports two general solu- all,,,=
tion method for internal and external conjugated
and of the equation for the solid
problems of convective heat transfer.

1. INTERNAL CONJUGATED PROBLEMS (4)

The general method developed for internal with the boundary conditions
heat transfer conjugated problems is based on
reducing the problem to a singular integral e,,I, = o = e,,, %mle - a = F,(r) (5)
equation for the unknown temperature of the
e,,I, = 1 = uo; $-,1, = 6 = hm (6)
surface in a flow [ll]. The method permits
exact solutions for the case of both steady- and At the solid-fluid interface the conjugation
unsteady-state heat transfer in laminar and conditions are
turbulent flows to be obtained.
e,l, = 1 = efmls = 1 = ur); (7)
The present paper is restricted (because of
the lack of space) to consideration of a steady-
state heat transfer problem with laminar forced
convection in circular and plane tubes for the Note, that the conditions at 5 + 00 in equations
developed Poiseuille velocity distribution with (2) and (5) derive from a solution of the problem
allowance for mechanical energy dissipation. stated at 5 + co.
1. Mathematically the problem reduces to a 2. Solutions (l)-(3) are easy to obtain by
solution of dimensionless equations for a fluid applying the superposition principle in the form

x exp[- L(5 - 5J dtl + f

II=1
%,K,,(tl)exp [- b,,5]; I (9)

where
P,, are the rootsof the ch~acte~st~c equation 3. The bo~d~~-v~~c problem solution for
the solid (4H6) is found by the generalized
BmPm)= 0
1Fl@nt, (‘10) Fourier sine transformation
which can be found, for example, in [137.

x exPC- WI (11)
sin iucda > ;
sh ol(b- 1)

Pe~or~~g subs~tutio~ of (9), (11) and (12) into transform of the unknown function x,(a)
the conjugation conditions (Q written for the
transforms, we obtain the equation for the CJ+&)
~~(~~- b,(t) 7 @dt = f+(t); (13)
i, r---f
where,

- am- 1eth ,,/(Q(8 - 1);

CQNVECTIVE HEAT TRANSFER

4. Equation (13) is a singular integral equation

cantoining the Cauchy formulation (14). To
solve it, we take advantage of the concept of
analy#ical extension to the ~Qmpfex regiw and
reduce equation (13) to the Riemann boundary-
value problem with discontinuity coefficients.
On introduction of the analytical function The nonhomogeneous boundary-value Rie-
marm problem (15) may be solved in a general
form, if the problem index is not negative. It is
easy to demonstrate that for the present case the
index is zero [ll] and solution (13) assumes the
where contour L is; a positive part of the real form
axis. Using Sokhatsky-PIemel’s formulae, we
arrive at

Substitution of (15) aBd (16) into f13) yields

W(t) = G(t)@-(t) + g(t); (17)
1052 A. V. LUIKOV, V. A, ALEKSASHENKO and A. A. ALEKSASHENKO

The desired temperatures /3,,,and 8,, are found a2h, a2h, o

-p+dy2= (21)
from equations (9). (11) and (12) by considering
(18) and the relationship in reference [ 12).
with the boundary conditions

x,(t) = (1~4 ;y
-
5x,(4 exp[ - ~1 sin4 du.
0 ah,1
h/
I
5 + 0 = ho0
5-m =O
at1
(22)

EXTERNAL CONJUGATED PROBLEMS $1, = o = he+ x(t); h,(,, _ 1 = h,. (23)

The general method developed for external At the body-fluid interface the conjugation
conjugated problems of convective heat transfer conditions are given as
is based upon series expansion with respect to
powers of the parameter of the generalized hi, = o = hsjy z 0;
Fourier sine transform of the temperature
profile for flow past a surface.
This method has permitted the authors to ?Y
(24)
ah,
=-y;O.
I
obtain the solutions for both steady- and Later on, the unknown function ~(5) is deter-
unsteady-state heat transfer problems for a plate mined from the second conjugation condition
in liquid or gas flows; to allow for radiation (24) for the flows.
effect [ 1l] as well as injection (or suction) ; in As is known, in case of a linear dependence
addition, by the method suggested the authors of viscosity upon temperature, the hydrodynamic
have succeeded to obtain new criterion K for a problem is autonomous (i.e. functions (p(c),
mutual thermal effect of the solid body and the (p’(l), cp”([) are known, this is the ordinary
liquid. Blasius profile [3]).
The present paper is restricted to considera- The solution of problem (19) and (20) for the
tion of a steady-state problem for a plate (length gas is to be sought in the form
L, thickness b, thermal conductivity A,) in a
longitudinal gas flow of constant velocity u, and K,5) = 1+ y M2,Ni)+ liK,O; (25)
temperature T, at infinity ; the internal surface
temperature is assumed to be constant (a similar where
solution is obtained for the case when the flow
is prescribed over the internal surface). s(c) = 20 7 [cp”(i,)]” d9: [~p”(i,)]~ - “d&i (26)
6 0
Gas velocity is assumed to be temperature-
dependent and to follow a linear law ; the Prandtl and Ii satisfies the boundary-value problem
number is considered to be constant.
1. Mathematically the problem reduces to a (27)
solution of dimensionless equations for the gas
q, = 0 = x(C); h(i-co= 0. (28)
a@ - 1)
!?? + $$ - 2acp’$!! = _ ____ The generalized Fourier sine transformation
d12 i X 4
x M2, [cp” K)12 (19) U,(E) = ~li~o 7 u(t) sin g< exp (- ~0 d5 (29)
b
with the boundary conditions
yields the transformed equations (27) and (28)
h[< = 0 = h,+ x(0: h[,+ m = 1 cw

1
a%, ah
and the eqli;irion for the solid F+Ocp2+2arp’([)
al [
r;,++ =o (30)
CONVECTIVE HEAT TRANSFER 1053

h,l,=, = x,(a); &o = 0. (31) 1~~= 0.5 (k - 2): k = 1,2.. . ; j = 1,2.. . ;

(k-1/4)
The solution of equations (30) and (31) is to be An-l-4nG-1-4n
found in a series form x A,_, C,_, + x
c
(2n + l)!
?I=1
h,([, a) = 1 Ck Zvk auk ; 02)
II P 4
where Z(c) satisfies the boundary-value problem -.G-4” k#4j
-1 (2n)! ’
Z”(5) + acpz (0 + 2o(p’ (1 + 1’) Z(5) = 0: (33) n=t
zI,=, = 1 ZI,,, = 0. (34) (I;- l/4)

&,-4nCk-l-4n
C,, vk are constants to be determined later on. %A,_, Ch-, + x
Here the transform of the unknown function is c (2n + l)!
n=l
w41
x,(a) = C Cvayk. (35)
- k = 4j
The transformed solution (21x23) with allow- ’
ance for (29) may be written as n=l

h _ ho0 WW + Y)
I
Applying I$ and Ck values obtained and invert-
fs x,(a) ing by the formula
a sh a

- (ho0- U sh(4 (I+ Y)+ (h,- hoobhccy

a sh a
apti-2rccl + 1)
n 5””
sin

[
c/J+ 1)n .
2 1 ’
(38)

the desired temperature profile in a boundary

(36) layer is determined from equation (35) as
Differentiation of equations (25) and (36) and NC,5) = 1 + IJk - lV81AIf,;K,
substitution of the results obtained into the
+ @b- kJ ; [cp”(t,)l”dMj [cp”K,)l”dW1
second conjugation condition (24) for flows 0

written for transforms yield the equation from + x z’_ l(O) (hb - h&l

which constants vk and C, are determined < [q”(~)/cp”(0)]“~-“‘5 + 0*25Z,.5(5)~-2’5

+0,25x Z;., (0) Z,(<)t- 3
xc Ck AkaYk +“‘5 sh a
k + O.~~X~Z;(O)Z;.~(O)Z~.~(/J~-~‘~
=(h,-hoo)cha+hoo-h,+a~C,a’kcha; +(7.542/16,/n)
k
x3z;.5(o)z;(o)z;.5(o)z3(~)~-4+ . . . .
where
(39)
r(r,. + 1) sin [( w,: + l)ni2]
4 = zv, (0)r( v,. + 15) sin [( vk + 1*5)n/2] With c = 0 the surface temperature in a gas flow
is obtained from equation (39) as
”

x = 0.5 5 ,/Re,.
Lf h, = hl~,,=hl,=,=h,+xZ’_,(O)(h,-h,)
< t-o’5 + 0.25<-2’5 + 0.25~z;.,(O)<-~
Performing expansion in power series of the
+ 0~25x2Z;(0)Z’I.5(O)~-3’5
functions entering into equation (37) and equat-
ing coefficients with like powers a give vI,and C, + (7*542/16Jn)
values determined by the formulae x32;.5(o)2;(o)2;.5(o)6-4 + ... . (40)
1054 A. V. LUIKOV, V. A. ALEKSASHENKO and A. A. ALEKSASHENKO

Note, that functions Z, in equations (39) Numerical calculation is performed with allow-
and (40) with I$ found are determined as the ance for (46) the constant C being determined at
boundary-value problem solution the end of the calculation at 5 = 0. [ll] reports
results obtained by the finite-difference method.
z;, + ocpz;, + 2a(p’(l + v,) Z, = 0
(41) Note that functions Z,, at Ye> 0 are non-
Z&o = 1 Z”I,([+m = 0. monotonous.
The exact solutions are obtained only for two It should be emphasized that solutions (39)
cases, when V, = - 1, V, = - 05. and (40) include the new criterion

z-,(0 = r [v’K,)l”d5, (f [VKJ”dT,)-’ :

x = 0.5$(Re,) = 0*5+” ,blT/,); (47)
(42) I
which comprises characteristics of both the
rp”K) a flow and the body. At 2 + 0 (that corresponds to
Z-,.,(0 = - . (43)
[ cp”(0)
1 an infinite heat conductivity of a solid or
With the other values of K the functions Zv, infinitesimal thickness of a plate b + 0) there is
may be obtained only numerically which re- obtained the known gas temperature profile in a
quires an asymptotic solution of equation (41) flow at constant temperature past a solid body
to be found first. The Blasius functions cp([) and
cp’([)entering into equation (41) are not expressed
analytically but for large [(c 2 35) they may be
represented as
(p(5) E 2([ - 0.86); (p’(l) = 2; r 2 3.5. (44)
Substitution of (44) into (41) produces the
equation for large c Generally speaking, the criterion x may take
on values from 0 to co (for a majority of real
“Yk ([) + 20([ - 0.86)&, (5)
situations x is of an order of unity); note that
+ 4cr(l + Vk)Z(C)= 0 the extreme case x + co corresponds to thermal
Z,‘r(@m = 0 (45) insulation of a body in a flow, and x + 0 is
characteristic of the case when use of Newton’s
the solution of which is iaw is possible with an error that may be found.
In all other cases (0 < x < 00) the heat transfer
((33.5
= coefficient may be only formally introduced.
Hence the Nusselt number with regard for
{exp [ - O.S(l- 086)2]}. (46) equation (40) is determined as

Nu(<) = O-5J(R eg < O-664Y(Pr) - xZ’_ 1(O)[0.25 Z;.,(0)t-2’5

+ 0.25 x Z;(O) Z;.,(O) t-3 + 0.25 x2 Z;.,(O) Z’,(O)z;.,(O) {-3’5
+ (7.5 d2/16 ,/x) x3 Z;.,(O) Z;(O) Z;.,(O) <-” + . . . ] >
x < 1 + xZ’_,(O) [<-o’5 + 0.25 r-2’5 + 025 xZ;.,(O)~-~
+ 0.25 x2 ZXO)Z;.,(O) t-3’5 + (7.5 ,/2/16 J ?I) x3 &.5(O) z;(o) z;.,(o) <-’ + .. *] -1.

(49)
CONVECTIVE HEAT TRANSFER 10.55

It is easy to see that at x + 0 the well-known From the graphs obtained it is seen that h,
Nusselt numbers are obtained from (49) and Nu may differ considerably from those
usually used (at x + 0). Thus, in case (a) at
Nu* = lirn~~ = 0,332 ,,@eJ . .$‘(I+). (50)
individua1 points NU values differ several times
2. Consider* two calculational examples by from the ordinary Nu * ; the same refers to the
formulae derived in [ 1l] temperature of the surface in a flow (since at
(a) Poor conducting material (As = 2.52 x+0, hti = hb). At the same time for a good
W/(mdeg): The plate is b = O.Olm thick. The conducting material h, and Nu values do not
flow characteristics are M, = 3 ; q = 280°K ; differ very much (by not more than 25 Der cent) I

‘Tb,= 223°K 11, = 9.54 x 10-6m2;s* 3 = 2.04 from the corresponding values of h$ and Nu*
x lo-’ W/mdeg; U,(l) = 900 m/s’; 1m~,(2) = (for the case x -+ 0).
400 mjs. In this case x(l) = 3.95; zc2) = 2.64. By
equations (40) and (49) there have been calculated
the temperature of the surface and the Nusselt
numbers. The results are presented in Figs. 1
and 2.
I I I I 1 I I I 1
RFATRENCES
1. H. SCHLICHTING,Grenzschicht-Theorie. Verlag. G.
Braun, Karlsruhe (1964).
2. V. S. AVDUEV~KY ef al, Heut Transfer F~~m~tals in
l-s:i~ ~~~_I Aircrufi und Rocket Engineering. Oborongiz, Moscow,
( 1960).
3. L. G. LOITSIANSKY, Luminur Boundury hyer. Fizmatgiz,
Moscow (1962).
4. E. R. VAN-DRIEST,Turbulent boundary layer in com-
0 IO 20 30 40 50 60 70 80 90
pressible fluids, J. Aeronaut. Sci. 18, 145-161 (1951).
E
5. D. R. CHAPMANand M. W. RUBE%N,Temperature and
FIG. 1. velocity profiles in the compressible laminar boundary
layer with arbitrary distribution of surface temperature,
J. Aeronaut Sci. l&547-565 (1949).
I
I ’ ’ ’ ’ I I I I
6. H. SCHLICHTING, Der WLrmeiibergang an einer linearen
langsangestrijmten ebene: Platte mit vergnderlicher
Wandttemperatur, Forsch. Geb. fngWe.s 17, l-7(1951).
7. A. V. LU~KOV,Heat Conduction Theory. Vyssh. Shkola,
Moscow f 1967).
8. V. R. SHELYAG,Calculation of temperature fields m a
laminar gas flow in flat slot channels, J. Engng Phys. 12
(Z), 227-234 (1967).
9. A. A. POMERANTSEV, Wall heating in supersonic flow,
J Engng Phys. No. 8,39-47 (1960).
2.0
10. V. A. ALEKSASHENKO, The steady conjugated problem
of heat transfer in fluid flow in a semi-infinite tube
I.0
including the effect of &cous dissipation, J.Engng
III I I I III I Phys. 14 (1), 100-108 (1968).
0 IO 20 30 40 50 60 70 80 90 11. V. A. ALEKSASHENKO, Anaiytical solution of some con-
t jugated problems of convective heat exchange, Thesis,
Minsk (1969).
FIG. 2.
12. J. R. MACDONALD and M. K. BRACHMAN. Linear-system
internal transform relations, Rev. Mod. Phys. 393422
(b) Well conducting material (Af = 15.6 f 1956).
W/(mdeg: In this case x(l) = 0.639. Calcula- 13. S. S. KUTATELADZE, F~~~ta~ of Heat Trun.$er
Theory. Mashgiz, Moscow-Leningrad (1962).
tional results similar to the first case are pre- 14. F. D. GAKHOV,Boundary-Value Problems. Fizmatgiz,
sented in Figs. 1 and 2. Moscow(1963).
1056 A. V. LUIKOV, V. A. ALEKSASHENKO and A. A. ALEKSASHENKO

Abstract-The present paper presents solution methods of convective heat transfer problems which take
into account heat propagation in the solid in contact with a moving fluid. The method is referred to as the
solution of conjugated problems.
In particular, the paper treats heat transfer in laminar fluid flow in circular and planar tubes with
allowance for the dissipation of mechanical energy. In addition, there are considered both steady-and
unsteady-state heat transfer problems for flow of a compressible fluid past a plate. In all cases heat transfer
in the fluid is discussed in relation to that in a solid wall.
On the basis of the analysis of the solution a new criterion is introduced which characterizes the effect of
thermophysical properties of the wall on heat transfer. A few examples are considered for illustration
purposes.

METHODES ANALYTIQUES POUR DES PROBLEMES CONJUGUES

DE CONVECTION THERMIQUE

R&me-Cet article concerne les methodes de solution de problemes de convection thermique en tenant
compte de la propagation thermique dans un solide en contact avec un Ccoulement de fluide. La methode
s’inspire de la recherche de la solution des problemes conjuguts. En particulier, l’article traite le transfert
thermique dans un tcoulement de fluide laminaire a I’inttrieur de tubes circulaires et plans avec admission
d’energie mecanique. De plus, on considbre a la fois les problemes de transfert thermique stationnaire et
instationnaire pour une plaque dans un ircoulement compressible. Dans tous les cas, le transfert de chaleur
dans le fluide est discute en relation avec celui dans une paroi solide. A partir de cette analyse, on introduit
un nouveau crittre qui caracterisel’effet des proprietes thermophysiques de la paroi sur le transfert thermique.
En illustration, on considbre quelques exemples.

ANALYTISCHE METHODEN BE1 KONVEKTIVER WARMEUBERTRAGUNG FUR

KONJUGIERTE PROBLEME

Zusammenfassung-Die vorliegende Arbeit berichtet iiber Lbsungsmethoden ftir konvektive Warmciiber-

tragungsprobleme unter Beriicksichtigung der Wlrmeausbreitung in einem festen K&per der in Kontakt
mit einer Flilssigkeitsstrijmung ist.
Im Einzelnen wird die Wlrmetibertragung in einem zylindrischen und cbenen Rohr bei laminarer
Striimung und bei Berticksichtigung der Reibung behandelt. Ausserdem wird sowohl das stationare wie
instationare Wlrmetibertragungsproblem ftir eine cbene Platte bei kompressibler Stromung untersucht.
In allen Fallen wird der Warmestrom in der Fltissigkeit mit dem Wandwlrmestrom des festen Kiirper\
verglichen und diskutiert.
Mit der analytischen Losung wird ein ncues Kriterium. das den Einfluss der thermo-physikalischen
Eigenschaften der Wand auf den Warmestrom charakterisiert, cingeftihrt. Das Problem wird an einigen
Beispielen erllutert.