Prediction of thermal distribution and fluid flow in the domain with multi-solid structures using Cubic-Interpolated Pseudo-Particle model.

A nanofluid is a suspension of very small solid particles in a continuous fluid with significant improvement of heat transfer characteristics in the main liquid. In general, in industrial equipment, the heat transfer rate can be improved with optimization of equipment including the domain structure and using the different types of nanofluids. Still, there is a big challenge to analyze the heat transfer and fluid circulation in the domain. Having nanofluids with experimental observation as using sensors and probes are destructive for the liquid stream and they are costly to observe the details of particles and the original fluid. Over the 20 years, different numerical methods have been implemented in the modeling of the heat and fluid distribution in industrial equipment containing nanofluids. Among all mathematical and numerical methods, Cubic-Interpolated Pseudo-Particle (CIP) model provides a strong potential in the prediction of the fluid structure and heat analysis, when there is a complex structure of thermal walls and high concentration of nanoparticles. However, this method is not frequently used by researchers in nanofluids analysis. In this study, the Cubic-Interpolated Pseudo-Particle model is applied to predict the flow in the square domain. different thermal walls (multi-solid structure) and hot cylindrical wall are specifically used to observe the fluid flow and heat distribution in the domain. Additionally, for a better understanding of the flow in the domain, different numbers of cylinders are used and also different amounts of nanofluid in the continuous fluid are added. The results show that adding more walls in the domain causes the change in the vortex structure. Furthermore, using nanofluid results in better heat transfer rate in the system. The CIP method is also a capable tool to predict the heat and fluid flow in the multi-solid structure domain.

1. Introduction

It is well perceived that although numerous industrial techniques are applied to increase heat transfer, the low amount of process fluid's thermal conductivity prevents great compactness and effectiveness of heat transfer equipments. A novel and promising procedure to improve the thermal conductivitiy of fuids is to apply small solid particles in the fuids [1].

To make heat transfer equipment less energy consuming, it is necessary to reach a massive increase in heat flux and diminish its dimensions. Heat transfer in fluids has a critical role in various applications such as chemical processes, power generation, reactor insolation, heating or cooling processes, microelectronics, membrane contactors, solar collectors, cooling of the electronic devices and home ventilation [2-10].

Such equipment must be designed in a way that it can provide the highest thermal efficiency and lowest size. In this regard, the small thermal conductivity of conventional coolant materials including water oil and ethylene glycol is a serious obstacle. For example, to achieve a certain thermal efficiency in a double tube heat exchanger, it is necessary to increase the heat exchanger length when a decrease occurs in the coolant’s thermal conductivity [11, 12]. Here, using nanofluids to enhance the heat transfer rate provides some advantages like size reduction in the fabricated heat transfer system, improvement in heat transfer, reduction in cost and weight of the fabricated thermal equipment, the least clogging, apparatus miniaturization, and microchannel cooling [13-16].

The low-quality heat analysis features of main liquids in comparison with most solids introduce a major problem to the efficiency of industrial equipment. Hence, numerous efforts have been made for the preparation of ultrafine solid materials with a thermal rate more significant than the commonly used liquids [17-19].

Suspending nanoparticles can be considered as a novel idea in improving the industrial equipment. Several materials such as polymeric, metallic, and non-metallic particles are added to fluids for making some slurries. Nevertheless, the typical slurries with particles in the size range of micrometers (μm) to millimeters (mm) would probably lead to some serious issues [20, 21]. The abrasion caused by these particles leads to erosion of pipelines and clogging of flow channels. Moreover, their suspensions frequently have to deal with rheological and instability issues because of the inclination of the particles to rapid settlement. Hence, despite the higher thermal conductivities of the slurries, they are not practical [22, 23].

A set of experimental works have indicated a considerable growth in heat transfer in a nanofluid containing aluminum or copper nanoparticles suspended in mineral oils or even water. Despite an increasing body of works conducted on heat transfer increase and nanofluids, many of them still have difficulties dealing with the macroscopic phenomenon occurring inside the nanofluid. Thus, they provide little information about the micro-scale enhancement mechanism of energy transport [24-26].

