Journal of Thermodynamics

Journal of Thermodynamics / 2014 / Article

Research Article | Open Access

Volume 2014 |Article ID 670129 | 10 pages | https://doi.org/10.1155/2014/670129

Thermohydraulic Performance of a Series of In-Line Noncircular Ducts in a Parallel Plate Channel

Academic Editor: Shou-Shing Hsieh
Received28 May 2014
Revised12 Aug 2014
Accepted20 Aug 2014
Published21 Sep 2014

Abstract

Heat transfer and fluid flow characteristics for two-dimensional laminar flow at low Reynolds number for five in-line ducts of various nonconventional cross-sections in a parallel plate channel are studied in this paper. The governing equations were solved using finite-volume method. Commercial CFD software, ANSYS Fluent 14.5, was used to solve this problem. A total of three different nonconventional, noncircular cross-section ducts and their characteristics are compared with those of circular cross-section ducts. Shape-2 ducts offered minimum flow resistance and maximum heat transfer rate most of the time. Shape-3 ducts at Re < 100 and Shape-2 ducts at Re > 100 can be considered to give out the optimum results.

1. Introduction

The heat transfer enhancement in most of the engineering applications is a never ending process. The need for better heat transfer rate and low flow resistance has led to extensive research in the field of heat exchangers. Higher heat transfer rate and low pumping power are desirable properties of a heat exchanger. The duct shape and its arrangement highly influence flow characteristics in a heat exchanger. Flow past cylinders, especially circular, flat, oval, and diamond arranged in a parallel plate channel, were extensively studied by Bahaidarah et al. [13]. They carried domain discretization in body-fitted coordinate system while the governing equations were solved using a finite-volume technique. Chhabra [4] studied bluff bodies of different shapes like circle, ellipse, square, semicircle, equilateral triangle, and square submerged in non-Newtonian fluids. Kundu et al. [5, 6] investigated fluid flow and heat transfer coefficient experimentally over a series of in-line circular cylinders in parallel plates using two different aspect ratios for intermediate range of Re 220 to 2800. Grannis and Sparrow [7] obtained numerical solutions for the fluid flow in a heat exchanger consisting of an array of diamond-shaped pin fins. Implementation of the model was accomplished using the finite element method. Tanda [8] performed experiments on fluid flow and heat transfer for a rectangular channel with arrays of diamond shaped elements. Both in-line and staggered fin arrays were considered in his study. Jeng [9] experimentally investigated pressure drop and heat transfer of an in-line diamond shaped pin-fin array in a rectangular duct. Terukazu et al. [10] studied the heat transfer characteristics and flow behaviours around an elliptic cylinder at high Reynolds number. Gera et al. [11] numerically investigated a two-dimensional unsteady flow past a square cylinder for the Reynolds number (Re) considered in the range of 50–250. The features of the flow past the square cylinder were observed using CFD. Olawore and Odesola [12] numerically investigated two-dimensional unsteady flow past a rectangular cylinder. The effect of vortical structure and pressure distribution around the section of rectangular cylinders are studied in their work. Chen et al. [13, 14] analysed flow and conjugate heat transfer in a high-performance finned oval tube heat exchanger element and calculated them for a thermally and hydrodynamically developing three-dimensional laminar flow. Computations were performed with a finite-volume method based on the SIMPLEC algorithm.

Zdravistch et al. [15] numerically predicted fluid flow and heat transfer around staggered and in-line tube banks and found out close agreement experimental test cases. Tahseen et al. [16] conducted numerical study of the two-dimensional forced convection heat transfer across three in-line flat tubes confined in a parallel plate channel, the flow under incompressible and steady-state conditions. They solved the system in the body-fitted coordinates (BFC) using the finite-volume method (FVM). Gautier et al. [17] proposed a new set of boundary conditions to improve the representation of the infinite flow domain.

