Hall Effects on Unsteady Free Convection in a Heated Vertical Channel in Presence of Heat Generation

Transcription

1 Applied athematics and Physics, 013, Vol. 1, No. 3, Available online at Science and Education Publishing DOI: /amp-1-3- Hall Effects on Unsteady Free Convection in a Heated Vertical Channel in Presence of Heat Generation Sanatan Das 1,*, Rabindra Nath Jana 1 Department of athematics, University of Gour Banga, alda, India Department of Applied athematics, Vidyasagar University, idnapore, India *Corresponding author: jana61171@yahoo.co.in Received September 11, 013; Revised September 4, 013; Accepted September 6, 013 Abstract Hall effects on an unsteady free convective flow of a viscous incompressible electrically conducting fluid between two heated vertical plates in the presence of a transverse applied magnetic field and heat generation have been studied. The governing equations are solved analytically using the Laplace transform technique. In order to get the physical insight into the problem, the effects of important parameters on temperature and velocity profiles are shown graphically and discussed in detail. Keywords: Hall currents, free convection, heat generation, Prandtl number, Grashof number, frequency parameter Cite This Article: Sanatan Das, and Rabindra Nath Jana, Hall Effects on Unsteady Free Convection in a Heated Vertical Channel in Presence of Heat Generation. Applied athematics and Physics 1, no. 3 (013): doi: /amp Introduction Heat transfer has emerged as a central discipline in contemporary engineering science. There has been growing interest in buoyancy-induced flows and the associated heat over the past three decades, because of the importance of these flows in many different areas, such as cooling of electronic equipment, pollution, materials processing, energy systems and safety in thermal processes. Several books, for instance, Turner [1], Jaluria [], Kakac et al. [3], Gebhart et al. [4] and Bejan and Kraus [5] may be consulted for detailed discussions on this subject. The study of heat generation or absorption in moving fluids is important in problems dealing with chemical reactions and those concerned with dissociating fluids. Heat generation effects may alter the temperature distribution and this in turn can affect the particle deposition rate in nuclear reactors, electronic chips and semi conductor wafers. Although exact modeling of internal heat generation or absorption is quite difficult, some simple mathematical models can be used to express its general behaviour for most physical situations. Heat generation or absorption can be assumed to be constant, space-dependent or temperature-dependent. The use of magnetic field that influences heat generation/absorption process in electrically conducting fluid flows has many engineering applications. For example, many metallurgical processes which involve cooling of continuous strips or filaments, which are drawn through a quiescent fluid. The properties of the final product depend to a great extent on the rate of cooling. The rate of cooling and therefore, the desired properties of the end product can be controlled by the use of electrically conducting fluids and the applications of the magnetic fields (Vajravelu and Hadjinicalaou [6]). Some processes in which heat is generated or absorbed may proceed in equipment and in the surrounding medium, for example, heat generation in the fuel elements of nuclear reactors due to the deceleration of nuclear fission fragments and slowing down of neutron fluxes, during a number of chemical reactions, etc. In the ionized gases, the current is not proportional to the applied potential except when the field is very weak in an ionized gas where the density is low and the magnetic field is very strong, the conductivity normal to the magnetic field is reduced due to the free spiraling of electrons and ions about the magnetic lines of force before suffering collisions and a current is induced in a direction normal to both electric and magnetic fields. This phenomenon, well known in the literature, is called the Hall effect. The study of hydromagnetic flows with Hall currents has important engineering applications in problems of magnetohydrodynamic generators and of Hall accelerators as well as in flight magnetohydrodynamics. It is well known that a number of astronomical bodies posses fluid interiors and magnetic fields. It is also important in the solar physics involved in the sunspot development, the solar cycle and the structure of magnetic stars. The unsteady hydromagnetic flow of a viscous incompressible electrically conducting fluid through a vertical channel is of considerable interest in the technical field due to its frequent occurrence in industrial and technological applications. The Hall effects on the flow of ionized gas between parallel plates under transverse magnetic field have been studied by Sato [7]. iyatake and Fujii [8] have discussed the free convection flow between vertical plates - one plate isothermally heated and other thermally insulated. Natural convection flow between vertical parallel plates- one plate with a uniform heat flux and the

