Engineering and Computer Science

Enhancement of Forced Convection Heat Transfer for SiO2 flow through Channel with Ribs at Constant Heat Flux

Khudheyer Ahmed F., Nawaf, Taha S.,
Article Date Published : 15 June 2018 | Page No.: EC-2018-35-43 | Google Scholar


Download data is not yet available.


Numerical investigation of forced convective flow in a 2-dimensional microchannel. This investigation is analyzed with nanoparticles SiO2 and water as a base fluid studying the influence of turbulence model inside multi geometrical channel (Triangular, Trapezoidal, Semi-circular, and Rectangular) by using "Finite Volume Method (FVM)". The heat flux is applied on the lower wall of channel and the upper is insulated. The diameter of nanoparticles is 20 nm. The Reynolds number ranges are from 10000 to 30000 for ratio of groove width (B) to channel height (H) was used 0.75. The volume fractions range is between 1-4%. Triangular channel score higher Nusselt number and lower friction factor than other cases against Reynolds number. When the volume fraction was increase, the Nusselt number increased and friction factor decreased, this gives 4% has the optimal properties.


Heat transfer enhancement is the main object of this paper. Internal forced convection was applied in many heat transfer systems and processes, which include heating and cooling of fluids flowing inside conduits. Heat exchangers is a device used to transfer heat between a solid object and a fluid, or between two or more fluids.

In pipes, when carry cold and hot fluids, such as oils, steam, refrigerants, and other chemicals. In automotive systems, it includes fuel lines, exhaust system, and cooling system. Internal forced convection is also in application of air ducts which is used for heating and air conditioning. In heat exchanger, convection heat transfer is enhanced in thermal systems by many engineering application by reducing size, weight, and cost.

In various chemical process plants like heat exchanger, most heat transfer frequently utilized mode is forced convection. Ribs, wires, or grooves are welded on the surface to get roughness elements to improve heat transfer [1, 2].

Numerical and experimental investigation of forced convective 2-D flow over grooved duct by Eiamsa and Promvonge [3]. 9 transverse grooves was set on the lower of incompressible air flow channel. They revealed that heat transfer can be enhanced by reverse re-circulation flow.

Groove geometry influence on turbulent heat transfer was studied by Abdulmajeed A. Ramadhan [4]. "Four different geometries of grooved tube channel were selected (triangular, circular trapezoidal, and rectangular)". The tube length was set on 180 cm, tube diameter was 5 cm, and two groove depth to tube diameter ratios (e/D) were chosen 0.1 and 0.2. The range of Reynolds number was chosen from 10000 to 20000. CFD code was used to investigate this study. It is that heat transfer was improved by 64.4%, maximum gain of 1.52 on thermal performance, and triangular grooves give highest thermal factor by as (e/D) of 0.1.

Fathinia 2012 [5] studied effect of types of fluid on 2-D grooves channel with turbulent of Reynolds number numerically. FVM, RNG k−ε simulation of model. Types of nanofluids are selected Al2O3, ZnO QUOTE _x0001_ QUOTE _x0001_ , and SiO2 QUOTE _x0001_ with respect to water was a base fluid. Volume fraction range was from 1 to 4% and diameter of nanoparticles was used from 20 nm QUOTE _x0001_ to 50 nm QUOTE _x0001_ . The result show that SiO2 QUOTE _x0001_ was best than other nanoparticles, and 20 nm is the best diameter and 4% volume fraction which improve heat transfer 114%.

Model Description

Physical Model

A "Computational fluid dynamic" CFD analysis of 2-D lower triangular ribs channel between parallel plates that illustrated in Figure 1 that show 8 ribs and 9 grooves across bottom of duct. The entire length of channel (L) is 1870 mm, and the height (H) is set to be H=40 mm. The width of groove (B) is set to be (3/4 H=30 mm). Figures 2, 3, and 4 show Trapezoidal, Circular, and Rectangular channels w.r.t diameter of circular ribs (D=2S= 40 mm). The entrance length (20H=800 mm) is regarded to create fully develop flow. The length of last rib is fixed by (10 H= 400 mm). The test section length (l) and rib height (S) is 670 mm and 20 mm respectively.

Rib length (s) is set to be (5/4 H=50 mm). The ratio of groove width to channel height is set to be (B/H=3/4).

Boundary Conditions

