# Computational Investigation on Fully Developed Periodic Laminar Flow Structure in Baffled Circular Tube with Various BR

Authors: Withada Jedsadaratanachai; Nuthvipa Jayranaiwachira; Pongjet Promvonge
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

This paper presents a 3D numerical analysis of fully developed periodic laminar flow in a circular tube fitted with 45° inclined baffles with inline arrangement. The computations are based on a finite volume method, and the SIMPLE algorithm has been implemented. The characteristics of fluid flow are presented for Reynolds number, Re = 100–1000, based on the hydraulic diameter (D) of the tube. The angled baffles were repeatedly inserted at the middle of the test tube with inline arrangement to generate vortex flows over the tested tube. Effects of different Reynolds numbers and blockage ratios (b/D, BR) with a single pitch ratio of 1 on flow structure in the tested tube were emphasized. The flows in baffled tube show periodic flow at x/D
≈ 2-3, and become a fully developed periodic flow profiles at x/D
≈ 6-7, depending on Re, BR and transverse plane positions. The computational results reveal that the higher of BR and closer position of turbulators, the faster of fully developed periodic flow profiles.


---

## Body

## 1. Introduction

The studies of fluid flow structure and heat transfer behavior are among the most important in the design of heat exchangers in applied engineering works. Many works studied flow structure and heat transfer by using turbulators, ribs fins grooves, or baffles. It is a matter of concern that the enhancement of heat transfer rate depends on effect of turbulators that can change flow structure and behavior in the tested section.Sripattanapipat and Promvonge [1] numerically studied the laminar periodic flow and thermal behaviors in a two-dimensional channel fitted with staggered diamond-shaped baffles. They found that the diamond baffle with half apex angle of 5–10° slightly performs better flow than flat baffle. Promvonge et al. [2] also examined numerically the laminar heat transfer in a square channel with 45° angled baffle placed on one wall. They reported that a single streamwise vortex flow appears and also induces impingement jets on the wall of the interbaffle cavity and the BTE sidewall. Moreover, Promvonge et al. [3, 4] also applied the numerical manner to investigate the laminar flow structure and thermal behaviors in a square channel with 30° or 45° inline baffles on two opposite walls. Two streamwise counterrotating vortex flows were created along the channel lead to increased heat transfer rate over the smooth channel while the maximum thermal enhancement factors of about 2.6 at BR = 0.20, PR = 1, and Re = 1000 and of around 4.0 at BR = 0.15, PR = 2, and Re = 2000 for using the 45° and 30° baffles were reported, respectively.Li and Braun [5] compared flow configuration and transport mechanisms in rectangular and cylindrical enclosures with and without baffles. They found that the presence of the baffle interrupted direct wall layer interactions and circumferential partitioning of the flow. Therefore, the transport of mass and energy was reduced in comparison to the case with no baffle. A numerical investigation of laminar forced convection in a three-dimensional channel with baffles for periodically fully developed flow and with a uniform heat flux in the top and bottom walls was conducted by Lopez et al. [6]. Guo and Anand [7] studied the three-dimensional heat transfer in a channel with a single baffle in the entrance region. Numerical studies for both solid and porous baffles in a two-dimensional channel for the turbulent flow [8] and for the laminar flow regimes [9, 10] were conducted and similar thermal performance results for both the solid and porous cases were reported. Ko and Anand [11] carried out an experiment for turbulent channel flow with porous baffles and found that the case with porous baffles presents a flow behavior as good as the one with solid baffles. Tsay et al. [12] investigated numerically by using baffles for enhancement of heat transfer in laminar channel flow over two heated blocks mounted on the lower plate.The concept of periodically fully developed flow was widely used. However, most investigations studied only effects of turbulators and their parameters. Therefore, the study on developing of periodic flow structure concept has rarely been reported. In this present work, the numerical computations for three-dimensional laminar periodic channel flows over a 45° inclined baffle are conducted to examine the changes in the flow structure and a way to develop to fully developed periodic flow profile.


## 2. Flow Descriptions

### 2.1. Baffle Geometry and Arrangement