2 46 Applied athematics and Physics other thermally insulated has been investigated by Tanaka et al. [9]. Hall effects on the hydromagnetic convective flow through a channel with conducting walls have been investigated by Dutta and Jana [10]. Joshi [11] has studied the transient effects in natural convection cooling of vertical parallel plates. Singh [1] have described the natural convection in unsteady Couette motion. Singh et al. [13] have studied the unsteady free convective flow between two vertical parallel plates. The natural convection in unsteady HD Couette flow with heat and mass transfers has been analyzed by Jha [14]. Narahari et al. [15] have studied the unsteady free convective flow between long vertical parallel plates with constant heat flux at one boundary. The unsteady free convective flow in a vertical channel due to symmetric heating have been described by Jha et al. [16]. Singh and Paul [17] have studied the unsteady natural convective between two vertical walls heated/cooled asymmetrically. Guria and Jana [18] have discussed the Hall effects on the hydromagnetic convective flow through a rotating channel under general wall conditions. Jha and Ajibade [19] have studied the unsteady free convective Couette flow of heat generating/absorbing fluid. Rajput and Sahu [0] have studied the unsteady free convection HD flow between two long vertical parallel plates with constant temperature and variable mass diffusion. Das et al. [1] have studied the radiation effects on free convection HD Couette flow started exponentially with variable wall temperature in presence of heat generation. The effect of radiation on transient natural convection flow between two vertical walls have been described by andal et al. []. Das et al. [3] have studied the radiation effects on unsteady HD free convective Couette flow of heat generation/absorbing fluid. Sarkar et al. [4] have studied an oscillatory HD free convective flow between two vertical walls in a rotating system. Kalita [5] has obtained an exact solution of an unsteady free convection HD flow and heat transfer between two heated vertical plates with heat source. The object of the present paper is to study the Hall effects on the unsteady HD free convective flow of a viscous incompressible electrically conducting fluid in a vertical channel. The problem is solved analytically for the velocity, temperature, rate of heat transfer and shear stress. The effects of various emerging parameters on the velocity and temperature distributions are shown and discussed with the help of graphs.. Formulation of the Problem and Its Solution Consider the unsteady magnetohydrodynamic flow of a viscous incompressible electrically conducting fluid between two infinitely long vertical parallel plates with variable temperature separated by a distance h. The flow is set up by the buoyancy force arising from the temperature gradient. Choose a Cartesian co-ordinates system with the x - axis along the vertical plates at z = h, the z - axis perpendicular to the plates, y -axis is normal to the zx -plane and the origin of the axes at the middle point between the plates [See Figure 1]. Initially, at time t 0, the fluid is rest and the temperature of the fluid and the plates are same as T h. At time t >0, the temperature of the plates at z = h is raised to n t Th + ( T0 Th) 1 e, where n is the frequency of the temperature oscillations and T 0 being the reference temperature. A uniform transverse magnetic field B 0 is applied perpendicular to the plates. We assume that the flow is laminar and the pressure gradient term in the momentum equation can be neglected. It is assumed that the effect of viscous and Joule dissipations are negligible. As the plates are infinitely long, the velocity field and temperature distribution are functions of z and t only. Figure 1. Geometry of the problem Under the usual Boussinesq approximation, the flow is governed by the Navier-Stokes equations u u B0 = ν + gβ ( T T ), h + jy t z ρ ρ v v B0 = ν j, x t z ρ T T cp = k + Q ( T T ), h t z (1) () (3) where ρ is the fluid density, T the fluid temperature, T h is the plate temperature, ν the kinematic viscosity, u and v are fluid velocity components, β the coefficient of thermal expansion, g the acceleration due to gravity, k the thermal conductivity, c p the specific heat at constant pressure and Q a constant. The initial and boundary conditions for the velocity and temperature distributions are t 0 : u = v = 0, T = Th for h z h, (4) nt t > 0 : u = 0, v = 0, T = Th + ( T0 Th)(1 e ) at z = h. The generalized Ohm's law, on taking Hall currents into account and neglecting ion-slip and thermo-electric effect, is (see Cowling [6]) ωτ e e j + ( j B) = σ ( E+ q B), (5) B0