Conventional tube/duct shapes have been investigated for fluid flow and heat transfer characteristics. In this paper, two-dimensional, steady state, laminar flow over nonconventional noncircular (Shape-1, Shape-2, and Shape-3) in-line ducts confined in a parallel plate channel is considered. The different geometries investigated are Circular, Shape-1, Shape-2, and Shape-3. Only representative cases are discussed in this paper.

2. Mathematical Formulation

The governing mass, momentum, and energy conservation equations [1] for steady incompressible flow of a Newtonian fluid are expressed as

Commercial CFD software, ANSYS Fluent 14.5, was used to solve the governing equations. The pressure-velocity coupling is done by semi-implicit method for pressure linked equations (SIMPLE) scheme.

The grid for all the cases is generated using the Automatic Method with all Quadrilateral elements in the Mesh Component System of the ANSYS Fluent 14.5 software. In the numerical simulation, the grid of 11,000 nodes is found to be fair enough for grid independent solution. Refining the grid proved to be unnecessary, resulting in utilizing more computer resources. Figure 2 shows the grid generated for all the configurations.

3. Geometric Configuration

A total of four different duct cross-sections, namely, Circular, Shape-1, Shape-2, and Shape-3 were analysed. The five in-line ducts were confined in a parallel plate channel. The distance between two consecutive ducts was kept constant. , , and ratio (shown in Figure 1) kept constant for all cross-sections of the ducts. The height of obstruction is kept constant; that is, the diameter () of circular cross-section duct is taken for all the other cross-sections (see Figure 1). The unobstructed length of the pipe before the first and the last module is kept as one module length and three module lengths, respectively, so that the flow is fully developed. It is also ensured that the flow is fully developed when it enters the channel. The ducts and the parallel plate channel walls were assumed to be of infinite extent in the -direction, that is, perpendicular to the paper. Hence, the flow could be considered as two-dimensional. Since, the geometry is symmetric along the -axis, only the lower half portion is considered for the computation purpose.

4. Boundary Conditions

Figure 3 represents the schematic of the boundary conditions for all the cross-sections used. A fully developed flow with velocity profile () and with temperature () was assigned for the flow at inlet. The temperature was kept different from that of the plate and duct walls. A no slip () boundary condition was given to the plate and duct surface. The surface temperatures of the plates and the ducts were taken to be constant (). At the symmetry, normal component of velocity and normal gradient of other velocity components were taken as zero. At the outflow boundary, there is no change in velocity across boundary.

5. Validation

The parallel plate channel problem corresponds to the rectangular duct with infinite aspect ratio. The calculated value of Nusselt number is 7.541. Comparing the calculated Nusselt number with that of the theoretical Nusselt number value of 7.54 [18] shows that our computed values are in good agreement with those of the theoretical values.

Using the geometric parameters and for circular cross-section, the normalized pressure drop (), and Nusselt number for third module of the parallel plate channel, when compared with the values of Kundu et al. [5], shows good agreement between both the values.

6. Results and Discussions

The numerical simulation was carried out for Reynolds number ranging from 25 to 350 and Prandtl number was taken as 0.74.

6.1. Onset of Recirculation

For any duct shape, no separation of flow was observed below Re = 30. Circular duct displayed the signs of recirculation at the lowest Re value (Re = 32) [1]. The Shape-1 ducts displayed the onset of recirculation at Re value of 48; Shape-2 ducts showed the signs of onset of recirculation at Re = 60, while the Shape-3 showed the signs of recirculation just under Re = 40.

As the values of Reynolds number are increased, the recirculation region starts increasing. At Re = 150, the length of the recirculation region is slightly smaller than the distance between the two consecutive ducts in the flow direction, while, at Re = 350, the length of the recirculation region is greater than the distance between the two consecutive ducts. Since there is no obstruction to the flow downstream of the last duct, the elliptic behaviour of the flow is observed as shown in Figures 58.

6.2. Fluid Flow Characteristics

