 Open Access
 Total Downloads : 226
 Authors : Chinyere Mbachu
 Paper ID : IJERTV4IS080445
 Volume & Issue : Volume 04, Issue 08 (August 2015)
 Published (First Online): 25082015
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Asymptotic Solutions and Range of Validity for Heat Transfer in PackedBed and Capillary Tubes
Chinyere Mbachu11 1Department of Chemical Engineering, Tennessee Technological University,
Prescott 214, 1020 Stadium Drive,
P. O. Box 5013, Cookeville, TN 38505, Tennessee, USA
Abstract: In this paper, an asymptotic analytical solution was used to evaluate convectivediffusive heat transfer in packed bed and capillary tubes. Asymptotic solution to the heat transfer problems found in packedbed systems with axial convective, radial diffusive with constant heat flux at the wall, and in capillary tubes were studied and analyzed from the results of the analysis reported by Bird et al. [1] and Arce and Trigatti, [2]. The role of asymptotic solution and its ranges of the validity of both for the packedbed and capillary tubes were assessed and ranges of validity reported.
Keywords: Asymptotic solution; heat transfer; convection; capillary tube; packedbed tube.
I INTRODUCTION
while diffusion dominates the transport in the radial direction of the flow [5]. The fluid is assumed to exchange heat at a constant rate without a volumetric source present within the bulk domain.
B Mathematical Model Equation for Capillary Tube
The basic equations for the system under study and the boundary conditions were properly derived from the general heat transfer equation given by Bird et al. [1]. For the capillary tube, the mathematical model for the steady state convective diffusion equation under the assumptions stated is given by
Capillary and packedbed tubes are used in
T 1
T
various industrial applications and are important to the
pC pVz r z
k r r r r
(1)
practical aspects of engineering processes such as reactors, heaters, evaporators, thermal, and energy storage units [1, 2]. Analysis of capillary tubes with laminar flow regime and packedbed tubes for Newtonian fluids under incompressible conditions was presented with focus on
Where
Vz r
r 2
Vmax 1
(2)
assessing the range of validity of an asymptotic solution for the heat transfer problem. A restriction on the asymptotic solutions is based on the fact that convectivediffusive transfer within the domain is controlled by a constant rate at the capillary walls. The assessment of the asymptotic solution is very important in the analysis of the design of reactors or heat exchangers [3]. Several studies have been done to correlate heat transfer by forced convection from circular tubes [4]. However, it appears that the validity ranges are not reported in the literature.
R
The boundary conditions for equation (2) and the parameters are introduced as follows:
at r = 0, T = finite or, alternatively, 0 (symmetry) (3)
r
II MODEL FORMULATION
A General Assumptions for the Mathematical Model
The geometry of the capillary and packedbed tubes are assumed to be cylindrical with a length that
at r = R, k T
r
at z = 0, T = T1
q0
(constant) (4)
(5)
satisfies the long channel assumption (R/L 1), see
By defining the following dimensionless variables [1]
figures 1 and 2. The energy equation for capillary
T T1 ;
r ; Z z
(6)
cylindrical coordinates is given by Bird et al. [1] and was described in section B. The fluid is assumed to be a Newtonian fluid and the flow is assumed to be pressure
q0 R / k
R c
p vz,max
R2 / k
driven, fully developed flow, and incompressible. Convection is assumed to dominate in the axial direction
The capillary equation of heat transfer then reduces to
2
1
at 0 ,
finite
(8)
(1
) Z
(7)
1
at ,
1
(9)
The boundary conditions in equations (1) – (5) in non dimensional terms are
at 0 , 0
(10)
An asymptotic solution for this model was explored in the sections below.
Figure 1: System description of a capillary tube.
C Mathematical Model Equation for PackedBed Tube
The convectivediffusive transport equation for the case of a Poiseuille flow in a packedbed tube as reported by Arce and Trigatti [2] is given by
Where v is the areaaveraged velocity of the fluid moving through the packedbed, and with the same assumptions stated above, the wall of the packedbed tube at Rt shows that the system (see figure 1) either gains or losses heat at a constant rate q0 and at
T k T
pm c pmv z
r
r r r
(11)
the center of the tube, symmetry is assumed by Arce and Trigatti [2]. The boundary conditions for equation 11 can be stated as
Then, by using eq. (12) and eq. (13), the following non dimensionless forms are as follows
T 1
0 ;

k T
q ; and T x 0 T
0
Z
(14)
r r 0
r r Rc
0
(12)
And if the dimensionless variables are given as
And the boundaries and entrance conditions are given by
r ;
z ; and
T T0
0;
1;
and
Z 0 0
pm
R m c
vR2 / k
R
q
k
0
(13)
0
1
(15)
An asymptotic solution for this differential model was obtained below.
Figure 2: System description of a packedbed tube.
III ASYMPTOTIC SOLUTIONS