3 Applied athematics and Physics 47 where j is the current density vector, q the velocity vector, B the magnetic field vector, E the electric field vector, ω e the cyclotron frequency, σ the electrical conductivity of the fluid and τ e the collision time of electron. We assume that the magnetic Reynolds number for the flow is small so that the induced magnetic field can be neglected. This assumption is justified since the magnetic Reynolds number is generally very small for partially ionized gases (Shereliff [7]). The solenoidal relation B =0 for the magnetic field gives Bz = B 0 = constant everywhere in the fluid where B (0,0, B0 ). Further, if ( jx, jy, j z) be the components of the current density j, then the equation of the conservation of the current density j =0 gives j z = constant. This constant is zero since j z = 0 at the wall which is electrically nonconducting. Thus j z = 0 everywhere in the flow. Since the induced magnetic field is neglected, the axwell's B equation E = becomes E =0 which gives t E x E y =0 and =0. This implies that E x = constant z z and E y = constant everywhere in the flow. We choose this constants equal to zero, i.e. Ex = E y =0. In view of the above assumption, the equation (5) gives jx + m jy = σ vb0, (6) j mj ub y x = σ 0, (7) where m = ωτ e eis the Hall parameter. Solving (6) and (7) for j x and j y, we have σ B0 jx = 1 ( v + mu + m ), (8) σ B0 jy = ( mv u). (9) 1+ m On the use of (8) and (9), the momentum equations (1) and () along x - and y -directions become u u σ B0 = ν + gβ ( T T ) ( ), h u mv (10) t y ρ(1 + m ) v v σ B0 = ν ( v + mu). t y ρ(1 + m ) Introducing non-dimensional variables h z ν t T Th ( u1, v1) = ( uv, ), η =, τ =, θ =, ν h h T0 Th equations (10), (11) and (3) become u1 u1 = + Grθ ( u1 mv1 τ η 1+ m ), (11) (1) (13) v1 v = 1 ( v 1+ mu1), τ η 1+ m θ θ Pr = + QPrθ, τ η (14) (15) σ B0 h where = is the magnetic parameter, ρν gβ ( T0 Th ) h ρc p ν Gr = the Grashof number, Pr = ν k Q ( T0 Th ) h the Prandtl number and Q = the heat k generation parameter. The initial and boundary conditions (4) become τ 0 : u1 = v1 = 0, θ = 0 for 1 η 1, n τ > 0 : u1 = 0, v1 = 0, θ = 1 e τ at η = 1, (16) nh where n = is the frequency parameter. ν Combining equations (13) and (14), we get where F F (1 + im) = + Gr θ F, τ 1+ m η (17) F = u 1 + iv 1 and i = 1. (18) Taking the Laplace transform of equations (17) and (15), we have d F ( a s) F Gr θ, dη + = (19) d θ Pr( Q s) θ 0, dη + = (0) where (1 + im) a =. (1) 1+ m The initial and boundary conditions for F( η, s) and θη (, s) are n F( 1, s) = 0, θ ( 1, s) =. () ss ( + n) Solutions of equations (19) and (0) subject to the boundary conditions () are given by n cos Pr( Q s) η θη (, s) =, ss ( + n) cos Pr( Q s) ngr F( η, s) = (1 Pr) s( s + b)( s + n) cos Pr( Q s) η cosh s+ aη. cos Pr( Q s) cosh s + a (3) (4)