From Figures 5, 6, 7, and 8, we can observe that the streamlines past the ducts 1–4 are similar in nature. Thus, the stream-wise velocity distributions at the module inlets (see Figure 9) and normalized pressure drop Figure 4 and module average Nusselt number for each module (Table 2) confirm the existence of periodically fully developed flow downstream of the second module. As the flow is periodically fully developed for the inner modules, it is sufficient to study one inner module so as to understand the essential of flow physics. In present study, we have selected module 3. Small variations for the first and the last modules are observed due to the end effects.

Figure 10 shows the normalized (dimensionless) pressure drop () as a function of Reynolds number. We can observe that, as the value of Re increases, there is descend in normalized pressure drop for all the ducts presented in this study.

The ratio of normalized pressure drop for noncircular duct module () to that for the circular duct module () is shown in Figure 11. The Shape-1 ducts present pressure drop almost similar to that of the Circular ducts due to their geometric similarity. Shape- 2 ducts present lowest pressure drop ratio (), which means that they offer lower resistance to flow than the circular ducts till Re = 300 resulting in lower pumping power requirement. Shape-3 ducts present lower values till Re = 150, but, after Re = 150, they increase significantly, resulting in more pumping power.

6.3. Heat Transfer Characteristics

Figure 12 presents the module average Nusselt number as a function of Re for module 3 for Shape-1, Shape-2, Shape-3, and Circular ducts Table 1. The Nusselt number for Re < 50 increases rapidly as the Re increases. Till Re = 50, there is no significant difference in Nusselt number values for different cross-section ducts, which implies that the duct shape has very little effect on heat transfer rate at Re < 50. At higher values of Re the recirculation of the flow takes place. These recirculation vortices do not transfer heat energy to the main stream of flow. Hence, their contribution to heat transfer is very less. At Re < 50, the Nu is almost the same for various cross-section ducts but, after Re = 50, the Nusselt number varies significantly with the cross-section of the duct. Circular cross-section ducts present the least Nu values till Re = 250, while Shape-2 cross-section ducts present the highest Nu values after Re = 250.


ReSecond moduleThird moduleFourth module

Kundu et al. [6]509.49.49.8
Present work9.519.489.48

Kundu et al. [6]20012.512.612.8
Present work12.6812.6312.64


ModuleCircularShape-1Shape-2Shape-3
Nu0/NuNu+Nu*/NuNu+Nu*/NuNu+Nu*

Re = 25
 18.0206.8828.0901.0097.0171.0201.0177.9850.9967.0321.0221.0237.9660.9937.0291.0211.024
 28.1186.8908.1771.0077.0151.0181.0168.0580.9937.0331.0211.0238.0610.9937.0391.0221.024
 38.1136.8838.1821.0097.0211.0201.0178.0450.9927.0241.0211.0238.0600.9947.0301.0211.024
 48.1046.8858.1741.0097.0091.0181.0158.0440.9937.0321.0211.0248.0580.9947.0351.0221.024
 58.1286.9088.1911.0086.9381.0041.0028.0840.9947.0671.0231.0258.0640.9927.0531.0211.024
Re = 50
 14.5389.4864.5250.9979.7991.0331.0344.4420.9799.8271.0361.0434.4760.9869.8321.0371.041
 24.5079.5174.4980.9989.7691.0271.0274.4190.9809.8121.0311.0384.4180.9809.8551.0361.042
 34.5109.4824.5010.9989.7971.0331.0344.4070.9779.7721.0311.0394.4340.9839.8141.0351.041
 44.5009.4794.4970.9999.7661.0301.0304.4060.9799.7771.0311.0394.4290.9849.8131.0351.041
 54.4969.5094.4820.9979.7941.0301.0314.4060.9809.7881.0291.0364.3940.9779.8511.0361.044
