<?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">AJCM</journal-id><journal-title-group><journal-title>American Journal of Computational Mathematics</journal-title></journal-title-group><issn pub-type="epub">2161-1203</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ajcm.2016.63026</article-id><article-id pub-id-type="publisher-id">AJCM-70667</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Numerical Investigation of a Mixed Convection Flow in a Lid-Driven Cavity
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Reyad</surname><given-names>Omari</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 Mathematics, Al-Balqa Applied University, Irbid University College, Irbid, Jordan</addr-line></aff><author-notes><corresp id="cor1">* E-mail:</corresp></author-notes><pub-date pub-type="epub"><day>04</day><month>07</month><year>2016</year></pub-date><volume>06</volume><issue>03</issue><fpage>251</fpage><lpage>258</lpage><history><date date-type="received"><day>12</day>	<month>May</month>	<year>2016</year></date><date date-type="rev-recd"><day>accepted</day>	<month>16</month>	<year>September</year>	</date><date date-type="accepted"><day>19</day>	<month>September</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>
 
 
  This study is devoted to the computational fluid dynamics (CFD) modeling of steady laminar mixed convection flow and heat transfer in lid driven cavity (10 ≤ &lt;i&gt;Re&lt;/i&gt;
   
  ≤ 1000). The ratio of the height to the width of the cavity is ranged over &lt;i&gt;H/L&lt;/i&gt; = 0.5 to 1.5. The governing equations are solved using commercial finite volume package FLUENT to visualize the nature of the flow and estimate the heat transfer inside the cavity for different aspect ratio. The simulation results are presented in terms of average Nusselt number of the hot wall, velocity profile, and temperature contours. It was found that the average Nusselt number inside the cavity is strongly governed by the aspect ratio as well as the Reynolds number. A parametric study is conducted to demonstrate the effect of aspect ratio on the flow and heat transfer characteristics. It is found that heat transfer enhancement was obtained by decreasing the aspect ratio and/or increasing the Reynolds number.
 
