Acessibilidade / Reportar erro

Integral transform solution of natural convection in a square cavity with volumetric heat generation

Abstract

The generalized integral transform technique (GITT) is employed to obtain a hybrid numerical-analytical solution of natural convection in a cavity with volumetric heat generation. The hybrid nature of this approach allows for the establishment of benchmark results in the solution of non-linear partial differential equation systems, including the coupled set of heat and fluid flow equations that govern the steady natural convection problem under consideration. Through performing the GITT, the resulting transformed ODE system is then numerically solved by making use of the subroutine DBVPFD from the IMSL Library. Therefore, numerical results under user prescribed accuracy are obtained for different values of Rayleigh numbers, and the convergence behavior of the proposed eigenfunction expansions is illustrated. Critical comparisons against solutions produced by ANSYS CFX 12.0 are then conducted, which demonstrate excellent agreement. Several sets of reference results for natural convection with volumetric heat generation in a bi-dimensional square cavity are also provided for future verification of numerical results obtained by other researchers.

Natural convection; Square cavity; Volumetric heat generation; Integral transform; Hybrid solution


PROCESS SYSTEMS ENGINEERING

Integral transform solution of natural convection in a square cavity with volumetric heat generation

C. AnI, II; C. B. VieiraI; J. Su1, * * To whom correspondence should be addressed

INuclear Engineering Program, COPPE, Phone: + (55) (21) 2562 8448, Fax: + (55) (21) 2562 8444, Universidade Federal do Rio de Janeiro, CP 68509, CEP: 21941-972, Rio de Janeiro - RJ, Brazil. E-mail: sujian@nuclear.ufrj.br

IIOffshore Oil/Gas Research Center, College of Mechanical and Transportation Engineering, China University of Petroleum (Beijing), Beijing 102249, China

ABSTRACT

The generalized integral transform technique (GITT) is employed to obtain a hybrid numerical-analytical solution of natural convection in a cavity with volumetric heat generation. The hybrid nature of this approach allows for the establishment of benchmark results in the solution of non-linear partial differential equation systems, including the coupled set of heat and fluid flow equations that govern the steady natural convection problem under consideration. Through performing the GITT, the resulting transformed ODE system is then numerically solved by making use of the subroutine DBVPFD from the IMSL Library. Therefore, numerical results under user prescribed accuracy are obtained for different values of Rayleigh numbers, and the convergence behavior of the proposed eigenfunction expansions is illustrated. Critical comparisons against solutions produced by ANSYS CFX 12.0 are then conducted, which demonstrate excellent agreement. Several sets of reference results for natural convection with volumetric heat generation in a bi-dimensional square cavity are also provided for future verification of numerical results obtained by other researchers.

Keywords: Natural convection; Square cavity; Volumetric heat generation; Integral transform; Hybrid solution.

INTRODUCTION

The investigations of natural convection flow and heat transfer processes, driven by buoyancy forces due to density differences caused by temperature variations in the fluid, have been motivated by their importance in many natural and industrial problems (Cheng, 2006, 2010; Kairi and Murthy, 2011). In all these investigations, considerable work has been done on the study of natural convection in volumetrically heated cavities due to applications in various technological areas such as nuclear reactor design (Baker et al., 1976a,b; Qiu et al., 2010), geophysics (Runcorn, 1962; Mckenzie et al., 1974) and astrophysics (Bethe, 1968; Tritton, 1975).