Re = 100
 13.20811.3403.1300.97611.7111.0331.0413.0380.94711.7111.0331.0523.1170.97211.7681.0381.048
 22.69011.2242.6921.00111.4221.0181.0172.6580.98811.4871.0231.0272.6200.97411.5701.0311.040
 32.72911.1342.7020.99011.4781.0311.0342.6310.96411.3841.0221.0352.6840.98411.4691.0301.036
 42.69711.1232.7011.00211.4121.0261.0252.6270.97411.3911.0241.0332.6640.98811.4701.0311.035
 52.66211.1722.6230.98511.4371.0241.0292.5650.96411.3831.0191.0322.5460.95611.5071.0301.046
Re = 150
 12.87112.5882.7790.96813.0091.0331.0452.6960.93912.9541.0291.0512.7740.96613.1101.0411.053
 22.02012.0702.0551.01712.2581.0161.0102.0441.01212.3451.0231.0191.9810.98112.3611.0241.031
 32.07111.9732.0580.99412.2811.0261.0282.0010.96612.1951.0191.0302.0680.99812.2531.0231.024
 42.02811.9372.0651.01812.2381.0251.0191.9920.98212.1871.0211.0272.0361.00412.2381.0251.024
 52.00811.9881.9680.98012.1881.0171.0241.9260.95912.1361.0121.0271.8960.94412.1991.0181.037
Re = 200
 12.68414.0532.6040.97014.5231.0331.0442.5320.94314.3891.0241.0442.6080.97214.7051.0461.056
 21.66712.6811.7231.03312.9671.0231.0111.7161.02913.0621.0301.0201.6430.98612.8821.0161.021
 31.70712.6301.7040.99812.8251.0151.0161.6660.97612.8521.0181.0261.7331.01512.8871.0201.015
 41.66412.5481.7261.03712.9401.0311.0191.6480.99012.7891.0191.0231.6931.01712.8331.0231.017
 51.66812.4351.6300.97812.5791.0121.0191.6060.96312.5921.0131.0261.5650.93812.5121.0061.028
Re = 250
 12.55115.7062.4910.97616.3221.0391.0482.4290.95216.1421.0281.0452.5020.98116.5821.0561.063
 21.45013.2251.5181.04713.6371.0311.0161.5061.03913.6501.0321.0191.4310.98713.3051.0061.011
 31.47413.1751.4751.00113.2071.0021.0021.4560.98813.4581.0211.0261.5181.03013.4201.0191.009
 41.43513.1251.5141.05513.6011.0361.0181.4290.99613.2951.0131.0141.4721.02613.2821.0121.003
 51.45412.5251.4250.98112.6011.0061.0131.4140.97312.7201.0161.0251.3640.93812.5401.0011.023
Re = 300
 12.45017.3522.4090.98318.2241.0501.0562.3550.96118.0241.0391.0532.4230.98918.4991.0661.070
 21.30213.7171.3771.05714.2091.0361.0171.3581.04314.1061.0281.0141.2850.98613.7030.9991.004
 31.31113.6251.3151.00313.5260.9930.9921.3121.00114.0011.0281.0271.3681.04313.8501.0171.002
 41.27713.6811.3661.07014.1981.0381.0151.2771.00013.7411.0041.0041.3161.03113.6000.9940.984
 51.30212.4351.2870.98912.5211.0071.0111.2840.98612.7131.0221.0271.2270.94212.5141.0061.026
Re = 350
 12.37018.9452.3470.99020.1241.0621.0662.2980.96919.9251.0521.0632.3580.99520.3791.0761.078
 21.19414.1491.2721.06614.6811.0381.0161.2451.04314.4711.0231.0081.1770.98614.0510.9930.998
 31.19014.0091.1961.00413.8120.9860.9851.2051.01214.4741.0331.0291.2561.05514.2171.0150.997
 41.16014.1851.2551.08214.7121.0371.0101.1631.00314.1270.9960.9951.1981.03313.8230.9750.964
 51.18612.2991.1861.00012.4551.0131.0131.1871.00112.6941.0321.0321.1240.94812.4831.0151.033

