Laminar natural convection of non-Newtonian nanofluids in a square enclosure with differentially heated side walls.
Today more than ever, ultrahigh-performance and controlled heat transfer plays an important role in the development of energy-efficient heat transfer systems and fluids which are required in many industries and commercial applications (1-3). However, conventional heat transfer fluids (e.g. water, oil or ethylene glycol) are inherently poor heat transfer fluids. For that, nano-sized particles dispersed in a base fluid, known as nanofluid (4), have been developed, used and researched extensively to enhance heat transfer. The presence of nanoparticles shows an unquestionable heat transfer enhancement in forced convection applications. But on the other hand, there is still a dispute on the role of nanoparticles on the heat transfer for the buoyancy driven flow.
Natural convection (i.e. flow caused by the temperature induced density variations) is one of the most extensively analysed configurations because of its fundamental importance as the "benchmark" problem to study convection effects (and compare as well as validate numerical techniques). In addition to the obvious academic interest, this configuration has relevance in various engineering systems (e.g. compact heat exchangers (5) and cooling systems (6), (7), selective laser melting process (8)). Although quite some various different configurations of the enclosure problem are possible, one of the most studied cases involves the two-dimensional square enclosure with differentially heated isothermal vertical walls and adiabatic horizontal walls. When the vertical walls are insulated to ensure adiabatic conditions and the lower horizontal wall held at the higher temperature then one has the classical Rayleigh--Benard configuration.
The present study analyses the case where the horizontal walls are adiabatic and the temperature difference driving the convection comes from the sidewalls as in the classic benchmark paper of de Vahl Davis (9). Although this configuration has been studied extensively for Newtonian fluids, a much smaller amount of information is available for the nanofluids. Most of the researchers have been focused on numerical modelling of buoyancy induced flows of Newtonian (mainly water-based) nanofluids. Recent numerical studies by Ternik et al. (10), Ternik and Rudolf (11), Oztop et al. (12) and Abu-Nada and Oztop (13) illustrated that nanoparticles suspended in Newtonian base fluid substantially enhance the heat transfer for any given Rayleigh number.
On the other hand, even less research work have been conducted if the rheological behaviour of nanofluid is more complex (i.e. the relationship between the shear stress and shear rate is non-linear). Natural convection heat transfer of non-Newtonian nanofluids (dispersion of [Al.sub.2][O.sub.3] and Ti[O.sub.2] nanoparticles in a 0,5 wt.% aqueous solution of carboxymethyl cellulose) in a vertical cylinder uniformly heated from below and cooled from top was investigated experimentally by Mahrood et al. (14). They showed that the heat transfer performance of nanofluids is significantly enhanced at low particle concentrations and that increasing nanoparticle concentration has a contrary effect on the heat transfer of nanofluids. Cheng (15) studied the free convection heat transfer over a truncated cone embedded in a porous medium saturated by a non-Newtonian power-law nanofluid with constant wall temperature and constant wall nanoparticle volume fraction and reported that increasing the thermophoresis parameter or the Brownian parameter tends to decrease the reduced Nusselt number. Forced convective heat transfer effect on the non-Newtonian nanofluid (containing [Al.sub.2][O.sub.3] and Xanthan aqueous solution as a liquid single phase) flow in the horizontal tube with constant heat flux was numerically investigated by Moraveji et al. (16). Their results showed that heat transfer coefficient and Nusselt number of non-Newtonian nanofluid increased with increasing concentration of Xanthan solution.
The above review of the existing literature shows that the problem of natural convection in a square enclosure filled with a non-Newtonian nanofluid is an issue still far from being completely solved. Framed in this general background, the purpose of the present study is to examine the effect of adding [Al.sub.2][O.sub.3], Au, Cu and Ti[O.sub.2] nanoparticles to the base fluid (which exhibits non-Newtonian viscous behaviour) on the heat and momentum transfer characteristics in a square cavity with differentially heated vertical and adiabatic horizontal walls over the range of the nanofluid Rayleigh number and volume fraction.
The rest of the paper is organised as follows. The necessary mathematical background and numerical details are presented in Section 2 which is followed by the grid refinement, numerical accuracy assessment and validation study (Section 3). Following this analysis, the results are presented and subsequently discussed (Section 4). The main findings are summarised and conclusions are drawn in the final section of this paper.
2. NUMERICAL MODELLING
The standard finite volume method (successfully used in numerical studies of isothermal (17) and non-isothermal (10), (11) generalized Newtonian fluid flows) is used to solve the coupled conservation equations of mass, momentum and energy. In the present framework, the second-order central differencing scheme is used for the diffusive terms and the second-order Higher upwind scheme for the convective terms. The convergence criteria were set to [10.sup.-9] for all the relative (scaled) residuals.
2.1 Governing equations
For the present study steady-state flow of an incompressible CMC based nanofluids is considered. It is assumed that the fluid phase and nanoparticles are in both, thermal and chemical equilibrium. Except for the density the properties of nanoparticles and fluid are taken to be constant. Table I presents the thermo-physical properties of CMC and Au nanoparticles at the reference temperature. It is further assumed that the Boussinesq approximation is valid for the buoyancy force.
Table I: Thermo-physical properties of nanofluids. [rho][kg / [c.sub.[rho]] k [W / [beta][1 / K] [m.sup.3]] [J / kgK] mK] 0,4 % CMC 997,1 4179 0,613 2,1 x [10.sup.-4] Au 19320 128,8 314,4 1,416 x [10.sup.-7] [Al.sub.2] 3970 765 40 8,5 x [10.sup.-6] [O.sub.3] Cu 8933 385 400 1,67 x [10.sup.-5] Ti[O.sub.2] 4250 686,2 8,9538 9,0 x [10.sup.-6]
The governing equations (mass, momentum and energy conservation) for a steady, laminar and incompressible flow are:
[[partial derivative][v.sub.i]]/[[partial derivative][x.sub.i]] = 0 (1)
[[rho].sub.nf][v.sub.j][[[partial derivative][v.sub.i]]/[[partial derivative][x.sub.j]]] - [[partial derivative]/[[partial derivative][x.sub.j]]][[[eta].sub.nf]([absolute value of [gamma]])[[[partial derivative][v.sub.i]]/[[partial derivative][x.sub.j]]]] = - [[partial derivative]p]/[[rho][x.sub.j]] + [([rho][beta]).sub.nf]g(T - [T.sub.c]) + [[partial derivative]/[[partial derivative][x.sub.j]]][[[eta].sub.nf]([absolute value of [gamma]])[[[partial derivative][v.sub.j]]/[[partial derivative][x.sub.i]]]] (2)
[([rho][c.sub.p]).sub.nf][v.sub.j][[partial derivative]T]/[[partial derivative][x.sub.j]] = [[partial derivative]/[[partial derivative][x.sub.j]]]([k.sub.nf][[partial derivative]T]/[[partial derivative][x.sub.j]]) (3)
In the momentum balance law (Eq. (2)) the constitutive (rheological) model is needed for the viscous function [[eta].sub.nf]([absolute value of [gamma]]). In the case of generalized Newtonian fluids, the form of the general relationship between the shear stress and shear rate is:
[[tau].sub.ij] = [[eta].sub.nf]([absolute value of [gamma]])[[gamma].sub.ij] (4)
where [absolute value of [gamma]] = [square root of ((1/2)[[SIGMA].sub.i][[SIGMA].sub.j][[gamma].sub.ij][[gamma].sub.ji])] is the II. invariant of the symmetrical rate of deformation tensor with Cartesian components [[gamma].sub.ij] = ([[partial derivative][v.sub.i]]/[[partial derivative][x.sub.j]]) + ([[partial derivative][v.sub.j]]/[[partial derivative][x.sub.i]]). Based on Eq. (4) the viscous function for the Power law (i.e. Ostwald-de Waele) rheological model is:
[[eta].sub.nf]([absolute value of [gamma]]) = K[[absolute value of [gamma]].sup.n -1] (5)
where K is the consistency and n is the dimensionless Power law index with experimentally determined values (for 0,4 % aqueous CMC solution) K = 0,048 [Pas.sup.n] and n = 0,882.
Relationships between properties of nanofluid (nf) to those of base fluid (bf) and pure solid (s) are given with the following empirical models (10), (11), (13):
* Density: [[rho].sub.nf] = (1 - [phi])[[rho].sub.bf] + [phi][[rho].sub.s]
* Dynamic viscosity: [[eta].sub.nf] = [[eta].sub.bf]/[(1 - [phi]).sup.2.5]
* Thermal expansion: [([rho][beta]).sub.nf] = (1 - [phi])[([rho][beta]).sub.bf] + [phi][([rho][beta]).sub.s]
* Heat capacitance: [([rho][c.sub.p]).sub.nf] = (1 - [phi])[([rho][c.sub.p]).sub.bf] + [phi][([rho][c.sub.p]).sub.s]
* Thermal conductivity: [k.sub.nf] = [k.sub.bf][[[k.sub.s] + 2[k.sub.bf] - 2[phi]([k.sub.bf] - [k.sub.s])]/[[k.sub.s] + 2[k.sub.bf] + [phi]([k.sub.bf] - [k.sub.s])]]
2.2 Geometry and boundary conditions
The simulation domain and the expected temperature distribution are shown schematically in Fig. 1. Two vertical walls of a square enclosure are kept at different constant temperatures ([T.sub.H] > [T.sub.C]), whereas the other boundaries are considered to be adiabatic in nature. Both velocity components (i.e. [v.sub.x] and [v.sub.y]) are identically zero on each boundary because of the no-slip condition and impenetrability of rigid boundaries. The temperatures for cold and hot vertical walls are specified ([T.sub.x = 0] = [T.sub.H], [T.sub.x = L] = [T.sub.C]). The adiabatic temperature boundary conditions for the horizontal insulated boundaries are given by [partial derivative]T/[partial derivative]y = 0 at y = 0 and y = L.
In the present study, the heat transfer characteristics are presented (and compared for the same value of the nanofluid Rayleigh number) in terms of the mean Nusselt number:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (6)
and relative enhancement of the heat transfer rate between the nanofluid (nf) and the base fluid (bf):
E = [[[Q.sub.nf] - [Q.sub.bf]]/[Q.sub.bf]] x 100% = [[[k.sub.nf]/[k.sub.bf]][[[bar.Nu].sub.nf]/[[bar.Nu].sub.bf]] - 1] x 100% (7)
where [[Q.sub.nf] = [[k.sub.nf][[bar.Nu].sub.nf]A[DELTA]T]/L and [[Q.sub.bf] = [[k.sub.bf][[bar.Nu].sub.bf]A[DELTA]T]/L are the heat transfer rates of nanofluid and base fluid. Here, the nanofluid Rayleigh number [Ra.sub.nf] represents the ratio of the strengths of buoyant thermal transports to thermal diffusion and for the non-Newtonian Power-law nanofluids is expressed as (18), (19):
[Ra.sub.nf] = [[[rho].sub.nf][([rho][c.sub.p]).sub.nf.sup.n]g[[beta].sub.nf][DELTA]T[L.sup.2n + 1]]/[k.sub.nf.sup.n]K (8)
In order to investigate the influence of solid particles volume fraction [phi] on the heat transfer characteristics (i.e. mean Nusselt number and heat transfer rate enhancement), the Prandtl number of CMC based Power-law nanofluids is defined (18), (19):
[Pr.sub.nf] = [K/[[rho].sub.nf]][(k/[rho][c.sub.p]).sub.nf.sup.n - 2][L.sup.2 - 2n] (9)
Using Eq. (9) we show (Fig. 2) that value of the nanofluid Prandtl number decreases as the volume fraction of nanoparticles in CMC base fluid increases and that for the investigated range of [phi] its values are [Pr.sub.nf] [much greater than] 1.
3. GRID REFINEMENT, NUMERCAL ACCURACY AND VALIDATION
The grid independence and numerical accuracy of the present results has been established on a detailed analysis using three different uniform meshes; MI(50 x 50), MII(100 x 100) and MIII(200 x 200). With each grid refinement the number of elements in particular direction is doubled and element size is halved. Such a procedure is useful for applying the Richardson's extrapolation technique which is a method for obtaining a higher-order estimate of the flow value (value at infinite grid) from a series of lower-order discrete values (20).
For a general primitive variable the grid-converged (i.e. extrapolated to the zero element size) value according to Richardson extrapolation is given as:
[[PHI].sub.ext] = [[PHI].sub.MIII] - ([[PHI].sub.MII] - [[PHI].sub.MIII])/([r.sup.p] - 1) (10)
where [[PHI].sub.MIII] is obtained on the finest grid and MII is the solution based on next level of coarse grid, r = 2 is the ratio between the coarse to fine grid spacing and p = 2 is the order of convergence.
The variation of the mean Nusselt number [bar.Nu] with the grid refinement is provided in tabulated form in Table II. The "percent" numerical error for the mean Nusselt number:
Error = [absolute value of ([[bar.Nu].sub.MIII] - [[bar.Nu].sub.ext])/[[bar.Nu].sub.ext]] x 100% (11)
as given in Table II is a quantification of the relative difference between the numerical predictions of [bar.Nu] on MIII and the extrapolated value Nuext obtained from Richardson's extrapolation technique. It can be seen that the differences with grid refinement are exceedingly small and the agreement between mesh MIII and extrapolated value is extremely good; the discretisation error for [bar.Nu] is less than 0,3 %.
Table II: Effect of mesh refinement upon the mean Nusselt number for CMC-Au nanofluid ([Ra.sub.nf] = [10.sup.5], [phi] = 10 %). MI MII MIII [[bar.Nu].sub.ext] Error [bar.Nu] 6,4705 6,2260 6,1725 6,1547 0,29
In addition to the aforementioned grid-dependency study, the present simulation results have also been compared against the results of Turan et al. (18) for natural convection of the shear-thinning Power-law fluid. The comparisons between the present numerical results (obtained with mesh MIII) with the benchmark values (summarised in Table III) are extremely good and entirely consistent with our grid-dependency studies.
Table III: Comparison of the present results with the benchmark values (n = 0,80, Pr = 1000) Ra = Ra = Ra = [10.sup.4] [10.sup.5] [10.sup.6] [bar.Nu] [bar.Nu] [bar.Nu] Present 3,3601 7,6432 16,3943 study Turan et 3,3595 7,6268 16,2643 al. (18)
All these results and comparison with existing numerical data from the literature gave sufficient confidence in the present numerical procedure allowing us to proceed with simulations over the whole range of [Ra.sub.nf] and [phi]. Accounting for numerical accuracy and number of elements mesh MIII was found to be a fair compromise and all results presented in continuation were obtained with this mesh.
4. RESULTS AND DISCUSSION
4.1 Temperature and velocity flow field
It is useful to inspect the distributions of dimensionless temperature [theta] = (T - [T.sub.c])/([T.sub.H] - [T.sub.c]) and velocity [v.sub.y] * = [v.sub.y]L/[[alpha].sub.nf] in order to understand the influences of [Ra.sub.nf] and [phi] on the heat transfer characteristics during natural convection of CMC-based nanofluids in the square enclosure. The distributions of [theta] and [v.sub.y] * along the horizontal mid-plane for CMC-Au nanofluid are shown in Fig. 3 for different values of [phi].
It is evident from Fig. 3 (left) that the distributions of [theta] become increasing non-linear for increasing values of [Ra.sub.nf]. This statement is supported by the data plotted in Fig. 3 (right) which demonstrates that the magnitude of the velocity component increases significantly with increasing [Ra.sub.nf]. For a given set of values of [phi] an increase in [Ra.sub.nf] gives rise to strengthening of buoyancy forces in comparison to viscous forces which can be seen from Fig. 3 (right) where the magnitude of [v.sub.y] * increases with increasing [Ra.sub.nf]. As the convective transport strengthens with increasing [Ra.sub.nf] the distribution of [theta] becomes significantly more non-linear with increasing [Ra.sub.nf] (see the profiles in Fig. 3 (left)). At [Ra.sub.nf] = [10.sup.3] the distribution of [theta] is completely linear and the vertical velocity component [v.sub.y] * is essentially negligible due to very weak flow as the effects of buoyancy forces are dominated by viscous effects. Under this circumstance, the heat transfer takes place principally by conduction across the square enclosure. The effects of buoyancy force strengthens relative to the viscous force with increasing [Ra.sub.nf], which in turn augments heat transfer by convection due to stronger buoyancy-driven flow with higher vertical velocity magnitude. This effect is clearly evident from Fig. 3, which indicates that the vertical velocity magnitude does indeed increase with increasing [Ra.sub.nf] and that the distribution of non-dimensional temperature becomes increasingly non-linear with strengthening of convective transport for higher values of [Ra.sub.nf].
In addition, it is interesting to notice that heat and momentum characteristics at given [Ra.sub.nf] are independent of the volume fraction [phi] (see Fig. 3) and the type of nanoparticles (see Fig. 4). This finding is a reflection of the nanofluid Prandtl number values considered in the present study. Its value decreases with increasing [phi] (as shown in Fig. 2) and for this range of [Pr.sub.nf] values (i.e. [Pr.sub.nf] [much greater than] 1) the hydrodynamic boundary layer thickness remains much greater than the thermal boundary layer thickness and as a result a change in [phi] (and therefore in [Pr.sub.nf]) principally modifies the relative balance between viscous and buoyancy forces so the heat transport in the thermal boundary layer gets only marginally affected.
Fig. 4 shows that the effects of convection strengthen with increasing [Ra.sub.nf] which is reflected in progressively curved isotherms for higher values of nanofluid Rayleigh number. It is clear from Fig. 4 that the thermal boundary layer thickness decreases with increasing [Ra.sub.nf]. As a results, the thinning of boundary layer thickness for larger values of [Ra.sub.nf] acts as to enhance the magnitude of wall heat flux for the vertical walls, which gives rise to an increase in the mean Nusselt number.
4.2 Mean Nusselt number
The variation of the mean Nusselt number [bar.Nu] with volume fraction [bar.Nu] and nanofluid Rayleigh number [Ra.sub.nf] for CMC-Au nanofluid is shown in Fig. 5 which indicates that [bar.Nu] increases monotonically with increasing [Ra.sub.nf]. Moreover, it is clear that the changes in the volume fraction [bar.Nu] do not affect the distributions of [theta] and [v.sub.y] * (see Figs. 3 and 4) and therefore [phi] (and consequently [Pr.sub.nf]) does not have a major influence on the value of [bar.Nu] in the range of volume fraction considered here. In addition, the same variations of [bar.Nu] with [phi] and [Ra.sub.nf] have been observed (but not shown due to the replication of results) for all CMC-based nanofluids (i.e. Au, [Al.sub.2][O.sub.3], Cu and Ti[O.sub.2]) considered in the present study.
Previous observations are consistent with earlier findings in the context of Newtonian as well as non-Newtonian Power-law fluids (18) which demonstrated negligible Pr dependence of [theta] and [bar.Nu] for Pr [much greater than] 1 (Fig. 5). For large values of Pr the hydrodynamic boundary layer thickness remains much greater than the thermal boundary layer thickness thus the transport characteristics are primarily driven by buoyancy and viscous forces which is reflected in the weak Prandtl number dependence of [bar.Nu] for large values of Pr.
4.3 Heat transfer rate enhancement
The relative enhancement of the heat transfer rate E (as defined by Eq. (7)) is plotted with respect to the Au nanoparticles volume fraction [phi] at two values (i.e. minimum and maximum) of the nanofluid Rayleigh number [Ra.sub.nf] (Fig. 7 (left)) and to the volume fraction of all nanoparticles investigated in the present study (Figure 7 (right)). Considering the whole range of the nanofluid Rayleigh number, the Figure 7 (left) illustrates that the enhancement of the heat transfer rate increases with increasing [phi], but is practically unaffected by the value of [Ra.sub.nf]. This is consistent with the fact that the ratio of the nanofluid and the base fluid mean Nusselt number [[bar.Nu].sub.nf]/[[bar.Nu].sub.bf] = 1 remains constant for all values of the nanoparticles volume fraction [phi] at given nanofluid Rayleigh number [Ra.sub.nf] (see Fig. 5). As a result, the heat transfer rate enhancement is affected only by the ratio of the nanofluid to the base fluid thermal conductivity [k.sub.nf]/[k.sub.bf], which increases with the increasing [phi] (Fig. 8).
At any given solid volume fraction [phi] the heat transfer rate enhancement depends on the thermo-physical properties (i.e. the ratio of the nanofluid to the base fluid thermal conductivity [k.sub.nf]/[k.sub.bf]) of the nanoparticles (Fig. 7, right). Nanoparticles having the higher ratio [k.sub.nf]/[k.sub.bf] (e.g. Au and Cu nanoparticles as shown in Fig. 8) yield greater enhancement of heat transfer rate.
In the present study, steady and laminar natural convection of non-Newtonian nanofluids in a square enclosure with differentially heated vertical walls have been studied by numerical means. The effects of nanofluid Rayleigh number ([10.sup.3] [less than or equal to] [phi] [less than or equal to] [Ra.sub.nf][10.sup.6]) and solid volume fraction (0 % [less than or equal to] [phi] [less than or equal to] 10%) on heat and momentum transport have been systematically investigated.
The influence of computational grid refinement on the present numerical predictions was studied throughout the examination of grid convergence for CMC-Au nanofluid at [Ra.sub.nf] = [10.sup.5] and [phi] = 10 %. By utilizing extremely fine meshes the resulting discretisation error for the mean Nusselt number is below 0,3 %.
Numerical method was validated for the case of non-Newtonian Power-law fluid natural convection in a square cavity for which the results are available in the open literature. Remarkable agreement of present results with the benchmark results yields sufficient confidence in present numerical procedure and results.
Highly accurate numerical results released some important points such as:
* For a given set of values of solid volume fraction [phi] an increase in the nanofluid Rayleigh number [Ra.sub.nf] gives rise to strengthening of buoyancy forces in comparison to the viscous forces.
* Heat and momentum characteristics at given nanofluid Rayleigh number [Ra.sub.nf] are independent of the volume fraction [phi] and the type of nanoparticles.
* The mean Nusselt number [bar.Nu] monotonically increases with increasing nanofluid Rayleigh number [Ra.sub.nf] but its variation (with [phi] and [Ra.sub.nf]) is the same for all CMC-based nanofluids considered herein. In addition, changes in solid volume fraction [phi] do not affect the value of [bar.Nu].
* Enhancement of the heat transfer rate E increases with increasing solid volume fraction [phi], but is practically unaffected by the value of nanofluid Rayleigh number [Ra.sub.nf].
* Nanoparticles having the higher ratio [k.sub.nf]/[k.sub.bf] yield greater enhancement of heat transfer rate.
The research leading to these results was carried out within the framework of a research project "Production technology of Au nano-particles" (L2-4212) and has received funding from the Slovenian Research Agency (ARRS).
(1.) Neslusan, M.; Mrkvica, I.; Cep, R.; Kozak, D.; Konderla, R. (2011). Deformations after heat treatment and their influence on cutting process, Technical Gazette, Vol. 18, No. 4, 601-608
(2.) Ghosh, A.; Chattopadhyaya, S.; Hloch, S. (2012). Prediction of weld bead parameters, transient temperature distribution & HAZ width of submerged arc welded structural steel plates, Technical Gazette, Vol. 19, No. 3, 617-620
(3.) Brezocnik, M.; Buchmeister, B.; Gusel, L. (2011). Evolutionary algorithm approaches to modeling of flow stress, Materials and manufacturing processes, Vol. 26, No. 3, 501-507, doi:10.1080/10426914.2010.523914
(4.) Choi, S. U. S. (1995). Enhancing thermal conductivity of fluids with nanoparticles, Siginer, D. A.; Wang, H. P. (Eds.), Developments Applications of Non-Newtonian Flows, FED-Vol. 231/MD-Vol. 66, ASME, New York, 99-105
(5.) Rek, Z.; Rudolf, M.; Zun, Z. (2012). Application of CFD Simulation in the Development of a New Generation Heating Oven, Strojniski vestnik--Journal of Mechanical Engineering, Vol. 58, No. 2, 134-144, doi:10.5545/sv-jme.2011.163
(6.) Venko, S.; Vidrih, B.; Pavlovic, E.; Medved, S. (2012). Enhanced Heat Transfer on Thermo Active Cooling Wall, Strojniski vestnik--Journal of Mechanical Engineering, Vol. 58, No. 11, 623-632, doi:10.5545/sv-jme.2012.436
(7.) Stritih, U.; Butala, V. (2011). Energy Savings in Building with a PCM Free Cooling System, Strojniski vestnik--Journal of Mechanical Engineering, Vol. 57, No. 2, 125-134, doi:10.5545/sv-jme.2010.066
(8.) Contuzzi, N.; Campanelli, S. L.; Ludovico, A. D. (2011). 3D finite element analysis in the selective laser melting process, International Journal of Simulation Modelling, Vol. 10, No. 3, 113-121, doi:10.2507/IJSIMM10(3)1.169
(9.) De Vahl Davis, G. (1983) Natural convection of air in a square cavity: a bench mark numerical solution, International Journal for Numerical Methods in Fluids, Vol. 3, No. 3, 249-264, doi:10.1002/fld.1650030305
(10.) Ternik, P.; Rudolf, R.; Zunic, Z. (2012). Numerical study of heat transfer enhancement of homogeneous water-Au nanofluid under natural convection, Materials and Technology, Vol. 46, No. 3, 257-261
(11.) Ternik, P.; Rudolf, R. (2012). Heat transfer enhancement for natural convection flow of water-based nanofluids in a square enclosure, International Journal of Simulation Modelling, Vol. 11, No. 1, 29-39, doi:10.2507/IJSIMM11(1)3.198
(12.) Abu-Nada, E.; Oztop, H. F. (2009). Effects of inclination angle on natural convection in enclosures filled with Cu-water nanofluid, International Journal of Heat and Fluid Flow, Vol. 30, No. 4, 669-678, doi:10.1016/j.ijheatfluidflow.2009.02.001
(13.) Oztop, H. F.; Abu-Nada, E.; Varol, Y.; Al-Salem, K. (2011). Computational analysis of non-isothermal temperature distribution on natural convection in nanofluid filled enclosures, Superlattices and Microstructures, Vol. 49, No. 4, 453-467, doi:10.1016/j.spmi.2011.01.002
(14.) Mahrood, M. R. K.; Etemad, S. G.; Bagheri, R. (2011). Free convection heat transfer of non Newtonian nanofluids under constant heat flux condition, International Communications in Heat and Mass Transfer, Vol. 38, No. 10, 1449-1454, doi:10.1016/j.icheatmasstransfer.2011.08.012
(15.) Cheng, C.-Y. (2012). Free convection of non-Newtonian nanofluids about a vertical truncated cone in a porous medium, International Communications in Heat and Mass Transfer, Vol. 39, No. 9, 1348-1353, doi:10.1016/j.icheatmasstransfer.2012.08.004
(16.) Moraveji, M. K.; Haddad, S. M. H.; Darabi, M. (2012). Modeling of forced convective heat transfer of a non-Newtonian nanofluid in the horizontal tube under constant heat flux with computational fluid dynamics, International Communications in Heat and Mass Transfer, Vol. 39, No. 7, 995-999, doi:10.1016/j.icheatmasstransfer.2011.05.003
(17.) Bilus, I.; Ternik, P.; Zunic, Z. (2011). Further contributions on the flow past a stationary and confined cylinder: Creeping and slowly moving flow of Power law fluids, Journal of Fluids and Structures, Vol. 27, No. 7, 1278-1295, doi:10.1016/j.jfluidstructs.2011.06.004
(18.) Turan, O.; Sachdeva, A.; Chakraborty, N.; Poole, R. J. (2011). Laminar natural convection of power-law fluids in a square enclosure with differentially heated side walls subjected to constant temperatures, Journal of Non-Newtonian Fluid Mechanics, Vol. 166, No. 17-18, 1049-1063, doi:10.1016/j.jnnfm.2011.06.003
(19.) Turan, O.; Sachdeva, A.; Poole, R. J.; Chakraborty, N. (2012). Laminar Natural Convection of Power-Law Fluids in a Square Enclosure With Differentially Heated Sidewalls Subjected to Constant Wall Heat Flux, Journal of Heat Transfer, Vol. 134, No. 12, 122504/1-15, doi:10.1115/1.4007123
(20.) Roache, P. J. (1994). Perspective: A method for uniform reporting of grid refinement studies, Journal of Fluids Engineering, Vol. 116, No. 3, 405-413, doi:10.1115/1.2910291
Ternik, P.* & Rudolf, R.**, ***
* Ternik Primoz--Private Researcher, Bresterniska ulica 163, 2354 Bresternica, Slovenia
** University of Maribor, Faculty of Mechanical Engineering, Smetanova 17, 2000 Maribor, Slovenia
*** Zlatarna Celje d.d., Kersnikova ulica19, 3000 Celje, Slovenia
E-Mail: email@example.com, firstname.lastname@example.org
|Printer friendly Cite/link Email Feedback|
|Author:||Ternik, P.; Rudolf, R.|
|Publication:||International Journal of Simulation Modelling|
|Date:||Mar 1, 2013|
|Previous Article:||Calendar of events.|
|Next Article:||Applying swarm intelligence to design the reconfigurable flow lines.|