All boundary conditions for steady, incompressible, 2-dimensional flow. Uniform heat flux (q''=500W/m) is applied on lower wall test section length. Mean inlet velocity is 10 m/s. It is assumed that the entering flow is at temperature QUOTE _x0001_ , and QUOTE _x0001_ is the ambient temperature and zero pressure gradient at the exit, and at wall is no slip condition.

Figure 1::

Figure 2::

Figure 3::

Figure 4::

Mathematical modeling

Governing equations

Governing equations under previous conditions of "Continuity equation, Navier-Stokes equations and Energy equation" are as [6]:

Continuity equation:

∂/(∂x_i ) (ρ u_i )=0 (1)

Momentum equation:

∂/(∂x_i ) (ρu_i u_j )=-∂P/(∂x_i )+∂/(∂x_j ) [μ((∂u_i)/(∂x_j )+(∂u_j)/(∂x_i ))] +∂/(∂x_i ) (-ρ (u_i^' u_j^' ) ̅ ) (2)

Energy equation:

(∂ )/(∂x_i ) (ρu_i T)=∂/(∂x_j ) ((Γ+Γ_t ) ∂T/(∂x_j )) (3)

Where QUOTE _x0001_ is "molecular thermal diffusivity" and QUOTE _x0001_ are "turbulent thermal diffusivity" have been given as:

Γ=μ/Pr, and Γ_t=μ_t/Pr_t (4)

In eq. (2), Reynolds stress (-ρ (u_i^' u_j^' ) ̅ ) can be modeled by k–ε turbulence model. "The Boussinesq hypothesis is employed to relate the Reynolds stress to the mean velocity gradient":"-ρ (u_i^' u_j^' ) ̅=μ_t ((∂u_i)/(∂x_j )+(∂u_j)/(∂x_i )) (5)

The turbulence viscosity ( QUOTE _x0001_ ) expression is computed from an appropriate turbulence model is given as:

μ_t=ρC_μ k^2/ε" (6)

The modeled TKE is written as:

∂/(∂x_i ) (ρku_i )=∂/(∂x_j ) [(μ+(∂μ_t)/σ_k ) ∂k/(∂x_j )]+G_k-ρε (7)

TKE dissipation rate, ε, is similarly given by:

∂/(∂x_i ) (ρεu_i )=∂/(∂x_j ) [(μ+(∂μ_t)/σ_ε ) ∂ε/(∂x_j )]+C_1ε ε/k G_k-C_2ε ρ ε^2/k (8)

Where (ρε) is TKE destruction rate, while G_k its generation rate given by:

"G_k=-ρ (u_i^' u_j^' ) ̅ (∂u_j)/(∂x_i )" (9)

The k–ε model constants can be chosen to be empirical constants [7] the turbulence transport equations is in Table 1

Table 1. Constants in k–ε turbulence model transport equations

C_1ε C_2ε σ_k σ_ε Pr_t
0.09 1.44 1.92 1.0 1.3 0.9

The governing equations can be solved by utilizing CFD code FLUENT 63 with the standard k–ε turbulence model and the SIMPLE algorithm The convergence criteria of 10 and 10 are assumed for the residuals of energy and all other variables, respectively

Solution procedure

Discretization of governing equations and turbulence model by using FVM and CFD was used to solve governing equations. For convective and diffusive terms, QUICK method was applied.

The setup is completed and number of iteration 250 is enough in Setup term.

Discretization of non-linear equations were implemented implicitly. SIMPLE method was chosen to estimate pressure field by velocity-pressure coupling algorithm.

Uniform velocity profile are applied at entrance region. The constant heat flux is set on bottom (test section) at same boundary conditions. Intensity of turbulence was kept 10% at entrance region. For the present case, friction factor ( QUOTE _x0001_ ) and Nusselt number ( QUOTE _x0001_ ) are benefit parameters.

The friction factor can be estimated by pressure drop ( QUOTE _x0001_ ) along test section length ( QUOTE _x0001_ ).

f=∆P/(1/2 ρu_av^2 ) D_h/L (10)

Where D_h=(4A_c)/P is the hydraulic diameter. The average Nusselt number may be obtained by:

"Nu=(h D_h)/k_f " (11)

3.3 Thermo physical properties of nanofluids

The following equations were used to estimate the thermo physical properties of nanofluids. Table 2 below gives the thermo-physical properties of water and nanoparticle SiO_2.

ρ (kg/m^3) Cp (J/kg.K) K (W/m.K) μ (kg/m.s)
water 998.2 4182 0.6 1.003*〖10〗^(-3)
2220 745 1.5 --