Generally, liquid molecular particles surround the solid nanoparticles suspended in a base liquid. These solid nanoparticles are always bombarded with the surrounding liquid particles. In addition, solid particles interact constantly. During the nanofluid preparation, some surface dispersants are applied to obtain stable nanoparticles with even distributions. These dispersants, in turn, lead to a set of complicated interactions between the liquid and solid particles. By exerting some internal and external forces, the suspended nanoparticles move irregularly even if the suspension is stationary at a microscopic scale [27, 28]. For instance, Brownian force put the nanoparticles in Brownian motion while gravitational forces lead to their sedimentation. The random motion of the nanoparticles may be accompanied by the aggregation of the particles. Here, the aggregation of the solid ones may occur because of the interaction of the nanoparticles. All nano characteristics play a key role in controlling the energy transport and flow processes in a given nanofluid [29-31]. Eastman et al conducted some experimental investigations about the influence of solid nanoparticles addition to fluid on the thermal conductiveity. They corroborated that tha addition of 5 vol.% CuO nanoparticles to water resulted in increasing its thermal conductivity by approximately 60% [32].

Many numerical methods have been used in the prediction of the nanofluids in industrial structures [33-36]. In the present study, the CIP method is employed as an efficient method for handling the heat transfer analysis in square shape domain with multi-solid structure with hot walls. Additionally, this method is used to model the nanofluid in the multi-structure domain.

2. Method

The square shape domain is used to simulate the liquid flow and heat distribution in this study. The right and left walls represent a hot and cold wall, while the top and bottom walls illustrate adiabatic walls. In this study different multi-solid walls are used for examining the CIP method in prediction of thermal distribution in the domain. The first structure is one cylindrical solid domain in the domain, and for the multi-structure domain, the triple cylinders are considered in simulations. At the beginning, the air is simulated in the domain, but after that water is mixed with Silicon Dioxide (TiO2) and add it into the domain with one object for observation of the thermal distribution and liquid stream.

To calculate the heat distribution and fluid flow, the CIP method is used in this study. Due to the flexibility of this method for nonlinear computing problems, this method splits the hyperbolic equation into two different stages including advection and non-advection. the voracity and energy equations in the two-dimensional cavity are calculated. The cavity domain has been also used in several studies to observe the liquid flow pattern and thermal distributions [37-40]. The vorticity and energy equations are written as follows:

Vorticity equation

Energy equation

The kinematics equation can be written as: where the thermal diffusivity is calculated based on:

The effective density of main liquid with ultrafine solid particles described as:

The details of all nanofluids and non-dimensional parameters can be found in [41].

3. Results and discussion

For validation of this paper, the Nu number in our study is compared with previous work. Table 1 shows a great agreement between the numerical results and existing data in the liteature for Pr = 0.7 and square shape domain. Additionally, the obtained results for Ra = 10000 are compared with literature and is percieved that our finding perfectly match with previous work (Fig 1). In this study, the analysis of the heating structure system and effect of nanofluids on the temperature distribution are investigated. To start the analysis, a hot wall at the right, a cold wall at the left and the adiabatic walls on the top and bottom are added. In addition to the walls, the hot cylinder is added in the domain. In general, having a square shape domain provides the normal fluid and heat distribution, while by adding hot cylinder in the domain, we can change this distribution which affects the optimization of heating equipment. Fig 2 shows that adding the hot cylinder causes a high-temperature region on the top and bottom right, and the fluid velocity in the domain appears with a large vortex at the left of the domain. Moreover, the maximum heat transfer appears at the left side of the cylinder. Fig 3 shows the temperature profiles along the y and x axis at the middle of the domain as well as the Nu number along with the hot wall. As the figure shows there is a low temperature at the bottom of the domain. As the position in y direction increases, the temperature increases up to the maximum value. By getting close to the cylinder, T is the maximum. However, after the cylinder area, almost near the adiabatic wall on the top, a small decrement in temperature is seen. Figs 4 and 5 show the velocity profile along the vertical and horizontal lines. As can be seen, due to having the vortex near the left wall, the positive and negative velocity representing the vortex area can be observed.

Table 1

Comparison of the average Nusselt number along the hot wall (Pr = 0.7).

Domain size	Ra	Present work	Ref. [42]	Ref. [43]
100 × 100	104	2.253	2.245	2.243
200 × 200	104	2.248	2.245	2.243