The system of interest is a horizontal circular tube with inclined baffles which are repeatedly inserted at the middle of the test tube as shown in Figure1(a). The circular tube is divided into 3 parts: entry region, test section, and exit. The details of the baffled tube and computation domain and its grid are shown in Figures 1(b) and 1(c), respectively. The air enters the tube at an inlet temperature, T
in, and flows over baffle where b is the baffle height, D, set to 0.05 m, is the tube diameter, and b/D is known as the blockage ratio, BR = 0.10–0.25. The distance between the inclined baffle cell or axial pitch,L, is set to L = D in which L/D is defined as the pitch spacing ratio, PR, and set to 1. To investigate an effect of Reynolds number, Re= 100–1000.(a) Circular tube geometry with 45° inclined baffle inserted (b) baffled tube details, and (c) computational domain.
(a)
(b)
(c)


### 2.2. Boundary Conditions

For baffled circular tube, a uniform air velocity is introduced at the inlet while a pressure outlet condition is applied at the exit. The physical properties of the air have been assumed to be constant at mean bulk temperature. Impermeable boundary and no-slip wall conditions have been implemented over the tube walls as well as the baffle.


## 2.1. Baffle Geometry and Arrangement

The system of interest is a horizontal circular tube with inclined baffles which are repeatedly inserted at the middle of the test tube as shown in Figure1(a). The circular tube is divided into 3 parts: entry region, test section, and exit. The details of the baffled tube and computation domain and its grid are shown in Figures 1(b) and 1(c), respectively. The air enters the tube at an inlet temperature, T
in, and flows over baffle where b is the baffle height, D, set to 0.05 m, is the tube diameter, and b/D is known as the blockage ratio, BR = 0.10–0.25. The distance between the inclined baffle cell or axial pitch,L, is set to L = D in which L/D is defined as the pitch spacing ratio, PR, and set to 1. To investigate an effect of Reynolds number, Re= 100–1000.(a) Circular tube geometry with 45° inclined baffle inserted (b) baffled tube details, and (c) computational domain.
(a)
(b)
(c)


## 2.2. Boundary Conditions

For baffled circular tube, a uniform air velocity is introduced at the inlet while a pressure outlet condition is applied at the exit. The physical properties of the air have been assumed to be constant at mean bulk temperature. Impermeable boundary and no-slip wall conditions have been implemented over the tube walls as well as the baffle.


## 3. Mathematical Foundation

The numerical model for fluid flow in a circular tube was developed under the following assumptions.(i)
Steady, three-dimensional, laminar, and incompressible fluid flows.
(ii)
Constant fluid properties.
(iii)
Ignoring body forces.Based on the above assumptions, the tube flow is governed by the continuity and the Navier-Stokes equations. In the Cartesian tensor system, these equations can be written as follows.Continuity equation:(1)
∂
∂
x
i
(
ρ
u
i
)
=
0
.Momentum equation:(2)
∂
(
ρ
u
i
u
j
)
∂
x
j
=
-
∂
p
∂
x
i
+
∂
∂
x
j
[
μ
(
∂
u
i
∂
x
j
+
∂
u
j
∂
x
i
)
]
.Apart from the energy equation discretized by the QUICK scheme, the governing equations were discretized by the power law scheme, decoupling with the SIMPLE algorithm and solved using a finite volume approach [13]. The solutions were considered to be converged when the normalized residual values were less than 10−5 for all variables.The Reynolds number is defined as(3)
Re
=
ρ
u
-
D
μ
.The friction factor,f, is computed by pressure drop, Δ
p, across the length of the periodic channel, L, as
(4)
f
=
(
Δ
p
/
L
)
D
(
1
/
2
)
ρ
u
-
2
.The variation inf  values for the 45° inclined baffle at BR = 0.20 and Re = 600 is less than 0.2% with increasing the number of cells from 500,000 to 800,000. Therefore, there is no such advantage in increasing the number of cells beyond this value and, thus, the grid system of 500,000 cells is adopted for the current computation.


## 4. Results and Discussion

The results of 45° inclined baffle with inline arrangement in circular tested tube are expressed in six parts, such as validation, fully developed periodic concept, effect of transverse plane position, flow structure, effect of Reynolds number, and effect of blockage ratio.


### 4.1. Validation