Thermo-physical properties of materials

Effective Density [8]

ρ_nf=(1-ϕ) 〖 ρ〗_f+ϕρ_np (12)

The effective heat capacity is given as [8]:

〖(ρC_p)〗_nf=(1-ϕ) 〖 (ρC_p)〗_f+ϕ〖(ρC_p)〗_np (13)

Where QUOTE _x0001_ and QUOTE _x0001_ are heat capacities of base fluid and Nanoparticles, respectively.

The effective thermal conductivity equation of nanofluid is written as [8]:

k_eff=k_static+k_Brwnian (14)

k_static=k_f [((k_np+2k_f )-2ϕ(k_f-k_np ))/((k_np+2k_f )+2ϕ(k_f-k_np ) )] (15)

k_Brwnian=5⨯〖10〗^4 βϕ〖 ρ〗_f C_(p,f) ((K T)/(2ρ_np R_np ))^0.5 f(T,ϕ ) (16)

Where QUOTE _x0001_ and QUOTE _x0001_ are the thermal conductivities of nanoparticles and base fluid respectively.

Boltzmann Constant:

k=1.3807⨯ 〖10〗^(-23) J/ K (17)

f(T,ϕ )=(2.8217⨯ 〖10〗^(-2) ϕ+3.917⨯ 〖10〗^(-2) )(T/T_o )+(-3.0669⨯〖10〗^(-2) ϕ-3.391123⨯〖10〗^(-2)) (18)

For 1% ≤ ϕ ≤ 4% and 300K, where T is the fluid temperature, and T0 is the reference temperature. The effective viscosity can be obtained using following mean empirical correlations [14]:

μ_eff=μ_f⨯1/((1-34.87(d_p/d_f )^(-0.3)⨯∅^1.03 ) ) (19)

d_f=0.1*[6M/(Nπρ_(f_o ) )]^(1/3) (20)

Where M is the molecular weight of base fluid, N is the Avogadro number, and QUOTE _x0001_ is the mass density of the based fluid calculated at temperature T0=293K.

The effective thermal expansion is expressed as [8]:

〖(ρβ)〗_nf=(1-ϕ) 〖 (ρβ)〗_f+ϕ(〖ρβ)〗_np (21)

For different particle materials the 𝛽 equations are listed in Table 3.

Types of particles Concentration
SiO 2 1.9526(100 ϕ)^(-1.4594) 1%≤ϕ≤10%

QUOTE _x0001_ Values for different particles

Numerical results

For Triangular ribs channel, number of Nodes and Elements are taken of 22902 and 16719 respectively. For Trapezoidal ribs channel, number of nodes and elements were 21123 and 14461 respectively, and for semi-circular 18770 nodes and 15107 elements, and for Rectangular ribs 25712 nodes and 18199 elements. Additional increase in Nodes and Elements is lower than 1.5% variation of average Nusselt number that is taken a criterion for grid independence to save computer memory and calculation time. Figures 5 a, b, c, and d show mesh of fluid domain of (test section) in micro-channels for different shapes.

Figure 5 a::

Figure 5 b::

Figure 5 c::

Figure 5 d::

4.2 Validating the results

In order to validate the results obtained from CFD code developed in present numerical study of Triangular, Trapezoidal, and Semi-circular ribs channel, the friction factor for turbulent convective heat of nanofluids in three geometric channel are calculated and compared with the previous experimental results of smooth channel and also with Rectangular ribs channel Fathinia et, al [ 5] as shown in Fig6 According to this figure, the results are in good agreement Moreover, the average Nusselt number for this three geometrical channels compared with smooth channel and Rectangular ribs channel Fathinia et,al [ 5] as shown in Figure 7 From Fig 7, it is found that the present results achieved better and higher Nusselt number than previous studyIn order to validate the results obtained from CFD code developed in present numerical study of Triangular, Trapezoidal, and Semi-circular ribs channel, the friction factor for turbulent convective heat of nanofluids in three geometric channel are calculated and compared with the previous experimental results of smooth channel and also with Rectangular ribs channel Fathinia et, al [5] as shown in Fig. 6 According to this figure, the results are in good agreement. Moreover, the average Nusselt number for this three geometrical channels compared with smooth channel and Rectangular ribs channel Fathinia et,al [5] as shown in Figure 7. From Fig. 7, it is found that the present results achieved better and higher Nusselt number than previous study.

