∂ ∂ ∂ ∂ ∂ ∂ ∂ + + = − + + − + − ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂

Because of the great importance of thermal instability in nature, in chemical processes, in separation processes, in industrial applications as well as in geophysical and astrophysical engineering, the effect of thermal diffusion on the combined MHD heat transfer in an unsteady flow past a continuously moving semi-infinite vertical porous plate which is subjected to constant heat has been investigated numerically under the action of strong applied magnetic field taking into account the induced magnetic field. This study is performed for cooling problem with lighter and heavier particles. Numerical solutions for the velocity field, induced magnetic field as well as temperature distribution are obtained for associated parameters using the explicit finite difference method. The obtained results are also discussed with the help of graphs to observe effects of various parameters on the above mentioned quantities.


Introduction
Magnetohydrodynamics (MHD) is currently undergoing a period of great enlargement and differentiation of subject matter.The science of magnetohydrodynamics was concerned with geophysical and astrophysical problems for a number of years.In recent years, the possible use of MHD is to affect a flow stream of an electrically conducting fluid for the purpose of thermal protection, braking, propulsion and control.From the point of applications, several investigators have made model studies on the effect of magnetic field on free convection flows.Some of them are Georgantopoulos [1], Nanousis et al. [2] and Raptis and Singh [3].Along with the effects of magnetic field, the effect of transpiration parameter, being an effective method of controlling the boundary layer has been considered by Kafoussias [4] and Singh [5].Singh et al. [6] studied magnetohydrodynamics heat and mass transfer flow of a viscous incompressible fluid past an infinite vertical porous plate under oscillatory suction velocity normal to the plate.The combined boundary layer heat and mass transfer of an electrically conducting fluid in MHD natural convection adjacent to a vertical surface was analyzed by Chen [7].The boundary layer behavior and mathematical configuration was previously discussed by Sakidas [8] [9].
All the above works are related to the stationary vertical surface.However, the flow past a continuously moving surface has many applications in manufacturing processes such as hot rolling, metal and plastic extrusion, continuous casting, and glass fiber and paper production.The heat transfer flow past a continuous moving plate with variable temperature was investigated by Soundalgekar and RamanaMurty [10].Al-Sanea [11] studied the steady laminar flow and heat transfer characteristics of a continuously moving vertical sheet of extruded material.An analytical study of the one dimensional steady combined heat and mass transfer by mixed convection flow of an incompressible electrically conducting viscous fluid past an electrically non-conducting continuously moving infinite vertical porous plate under the action of strong magnetic field with constant suction velocity, constant heat and mass fluxes was done by Chaudhary and Sharma [12].The level of concentration of foreign mass was assumed very low in this study so that the thermal and mass diffusion were neglected.They used perturbation technique to obtain the solution and did not show the complete results due to some analytical restrictions.Along with these studies, the effect of thermal diffusion on MHD free convection and mass transfer flows have also been considered by many investigators due to its important role particularly in isotope separation and in mixtures between gases with very light molecular weight (H 2 ,He,) and medium molecular weight (N 2 , air) (Eckert and Drake, [13]).Considering these aspects, model studies were carried out by many investigators of who the names of Kafoussias [14], Nanousis [15], Sattar and Alam [16] and Alam et al. [17] are worth mentioning.Recently both Soret and Dufour effects on steady mixed convection flow past a semi-infinite vertical porous flat plate in a porous medium with variable suction was investigated by Alam and Rahman [18].Quite recently Alam et al. [19] studied Numerical Study of the Combined Free-Forced Convection and Mass Transfer Flow Past a Vertical Porous Medium with Heat Generation and Thermal Diffusion.Very recently, a number of studies of unsteady heat transfer boundary layer flow were reported in the literature [20]- [26].However, the effect of induced magnetic is still not getting abundant attraction to the researchers.
We investigated numerically the transient heat transfer of viscous incompressible boundary layer fluid flow through a porous plate with induced magnetic field taking into account the strong magnetic field with constant heat fluxes.In this study, the viscous dissipation and joule heating terms in the energy equation have been considered for high speed flows.The governing equations of the problem contain a system of partial differential equations which are transformed by usual transformation into a non-dimensional system of partial coupled non-linear differential equations.The obtained non-similar partial differential equations solved numerically by finite difference method [27].The results of this study will be discussed for the different values of the wellknown parameters and shown graphically.