Verification for the friction factor of the smooth circular tube with no baffle is initially performed by a comparison with the values from previous correlations under a similar operating condition as shown in Figure2. The present numerical smooth duct result shows the excellent agreement with correlation solutions obtaining from the open literature [14] for friction factor, within ±1.5% maximum deviations.Figure 2
Verification of friction factor for smooth circular tube with no baffle.


### 4.2. Fully Developed Periodic Concept

Understanding of a fully developed periodic profile condition in the baffled tube is needed before discussing the results. The fully developed periodic flow conditions in the tube mounted repeatedly with inline 45° inclined baffles can be displayed by considering the axialu/u
0 distribution as depicted in Figures 4 and 5. The u/u
0 distributions of a baffled test tube are presented at different transverse plane positions as Figure 3.Figure 3
YZ sxis of baffled tube.Figure 4
Velocity profile atz/D = 0 with various y/D value for BR = 0.20 and Re = 600.Figure 5
Velocity profile aty/D = 0.25 with various z/D value for BR = 0.20 and Re = 600.The profile of velocity is explained into 2 forms. First, the velocity profile that structures similarity to other modules but the values are not equal as shown in Figures4 and 5 at x/D = 3-4 and x/D = 2-3, respectively; this profile is called “periodic flow profile”. Second, the profile that equal on both structures and values is called “fully developed periodic flow profile” as Figures 4 and 5 when the x/D ≥ 7.


### 4.3. Effect of Transverse Plane Position

The effect of transverse plane position is expressed into 2 axes, such asy/D and z/D, as seen in Figures 4 and 5, respectively.Figure4 displays the axial u/u
0 profiles along the baffled tube with BR = 0.2 and Re = 600 at location z/D= 0 with different y/D values of 0.02, 0.25 and 0.48. It can be seen in the figure that the u/u
0 profiles for all cases become periodic at the 2nd module and tend to increase to a fully developed periodic flow at about the 6th-7th module or at x/D
≈ 6-7. In addition, theu/u
0 profile at the location y/D= 0.02 close to baffle becomes a fully developed periodic flow profile slightly faster than that at y/D = 0.25 and 0.48.Figure5 shows the axial distribution of u/u
0 for y/D = 0.25 with z/D values of −0.25, 0 and 0.25. In this figure, it is apparent that the u/u
0 profiles for different z/D locations at a fixed location of y/D = 0.25 have no effect for fully developed periodic flow regions while those for different y/D locations at a fixed point of z/D = 0.25 show different locations to obtain the fully developed periodic flow regime. In closer inspection, the flow in the region closer to the baffle becomes a fully developed periodic profile faster than that farther away from the baffle.Therefore, the concept of fully developed periodic flow profile can be applied efficiently to laminar tube flow through baffles if the test tube is sufficiently long (x/D > 10 is preferable, in general). It can be observed that the mean u/u
0 value is not much different for both the periodic regions. Further, with considering both convergent time and solution precision, only a fully developed periodic flow tube model is employed in the computation.


### 4.4. Flow Structure

