International Scholarly Research Notices
Volume 2014 (2014), Article ID 256898, 11 pages
http://dx.doi.org/10.1155/2014/256898

Research Article

MHD Two-Fluid Flow and Heat Transfer between Two Inclined Parallel Plates in a Rotating System

P. Sri Ramachandra Murty and G. Balaji Prakash

Department of Mathematics, GIT, GITAM University, Visakhapatnam 530 045, India

Received 17 April 2014; Revised 18 November 2014; Accepted 23 November 2014; Published 17 December 2014

Academic Editor: Majid Ghassemi

Copyright © 2014 P. Sri Ramachandra Murty and G. Balaji Prakash. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Two-phase magnetohydrodynamic convective flow of electrically conducting fluid through an inclined channel is studied under the action of a constant transverse magnetic field in a rotating system. The fluids in the two phases are steady, incompressible, laminar, immiscible, and electrically conducting, having different densities, viscosities, and thermal and electrical conductivities. The transport properties of both the fluids are assumed constant. The bounding infinite inclined parallel plates are maintained at different constant temperatures, making an angle with the horizontal. Approximate solutions for velocity and temperature distributions are obtained by using a straightforward regular perturbation technique. An in-depth study has been done on the effects of rotation parameter, Hartmann number, inclination angle, the ratio of electrical conductivities, and viscosities of two fluids on the flow. It is observed that the effect of increasing rotation is to decrease the primary velocity. Further it is noticed that as the rotation increases, the secondary velocity increases for smaller rotation, while for larger rotation it decreases. It is also found that the temperature distribution decreases as the rotation increases.

1. Introduction

The study of convective Hartmann flow and heat transfer between two parallel plates is receiving considerable interest in the current literature. Shail [1] studied the problem of Hartmann flow of a conducting fluid in a horizontal channel of insulated plates with a layer of nonconducting fluid overlying a conducting fluid in a two-fluid flow. Rudraiah et al. [2] and Umavathi [3] have done a detailed analysis on free and forced convective heat transfer to an electrically conducting fluid in a channel. Lohrasbi and Sahai [4], Malashetty and Leela [5, 6], and Raju and Murty [7] have studied the Hartmann flow characteristics of two fluids in a horizontal channel. Seth et al. [8] have presented the detailed analysis on Hartmann flow in a rotating system in the presence of inclined magnetic field with Hall effects. Chauhan and Rastogi [9, 10] have considered Hall current and heat transfer effects on MHD flow and MHD Couette flow in a channel partially filled with a porous medium in a rotating system.

Numerous publications dealing with both the experimental and the theoretical aspects of the two-phase flow systems with or without considering the heat transfer problems associated with MHD power generators, MHD devices, and thermonuclear power generations have appeared in the literature. Abdul Mateen [11, 12] has studied the magnetohydrodynamic flow and transient magnetohydrodynamic flow of two immiscible fluids through a horizontal channel. Raju and Nagavalli [13] have studied the MHD two-layered unsteady flow and heat transfer through a horizontal channel in the presence of applied magnetic and electric fields in a rotating system. Recently, Murty and Linga Raju [14] have investigated magnetohydrodynamic two-phase flow and heat transfer between two parallel porous walls in a rotating system.

Basically, the inclined geometry has enormous applications in the heat transfer technology like solar collector. Though such problems with inclined geometry are close to realistic practical situations, much attention has not been given to handling them except the studies by Malashetty and Umavathi [15], Malashetty et al. [16], S. Daniel and Y. S. Daniel [17], and Murty and Prakash [18]. Therefore, in this paper, we have studied two-phase magnetohydrodynamic convective flow of electrically conducting fluid through an inclined channel under the action of a constant transverse magnetic field when rotated by an angular velocity about an axis perpendicular to the plates.

2. Mathematical Formulation

A steady laminar and fully developed two-phase magnetohydrodynamic convective flow driven by a common constant pressure gradient and temperature gradient has been considered in the presence of a constant magnetic field applied transversely to the direction of the flow. The physical configuration is shown in Figure 1, which consists of two infinite inclined parallel plates maintained at different constant temperatures, extending in the and directions, making an angle with the horizontal. The regions and are occupied by two different electrically conducting incompressible fluids having density , viscosity , electrical conductivity , and thermal conductivity . The whole system is rotated with an angular velocity in a counterclockwise direction about -axis perpendicular to the plates. The transport properties of both the fluids are assumed constant. The suffix   () represents the values for phases I and II, respectively.