4 48 Applied athematics and Physics The inverse Laplace transforms of equations (3) and (4) give the solution for the temperature distribution and the velocity field as cos PrQη cos Pr( Q + n) η (, ) = cos PrQ cos Pr( Q + n) (5) k s n (k 1)( 1) e 1τ π + π + cos(k + 1) η, Pr k=0 s1( s1+ n) θητ Gr cos PrQη cosh( α+ iβη ) F( ητ, ) = PrQ + a cos PrQ cosh( α + iβ) cos Pr( Q + n) η nτ Gre cos Pr( Q + n) + PrQ + a n(1 Pr) cosh( α1+ iβ1) η cosh( α1+ iβ1) (6) nπ Gr k + (k + 1)( 1) 1 Pr k=0 s e 1τ Pr s1( s1+ b)( s1+ n) π s e τ cos(k + 1) η, s( s + b)( s + n) where αβ, = ( 1 m ) 1, + ± (1 + m ) 1 1 n m α1, β1 = ± n, m 1 m m π = + (7) s1 Q (k 1), 4Pr (1 + im) π PrQ + a s = + (k+ 1) and b =. 1+ m 4 1 Pr In the absence of Hall currents ( m = 0 ), the equation (6) reduces to u1 ( ητ, ) Gr cos PrQη cosh = PrQ + cos PrQ cosh nτ Gre + PrQ + n(1 Pr) cos Pr( Q + n) η cos n cos Pr( Q + n) cos n nπ Gr k + (k + 1)( 1) 1 Pr k =0 η η π Q (k+ 1) τ 4Pr e [ π (k + 1) π (k + 1) Pr Q + n Q 4Pr 4Pr π (k + 1) + Q (1 Pr) 4Pr π + (k+ 1) τ 4 e ] (8) π (k + 1) π (k + 1) + n π (k + 1) Pr( Q + ) (1 Pr) 4 π cos(k + 1) η, v1 ( ητ, ) = 0. (9) Equation (8) is not identical with the equation (15) of Kalita [5] due to the mathematical error done by Kalita. It is also seen from figures of Kalita [5] that the temperature and velocity boundary conditions do not satisfy. 3. Results and Discussion We have presented the non-dimensional velocity components and temperature distribution for several values of magnetic parameter, Hall parameter m, heat generation parameter Q, Prandtl number Pr, frequency parameter n, Grashof number Gr and time τ against η in Figures -17. It is seen from Figure that the primary velocity u 1 decreases while the magnitude of the secondary velocity v 1 increases with an increase in magnetic parameter. This agrees with the expectations, since the magnetic field exerts a retarding force on the free convective flow. Figure 3 reveals that the primary velocity u 1 decreases whereas the magnitude of the secondary velocity v 1 increases with an increase in Hall parameter m. Figure 4 and Figure 5 show that the primary velocity u 1 and the magnitude of the secondary velocity v 1 increase with an increase in either heat generation parameter Q or Grashof number Gr. Physically, the heat generation (thermal source) has the tendency to make more heating the fluid. This behaviour is seen from Figure 4 in which the velocity components increase as Q increases. An increase in Grashof number leads to an increase in velocity, this is because, increase in Grashof number means more heating and less density. It is illustrated from Figure 6 that the primary velocity u 1 and the magnitude of the secondary velocity v 1 decrease with an increase in Prandtl number Pr. Figure 7 and Figure 8 show that both the primary velocity u 1 and the magnitude

5 Applied athematics and Physics 49 of the secondary velocity v 1 increase with an increase in either frequency parameter n or time τ. It is reveal form Figure 9 that the steady state reaches when time progresses. It is seen from Figure 10 that the fluid temperature θ increases with an increase in heat generation parameter Q. This result agrees with expectations, as Q increases, the fluid is more heated and hence the fluid temperature increases. It is observed from Figure 11 that the fluid temperature θ decreases with an increase in Prandtl number Pr. Smaller Pr fluids possess higher thermal conductivities so that heat can diffuse away from the plate surface faster than for higher Pr fluids (thicker boundary layers). Our computations show that a rise in Pr depresses the temperature function, a result consistent with numerous other studies. Figure 1 and Figure 13 show that the fluid temperature θ increases with an increase in either frequency parameter n or time τ. Figure. Velocities u 1 and v 1 for different when m = 0.5, Q = 0.5, Pr = 0.05, n =, Gr = 5 and τ = 0. Figure 3. Velocities u 1 and v 1 for different m when = 5, Q = 0.5, Pr = 0.05, n =, Gr = 5 and τ = 0.

6 50 Applied athematics and Physics Figure 4. Velocities u 1 and v 1 for different Q when = 5, m = 0.5, Pr = 0.05, n =, Gr = 5 and τ = 0. Figure 5. Velocities u 1 and v 1 for different Gr when = 5, Q = 0.5, Pr = 0.05, n =, m = 0.5 and τ = 0.

7 Applied athematics and Physics 51 Figure 6. Velocities u 1 and 1 v for different Pr when = 5, Q = 0.5, m = 0.5, n =, Gr = 5 and τ = 0. Figure 7. Velocities u 1 and v 1 for different n when = 5, Q = 0.5, m = 0.5, Pr = 0.05, Gr = 5 and τ = 0.

8 5 Applied athematics and Physics Figure 8. Velocities u 1 and v 1 for different τ when =5, Q = 0.5, m = 0.5, n =, Gr = 5 and Pr = 0.05 Figure 9. Velocities u 1 and v 1 for different η when = 5, Q = 0.5, m = 0.5, n =, Gr = 5 and Pr = 0.05

9 Applied athematics and Physics 53 Figure 10. Temperature θ for different Q when Pr = 0.05, n = and τ = 0. Figure 11. Temperature θ for different Pr when Q = 0.5, n = and τ = 0.

