<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">JECTC</journal-id><journal-title-group><journal-title>Journal of Electronics Cooling and Thermal Control</journal-title></journal-title-group><issn pub-type="epub">2162-6162</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jectc.2016.63010</article-id><article-id pub-id-type="publisher-id">JECTC-70463</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Engineering</subject><subject> Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Numerical Study of Fluid Flow over Bundle Tubes
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Majid</surname><given-names>Almas</given-names></name><xref ref-type="aff" rid="aff1"><sub>1</sub></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>Department of Mechanical and Materials Engineering, Florida International University, Miami, FL, USA</addr-line></aff><author-notes><corresp id="cor1">* E-mail:</corresp></author-notes><pub-date pub-type="epub"><day>08</day><month>09</month><year>2016</year></pub-date><volume>06</volume><issue>03</issue><fpage>109</fpage><lpage>119</lpage><history><date date-type="received"><day>June</day>	<month>6,</month>	<year>2016</year></date><date date-type="rev-recd"><day>Accepted:</day>	<month>September</month>	<year>5,</year>	</date><date date-type="accepted"><day>September</day>	<month>8,</month>	<year>2016</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  In this paper, a bank of tubes containing a flowing fluid which is immersed in a cross flow second medium of fluid with different temperature has been studied numerically using computational fluid dynamics. Laminar steady flow with a low Reynolds number has been studied in this work. Inlet mass flow rate and the bulk temperature are known and numerical method has been implemented to study the convective heat transfer to investigate the temperature and flow fields. Effects of different inlet bulk temperatures and mass flow rates have been investigated on temperature and pressure variations.
 
