Electrohydrodynamic Nanofluid Hydrothermal Treatment in an Enclosure with Sinusoidal Upper Wall

The influence of non-uniform electric filed on Fe3O4-Ethylene glycol nanofluid hydrothermal treatment in an enclosure with sinusoidal upper and moving lower walls is investigated in this study. Control Volume based Finite Element Method (CVFEM) is utilized to simulate in the presented model. Numerical investigation are conducted for the sundry parameters such as Reynolds number; nanoparticle volume fraction and supplied. Results show that supplied voltage can change the flow shape. Coulomb force causes isotherms denser near the moving wall. Heat transfer rises with augment of supplied voltage and Reynolds number. Effect of electric filed on heat transfer is more pronounced at low Reynolds number. Finally, a comparison with the existing literature is also made.


Introduction
Unlike free or forced convection, Electrohydrodynamic convection has not yet been well described.In the previous works, an experimental method is applied to investigate electric field effect.Also, there are few papers in which authors considered the combined effect of nanofluid and electric field for enhancement in heat transfer.Furthermore, there is no paper in which Control Volume based Finite

OPEN ACCESS
Element Method applied for Electrohydrodynamic nanofluid flow and heat transfer.Replacing a smooth wall with a wavy-wall is another commonly used passive method for enhancing the heat transfer in applications where space restrictions are stringent.A wavy-walled geometry is easy to generate and can provide significant heat transfer augmentation if operated in an appropriate Reynolds number range.The corrugations in the wavy-wall act as turbulence generators with the sole purpose of improving the heat transfer.The wavy surface structure adds a degree of complexity to the flow by inducing streamline curvature, flow separation and flow reattachment, thus leading to complex flow situations.The wavy wall provides considerable enhancement of heat transfer with a low pressure drop penalties, thus posing a potential to be considered in several applications such as turbine blade cooling.The target of the current article is to investigate Electrohydrodynamic nanofluid flow and force convective heat transfer via CVFEM.
One of the important active methods for improving rate of heat transfer is utilizing an electric field.The influence of electric field on buoyancy-induced flow was studied by Shu and Lai [1].In their study, electrodynamics equations are separated from fluid dynamics equations.The influence of electrode arrangements on rate of heat transfer was studied by Kasayapanand et al. [2].The influence of electric field on fluid flow over a plate has been investigated by Velkoff and Godfrey [3]. Lee et al. [4] changed the fluid flow style in a channel by adding electric field.Hsiao [5] studied the heat and mass transfer of a steady laminar boundary-layer flow of an electrically magnetic conducting fluid of second-grade subject to suction and to a transverse uniform magnetic and electric field past a semi-infinite stretching sheet.
In order to improve the thermo physical characteristics of base fluid, nano scale metal can be added in the base fluid.In which way, thermal conductivity of fluid rises and in turn heat transfer can be enhanced.Khanafer et al. [6] used Finite Volume Method to simulate nanofluid hydrothermal improvement.They indicated that adding nanoparticles leads to rise in Nusselt number.Rahman et al. [7] investigated nanofluid hydrothermal treatment in a tilted enclosure with moving wall.They proved that temperature enhances with rise of nanofluid volume fraction.Combined nanofluid free and force convection heat transfer in an enclosure with moving wall has been studied by Rahman et al. [8].They verified that rate of improvement is more pronounced at higher tilt angle.Sheikholeslami and Ellahi [9] studied three dimensional mesoscopic simulation of magnetic field effect on natural convection of nanofluid.They found that thermal boundary layer thickness increases with increase of Lorentz forces.Energy conversion for conjugate conduction, convection and radiation analysis has been performed by Hsiao [10] for free convection heat transfer.Their results were shown that the free convection effect will be produced a larger heat transfer effect better than the forced convection.Sheikholeslami et al. [11] studied the effect of thermal radiation on magnetohydrodynamics nanofluid flow and heat transfer.They found that found that concentration boundary layer thickness decreases with the increase of radiation parameter.Boutina and Bessaïh [12] presented the laminar mixed convection air-cooling of two identical heat sources, simulating electronic components, mounted in an inclined channel.Their results indicated that the Reynolds number, the inclination angle of the channel, the dimensions of the heat sources, and the spacing between them have considerable effects on the improvement of the heat transfer inside the channel.Asymmetric laminar flow and heat transfer of nanofluid between contracting rotating disks was investigated by Hatami et al. [13].Their results indicated that temperature profile becomes more flat near the middle of two disks with the increase of injection but the opposite trend is observed with increase of expansion ratio.Some relevant studies on the topic can be seen from the list of references [14][15][16][17][18][19][20][21][22].
Due to the importance of heat exchangers in many engineering applications, the subject of potential heat transfer enhancement in these devices has received substantial attention in research and practice.In the field of heat exchanger performance, it is known that utilization of wavy walls (or corrugations) can potentially lead to better heat transfer performance due to enhanced mixing of the fluid by the wavy surfaces [23,24].
Control Volume based Finite Element Method (CVFEM) is a new numerical method in which benefits of both finite element and finite volume methods are combined.This method can be used for simulation of complex geometries with multi-physics ( [25]).Sheikholeslami et al. [26] utilized an active method (magnetic field) to find the impact of Hartmann number on natural convection of nanofluid.They proved that the effect of Lorentz forces is greater for high Relight number.Ferrofluid heat transfer treatment in the presence of variable magnetic field has been studied by Sheikholeslami and Rashidi [27].They found that Nusselt number increases by considering magnetic field dependent viscosity.The lower wall has the velocity of Lid U and others are stationary.The lower wall has constant temperature 1