Fig 1

Comparison of the temperature at the mid-section (Y = 0.5) of the cavity by Khanafer et al. [2] (left) and the present results.

Fig 2

Temperature and velocity profiles for the hot object in the square domain.

Fig 3

Temperature and distribution at x = 0.5 and y = 0.5 for the hot object in the square domain.

Fig 4

Velocity profile in vertical direction for the hot object in the square domain.

Fig 5

Velocity profile in horizontal direction for the hot object in the square domain.

To explore the effect of the hot cylinder in this study, more cylinders are added in the domain. Adding more solid structure causes a very hot region in the domain and small vortex structure between objects (Fig 6 and Fig 7). Figs 8 and 9 illustrate the velocity profiles along the vertical and horizontal lines. The results show the small shape of the vortex in the domain due to having the multi-solid structure in the cavity.

Fig 6

Temperature and velocity profiles for three hot objects in the square domain.

Fig 7

Temperature and distribution at x = 0.5 and y = 0.5 for three hot objects in the square domain.

Fig 8

Velocity profile in vertical direction for three hot objects in the square domain.

Fig 9

Temperature and distribution at x = 0.5 and y = 0.5 for three hot objects in the square domain.

After adding nanofluid into the system with only one cylindrical object, the fluid flow and heat distribution for the different amount of nanoparticles are able to be observed. Adding more particles in the main liquid causes higher heat transfer rate in the system. Fig 10 shows the effect of the amount of nanofluids on the thermal distribution. As the graph shows, by adding more nanoparticles into the system,more heat transfer rate is observed. Fig 11 shows the effect of nanofluids on the Nu number. The higher amount of nanoparticles causes more uniform distribution in the square shape domain. In addition, the high amount of nanofluids has an effect on the vortex structure. The nanofluids result in smaller vortex domain compared with the low amount of nanoparticles (Fig 12).

Fig 10

Effect of different amount of nanofluids on the thermal distribution.

Fig 11

Effect of different amount of nanofluids on Nu number.

Fig 12

Effect of different amount of nanofluids on the liquid pattern.

Fig 13 shows the temperature distribution along the hot wall. The results show that the temperature is uniformly distributed along this wall. The amount of nanoparticles also decrease the velocity of the main liquids and generates very small vortex shape with lower maximum and minimum velocity compared with zero nanoparticle system (Fig 14).

Fig 13

Effect of different amount of nanofluids on the temperature distribution.

Fig 14

Effect of different amount of nanofluids on the velocity distribution.

4. Conclusions

Nanofluids have a significant effect on the improvement of thermal equipment in industries. This package of fluid-structure enables engineers and researchers to design and optimize a high heat transfer mechanism in the different processes, particularly heat exchanging systems. There is a need to use capable numerical methods in simulation of multi-solid walls without the usual limitation of numerical methods such as computational resources and stability of numerical methods in the complex multi-structure body. This study uses the CIP method for the prediction of multi-solid walls in the present of nanofluids in the system. The results show that the CIP method can accurately predict the nanofluid pattern and thermal distribution in the domain. Using multi-solid walls causes controlling the liquid vortex and change the thermal distribution. Additionally, adding nanoparticles in the main liquid results in higher heat transfer rates throughout the domain which increases the overall efficiency of the system. The main limitation of the current model is computing time in a large domain of simulation and coupling mechanistic models to CIP computing nodes. Additionally, calculation of turbulence spectrum, eddy structure and velocity distribution for high Reynolds number can be challenging in the domain, particularly when forcing flow such as shear flow is a main forcing scheme. To solve the computing time, the model can be written in the parallel cluster nodes, and with correct transferring information through boundaries, the computational time significantly decreases. The main suggestion for coupling mechanistic models (such as chemical reactions and deformation of lagrangian particles) in every calculation, CIP computing nodes can transfer fluid flow and thermal information into the mechanistic model. Then in every physical time we can see rate of chemical reactions or particle deformation and dynamics. For calculation of turbulence flow, Large Eddy Simulation (LES) is suggested to resolve more eddy and vortex structure in the domain. However, this method of calculation can be time consuming for large domain.