Figure 6::

Figure 7::

Results and discussion

The simulations present five types of channel (Triangular, Trapezoidal, Semi-circular, Rectangular ribs, and Smooth duct). Figure 6 above gives the results that friction factor decrease with increasing of Reynolds number values in different shape of channel, the results above give Smooth channel has lower and the Triangular ribs score lower than other ribs channel shapes. Figure 7 gives the Nusselt number increase with increasing Reynolds number, the results indicate that Triangular ribs channel have a higher Nusselt number than other case.

Triangular ribs channel used range of nanoparticles SiO2 volume fraction of QUOTE _x0001_ w.r.t nanoparticles diameter dp= 20 nm and water is considered as a base fluid of nanoparticles. The results of friction factor and Nusselt number vs Reynolds number compare of multi-concentrations of nanofluid with water only. Figure 8 below gives friction factor decrease with increasing of volume fraction against Reynolds number.

Figure 8::

Figure 9::

The velocity contour and vector shown in Figure 10 a, and b notice that the velocity is decreasing when it was close to ribs in the test section because of the entropy and it is reach to 13 m/s when it was close to upper channel wall. Pressure income is fully developed because of the first rectangular rib (20H).

Figure 10 a::

Figure 10 b::


The investigation is analyzed of nanoparticles SiO2 and water as a base fluid by using "Finite Volume Method (FVM)" with 20 nm as diameter of nanoparticles. The heat flux is set on lower (Test section) of the channel and the upper wall is insulated. Five geometrical rib groove channel are studied (Triangular, Trapezoidal, Semi-circular, Rectangular, and Smooth duct) The Reynolds number range is from 10000 to 30000. The groove width to channel-height to ratio (B/H) of 0.75. The results show that triangular ribs channel is the optimum selection to enhance heat transfer and fluid flow characteristics, rather than other cases. The increase of volume fraction of nanoparticles gives increase in Nusselt number and decrease in friction factor (4% is optimum choice of volume fraction).

Symbol Quantity
Ac Surface Area, m2 t Temperature, K
B Distance between grooves, mm uav Mean velocity, m/s
C? Turbulence model constant u?ij Fluctuation velocity components
Dh Hydraulic diameter of channel, mm
f Friction factor CFD Computational Fluid Dynamics
Gk TKE generation rate QUICK Quadratic Upstream Interpolation for
h Convective heat transfer coefficient Convective Kinetics scheme
H Channel height, mm RNG Random Number Generation
k Turbulent kinetic energy SIMPLE Semi Implicit Method for Pressure
kf Thermal conductivity Linked Equations
l Length of test section, mm TKE Turbulent Kinetic Energy
L Length of channel, mm Greek symbols
Nu Nusselt number ? Thermal diffusivity
p Parameter, mm ? Turbulent dissipation rate
P Static pressure, Pa ? Kinematic viscosity
Pr Prandtl number ?t Eddy viscosity
q" Heat flux, W/ m2 ? Density, kg/m3
Re Reynolds number ?ij Reynolds stress
s Rib land/length, mm ? Wall shear stress
S Rib height, mm ? Turbulent specific dissipation rate

7. Nomenclature


  1. B. E. Launder and D. B. Spalding, Mathematical Models of Turbulence. 169 S. m. Abb. London/New York 1972. Academic Press. Preis geb. $ 7.50 Szablewski W. ZAMM - Zeitschrift fĂźr Angewandte Mathematik und Mechanik.1973;:424-424. CrossRef Google Scholar
  2. The effect of nanofluids flow on mixed convection heat transfer over microscale backward-facing step Kherbeet ASh, Mohammed HA, Salman BH. International Journal of Heat and Mass Transfer.2012-oct;:5870-5881. CrossRef Google Scholar

Author's Affiliation

Copyrights & License

International Journal of Scientific Research and Management, 2018.
Creative Commons License

This work is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License.

Article Details

Issue: Vol. 6 No. 06 (2018)
Page No.: EC-2018-35-43
Section: Engineering and Computer Science

How to Cite

F., K. A., & Taha S., N. (2018). Enhancement of Forced Convection Heat Transfer for SiO2 flow through Channel with Ribs at Constant Heat Flux. International Journal of Scientific Research and Management, 6(06), EC-2018.

Download Citation

  • HTML Viewed - 133 Times
  • PDF Downloaded - 60 Times
  • XML Downloaded - 18 Times

Most read articles by the same author(s)