10 54 Applied athematics and Physics Figure 1. Temperature θ for different n for Pr = 0.05, Q = 0.5 and τ = 0. θ Figure 13. Temperature θ for different τ for Pr = 0.05, n = and Q = 0.5 The rate of heat transfer at the plates η = 1 is θ (1, τ) = and is given by η η =1 ' θ (1, τ) = PrQ tan PrQ Pr( Q + n)tan Pr( Q + n) s1τ nπ (k+ 1) e +, Pr k=0 s1( s1+ n) (30)

11 Applied athematics and Physics 55 ' Table 1. Rate of heat transfer θ (1, τ ) when n = and τ = 0.5 Q \ Pr Numerical results of the rate of heat transfer at the plate η = 1 is θ (1, τ) which is presented in Tables 1- for several values of heat generation parameter, Prandtl number Pr, frequency parameter n and time τ. Table 1 shows that the rate of heat transfer θ (1, τ) increases with an increase in Prandtl number Pr. This may be explained by the fact that frictional forces become dominant with increasing values of Pr and hence yield greater heat transfer rates. On the other hand, the rate of heat transfer θ (1, τ) decreases with an increase in heat generation parameter Q. Table shows that the rate of heat transfer θ (1, τ) increases when time τ progresses. The rate of heat transfer θ (1, τ) increases with an increase in frequency parameter n. The negative value of θ (1, τ) physically explains that there is heat flow from the plate to the fluid. ' Table. Rate of heat transfer θ (1, τ ) when Q = 0.5 and Pr = 0.05 n \ τ The non-dimensional shear stresses at the plate η =1 is given by F τx + iτy = η η =0 Gr = PrQ tan PrQ + ( α + iβ) tanh( α + iβ) PrQ + a nτ Gre + Pr( Q + n)cos Pr( Q + n) PrQ + a n(1 Pr) + ( α + iβ )tanh( α + iβ ) π Gr n (k + 1) (1 Pr) k=0 ] s1τ sτ e e, Prs1( s1+ b)( s1+ n) s( s + b)( s + n) where α, β, α 1, β 1, s 1 and s are given by (31). (31) Figure 14. Shear stresses τ x and τ y for different Numerical results of the non-dimensional shear stresses τ x and τ y at the plate ( η = 1) due to the primary and the secondary flows are plotted in Figures against Hall parameter m for several values of magnetic parameter when Gr = 5, Q = 0.5, Pr = 0.05, n = and τ = 0., heat generation parameter Q, Prandtl number Pr, Grashof number Gr, frequency parameter n and time τ. Figure 14 shows that the magnitude of the shear stress τ x due to the primary flow decreases while the shear stress

12 56 Applied athematics and Physics τ y at the plate η = 1 due to the secondary flow increases with an increase in magnetic parameter. It is seen from Figure 15 that the magnitude of the shear stress τ x and the shear stress τ y increases with an increase in heat generation parameter Q. Figure 16 that the magnitude of the shear stress τ x and the shear stress τ y decreases with an increase in Prandtl number Pr. It is seen from Figures that the magnitude of the shear stress τ x and the shear stress τ y increases with an increase in either Grashof number Gr or frequency parameter n or time τ. These results are in agrement with the fact that the velocity increases with an increase in Gr or n or τ. Figure 15. Shear stresses τ x and τ y for different Q when =5, Gr = 5, Pr = 0.05, n = and τ = 0. Figure 16. Shear stresses τ x and τ y for different Pr when =5, Gr =5, Q = 0.5, n =

13 Applied athematics and Physics 57 Figure 17. Shear stresses τ x and τ y for different Gr when =5, Q = 0.5, Pr = 0.05, n = and τ = 0. Figure 18. Shear stresses τ x and τ y for different n when =5, Gr =5, Pr = 0.05, Q = 0.5 and τ = 0.