Figure 1: Physical configuration.

With these assumptions, the governing equations of motion and energy for Boussinesq fluids as in Malashetty et al. [16] for both phases are

where and are the primary and secondary velocity components along and directions, respectively, is the temperature, is the coefficient of thermal expansion, and is the acceleration due to gravity. The fluid and the thermometric boundary conditions are unchanged by the addition of electromagnetic field. The no-slip condition requires that the velocity must be vanishing at the wall. In addition, the fluid velocity, shear stress, temperature, and heat flux must be continuous across the interface.

The boundary and interface conditions are

In making these equations dimensionless, the following transformations are used:

With the above nondimensional quantities, the governing equations (1) become

where ,  ,  ,  ,   , and , , , , and are all equal to 1 for phase I.

The nondimensional forms of the boundary and interface conditions (2) become

The asterisks have been dropped for simplicity. Further writing and , (4) can be written in complex form as

which are to be solved subject to the boundary and interface conditions:

3. Solutions

The governing equations of momentum (6) along with the energy equations (7) are to be solved subject to the boundary and interface conditions (8) for the velocity and temperature distributions. Here, we consider the Eckert number very small. Hence, the product PrEc () is very small and can be exploited to use the regular perturbation method. The solutions are assumed in the form

where are solutions for the case equal to zero. are perturbed quantities relating to , respectively. Substituting the above solution in (6) and (7) and equating the coefficients of similar powers of to zero, we get the zeroth and the first order equations as follows.

Zeroth order equations:

First order equations:

The corresponding Boundary conditions (8) reduce to

We observe that (10) and (11) are linear and coupled and therefore can be solved exactly. Here, we consider ,  ,  , and .

Solutions of the zeroth order equations (10) using boundary conditions (12) are

Solutions of the first order equations (11) using boundary conditions (13) are

The constants appearing in (15) and (16) are not given for the sake of brevity. Since the problem contains too many nondimensional parameters, for the sake of conciseness, we fix , , , , and . In the figures, all the other parameters except the varying one are chosen from the set (, Gr, , , , , ) = (2.0, 5.0, 30°, 0.5, 1.0, 2.0, 2.0).

4. Results and Discussion

Two-phase magnetohydrodynamic convective flow between two infinite inclined parallel plates in a rotating system is studied analytically. The resulting differential equations are solved using perturbation method to obtain approximate solutions for temperature distribution and primary and secondary velocity distributions. Here, we note that when , that is, in the absence of rotation, these results are in agreement with that of Malashetty et al. [16].

The effect of rotation parameter on primary velocity and secondary velocity is shown in Figures 2 and 3, respectively. From Figure 2, it is observed that the primary velocity decreases with the increase in the rotation parameter. The rotation parameter defines the relative magnitude of the Coriolis force and the viscous force in the regime. As the high magnitude Coriolis forces oppose the buoyancy force, the velocity will be decreased. From Figure 3, it is concluded that as the rotation parameter increases in (0, 2), the secondary velocity also increases but outside this range as increases, it decreases. Therefore, by increasing the rotation parameter , the secondary flow becomes oscillatory.

Figure 2: Primary velocity profiles for different values of rotation parameter .

Figure 3: Secondary velocity profiles for different values of rotation parameter .

Figures 4 and 5 show the effect of the angle of inclination on primary velocity and secondary velocity , respectively. As the angle of inclination increases, both the primary and secondary velocities increase because the magnitude of the buoyancy force increases with increase in the inclination angle. The effect of the ratio of the viscosities on primary and secondary velocities is shown in Figures 6 and 7, respectively. The smaller the value of the viscosity of the fluid in the lower phase compared to the fluid in the upper phase, the larger the primary and secondary flow fields.

Figure 4: Primary velocity profiles for different values of inclination angle .

Figure 5: Secondary velocity profiles for different values of inclination angle .

Figure 6: Primary velocity profiles for different values of ratio of viscosities .

Figure 7: Secondary velocity profiles for different values of ratio of viscosities .

