Introduction

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. 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 Al2O3/water, Al2O3/ethylene-glycol, CuO/water, and CuO/ethylene-glycol nanofluids.

  2. 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 Al2O3/water and TiO2/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. [1012] investigated the convective heat transfer of Al2O3/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 Al2O3 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 Al2O3 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 Al2O3/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 Al2O3/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. 3.

    Boiling heat transfer: The boiling process of nanofluids was investigated experimentally by several researchers [1821]. Das et al. [18] observed the nanofluids' boiling performance deterioration. Soltani et al. [19], through experimental measurements of boiling heat transfer characteristics of Al2O3/water and SnO2/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 Al2O3/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 Al2O3/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 1.

Figure 1
figure 1

Geometry of a triangular duct.

The assumptions of the model presented in this article can be summarized as follows:

  1. 1.

    Fully developed, steady-state, laminar flow.

  2. 2.

    Constant wall temperature.

  3. 3.

    Neglecting the axial diffusion terms in the equations of conservation of momentum and energy.

  4. 4.

    Neglecting x, y direction convective terms.

  5. 5.

    Equilateral triangular duct considered for investigation.

  6. 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:

x ( μ u x ) + y ( μ u y ) = d p d z
(1)

The dimensionless parameters are defined as follows:

Z = z 2 a
(2)
Y = y 2 b
(3)
X = x 2 a
(4)
θ = T T w T i T w
(5)
u * = μ u ( d p d Z ) D h 2
(6)

The dimensionless momentum equation can be written as [32]

2 u * X 2 + ( a b ) 2 2 u * Y 2 + 4 a 2 D h 2 = 0
(7)

Consequently, the dimensionless velocity profile for equilateral triangular duct defining U = u u m is calculated as [31]

U = 15 b 2 [ b 2 Y 3 + 3 a 2 Y X + ( b 2 Y 2 + a 2 X 2 ) ( 4 27 ) b 2 ]
(8)

On the basis of the above assumptions, the energy equation for constant property flow is defined as

k eff ρ C p ( x ( T x ) + y ( T y ) ) = u d T d z
(9)

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]:

k eff = k nf + k d
(10)

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[3739]:

k d = c ( ρ C p ) nf u m v d p a
(11)

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

2 a U θ Z = [ k nf ( ρ C p ) nf u m + c ( ρ C p ) nf v d p a ( ρ C p ) nf ] 2 θ X 2 + ( a b ) 2 [ k nf ( ρ C p ) nf u m + c ( ρ C p ) nf v d p a ( ρ C p ) nf ] 2 θ Y 2
(12)

where Peclet number can be used for simplifying the equation

P e nf = 2 a u m ( ρ C p ) nf k nf
(13)

Consequently, the temperature distribution equation takes the form:

2 a U θ Z = [ 2 a P e nf + c v d p a ] 2 θ X 2 + ( a b ) 2 [ 2 a P e nf + c v d p a ] 2 θ Y 2
(14)

Dimensionless bulk temperature is defined in the following form [40]:

θ b = u * θ d A u * d A
(15)

An energy balance in a control volume in the duct will give the equation for the estimation of the local Nusselt number as follows:

N u = ( θ X ) X = 1 θ b
(16)

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]:

( ρ C p ) nf = ( 1 ν ) ( ρ C p ) f + ν ( ρ C p ) s
(17)

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, 4245], 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:

k nf = [ k s + 2 k w + 2 ( k s k w ) ( 1 + β ) 3 ν k s + 2 k w ( k s k w ) ( 1 + β ) 3 ν ] k w
(18)

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 2.

Figure 2
figure 2

Hexagonal computational cell and finite differencing along a characteristic.

Equation 14 was discretized using central differencing for ( 2 θ Y 2 ) and ( 2 θ X 2 ) and backward differencing for ( θ Z )

Boundary conditions in dimensionless form read as follows:

θ ( i , j , 1 ) = 1
(19)
θ ( 1 , j , k ) = 0
(20)
θ ( i , 1 , k ) = 0
(21)
θ ( i , m , k ) = 0
(22)
θ x ( m + 1 2 , j , k ) = 0
(23)

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 3 displays the comparison of Nusselt number values computed by Shah and London [31] with the computed values from the present simulations.

Figure 3
figure 3

Comparison between model predictions and results defined by Shah and London [31].

Results and discussion

In this section, the effect of nanoparticle's diameter, nanoparticle's concentration, and Reynolds number on heat transfer performance of Al2O3/water nanofluid is investigated.

