1. Introduction
Labropulu [1], elegantly stated that the movement of a fluid (liquid or gas) through a magnetic field may induce an electric field, and consequently a current. The current then interacts with the magnetic field to produce a force that acts on the flowing fluid. It has long been recognized that these interactions, and magneto-hydro-dynamic (MHD) flows in general, occur in nature and in man-made media and devices [1], and find applications in both natural settings, such as in plasma dynamics and in astrophysics [2], and in engineering applications, such as in power generation [3] [4].
Importance of these and many other applications have motivated various researchers to propose methods of solutions that are either specialized or approximate to provide satisfactory solutions to the inherently complex governing equations that lack complete analytical solutions. Many methods have been proposed and implemented, and have been reported in the literature (cf. [5] - [10] and the references therein).
Occurrence of MHD flows in the earth interior and in other porous materials underscores the need for studies of MHD flows in porous media. These flows enjoy a spectrum of engineering applications that have been investigated over the past few decades (cf. [4] and the references therein). Equations governing MHD flow through porous media have been thoroughly and elegantly analyzed by Geindreau and Auriault [4], who also presented a literature review of both the theoretical work, and the experiments that have been conducted in field.
Of particular interest to the current work is plane transverse MHD flow through a porous structure, wherein the magnetic vector field is considered normal to the plane of the flow, shown in Figure 1.
This type of flow might be of interest in some applications including lubrication theory where the magnetic field might have an impact on the load carrying capacity in mechanisms that include porous lining [11] [12] [13]. In the analysis to follow, we consider the steady, incompressible MHD flow through a porous structure with constant permeability. We parallel closely the elegant analysis and approach provided by Labropulu [1], in her study of Riabouchinsky flows of ordinary and second grade MHD fluids in free space in that we state, define and formulate the problem in a similar manner and assume that streamfunction is linear in one of the space variables. However, we present a method of solution that is based on the inverse approach introduced in by Hamdan [14] and received success in the solution of flow problems in porous media [15]. We provide expressions for the streamfunction, velocity components, vorticity, pressure functions and magnetic field.
While we consider in this work fluid flow in unconfined domains, solutions might be of interest in lubrication theory with porous lining under the influence of an electromagnetic field.
2. Governing Equations
Steady, magneto-hydro-dynamic (MHD) flow of an incompressible fluid, through a porous sediment, is governed by the following equations:
(1)
(2)
Figure 1. Representative sketch of flow direction.
(3)
(4)
wherein
is the magnetic permeability,
is the velocity vector field,
is the magnetic field, p is the pressure function,
is the fluid density,
is the base fluid viscosity coefficient,
is the effective viscosity of the fluid in the porous medium, k is the (constant) permeability,
is the electrical conductivity, and
is the gradient operator. We point out that the electric current density is given by
.
Employing the identity
(5)
Equation (3) takes the form
(6)
Now, governing Equations (1), (4), and (6) represent a determinate system of seven scalar equations in the seven unknown functions
,
and p, while Equation (2) represents an additional condition on
that signifies the absence of magnetic poles in the flow field.
Assuming the flow is in two space dimensions, (x, y), then the velocity vector field is given by
(7)
Assuming that that the magnetic field acts in a direction perpendicular to the plane of the flow, we take
(8)
Equations (2) and (8) imply that
and
, where subscript notation denotes partial differentiation.
Using (7) and (8) in (1), (4) and (6), we obtain the following governing equations in components’ form:
(9)
(10)
(11)
(12)
where
is the Laplacian in two dimensions, and
(13)
Equations (9), (10), (11), (12) and (13) represent five scalar equations in the five unknowns
and
(the magnetic pressure function). These equations govern the two-dimensional flow of a finitely conducting fluid (
is finite) in a constant permeability porous material in the presence of a magnetic field.
If the fluid is infinitely conducting then
,
and
. Diffusion Equation (12) thus reduces to:
(14)
3. Vorticity-Streamfunction and Generalized Pressure Formulation
For three dimensional flow, the vorticity vector is the curl of the velocity vector. This definition reduces to the following expression for vorticity,
, of a fluid in two space dimensions
(15)
Defining the generalized pressure function
(16)
then Equations (10) and (11) are replaced, respectively, by
(17)
(18)
Now, equation of continuity (9) implies the existence of a streamfunction
such that
(19)
and
(20)
With the introduction of the streamfunction, the equation of continuity (9) is automatically satisfied, and the vorticity definition (15) becomes
(21)
Using (19) and (20) in (17), (18), (12) and (14), we obtain, respectively
(22)
(23)
(24)
(25)
Once Equations (21), (22), (23) and (24), or (25), are solved for
and H as functions of x and y, velocity components can be computed from (19) and (20), then
can be obtained from (16) and finally p from (13).
4. Compatibility Equation and the Form of Streamfunction
The generalized pressure function
satisfying (22) and (23) is assumed to be a continuous function with the integrability condition,
. Integrability condition leads us to the following compatibility equation, derived by differentiating (22) with respect to y and (23) with respect to x, and equating the resulting derivatives:
(26)
or
(27)
where
.
The streamfunction
is a solution to (26) or (27). Its determination is a formidable task at present. We therefore assume that the streamfunction is linear in one of the space dimensions. We thus let
(28)
Using (28) in (26), and equating coefficients of similar powers of y, we obtain the following ordinary differential equations for
and
, wherein prime notation denotes ordinary differentiation
(29)
(30)
Integrating (29) and (30) once with respect to x, we obtain, respectively:
(31)
(32)
where
and
are arbitrary constants, whose values will be discussed following solution of pressure function, below. Once Equations (31) and (32) are solved for
and
, the streamfunction given by (28) can then be obtained.
Using (28) in (19) and (20), we obtain the following expressions for the velocity components:
(33)
and
(34)
Equations (21) and (28) give the following expression for vorticity:
(35)
Using (28) in (22) and (23) we obtain, respectively;
(36)
(37)
From (31) and (32), we can write
(38)
(39)
Using (38) and (39) in (36) and (37), followed by integrating (36) and (37), we obtain the following expression for the generalized pressure function
(40)
where
is an arbitrary constant.
Using (28) in the diffusion Equation (24), we obtain
(41)
If the fluid is infinitely conducting, then (40) reduces to
(42)
Using (33), (34) and (40) in (16), we obtain
(43)
while using (43) in (13) gives
(44)
5. Determination of the Streamfunction
In the absence of a general solution to Equations (31) and (32), the functions
and
have customarily be determined in accordance with the following procedure that received attention in the study of Navier-Stokes flow and in MHD flow [1]:
1) Find a particular solution for
satisfying (31).
2) Substitute the solution for
, found in step 1), in (32) and find a general solution for
.
3) Once
and
are found, the streamfunction
can be calculated using (28).
This procedure has some drawbacks that have been discussed in details in [14] and [15], and include the difficulty of finding particular solutions to (31) for general flows, and the presence of many arbitrary constants in the function
. These and other drawbacks can be overcome with the following modified, inverse procedure.
1) Assume the form of the function
.
2) Substitute the assumed form of
in (32). This produces an ordinary differential equation for
.
3) Solve the equation that arises in step 2) for
. Solution obtained in this step will involve some arbitrary constants and the solution must satisfy (31).
4) Substitute
that is obtained in step 3) in (31) and determine the arbitrary constants.
We apply the modified algorithm to obtain a solution to Equations (31) and (32). In order to parallel solutions obtained for flow through porous media without magnetic effects, discussed in [15], we assume the following form of
:
(45)
wherein
and
are known parameters.
Using (45) in (31) and (32) we obtain the following first order ordinary differential equation for
:
(46)
whose solution takes the form
(47)
where
is an arbitrary constant.
Using (45) and (47) in (28), we obtain the following expression for the streamfunction:
(48)
Using (48), velocity components u and v, given by Equations (33) and (34), respectively, and vorticity, given by Equation (35), take the following forms, respectively
(49)
(50)
(51)
The generalized pressure function
given by Equation (40), and the pressures
and p, given by Equations (43) and (44), respectively, take the following forms, respectively
(52)
where
is an arbitrary constant.
(53)
(54)
Equation (54) requires the determination of
in order to determine the pressure. To accomplish this, we consider the cases of infinitely conducting fluid and finitely conducting fluid, discussed in what follows.
If the fluid is infinitely conducting then the diffusion equation is given by Equation (25), or, equivalently, (42). In this case if
or if
, where
is an arbitrary function of
, then (25) is satisfied. For example, the following function satisfies (42):
(55)
For finitely conducting fluid, we follow the procedure introduced by Labropulu [1] which leads us to a solution to the diffusion Equation (41) of the form
(56)
where
(57)
and the functions
and
in (56) must be such that
(58)
or equivalently
(59)
Using (45) and (47), we write (59) as
(60)
where
(61)
and B is an arbitrary constant.
For the sake of illustration, if we take
in (56), then
, where
satisfies
(62)
whose solution is given by
(63)
where
and
are arbitrary constant.
6. Determination of the Arbitrary Constants
In the above analysis, we derived expressions for the flow quantities under the assumption that the streamfunction is linear in one of the space variables. Solutions obtained are given by Equations (48) to (54), and (55) or (56), and involve seven arbitrary constants,
to
.
The expression for
given by Equation (47), namely
(47)
contains two arbitrary constants,
and
. This solution for
must satisfy (29). We can use (29) as a condition to help us determine the arbitrary constants
and
.
Upon substituting (47) into (29), and simplifying, we obtain the following condition on
and
(64)
Condition (64) implies that either
or
.
Case 1: If
and
When
then the pressure is unbounded, unless the fluid extent is finite. This leads to the following observation.
Observation (1). In the MHD flow at hand, if the fluid is of finite extent then the flow variables are given by
(65)
(66)
(67)
(68)
(69)
(70)
(71)
where
(72)
(73)
and
and
satisfy
(74)
If in the above solutions,
, we recover solution to flow through a free-space domain, that is a domain without the presence of a porous matrix.
Case 2: If
For fluid of infinite extend, we take
so that the pressure is bounded. Taking
as well leads to the following observation.
Observation (2). In the MHD flow at hand, if the fluid is of infinite extent and
, then the flow variables are given by
(75)
(76)
(77)
(78)
(79)
(80)
(81)
where
(82)
(83)
and
and
satisfy
(84)
If in the above solutions,
, we recover solution to flow through a free-space domain, that is a domain without the presence of a porous matrix.
Case 3: If
and
For fluid of infinite extent, coefficients of y and y2 must vanish in order for the pressure function to be bounded. We therefore take
and state the following observation.
Observation (3). In the MHD flow at hand, if the fluid is of infinite extent and
, then the flow variables are given by
(85)
(86)
(87)
(88)
(89)
(90)
(91)
(92)
(93)
and B is an arbitrary constant, and the functions
and
satisfy
(94)
If in the above solutions,
, we recover solution to flow through a free-space domain, that is a domain without the presence of a porous matrix.
In all of the cases discussed above, all remaining arbitrary constants must be determined using physical conditions on the flow variables.
7. Potential Flow
Using the presented method of solution, if we choose
as a linear function in the variable x, namely
(95)
where C and D are known parameters, Equations (31) and (32) yield
(96)
and Equation (28) gives the following expression for the streamfunction:
(97)
Velocity components take the form:
(98)
(99)
and the vorticity vanishes since
(100)
If in the above solutions,
, we recover solution to flow through a free-space domain, that is a domain without the presence of a porous matrix.
8. Conclusion
In this work, we considered plane transverse MHD flow through a porous structure of contant permeability. In order to solve the governing equations, we assumed the streamfunction to be linear in one of the space variables, namely
. We derived a compatibility equation that resulted in two coupled ordinary differential equations the solution to which comprises solution for the streamfunction of the flow. In order to avoid the presence of a large number of arbitrary constants, we implemented an inverse approach to obtain the solution. Advantage to this implemented approach is avoidance of Riabouchinsky assumption that
must be taken as zero. We have obtained various forms of the flow variables depending on three observations made about some of the arbitrary constants. Solutions obtained are for both cases of finitely conducting and infinitely conducting fluids, and produce same forms of solutions obtained by Labropulu [1], when permeability goes to infinity. It is also worth noting that solutions for the streamfunction, velocity and viscosity are dependent on permeability of the medium and are the same results obtained by Kamel and Hamdan [15].
Acknowledgements
The authors wish to express their sincere thanks to the referees of this work who constructive suggestions lead to enhancing this work.
Nomenclature
: arbitrary constants
: functions of x
: functions of y
R: arbitrary function of the streamfunction
: arbitrary parameters
: base fluid viscosity coefficient
: effective viscosity
: magnetic permeability
: space variables
: velocity vector field
u: tangential velocity component
v: normal velocity component
h: generalized pressure function
: magnetic field
H: magnetic field third component
p: pressure function
: magnetic pressure function
: fluid density
k: permeability
: electrical conductivity
: vorticity vector third component
: streamfunction
: gradient operator
: Laplacian operator
: curl operator