</p></abstract><kwd-group><kwd>Heat Transfer</kwd><kwd> Bundle Tube</kwd><kwd> Mass Flow Rate</kwd><kwd> Numerical Solution</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Periodic heat flows have broad industrial and engineering applications in steam generation in a boiler, condenser or air cooling in the coil of an air conditioner. Forced convection was broadly employed in different kinds of heat exchanger, and in many other applications in which the convective heat transfer has been studied and discussed. Theoretically calculation of heat transfer in a fully developed laminar flow in smooth pipe with uniform heat flux at the wall showed that the Nusselt number was constant [<xref ref-type="bibr" rid="scirp.70463-ref1">1</xref>] . A correlation for heat transfer forced convection has been proposed by Colburn [<xref ref-type="bibr" rid="scirp.70463-ref2">2</xref>] in which the comparison has been performed with fluid friction. The heat transfer of liquids in tubes has been studied experimentally by Sieder and Tate [<xref ref-type="bibr" rid="scirp.70463-ref3">3</xref>] . They derived the heat transfer correlation with viscosity ratio caused by heat transfer temperature difference. For the optimal design and operation of dry cooling systems to make clear the heat and mass transfer characteristics of air-cooled heat exchangers the demand on the performance of heat exchangers and the need to enhance their performance have been increased and this matter has been discussed by Gough [<xref ref-type="bibr" rid="scirp.70463-ref4">4</xref>] . In power plants where there is a shortage of water finned bundle tubes are frequently used to the air-cooled heat exchangers of the dry cooling system. The total thermal resistance for these kinds of heat exchangers includes three parts: the air-side convective thermal resistance, the wall conductive thermal resistance and the liquid-side convective thermal resistance. A lot of numerical and experimental studies have been done for improving the heat transfer performance of compact heat exchangers using slotted fins [<xref ref-type="bibr" rid="scirp.70463-ref5">5</xref>] and [<xref ref-type="bibr" rid="scirp.70463-ref6">6</xref>] . Ribatski and Jacobi [<xref ref-type="bibr" rid="scirp.70463-ref7">7</xref>] published a review of studies on tube bundle. In this study, a bank of tubes containing a flowing fluid at one temperature that is immersed in a second fluid in cross flow at a different temperature has been studied numerically. Effects of different inlet bulk temperatures and mass flow rate have been investigated on temperature and pressure variation.</p></sec><sec id="s2"><title>2. Problem Definition</title><p><xref ref-type="fig" rid="fig1">Figure 1</xref> shows the schematic of the 2 dimensional section of a tube bank which has been studied in this paper. The bank of tubes with a diameter of 1 cm is staggered across the cross-fluid flow. They are arranged with uniform space. In the x and y-directions their centers are separated by a distance of 2 and 1cm respectively. The depth of the tube bank is 1 m. Because of symmetry of tube bank geometry, only a portion of the domain has been considered for simulation. The computational domain and the grid generation are shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>. The inlet mass rate is 0.05 kg/s at the inlet</p><fig-group id="fig1"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Geometry and grid generation of the problem.</title></caption><fig id ="fig1_1"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1520084x2.png"/></fig><fig id ="fig1_2"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1520084x3.png"/></fig></fig-group><p>boundary of the periodic module. The temperature of the tube wall (T<sub>wall</sub>) is 400 K and the bulk temperature of the cross flow water (T<sub>bulk</sub>) is 300 K. The physical water properties are also mentioned in <xref ref-type="fig" rid="fig1">Figure 1</xref>. A hybrid mesh has been used for the simulation. For the regions surrounding tube walls Quadrilateral cells and for the rest of the domain triangular cells have been used. Second order upwind has been used for both momentum and energy equations. The convergence criteria for continuity, momentum and energy are 10e−3, 10e−3 and 10e−6 respectively.</p></sec><sec id="s3"><title>3. Results and Discussions</title><p><xref ref-type="fig" rid="fig2">Figure 2</xref> shows the contours of pressure inlet bulk temperature of 300 K and mass flow rate of 0.05 kg/s. As shown in the figure, pressure has its maximum value where the inlet mass flow hit the front face of each tube and this value reduces as the flow moves downstream of the tube. The maximum value at the front side is 0.08 Pa and this value is at its lowest at the top and bottom of the tubes as seen in the paper and is around 0.006 Pa behind each tube.</p><p><xref ref-type="fig" rid="fig3">Figure 3</xref> depicts the contours of temperature at the inlet bulk temperature of 300 K</p><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> Pressure contour for inlet bulk temperature of 300 K and mass flow rate of 0.05 kg/s</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1520084x4.png"/></fig><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> Temperature contour for inlet bulk temperature of 300 K and mass flow rate of 0.05 kg/s</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1520084x5.png"/></fig><p>and mass flow rate of 0.05 kg/s. As shown in the figure, temperature has its maximum value as the initial surface temperature of each tube is 400 K which more than the inlet bulk temperature, 300 K. As the inlet flow passes over the tube the temperature of the front tubes side decreases to around 360 K. There is a pass between the tubes which has its lowest temperature value around 270 K. As can be seen the value of temperature behind the tubes reduces to 360 K.</p><p><xref ref-type="fig" rid="fig4">Figure 4</xref> displays the contours of temperature at the inlet bulk temperature of 500 K and mass flow rate of 0.05 kg/s. As shown in the figure, the initial surface temperature of each tube is 400 K which less than the inlet bulk temperature, 470 K. As the inlet flow passes over the tube the temperature of the front tubes side increases to around 360 K. There is a pass between the tubes which has its highest temperature value around 520 K. As can be seen the value of temperature behind the tubes increases to 450 K.</p><p><xref ref-type="fig" rid="fig5">Figure 5</xref> displays the contours of temperature at the inlet bulk temperature of 700 K and mass flow rate of 0.05 kg/s. In this figure the inlet bulk temperature increases even more to 700 to see how the flow field temperature increases over the tube bundles. As</p><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> Temperature contour for inlet bulk temperature of 500 K and mass flow rate of 0.05 kg/s</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1520084x6.png"/></fig><p>shown in the figure, the initial surface temperature of each tube is 400 K which less than the inlet bulk temperature, 700 K. As the inlet flow passes over the tube the temperature of the front tubes side increases to around 620 K. There is a pass between the tubes which has its highest temperature value around 760 K. As can be seen the value of temperature behind the tubes increases to 550 K.</p><p><xref ref-type="fig" rid="fig6">Figure 6</xref> demonstrates the temperature profiles at three different locations of 0.01, 0.02 and 0.03 m for the inlet bulk temperature of 300 K and mass flow rate of 0.05 kg/s. As seen in the figure at location x = 0.01 m the temperature value drops from 330 K to</p><fig id="fig5"  position="float"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Temperature contour for inlet bulk temperature of 700 K and mass flow rate of 0.05 kg/s</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1520084x7.png"/></fig><p>280 K and thereafter in increases sharply to the surface temperature of the tube, 400 K. At location x = 0.02 m, the temperature decreases from 320 K to around 280 and then increases to about 360 K and slightly reduces to 350 K. At location x = 0.03 m the temperature dramatically drops from tube surface temperature of 400 K to around 280 K and then increases to 340 K.</p><p><xref ref-type="fig" rid="fig7">Figure 7</xref> demonstrates the temperature profiles at three different locations of 0.01, 0.02 and 0.03 m for the inlet bulk temperature of 500 K and mass flow rate of 0.05 kg/s. As seen in the figure at location x = 0.01 m the temperature value increases from 470 K to 520 K and quickly drops sharply to the tube surface temperature of 400 K. At location</p><fig id="fig6"  position="float"><label><xref ref-type="fig" rid="fig6">Figure 6</xref></label><caption><title> Temperature profiles for inlet bulk temperature of 300 K and mass flow rate of 0.05 kg/s at locations 0.01, 0.02 and 0.03 m</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1520084x8.png"/></fig><fig id="fig7"  position="float"><label><xref ref-type="fig" rid="fig7">Figure 7</xref></label><caption><title> Temperature profiles for inlet bulk temperature of 500 K and mass flow rate of 0.05 kg/s at locations 0.01, 0.02 and 0.03 m</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1520084x9.png"/></fig><p>x = 0.02 m, the temperature increases from 480 K to almost 520 and then decreases to about 440 K and stays constant around this temperature. At location x = 0.03 m the temperature dramatically grows from tube surface temperature of 400 K to around 510 K and then reduces to 460 K.</p><p><xref ref-type="fig" rid="fig8">Figure 8</xref> shows the temperature profiles at three different locations of 0.01, 0.02 and 0.03 m for the inlet bulk temperature of 700 K and mass flow rate of 0.05 kg/s. As seen in the figure at location x = 0.01 m the temperature value increases from 600 K to 760 K and quickly drops sharply to the tube surface temperature of 400 K. At location x = 0.02 m, the temperature increases from 640 K to almost 760 and then decreases to about 540 K and stays constant around this temperature. At location x = 0.03 m the temperature dramatically grows from tube surface temperature of 400 K to around 750 K and then decreases to 600 K.</p><p>Here in <xref ref-type="fig" rid="fig9">Figure 9</xref> and <xref ref-type="fig" rid="fig1">Figure 1</xref>0, the effects of inlet mass flow rates have been studied on flow field over the bundle tubes. <xref ref-type="fig" rid="fig9">Figure 9</xref> shows the Pressure contour for inlet bulk temperature of 300 K and mass flow rates of 0.05, 0.09 and 0.13 kg/s. As can be seen through the figures the value of pressure increases as the inlet mass flow rate grows. At mass flow rate of 0.05 kg/s the static pressure is 0.082 Pa and this value increases to 0.25 and 0.54 for in let mass flow rates of 0.09 and 0.13 kg/s respectively. As also can be seen from the figures the pressure contours at the top and bottom of the tube become smoother with less sharps edges and also more evenly distributed.</p><fig id="fig8"  position="float"><label><xref ref-type="fig" rid="fig8">Figure 8</xref></label><caption><title> Temperature profiles for inlet bulk temperature of 700 K and mass flow rate of 0.05 kg/s at locations 0.01, 0.02 and 0.03 m</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1520084x10.png"/></fig><fig id="fig9"  position="float"><label><xref ref-type="fig" rid="fig9">Figure 9</xref></label><caption><title> Pressure contour for inlet bulk temperature of 300 K and mass flow rate of 0.05 kg/s at locations 0.05, 0.09 and 0.13 m</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1520084x11.png"/></fig><fig id="fig10"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>0</label><caption><title> Temperature contour for inlet bulk temperature of 300 K and mass flow rate of 0.05 kg/s at locations 0.05, 0.09 and 0.13 m</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1520084x12.png"/></fig><p><xref ref-type="fig" rid="fig1">Figure 1</xref>0 depicts the temperature contour for inlet bulk temperature of 300 K and mass flow rates of 0.05, 0.09 and 0.13 kg/s. As can be seen through the figures the value of surface tube temperature does not change but the flow field around the tubes changes and becomes stronger. The temperature contours become sharper behind the tubes and the temperature boundary layer gets thinner as the inlet mass flow rates increases form 0.05, 0.09 to 0.13 kg/s.</p></sec><sec id="s4"><title>4. Conclusion</title><p>Effects of different inlet bulk temperatures and mass flow rate have been investigated on temperature and pressure variation. As numerical results show, the value of pressure increases as the inlet mass flow rate grows and the temperature contours become sharper behind the tubes and the temperature boundary layer gets thinner as the inlet mass flow rates increase. The surface temperature reduces as the inlet bulk temperature passes over the bundle tubes with lower temperature and this value increases as the inlet flow temperature increases to the higher value than the surface temperature.</p></sec><sec id="s5"><title>Acknowledgements</title><p>The author would like to thank both the Saudi Arabian Cultural Mission in Washington DC and King Abdulaziz University in Jeddah, KSA for their support.</p></sec><sec id="s6"><title>Cite this paper</title><p>Almas, M. (2016) Numerical Study of Fluid Flow over Bundle Tubes. Journal of Electronics Cooling and Thermal Control, 6, 109-119. http://dx.doi.org/10.4236/jectc.2016.63010</p></sec></body><back><ref-list><title>References</title><ref id="scirp.70463-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Shan, R.K. and Bhatti, M.S. (1987) Laminar Convective Heat Transfer in Ducts, Hand Book of Single-Phase Convective Heat Transfer. Wiley Interscience, New York</mixed-citation></ref><ref id="scirp.70463-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Colburn, A.P. (1933) A Method of Correlating Forced Convection Heat Transfer Data and A Comparison with Fluid Friction. International Journal of Heat and Mass Transfer, 7, 1359-1384. http://dx.doi.org/10.1016/0017-9310(64)90125-5</mixed-citation></ref><ref id="scirp.70463-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Sieder, E.N. and Tate, G.E. (1936) Heat Transfer and Pressure Drop of Liquids in Tubes. Industrial &amp; Engineering Chemistry Research, 28, 1429-1435. http://dx.doi.org/10.1021/ie50324a027</mixed-citation></ref><ref id="scirp.70463-ref4"><label>4</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Gough</surname><given-names> M.J. </given-names></name>,<etal>et al</etal>. (<year>2011</year>)<article-title>Process Heat Transfer Enhancement to Upgrade Performance, Throughput and Reduced Energy Use</article-title><source> Chemical Engineering Transactions</source><volume> 29</volume>,<fpage> 1</fpage>-<lpage>6</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.70463-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Wang, Y., He, Y.L., Mei, D.H. and Tao, W.Q. (2011) Optimization Design of Slotted Fin by Numerical Simulation Coupled with Genetic Algorithm. Applied Energy, 88, 4441-4450. http://dx.doi.org/10.1016/j.apenergy.2011.05.030</mixed-citation></ref><ref id="scirp.70463-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Yayla, S. (2013) Flow Characteristic of Staggered Multiple Slotted Tubes in the Passage of a Fin Tube Heat Exchanger. Strojniski Vestnik-Journal of Mechanical Engineering, 59, 462-472. http://dx.doi.org/10.5545/sv-jme.2012.902</mixed-citation></ref><ref id="scirp.70463-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Ribatski, G. and Jacobi, A.M. (2005) Falling Film Evaporation on Horizontal Tubes—A Critical Review. International Journal of Refrigeration, 28, 635-653. http://dx.doi.org/10.1016/j.ijrefrig.2004.12.002</mixed-citation></ref></ref-list></back></article>