Figure 4 shows the average Nusselt numbers versus Re for Al2O3/water nanofluid (with 1.0% volume concentration of 10-nm Al2O3 nanoparticles) and pure water. As shown in Figure 4, the slope of Nu versus Re is greater for Al2O3/water compared with pure water, which means a considerable enhancement of heat transfer due the addition of nanoparticles to the base fluid. For example, at Re = 1500, Nusselt number of water is increased from 3.47 to 4.22 with the addition of Al2O3 nanoparticles. It is known that the addition of Al2O3 nanoparticles will increase the thermal conductivity of the working fluid and hence the heat transfer capability [512]. Besides, the nanoparticles with dispersion effect and Brownian motion hit the tube wall and absorb heat, and then mix back with the bulk of the fluid to cause a better heat transfer. The presence of nanoparticles inside the fluid causes the collision between the heating surface and the particles, thereby producing higher heat transfer coefficients. This means that the addition of nanoparticles to fluid changes the flow structure so that besides the increase in thermal conductivity, dispersion and fluctuation of nanoparticles, especially near the tube wall, lead to the increase in the energy exchange rates and augment the heat transfer rate between the fluid and the tube wall [22, 23]. Moreover, the local Nu in fluid flow inside channel is related to the thickness of the thermal boundary layer, and a decrease in thermal boundary-layer thickness increases local Nu. One of the possible mechanisms responsible for the exhibition of the thermal boundary-layer thickness decrement by nanofluid is the migration of the nanoparticles due to shear action, Brownian motion, and the viscosity gradient in the cross section of the channel [47].

Figure 4
figure 4

Comparison between nanofluid and pure fluid heat transfer.

Figure 5 shows the plots of the average Nusselt number versus Re at various concentrations of Al2O3 for 10-50 nm nanoparticles. This figure indicates that the average Nusselt number increases with the concentration of the nanoparticles, and better enhancement is seen at higher Reynolds numbers. For example, at dp = 10 and Re = 400, by increasing nanoparticle's concentration from 0.01 to 0.04, the average Nusselt number increases from 2.588 to 3.345, or at higher Reynolds number (Re = 2050), the Nusselt number changes from 4.89 to 6.02. The average Nusselt number at the same diameter increases according to Reynolds number. The results illustrate that by increasing Reynolds number from 500 to 2070 at dp = 30 nm, for 0.01 volume concentration of Al2O3/water nanofluid, the average Nusselt number increases from 2.57 to 4.53.

Figure 5
figure 5

The influence of Al 2 O 3 nanoparticles' volume concentration on the Nusselt number over a range of Reynolds numbers with diameter of nanofluids in the range of 10-50 nm.

During the nanofluid flow through the channel, migration of the nanoparticles and clustering due to non-uniform shear rate across the channel's cross section affect the heat transfer performance. Taking into account the increase in thermal conductivity of the nanofluid, other factors such as chaotic movement of nanoparticles, Brownian motion, and particles' migration must also be considered in the interpretation of heat transfer performance of nanofluids [35]. An increase in the volume fraction of the nanoparticles intensifies the interaction and collision of the nanoparticles. Also, diffusion and relative movement of these particles near the channel walls lead to the rapid heat transfer from the walls to the nanofluid. In other words, increasing the concentration of the nanoparticles intensifies the mechanisms responsible for the enhanced heat transfer. Moreover, at high flow rates, the dispersion effects and chaotic movement of the nanoparticles intensify the mixing fluctuations and change the temperature profile to a flatter profile similar to turbulent flow and cause an increase in the heat-transfer coefficient. At low flow rates, clustering and agglomeration of nanoparticles may exist in the nanofluid flow, and therefore, at a low Re, a lower heat transfer enhancement can be observed.

Figure 6 displays the effect of nanoparticle's diameter on the Nusselt number for Al2O3 nanofluids of constant volume concentrations. It can be seen that the average Nusselt number increases with the decreasing size of nanoparticles at the same concentration, particularly at high concentrations. For example, by increasing the size of the nanoparticles from 10 to 50 nm in 0.02 concentration at Re = 400, the average Nusselt numbers decrease from 2.845 to 2.376. Also, at Reynolds number 2050 in 0.02 volume concentration, increasing the nanoparticle's size from 10 to 50 nm leads to a decrease in Nu from 5.273 to 4.406. Similar kind of enhancement in Nusselt number with smaller particle size was observed from the experiments conducted by Zeinali Heris et al. [10].

Figure 6
figure 6