Capillary Tube Model Equation
The fact that this solution is valid for case z 0 implies that the condition at z = 0 needs to be replaced by a global balance as stated by Bird et al. [1]. This is the reason to
use , Z instead of simply , Z .
2 R
For distances that are away (z 0) from the entrance of the capillary tube, the temperature profile can be approximated by a linear relation in the zdirection as reported by Bird et al. [1]. This assumption yields the mathematical formulation of the asymptotic solution and in terms of the nomenclature used above as follows
2Rzq0 pC p T T1 vz rdrd
0 0
Now, in dimensionless form, this equation leads to
1
(17)
0
, Z C Z
(16)
Z , Z 1 2 d 0
(18)
As stated by the Bird et al. [1] function in equation
(16) must satisfy
Therefore, the final solution for the asymptotic solution for the cylindrical tube is obtained (see Figure 3) and given by
1 d d 2
, Z 4Z 2 1 4 7 (29)
C 1
d d o
(19)
4 24
B PackedBed Tube Model Equation
And by applying the boundary conditions
0
at 0 , = finite and (20)
Similarly to the case of a capillary tube, the asymptotic solution for the Poiseuille flow with convectivediffusive heat transfer in a packedbed tube can be written as
at 1 , d 1
dr
(21)
, C
(30)
By solving equation (19), the following result is obtained
In order to determine function , this must satisfy the
2
4
following differential equation
C0 C1 ln C2
(22)
d d
4 16
d d c0
(31)
Where C1
an
C2 are constants of integration that by
applying the boundary condition given by equations (20) and (21) can be determined as
As well as the boundary conditions
C1 0
C 7
2 24
(23)
(24)
d
d 0
0 ; and
d
d 1
1
(32)
The equation (22) leads to the following general solution for
The solution is reported by Arce and Trigatti [2] and only a summary of the key steps for the convenience of the reader was done. As before, see Equation (17), the condition at the inlet of the packedbed must be replaced
2
4 7
by an integral condition valid also for the function
C0
4
16 24
(25)
, Z . This equation can be derived using an integral energy balance in the system.
2
R
[Energy due to transport by convection]
Substituting equation (25) into Equation (16) and by applying the asymptotic solution, then the temperature
0
0 rdrdpm c pmV (T T0 )
(33)
profile is given by
2 4
[Energy loss at the wall of the tube] wRxq0(34)
0
0
, Z C Z C C ln C
1
4 16 1 2
Now, by using Equation (33) and Equation (34) and invoking the energy balance in the system gives
(26)
' ' '
Or, alternatively,
d ,
0
(35)
2
4 7
After using the nondimensional variables, the solution for
leads to
, Z C0 Z C0
4
16 24
(27)
c
o 2 c ln c
(36)
4
1 2
Constant C0
C0 4
can now be obtained to yield
(28)
And after using the boundary conditions in Equation (32), function is given by
1 1
4
2 2
(37)
, Z changes sign within the domain of validity. This change of sign is a result that the asymptotic solution is no longer valid; therefore, the condition determines the
The constant C0
give
can be obtained from Equation (35) to
values from which meaningful results.
, Z
yields physically
C0 2
(38)