Figure 13 presents the heat transfer enhancement ratio (Nu+) as a function of Re for module 3 for Shape-1, Shape-2, and Shape-3 ducts. Nu+ is the ratio of Nusselt number for noncircular duct module (Nu) to that for the circular duct module (). It represents how effective the heat transfer rate will be using the noncircular ducts as compared to the circular ducts. For Re < 50, all the ducts show similar Nu+ values. Shape-2 ducts show least Nu+ values at Re < 200 but, after Re = 200, they rise up significantly giving highest Nu+ values as compared to their counter parts. Shape-1 ducts show constant drop in the Nu+ values and, at Re > 250, they have heat transfer rate less than the circular ducts.

Figure 14 presents the heat transfer performance ratio () as a function of Re for module 3 for Shape-1, Shape-2, and Shape-3 ducts. represents heat transfer enhancement per unit increase in pumping power. The Shape-1 ducts performed poorly as compared to Circular ducts at Re > 250. Shape-3 ducts presented good results at Re < 100 but then they presented a constant drop in the values, while Shape-2 ducts displayed promising results after Re = 100.

7. Conclusion

In an environment where the heat transfer rate and pumping power is critical, Shape-2 ducts show promising results as compared to their counterparts. The duct shape has very little influence on the heat transfer rate at low Reynolds number. For optimal design of a heat exchanger for heat transfer enhancement per unit increase in pumping power, Reynolds number and cross-section of the duct must be carefully chosen. So, for heat exchanger, Shape-3 ducts at Re < 100 and Shape-2 ducts at Re > 100 can be considered to give out the optimum results.

Nomenclature

Area of heat transfer surfaces in a module (m2)
Specific heat (J/kg K)
Blockage height of a tube inside the channel (m)
Module friction factor
Module friction factor for circular tubes
Height of the channel (m)
Module length (m)
Thermal conductivity (W/m·K)
Length of the tube cross-section along the flow direction (m)
Module average Nusselt number
Module average Nusselt number for the reference case of circular tubes
Heat transfer enhancement ratio ()
Heat transfer performance ratio
Mass flow rate (kg/s)
Pressure (Pa)
Reynolds number (VH/)
Temperature (K)
Average of module inlet bulk temperature and module outlet bulk temperature (K)
Temperature of heat transfer surfaces (K)
Velocity at the channel inlet (m/s)
Velocity field
Vertical distance from the channel bottom (m)
Pressure drop across a module (Pa)
Normalized pressure drop ()
Change in bulk temperature across a module (K)
Dynamic viscosity (N·s/m2)
Density
Del operator.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to present their sincere gratitude towards the Faculty of Mechanical Engineering in Symbiosis Institute of Technology, Pune, and Rajarshi Shahu College of Engineering, Pune.