Effect of nanoparticle's diameter on the Nusselt number for. (a) 1% volume concentration, (b) 2% volume concentration, (c) 3% volume concentration, and (d) 4% volume concentration of Al2O3 nanofluids.

Since heat transfer between the nanoparticles and the fluid takes place at the particle-fluid interface, the ratio of the surface area of nanoparticles to their volume is the most important factor for heat transfer enhancement by nanofluids. On decreasing the size of nanoparticles, the ratio of surface area-to-volume of nanoparticles increases, allowing them to absorb and transfer heat more efficiently. Moreover, for the particles with very small diameter, the particles' distribution is fairly uniform; on the other hand, by increasing the nanoparticles' mean diameter, non-uniformity on the particles' distribution becomes more important, and the nanoparticles' concentration becomes higher at the vicinity of wall for which the viscous forces are important; this causes a decrement in the heat transfer enhancement while using the nanofluid with large particle's size.

For laminar flow, the heat transfer coefficient is mainly proportional to the fluid's thermal conductivity. Also, as already mentioned, the main effects of nanoparticles inside the fluid considering the Brownian motion and fluctuation are the change in the flow structure of fluid to semi-turbulence regime, and flattening of the transverse temperature gradient in the bulk of the fluid [9], and hence enhancing convective heat transfer of the nanofluid. However, in the turbulent regime, this mechanism is not dominant, and the thermal conductivity increment is the only factor for heat transfer enhancement in turbulent flow. For the turbulent flow, as the heat transfer coefficient depends to a smaller degree on the thermal conductivity of the fluid, the effect of thermal conductivity becomes less pronounced, since, in the presence of the nanoparticles, the aforementioned heat transfer enhancement may be decreased for turbulent flow.

A detailed review of the literature revealed that it might be beneficial to avail the advantage (i.e., lower pressure drop) of using triangular cross-sectional ducts in thermal engineering systems. However, heating exchange rates will decrease through such conduits. On the contrary, the results of this preliminary numerical study revealed that this could be compensated by using nanofluids in these systems, and so this will enhance the heat transfer rates. The use of nano-sized solid particles additives suspended into the base fluid (nanofluids) is a technique recommended for the enhancement of heat transfer in the triangular ducts. There are very few correlations available to exactly predict the heat transfer performance of nanofluids, as well as correlations which include the effects of solid particles' concentration, shape, size, dispersion, and nanoparticles' random movement are not suffice. Therefore, further research on convective heat transfer of nanofluids, and more theoretical and experimental research studies are needed to clearly understand and accurately predict their hydrodynamic and thermal characteristics especially in triangular ducts.

Conclusions

In this article, the laminar flow-forced convection of Al2O3/water nanofluid in a triangular duct is studied numerically. The results indicate that the addition of nanoparticles to base fluid, besides the thermal conductivity increment, affects the structure of the flow field and leads to heat transfer enhancement, because of the dispersion and random movement of nanoparticles inside the fluid. The results obtained by the numerical solutions show that decreasing the nanoparticle's size increases Nusselt number at a specific concentration, and increasing the nanoparticles' concentration increases Nusselt number at constant particle size. The results obtained in this preliminary study indicate that, in the case of using triangular cross-sectional ducts in thermal engineering systems, because of their low pressure drop, the decrement of heating exchange rate could be compensated by the use of nanofluids in these systems. Consequently, the flow of the nanofluids through triangular conduits has both the benefits of low pressure drop and high heat transfer rate.

List of symbols

dp: Nanoparticles diameter (nm)

k d: Dispersion Thermal conductivity

kw: Thermal conductivity of water (W/mK)

ks: Thermal conductivity of solid nanoparticles (W/mK)

knf: Thermal conductivity of nanofluid (W/mK)

keff: Effective thermal conductivity (W/mK)

Cpf: Specific heat of fluid (kJ/kg K)

Cpnf: Specific heat of nanofluid (kJ/kg K)

Cps: Specific heat of nanoparticles (kJ/kg K)

Penf: Peclet number of nanofluid

T: Nanofluid local temperature (K)

Ti: Inlet temperature of nanofluid (K)

Tw: Triangle wall temperature (K)

u: Local axial velocity(m/s)

U = u u m

: dimensionless velocity profileu m : Average axial velocity (m/s)

u*: Dimensionless velocity

Greek letters

θ: Dimensionless temperature

θb: Dimensionless bulk temperature

ρs: Nanoparticles density (kg/m3)

ρf: Fluid density (kg/m3)

ρnf: Nanofluid density (kg/m3)

μ: Dynamic viscosity (kg/ms)

v: Volume fraction