The increase of heat transfer coefficient is one of the most important technical aims for industry and researches. Also, the decrease in the pressure drop for systems that generate high fluid pressure drop is very noticeable. The aim, therefore, for achieving the optimization of heat exchangers must be always to increase the heat transfer, and simultaneously minimize the increase in the pressure drop [
1]. Increased efforts are being directed to produce heat exchangers with higher efficiency to achieve savings of energy, material, and labor [
2]. Improvements in heat transfer augmentation depend on performance and manufacturing cost [
3]. Consequently, there is an increased need for utilization of a variety of duct geometries for heat transfer applications with forced convection and internal flow [
2]. Because of the size and volume constraints in applications, such as aerospace, nuclear, biomedical engineering, and electronics, the utilization of non-circular flow passage geometries may be required, particularly, in respect of compact heat exchangers [
2]. Consequently, duct with non-circular cross section (triangular) is used in this study because of its low pressure drop, but it causes decrement of heat transfer. Therefore, for compensating this decrement, nanofluid, instead of pure fluid, was used in this study because of the former's potential to increase the heat transfer of the system. Nanofluids are created by dispersing nanometer-sized particles (<100 nm) in a base fluid such as water, ethylene-glycol, or propylene-glycol [
4].
Choi [
5] was the first person to have created fluids containing suspension of nanometer-sized particles which are called the nanofluids and disclosed their significant thermal properties through the measurement of the convective heat transfer coefficient of those fluids. Various benefits of the application of nanofluids, such as improved heat transfer, size reduction of the heat transfer system, minimal clogging, microchannel cooling, and miniaturization of systems, were achieved in his study. Since then investigations have been continued at three phases as described below:
1. Conductive heat transfer was investigated in studies by many researchers, e.g., Lee et al. [
6], which include measurement of conductive heat transfer coefficients of Al
2O
3/water, Al
2O
3/ethylene-glycol, CuO/water, and CuO/ethylene-glycol nanofluids.
2. Convective heat transfer was also studied in some published articles. For example, Pak and Cho [
7] investigated convective heat transfer in the turbulent flow regime using Al
2O
3/water and TiO
2/water nanofluids, and found that the Nusselt number of the nanofluids increased with increasing volume fraction of the suspended nanoparticles, and the increasing Reynolds number. Lee and Choi [
8] studied convective heat transfer of laminar flows of an unspecified nanofluid in microchannels, and observed a reduction in thermal resistance by a factor of 2. Nanofluids were also observed to have the ability to dissipate a heat power three times more than pure water could do. Xuan and Li [
9] measured convective heat transfer coefficient of Cu/water nanofluids, and found substantial heat transfer enhancement. For a given Reynolds number, heat transfer coefficient of nanofluids containing 1% volume Cu nanoparticles was shown to be approximately 12% higher than that of pure water.
Zienali Heris et al. [
10-
12] investigated the convective heat transfer of Al
2O
3/water and CuO/water nanofluids in circular tubes, and observed that the heat transfer coefficient was enhanced by increasing the concentration of nanoparticles in the nanofluids. However, the 20-nm Al
2O
3 nanoparticles showed an improved heat transfer performance compared with the 50-nm CuO nanoparticles, especially at high concentrations.
Maiga et al. [
13] studied numerically the heat transfer enhancement in turbulent tube flow using Al
2O
3 nanoparticles suspension. Their results showed that the inclusion of nanoparticles into the base fluid produced an augmentation of the heat transfer coefficient which has been found to increase appreciably with an increase in the concentration of particles.
Akbari and Behzadmehr [
14] investigated the developing laminar-mixed convection flow of a nanofluid consisting of Al
2O
3/water in a horizontal tube and hypothesized that the nanoparticles' concentration did not have any significant effect on the secondary flow pattern and the axial velocity.
Das and Ohal [
15] studied numerically the behavior of nanofluids inside a partially heated and partially cooled square cavity to gain insight into heat transfer and flow processes induced by a nanofluid. Ben Mansour et al. [
16] investigated the conjugate problem of developing laminar-mixed convection flow and heat transfer of Al
2O
3/water nanofluid inside an inclined tube subjected to a uniform wall heat flux. The Cu/water nanofluid-forced convective heat transfer performance in a circular tube was experimentally investigated by Zeinali Heris et al. [
17]. Based on their experimental results, they observed that the heat transfer coefficient was influenced by Peclet number, as well as by Cu nanoparticles' volume concentrations. They also stated that there was an optimum concentration for Cu nanoparticles in water, in which improved enhancement for heat transfer can be found [
17].
3. Boiling heat transfer: The boiling process of nanofluids was investigated experimentally by several researchers [
18-
21]. Das et al. [
18] observed the nanofluids' boiling performance deterioration. Soltani et al. [
19], through experimental measurements of boiling heat transfer characteristics of Al
2O
3/water and SnO
2/water Newtonian nanofluids, showed that nanofluids possess noticeably higher boiling heat transfer coefficients than those of the base fluid.
Bang and Chang [
21] studied boiling heat transfer characteristics of nanofluids with alumina nanoparticles suspended in water. They found that the addition of alumina nanoparticles caused a decrease of pool nucleate boiling heat transfer.
Noie et al. [
22] investigated heat transfer enhancement using Al
2O
3/water nanofluids in a two-phase closed thermosyphon. Experimental results showed that for different input powers, the efficiency of the TPCT increases up to 14.7% when Al
2O
3/water nanofluid was used instead of pure water. The comparison between heat transfer enhancements using metallic and oxide nanoparticles was done by Hamed Mosavian et al. [
23] whose results indicated the enhancement of heat transfer with increasing nanoparticles. Based on their experimental results, metallic nanoparticles showed better enhancement of heat transfer coefficient in comparison with oxide particles.
There are many passive cases regarding nanofluids that are still unrecognized. The most of the investigations were regarding heat transfer in circular ducts, and there is no report regarding ducts with a triangular cross section which causes lower pressure drop than the other forms of ducts. Different criteria for selecting and optimizing the heat exchanger passage geometries were outlined by Bergles [
24]. Kays and London [
25] showed that a compact heat exchanger, with a triangular cross-sectional internal flow passage, had a high ratio of heat transfer area to flow-passage volume. As far back as the late 1950s and early 1960s, Eckert et al. [
26], Sparrow [
27], Sparrow and Haji-Sheikh [
28], and Schmidt and Newell [
29] used approximate solution methods to study the pressure drop and convective heat transfer in fully developed laminar flow in ducts with cross sections having an equilateral or isosceles triangular section. Shah [
30] and Shah and London [
31] studied the heat transfer characteristics of laminar flow in a wide variety of channel shapes, including equilateral triangular, equilateral triangular with rounded corners, isosceles triangular, right triangular, and arbitrary triangular cross-sectional ducts, for an extensive range of thermal boundary conditions. Recently, Zhang [
32] has reported Nusselt numbers for laminar hydrodynamically fully developed and thermally developing flow for a uniform wall temperature condition in isosceles triangular ducts with apex angles ranging from 30° to 120°. Gupta et al. [
33] studied fully developed laminar flow and heat transfer in equilateral triangular cross-sectional ducts following serpentine and trapezoidal path. Kuznetsov et al. [
34] studied the effects of thermal dispersion and turbulence in forced convection in a composite parallel-plate channel and stated: "Although the flow in the porous region remains laminar, thermal dispersion may have a dramatic impact on heat transfer in the channel."
The aim of this article is to study the laminar-forced convective flow heat transfer of nanofluid in a triangular duct with constant wall temperature using the dispersion model.
Mathematical modeling
Laminar flow-forced convection of Al2O3/water nanofluid in a triangular duct is studied numerically. The duct configurations and coordinate system are shown in Figure .
The assumptions of the model presented in this article can be summarized as follows:
1. Fully developed, steady-state, laminar flow.
2. Constant wall temperature.
3. Neglecting the axial diffusion terms in the equations of conservation of momentum and energy.
4. Neglecting x, y direction convective terms.
5. Equilateral triangular duct considered for investigation.
6. Neglecting viscous dissipations in Cartesian coordinates
For the hydrodynamically developed and thermally developing flow, there is only one nonzero component of velocity (u), and the constitutive equations of motion reduce to a single nonlinear partial differential equation of the form:
The dimensionless parameters are defined as follows:
The dimensionless momentum equation can be written as [
32]
Consequently, the dimensionless velocity profile for equilateral triangular duct defining