Figures 8 and 9 represent the effect of the ratio of heights on primary and secondary velocities, respectively. The smaller the height of the upper phase compared to the lower phase, the larger the primary as well as secondary flow fields. The effect of the ratio of the electrical conductivity on primary velocity and secondary velocity is shown in Figures 10 and 11, respectively. We have observed that as the ratio of the electrical conductivity increases, the primary velocity increases but the secondary velocity decreases.

Figure 8: Primary velocity profiles for different values of ratio of heights .

Figure 9: Secondary velocity profiles for different values of ratio of heights .

Figure 10: Primary velocity profiles for different values of ratio of electrical conductivities .

Figure 11: Secondary velocity profiles for different values of ratio of electrical conductivities .

Figures 12 and 13 show the effect of Hartmann number on primary and secondary velocities, respectively. The effect of increasing Hartmann number is to decrease both the primary and secondary velocities. This is because an increase in applied magnetic field strength causes greater interaction between the fluid motion and the magnetic field, therefore, an increase in the Lorentz force. Since this force opposes the buoyancy force, both the velocities will be decreased. The effect of Grashof number Gr on primary and secondary velocities is shown in Figures 14 and 15, respectively. We find that an increase in the value of Grashof number Gr increases both the primary and secondary velocities.

Figure 12: Primary velocity profiles for different values of Hartmann number .

Figure 13: Secondary velocity profiles for different values of Hartmann number .

Figure 14: Primary velocity profiles for different values of Grashof number Gr.

Figure 15: Secondary velocity profiles for different values of Grashof number Gr.

The effect of rotation parameter on temperature can be seen in Figure 16. From the figure, it is evident that the temperature decreases with the increase in the rotation parameter . The rotation parameter    defines the relative magnitude of the Coriolis force and the viscous force in the regime. As the high magnitude Coriolis forces oppose the buoyancy force, the velocity will be decreased leading to a reduction in the viscous and Joule dissipation and so to a reduction in the temperature.

Figure 16: Temperature profiles for different values of rotation parameter .

Figure 17 shows the effect of the angle of inclination on temperature . As the angle of inclination increases, the temperature also increases because the magnitude of the buoyancy force increases with the increase in the inclination angle. The effect of the ratio of viscosities on the temperature is shown in Figure 18. From the figure, it is observed that the less viscous fluid in the lower phase adds the heat transfer. Figure 19 exhibits the effect of the ratio of heights on the temperature . From the figure, it is noticed that the smaller the height of the upper phase compared to the lower phase, the larger the magnitude of the temperature. Figure 20 shows the effect of the ratio of electrical conductivities on the temperature . From the figure, it is concluded that as the ratio of electrical conductivities increases, the temperature also increases.

Figure 17: Temperature profiles for different values of inclination angle .

Figure 18: Temperature profiles for different values of ratio of viscosities .

Figure 19: Temperature profiles for different values of ratio of heights .

Figure 20: Temperature profiles for different values of ratio of electrical conductivities .

Figure 21 represents the effect of Hartmann number on the temperature . From the figure, it is clear that the effect of increasing is to decrease the temperature. This is because an increase in the applied magnetic field strength causes greater interaction between the fluid motion and the magnetic field, therefore, an increase in the Lorentz force. Since this force opposes the buoyancy force, the temperature will be decreased. The effect of Grashof number Gr on temperature is shown in Figure 22; we observe from this figure that an increase in the value of Grashof number Gr increases the temperature .

Figure 21: Temperature profiles for different values of Hartmann number .

Figure 22: Temperature profiles for different values of Grashof number Gr.

Nomenclature

: | Magnetic field strength
: | Ratio of the coefficients of thermal expansion
: | Specific heat at constant pressure
Ec: | Eckert number
: | Acceleration due to gravity
: | Ratio of the heights of the two phases
Gr: | Grashof number
: | Ratio of the thermal conductivities
: | Thermal conductivities of phases I and II, respectively
: | Hartmann number
: | Ratio of the viscosities
: | Ratio of the densities
: | Nondimensional pressure gradient
: | Prandtl number
: | Reynolds number
: | Rotation parameter
: | Ratio of the electrical conductivities
: | Temperature
: | Temperature of the boundaries
: | Primary velocity
: | Secondary velocity
: | Average velocity
: | Space coordinates.

Greek Symbols

: | Coefficient of thermal expansion
: | Electrical conductivity
: | Angle of inclination
: | Density
: | Kinematic viscosity
: | Viscosity
: | Product of Prandtl number and Eckert number
: | Difference in temperature
: | Nondimensional temperature
: | Angular velocity.

