Induced magnetic field and viscous dissipation on flows of two immiscible fluids in a rectangular channel.

The unsteady, magneto-hydrodynamic generalized Couette flows of two immiscible fluids in a rectangular channel with isothermal walls under the influence of an inclined magnetic field and an axial electric field have been investigated. Both fluids are considered electrically conducting and the solid boundaries are electrically insulated. Approximate analytical solutions for the velocity, induced magnetic, and temperature fields have been determined using the Laplace transform method along with the numerical Stehfest's algorithm for the inversion of the Laplace transforms. Also, for the nonlinear differential equation of energy, a numerical scheme based on the finite differences has been developed. A particular case has been numerically and graphically studied to show the evolution of the fluid velocity, induced magnetic field, and viscous dissipation in both flow regions.

Introduction

The dynamics that occur in the motion of liquid–solid, liquid–liquid and liquid–gas environments are a rich source of interdisciplinary research, for example, in nanotechnology, the nuclear industry, as well as when cooling devices in microelectronics. The overall properties and applicability of multiphase systems are totally dependent on interface shapes, which influence and are influenced by each phase's flow area, where the heat fluxes are not equal since the transition in interfacial energy is considered[1,2].

Magneto-hydrodynamic (MHD) devices, generators, accelerators, and flow meters use the flow and heat transfer of electrically conducting fluids in channels and circular pipes under the influence of a transverse magnetic field, and have uses in nuclear reactors, filtration, and heat exchangers, etc. Decades prior, there was an involvement in the influence of the outer magnetic field on heat-physical mechanisms. Blum et al.[3] studied heat and mass transfer in the presence of a magnetic field and that was one of the first works. Many researchers have investigated the flow and heat transfer of a viscous, incompressible, electrically conducting fluid, between two infinite, parallel, insulating plates[4-6]. Convective heat transfer in channels has also become a common research subject in recent decades due to its applications in solar energy, gas cooled reactor protection, and crystal growth in liquids, among other things. The first research into a two-phase liquid metal magneto-fluid-mechanics generator was started by Thome[7].In this study have been analyzed effects of the magnetic field on distribution of gas in the field direction, slip ratio, and two-phase pressure drop, by using streams of sodium–potassium alloy and nitrogen that were mixed and pumped through a vertical, rectangular channel with a transverse magnetic field applied perpendicularly to the long side of the cross section. Postlethwaite and Sluyter[8] gave a summary of the heat transfer issues that MHD generators face. Lohrasbi and Sahai[9] investigated MHD two-phase flow and heat transfer in a horizontal parallel-plate channel and identified analytical solutions for velocity and temperature profiles where only one of the fluids is electrically conducting. Closed-form solutions for two-phase flow and heat transfer in a horizontal channel with both phases electrically conducting were mentioned by Malashetty and Leela[10,11].

Malashetty and Umavathi[12] recently investigated two-phase MHD flow and heat transfer in an inclined channel with buoyancy effects for the case where only one of the phases is electrically conducting. The studies mentioned above are helpful in determining the impact of slag layers on the heat transfer aspects of coal fired MHD generators.

Sai et al.[13] investigated the flow of two immiscible fluids in a channel with one permeable wall caused by an oscillatory pressure gradient. They found slip conditions on the channel wall but assumed continuous velocities and shear stresses at the fluid–fluid interface. The authors studied the flow of water and mercury under the influence of a sinusoidal time-dependent pressure gradient as a particular case. By introducing cross terms to the generalized Darcy law, Pasquier et al.[14] developed a statistical model of creeping flows that provides an explicit coupling term in both phases. Using an extension of the Buckley-Leverett theorem, the effect of the additional terminology on saturation profiles and pressure drops was studied. Nicodijevic et al.[15] have studied magneto-hydrodynamic Couette flows of two immiscible, electrically conducting Newtonian fluids in a horizontal channel. The flow is influenced by an electric field and an inclined magnetic field while, the channel walls are isothermal and insulated. The induced magnetic field has been also determined. Other interesting results of immiscible fluids flow can be found in[16-18].

Thermal management systems/devices must be able to withstand a variety of thermal boundary conditions due to their vast range of applications in various sectors. Uniform heat flux, uniform wall temperature, and insulated wall boundary conditions are common examples. Numerous studies on studying the convective heat transfer process under such boundary circumstances have been published in the literature. However, depending on the flow environment, an asymmetric thermal boundary condition may be required to investigate the system's underlying thermal transport properties. We'd like to point out that there are a few studies in the literature that analyze the thermo-hydrodynamics of both Newtonian and non-Newtonian fluids while taking the aforementioned thermal boundary conditions into account[19-21].

