Plane Transverse MHD Flow through Porous Media

Plane, transverse MHD flow through a porous structure is considered in this work. Solution to the governing equations is obtained using an inverse method in which the streamfunction of the flow is considered linear in one of the space variables. Expressions for flow quantities are obtained for finitely conducting and infinitely conducting fluids.


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 pro-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.

Governing Equations
Steady, magneto-hydro-dynamic (MHD) flow of an incompressible fluid, through a porous sediment, is governed by the following equations: wherein m µ is the magnetic permeability, v is the velocity vector field, H is the magnetic field, p is the pressure function, ρ is the fluid density, µ is the base fluid viscosity coefficient, e µ 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 ∇ × H .
Employing the identity Now, governing Equations (1), (4), and (6) represent a determinate system of seven scalar equations in the seven unknown functions v , H and p, while Equation (2) represents an additional condition on H 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 ( ) 0, 0, H = H (8) Equations (2) and (8)  Using (7) and (8) in (1), (4) and (6), we obtain the following governing equations in components' form:

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 x y v u = − ω (15) Defining the generalized pressure function ( ) then Equations (10) and (11) are replaced, respectively, by Now, equation of continuity (9) implies the existence of a streamfunction and With the introduction of the streamfunction, the equation of continuity (9) is automatically satisfied, and the vorticity definition (15) becomes Using (19) and (20) in (17), (18), (12) and (14), we obtain, respectively

Compatibility Equation and the Form of Streamfunction
The generalized pressure function ( ) , h x y satisfying (22) and (23) 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 Using (28) in (26), and equating coefficients of similar powers of y, we obtain the following ordinary differential equations for ( ) f x and ( ) g x , wherein prime notation denotes ordinary differentiation Integrating (29) and (30) once with respect to x, we obtain, respectively: where 1 C and 2 C are arbitrary constants, whose values will be discussed following solution of pressure function, below. Once Equations (31) and (32) are solved for ( ) f x and ( ) g x , the streamfunction given by (28) can then be obtained.

Determination of the Streamfunction
In the absence of a general solution to Equations 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 ( ) g x . These and other drawbacks can be overcome with the following modified, inverse procedure. 1) Assume the form of the function ( ) g x .
2) Substitute the assumed form of ( ) g x in (32). This produces an ordinary differential equation for 3) Solve the equation that arises in step 2) for ( ) f x . Solution obtained in this step will involve some arbitrary constants and the solution must satisfy (31).

4) Substitute
( ) f x 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 ( ) wherein , α β and γ are known parameters.
Using (45) in (31) and (32) we obtain the following first order ordinary diffe- whose solution takes the form ( ) where 4 C is an arbitrary constant.
Using (45) and (47) in (28), we obtain the following expression for the streamfunction: Using (48), velocity components u and v, given by Equations (33) and (34), respectively, and vorticity, given by Equation (35), take the following forms, re- The generalized pressure function ( ) x e e y h C C y C y C k where 5 C is an arbitrary constant.
and B is an arbitrary constant.
For the sake of illustration, if we take whose solution is given by where 6 C and 7 C are arbitrary constant.

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. Solu- e contains two arbitrary constants, 2 C and 4 C . This solution for ( ) f x must satisfy (29). We can use (29) as a condition to help us determine the arbitrary constants 2 C and 4 C .
and ( ) F y and ( ) If in the above solutions, k → ∞ , we recover solution to flow through a free-space domain, that is a domain without the presence of a porous matrix. x y y C y k  In all of the cases discussed above, all remaining arbitrary constants must be determined using physical conditions on the flow variables.

Potential Flow
Using the presented method of solution, if we choose ( ) g x as a linear function in the variable x, namely ( ) where C and D are known parameters, Equations (31) and (32) yield ( ) and Equation (28) gives the following expression for the streamfunction: Velocity components take the form: If in the above solutions, k → ∞ , we recover solution to flow through a free-space domain, that is a domain without the presence of a porous matrix.

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 ( ) ( ) yf x g x = + ψ . 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 ( ) g x 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