Mathematical Formulations
From the basis of studying Magneto Fluid Dynamics (MFD), for a boundary layer [8] [9] heat transfer by mixed convection flow of an electrically conducting viscous fluid that approximately grossly neutral with induced magnetic field, thermal diffusion, constant heat and mass fluxes, the generalized continuity equation, momentum equation, magnetic induction equation, energy equation and species equation together with the Ohm's law and Maxwell's equations are as furnished below, 1 Initially we consider that the plate as well as the fluid is at the same temperature ( ) everywhere in the fluid is same.Also it is assumed that the fluid and the plate is at rest initially in its own plane and instantaneously at time 0 t > , the temperature of the plate is raised to , which is there after maintained constant, where w T is the temperature at the wall and T ∞ is the temperature of the species far away from the plate.
The x-component of the momentum equation reduces to the boundary layer equation if the only contribution to the body force is made by gravity, the body force per unit volume, x F g = − , where g is the local acceleration due to gravity.The physical model of this study is provided in Figure 1.
There is no body force in the y-direction i.e.
which implies that ( ) Hence the x-component pressure gradient at any point in the boundary layer must equal to the pressure gradient in the quiescent region outside the boundary layer.However, in this region, 0 u v = = .Therefore the x-component pressure gradient of the momentum equation, P g x , where ρ ∞ is the density of the surrounding fluid at the temperature T ∞ .
( ) With the help of the Equation ( 7), the momentum Equations ( 2) and (3) become ( ) The magnetic Reynolds number of the flow is not taken to the small enough so that the induced magnetic field ( ) , , 0 x y H H is not negligible.The divergence equation of Maxwell's equation for the magnetic field gives, 0 For considering the direction of propagation is along the y-axis and H does not have any variation along the y-axis, which implies that 0 ( ) is the constant induced magnetic field.Hence the Equations ( 1) to ( 7) become ( ) By neglecting the small order terms from the Equations ( 10) to (15) we get Continuity equation Momentum equation ( ) Magnetic induction equation Energy equation and the corresponding initial and boundary conditions for the problem are everywhere 0, 0, 0, 0, where x, y are Cartesian coordinate system; u, v are x, y component of flow velocity respectively; g is the local acceleration due to gravity; β is the thermal expansion coefficient; υ is the kinematic viscosity; e µ is the magnetic permeability; ρ is the density of the fluid; 0 H is the constant induced magnetic field; x H be the x-component induced magnetic field; σ is the electrical conductivity; κ is thermal conductivity; p c is the specific heat at the constant pressure and w H is the induced magnetic field at the wall.Since the solutions of the governing Equations ( 16) to (19) under the initial (20) and boundary (21) conditions will be based on the finite difference method it is required to make the said equations dimensionless [27].After introducing the dimensionless quantities and proper simplifications we obtain the following nonlinear coupled partial differential equations in terms of dimensionless variables, Also the associated initial (20) and boundary (21) conditions becomes 0, 0, 0, 0, 0

Important Physical Parameters
( )

Numerical Solutions
Now we attempt to solve the governing second order nonlinear coupled dimensionless partial differential equations with the associated initial and boundary conditions.For solving a transient free convection flow with mass transfer past a semi-infinite plate, Carnahan et al.
[25] used the explicit finite difference method which is conditionally stable.On the contrary, the same problem was studied by Soundalgekar and Ramana Murty [10] by an implicit finite difference method which is unconditionally stable.The only difference between the two methods is that the implicit method being unconditionally stable is less expansive from the point of view of computer time.However, these two methods respectively employed by Carnahan et al.
From the concept of the above discussion, for simplicity the explicit finite difference method has been used to solve Equations ( 22) to (25) subject to the conditions given by ( 26) and (27).To obtain the difference equations the region of the flow is divided into a grid or mesh of lines parallel to X and Y axes where X-axis is taken along the plate and Y-axis is normal to the plate.
Let U ′ , V ′ , x H ′ and T ′ denote the values of U, V, x H and T at the end of a time-step respectively.Using the explicit finite difference approximation [27] we have, , Also the related boundary conditions are 0 0 0 0 , , , , 0, 0, 0, 0 Here the subscripts i and j designate the grid points with x and y coordinates respectively and the superscript n represents a value of time, n τ τ = ∆ where 0,1, 2, .n =  From the initial condition (4.5.5), the values of U, x H and T are known at 0 τ = .During any one time-step, the coefficients , i j U and , i j V appearing in Equations (29) to (31) are treated as constants.Then at the end of any time-step τ ∆ , the new temperature T ′ , the new velocity U ′ , the new induced magnetic field x H ′ and V ′ at all interior nodal points may be obtained by successive applications of Equations ( 31), ( 30), ( 29) and ( 28) respectively.This process is repeated in time and provided the time-step is sufficiently small, U, V, x H and T should eventually converge to values which approximate the steady-state solution of Equations ( 22) to (25).These converged solutions are shown graphically in Figures 2-13.