The flow structure in the tube repeatedly mounted with the inline 45° inclined baffles can be displayed by considering the streamline in transverse planes for Re = 600 and BR = 0.20 as shown in Figures6–9. By partially describing the flow structure into 2 parts, the developing flow structure and the fully developed periodic flow profile.Figure 6
Streamline in transverse plane for BR = 0.20 at Re = 600.Figure 7
Streamline in transverse planes m1 to m8 for BR = 0.20 at Re = 600.Figure 8
Streamline in transverse plane for a module at BR = 0.20 at Re = 600.Figure 9
Streamline in transverse planes A1 to A5 for a module at BR = 0.20 at Re = 600.The concept of fully developed periodic flow profile is shown in Figures6 and 7 for overall streamline planes and the details of planes m1 to m8, respectively. Streamline in transverse planes of each module is plotted at a similar position, in the middle of every element. In Figure 6 the flow profiles are found to be similar as planes m4 to m14. The detail of streamline in transverse planes is shown in Figure 7 for only m1 to m8; causes m9 to m14 give similar structure as depicted in Figure 6. It is visible that the plane m1and the two counterrotating vortex pairs caused by the baffle appear on the lower and upper parts of the module. The core of pair vortex appears at the left of the plane from the middle of plane (z/D = 0). The next module of plane m2 and the two counterrotating vortex appears similarly to m1 but the position of vortex core shows a difference as seen at the right of the plane. Planes m3 to m8 show similarity on both flow and vortex core position, at the middle of the plane and on the baffle position. In addition, the result from the streamline in the transverse plane plotted can be describing flow structure and periodic flow profile but the fully developed periodic flow profile must be explained by the values of velocity.The change of flow structure in fully developed periodic flow profile (Figure8) from the start to the end of module at x/D = 9-10 is presented and the details of each plane, A1 to A5, are shown in Figure 9. Plane A1, at the middle of the baffle point, two counterrotating vortexes from left to right, and vortex core appear in the middle of the plane over the baffle; after that, the vortex core moves along the baffle to plane A2 and appear the new pairs of vortex on the left of the plane. Plane A3 is the end of baffle at the right plane and start point of baffle at the left plane; vortex core before gradually faded and the new one will be appear. Plane A4, the new core of vortex, moves along the baffle as well as next plane, A5. The flow structure presents repeatedly from plane A1 to A5 in the fully developed periodic flow profile region.


### 4.5. Effect of Re

The effect or Re is described by consideringu/u
0 profile versus x/D at different Re as depicted in Figure 10; the inclined baffles with BR = 0.20 at z/D = 0 and y/D = 0.25 for Re in range from 100 to 1000 are presented. In this figure, the flow structure along the x-axis gives similarity structure for all Re values. The flow characteristic begins periodic flow at x/D
≈ 2-3 and becomes a fully developed periodic flow profile at x/D
≈ 6-7 for all cases. This means that the Re is not effective for the developing of fully developed periodic flow structure.Figure 10
Velocity profile atz/D = 0 and y/D = 0.25 with various Re values for BR = 0.20.


### 4.6. Effect of BR

The effect of BR can be displayed by considering the velocity profile,u
i/u
0, as depicted in Figures 11(a), 11(b), 11(c), and 11(d) for BR = 0.10, 0.15, 0.20, and 0.25 at Re = 600, y/D = 0.25, and z/D = −0.25, respectively.Velocity profile aty/D = 0.25, z/D = −0.25 and Re = 600 for (a) BR = 0.10, (b) BR = 0.15, (c) BR = 0.20, and (d) BR = 0.25.
(a)
(b)
(c)
(d)In Figure11(a), the v/u
0 profile shows periodic profile at x/D
≈ 2-3 and becomes fully developed periodic flow profile at x/D
≈ 9-10. The u/u
0 profile provides periodic profile at x/D
≈ 5-6 but does not become fully developed periodic flow structure. The w/u
0 profile shows periodic profile and fully developed periodic flow profile at x/D
≈ 2-3 and x/D
≈ 11-12, respectively. In addition, the flow structure of BR = 0.10 does not seem to be going to be fully developed periodic flow profile in range of 14 modules tested. Fully developed periodic flow concept for BR ≤ 0.10 will appear when increasing the module of test section or decreasing the distance between baffle, PR < 1.In Figures11(b), 11(c), and 11(d), BR = 0.15–0.25; u/u
0 profile presents periodic at x/D
≈ 2-3 and becomes fully developed periodic at x/D
≈ 6-7 when v/u
0 and w/u
0 develop into fully developed periodic flow profile slightly faster than x-velocity around a half to one module. This can be attributed that the fully developed periodic flow profiles perform faster when increasing BR.


## 4.1. Validation

Verification for the friction factor of the smooth circular tube with no baffle is initially performed by a comparison with the values from previous correlations under a similar operating condition as shown in Figure2. The present numerical smooth duct result shows the excellent agreement with correlation solutions obtaining from the open literature [14] for friction factor, within ±1.5% maximum deviations.Figure 2
Verification of friction factor for smooth circular tube with no baffle.


## 4.2. Fully Developed Periodic Concept