Many experimental, analytical and numerical approaches have been employed to investigate the characteristics of flow and heat transfer in enclosures with volumetric heat generation. Lee and Goldstein (1988) used a Mach-Zehnder interferometer to investigate experimentally the temperature distribution and the heat transfer rates within an inclined square enclosure containing fluid with internal energy sources bounded by four rigid planes of constant equal temperature. The experimental results of Lee and Goldstein (1988) were later employed to compare with the numerical solutions obtained by May (1991). Based on a finite-difference procedure, Acharya and Goldstein (1985) and Rahman and Sharif (2003) separately studied natural convection in an inclined and vertical square enclosure containing internal energy sources and subjected to external heating, and analyzed the flow patterns for different Rayleigh numbers (both internal and external). To verify the steady-state experimental observations of Kawara et al. (1990), Fusegi et al. (1992) solved numerically natural convection in a differentially heated rectangular cavity with internal heat generation by a control volume-based finite difference technique. Using a finite-volume technique, Di Piazza and Ciofalo (2000) performed direct numerical two-dimensional simulations of low-Prandtl number free convection in a volumetrically heated rectangular enclosure with aspect ratio (AR = 4), adiabatic top/ bottom walls and isothermal side walls. Different flow regimes (steady-state, periodic and chaotic) were obtained based on the values of the Grashof number. Extending previous work, Arcidiacono et al. (2001) and Arcidiacono and Ciofalo (2001) predicted the flow and temperature fields at different Grashof numbers in a square cavity (AR = 1), and subsequently in a shallow cavity (AR = 0.25). By expressing the governing equations in terms of streamfunction, vorticity and temperature formulation, Joshi et al. (2006) conducted a semi-analytical study of natural convection in cavities of different aspect ratios with uniform volumetric heat generation, considering two different boundary conditions, viz., all isothermal walls and only adiabatic horizontal walls. Daniels and Jones (1998) used matched asymptotic expansions method to analyze natural convection in a shallow rectangular cavity due to internal heat gene-ration, where the nonlinear convective effects were studied. Moreover, Qiu et al. (2010) investigated numerically the local Nusselt number, the flow and temperature field at different conditions of two-dimensional natural convection with non-uniform heat generation in a confined enclosure.

Recently, An and Su (2011) employed a hybrid numerical-analytical approach, known as GITT, to study the dynamic response of clamped axially moving beams, where excellent convergence behavior of the integral transform solution was shown by comparing the vibration displacement of different points along the beam length. Indeed, past studies have shown that GITT is a powerful and applicable method for solving diffusion and convection-diffusion problems in heat and fluid flow (Cotta, 1993; Cotta and Mikhailov, 1997; Cotta, 1998). The most interesting feature of this technique is the automatic and straightforward global error control procedure, which makes it particularly suitable for benchmarking pur- poses, and the only mild increase in overall computational effort with increasing number of independent variables. Leal et al. (1999) examined the convergence characteristics of the GITT solution for the steady laminar natural convection of a Newtonian fluid inside rectangular enclosures differentially heated at the vertical walls and demonstrated the usefulness of this integral transform method by performing critical comparisons against previously reported benchmark solutions. As a continuation of their research, Leal et al. (2000) studied the behavior of fluid flow and heat transfer for transient laminar natural convection with variable physical properties inside cavities, and presented the advantages of utilizing the false transient technique for obtaining steady solutions with the GITT method.

The objective of this work is to extend the GITT method to study natural convection in a square cavity with uniform volumetric heat generation, which is one of the key issues concerned in nuclear safety. Since there is no approximation involved in the analytical derivation of the GITT approach, and the problem solution is represented by the sum of a rapidly-converging in finite series, the GITT solutions with the desired precision may serve as benchmarks for testing other numerical methods. Two different boundary conditions are analyzed for the cavity: (a) horizontal walls adiabatic and vertical walls isothermal, (b) all walls isothermal. In the next section, the mathematical formulation of the problem of natural convection in a volumetrically heated cavity is presented. Then, the hybrid numerical-analytical solution is obtained by carrying out integral transform, and numerical results with automatic global accuracy control are presented for different boundary conditions. A comparison against the simulation results of square cavities obtained by ANSYS CFX 12.0 is then performed to assess the accuracy and convergence of the present approach.

MATHEMATICAL FORMULATION