Subscript

:
Value for phase.

Conflict of Interests

The authors have declared that no conflict of interests exists.

References

R. Shail, “On laminar two-phase flows in magnetohydrodynamics,” International Journal of Engineering Science, vol. 11, no. 10, pp. 1103–1108, 1973.

N. Rudraiah, V. Kumudini, and W. Unno, “Theory of nonlinear magnetoconvection and its application to solar convection problems. I, II,” Publications of the Astronomical Society of Japan, vol. 37, no. 2, pp. 183–233, 1985.

J. C. Umavathi, “A note on magneto convection in a vertical enclosure,” International Journal of Non-Linear Mechanics, vol. 31, no. 3, pp. 371–376, 1996.

J. Lohrasbi and V. Sahai, “Hall current and heat transfer effects on MHD flow in a channel partially filled with a porous medium in a rotating system,” Applied Scientific Research, vol. 45, pp. 53–66, 1988.

M. S. Malashetty and V. Leela, “Magnetohydrodynamic heat transfer in two fluid flow,” in Proceedings of the ASME/AIChE 27th National Heat Transfer conference and Exposition, pp. 28–31, July 1991.

M. S. Malashetty and V. Leela, “Magnetohydrodynamic heat transfer in two phase flow,” International Journal of Engineering Science, vol. 30, no. 3, pp. 371–377, 1992.

T. L. Raju and P. S. R. Murty, “Hydromagnetic two phase flow and heat transfer through two parallel plates in a rotating system,” The journal of the Indian Academy of Mathematics, vol. 2, pp. 343–360, 2006.

G. S. Seth, R. Nandkeolyar, and M. S. Ansari, “Hartmann flow in a rotating system in the presence of inclined magnetic field with hall effects,” Tamkang Journal of Science and Engineering, vol. 13, no. 3, pp. 243–252, 2010.

D. S. Chauhan and P. Rastogi, “Heat transfer effects on rotating MHD flow in a channel partially filled with a porous medium in a rotating system,” Turkish Journal of Engineering and Environmental Sciences, vol. 33, pp. 167–184, 2009.

D. S. Chauhan and P. Rastogi, “Heat transfer effects on rotating MHD couette flow in a channel partially filled by a porous medium with hall current,” Journal of Applied Science & Engineering, vol. 15, no. 3, pp. 281–290, 2012.

Abdul Mateen, “Magnetohydrodynamic flow and heat transfer of two immiscible fluids through a horizontal channel,” International Journal of Current Engineering and Technology, vol. 3, pp. 1952–1956, 2013.

Abdul Mateen, “Transient magnetohydrodynamic flow of two immiscible fluids through a horizontal channel,” International Journal of Engineering Research, vol. 3, no. 1, pp. 13–17, 2014.

T. L. Raju and M. Nagavalli, “MHD two-layered unsteady fluid flow and heat transfer through a horizontal channel between parallel plates in a rotating system,” International Journal of Applied mechanics and Engineering, vol. 19, no. 1, pp. 97–121, 2014.

P. S. R. Murty and T. Linga Raju, “MHD two-phase flow and heat transfer between two parallel porous walls in a rotating system,” Science Domain International, vol. 4, no. 13, pp. 1894–1907, 2014.

M. S. Malashetty and J. C. Umavathi, “Two-phase magnetohydrodynamic flow and heat transfer in an inclined channel,” International Journal of Multiphase Flow, vol. 23, no. 3, pp. 545–560, 1997.

M. S. Malashetty, J. C. Umavathi, and J. Prathap Kumar, “Convective magnetohydrodynamic two fluid flow and heat transfer in an inclined channel,” Heat and Mass Transfer, vol. 37, no. 2-3, pp. 259–264, 2001.

S. Daniel and Y. S. Daniel, “Convective flow two immiscible fluids and heat transfer with porous along an inclined channel with pressure gradient,” Research Inventy: International Journal of Engineering and Science, vol. 2, no. 4, pp. 12–18, 2013.

P. S. R. Murty and G. B. Prakash, “Heat transfer aspects on rotating MHD two-phase convective flow through an inclined channel in the presence of electric field,” Physical Science International Journal, Science Domain International, vol. 4, pp. 1260–1279, 2014.