Case 2 PackedBed Tube
For the packedbed tube, the system is assumed to
Finally, the asymptotic function for , z is given by
1
be cooling. From the physical point of view, under Poiseuille flow conditions and in order to establish the ranges of validity for the packedbed tube, the function
, , it is seen that when Z takes smaller values, the
, 2
1 2 2
4
(39)
function ,
tends to cross the axis
, 0 as it should be since the validity of
For the case of the packedbed tube in cylindrical coordinates is shown in Figure 4.
IV RANGE OF VALIDITY OF ASYMPTOTIC
, is for large. A condition for the validity can be obtained by simply stating that
SOLUTION
A Case 1 Capillary Tube
The ranges of validity of the asymptotic solution in a capillary tube can be obtained using Equation (29) and
, 0
With this, Equation (39) can be written as
1 1 2 2
(42)
(43)
after imposing the condition following equation is obtained
, Z 0
, the
8
This equation is used to plot the graph as seen in
Z*
*2
4
1 *4
16
7
for 0 1
96
(40)
Figure 4. From this figure 4, it is possible to see the following: first for 0 1 and Z 0.125 , the
solution
0 and belongs to a physically meaningful
From the physical point of view, this system describes a heating fluid under Poiseuille flow conditions. This implies that according to the definition of used in here, the proper domain for this solution should be 0. Now, values that are located within the domain of 0 (see Figure 3) are not physically meaningful. Therefore, the boundary between the two domains for can be calculated from the condition 0 or alternatively
domain. However, for values of Z 0.125, the asymptotic solution is still valid (partially) for values of
0.5 1 . In other words, if the full domain of
0 1 is desired to be considered, the Z is restricted
to Z 0.125 in order for to be within the range of validity. On the other hand, if the full range of Z
0 Z is considered, then the partial values of
1 *2
1 *4 7
0.5 1
must be used in order for
to be
Z*
4 16
96
(41)
within the region of validity.
The plot of this boundary is shown in Figure 3 and it is observed that the upper region for the boundary between the two domains for is the physical meaningful domain. From Figure 3, this is clearly when Z 0.073
and 0 1 . In other words, for all values
V NUMERICAL ILLUSTRATIONS
Based on the prediction suggested by the validity of the capillary tube, Figure 5 is presented. All values of Z0.073 (i.e., Z=1, Z= 0.075, Z=0.5) comply with the
solution , being within the physically
0 1 , the smallest value of Z , when the solution
meaningful domain. For example, for Z= 0.25, Z= 0.1, Z=
Z 0.073
Z 0.05 and Z= 0.01, the values of , fall outside the
, is valid is . Moreover, if values of are
such that Z 0Z Z 0.073, the solution is only valid for * 1 ; where * = 0.29. The reason
behind this condition as indicated above is the fact that
region of physically meaningful domain. Similarly, Figure 6 illustrates the results for the case of the packedbed tube for example, Z=0.25, Z=0.5, Z=0.75, and Z=1 all have
values of , being inside the physically
meaningful domain for this case. However, the values for
Z=0.1, Z=0.05, and Z=0.01 do not belong to this region. As a final remark, the results of these figures effectively
illustrate the useful predictions of the ranges of validity.
Figure 3: Sketch of the range of validity of asymptotic solution in a cylindrical capillary tube.
(No physically meaningful
(Physically meaningful
Figure 4: Sketch of the range of validity of asymptotic solution in a packedbed tube.
Figure 5: Heating process in a capillary tube.
Figure 6: Cooling process in a packedbed tube.
VI SUMMARY AND CONCLUDING REMARKS
Differential models for two cases of heat transfer (i.e. a capillary system under heating conditions and a packedbed system under cooling conditions) were developed and their boundary conditioned formulated. Asymptotic solutions for both models were proposed by following the suggestions of Bird et al. [1] that are consistent with a case of constant rate of cooling (or heating) at the external walls [6]. Moreover, to the best of
the researchers knowledge, the ranges of validity of both asymptotic solutions were determined and illustrations of the behavior of the systems (within the ranges) were presented. The analysis suggests that for these cases, the asymptotic solutions present limitations for both independent variables as opposed to only one as suggested before. Finally, these solutions and the range of validity are useful to improve the understanding of the behavior of the systems described.
ACKNOWLEDGMENT
Financial support of the Schlumberger Foundation Faculty for the Future Fellowship is greatly acknowledged. Thanks to Offices of Research and Graduate Studies, Tennessee Technological University, Cookeville, Tenessee, USA.
c pm
c0
c1 , c2 DH
NOMENCLATURE
Heat capacity of the cooling/heating fluid Proportionality constant
Integration constants
Diameter
REFERENCES

R. B., Bird, W. E. Stewart, and E. N. Lightfoot, Transport Phenomena, 2nd ed., Wiley, New York, 2002, pp. 311316.

P. E. Arce and I. M. B. Trigatti, Convectivediffusive heat transfer in packedbed with constant heat flux at the wall: An asymptotic solution for large, International Communications in Heat
and Mass Transfer. 21(3), 1994: 435446.

M. Mansour, A. Sheboul, and A. N. Hussein, An analytical solution for diffusion and reaction in a laminar flow tubular reactor, International Communications in Heat and Mass Transfer. 16, 1989: 603608.

S. M. Ross, A mathematical model of mass transport in a long permeable tube with radial convection, Journal of Fluid Mechanics. 63, 1974: 157175.

G. R. Ahmed and M. M. Yovanovich, Analytical method for forced convection from flat plates, circular cylinders, and spheres, Journal of Thermophysics and Heat Transfer. 9 (3), 1995: 516523.

P. E. Arce, A. E. Cassano, and H.A. Irazoqui, The tubular reactor with laminar flow regime: An integral equation approach homogeneous reaction with arbitrary kinetics. Computers and Chemical Engineering. 12, 1988: 11031113.
h Convective heat transfer coefficient
H Width of the packedbed

Thermal diffusivity
kg Global mass constant

Length
3
, m Density m Convective heat transfer coefficient
q0 Heat flux at the wall of the system
r Radial position of the packedbed tube R Radius of the packedbed tube
T Temperature profile
T0 Nondimensional temperature profile v, V Area Averaged Velocity of the fluid
z,max Maximum Velocity in the zdirection
vz (r) Average velocity in the z direction
Greek symbols
Dynamic viscosity (Pa s)
Psi
Pi
Chi
Theta
Gradient
Radial component
Nondimensional axial position