Problem Definition
T and the temperature of other walls is 0 T .Also the retain boundary conditions are depicted in Figure 1a.The formula of the upper sinusoidal wall is:

Mathematical Model
In order to simulate nanofluid hydrothermal treatment in existence of electric field, we should combined equations of electric fields with those of hydrothermal.The formulas of electric field are: . J 0 q t There exist two model for charge distribution: (1) conductivity model [28] and (2) mobility model [29].In first model, electro-convection rely on temperature gradient.But in second model, electro-convection is independent of temperature gradient in the liquid.In the case of free charge origination, second model is more acceptable according to experimental results.Electric current density can be defined as [30]: Where σ E → is ionic mobility, D q ∇ is diffusion [31], qV → is convection.
According to Equations ( 4) and ( 5), the equation for electric charge density can be obtained as follow: According to [30] the diffusion term can be taken negligible.Also D q ∇ in Equation ( 5) can be taken negligible and σ bq = [31].So Equation ( 5) can be considered as: In presence of electric field Coulomb forces should be added to momentum equation and Joule heating effect should be added in energy equation.So we have: and n f k are defined as [6]: (1 ) ( ) ( ) ( ) (1 ) ( ) (1 ) The thermo physical properties of the working fluid are given in Table1.Non-dimensional parameters are introduced as follows: , , In order to reach clear formulation, over bar will be deleted in next equations.So, the governing equations can be considered as follows: . 0 The formulas of vorticity and stream function are: , , , , Continuity equation has been satisfied by the stream function.By eliminating pressure between x-momentum and y-momentum, the vorticity equation can be obtained.

Numerical Method
In order to simulate this problem CVFEM is applied.We used triangular element with linear interpolation through the elements.Shape of control volume is depicted in Figure 1b.We developed FORTRAN code to solve governing equations.More information about this method can be found in [26].

Code Verification and Mesh Independency
In order to reach the mesh independent results, various meshes were examined for case of 6000 Here Γ is variable ( , , Ψ Ω Θ ).In order to verify the current results, we compared them with those of obtained in [32] (see Figure 2).Furthermore, free convection of nanofluid in a cavity was validated as shown in Figure 3.These figures proved that the current code has good accuracy.

Results and Discussion
Electric filed effect on hydrothermal behavior of nanofluid in an enclosure with moving lower and sinusoidal upper walls is presented.Table 2  The influence of Reynolds number and supplied voltage on streamlines and isotherm are shown in Figures 4-6.At Re = 3000, one main eddy and two very small eddied at upper corners of the enclosure exist in streamline.As electric filed applied the main eddy turn in to two smaller ones.Also isotherm becomes denser near the hot wall due to existence these eddies.As Reynolds number increases up to 6000, the upper right eddy becomes stronger and isotherms become denser near the bottom wall.
Since convective heat transfer is helped more successfully at a greater flow rate, Nusselt number is improved at great Reynolds numbers.So isotherm becomes more distorted at greater values of Reynolds number.As electric filed increases, the main cell convert to two eddies which are rotates opposite direction.Nusselt umber is an increasing function of supplied voltage.In absence of electric filed, Nusselt number for Re = 6000 is 1.117 times higher than that of Re = 3000 while in presence of electric filed ( ) Nusselt number for Re = 6000 is 2.298385 times lower than that of Re = 3000.Also it can be concluded that Nusselt number at 10 Δϕ = for Re = 3000, 4500 and 6000 are times higher than those of obtained at 0 Δϕ = .This observation confirms that impact of electric filed is more marked for lower Reynolds number.

Conclusions
Fe3O4-Ethylene glycol nanofluid hydrothermal treatment in existence of electric field in lid driven cavity with sinusoidal upper wall is presented.The numerical procedure is conduct via Control Volume based Finite Element Method.We reduced continuity and momentum equations to vorticity stream function.We solved the coupled equations together by considering electric filed effect.Numerical outcomes are obtained for various values of supplied voltage, nanofluid volume fraction and Reynolds number.Results show that supplied voltage causes the flow form to change.By applying electric field isotherms becomes denser near the lid wall and this observation is more market for stronger electric filed.In absence of electric filed, Nusselt number for Re = 6000 is 1.117 times higher than that of Re = 3000 while in presence of electric filed ( ) , Nusselt number for Re = 6000 is 2.298385 times lower than that of Re = 3000.Also, it can be concluded that Nusselt number at 10 Δϕ = for Re = 3000, 4500 and 6000 are times higher than those of obtained at 0 Δϕ = .Rate of augmentation in heat transfer rises with increase of strength of electric filed.Also, Nusselt number enhances with rise of Reynolds number.Improvement in Nusselt number due to existence of eclectic field becomes stronger for low Reynolds numbers.

Figure 1 Figure 1 .
Figure 1 illustrates the physical geometry along with the important parameters and mesh of the enclosure.

Figure 2 Figure 2 .
Figure 2 depicts the distribution of electric density for different Reynolds number and supplied voltage.4kV Δϕ = 10kV Δϕ = lid wall can be obtained as:

Figure 3 .
Figure 3. (a) Comparison of the local Nusselt number over the lid wall between the present results and Moallemi and Jang [30] at Re = 500, Ri = 0.4, and Pr = 1; (b) Comparison of average Nusselt number between the present results and numerical results by Khanafer et al. [6] 4 10 Gr = , 0.1 φ = and

Figure 8 .
Figure 8. Effects of Reynolds number and supplied voltage on average Nusselt number.

Table 2 .
illustrates the properties of Ethylene glycol and Fe3O4.Thermo physical properties of Ethylene glycol and nanoparticles.