The purpose of this article it to investigate the unsteady, magneto-hydrodynamic generalized Couette flows of two immiscible fluids in a rectangular channel with isothermal walls under the influence of an inclined magnetic field and an axial electric field have been investigated. Both fluids are considered electrically conducting and the solid boundaries are electrically insulated. The problem is formulated into the dimensionless form and, at the interface the velocity, shear stress, induced magnetic and temperature fields are considered continuous functions. The bottom wall of channel is fixed, while the upper wall is moving with a given time-dependent velocity. The nonslip conditions on the solid boundaries are also considered. Approximate analytical solutions for the velocity, induced magnetic, and temperature fields have been determined using the Laplace transform method along with the numerical Stehfest's algorithm for the inversion of the Laplace transforms. Also, for the nonlinear differential equation of energy, a numerical scheme based on the finite differences has been developed. A particular case characterized by a time-exponential velocity of the upper wall, has been numerically and graphically studied to show the evolution of the fluid velocity, induced magnetic field, and viscous dissipation in both flow regions.

Problem formulation

The physical model shown in Fig. 1 consists of two infinite, parallel plates extending in the x-direction, and z-direction. The y-axis is perpendicular to the channel walls. The fluids in the two domains are assumed immiscible and electrically conducting. The flow is incompressible, unsteady, one-dimensional, and fully developed.

Figure 1

Physical model. Coordinate system.

The region is occupied by a fluid of viscosity , density , electrical conductivity , thermal conductivity , specific heat capacity , and the magnetic permittivity .

The region is field by a different fluid of viscosity , density , electrical conductivity , thermal conductivity , specific heat capacity , and the magnetic permittivity .

It is supposed that fluids have laminar, one-dimensional flow with the velocity fields.where is the unit vector along the x-direction.

A uniform magnetic field of the strength is applied such that the angle between vectors is equal to therefore, the applied magnetics vector field is

Due to the fluid motion, a magnetic field is induced along of the lines of motion, namelyis induced along of the lines of motion.

The total magnetic field is given by

The walls of the rectangular channel are kept at the constant temperature respectively and the plates are electrically insulated.

Also, a time-dependent electric fieldis applied to the fluids.

In the above hypotheses, the continuity equation is satisfied by the velocity fields in (1).

The electro-magneto-hydrodynamic fluids flow problem is mathematically described by[22-24]:

The momentum equation

The constitutive equation

The magnetic induction equation

The energy equationwhere, the dissipation function is given byand the current density vector due to the magnetic field and electric field is

The magnetic body force becomes

In the case of , the above equations reduce to

In this paper, we shall study the problem for zero pressure gradient in the x-direction namely . Also, the Eq. (14) that gives the pressure gradient in the y-direction will not be discussed, because it is clear that if are known the variation of the pressure p is easily determined.

The flow, induced magnetic and temperature field of the fluid in two regions are determined by the equations

Along with Eqs. (18)–(21), the following initial, boundary and interface conditions are considered:

In the present problem, it has assumed that the two plates are maintained at constant temperatures. Under the flow's assumptions, the temperature fields can be considered independent of the longitudinal variable x. The term involving in the energy Eq. (21) drops out for such conditions.

We introduce the non-dimensional entities

Equations (18)–(21) become (the notation “” is neglected)

The non-dimensional boundary conditions are

Velocity field distribution

To solve Eqs. (33) and (34) along with initial-boundary conditions (37) and (38), we use the Laplace transform method. Applying the Laplace transform to Eqs. (33) and (34), we obtainwhere, denotes the Laplace transform of function .

Eliminating between Eqs. (41) and (42), we obtain the following equation for the Laplace transform of the velocity :

Functions have to satisfy the boundary and interface conditions

The general solution of Eq. (43) iswhere

Using (44), the following expressions of the functions are found:

Although the inverse Laplace transform of the function (45) can be obtained using the residue theorem, the analytical expression of the is complicated and difficult to be used for numerical evaluation. For this reason, we prefer to obtain the numerical values of the inverse Laplace transform of using Stehfest’s numerical algorithm[25,26]

According to the Stehfest’s algorithm, the approximate values of function are given bywhere

In the above relation is an even, positive integer number and denotes the integer part of the real number .

The induced magnetic field

The magnetic field induction is given by Eq. (35) along with conditions (39). Applying the Laplace transform to Eq. (35), using initial conditions (37)2 and the velocity expression (45), we obtain for the Laplace transform of the differential equationalong with the boundary and interface conditions

The general solution of Eq. (50) is given bywhere

We introduce the notationswe have

Using the boundary and interface conditions (51), we obtain the algebraic system for the unknown functions whose solution iswherewith

Numerical values of the induced magnetic fields (52) shall be determined with the numerical approximate formula (48) updated for .

The temperature field

Approximate analytical solution

Let’s introduce notationwhere is given by Eq. (45).

Temperature fields are solutions of the equationsalong with the initial-boundary and interface conditions

Applying the Laplace transform to Eq. (65), we obtain the transformed equation

The general solution of Eq. (70) iswhere , have to satisfy conditions

From Eq. (72), we haverespectively,

Replacing (74) into (71), we obtainwhere,

Unknown functions , shall be determined using the boundary and interface conditions

After direct calculations, we obtainwhere,