Results and Discussion
For the purpose of discussing the results of the problem, the approximate solutions are obtained for various parameters with small values of Eckert number.In order to analyze the physical situation of the model, we have computed the steady state numerical values of the non-dimensional velocity U, induced magnetic field             To get the steady state solutions, the computations have been carried out up to 80 τ = .We observed that the results of the computations, however, show no visible changes after 40 τ = .Thus the solution for 40 τ = are essentially steady state solutions.
Because of the special importance of cooling problem in nuclear engineering in connection with the cooling of reactors, the value of the Grashof number for heat transfer is taken to be positive ( ) The velocity profiles are shown on Figures 2-5.From these Figures, we see that, as the values of Grashof number ( r G ) increases the velocity increases.On the other hand, in case of the time τ increases the velocity in- crease up to τ = 40.But at τ = 50 the velocity remain steady state as the time being increases.
From the Figures 6-9 it is observed that the induced magnetic field decreases as the value of the Grashof number ( r G ) increases.Also induced magnetic field profile increases as time increases.But at τ = 50, the in- duced magnetic field also remain steady state as on the time goes.
Finally the influence of Grashof number ( r G ) on temperature profile is shown in Figures 10-13.It is interes- tingly observed that the temperature profile increases as the value of r G increases.In these case the tempera- ture profile for different values of r G initially increases and later on it is decreases as the value of r G in- creases.Beside this, for the time τ increases the temperature profile also increase up to τ = 40.But after that the temperature remains steady state as the time increases on the same manner.

Conclusions
An unsteady heat transfer flow through an electrically conducting incompressible viscous fluid past an electrically non-conducting continuously moving semi-infinite vertical plate under the action of strong magnetic field taking into account the induced magnetic field constant heat is investigated in this work.The resulting governing system of dimensionless coupled non-linear partial differential equations is numerically solved by an explicit finite difference method.The results are discussed for different values of important parameters as Magnetic Force Number, Magnetic diffusivity numbers, Grashof number, Prandtl number and Eckert number.Some of the important findings obtained from the graphical representation of the results are listed herewith;  The velocity & temperature profiles increases, when the value of r G increases. The magnetic induction decreases with the increases of r G .As the basis for many scientific and engineering applications for studying more complex vertical problems involving the flow of electrically conducting fluids, it is hoped that the present investigation of the study of applied physics of flow over a vertical surface can be utilized.In the migration of underground water or oil as well as in the filtration and water purification processes, the findings may be useful for study of movement of oil or gas and water through the reservoir of an oil or gas field.The results of the problem are also of great interest in geophysics and astrophysics in the study of interaction of the geomagnetic field with the fluid in geothermal region.

Figure 1 .
Figure 1.Physical configuration and coordinate system.
x H and temperature T within the boundary layer for different values of magnetic parameter (M), magnetic diffusivity ( ) m P , Grashof number ( ) r G , Prandtl number ( ) r P and Eckert number ( ) c E .

Figure 2 .
Figure 2. Velocity profiles for different values of r G where, 0.01 M = , 1.00 m P = ,

Figure 3 .
Figure 3. Velocity profiles for different values of r G where, 0.01 M = , 1.00 m P = ,

Figure 4 .
Figure 4. Velocity profiles for different values of r G where, 0.01 M = , 1.00 m P = ,

Figure 5 .
Figure 5. Velocity profiles for different values of r G where, 0.01 M = , 1.00 m P = ,

Figure 6 .
Figure 6.Induced magnetic fields for different values of r G where, 0.01 M = , 1.00 m P = , 0.71 r P = , and 0.01 c E = at time 10 τ = .

Figure 7 .
Figure 7. Induced magnetic fields for different values of r G where, 0.01 M = , 1.00 m P = , 0.71 r P = , and 0.01 c E = at time 30 τ = .

Figure 8 .
Figure 8. Induced magnetic fields for different values of r G where, 0.01 M = , 1.00 m P = , 0.71 r P = , and 0.01 c E = at time 50 τ = .

Figure 9 .
Figure 9. Induced magnetic fields for different values of r G where, 0.01 M = , 1.00 m P = , 0.71 r P = , and 0.01 c E = at time 60 τ = .

Figure 10 .
Figure 10.Temperature profiles for different values of r G where, 0.01 M = , 1.00 m P = , 0.71 r P = , and 0.01 c E = at time 10 τ = .

Figure 11 .
Figure 11.Temperature profiles for different values of r G where, 0.01 M = , 1.00 m P = , 0.71 r P = , and 0.01 c E = at time 30 τ = .

Figure 12 .
Figure 12.Temperature profiles for different values of r G where, 0.01 M = , 1.00 m P = , 0.71 r P = , and 0.01 c E = at time 50 τ = .

Figure 13 .
Figure 13.Temperature profiles for different values of r G where, 0.01 M = , 1.00 m P = , 0.71 r P = , and 0.01 c E = at time 60 τ = .
2.0.Since the most important fluids are atmospheric air, salt water and water so the results are limited to 0for water at 20˚C).However the values of another parameters M, m P and c E are chosen arbitrarily as M = 0.010, 05 and 0.09 (for Velocity profiles).Along with the obtained steady state solutions, the flow behaviors in case of cooling problem are discussed graphically.The profiles of the transient velocity, induced magnetic field and temperature versus Y are illustrated in Figures 2-13.