Understanding of a fully developed periodic profile condition in the baffled tube is needed before discussing the results. The fully developed periodic flow conditions in the tube mounted repeatedly with inline 45° inclined baffles can be displayed by considering the axialu/u
0 distribution as depicted in Figures 4 and 5. The u/u
0 distributions of a baffled test tube are presented at different transverse plane positions as Figure 3.Figure 3
YZ sxis of baffled tube.Figure 4
Velocity profile atz/D = 0 with various y/D value for BR = 0.20 and Re = 600.Figure 5
Velocity profile aty/D = 0.25 with various z/D value for BR = 0.20 and Re = 600.The profile of velocity is explained into 2 forms. First, the velocity profile that structures similarity to other modules but the values are not equal as shown in Figures4 and 5 at x/D = 3-4 and x/D = 2-3, respectively; this profile is called “periodic flow profile”. Second, the profile that equal on both structures and values is called “fully developed periodic flow profile” as Figures 4 and 5 when the x/D ≥ 7.


## 4.3. Effect of Transverse Plane Position

The effect of transverse plane position is expressed into 2 axes, such asy/D and z/D, as seen in Figures 4 and 5, respectively.Figure4 displays the axial u/u
0 profiles along the baffled tube with BR = 0.2 and Re = 600 at location z/D= 0 with different y/D values of 0.02, 0.25 and 0.48. It can be seen in the figure that the u/u
0 profiles for all cases become periodic at the 2nd module and tend to increase to a fully developed periodic flow at about the 6th-7th module or at x/D
≈ 6-7. In addition, theu/u
0 profile at the location y/D= 0.02 close to baffle becomes a fully developed periodic flow profile slightly faster than that at y/D = 0.25 and 0.48.Figure5 shows the axial distribution of u/u
0 for y/D = 0.25 with z/D values of −0.25, 0 and 0.25. In this figure, it is apparent that the u/u
0 profiles for different z/D locations at a fixed location of y/D = 0.25 have no effect for fully developed periodic flow regions while those for different y/D locations at a fixed point of z/D = 0.25 show different locations to obtain the fully developed periodic flow regime. In closer inspection, the flow in the region closer to the baffle becomes a fully developed periodic profile faster than that farther away from the baffle.Therefore, the concept of fully developed periodic flow profile can be applied efficiently to laminar tube flow through baffles if the test tube is sufficiently long (x/D > 10 is preferable, in general). It can be observed that the mean u/u
0 value is not much different for both the periodic regions. Further, with considering both convergent time and solution precision, only a fully developed periodic flow tube model is employed in the computation.


## 4.4. Flow Structure

The flow structure in the tube repeatedly mounted with the inline 45° inclined baffles can be displayed by considering the streamline in transverse planes for Re = 600 and BR = 0.20 as shown in Figures6–9. By partially describing the flow structure into 2 parts, the developing flow structure and the fully developed periodic flow profile.Figure 6
Streamline in transverse plane for BR = 0.20 at Re = 600.Figure 7
Streamline in transverse planes m1 to m8 for BR = 0.20 at Re = 600.Figure 8
Streamline in transverse plane for a module at BR = 0.20 at Re = 600.Figure 9
Streamline in transverse planes A1 to A5 for a module at BR = 0.20 at Re = 600.The concept of fully developed periodic flow profile is shown in Figures6 and 7 for overall streamline planes and the details of planes m1 to m8, respectively. Streamline in transverse planes of each module is plotted at a similar position, in the middle of every element. In Figure 6 the flow profiles are found to be similar as planes m4 to m14. The detail of streamline in transverse planes is shown in Figure 7 for only m1 to m8; causes m9 to m14 give similar structure as depicted in Figure 6. It is visible that the plane m1and the two counterrotating vortex pairs caused by the baffle appear on the lower and upper parts of the module. The core of pair vortex appears at the left of the plane from the middle of plane (z/D = 0). The next module of plane m2 and the two counterrotating vortex appears similarly to m1 but the position of vortex core shows a difference as seen at the right of the plane. Planes m3 to m8 show similarity on both flow and vortex core position, at the middle of the plane and on the baffle position. In addition, the result from the streamline in the transverse plane plotted can be describing flow structure and periodic flow profile but the fully developed periodic flow profile must be explained by the values of velocity.The change of flow structure in fully developed periodic flow profile (Figure8) from the start to the end of module at x/D = 9-10 is presented and the details of each plane, A1 to A5, are shown in Figure 9. Plane A1, at the middle of the baffle point, two counterrotating vortexes from left to right, and vortex core appear in the middle of the plane over the baffle; after that, the vortex core moves along the baffle to plane A2 and appear the new pairs of vortex on the left of the plane. Plane A3 is the end of baffle at the right plane and start point of baffle at the left plane; vortex core before gradually faded and the new one will be appear. Plane A4, the new core of vortex, moves along the baffle as well as next plane, A5. The flow structure presents repeatedly from plane A1 to A5 in the fully developed periodic flow profile region.