</p></abstract><kwd-group><kwd>CFD Simulation</kwd><kwd> Mixed Convection</kwd><kwd> Lid Driven Cavity</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>One of the major problems in the thermos-fluids is the free and forced convection flow and heat transfer in different aspect of cavity. The lid-driven cavity flow is the motion of a fluid inside a rectangular cavity created by a constant translational velocity of one side while the other sides remain at rest. Fluid and heat flow behaviors inside lid driven cavities have been the subject of extensive computational and experimental studies over the past years. It has many applications in engineering, material processing, dynamics of lakes, metal casting and galvanizing.</p><p>The influence of a lid movement on the fluid flow and thermal structures in a lid-driven cavity was conducted by [<xref ref-type="bibr" rid="scirp.70667-ref1">1</xref>] . They concluded that the maximum local heat transfer rate occurs at the starting area of the sliding lid and decreases along the sliding lid. Mixed convection and conduction heat transfer in open cavities investigated numerically by [<xref ref-type="bibr" rid="scirp.70667-ref2">2</xref>] . They found that heat transfer across the cavity is enhanced when the cavity aspect ratio is increased. A combined forced and natural convection in a lid-driven cavity is investigated experimentally by [<xref ref-type="bibr" rid="scirp.70667-ref3">3</xref>] . The case of transient laminar mixed convention flow and heat transfer of an electrically conducting and heat generating in a lid driven cavity was investigated by [<xref ref-type="bibr" rid="scirp.70667-ref4">4</xref>] . He concluded that the existence of the internal heat generation decreases the average Nusselt number significantly. The unsteady free and forced convection in partially divided cavity were investigated by [<xref ref-type="bibr" rid="scirp.70667-ref5">5</xref>] . They found that the time to reach the steady state depends on the structure of divider.</p><p>The cubic interpolated pseudo particle method and verified their results with the shear driven flow in square cavities were considered by [<xref ref-type="bibr" rid="scirp.70667-ref6">6</xref>] . Li et al. [<xref ref-type="bibr" rid="scirp.70667-ref7">7</xref>] applied the new version of multiple relaxation time lattices Boltzmann method to investigate the fluid flow in deep cavity. Manca et al. [<xref ref-type="bibr" rid="scirp.70667-ref8">8</xref>] showed a numerical analysis of laminar mixed convection in an open cavity with a heated wall bounded by a horizontally insulated plate. Their results were reported for Reynolds numbers in the range of 100 to 1000 and aspect ratio from 0.1 to 1.5. They presented that the maximum decrease in temperature was occurred at higher Reynolds number. They found that the effect of the aspect ratio was played a significant role on streamlines and isotherm patterns for different heating configurations. Their investigations also indicate that opposing forced flow configuration has the highest thermal performance, in terms of both maximum temperature and average Nusselt number.</p><p>Gau and Sharif [<xref ref-type="bibr" rid="scirp.70667-ref9">9</xref>] investigated the mixed convection in rectangular cavities at various aspect ratios with moving isothermal side walls and constant wall flux heat source on the bottom wall. Numerical simulation of unsteady mixed convection in a driven cavity using an externally excited sliding lid is conducted by Khanafer et al. [<xref ref-type="bibr" rid="scirp.70667-ref10">10</xref>] . They observed that, Re and Gr would either enhance or retard the energy transport process and drag force behavior depending on the conduct of the velocity cycle. Recently, nanofluids, which are considered very promising in heat transfer enhancement, were used in several studies. Khanafer et al. [<xref ref-type="bibr" rid="scirp.70667-ref11">11</xref>] investigated the problem of buoyancy-driven heat transfer enhancement of nanofluids in a two-dimensional enclosure. Jang and Choi [<xref ref-type="bibr" rid="scirp.70667-ref12">12</xref>] investigated free convection in a rectangular cavity (Benard convection) with nanofluids. Jou and Tzeng [<xref ref-type="bibr" rid="scirp.70667-ref13">13</xref>] present a numerical study of the heat transfer performance of nanofluids inside two dimensional rectangular cavity.</p><p>Wang et al. [<xref ref-type="bibr" rid="scirp.70667-ref14">14</xref>] studied free convection heat transfer in horizontal and vertical rectangular cavities filled with nanofluids. Santra et al. [<xref ref-type="bibr" rid="scirp.70667-ref15">15</xref>] investigated heat transfer characteristics of copper water nanofluid in a differentially heated square cavity with different viscosity models.</p><p>Many researchers, investigate the mixed convection heat transfer characteristics of a square cavity over the years, however, there seems to be lack studies relates the numerical investigation covered the effect of aspect ratio. Therefore, the main objective of this study is to study the effect of aspect ratio (0.5 ≤ AR ≤ 1.5) on fluid and heat transfer characteristic for laminar, steady and mixed convection using commercial finite volume package FLUENT. Also, the present study focuses on flow visualizations on moderate Reynolds number range varying the aspect ratio of driven cavity.</p></sec><sec id="s2"><title>2. Theoretical Formulation</title><p>The fluid flow and heat transfer in open lid driven cavity can be simulated by a set of mass, momentum and energy conservation equations. The flow is assumed to be two-dimensional, laminar, incompressible, Newtonian and neglecting the viscous effect. The dimensionless governing non-linear partial differential equations can be written as follows:</p><disp-formula id="scirp.70667-formula86"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/7-1100525x6.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.70667-formula87"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/7-1100525x7.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.70667-formula88"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/7-1100525x8.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.70667-formula89"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/7-1100525x9.png"  xlink:type="simple"/></disp-formula><p>The dimensional variables are as follows:</p><disp-formula id="scirp.70667-formula90"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/7-1100525x10.png"  xlink:type="simple"/></disp-formula><p>The L is the reference length dimension (width of the cavity along lower lid), while U<sub>∞</sub> is the reference velocity dimension. The fluid property ν, refers to the kinematic viscosity. The Reynold’s number, Re, is the ratio of inertial to viscous forces, which influences the fluid flow features within the cavity. No-slip velocity boundary condition (u = v = 0) is applied on all the walls, except the top lid. On the top lid (U = 1 and V = 0) is applied. The bottom boundary of the domain is modeled as wall .The boundary conditions which describing the current simulated computational domain as well as the surface boundary layer is depicted in <xref ref-type="fig" rid="fig1">Figure 1</xref>.</p><disp-formula id="scirp.70667-formula91"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/7-1100525x11.png"  xlink:type="simple"/></disp-formula><p>A finite volume method is employed using commercial software FLUENT 6.2 to solve the governing equations together with the specified boundary conditions. The velocities and pressure are coupled by the SIMPLE algorithm. More cells are constructed near the surface of the cavity to compensate the high velocity gradient in the boundary layer region of the viscous flow. A commercial software GAMBIT is used for grid generation. A second order upwind scheme is used for the convection. Here in this study, following [<xref ref-type="bibr" rid="scirp.70667-ref10">10</xref>] , we define the Nusselt number, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-1100525x12.png" xlink:type="simple"/></inline-formula>as follows:</p><disp-formula id="scirp.70667-formula92"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/7-1100525x13.png"  xlink:type="simple"/></disp-formula><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-1100525x14.png" xlink:type="simple"/></inline-formula>is the fluid thermal conductivity, Y is appropriate cavity length and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/7-1100525x15.png" xlink:type="simple"/></inline-formula> is non dimnial temperature. The grid independence is achieved by comparing the results of the different grid cell size. It was found that 65,000 cells are satisfactory, and any increase beyond this size would lead to an insignificant change in the resulting solution.</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Schematic of the square lid-driven cavity</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/7-1100525x16.png"/></fig></sec><sec id="s3"><title>3. Results and Discussion</title><p>Simulation results for laminar mixed convection in square cavity (AR = 1) are compared to experimental data to verify the validity of the CFD simulation solution. Then, the code is employed to investigate the mixed convection fluid flow and heat transfer inside the different aspect ratio of cavity. <xref ref-type="fig" rid="fig2">Figure 2</xref> and <xref ref-type="fig" rid="fig3">Figure 3</xref> show the U and V-velocity profiles at the vertical and horizontal mid-plane of the square cavity for (AR = 1). As can be seen there is an excellent agreement in the velocity profile between CFD simulations in this study and the experimental measured by [<xref ref-type="bibr" rid="scirp.70667-ref16">16</xref>] . It has been noticed that the heating of the surface of the cavity has no effect on the velocity profile.</p><p><xref ref-type="fig" rid="fig4">Figure 4</xref> shows the influence of the aspect ratio on the local Nusselt number for Re = 100. It is worth indicating that all the simulation results of local Nusselt number and temperature profile presented in this work correspond to the cold lid. It is noticeable, that the local Nusselt number decreases as x increase and increases as the</p><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> Computed of U-velocity profiles along a vertical line passing through the geometric centre of the cavity (at X = 0.5) at various aspect ratios, Re = 1000</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/7-1100525x17.png"/></fig><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> Computed of V-velocity profiles along a horizontal line passing through the geometric centre of the cavity (at Y = 0.5) at various aspect ratios, Re = 1000</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/7-1100525x18.png"/></fig><p>aspect ratio decreases for Re = 100. Decreasing the aspect ratio means that the cold lid would be close to the hot wall and sense the heat significantly. Therefore, the heat transfer enhanced near the cold lid.</p><p><xref ref-type="fig" rid="fig5">Figure 5</xref> shows the influence of the aspect ratio on the profile temperature for Re = 100 to the right wall. It is also noticeable, the profile temperature decreases as y increase and increases as the aspect ratio increases for Re = 100. For the same reason, increasing the aspect ratio means that the right wall would be close to the hot wall and sense the heat significantly. Therefore, the heat transfer enhanced near the right wall.</p><p><xref ref-type="fig" rid="fig6">Figure 6</xref> shows the influence of the aspect ratio on the profile temperature for Re = 100 to the left wall. It is also noticeable, the profile temperature decreases as y increase and increases as the aspect ratio increases for Re = 100. For the same reason, increasing the aspect ratio means that the right wall would be close to the hot wall and sense the heat significantly. Therefore, the heat transfer enhanced near the right wall.</p><p>The contours of static temperature that has the same temperature difference and various aspect ratio are shown in <xref ref-type="fig" rid="fig7">Figure 7</xref>. The mixed convection cases show different pattern. In order to visualize the overall thermal</p><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> Distribution of local Nusselt number at hot lid for various aspect ratios, Re = 100</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/7-1100525x19.png"/></fig><fig id="fig5"  position="float"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Typical profile of temperature θ distribution at hot lid for various aspect ratios, Re = 100</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/7-1100525x20.png"/></fig><p>patterns, the isotherm contours plots from the computational data for different aspect ratio are at Reynolds number = 100 is shown in <xref ref-type="fig" rid="fig8">Figure 8</xref>. For smaller aspect ratio the centre of the primary eddy moves right with</p><fig id="fig6"  position="float"><label><xref ref-type="fig" rid="fig6">Figure 6</xref></label><caption><title> Typical profile of temperature θ distribution at left wall for various aspect ratios, Re = 100</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/7-1100525x21.png"/></fig><fig id="fig7"  position="float"><label><xref ref-type="fig" rid="fig7">Figure 7</xref></label><caption><title> Typical profile of temperature θ distribution at mid of cavity for various aspect ratios, Re = 100</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/7-1100525x22.png"/></fig><fig id="fig8"  position="float"><label><xref ref-type="fig" rid="fig8">Figure 8</xref></label><caption><title> Contours of static temperature in (K) inside cavity for various aspect ratios, Re = 100</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/7-1100525x23.png"/></fig><p>respect to right wall cavity. Further, an energy rate increases due to circulation movements of fluid molecules which enhance the thermal dispersion especially for smaller cavity aspect ratio.</p></sec><sec id="s4"><title>Cite this paper</title><p>Reyad Omari, (2016) Numerical Investigation of a Mixed Convection Flow in a Lid-Driven Cavity. American Journal of Computational Mathematics,06,251-258. doi: 10.4236/ajcm.2016.63026</p></sec></body><back><ref-list><title>References</title><ref id="scirp.70667-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Bruneau, Ch.-H. and Saad, M. (2006) The 2D Lid-Driven Cavity Problem Revisited. Computers &amp; Fluids, 35, 326-348. http://dx.doi.org/10.1016/j.compfluid.2004.12.004</mixed-citation></ref><ref id="scirp.70667-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Batchelor, G.K. (1956) On Steady Laminar Flow with Closed Streamlines at Large Reynolds Numbers. Journal of Fluid Mechanics, 1, 177-190. http://dx.doi.org/10.1017/S0022112056000123</mixed-citation></ref><ref id="scirp.70667-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Wang, X.-Q., Mujumdar, A.S. and Yap, C. (2006) Free Convection Heat Transfer in Horizontal and Vertical Rectangular Cavities Filled with Nanofluids. In: International Heat Transfer Conference IHTC-13, Sydney, 1-12. http://dx.doi.org/10.1615/ihtc13.p7.410</mixed-citation></ref><ref id="scirp.70667-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Jou, R.Y. and Tzeng, S.C. (2006) Numerical Research of Nature Convective Heat Transfer Enhancement Filled with Nanofluids in Rectangular Enclosures. International Journal of Heat and Mass Transfer, 33, 727-736. http://dx.doi.org/10.1016/j.icheatmasstransfer.2006.02.016</mixed-citation></ref><ref id="scirp.70667-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Jang, S.P. and Choi, S.U.S. (2004) Free Convection in a Rectangular Cavity Benard Convection with Nanofluids. in: Proceedings of the IMECE, Anaheim, California, USA,. http://dx.doi.org/10.1115/imece2004-61054</mixed-citation></ref><ref id="scirp.70667-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Khanafer, K., Vafai, K. and Lightstone, M. (2003) Buoyancy Driven Heat Transfer Enhancement in a Two-Dimensional Enclosure Utilizing Nanofluids. International Journal of Heat and Mass Transfer, 46, 3639-3653. http://dx.doi.org/10.1016/S0017-9310(03)00156-X</mixed-citation></ref><ref id="scirp.70667-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Khanafer, K.M., Al-Amiri, A.M. and Pop, I. (2007) Numerical Simulation of Unsteady Mixed Convection in a Driven Cavity Using an Externally Excited Sliding Lid. European Journal of Mechanics—B/Fluids, 26, 669-687. http://dx.doi.org/10.1016/j.euromechflu.2006.06.006</mixed-citation></ref><ref id="scirp.70667-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Gau, G. and Sharif, M.A.R. (2004) Mixed Convection in Rectangular Cavities at Various Aspect Ratios with Moving Isothermal Side Walls and Constant Flux Heat Source on the Bottom Wall. International Journal of Thermal Sciences, 43, 465-475.</mixed-citation></ref><ref id="scirp.70667-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Manca, O., Nardini, S., Khanafer, K. and Vafai, K. (2003) Effect of Heated Wall Position on Mixed Convection in a Channel with an Open Cavity. Numerical Heat Transfer, Part A: Applications, 43, 259-282.</mixed-citation></ref><ref id="scirp.70667-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Li, S.L., Chen, Y.C. and Lin, C.A. (2011) Multi Relaxation Time Lattice Boltzmann Simulations of Deep Lid Driven Cavity Flows at Different Aspect Ratios. Computers &amp; Fluids, 45, 233-240. http://dx.doi.org/10.1016/j.compfluid.2010.12.012</mixed-citation></ref><ref id="scirp.70667-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Sidik, N.A.C. and Attarzadeh, S.M.R. (2011) An Accurate Numerical Prediction of Solid Particle Fluid Flow in a Lid-Driven Cavity. International Journal of Mechanics, 5, 123-128.</mixed-citation></ref><ref id="scirp.70667-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">How, S.P. and Hsu, T.H. (1998) Transient Mixed in a Partially Divided Enclosure. Computers &amp; Mathematics with Applications, 36, 95-115. http://dx.doi.org/10.1016/S0898-1221(98)00186-2</mixed-citation></ref><ref id="scirp.70667-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Chamkha, A.J. (2002) Hydromagnetic Combined Convection Flow in a Vertical Lid-Driven Cavity with Internal Heat Generation or Absorption. Numerical Heat Transfer, Part A: Applications, 41, 529-546. http://dx.doi.org/10.1080/104077802753570356</mixed-citation></ref><ref id="scirp.70667-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Prasad, A.K. and Koseff, J.R. (1996) Combined Forced and Natural Convection Heat Transfer in a Deep Lid Driven Cavity Flow. International Journal of Heat and Fluid Flow, 17, 460-467. http://dx.doi.org/10.1016/0142-727X(96)00054-9</mixed-citation></ref><ref id="scirp.70667-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Chenak, A.K., Mbaye, M., Vasseur, P. and Bilgen, E. (1995) Mixed Convection and Conduction Heat Transfer in Open Cavities. Heat and Mass Transfer, 30, 229-235. http://dx.doi.org/10.1007/BF01602767</mixed-citation></ref><ref id="scirp.70667-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Mohamad, A.A. and Viskanta, R. (1995) Flow and Heat Transfer in a Lid-Driven Cavity Filled with a Stably Stratified Fluid. Applied Mathematical Modelling, 19, 465-472. http://dx.doi.org/10.1016/0307-904X(95)00030-N</mixed-citation></ref></ref-list></back></article>