Consider a square cavity of side H containing a fluid with uniform volumetric heat generation rate q'''. Two different boundary conditions are considered in the present study: (a) the horizontal walls are adiabatic and the vertical walls are isothermal, maintained at a constant temperature T0, and (b) all walls are isothermal, maintained at a contant temperature T0, as illustrated in Fig. 1.


The governing equations for two-dimensional steady-state natural convection, under the Boussinesq approximation, are written as:

where x and y are respectively the horizontal and vertical spatial coordinates, u and v the velocity components, p the pressure, T the temperature, g the gravitational acceleration, ρ the density, β the thermal expansion coefficient, cp the specific heat, and k the thermal conductivity of the fluid.

The boundary conditions for the two problems considered are respectively given by:

For two-dimensional incompressible flow, the stream-function ψ* can be defined, which satisfies identically the continuity equation:

The governing Equations (2a), (2b), (3) can be combined into the stream-function only formulation as follows:

where x* and y* identify the dimensional spatial variables, v the kinematic viscosity and α the thermal diffusivity. The dimensionless governing Equations (7), (8) are written as follows:

respectively.

INTEGRAL TRANSFORM SOLUTION

The main procedures of the GITT-type integral transformation process are: (a) to choose auxiliary eigenvalue problems from homogeneous versions of the original problem; (b) to obtain eigenvalues and eigenfunctions which satisfy the orthogonality property; (c) to define the integral transform pair, including the transform formula and the inverse formula; (d) to transform the original problem with transform formula to obtain a set of ordinary differential equations with the calculable coefficients; (e) to solve the ordinary differential equation system by mathematical tools software or codes, such as Mathematica, IMSL Library, etc.; and (f) to obtain the original problem solutions using the inverse formula.

According to the principle of the generalized integral transform technique, the auxiliary eigenvalue problem needs to be chosen for the dimensionless governing Equations (14), (15) with the homogenous boundary conditions (16 or 17). The 'x' direction is selected to be eliminated through integral transformation, and the eigenvalue problem of the biharmonic-type, previously studied in Leal et al. (1999), is adopted for the streamfunction representation as follows:

where Xi(x) and µi are, respectively, the eigenfunctions and eigenvalues of problem (18). The eigenfunctions satisfy the following orthogonality property

with for and for The norm, or normalization integral, is written as

The related eigenfunctions of problem (18) are given by:

where the eigenvalues are obtained from the transcendental equation:

and the normalization integral is evaluated as

Therefore, the normalized eigenfunction coincides, in this case, with the original eigenfunction itself, i.e.

For the temperature expansion, the classical second-order diffusion operator yields a Sturm-Liouville type problem, readily solved with the appropriate boundary conditions of the first kind at the lateral walls, to yield the eigenfunctions and related quantities:

The solution methodology proceeds towards the proposition of the integral transform pair for the potentials, the integral transformation itself and the inversion formula. For the streamfunction field:

The integral transformation process is now employed through operation of (14) with , to find the transformed streamfunction system:

where the related coefficients are obtained analytically through symbolic computation packages, such as Mathematica (Wolfram, 2003), and automatically generated in Fortran form through their definitions:

Similarly, Eq. (15) is operated on with to yield the transformed temperature problem:

where the associated coefficients, symbolically computed and automatically generated within the code, are obtained from

The resulting coupled infinite system of ordinary differential equations with boundary condition at two points for the transformed potentials is then described by Equations (31) and (33), together with the integral transformed boundary conditions:

respectively.

For computational purposes, the expansions are truncated at sufficiently large orders, and , respectively. This truncated system is then solved by a computer program developed in Fortran 90, based on the use of the subroutine DBVPFD from the IMSL Library (IMSL, 2003) with automatic control of the local relative error ( is selected for this problem). For this purpose, the system composed of (31) and (33) is first rewritten as a first order ODE system, i.e.

Once the transformed potentials, and , have been numerically evaluated under controlled accuracy, the inversion formulae (29b) and (30b) are recalled to provide explicit analytical expressions, in the direction, for the original potentials and . The velocity components, and , are analytically expressed from the inversion formula (29b) and the streamfunction definition, in the form:

RESULTS AND DISCUSSION

We now present numerical results for natural convection in a cavity with volumetric heat generation. Two different boundary conditions are considered: i) adiabatic horizontal walls and isothermal vertical walls, ii) all walls isothermal. A Prandtl number of is employed in this study. The solutions of the systems (31, 33) are obtained with and to analyze the convergence behavior, then the accuracy is verified by the numerical results obtained with ANSYS CFX 12.0, where a computational grid is used.

The streamline and isotherm contour plots for an internally heated square enclosure with adiabatic horizontal walls and isothermal vertical walls at are illustrated in Fig. 2. The flow field takes on a pattern possessing bilateral symmetry, where two counter-rotating convective rolls appear. The dimensionless temperature is maximum at the top of the vertical centerline and decreases near the vertical walls, which explains why the flows rise in the central region of the cavity, divide close to the top adiabatic wall, then move downward separately along the vertical walls.


The convergence behavior of the GITT solution is examined by computing the dimensionless temperature and dimensionless vertical velocity with different truncation orders and at the positions , and , as presented in Table 1. For the dimensionless temperature, it can be observed that convergence is essentially achieved with reasonably low truncation orders and .

For a full convergence to five significant digits, more terms (e.g., and are required. As would be expected, convergence rates of the vertical velocity areslower than those of the temperature, which is due to the fact that the velocity is proportional to the first derivative of the streamfunction (Eq. (6b)), while the temperature is defined by the original eigenfunction expansions (Eq. (30b)). This behavior becomes more apparent as increases, but a full convergence to four significant digits still can be obtained for most cases. The results are in excellent agreement with the ones calculated by ANSYS CFX 12.0, which verifies the accuracy of the GITT solution.

The highest dimensionless temperatures, the corresponding positions and the spatially averaged dimensionless temperatures for a cavity with adiabatic horizontal walls at several values of are summarized in Table 2. The GITT yields good agreement with the solution obtained by ANSYS CFX 12.0, with coincidence of four digits for and three digits for . The peak temperature increases and the average temperature decreases monotonically with ; at , the value of is nearly 1.7% lower than that at , which indicates that convection contributes more to heat transport as increases. The variations of and with are plotted in Fig. 3 on a linear-log scale, where the GITT solution curves are constructed by the corresponding values computed for every increment of 200 in Rayleigh numbers.


Figs. 4 and 5 show the variation of the dimensionless temperature and the profile of the dimensionless vertical velocity along the horizontal centerline at , respectively. Fig. 6 reports the profile of the dimensionless horizontal velocity along at . The results are compared with numerical simulations carried out by ANSYS CFX 12.0, for which a close agreement is observed.




The steady-state dimensionless stream function and temperature for a cavity with isothermal walls at are shown in Fig. 7. The streamline exhibits flow behavior similar to that obtained for a cavity with adiabatic horizontal walls; only one circulation cell forms in each half of the enclosure, with a perfect bilateral symmetry. It can be noted that the streamfunction values when all walls are isothermal are approximately half of the corresponding values


for the case with adiabatic horizontal walls, which has been also captured by Joshi et al. (2006). With an approximate up-down symmetry, the isotherm plot indicates that the maximum dimensionless temperature occurs nearly at the center of the cavity and the lowest temperature is located along the side walls.

An analysis similar to the above is performed on the convergence behavior of the dimensionless temperature and dimensionless vertical velocity at the center point of a cavity with isothermal walls, as shown in Table 3. In this case, the Rayleigh numbers are set equal to and respectively. For reference purposes, parallel computations are implemented with ANSYS CFX 12.0. The excellent convergence rates of both and give agreement to four significant digits, with values of and

The values of and for a cavity with isothermal walls are listed in Table 4 to corroborate the accuracy of the GITT solution, irrespective of the magnitude of . Compared to the case with adiabatic horizontal walls, a different behavior is found for the maximum temperature, which declines slightly and the corresponding location moves upward continuously as the value of increases. The value of at is about 0.32% lower than that at

Figures 8 and 9 depict the variation of the dimensionless temperature and the profile of the dimensionless vertical velocity along the horizontal centerline at for a cavity with isothermal walls, respectively. Figure 10 reports the profile of the dimensionless horizontal velocity along at . The plots are qualitatively similar to the ones obtained from a cavity with adiabatic horizontal walls at It is observed that the direction of the horizontal velocity changes its sign at around while the vertical velocity is positive throughout the central region of the cavity and negative near the boundaries.




CONCLUSIONS

The generalized integral transform technique (GITT) has been shown in this work to be an adequate approach for the analysis of natural convection in a square cavity containing internal energy sources, providing an accurate numerical-analytical solution for the temperature and velocity components. The numerical results reported are in excellent agreement with the solution obtained by ANSYS CFX 12.0, which supports the validity of the present computations. This approach can be either employed for benchmarking purposes, yielding sets of reference orders and exceptional computational performances. It is important to note that the methodology developed is applicable to any boundary conditions, and may also be applied to the analysis of natural convection in a rectangular cavity with different aspect ratios or an inclined cavity.

ACKNOWLEDGMENTS

The authors gratefully acknowledge the financial support provided by CNPq, CAPES and FAPERJ of Brazil for their research work. Chen An also would like to acknowledge the financial support provided by the China Scholarship Council.

REFERENCES

Acharya, S., Goldstein, R. J., Natural-convection in an externally heated vertical or inclined square box containing internal energy-sources. Journal of Heat Transfer-Transactions of the ASME, 107(4), 55-865 (1985).

An, C., Su, J., Dynamic response of clamped axially moving beams: Integral transform solution. Applied Mathematics and Computation, 218(2), 249-259 (2011).

Arcidiacono, S., Ciofalo, M., Low-Prandtl number natural convection in volumetrically heated rectangular enclosures III. Shallow cavity, AR = 0.25. International Journal of Heat and Mass Transfer, 44(16), 3053-3065 (2001).

Arcidiacono, S., Di Piazza, I., Ciofalo, M., Low-Prandtl number natural convection in volumetrically heated rectangular enclosures II. Square cavity, AR = 1. International Journal of Heat and Mass Transfer, 44(3), 537-550 (2001).

Baker, L., Faw, R. E., Kulacki, F. A., Postaccident heat removal - Part I: Heat transfer within an internally heated, nonboiling liquid layer. Nuclear Science and Engineering, 61(2), 222-230 (1976a).

Baker, L., Faw, R. E., Kulacki, F. A., Postaccident heat removal - Part II: Heat transfer from an internally heated liquid to a melting solid. Nuclear Science and Engineering, 61(2), 231-238 (1976b).

Bethe, H. A., Energy production in stars. Science, 161(3841), 541-547 (1968).

Cheng, C. Y., Non-Darcy natural convection heat and mass transfer from a vertical wavy surface in saturated porous media. Applied Mathematics and Computation, 182(2), 1488-1500 (2006).

Cheng, C. Y., Natural convection boundary layer flow of fluid with temperature-dependent viscosity from a horizontal elliptical cylinder with constant surface heat flux. Applied Mathematics and Computation 217(1), 83-91 (2010).

Cotta, R. M., Integral Transforms in Computational Heat and Fluid Flow. CRC Press, Boca Raton, FL (1993).

Cotta, R. M., The Integral Transform Method in Thermal and Fluids Science and Engineering. Begell House , New York (1998).

Cotta, R. M., Mikhailov, M. D., Heat Conduction - Lumped Analysis, Integral Transforms, Symbolic Computation. Ed. John Wiley & Sons, Chichester, England (1997).

Daniels, P. G., Jones, O. K., Convection in a shallow rectangular cavity due to internal heat generation. International Journal of Heat and Mass Transfer, 41(23), 3979-3987 (1998).

Di Piazza, I., Ciofalo, M., Low-Prandtl number natural convection in volumetrically heated rectangular enclosures I. Slender cavity, AR = 4. International Journal of Heat and Mass Transfer, 43(17), 3027-3051 (2000).

Fusegi, T., Hyun, J. M., Kuwahara, K., Numerical study of natural-convection in a differentially heated cavity with internal heat-generation - Effects of the aspect ratio. Journal of Heat Transfer-Transactions of the ASME, 114(3), 773-777 (1992).

IMSL, IMSL Fortran Library version 5.0. MATH/ LIBRARY Volume 2. Visual Numerics, Inc., Houston, TX (2003).

Joshi, M. V., Gaitonde, U. N., Mitra, S. K., Analytical study of natural convection in a cavity with volumetric heat generation. Journal of Heat Transfer-Transactions of the ASME, 128(2), 176-182 (2006).

Kairi, R. R., Murthy, P. V. S. N., Effect of viscous dissipation on natural convection heat and mass transfer from vertical cone in a non-newtonian fluid saturated non-darcy porous medium. Applied Mathematics and Computation, 217(20), 8100-8114 (2011).

Kawara, Z., Kishiguchi, I., Aoki, N., Michiyoshi, I., Natural convection in a vertical fluid layer with internal heating. In: Proceedings of the 27th National Heat Transfer Symposium of Japan. Vol. II. pp. 115-117 (1990).

Leal, M. A., Machado, H. A., Cotta, R. M., Integral transform solutions of transient natural convection in enclosures with variable fluid properties. International Journal of Heat and Mass Transfer, 43(21), 3977-3990 (2000).

Leal, M. A., Pérez-Guerrero, J. S., Cotta, R. M., Natural convection inside two-dimensional cavities: The integral transform method. Communications in Numerical Methods in Engineering, 15, 113-125 (1999).

Lee, J. H., Goldstein, R. J., An experimental-study on natural-convection heat-transfer in an inclined square enclosure containing internal energy-sources. Journal of Heat Transfer-Transactions of the ASME, 110(2), 345-349 (1988).

May, H. O., A numerical study on natural-convection in an inclined square enclosure containing internal heat-sources. International Journal of Heat and Mass Transfer, 34(4-5), 919-928 (1991).

Mckenzie, D. P., Roberts, J. M., Weiss, N. O., Convection in earth's mantle - Towards a numerical-simulation. Journal of Fluid Mechanics, 62, 465-538 (1974).

Qiu, S. Z., Qian, L. B., Zhang, D. L., Su, G. H., Tian, W. X., Numerical research on natural convection in molten salt reactor with non-uniformly distributed volumetric heat generation. Nuclear Engineering and Design, 240(4), 796-806 (2010).

Rahman, M., Sharif, M. A. R., Numerical study of laminar natural convection in inclined rectangular enclosures of various aspect ratios. Numerical Heat Transfer Part A-Applications 44 (4), 355-373 (2003).

Runcorn, S. K., Convection currents in earth's mantle. Nature, 195(4848), 1248-1249 (1962).

Tritton, D. J., Internally heated convection in atmosphere of venus and in laboratory. Nature, 257(5522), 110-111 (1975).

Wolfram, S., The Mathematica Book. 5th Ed. Wolfram Media/Cambridge University Press, Champaign, IL, USA (2003).

(Submitted: January 4, 2012 ; Revised: October 9, 2012 ; Accepted: October 10, 2012)

  • Acharya, S., Goldstein, R. J., Natural-convection in an externally heated vertical or inclined square box containing internal energy-sources. Journal of Heat Transfer-Transactions of the ASME, 107(4), 55-865 (1985).
  • An, C., Su, J., Dynamic response of clamped axially moving beams: Integral transform solution. Applied Mathematics and Computation, 218(2), 249-259 (2011).
  • Arcidiacono, S., Ciofalo, M., Low-Prandtl number natural convection in volumetrically heated rectangular enclosures III. Shallow cavity, AR = 0.25. International Journal of Heat and Mass Transfer, 44(16), 3053-3065 (2001).
  • Arcidiacono, S., Di Piazza, I., Ciofalo, M., Low-Prandtl number natural convection in volumetrically heated rectangular enclosures II. Square cavity, AR = 1. International Journal of Heat and Mass Transfer, 44(3), 537-550 (2001).
  • Baker, L., Faw, R. E., Kulacki, F. A., Postaccident heat removal - Part I: Heat transfer within an internally heated, nonboiling liquid layer. Nuclear Science and Engineering, 61(2), 222-230 (1976a).
  • Baker, L., Faw, R. E., Kulacki, F. A., Postaccident heat removal - Part II: Heat transfer from an internally heated liquid to a melting solid. Nuclear Science and Engineering, 61(2), 231-238 (1976b).
  • Bethe, H. A., Energy production in stars. Science, 161(3841), 541-547 (1968).
  • Cheng, C. Y., Non-Darcy natural convection heat and mass transfer from a vertical wavy surface in saturated porous media. Applied Mathematics and Computation, 182(2), 1488-1500 (2006).
  • Cheng, C. Y., Natural convection boundary layer flow of fluid with temperature-dependent viscosity from a horizontal elliptical cylinder with constant surface heat flux. Applied Mathematics and Computation 217(1), 83-91 (2010).
  • Cotta, R. M., Integral Transforms in Computational Heat and Fluid Flow. CRC Press, Boca Raton, FL (1993).
  • Cotta, R. M., The Integral Transform Method in Thermal and Fluids Science and Engineering. Begell House , New York (1998).
  • Cotta, R. M., Mikhailov, M. D., Heat Conduction - Lumped Analysis, Integral Transforms, Symbolic Computation. Ed. John Wiley & Sons, Chichester, England (1997).
  • Daniels, P. G., Jones, O. K., Convection in a shallow rectangular cavity due to internal heat generation. International Journal of Heat and Mass Transfer, 41(23), 3979-3987 (1998).
  • Di Piazza, I., Ciofalo, M., Low-Prandtl number natural convection in volumetrically heated rectangular enclosures I. Slender cavity, AR = 4. International Journal of Heat and Mass Transfer, 43(17), 3027-3051 (2000).
  • Fusegi, T., Hyun, J. M., Kuwahara, K., Numerical study of natural-convection in a differentially heated cavity with internal heat-generation - Effects of the aspect ratio. Journal of Heat Transfer-Transactions of the ASME, 114(3), 773-777 (1992).
  • IMSL, IMSL Fortran Library version 5.0. MATH/ LIBRARY Volume 2. Visual Numerics, Inc., Houston, TX (2003).
  • Joshi, M. V., Gaitonde, U. N., Mitra, S. K., Analytical study of natural convection in a cavity with volumetric heat generation. Journal of Heat Transfer-Transactions of the ASME, 128(2), 176-182 (2006).
  • Kairi, R. R., Murthy, P. V. S. N., Effect of viscous dissipation on natural convection heat and mass transfer from vertical cone in a non-newtonian fluid saturated non-darcy porous medium. Applied Mathematics and Computation, 217(20), 8100-8114 (2011).
  • Kawara, Z., Kishiguchi, I., Aoki, N., Michiyoshi, I., Natural convection in a vertical fluid layer with internal heating. In: Proceedings of the 27th National Heat Transfer Symposium of Japan. Vol. II. pp. 115-117 (1990).
  • Leal, M. A., Machado, H. A., Cotta, R. M., Integral transform solutions of transient natural convection in enclosures with variable fluid properties. International Journal of Heat and Mass Transfer, 43(21), 3977-3990 (2000).
  • Leal, M. A., Pérez-Guerrero, J. S., Cotta, R. M., Natural convection inside two-dimensional cavities: The integral transform method. Communications in Numerical Methods in Engineering, 15, 113-125 (1999).
  • Lee, J. H., Goldstein, R. J., An experimental-study on natural-convection heat-transfer in an inclined square enclosure containing internal energy-sources. Journal of Heat Transfer-Transactions of the ASME, 110(2), 345-349 (1988).
  • May, H. O., A numerical study on natural-convection in an inclined square enclosure containing internal heat-sources. International Journal of Heat and Mass Transfer, 34(4-5), 919-928 (1991).
  • Mckenzie, D. P., Roberts, J. M., Weiss, N. O., Convection in earth's mantle - Towards a numerical-simulation. Journal of Fluid Mechanics, 62, 465-538 (1974).
  • Qiu, S. Z., Qian, L. B., Zhang, D. L., Su, G. H., Tian, W. X., Numerical research on natural convection in molten salt reactor with non-uniformly distributed volumetric heat generation. Nuclear Engineering and Design, 240(4), 796-806 (2010).
  • Rahman, M., Sharif, M. A. R., Numerical study of laminar natural convection in inclined rectangular enclosures of various aspect ratios. Numerical Heat Transfer Part A-Applications 44 (4), 355-373 (2003).
  • Runcorn, S. K., Convection currents in earth's mantle. Nature, 195(4848), 1248-1249 (1962).
  • Wolfram, S., The Mathematica Book. 5th Ed. Wolfram Media/Cambridge University Press, Champaign, IL, USA (2003).
  • *
    To whom correspondence should be addressed
  • Publication Dates

    • Publication in this collection
      10 Jan 2014
    • Date of issue
      Dec 2013

    History

    • Received
      04 Jan 2012
    • Accepted
      10 Oct 2012
    • Reviewed
      09 Oct 2012
    Brazilian Society of Chemical Engineering Rua Líbero Badaró, 152 , 11. and., 01008-903 São Paulo SP Brazil, Tel.: +55 11 3107-8747, Fax.: +55 11 3104-4649, Fax: +55 11 3104-4649 - São Paulo - SP - Brazil
    E-mail: rgiudici@usp.br