14 58 Applied athematics and Physics Figure 19. Shear stresses τ x and 4. Conclusion τ y for different timeτ when Hall effects on an unsteady free convective flow of a viscous incompressible electrically conducting fluid between two vertical heated plates in the presence of transverse applied magnetic field and heat generation have been studied. It is found that magnetic field retards the fluid motion due to the opposing Lorentz force generated by the magnetic field. Hall parameter increases the velocity components. An increase in either heat generation parameter or Grashof number or frequency parameter leads to increase the primary velocity and the magnitude of the secondary velocity. The fluid temperature increases with an increase in either heat generation parameter or frequency parameter. Further, the Hall parameter or heat generation parameter or Grashof number or frequency parameter increases the shear stresses at the right plate. Prandtl number or frequency parameter tends to enhance the heat transfer efficiency. This model is useful to gain a deeper knowledge of the various industrial processes. Acknowledgements The authors would like to express thanks to the anonymous referees for their valuable suggestions. References [1] Turner, J. S., Buoyancy Effects in Fluids, Cambridge University Press, Cambridge, [] Jaluria, Y., Natural Convection Heat and ass Transfer, Pergamon Press, Oxford, =5, Gr =5, Pr = 0.05, n = and Q = 0.5 [3] Kakac, S., Aung, W. and Viskanta, R., eds., Natural Convection: Fundamentals and Applications, Hemisphere Publishing, New York, [4] Gebhart, B., Jaluria, Y., ahajan, R. L., and Sammakia, B., Buoyancy Induced Flows and Transport, Hemisphere Publishing, New York, [5] Bejan, A. and Kraus, A.D., Heat transfer handbook, John Wiley & Sons, Inc., Hoboken, New Jersey, 003. [6] Vajravelu, K. and Hadjinicolaou, A., (1997). Convective heat transfer in an electrically conducting fluid at a stretching surface with uniform free stream. Int. J. Eng. Sci [7] Sato, H., The Hall effects in the viscous flow of ionized gas between parallel plates under transverse magnetic field. J. Phys. Soc. Japn [8] iyatake, O. and Fujii, T., Free convection heat transfer between vertical plates - one plate isothermally heated and other thermally insulated. Heat Transfer Japan Res [9] Tanaka, H., iyatake, O., Fujii, T. and Fujii,., Natural convection heat transfer between vertical parallel plates- one plate with a uniform heat flux and the other thermally insulated, Heat Transfer Japan Res [10] Datta, N. and Jana, R. N., Hall effects on hydromagnetic convective flow through a channel with conducting walls. Int. J. Eng. Sci. 15, [11] Joshi, H.., Transient effects in natural convection cooling of vertical parallel plates. Int. Comm. Heat and ass Transfer [1] Singh, A. K., Natural convection in unsteady Couette motion. Defense Science Journal 38(1) [13] Singh, A. K., Gholami, H. R. and Soundalgekar, V.., Transient free convection flow between two vertical parallel plates. Heat and ass Transfer [14] Jha, B. K., Natural convection in unsteady HD Couette flow. Heat and ass Transfer [15] Narahari,. Sreenadh, S. and Soundalgekar, V.., Transient free convection flow between long vertical parallel plates with constant heat flux at one boundary. J. Thermophysics and Aeromechanics. 9() [16] Jha, B. K., Singh, A. K. and Takhar, H. S., Transient free convection flow in a vertical channel due to symmetric heating, Int. J. Appl. ech. Eng., 8(3),

15 Applied athematics and Physics 59 [17] Singh, A. K. and Paul, T., Transient natural convection between two vertical walls heated/cooled asymetrically. Int. J. Appl. ech. Eng. 11(1) [18] Guria,. and Jana, R. N., Hall effects on the hydromagnetic convective flow through a rotating channel under general wall conditions. agnetohydrodynamics. 43(3) [19] Jha, B. K. and Ajibade, A. O., Unsteady free convective Couette flow of heat generating/absorbing fluid. Int. J. Energy and Tech. (1) [0] Rajput, U.S. and Sahu, P.K., Transient free convection HD flow between two long vertical parallel plates with constant temperature and variable mass diffusion. Int. J. ath. Analysis. 5(34) [1] Das, S., Sarkar, B. C. and Jana, R. N., Radiation effects on free convection HD Couette flow started exponentially with variable wall temperature in presence of heat generation. Open J. Fluid Dynamics [] andal, C., Das, S. and Jana, R. N., Effect of radiation on transient natural convection flow between two vertical walls. Int. J. Appl. Inf. Systems. () [3] Das, S., Guchhait, S. K. and Jana, R. N., Radiation effects on unsteady HD free convective Couette flow of heat generation/absorbing fluid. Int. J. Com. Appl. 39(3) [4] Sarkar, B. C., Das, S. and Jana, R. N., Oscillatory HD free convective flow between two vertical walls in a rotating system. Advance in Applied Science Reasearch. 3(5) [5] Kalita, B., Unsteady free convection HD flow and heat transfer between two heated vertical paltes with heat source: an exact solution, J. Applied athematics & Bioinformatics, (3) [6] Cowling, T. G., agnetohydrodynamics, New York, Intersscience, [7] Shereliff, J.A., A text book of magnetohydrodynamic-1st ed. London, Pergamon Press, 1965.