References

  1. H. M. S. Bahaidarah, M. Ijaz, and N. K. Anand, “Numerical study of fluid flow and heat transfer over a series of in-line noncircular tubes confined in a parallel-plate channel,” Numerical Heat Transfer Part B: Fundamentals, vol. 50, no. 2, pp. 97–119, 2006. View at: Publisher Site | Google Scholar
  2. H. M. S. Bahaidarah, N. K. Anand, and H. C. Chen, “A numerical study of fluid flow and heat transfer over a bank of flat tubes,” Numerical Heat Transfer A, vol. 48, no. 4, pp. 359–385, 2005. View at: Publisher Site | Google Scholar
  3. H. M. S. Bahaidarah, A numerical study of heat and momentum transfer over a bank of flat tubes [Ph.D. thesis], Texas A&M University, College Station, Tex, USA, 2004.
  4. R. P. Chhabra, “Fluid flow and heat transfer from circular and noncircular cylinders submerged in non-newtonian liquids,” Advances in Heat Transfer, vol. 43, pp. 289–417, 2011. View at: Publisher Site | Google Scholar
  5. D. Kundu, A. Haji-Sheikh, and D. Y. S. Lou, “Heat transfer predictions in cross flow over cylinders between two parallel plates,” Numerical Heat Transfer; Part A: Applications, vol. 19, no. 3, pp. 361–377, 1991. View at: Publisher Site | Google Scholar
  6. D. Kundu, A. Haji-Sheikh, and D. Y. S. Lou, “Heat transfer in crossflow over cylinders between two parallel plates,” Journal of Heat Transfer, vol. 114, no. 3, pp. 558–564, 1992. View at: Publisher Site | Google Scholar
  7. V. B. Grannis and E. M. Sparrow, “Numerical simulation of fluid flow through an array of diamond-shaped pin fins,” Numerical Heat Transfer A: Applications, vol. 19, no. 4, pp. 381–403, 1991. View at: Publisher Site | Google Scholar
  8. G. Tanda, “Heat transfer and pressure drop in a rectangular channel with diamond-shaped elements,” International Journal of Heat and Mass Transfer, vol. 44, no. 18, pp. 3529–3541, 2001. View at: Publisher Site | Google Scholar
  9. T.-M. Jeng, “Thermal performance of in-line diamond-shaped pin fins in a rectangular duct,” International Communications in Heat and Mass Transfer, vol. 33, no. 9, pp. 1139–1146, 2006. View at: Publisher Site | Google Scholar
  10. O. Terukazu, N. Hideya, and T. Yukiyasu, “Heat transfer and flow around an elliptic cylinder,” International Journal of Heat and Mass Transfer, vol. 27, no. 10, pp. 1771–1779, 1984. View at: Publisher Site | Google Scholar
  11. B. Gera, P. K. Sharma, and R. K. Singh, “CFD analysis of 2D unsteady flow around a square cylinder,” International Journal of Applied Engineering Research, vol. 1, no. 3, pp. 602–610, 2010. View at: Google Scholar
  12. A. S. Olawore and I. F. Odesola, “2D flow around a rectangular cylinder: a computational study,” International Journal of Science and Technology, vol. 2, no. 1, pp. 1–26, 2013. View at: Google Scholar
  13. Y. Chen, M. Fiebig, and N. K. Mitra, “Conjugate heat transfer of a finned oval tube part A: flow patterns,” Numerical Heat Transfer A: Applications, vol. 33, no. 4, pp. 371–385, 1998. View at: Publisher Site | Google Scholar
  14. Y. Chen, M. Fiebig, and N. K. Mitra, “Conjugate heat transfer of a finned oval tube part B: heat transfer behaviors,” Numerical Heat Transfer A, vol. 33, no. 4, pp. 387–401, 1998. View at: Publisher Site | Google Scholar
  15. F. Zdravistch, C. A. Fletcher, and M. Behnia, “Numerical laminar and turbulent fluid flow and heat transfer predictions in tube banks,” International Journal of Numerical Methods for Heat and Fluid Flow, vol. 5, no. 8, pp. 717–733, 1995. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  16. T. A. Tahseen, M. Ishak, and M. M. Rahman, “A numerical study of forced convection heat transfer over a series of flat tubes between parallel plates,” Journal of Mechanical Engineering and Sciences, vol. 3, pp. 271–280, 2012. View at: Google Scholar
  17. R. Gautier, D. Biau, and E. Lamballais, “A reference solution of the flow over a circular cylinder at Re=40,” Computers & Fluids, vol. 75, pp. 103–111, 2013. View at: Publisher Site | Google Scholar
  18. Y. A. Cengel and A. J. Ghajar, Heat and Mass Transfer: Fundamentals and Applications, Tata McGraw-Hill, 4th edition, 2011.

Copyright © 2014 Siddharth D. Mhaske et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

785 Views | 384 Downloads | 0 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at help@hindawi.com to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19.