Numerical values of the temperature could be determined with the numerical approximate formula (48) updated for . However, the expressions of the temperature fields (75) are quite complicated and also contain functions that are defined by means of some definite integrals. For this reason, the numerical inversion with the Stehfest's algorithm becomes difficult. Therefore, in the following, we will elaborate a numerical scheme adequate to the integration of the differential Eq. (65) with the initial-boundary conditions (66)–(69).

Numerical scheme

Temperature fields are solutions of the equationsalong with the initial, boundary and interface conditions

We consider the step size for the time t, step size for the spatial coordinate , and the sequences . Using the approximate formulas for the derivatives,

Equation (82)1 is written as

where,

The initial and boundary conditions (83) and (84) lead to the following relations:

For the interval , we consider the sequence

From the interface condition (85)1, we obtainwhile, from the condition (85)2, we have

Equation (82)2 is written aswhere,

Using the boundary condition (84)2, we obtain

Now, the above relations, determine the values of the temperature fields.

Indeed, using Eqs. (93) and (95), the values are determined.

The values of the temperature field in zone are determined using Eq. (87) together with Eq. (92).

Numerical results and discussion

Flows inside channels, such as magneto-hydrodynamic, electro-hydrodynamic, and electrokinetic flows are important in many microscale applications. The generation of such flows is often done with the help of external fields (mechanical, thermal, magnetic, electric) applied to the fluids. In this section, the generalized Couette flow studied in previous paragraphs is numerically and graphically simulated in the case when the motion of the upper plate is given by the function .

For the non-dimensional parameters that characterize the studied fluids and fields, the following numerical values were used: .

Figure 2 shows the profiles of fluid velocities in both flow regions for different values of the time t. As expected, the fluid flows more slowly in the region bounded by the lower wall. This behavior is due to the fact that the lower wall of the canal is fixed, so the viscosity force will slow down the fluid motion. In the region of the canal bounded by the upper wall, the fluid moves at a higher velocity due to the movement of the wall with the velocity given by the function . It is observed in Fig. 2 that for large values of the time t, fluid velocity on the upper wall becomes almost constant, because .

Figure 2

The profiles of velocities versus y for different values of the time t.

This property is most clearly evident from Fig. 3 which shows the profiles of velocities as a function of time t for three positions of the channel, namely, .The curves in Fig. 2, show that for fluid velocities tend to become constant. Also, Figs. 2 and 3 highlight that fluid velocity satisfy the interface condition .

Figure 3

Time-evolution of velocities in three different positions of channel.

Figures 4 and 5 show profiles of the ratio of induced magnetic field, and the external applied magnetic field, versus spatial coordinate y, respectively versus time t. Note that in the position of the channel, where the velocities are equal, the variation of the induced magnetic field is slower. It can be seen that in this position the induced magnetic field has a slow increase and after the value , the values of the induced magnetic field are almost constant. In the first flow zone, in which the velocity gradient has large variations, the induced magnetic field also has significant variations. These variations are attenuated in the second flow area where the speed gradient variations are smaller. Also, we note that for , the values of induced magnetic field become almost constant for any position of the channel. This behavior is in accordance with the fluid motion whose velocity tends to constant values for large values of the time t.

Figure 4

The profiles of induced magnetic fields versus y for different values of the time t.

Figure 5

Time-evolution of induced magnetic fields in three different positions of channel.

Figures 6 and 7 are plotted to show the variation in space and time of the dissipative functions . It is known that viscous dissipation acts as a source of energy in the fluid flow and it affects the temperature distribution. The shear stress within the fluid layer induced by the motion of the upper plate significantly influences the energy dissipation. It is observed in Figs. 6 and 7 that for a short time interval, the dissipative function in the region bounded by the lower wall is higher than the dissipative function in the other region. After this instant, the dissipative function in the region bounded by the upper plate becomes higher. This fact is generating by the evolution in time of the velocity in this region.

Figure 6

The profiles of dissipative functions versus y for different values of the time t.

Figure 7

Time-evolution of the dissipative functions in three positions of channel.

Conclusions

Unsteady, magneto-hydrodynamic generalized Couette flows of two immiscible fluids in a rectangular channel have been studied.

The flow is influenced by an inclined magnetic field and an axial time-oscillating electric field. Both fluids are considered electrically conducting and the solid boundaries are electrically insulated.

The problem is formulated into the dimensionless form and at the interface the velocity, shear stress, induced magnetic and temperature fields are considered continuous functions. Both walls are isothermal. The bottom wall of channel is fixed, while the upper wall is moving with a given time-dependent velocity. The nonslip conditions on the solid boundaries are also considered.

Approximate analytical solutions for the velocity, induced magnetic, and temperature fields have been determined using the Laplace transform method along with the numerical Stehfest's algorithm for the inversion of the Laplace transforms.

A numerical scheme based on the finite differences has been developed for the non-linear energy equation.

A particular case characterized by a time-exponential velocity of the upper wall, has been numerically and graphically studied to show the evolution of the fluid velocity, induced magnetic field, and viscous dissipation in both flow regions.