Numerical Solution of Fluid Flow in Horizontal Tube under Effects of Radiation Field ()
1. Introduction
The flow and heat transfer of electrically conducting fluid in channels and circular tubes under the effect of a transverse magnetic field occurs in magnetohydrodynamic (MHD) generators, pumps and flow meters and has applications in nuclear reactors, filtration, geothermal systems and others [1].
In 2009, Singha studied an analytical solution to the problem of magnetohydrodynamic (MHD) free convection flow of an electrically conducting fluid between two heated parallel plates in the presence of an induced magnetic field [2].
In 2010, Cooky C., and Omubo-Pepple V., investigated the combined effects of radiative heat transfer and transverse magnetic field on steady flow of an electrically conducting optically thin fluid through a horizontal channel filled with porous medium and non-uniform temperature at walls; they concluded that increasing magnetic field, radiation and porosity parameters reduced the velocity and temperature profiles, and the shear stress at the wall white increasing radiation parameters causes an increase in the magnitude of the rate of heat transfer [3].
In 2011, Dalia S., and Ioan P., have studied the entropy generation production for a problem of mixed convection flow of a fluid saturated porous medium through an inclined channel with uniform heat walls; he used the analytical results obtained for the velocity and temperature profiles to obtain the entropy production [4].
In 2012, Sudershan S., Reddy R., and Reddy J., have discussed the effect of radiation and chemical reaction on a free convection MHD flow through a porous medium boundary by vertical surface [5].
In 2015, Mukhopadhyay S. and Mandal I., studied the numerical solutions for steady MHD mixed convection boundary layer flow and heat transfer over a porous plate in the presence of the velocity and thermal slip; we found the effect of the magnetic field of non-compressible liquid on increasing fluid velocity, which in turn leads to a decrease in the temperature of the liquid [6].
In 2017, Kumari K. and Gayal M. presented the effect of mass transfer, viscous dissipative and suction parameter on two dimensional steady hydromagnetic viscous fluid flow between two parallel plates in the presence of thermal radiation and we found that velocity, temperature and concentration decrease with increase suction parameter and Reynolds number in addition to the relation among the different quality physicals [7]. The aim of this work, is to study the numerical solution of the equations from the heat transfer, diffusion and motion in a porous tube with the presences of magnetic field and radiation has been investigated. It’s found that the parameters Pr, Sc, R and Gr have a significant effect on the solution of the equation.
2. Formulation of the Problem
Consider the unsteady flow of a dissipative fluid passing through a long horizontal tube with porous walls, if we assume that the temperature at the two horizontal walls is constant so that
represents the lower (upper) walls temperature of horizontal tube with
. Let u, v be the velocity components in directions x and y respectively, we assume that all the components in z direction vanish. The magnetic field is applied perpendicular to the x-axis which in the direction of flow and apposite to the y direction as it is illustrated by Figure 1:
The governing continuity, momentum, energy and concentration equations are:
(2.1)
(2.2)
(2.3)
Derive Equation (2.2) with respect to y, and Equation (2.3) with respect to x, and subtract the resulting equations we get:
(2.4)
And this called the general motion equation.
(2.5)
By using the Roseland approximations consider the radiative heat flux for optically thick fluid is given by [8]:
(2.6)
where
is the Stefan-Boltzmann constant and
is the mean absorption coefficient. Assume that the difference in temperature within the flow is sufficiently small such as that
can be expressed as a linear function of the temperature, we expand
in a Taylor’s series about
and neglected higher order terms, thus [9].
(2.7)
Hence the equation of energy Equation (2.5) becomes:
(2.8)
(2.9)
where u, v are the velocity component in x and y directions respectively, t is the time and
are the temperature, concentration, density, kinematic viscosity, pressure, electrical conductivity, Boltzmann number, gravitational acceleration, permeability of medium, thermal conductivity, specific heat with constant pressure, radiation flux respectively.
The following boundary conditions of the Problem given by:
(2.10)
Let us introduce the following similarity transformation [3]
(2.11)
And non-dimensional parameters:
(2.12)
where M is Hartmann number,
is Thermal diffusion, R is Radiation coefficient, N is new physical quantity, Sc is Schmidt number, Da is Darcy number, Pr is Prandtl number, Gr is Gratshof number for heat transfer,
is Gratshof number for mass transfer,
is dispersion parameter. The above Equations (2.11) and (2.12) reduce the Equations (2.1), (2.4), (2.8) and (2.9) into the following system of non-dimensional equations:
(2.13)
(2.14)
(2.15)
(2.16)
The boundary conditions (2.10) in the non-dimensional form become:
(2.17)
3. Method of Solution
In order to solve the system of Equations (2.13)-(2.16) with the boundary conditions (2.17), we resort to ADI finite difference method [10], and to achieve this, we have to start with the last Equation (2.16) then Equation (2.15) and finally Equation (2.14) as following:
3.1. Solution of Diffusion Equation
(3.1.1)
(3.1.2)
Equations (3.1.1) and (3.1.2) can be reduced to give:
(3.1.3)
where
(3.1.4)
Followed by:
(3.15)
where
(3.1.6)
where
3.2. Solution of Heat Equation
(3.2.1)
(3.2.2)
Equations (3.2.1) and (3.2.2) can be reduced to give:
(3.2.3)
where
followed by
(3.2.5)
where
(3.2.6)
where
3.3. Solution of the General Motion Equation
(3.3.1)
(3.3.2)
Equations (3.3.1) and (3.3.2) can be reduced to give:
(3.3.3)
where
(3.3.4)
followed by
(3.3.5)
where
(3.3.6)
where
With boundary conditions:
(3.3.7)
The coefficients U, V are treated as constants during any one time-step of the computation, each of the equations (diffusion, heat, motion) creating a tridiagonal system which are solved by using Gauss elimination method, all are given in [5].
4. Conclusions
We present in this section some of the results obtained from the computation done on Equations (2.14), (2.15) and (2.16) for different points in the region of solution and this is because of the appearance of effect on this equations, and we find the following results (Figures 2-11):
1) With the increase of the radiation parameter, the faster steady state will be reached.
2) In energy equation, when
and
, and the time-step increases, the faster steady state will be reached.
3) With the decrease of the Schmidt number Sc, the faster steady state will be reached.
4) In Darcy equation, when
and
, and the time-step increases, the faster steady state will be reached.
5) With the decrease of the Gratshof number Gr, the faster steady state will be reached.
6) With the decrease of the Prandtl number Pr, the faster steady state will be reached.
Figure 3. The effect of Pr on the temperature.
Figure 4. The effect of R on the temperature.
Figure 5. The effect of Gr on the temperature.
Figure 6. The effect of the time-step on the temperature.
Figure 7. The effect of Sc on the Darcy.
Figure 8. The effect of Gr on the Darcy.
Figure 9. The effect of the time-step on Darcy.
Figure 11. The effect of the time-step on the motion.
7) In motion equation, when
, and the time-step increases, the faster steady state will be reached.
Acknowledgements
The authors are grateful to Assist. Prof. Dr. Ala’a Abdul-Raheem Ahmed Hammodat for his valuable remarks.