Abstract
The magnetohydrodynamic (MHD) Couette flow of two immiscible fluids in a horizontal channel with isothermal walls in the presence of an applied electric and inclined magnetic field has been investigated in the paper. Both fluids are electrically conducting, while the channel plates are electrically insulated. The general equations that describe the discussed problem under the adopted assumptions are reduced to ordinary differential equations, and closed-form solutions are obtained in both fluid regions of the channel. Separate solutions with appropriate boundary conditions for each fluid have been obtained, and these solutions have been matched at the interface using suitable matching conditions. The analytical results for various values of the Hartmann number, the angle of magnetic field inclination, loading parameter, and the ratio of fluid heights have been presented graphically to show their effect on the flow and heat transfer characteristics.
1. Introduction
The flow and heat transfer of electrically conducting fluids in channels and circular pipes under the effect of a transverse magnetic field occurs in magnetohydrodynamic (MHD) generators, pumps, accelerators, and flowmeters and have applications in nuclear reactors, filtration, geothermal systems, and others.
The interest in the outer magnetic field effect on heat-physical processes appeared seventy years ago. Blum et al. [1] carried out one of the first works in the field of mass and heat transfer in the presence of a magnetic field. The flow and heat transfer of a viscous incompressible electrically conducting fluid between two infinite parallel insulating plates have been studied by many researchers [2–6] due to its important applications in the further development of MHD technology. Also convective heat transfer in channels has been an important research topic for the last few decades because of its applications in solar technology, safety aspects of gas cooled reactors and crystal growth in liquids, and so forth.
Yang and Yu [7] studied the problem of convective magnetohydrodynamic channel flow between two parallel plates subjected simultaneously to an axial temperature gradient and a pressure gradient numerically. The problem of an unsteady two-dimensional flow of a viscous incompressible and electrically conducting fluid between two parallel plates in the presence of a uniform transverse magnetic field has been analyzed by Bodosa and Borkakati [8] for the case of isothermal plates and one isothermal and other adiabatic. The MHD fully developed flow and heat transfer of an electrically conducting fluid between two parallel plates with temperature-dependent viscosity is studied in [9, 10]. An analytical solution to the problem of steady and unsteady hydromagnetic flow of viscous incompressible electrically conducting fluid under the influence of constant and periodic pressure gradient in presence of inclined magnetic field has been obtained exactly by Ghosh [11]. Borkakati and Chakrabarty [12] investigated the unsteady free convection MHD flow between two heated vertical parallel plates in induced magnetic field. Analytical investigation of laminar heat convection in a Couette-Poiseuille flow between two parallel plates with a simultaneous pressure gradient and an axial movement of the upper plate was carried out by Aydin and Avci [13]. Recently, Singha [14] gave an analytical solution to the problem of MHD free convective flow of an electrically conducting fluid between two heated parallel plates in the presence of an induced magnetic field.
All the mentioned studies pertain to a single-fluid model. Most of the problems relating to the petroleum industry, geophysics, plasma physics, magneto-fluid dynamics, and so forth involve multifluid flow situations. Hartmann flow of a conducting fluid and a non-conducting fluid layer contained in a channel has been studied by Shail [15]. His results predicted that an increase of the order 30% can be achieved in the flow rate for suitable ratios of heights and viscosities of the two fluids. Lohrasbi and Sahai [16] studied two-phase MHD flow and heat transfer in a parallel plate channel with the fluid in one phase being conducting. These studies are expected to be useful in understanding the effect of the presence of a slag layer on heat transfer characteristics of a coal-fired MHD generator.
There have been some experimental and analytical studies on hydrodynamic aspects of the two-fluid flow reported in the recent literature. Following the ideas of Alireza and Sahai [17], Malashetty et al. [18, 19] have studied the two fluid MHD flow and heat transfer in an inclined channel, and flow in an inclined channel containing porous and fluid layer. Umavathi et al. [20, 21] have presented analytical solutions of an oscillatory Hartmann two-fluid flow and heat transfer in a horizontal channel and an unsteady two-fluid flow and heat transfer in a horizontal channel. Recently, Umavathi et al. [22] have analysed the magnetohydrodynamic Poiseuille-Couette flow and heat transfer of two immiscible fluids between inclined parallel plates.
Recent studies show that magnetohydrodynamic (MHD) flows can also be a viable option for transporting weakly conducting fluids in microscale systems, such as flow inside the microchannel networks of a lab-on-a-chip device [23, 24]. In microfluidic devices, multiple fluids may be transported through a channel for various reasons. For example, increase in mobility of a fluid may be achieved by stratification of a highly mobile fluid or mixing of two or more fluids in transit may be designed for emulsification or heat and mass transfer applications. In that regard, magnetic field-driven micropumps are in increasing demand due to their long-term reliability in generating flow, absence of moving parts, low power requirement, flow reversibility, feasibility of buffer solution manipulation, and mixing efficiency [25, 26].
MHD flows inside channels can be propelled in many different ways, for example, in electromagnetohydrodynamics (EMHDs) axial flow along a channel is generated by the interaction between the magnetic field and an electric field acting normal to it. Regardless of the purpose of a multifluid EMHD flow, it is important to understand the dynamics of interfaces between the fluids and its effect on the transport characteristics of the system. Keeping in view the wide area of practical importance of multifluid flows as mentioned above, it is the objective of the present study to investigate the MHD Couette flow and heat transfer of two immiscible fluids in a parallel-plate channel in the presence of applied electric and inclined magnetic fields.
2. Mathematical Model
As mentioned in the introduction, the problem of the EMHD Couette two fluid flow has been considered in this paper. The fluids in the two regions have been assumed immiscible and incompressible, and the flow has been steady, one-dimensional, and fully developed. Furthermore, the two fluids have different kinematic viscosities and and densities and . The transport properties of the two fluids have been taken to be constant. The analytical solutions for velocities, magnetic field, and temperature distributions have been obtained and computed for different values of the characteristic parameters. The physical model, shown in Figure 1, consists of two infinite parallel plates extending in the and -direction. The upper plate moves with constant velocity in longitudinal direction. The region I has been occupied by a fluid of viscosity , electrical conductivity , thermal conductivity , and specific heat capacity , and the region II has been filled by a layer of different fluid of viscosity , thermal conductivity , specific heat capacity , and electrical conductivity .
A uniform magnetic field of the strength has been applied in the direction making an angle to the vertical line and due to the fact that the fluid motion magnetic field of the strength has been induced along the lines of motion.
The fluid velocity, treating the problem as a monodimensional, and the magnetic field distributions for the case of inclined and induced magnetic field [8, 11, 27, 28] are where is magnetic field vector and . The upper and lower plates have been kept at the two constant temperatures and , respectively, and the plates are electrically insulated. The described MHD two fluid flow problem is mathematically presented with a continuity equation: momentum equation: general magnetic induction equation: and an energy equation: where: In previous equations and in following boundary conditions, used symbols are common for the theory of MHD flows: -time, -specific heat capacity, -velocity in longitudinal direction, -thermodynamic temperature, -magnetic permeability and -dissipation function. The third term on the right hand side of (2.3) is the magnetic body force, and is the current density vector due to the magnetic field and electric field defined by where is the vector of the applied electric field.
Finally the continuity, momentum, and induction equation written in the classic quasi-static low magnetic Reynolds number approximation [29, 30] takes the following form: where
The fluid and thermal boundary conditions have been unchanged by the addition of electromagnetic fields. The no-slip conditions require that the fluid velocities are equal to the plate’s velocities, and boundary conditions on temperature are isothermal conditions. In addition, the fluid velocity, sheer stress, induced magnetic field, induced magnetic flux (induced currents at the interface of conductors [29]), temperature, and heat flux must be continuous across the interface . Equations which represent these conditions are
3. Velocity and Magnetic Field Distribution
The governing equation for the velocity in regions I and II can be written as: where suffix () represents the values for regions I and II, respectively. The equation for the magnetic field induction in the regions I and II can be written as It is convenient to transform (3.1) and (3.2) to a nondimensional form. The following transformations have been used: With the above nondimensional quantities, the governing equations become The nondimensional form of the boundary and interface conditions (2.12) to (2.17) becomes The solutions of (3.4) with boundary and interface conditions have the following forms: where
4. Temperature Distribution
Once the velocity distributions were known, the temperature distributions for the two regions have been determined by solving the energy equation subject to the appropriate boundary and interface conditions (2.18)–(2.20). In the present problem, it has been assumed that the two walls have been maintained at constant temperatures. The term involving in the energy equation (2.10) drops out for such a condition. The governing equation for the temperatures and in region I and II is then given by In order to nondimensionalize previous equation, the following transformations have been used beside the already introduced (3.3): With the above, nondimensional quantities (4.1) for regions I and II becomes: where In the nondimensional form, the boundary conditions for temperature and heat flux at the interface becomes The solution of (4.3) with boundary and interface conditions has the following form: where
5. Results and Discussion
Recent technological trends show that the use of external fields to generate the flow inside channels, such as electrohydrodynamic, MHD, and electrokinetic flows, can be more advantageous in many microscale applications. In order to show the results of the considered MHD Couette flow problem graphically, two fluids important for technical practice (selected for the development of MHD pump under the project TR35016) have been chosen, and the parameters , and take the values of 0.677; 0.0647 and 0.025, respectively. Fluids Prandtl number is and , while Eckert number is equal to and for all the results given in Figures 2 to 13 and except in Figure 14 where it takes different values. The part of obtained results has been presented graphically in Figures 2 to 13 to elucidate the significant features of the hydrodynamic and thermal state of the flow.
Figures 2 to 4 show the effect of the magnetic field inclination angle on the distribution of velocity, temperature, and the ratio of the applied and induced magnetic field.
Figure 2 shows the effect of the angle of inclination on velocity which predicts that the velocity increases as the inclination angle increases. These results are expected because the application of a transverse magnetic field normal to the flow direction has a tendency to create a drag-like Lorentz force which has a decreasing effect on the flow velocity. Dimensionless temperature in function of angle of inclination of applied magnetic field is shown in Figure 3. In region II containing higher electrical conductivity fluid, the viscous heating is less pronounced and the influence of applied magnetic field is more expressed.
It can be seen from Figures 2 and 3 that the magnetic field flattens out the velocity and temperature profiles and reduces the flow energy transformation as the inclination angle decreases.
Figure 4 shows that the ratio of an induced and externally imposed magnetic field increases as the inclination angle of an applied field increases, for negative values of . This ratio has tendency to change the sign while decreases and have positive values.
Figures 5 to 7 depict the effect of the Hartmann number, while the electric loading factor is equal to zero (so-called short-circuited case). The influence of the Hartmann number on the velocity profiles was more pronounced in the channel region II containing the fluid with greater electrical conductivity. Figure 5 illustrates the effect of the Hartmann number on the velocity field. It was found that for large values of Hartmann number, flow can be almost completely stopped in the region II, while in region I velocity decrease is significant.
The effect of increasing the Hartmann number on temperature profiles (Figure 6) in both of the parallel-plate channel regions was in equalizing the fluid temperatures.
The influence of the Hartmann number had quite similar effect on the ratio of induced and externally applied magnetic field as shown in Figure 7.
The influence of the induced magnetic field in the considered case is not so important, but in similar flow problems where the transversal component velocity is present, the knowledge of the imposed and induced field ratio can have great significance.
Of particular significance is the analysis when the loading factor is different from zero (value of loading factor defines the system as generator, flowmeter, or pump). Figure 8 illustrates that with the increase in the absolute value of loading factor the temperature in both regions increases. In region I, viscous heating decreases while Joule heating increases, and, in region II, viscous heating increases near the lower plate and towards the middle of the channel Joule heating is more pronounced.
Figure 9 shows the effect of the loading factor on velocity, which predicts the possibility to change the flow direction. For negative values, the flow rate increases.
The obtained results show that different values of the inclination angle, the Hartmann number, and the loading factor are a convenient control method for heat and mass transfer processes.
The ratio of an induced and externally imposed magnetic field had a considerable change when the loading parameter was different from zero, especially in region II.
Figure 10 also shows a direction change of the induced field in some areas of regions I and II. This property can be used together with the change of parameters , , and in order to obtain a precise flow and heat transfer process control.
The effect of the ratio of heights of the two regions on the velocity field is shown in Figure 11. It is interesting to note that decreasing of flattens out velocity profiles and for small values, even change curves shape.
The effect of ratio of the heights of the two regions on temperature field is same as its effect on velocity field, which is evident from Figure 12. It is found that the effect of decreasing is to decrease the temperature field. It is also interesting to note that for small , the ratio of induced and externally imposed magnetic field become negligible small.
Figure 14 demonstrates the temperature distribution for different values of Eckert number Ec. It is observed that increasing values of Eckert number is to increase the temperature distribution in the flow region. Increase in Eckert number enhances the temperature because the heat energy is stored in the liquid due to the frictional heating. Temperature change is more pronounced in the region 1, while in region 2 a linear change is observed.
6. Conclusion
The problem of MHD Couette flow and heat transfer of two immiscible fluids in a horizontal parallel-plate channel in the presence of applied electric and inclined magnetic fields was investigated analytically. Both fluids were assumed to be Newtonian, electrically conducting, and have constant physical properties. Separate closed form solutions for velocity, temperature, and magnetic induction of each fluid were obtained taking into consideration suitable interface matching conditions and boundary conditions. The results were numerically evaluated and presented graphically for two fluids important for technical practice. Only part of the results are presented for various values of the magnetic field inclination angle, Hartmann number, loading parameter, and ratio of fluid heights in region I and II.
Furthermore, it was concluded that the flow and heat transfer aspects of two immiscible fluids in a horizontal channel with insulating walls can be controlled by considering different fluids having different viscosities and conductivities and also by varying the heights of regions. The obtained results show also that different values of the inclination angle, the Hartmann number, and the loading factor are a convenient control method for heat and mass transfer processes.
Acknowledgments
This paper is supported by the Serbian Ministry of Sciences and Technological development (Project no. TR 35016; Research of MHD flow in the channels, around the bodies and application in the development of the MHD pump). The authors wish to thank the reviewer for his careful, unbiased, and constructive suggestions that significantly improved the quality of this paper.