is calculated as [
31]
On the basis of the above assumptions, the energy equation for constant property flow is defined as
where
keff is the effective thermal conductivity of the nanofluid. There are different approaches that are discussed in the literature for investigating the heat transfer of nanofluids. In the first approach (homogeneous model), the flow and energy equations of the base fluid are not affected by the presence of the nanoparticles. Under this assumption, both the fluid phase and solid nanoparticles flow at the same velocity and are in thermal equilibrium [
35,
36].
In the second approach, for the contribution of hydrodynamic dispersion and irregular movement of the nanoparticles, modified homogeneous model or dispersion model is adopted [
35,
36].
In this investigation, the dispersion model, in which the effect of random movement of nanoparticles inside the liquid is considered as excess terms in the heat transfer equation, is solved. The effective thermal conductivity of the nanofluid may take the following form [
35]:
where
kd is the dispersion thermal conductivity. With respect to the similarity between diffusion in porous media and nanofluid flow, the following formula has been proposed to calculate
kd [
37-
39]:
where (
c) is an unknown constant, and should be determined by matching experimental data. It depends on the diameter of the nanoparticles and flowing surface geometry. The comparison of the measured values of the nanofluid thermal conductivity with the calculated values from the proposed models indicates that the thermal dispersion is the main mechanism for enhancing fluid thermal conductivity inside channels filled with nanofluids under convective conditions. In fact, Equation 11 is a first approximation considering the dispersive effects of nanoparticles on the thermal conductivity of the nanofluid flowing through channels. According to the study of Khaled and Vafai [
40] in which heat transfer of nanofluid flow in a channel was investigated, the range of the value of c was chosen to be from 0 to 0.4. Comparing this study (triangular duct) with the channel flow in the Khaled and Vafai's investigation, the value of
c = 0.3 is used in this study. In order to examine the exact value of constant (
c), further experimental and numerical investigations are needed.
Finally, the energy equation for laminar flow in an equilateral triangular duct is
where Peclet number can be used for simplifying the equation
Consequently, the temperature distribution equation takes the form:
Dimensionless bulk temperature is defined in the following form [
40]:
An energy balance in a control volume in the duct will give the equation for the estimation of the local Nusselt number as follows:
Thermophysical properties of nanofluids
The thermophysical properties of nanofluid in the Equations 13-14 were calculated from nanoparticles and water properties using the following correlation at the mean bulk temperature [
5,
35,
41]:
Thermal conductivity is the most important parameter indicating the enhancement potential of the nanofluids. Based on the studies carried out to evaluate the thermal conductivity of the nanofluids [
5,
6,
42-
45], the theoretical models cannot predict the thermal conductivity of nanofluids. In the absence of experimental data, Yu and Choi's correlation [
44] used for the determination of the nanofluid's effective thermal conductivity is as follows:
where β is the ratio of the nanolayer thickness to the original particle radius, and β = 0.1 was used to calculate the nanofluid's effective thermal conductivity.
Validation of the simulation
In this article, the finite difference method is used for numerical solution. The discretization in the physical space (x, y) is performed by dividing the flow domain in equal triangular elements. The grid is constructed by drawing inside the triangular cross section three groups of parallel lines. The lines of each group are equally distanced and parallel to one of the three sides of the triangle.
The benefit of such a discretization is obvious. Since the boundaries of the computational domain are identical to the boundaries of the triangular cross sections of the channel, this method provides good accuracy in the numerical solution [
45,
46]. Hexagonal computational cell and finite differencing along a characteristic is shown in Figure .
Equation 14 was discretized using central differencing for

and

and backward differencing for

Boundary conditions in dimensionless form read as follows:
where m is gride number in x and y directions. The grid used in the present analysis is 560 × 560 × 100 nonuniform one with highly packed grid points in the vicinity of the tube wall and especially in the entrance region (560 in x, y direction; 100 in z direction). In order to ensure grid independence, the solution is tested for the 800 × 800 × 180 and the 500 × 500 × 90 configurations, with all of the latter giving similar values. Therefore, 560 × 560 × 100 configuration was accepted as the optimal grid size.
The duct is an equilateral triangle with a length of each side being 100 cm. The flow is laminar with the Reynolds number ranging from 100 to 2100. Because of the absence of experimental data for nanofluid in triangular ducts, in order to validate the computational model, the numerical results were compared with the theoretical data available for the conventional fluids in triangular duct as proposed by Shah and London [
31]. Figure displays the comparison of Nusselt number values computed by Shah and London [
31] with the computed values from the present simulations.