This paper presents the results of exact solutions and numerical simulations of strongly-conductive and weakly-conductive magnetic fluid flows. The equations of magnetohydrodynamic (MHD) flows with different conductivity coefficients, which are independent of viscosity of fluids, are investigated in a horizontal rectangular channel under a magnetic field. The exact solutions are derived and the contours of exact solutions of the flow for magnetic induction modes are compared with numerical solutions. Also, two classes of variational functions on the flow and magnetic induction are discussed for different conductivity coefficients through the derived numerical solutions. The known results of the phenomenology of magnetohydrodynamics in a square channel with two perfectly conducting Hartmann-walls are just special cases of our results of magnetic fluid.
The first classical study of electro-magnetic channel flow was carried out by Hartmann in the 1930s [
In magnetic fluids, the fluid dynamic phenolmena with magnetic induction create new difficulties for the solution of the problems under consideration. The classical Hartmann flow can be further generalized to include arbitrary electric energy extraction from or addition to the flow. In general, classical MHD flows are dealt with using the exact solution of the Couette flow which is presented when the magnetic Prandtl number is unity [
The exact solutions of appropriately simplified physical problems provide estimates for the approximate solutions of complex problems. In view of its physical importance, the flow in a channel with a considerable length, rectangular, two-dimensional, and unidirectional cross section, which is assumed steady, pressure-driven of an incompressible Newtonian liquid, is the simplest case to be considered. In such a flow, taking into account the symmetrical planes and and an exact solution is obtained by using the separation of variables. The solution indicates that, when the width-to-height ratio increases, the velocity contours become flatter away from the two vertical walls and that the flow away from the two walls is approximately one-dimensional (the dependence of on is weak) [
If all walls are electrically insulating, , Shercliff (1953) has investigated principle sketch of the phenomenology of Magnetohydrodynamics (MHD) channel flow of rectangular cross-section with Hartmann walls and side walls [
Other flow configurations in basic MHD may include Hele-Shaw cells. Wen et al. [6,7] were motivated to visualize the macroscopic magnetic flow fields in a square Hele-Shaw cell with shadow graphs for the first time, taking advantage of its small thickness and corresponding short optical depth. Examples of applications of MHD include the chemical distillatory processes, design of heat exchangers, channel type solar energy collectors and thermo-protection systems. Hence, the effects of combined magnetic forces due to the variations of magnetic fields on the laminar flow in horizontal rectangular channels are important in practice [8–10].
In the present study, we consider the characteristics of magnetic fluids in a horizontal rectangular channel under the magnetic fields and use the flow equations with a conductivity coefficient. The exact solutions of the strongly-conductive and weakly-conductive magnetic fluids are considered using the series expansion technique in order to obtain the relationship between the flow and magnetic induction. Also, a quadratic function on flow and magnetic induction is studied to verify the characteristic of flow field using the obtained solutions.
The configuration of the flow geometry is illustrated in
The Maxwell’s equations in their usual form
, , (3)
with the relation equations and the Ohm’s law given as
where is the permeability of free space, is the magnetic susceptibility (H/m), is the conductivity, is the magnetic induction, is the magnetic field (A/m), and is the magnetization (A/m).
We choose the axis such that the velocity vector of the fluid is and from the continuity Equation (1), we have. We also choose =, where is a constant representing magnetic induction.
Applying the Maxwell’s equation and, we have. To simplify our presentations, the following assumptions are made for related variables:
, (6)
, (7)
, (9)
And, Equation (3) is satisfied. As there is no excess charge in the fluid, then, by using (5), is obtained as follows
The magnetic fluid boundary conditions considered here are
at,
at,
We shall also assume that all quantities are independent of time, that is to say, the fluid we consider here is in a steady state.
The magnetic fluid is called strongly-conductive if the term appears [
(12)
where. Note that Sutton and Sherman [
For Hartmann flow, it is feasible to replace by for a simple model, then =
where.
The axial pressure gradient is taken to be if the gravitational field is neglected, and, where is the viscosity of magnetic fluid. Combining Equations (12)–(14) yields
Let
, (17)
then Equations (15) and (16) are reduced to
where the Hartmann number is defined as.
The solution for is obtained by expressing over the range as a cosine Fourier series,
where is a constant. The solution for is then written
(21)
where and are given as
It should be pointed out that the solution for is just the same function as in which and are displaced by and, respectively, which are given by
(23)
If the fluid is weakly-conductive and the field is not time-dependent, the term will disappear as shown in equations (104) [
Replacing with, we have = as well as in Subsection 2.1. Then, as the axial pressure gradient is taken to be, where
where is a constant. Let
, (28)
Thus, the solution of Equation (28) for as in (17) is also obtained
where and are as follows
Note that the solution for is the same function as in which and are given as follows
For a weak-conductive fluid, its solution is simply the solution of the conductive fluid with a conductive coefficient (for). We find that which is obviously independent of the fluid viscosity.
Here we only consider a unidirectional two dimensional flow without a magnetic field, so that, Equations (18) and (19) are reduced to
The boundary conditions are as follows
at (33)
at (34)
at,(35)
In order to obtain an exact solution of Equation (32), comparing with our above results, we have
where is a constant. The problem consisting Equation (32) and its conditions are solved similarly using the separation of variables, which has the solution as follows [
where
, (38)
In the present study, the flow fields and their associated functions are presented in the flow region with,. Since ranges from 10 to 100 in most practical problems, the initial magnetic induction is taken to be (kg.s-2.A-1), (H/m), and the constant.
In Figures 4(a) and 4(b), the development of the velocity profile in and directions are shown for various values of the conductivity coefficient. For the symmetry, we only consider two cases: (a),; and (b),. Several interesting observations are readily made from the results. The cooperate process of and is shown in the above analysis. In order to clearly show the self-governed process of and, the contours of the velocity versus coordinate, and the velocity versus coordinate are given. It is clear that the velocity gradients increase quickly near the boundary walls and as is increased. On the other hand, the exact solutions are multiple hyper-cosine functions of, and cosine functions of. Therefore, the velocity gra-
dient is larger near the boundary walls than that near the boundary walls for a given.
For a given Hartmann number, comparing
change with different. Differently, in MHD, Shercliff and Hunt considered the induction equation and used only the linear constitutive equation with Hartmann number. They gave velocity profile and current paths for different Hartmann number.
For many years the variation techniques have been effectively applied to problems in the theory of elasticity. However, they are rarely used in fluid dynamic problems. The great utility in elasticity problems are due to the fact that they can be conveniently applied to linear problems. This, of course, explains why they are not frequently used in fluid dynamics since most such problems are nonlinear [10,11].
For the conductive fluid problems of the type being considered here, we recall the governing Equations (15) and (16) which are linear for and and the variation technique may be tried. Firstly, consider the following integral [
where is some given function of, ,. Clearly, the value of the integral depends on the choice of the functions, and. Now, let us pose the following problem: to obtain functions, and to minimize the value of. As is well known from variational calculus, the necessary conditions that, and for minimized are the Euler equations:
Simply, let the parameters be fixed at, then. According to and in the Equation (17), we only considered the function of, and in order to minimize the special function as follows
where and are considered as the function of, and, , , and.
Let, , then. The expressions of and are called the velocity decompositions of the magnetic induction and the flow field with a variable coefficient of the flow and the magnetic induction, where. From
It is also easy to know that has the same variation characteristic as the following function
(44)
The differential of on is given by
(45)
It is obvious that is a nonlinear function of, and is very complex to study the variation characteristics of by using the method of mathematical analysis.
With Equation (45), the variation characteristics of the function is determined by flow and magnetic induction as a function of.
Based on the above analysis of Equation (17), let, with, we call the velocity of magnetic fluid flow, which is equivalence to the velocity of magnetic fluid flow evoked by magnetic force. A total energy function is defined by
=+ (46)
where is the kinetic energy, is the magnetic energy.
By the calculus of variations, we have
=+ (47)
As and are nonlinear functions of, it is very complex to study the variational characteristics of by using the method of mathematical analysis.
As seen in Figures 2(a) and 2(c), in a region of (0,0), the contour of the function is more similar to that of the flow for. Furthermore, the gradient of flow which is larger than that of magnetic induction, the distribution value of is largely affected depending on the flow. On the contrary, for, the gradient of magnetic induction is larger than that of flow as is largely affected depending upon the gradient of the magnetic induction.
It is observed that the difference of flow and magnetic induction are almost the same for in a region of (0,0), where the distribution of the function is determined by the gradients of both the magnetic induction and the flow in this region. Furthermore, near and for any, the difference of flow is acute singularly and the function is also changed singularly. It is noted that has only one limit point for, and F has two limit points for.
As seen in
1) For magnetic fluid of this work and Magnetohydrodynamics (MHD) of Reference [
tions are different, the Hartmann numbers are also different. For magnetic fluid, conductivity coefficient is an important coefficient to analyze the flow and the current. For MHD, Hartmann number is the controlling coefficient, Shercliff and Hunt studied velocity profile and current paths for different Hartmann number [3–5].
2) For conductivity coefficient, the velocity contours for steady unidirectional flow is shown in a rectangular channel with a cosine Fourier series function of. Our result is in agreement with published findings.
3) A velocity decomposition and composition function, and a total energy variational function, on the flow and magnetic induction are considered. The variational characteristics of are analyzed only using the characteristics of the resultant flow field and the magnetic induction, and the number of its limit points changes as changes. It is shown in numerical simulations that the gradient of total energy is affected by the kinetic energy and the magnetic energy as changes.
4) Theoretically, the strongly-conductive and weaklyconductive magnetic fluid flows are studied on different conductivity coefficients which are independent of fluid viscosity in a horizontal rectangular channel.