Abstract
A numerical investigation is carried out to analyze the impacts of internal heat source size, solid concentration of nanoparticles, magnetic field, and Richardson number on flow characteristics in an oppositely directed lid-driven wavy-shaped enclosure. The left and right vertical walls of the enclosure are cooled isothermally and moving with fixed velocity in upward and downward directions, respectively. The bottom wall is wavy shaped and isothermally cooled as the vertical walls while the top wall is kept adiabatic. A rectangular heater is placed horizontally in the center of the cavity. The physical problems are characterized by 2D governing partial differential equations accompanying proper boundary conditions and are discretized using Galerkin’s finite element formulation. The study is executed by analyzing different ranges of geometrical and physical parameters, namely, internal heat source length , solid concentration of nanoparticles , Hartmann’s number , and Richardson’s number . The results indicate that the overall heat transfer rate declines with the increasing length of internal heat source. The presence and rising values of solid concentration of nanoparticles cause the augmentation of heat transfer whereas the magnetic field has a negative influence and the Richardson number has a positive influence on heat transfer.
1. Introduction
The application of heat transfer mechanism in various shaped enclosures by different types of convection is swelling day by day. Besides cooling of electronic devices and heating and cooling buildings, glass industry, etc., it has many applications such as nuclear reactors, thermal storage, food processing, and solar collectors. The less thermal conductivity of conventional fluids such as water, mineral oils, and ethylene glycol is the main barrier to thermal systems enhancing effectiveness. To circumvent this restrain, researchers are continuously trying to develop some fluids that have advanced properties of heat transfer by appending metallic nanoparticles with base fluid.
Abu-Nada and Chamkha [1] examined the mixed convective CuO-water nanofluid flow in a moving lid sinusoidal enclosure. They reported that wavy wall effects strongly on heat transfer. They have shown that avg enhances in case of a wavy wall and reduces when the wall is smooth. Ali et al. [2] examined conjugate combined convection of nanofluid in an alternated baffles equipped horizontal channel. Alsabery et al. [3] investigated nanofluid flow and energy dissipation using rotating cylinders. Alsabery et al. [4] studied on nanofluid using a curvy shaped enclosure for investigation of the effect of heat source and amplitude on free convection. Alsabery et al. [5] investigated combined convection nanofluid flow using entropy generation and rotating cylinder. Armaghani et al. [6] used L-shaped enclosure to study MHD combined convection equipped with heat source/sink and hybrid nanofluid. They remarked that at solid volume fraction 0.1, entropy creation for sink heat generation is 20% less than that of source heat generation. They have also found that local for Cu-water is 29% more than that of pure water. Basak and Chamkha [7] studied nanofluid flow of free convection using square enclosure of different temperatures. Their report shows that, at , average is higher for all nanofluids and a larger enhancement of heat transfer exhibit for Cu-water and alumina-water nanofluid. Bensouici and Boudebous [8] focused on combined convection for various nanofluids cramped in a moving lid cavity furnished with a heat source. Their study shows that the nanofluids with high thermal conductivity materials transfer more heat than low thermal conductivity materials. Chamkha et al. [9] used a moving lid cavity to analyze MHD combined convection of nanofluid to determine the impacts of entropy generation, heat sink, and source. Choi and Eastman [10] were pioneers to hoist the thermal conductivity of water by appending high thermal conductivity metallic nanoparticles. Ghaneifar et al. [11] used horizontal channel equipped with Al2O3 nanofluid and double heat sources to study combined convection fluid flow. They reported that the Reynolds number has a negative impact on heat transfer. They also showed that heat transfer rate in a channel can be controlled by changing the distance between two heat sources. Ismael et al. [12] used magnetic field and heat flux for examining combined convection of nanofluid within partial slip enclosure. Kakac and Pramuanjaroenkij [13] used nanofluid to examine the enhancement of convective heat transfer. They have summarized the significant published articles which were conducted on nanofluids for the augmentation of forced convection heat transfer. Their literature investigation shows that suspending nanoparticles with pure fluids enhances the heat transfer ability of base fluids. Kalteh et al. [14] analyzed combined convection using a moving lid enclosure equipped with a heat source of triangular shaped. They reported that the mean Nusselt number upsurges with the enhancement of solid concentration of nanoparticles and downfalls with the increasing value of nanoparticle diameter. They also found that the mean Nusselt number is highest for Ag nanoparticles and lowest for TiO2 nanoparticles. Kasaeipoor et al. [15] examined nanofluid using a T-shaped enclosure to show the effect of magnetic field. Keya et al. [16] used a double-pipe heat exchanger in a moving lid enclosure to investigate conjugate combined convection flow. Mahalakshmi et al. [17] focused on MHD mixed convective fluid flow in a moving lid enclosure with a center heater for Ag-water nanofluid. They reported that higher center heater length and solid concentration of nanoparticles boost up the mean Nusselt number. They also found that the augmentation of reduces the average . Nasrin [18] analyzed MHD mixed convective fluid flow in a moving lid sinusoidal wavy walled cavity with a central heat-conducting body. She reported that average raises with the elevation of the Prandtl number and with the degradation of . She also showed that thermal boundary layer increases with the increasing value of the Richardson number. Nasrin et al. [19] investigated the hybrid and mono nanofluid using moving lid and wavy shaped enclosure. Their finding shows that the augmentation of mean is 14.7% more in the wavy enclosure than that of flat enclosure. They also found that amongst various types of hybrid nanofluid, heat transfer rate is higher in water-Cu-CNT hybrid nanofluid. Rahman et al. [20] used a circular cylinder as a heat source within rectangular enclosure to study combined convection. They showed that mean and fluid temperature in the cavity change nonmonotonically with the position of the heat-conducting cylinder in the cavity. They also showed that the maximum heat transfer rate is found for a cavity aspect ratio of 0.5 and the fluid average temperature is the lowest for a cavity aspect ratio of 2.0. Rahman et al. [21] studied nanofluid using a moving lid cavity to analyze the intensification of heat transfer. They reported the volume fraction of nanoparticle as a high-level control parameter because they noticed, under the same , nanofluid transfers more heat than the original fluid. They also reported that mean upsurges significantly by raising the Richardson number and maximum is found in the combined convection regime. Rashad et al. [22] used a moving lid and partly slip enclosure to investigate the effect of heated source and sink on MHD combined convection nanofluid flow. Roy and Basak [23] analyzed natural convection flow in an enclosure with nonuniformly heated walls. They reported that uniformly heated wall transfers more heat than nonuniformly heated wall. Saha et al. [24] used a wavy walled moving lid enclosure to investigate combined convection fluid flow. Taghikhani [25] examined MHD mixed convective nanofluid flow using moving lid enclosure furnished with Joule heating and two heat sources in wavy pattern. They reported that the increasing value of the solid volume fraction and Reynolds number significantly enhances average which boosts up convective cooling. They also showed that augmentation of the Hartmann number reduces the heat transfer rate and augments conduction heat transfer. Tiwari and Das [26] investigated the augmentation of heat transfer of a nanofluid in a double lid-driven nonuniformly heated enclosure. Zahan et al. [27] studied MHD combined convection of hybrid nanofluid subject to Joule heating using wavy walled moving lid enclosure. They reported that with the increasing value of from 0 to 50, heat transfer rate decreases by 12% and it increases by 16.89% at the undulation number three than that of flat surface.
By going through the above literature, the authors acknowledge that no investigation on mixed convection visualization with arrow surface approach into an oppositely directed lid-driven wavy walled cavity furnished with nanofluid and a central heat source of different sizes has been published. Thus, this paper is aimed at investigating the impacts of heat source length, solid concentration of nanoparticles, and magnetic field on flow characteristics under the above-mentioned circumstances for 2D laminar flow.
2. Physical Model
The treated problem is investigated inside an enclosure with height and width where filled with Cu-water nanofluid. The schematic diagram together with the rectangular Cartesian coordinate system is shown in Figure 1. The top wall is adiabatic whereas the vertical side walls and the curvy bottom wall of the cavity are maintained with same cold temperature . A rectangular heat source with temperature is positioned horizontally in the middle of the cavity. The left and right vertical walls are assumed to move in the upward and downward directions with velocity and , respectively. In standard temperature, the thermophysical properties of copper and water are given in Table 1. is the applied magnetic field along the negative -axis normal to the vertical wall, and acceleration of gravity is acting in the downward direction. The nanofluid inside the cavity is considered as laminar, incompressible, and Newtonian. The equilibrium temperature is maintained for water and copper nanoparticles.
3. Mathematical Formulation
Based on the above assumptions, the conservation equations for steady two-dimensional MHD mixed convection flow are given in a dimensional form as follows [2]:
Energy equation for rectangular center heater:
The associated conditions of heat source and boundaries are as follows:
Density of nanofluid at standard temperature is as Ali et al. [2]:
The other properties of nanofluid are given as follows:
Heat capacity as Alsabery et al. [5]: Coefficient of thermal expansion as Alsabery et al. [5]: Thermal diffusivity as Ali et al. [2]: Dynamic viscosity initiated by Brinkman as Rashad et al. [22]: Thermal conductivity introduced by Maxwell as Nasrin [18]:
To convert the governing equations (1) to (5) into nondimensional form, the following mentioned nondimensional parameters are used: where is the velocity of the vertical walls of the cavity and the nondimensional parameters , , , , and are the Grashof, Prandtl, Reynolds, Richardson, and Hartmann numbers, respectively.
Nondimensional forms of equations (1) to (5) are as follows:
Energy equation for rectangular center heater:
The nondimensional forms of the associated conditions of heat source and boundaries (6) are given as follows:
Local along vertical side walls is given as Mahalakshmi et al. [17]:
After integration (equation (17)), the average is obtained as Mahalakshmi et al. [17]:
For the wavy bottom surface, where is the spaces along -axis or -axis normal to the wavy wall. The nondimensional temperature gradient is given as where is the length of wavy wall.
4. Numerical Procedure
The governing equations have been simulated using finite element technique of Galerkin’s method [28, 29]. Finite element method is useful to solve the problems with complicated geometries where analytical solutions cannot be obtained. It divides up a very complicated problem into small elements that can be solved in relation to each other. This method has become an integral part of the design and development of numerous engineering systems. In this simulation process, the whole domain has been subdivided into definite number of triangular meshes. After that, the conservation equations (14) and (15) accompanying boundary conditions (16) are converted into a system of integral equations with the help of Galerkin’s method. By using the Newton-Raphson approach, these equations in nonlinear form are converted into linear form, and then, triangular factorization process is used for solving these linear algebraic equations to find the values of velocity components , pressure , and temperature . The convergence technique for the solution benchmark is demarcated as with as the iteration number and as . Figure 2 expresses the generated mesh of the computational domain for simulation. It consists of 25,657 domain elements and 829 boundary elements.
5. Code Verification
To verify the present solution procedure, it has been compared with the study of mixed convection using moving lid enclosure equipped with center heater and nanofluid of Mahalakshmi et al. [17]. The comparison has been concurred between current study and those of Mahalakshmi et al. demonstrated for center heater length and various Hartmann numbers at in Figures 3(a) and 3(b). These figures are almost concurrent. The comparison gives a very close propinquity with other published works. The verification with other research works boosts up the credence for graphical and numerical results of the current work.
(a)
(b)
6. Results and Discussion
MHD mixed convection in an oppositely directed lid-driven cavity filled with Cu-water nanofluid and furnished with wavy bottom wall and internal heat source is considered in this study. The nondimensional parameters characterizing the flow behavior and thermal state within the cavity are as follows: the Reynolds number is considered constant as 100, Grashof number , Richardson number , and Hartmann number . The length of internal heat source is considered as and is positioned horizontally. The whole investigation is conducted taking Cu-water nanofluid as functional fluid within the enclosure with the Prandtl number and volume fraction of nanoparticle . The flow characteristics and thermal distribution inside the cavity are investigated through streamlines and isothermal lines whereas heat transfer rate is analyzed through the average Nusselt number in three different regions, namely, forced, mixed, and free convection regions.
6.1. Effect of Center Heater Length
The effect of length on fluid flow at various Richardson numbers for different length is outlined in Figure 4 through streamlines. The flow characteristics and thermal state within the cavity are governed by forced convection in the region , by combined convection in the region , and by free convection in the region . The criterions like solid concentration of nanoparticles and are kept constant as and in pursuance of observing the unique power of heater size on fluid flow and temperature field. As seen from the figure, at , the flow within the enclosure is largely generated by forced convection which is created by the shear force emanated from the upward motion of the left wall and downward motion of the right wall. Since the shear stresses induced by the left and right walls of the enclosure are equal and opposite, streamline contours are permeated almost symmetrically by a primary clockwise rotating unicellular vortex around the center heater which is clearly seen in the first column (at ) of Figure 4. With the increasing length of , vortex pattern remains almost similar except the center of the vortices is enlarged as the takes more space in the center of the enclosure and the streamlines are shrunk near the vertical walls. Few tiny vortices are also created in the center of the vortex. With raising the Richardson number, to , where both the shear stress and buoyancy force are equally dominant, streamlines have started to change its pattern irrespective of length. And further increase of , at , the shear effect is overwhelmed by the buoyancy effect. In this region, due to the variation of temperature between heated center heater and cold side walls of the enclosure, the fluid near the center heater gets warmer and denser which flows upward and the cold fluid of higher density goes downward. Furthermore, the adiabatic wall resists the air entering into the cavity which creates the circulation of flow, and this continuous process boosts up the fluid flow and rises the number of vortices. The figure evidently expresses that the natural upward flow creates two clockwise vortices in the vertical sides of the cavity and an anticlockwise recirculation cell in the cavity center. By augmentation of length, the vortex of the middle part becomes horizontal from vertical due to the reduction of flow area. Figure 5 exhibits the arrow surface of fluid flow at different with , and . The clockwise and anticlockwise directions of the vortices are clearly seen from these figures.
The corresponding temperature distribution which is highly influenced by convection is depicted in Figure 6. As seen from the figures, the isotherms are denser at the cavity center as the exists there. In forced convection regime for all values of length, the isotherms are distributed around the which is the manifestation of sole conduction heat transfer. Similar trend continues in combined convective region with little distortion in isothermal lines, but in the natural convective region , a significant variation is discerned. In this region, the isothermal lines are raising in the upward direction creating a plume shape illustrating that more heat is being transferred due to convection. Figure 7 depicts the consequence of length on convection through the average Nusselt number in three different regimes with and . As seen from the figure, with the increasing length of , average decreases and vice versa at any constant value of . Figure 6 clearly demonstrates that by increasing the length, the isotherms occupy larger area around the transferring more heat by conduction. Enhancing conduction reduces convective heat transfer accordingly. The analyzed data of Figure 7 depicts that the mean Nusselt number is highest at , and by increasing the length to , it decreases by 30% and further increase of length to 0.6 average decreases by 4% only. To envisage the consequence of on convection, keeping the length constant, data analysis reports that heat transfer rate upsurges by 1% and 11% in the mixed and natural convection regimes, respectively, compared to that of forced convection regime at , , and . Overall heat transfer rate is maximum in the natural convection regime for and minimum in the forced convection regime for .
Figure 8 depicts the discrepancy of average thermal state with the Richardson number for different length at and . From the figure, it is clear that length has a positive relation with temperature increase. It is very logical because with the increasing value of perimeter, more heat will be disbursed and fluid temperature will increase accordingly.
6.2. Effect of Magnetic Field and Solid Concentration of Nanoparticle
Figure 9 presents the consequence of the Hartmann number on fluid flow at various volume fraction of nanoparticles through streamlines. In order to perceive the solitary effect of magnetic field and volume fraction of nanoparticles on fluid flow and temperature field, other parameters like the length and Richardson number are kept constant as and . It is evident from the figure that, at , the consequence of magnetic field is very insignificant and the streamlines are portrayed by a primary clockwise rotating unicellular vortex within the enclosure which is mostly produced by the mutual consequence of shear stress and buoyancy effect. The core of the recirculating vortex becomes anemic with swelling the value of from 10 to 30 then to 50 gradually. Since the Lorentz force generated by magnetic field resists the buoyancy force, the magnetic field and fluid velocity take action in opposite direction and the augmentation of slows down the locomotion of fluid in the enclosure which declines the strength of the circulation. Furthermore, no major change is noticed in the pattern of streamlines by escalating the value of solid concentration of nanoparticles. Figure 10 exhibits the arrow surface of streamlines at different with , , and . The clockwise and anticlockwise directions of the vortices are clearly seen from these figures.
The corresponding impact of the Hartmann number and solid concentration of nanoparticles on the thermal state is demonstrated in Figure 11. Some heat is generated within fluid of cavity by the internal friction which is created by the magnetic field. At , the isotherms are symmetrically spread around the center heater occupying the whole enclosure, but with the increasing value of , the isotherms are more scattered.
In order to present the heat transfer rate as a function of , the mean Nusselt number at the hot surface for the above-mentioned diverse solid concentration of nanoparticles is depicted in Figure 12. Here, the center heater length and Richardson number are considered constant as and . As shown in the figure, the heat transfer rate upsurges gradually by the escalating value of solid concentration of nanoparticles. The figure also expresses that the augmentation of the Hartmann number reduces the mean Nusselt number for all values of solid concentration of nanoparticles. The reason is that, by addition of Cu nanoparticles with base fluid water, thermal conductivity of fluid enhances and the increasing value of Cu nanoparticles accelerate the elevation of this thermal conductivity which boosts up the thermal transport of nanofluid inside the enclosure. The analyzed data of Figure 12 reveals that in the combined convection regime, the addition of 3%, 6%, and 9% of solid nanoparticles escalates the heat transfer rate by 8%, 16%, and 24%, respectively, at a constant magnetic field . Heat transfer rate is maximum at with and minimum at with .
7. Conclusion
MHD mixed convective Cu-water nanofluid flow has been studied in a wavy walled lid-driven enclosure. Based on the above discussion, the summary of the findings for varying several pertinent parameters is as follows: (i)The minutest center heater gives the maximum heat transfer rate in the natural convection regime(ii)Magnetic field controls heat transfer rate robustly, and it is maximum at . For base fluid, heat transfer rate is decreased by 23%, and for nanofluid, it is decreased by 25% with the increasing value of Ha from 10 to 70 at (iii)Solid concentration of nanoparticles influences strongly the flow behavior. In the mixed convection regime, the dissolved nanoparticles boost up the heat transport system and the nanofluid achieves higher average in comparison with pure base fluid. Average is upsurged by 8%, 16%, and 24% for nanofluid with solid volume concentrations 3%, 6%, and 9%, respectively, in comparison with that of pure base fluid under the same magnetic field (iv)Heat transfer rate muscularly depends on the Richardson number. It is minimum at low Richardson number . It enhances slightly in the combined convection region but vigorously in the free convection region for all center heater length, and it is maximum at (v)Center heater length has a strong influence on thermal state of fluid inside the cavity. The increasing length of center heater enhances average temperature in all the regions from forced to natural convection
Abbreviations
| : | Cold |
| : | Hot |
| : | Center heater (internal heat source) |
| : | Solid concentration of nanoparticles. |
| : | Specific heat |
| : | Gravitational acceleration (ms-2) |
| : | Convective heat transfer coefficient (Wm-2K-1) |
| : | Thermal conductivity (Wm-1K-1) |
| : | Temperature (K) |
| : | Temperature difference (K) |
| : | Velocity components (ms-1) |
| : | Nondimensional velocity components |
| : | Cartesian coordinates |
| : | Nondimensional Cartesian coordinates |
| : | Magnetic induction (Wb/m2) |
| : | Grashof number |
| : | Hartmann number |
| : | Reynolds number |
| : | Richardson number |
| : | Nusselt number |
| : | Prandtl number |
| : | Pressure (Nm-2) |
| : | Nondimensional pressure |
| : | Width of the cavity (m) |
| : | Height of the cavity (m). |
| : | Thermal diffusivity (m2s-1) |
| : | Thermal expansion coefficient (K-1) |
| : | Kinematic viscosity (m2s-1) |
| : | Electrical conductivity (Sm-1) |
| : | Dynamic viscosity (m2s-1) |
| : | Nondimensional temperature |
| : | Density (kgm-3). |
| : | Cold |
| : | Hot |
| : | Base fluid, water |
| : | Nanofluid. |
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.