## 4.5. Effect of Re

The effect or Re is described by consideringu/u
0 profile versus x/D at different Re as depicted in Figure 10; the inclined baffles with BR = 0.20 at z/D = 0 and y/D = 0.25 for Re in range from 100 to 1000 are presented. In this figure, the flow structure along the x-axis gives similarity structure for all Re values. The flow characteristic begins periodic flow at x/D
≈ 2-3 and becomes a fully developed periodic flow profile at x/D
≈ 6-7 for all cases. This means that the Re is not effective for the developing of fully developed periodic flow structure.Figure 10
Velocity profile atz/D = 0 and y/D = 0.25 with various Re values for BR = 0.20.


## 4.6. Effect of BR

The effect of BR can be displayed by considering the velocity profile,u
i/u
0, as depicted in Figures 11(a), 11(b), 11(c), and 11(d) for BR = 0.10, 0.15, 0.20, and 0.25 at Re = 600, y/D = 0.25, and z/D = −0.25, respectively.Velocity profile aty/D = 0.25, z/D = −0.25 and Re = 600 for (a) BR = 0.10, (b) BR = 0.15, (c) BR = 0.20, and (d) BR = 0.25.
(a)
(b)
(c)
(d)In Figure11(a), the v/u
0 profile shows periodic profile at x/D
≈ 2-3 and becomes fully developed periodic flow profile at x/D
≈ 9-10. The u/u
0 profile provides periodic profile at x/D
≈ 5-6 but does not become fully developed periodic flow structure. The w/u
0 profile shows periodic profile and fully developed periodic flow profile at x/D
≈ 2-3 and x/D
≈ 11-12, respectively. In addition, the flow structure of BR = 0.10 does not seem to be going to be fully developed periodic flow profile in range of 14 modules tested. Fully developed periodic flow concept for BR ≤ 0.10 will appear when increasing the module of test section or decreasing the distance between baffle, PR < 1.In Figures11(b), 11(c), and 11(d), BR = 0.15–0.25; u/u
0 profile presents periodic at x/D
≈ 2-3 and becomes fully developed periodic at x/D
≈ 6-7 when v/u
0 and w/u
0 develop into fully developed periodic flow profile slightly faster than x-velocity around a half to one module. This can be attributed that the fully developed periodic flow profiles perform faster when increasing BR.


## 5. Conclusions

3D laminar fully developed periodic flow characteristics in a circular test tube mounted repeatedly with inline 45° inclined baffles inserted in the middle of the test tube for 14 modules are numerically investigated.The baffled tube flow shows periodic flow profile atx/D
≈ 2-3 and becomes fully developed periodic flow profiles at x/D
≈ 6-7 downstream of the tested tube inlet for BR = 0.15–0.25, PR = 1, and Re = 100–1000 when BR = 0.10 does not become fully developed periodic flow. The fully developed periodic flow profile concept depends on BR and transverse plane position. the increasing BR and close to the baffle regimes are the factors for occurring the fully developed periodic flow profile. The Re in range of study has no effect for fully developed periodic flow profile.Thus, the fully developed periodic flow structure and behavior concept in tested tube at BR ≥ 0.15, PR = 1, and Re = 100–1000 are recommended for use instead of the full tube and all of systems in order to save more resources and time in numerical investigation. In case of BR ≤ 0.10, the fully developed periodic flow profile concept may be used for PR < 1 and the tested section is higher than 14 modules.



---

*Source: 735905-2014-03-05.xml*