Table of Contents
Advances in Numerical Analysis
VolumeΒ 2011, Article IDΒ 304124, 17 pages
Research Article

Unsteady Magnetohydrodynamic Heat Transfer in a Semi-Infinite Porous Medium with Thermal Radiation Flux: Analytical and Numerical Study

1Aerospace Engineering, Department of Engineering and Mathematics, Sheaf Building, Sheffield Hallam University, Sheffield S1 1WB, UK
2Thermal Engineering and Fluids Department, Technical University of Cartagena, Campus Muralla del Mar, 30202 Cartagena, Spain
3Magnetohydrodynamics, Applied Mathematics Program, Department of Mathematics, Narajole Raj College, P.O. Narajole, Midnapore 721 211, WB, India
4Institute for Advanced Studies, Tehran 14456-63543, Iran

Received 12 December 2010; Accepted 20 January 2011

Academic Editor: William J.Β Layton

Copyright Β© 2011 O. Anwar BΓ©g et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


The unsteady, buoyancy-induced, hydromagnetic, thermal convection flow in a semi-infinite porous regime adjacent to an infinite hot vertical plate moving with constant velocity, is studied in the presence of significant thermal radiation. The momentum and energy conservation equations are normalized and then solved using both the Laplace transform technique and Network Numerical Simulation. Excellent agreement is obtained between both analytical and numerical methods. An increase in Hartmann number (𝑀2) strongly decelerates the flow and for very high strength magnetic fields (𝑀2=20), the flow is reversed after a short time interval. The classical velocity overshoot is also detected close to the plate surface for low to intermediate values of 𝑀2 at both small and large times; however this overshoot vanishes for larger strengths of the transverse magnetic field (𝑀2=10). An increase in radiation-conduction parameter (πΎπ‘Ÿ) significantly increases temperature throughout the porous regime at both small and larger times, adjacent to the plate, but decreases the shear stress magnitudes at the plate. Temperature gradient is reduced at the plate surface for all times, with a rise in radiation-conduction parameter (πΎπ‘Ÿ). Shear stress is reduced considerably with an increase in Darcian drag parameter (𝐾𝑝).

1. Introduction

Transient MHD (magnetohydrodynamic) flows with and without heat transfer in electrically-conducting fluids have attracted substantial interest in the context of metallurgical fluid dynamics, re-entry aerothermodynamics, astronautics, geophysics, nuclear engineering, and applied mathematics. An early study was presented by Carrier and Greenspan [1] who considered unsteady hydromagnetic flows past a semi-infinite flat plate moving impulsively in its own plane. Gupta [2] considered unsteady magneto-convection under buoyancy forces. Singer [3] further assessed the unsteady free convection heat transfer with magnetohydrodynamic effects in a channel regime. Pop [4] reported on transient buoyancy-driven convective hydromagnetics from a vertical surface. Yu and Yang [5] investigated the influence of channel wall conductance on hydromagnetic convection. Rao [6] analyzed the unsteady magnetohydrodynamic convection heat transfer past an infinite plane. Further excellent studies of unsteady free convection magnetohydrodynamic flows were reported by Antimirov and Kolyshkin [7] for a vertical pipe and Rajaram and Yu for a parallel-plate channel [8]. Tokis [9] used Laplace transforms to analyze the three-dimensional free-convection hydromagnetic flow near an infinite vertical plate moving in a rotating fluid when the plate temperature undergoes a thermal transient. The influence of oscillatory pressure gradient on transient rotating hydromagnetic flow was considered by Ghosh [10]. Other transient MHD studies include the papers by Sacheti et al. [11], Attia [12] who included viscosity variation effects, Al-Nimr and Alkam [13] who considered open-ended vertical annuli, and Takhar et al. [14] who employed a numerical method to study flat-plate magnetohydrodynamic unsteady convection flow. Eswara et al. [15] examined the transient laminar magnetohydrodynamic convection in a cone due to a point sink, with the free stream velocity varying continuously with time and also for the case of an impulsive change either in the strength of the point sink or in the wall temperature. They showed numerically that magnetic field increases the skin friction but decreases heat transfer and that the transient nature of the convection flow is active for short durations with suction present and greater times with injection. Chamkha [16] has analyzed the unsteady MHD free three-dimensional convection over an inclined permeable surface with heat generation/absorption. Jha [17] presented exact solutions for transient-free convection MHD Couette channel flow with impulsive motion of one of the plates discussed. More recent communications on unsteady hydromagnetic heat transfer flows include the articles by Seddeek [18] incorporating variable viscosity effects, Zakaria [19] who considered a polar fluid, and Ghosh and Pop [20] who included Hall currents [21]. Zueco presented network simulation solutions for the transient natural convection MHD flow with viscous heating effects. BΓ©g et al. [22] studied the free convective MHD flow from a spinning sphere with impulsive motion using the Blottner difference method. Duwairi et al. [23] analyzed the unsteady MHD natural convection for the non-Boussinesq case that is, using a nonlinear density relationship for water at low temperatures. Bozkaya and Tezer-sezgin [24] have presented boundary element numerical solutions for transient magnetohydrodynamic flow in a rectangular duct with insulating walls, showing that 𝑛 increase in Hartmann number causes the formation of boundary layers for both the velocity and the induced magnetic field with the velocity becoming uniform at duct centre. With increasing magnetic field, the time for reaching steady-state solution is also reduced. In many industrial applications, hydromagnetic flows also occur at very high temperatures in which thermal radiation effects become significant. The vast majority of radiation-convection flows have utilized algebraic flux approximations to simplify the general equations of radiative transfer [25]. The most popular of these simplifications remains the Rosseland diffusion approximation which has been employed by for example, Ali et al. [26] and later by Hossain et al. [27]. Radiation magnetohydrodynamic convection flows are also important in astrophysical and geophysical regimes. Raptis and Massalas [28] considered induced magnetic field effects in their study of unsteady hydromagnetic-radiative free convection. Abd El-Naby et al. [29] numerically studied magnetohydrodynamic (MHD) transient natural convection-radiation boundary layer flow with variable surface temperature, showing that velocity, temperature, and skin friction are enhanced with a rise in radiation parameter increases, whereas Nusselt number is reduced. Ogulu and Motsa [30] studied the transient radiation-magnetohydrodynamic Couette flow in a channel with wall conductance variation, showing that heat transfer rates are boosted with a reduction in Hartmann number when the wall is nonconducting, but conversely increase with a rise in values of the magnetic Reynolds number at a given value of electrical conductance. Abd-El Aziz [31] studied the thermal radiation flux effects on unsteady MHD micropolar fluid convection. Ogulu and Prakash [32] obtained analytical solutions for variable suction and radiation effects on dissipative-free convective, optically-thin, magnetohydrodynamic flow using a differential approximation to describe the radiative flux. More recent studies involving thermal radiation and transient hydromagnetic convection include the analyses by Prasad et al. [33] which included species transfer and Zueco [34] who also considered viscous heating. In numerous geophysical and metallurgical flows, porous media may also arise. Classically the Darcian model is used to simulate the bulk effects of porous materials on flow dynamics and is valid for Reynolds numbers based on the pore radius, up to approximately 10. Chamkha [35] studied the transient-free convection magnetohydrodynamic boundary layer flow in a fluid-saturated porous medium channel, and later [36] extended this study to consider the influence of temperature-dependent properties and inertial effects on the convection regime. BΓ©g et al. [37] presented perturbation solutions for the transient oscillatory hydromagnetic convection in a Darcian porous media with a heat source present. Chaudhary and Jain [38] studied the influence of oscillating temperature on magnetohydrodynamic convection heat transfer past a vertical plane in a Darcian porous medium. In the present study, we shall consider the transient radiation-convection magnetohydrodynamic flow past a vertical plate adjacent to a Darcian regime. Both analytical and numerical solutions have been obtained.

2. Mathematical Model

Consider the unsteady free convective hydromagnetic flow of an incompressible, viscous incompressible, electrically-conducting fluid along an infinite hot vertical plate moving with constant velocity, adjacent to a saturated porous regime. The physical scenario is shown in Figure 1. The π‘₯ξ…ž-axis is oriented along the plate from the leading edge in the vertically upward direction, and the π‘¦ξ…ž-axis is perpendicular to this. A uniform magnetic field, 𝐡0, is applied parallel to the π‘¦ξ…ž-axis that is, transversely to the plate. Thermal radiation acts as a unidirectional flux in the π‘¦ξ…ž-direction. The fluid is gray and absorbing-emitting but non-scattering, and the magnetic Reynolds number is assumed to be small enough to neglect induced magnetic field effects. Hall and ionslip current effects are also neglected. The electromotive force generated by a magnetic field is a function of the speed of the fluid and the magnetic field strength. Following Shercliff [39], we define the electrical field intensity, 𝐸, using Maxwell’s equation: βˆ‡Γ—πΈ=βˆ’πœ•π΅πœ•π‘‘ξ…ž,(2.1) where π‘‘ξ…ž is time. The magnetic flux density, 𝐡, is defined as follows:𝐡=πœ‡πΈH,(2.2) where 𝐻 is the magnetic field strength and πœ‡π‘’ is the magnetic permeability. The generalized Ohm’s law defines the total current flow as follows:𝐽=𝜎(𝐸+𝑉×𝐡),(2.3) where 𝑉 is the velocity vector and 𝜎 is the electrical conductivity of the fluid. The electromagnetic retarding force, 𝐹magnetic, to be incorporated into the momentum conservation equation, then takes the form: 𝐹magnetic=𝐽×𝐡=𝜎(𝑉×𝐡)×𝐡.(2.4) Incorporating this magnetic retarding force in the momentum boundary layer equation, the appropriate conservation equations, under the Boussinesq approximation for the flow under the above assumptions, may be expressed asπœ•π‘’ξ…žπœ•π‘₯ξ…ž+πœ•π‘£ξ…žπœ•π‘¦ξ…ž=0,πœ•π‘’ξ…žπœ•π‘‘ξ…žπœ•=𝜈2π‘’ξ…žπœ•π‘¦ξ…ž2ξ€·T+π‘”π›½ξ…žβˆ’Tξ…žβˆžξ€Έβˆ’πœˆπ‘’ξ…žπΎβˆ’πœŽπ΅20πœŒπ‘’ξ…ž,(2.5)πœ•Tξ…žπœ•π‘‘ξ…ž=π‘˜1πœŒπΆπ‘πœ•2Tξ…žπœ•π‘¦ξ…ž2βˆ’1πœŒπΆπ‘πœ•π‘žπ‘Ÿπœ•π‘¦ξ…ž.(2.6) The appropriate boundary conditions at the wall and in the free stream are π‘’ξ…ž=0,Tξ…ž=Tξ…žβˆžforπ‘¦ξ…žβ‰₯0,π‘‘ξ…žπ‘’β‰€0,ξ…ž=π‘ˆ,Tξ…ž=Tξ…žπ‘€forπ‘¦ξ…ž=0,π‘‘ξ…žπ‘’>0,ξ…ž=0,Tξ…žβŸΆTξ…žβˆžforπ‘¦ξ…žβŸΆβˆž,(2.7) where π‘’ξ…ž is the velocity along the plate, π‘£ξ…ž is velocity normal to the plate, 𝜈 is the kinematic viscosity of the conducting fluid, 𝑔 denotes gravity, 𝛽 is the coefficient of thermal expansion, Tξ…ž is fluid temperature, Tξ…žβˆž is free stream temperature, π‘˜1 is thermal conductivity of the fluid, 𝐢𝑝 is the specific heat capacity, 𝜌 is the fluid density, π‘žπ‘Ÿ is radiative heat flux, 𝐾 is the permeability of the porous regime, 𝐡0 the component of magnetic field in the π‘₯ξ…ž-direction, Tξ…žπ‘€ is the plate temperature (isothermal), and π‘ˆ is the velocity of the plate. Following Isachenko et al. [40], we employ a diffusion-type radiation heat transfer approximation, namely,π‘žπ‘Ÿ=βˆ’4πœŽβˆ—3π‘˜βˆ—πœ•Tξ…ž4πœ•π‘¦ξ…ž,(2.8) where πœŽβˆ— and π‘˜βˆ— are respectively the Stefan-Boltzmann constant and the spectral mean absorption coefficient of the saturated medium. Assuming that the temperature differences within the saturated porous regime are sufficiently small such that Tξ…ž4 may be expressed as a linear function of the temperature, a power-series expansion of Tξ…ž4 about Tξ…žβˆž, neglecting higher order terms leads to Tξ…ž4=4Tβˆžξ…ž3Tξ…žβˆ’3Tβˆžξ…ž4.(2.9) Implementing (2.8) and (2.9) in (2.6), we arrive at the modified energy conservation equation: πœ•Tξ…žπœ•π‘‘ξ…ž=π‘˜1πœŒπΆπ‘πœ•2Tξ…žπœ•π‘¦ξ…ž2+1πœŒπΆπ‘4πœŽβˆ—3π‘˜βˆ—πœ•2Tξ…ž4πœ•π‘¦ξ…ž2.(2.10) In order to render solutions to the boundary value problem described by (2.5) and (2.10) subject to the spatial and temporal conditions specified in (2.7), we introduce a group of nondimensional transformations, defined as follows: 𝑒𝑒=ξ…žπ‘ˆπ‘¦,𝑦=ξ…žπ‘ˆπœˆπ‘‘,𝑑=ξ…žπ‘ˆ2𝜈T,πœƒ=ξ…žβˆ’Tξ…žβˆžTξ…žπ‘€βˆ’Tξ…žβˆž,π‘ƒπ‘Ÿ=πœŒπœˆπΆπ‘π‘˜1,πΊπ‘Ÿ=ξ€·Tπ‘”π›½πœˆξ…žπ‘€βˆ’Tξ…žβˆžξ€Έπ‘ˆ3,πΎπ‘Ÿ=16πœŽβˆ—Tβˆžξ…ž33π‘˜βˆ—π‘˜1,𝐾2𝑝=π‘”πœˆπΎπ‘ˆ2,𝑀2=𝜎𝐡20πœˆπœŒπ‘ˆ2,(2.11) where 𝑒 is dimensionless π‘₯ξ…ž-direction velocity, 𝑦 is dimensionless coordinate normal to the plane, 𝑑 is dimensionless time, πœƒ is dimensionless temperature, π‘ƒπ‘Ÿ is Prandtl number, πΊπ‘Ÿ is Grashof number, πΎπ‘Ÿ is the radiation-conduction parameter, 𝐾2𝑝 is the Darcian drag force coefficient (inverse permeability parameter), and 𝑀2 is the Hartmann magnetohydrodynamic parameter. The transformed equations are thereby reduced to the following pair of coupled, second-order partial differential equations:πœ•π‘’=πœ•πœ•π‘‘2π‘’πœ•π‘¦2+πΊπ‘Ÿπœƒβˆ’πΎ2π‘π‘’βˆ’π‘€2𝑒,(2.12)ξ€·1+πΎπ‘Ÿξ€Έπœ•2πœƒπœ•π‘¦2βˆ’π‘ƒπ‘Ÿπœ•πœƒπœ•π‘‘=0.(2.13)

Figure 1: Physical regime and coordinate system.

The corresponding transformed boundary conditions become𝑒=0,πœƒ=0for𝑦β‰₯0,𝑑≀0,𝑒=1,πœƒ=1for𝑦=0,𝑑>0,𝑒=0,πœƒβŸΆ0forπ‘¦βŸΆβˆž.(2.14) We note that the optically-thick radiative approximation is valid for relatively low values of the parameter, πΎπ‘Ÿ. The electrically nonconducting version (i.e., with 𝑀2=0) of (2.12) has recently been studied by Ghosh and BΓ©g [41] where extensive computations were provided of the influence of thermal radiation on the flow field. In the present study, we shall consider the supplementary influence of transverse magnetic field for the case where the fluid is saturated with air for which the Prandtl number is assumed to take the value 0.7.

3. Analytical Solution

The Laplace transform technique is now employed to generate closed-form solutions for the coupled, linear partial differential equations (2.12) and (2.13) subject to the boundary conditions (2.14). The solutions for the transient velocity (𝑒) and transient temperature (πœƒ) take the following expressions:1𝑒(𝑦,𝑑)=2𝐺1βˆ’π‘Ÿπ·2exp(βˆ’π·π‘¦)erfcπ‘¦βˆ’2𝐷𝑑2βˆšπ‘‘ξƒͺ+exp(𝐷𝑦)erfc𝑦+2𝐷𝑑2βˆšπ‘‘+𝐺ξƒͺξƒ­π‘Ÿπ·2𝑦erfc2ξƒŽπ‘ƒπ‘Ÿξ€·1+πΎπ‘Ÿξ€Έπ‘‘ξƒͺ+πΊπ‘Ÿπ·2𝑦2βˆšπœ‹ξ‚€βˆšξ‚π‘‘2βˆ’1βˆ’1/22Γ—ξƒ¬ξ‚»βˆ’ξ‚΅π‘¦exp24𝑑+𝐷2π‘‘βˆ’ξƒŽξ‚Άξ‚Όπ‘ƒπ‘Ÿ1+πΎπ‘Ÿξ‚΅βˆ’π‘¦exp2𝑃4π‘‘π‘Ÿ1+πΎπ‘Ÿξ‚Άξƒ­,ξƒ©π‘¦πœƒ(𝑦,𝑑)=erfc2ξƒŽπ‘ƒπ‘Ÿξ€·1+πΎπ‘Ÿξ€Έπ‘‘ξƒͺ.(3.1) The spatial gradients of these functions provide expressions for the dimensionless shear stress (i.e., related to the skin friction) and temperature gradient (i.e., related to the Nusselt number) at the plate surface: πœ•π‘’||||πœ•π‘¦π‘¦=0=12𝐺1βˆ’π‘Ÿπ·2ξ‚Άξƒ¬ξ‚€βˆš(βˆ’π·)erfcβˆ’π·π‘‘ξ‚ξ‚€π·βˆš+(𝐷)erfcπ‘‘ξ‚βˆ’2βˆšξ€·π‘‘πœ‹expβˆ’π·2π‘‘ξ€Έξƒ­βˆ’πΊπ‘Ÿπ·21βˆšπœ‹ξƒ¬ξƒŽπ‘ƒπ‘Ÿξ€·1+πΎπ‘Ÿξ€Έπ‘‘ξƒ­+πΊπ‘Ÿπ·214βˆšπœ‹ξ‚€βˆšξ‚π‘‘2βˆ’1βˆ’12expβˆ’π·2π‘‘ξ€Ύβˆ’ξƒŽπ‘ƒπ‘Ÿ1+πΎπ‘Ÿξƒ­,πœ•πœƒ||||πœ•π‘¦π‘¦=01=βˆ’βˆšπœ‹ξƒ©ξƒŽπ‘ƒπ‘Ÿξ€·1+πΎπ‘Ÿξ€Έπ‘‘ξƒͺ(3.2) in which 𝐷2=𝑀2+𝐾2𝑝.

4. Network Numerical Solution

Numerical solutions to the two-point boundary value problem have also been obtained with the Network Simulation Method (NSM). This powerful and robust computational method has been employed extensively by the authors in a wide spectrum of both linear and nonlinear steady and unsteady magnetohydrodynamic and thermal convection flows. Zueco [42] recently studied the periodic temperature variation effect on thermal convection in a horizontal channel. Other very recent studies employing NSM include the works of BΓ©g et al. [43] considered unsteady rotating Couette flow in a porous medium channel, [44] analyzed the magnetohydrdoynamic rotating flow in a Darcian channel with dissipation effects and Hall/ionslip currents, and [45] examined the effects of thermal stratification and non-Darcian drag on natural convection boundary layers in a porous regime. In the NSM technique, a second-order central-difference scheme is utilized to discretize the momentum and energy conservation equations and the resulting system of finite difference equations are solved employing the Pspice program [46]. A network model is subsequently designed, with component equations which are formally equivalent to the discretized ones. The electrical analogy relates the electrical current (𝐽) with the velocity flux (πœ•π‘’/πœ•π‘¦) and temperature flux (πœ•πœƒ/πœ•π‘¦), while the electrical potential (Ξ¦) is equivalent to the velocity, 𝑒 and temperature πœƒ. A number of networks are connected in series to make up the whole medium. After experimenting with a few sets of mesh sizes, a region of integration of 200 cells has been selected. Boundary conditions are subsequently added by means of special electrical devices (current or voltage control-sources) that is, resistors, capacitors, and so forth. Once the complete network model is designed, the Pspice code is employed for the numerical simulations. This code is designated the β€œelectric circuits simulator”. Using the Fourier Law, the spatial discretization of (2.12) and (2.13) gives □𝑦𝑑𝑒𝑖=ξ€·π‘’π‘‘π‘‘π‘–βˆ’β–΅π‘¦βˆ’π‘’π‘–ξ€Έβˆ’ξ€·π‘’(□𝑦/2)π‘–βˆ’π‘’π‘–+□𝑦(□𝑦/2)+πΊπ‘Ÿπœƒπ‘–β–΅π‘¦βˆ’β–΅π‘¦πΎ2π‘π‘’π‘–βˆ’β–΅π‘¦π‘€2𝑒𝑖,β–΅π‘¦π‘ƒπ‘Ÿξ€·1+πΎπ‘Ÿξ€Έπ‘‘πœƒπ‘–=ξ€·πœƒπ‘‘π‘‘π‘–βˆ’β–΅π‘¦βˆ’πœƒπ‘–ξ€Έβˆ’ξ€·πœƒ(□𝑦/2)π‘–βˆ’πœƒπ‘–+□𝑦.(□𝑦/2)(4.1) The electrical analogy is applied to (4.1) together with Kirchhoff’s law for the currents. To implement the boundary conditions at 𝑦=0 and 𝑦=1, constant voltage sources are employed for both conditions. The principal advantage of the NSM approach is that it avoids the necessity in traditional numerical difference schemes of manipulation of difference equations and the specified constraints concerning the convergence of numerical solutions. For example, the time-step used in transient problems, which is required for convergence is not a prerequisite as Pspice achieves this via sophisticated numerical algorithms largely analagous to those intrinsic to the standard difference numerical solvers, as described by Nagel [47]. Design of the model does require a comprehensive appreciation of electrical circuit theory. Momentum balance β€œcurrents” are defined systematically for each of the discretized equations and errors can be quantified in terms of the quantity of control volumes. The network model is shown in Figure 2 for the momentum equation (2.12) and Figure 3 for the energy equation (2.13).

Figure 2: Electronic network model for the momentum equation (2.12).
Figure 3: Electronic network model for the energy equation (2.13).

5. Results and Discussion

We have obtained a comprehensive range of solutions to the transformed conservation equations. To test the validity of our numerical NSM computations, we have compared the velocity and shear stress distributions in Tables 1 and 2 with the Laplace transform solutions. Very good correlation is apparent. In all computations the key thermophysical parameters have been prescribed as follows, unless otherwise stated: πΊπ‘Ÿ=10, πΎπ‘Ÿ=2.0, 𝐾𝑝=1.0, 𝑀2=5.0, π‘ƒπ‘Ÿ=0.7, and 𝑑=0.25, corresponding to free convection of air in a highly porous regime with strong magnetic field and high thermal radiation flux at intermediate time. Both Tables 1 and 2 correspond to distributions computed a short time after the initiation of motion that is, at 𝑑=0.2. In Table 1, we observe that an increase in Hartmann number (𝑀2) from 5.0 through 8 to 10 (strong magnetic flux density) causes a significant decrease in the flow velocity, 𝑒 with distance normal to the plate surface into the boundary layer (𝑦). This trend is consistent with many classical studies on magneto-convection showing that the hydromagnetic body force retards the flow that is, decelerates the fluid causing a thinning in the boundary layer thickness. Very high Hartmann numbers (i.e., 𝑀2≫1) are usually associated with the formation of a Hartmann boundary layer [39]. In Table 2, the shear stress, πœ•π‘’/πœ•π‘¦|𝑦=0, is found to be decreased significantly with an increase in Hartmann number (𝑀2) from 5.0 through 8.0, 10.0 to 12.0, for all values of the radiation-conduction parameter, πΎπ‘Ÿ. In all cases, the shear stresses are negative since the high values of Hartmann number, 𝑀2, retard the flow in the boundary layer to such an extent that reversal of the flow is caused. This result is significant in the design of, for example, MHD generators since a critical magnetic flux density may be applied (i.e., Hartmann number) to reverse the flow dynamics during operation. The change in shear stresses with an increase in radiation-conduction parameter indicates that an increase in thermal radiation has a positive effect on the flow that is, reduces the degree of flow reversal. For example, for 𝑀2=5.0, shear stress, πœ•π‘’/πœ•π‘¦|𝑦=0increases (for the NSM solutions) from βˆ’2.04602 for πΎπ‘Ÿ=0.5 (thermal conduction exceeds thermal radiation) to βˆ’2.0252 for πΎπ‘Ÿ=1.0 (for which the thermal radiation and thermal conduction mode contributions are approximately the same), to βˆ’1.9960 for πΎπ‘Ÿ=2.0 (thermal radiation exceeds thermal conduction). Therefore, for very high-strength magnetic field operating conditions, thermal radiation mitigates to some extent flow reversal effects. In Figure 4, the dimensionless shear stress profiles in time for different Hartmann numbers (𝑀2) are illustrated. A strong decrease is observed in shear stress from 𝑑=0 to 𝑑=0.05, after which profiles, although they continue to decrease with increasing 𝑀2 values, tend for 𝑀2=0, 5, 10, and 15, to the steady state. For these values of Hartmann number, profiles are always positive indicating that flow reversal does not occur. Comparing with the trends in Table 2, we can note that while the values of πΎπ‘Ÿ, 𝐾𝑝, and π‘ƒπ‘Ÿ are identical for the third column, the Grashof number (πΊπ‘Ÿ) is lower in Table 2 at 2.0, compared with the value of 10.0 in Figure 4. As such the flow is more strongly assisted by buoyancy forces in Figure 4 which prevents the reversal of flow for all the profiles, with the exception of 𝑀2=20 (very strong transverse magnetic field) which becomes negative for 𝑑>0.4, Thermal buoyancy force, +πΊπ‘Ÿπœƒ, is directly proportional to the Grashof free convection parameter (πΊπ‘Ÿ) and therefore would appear to assist the flow, whereas magnetic field inhibits flow acceleration in the regime. In Figure 5, temperature gradient profiles in time for the influence of radiation-conduction parameter (πΎπ‘Ÿ) are presented. In all cases, profiles are a maximum initially at the isothermal plate and decay quickly from the wall with time. An increase in πΎπ‘Ÿ from 0.1 through 0.5, 1, 1.5 to 2.0 is seen to markedly reduce heat transfer gradient especially at shorter times (0<𝑑<0.2); with further elapse of time all profiles converge that is, radiation effects are negligible for large times. Increasing πΎπ‘Ÿ implies a greater augmentation of heat transfer by thermal radiation which will serve to increase fluid temperatures in the regime; the spatial heat transfer rate πœ•πœƒ/πœ•π‘¦|𝑦=0 that is, temperature gradient at the wall will therefore be reduced as greater thermal energy (heat) will be imparted to the fluid-saturated regime raising temperatures within the porous regime. In Figure 6, the influence of the Darcian drag force parameter, 𝐾𝑝, on the time evolution of shear stress profiles, is depicted. As 𝐾𝑝 increases from 0.1 through 0.5, 1, 2 to 5, a very large escalation in Darcian drag force is caused, as expressed in (2.12) in the linear term, βˆ’πΎ2𝑝, which decelerates the flow and reduces the shear stress at the plate. Steady-state values are achieved faster with lower Darcian drag (𝐾𝑝=0.1) than with higher Darcian drag (𝐾𝑝=5). The effect of Grashof number on velocity gradient (i.e., shear stress) through time is presented in Figure 7. Increasing πΊπ‘Ÿ for the case of very low Darcian drag (i.e., highly permeable medium, 𝐾𝑝=0.1) strongly increases shear stress values at the wall that is, accelerates the flow over time. We note that values become negative for very low πΊπ‘Ÿ values since the magnetic impedance force (𝑀2=5.0) will dominate and have a greater inhibiting influence with low buoyancy that is, flow reversal accompanies lower thermal buoyancy forces for higher permeability regimes. In Figure 8, the effect of the Prandtl number (π‘ƒπ‘Ÿ) on the temporal shear stress distribution is shown again for πΊπ‘Ÿ=10.0, πΎπ‘Ÿ=2.0, 𝑀2=5.0 but with 𝐾𝑝=1.0. Increasing π‘ƒπ‘Ÿ strongly boosts the flow and increases shear stress profile values which remain positive for small times; however, with increasing elapse of time shear stress, values become negative indicating backflow occurs at the plate. For lower π‘ƒπ‘Ÿ values (0.7, 0.1), negative values are attained more quickly that is, backflow takes place quicker. In Figure 9, an increase in π‘ƒπ‘Ÿ is observed to the enhance temperature gradient. Prandtl number controls the relative thickness of the momentum and thermal boundary layers. When π‘ƒπ‘Ÿ is of low value, heat diffusion exceeds momentum diffusion. For π‘ƒπ‘Ÿ<1, the thickness of the thermal boundary layer therefore exceeds the thickness of the velocity boundary layer that is, temperatures will be greater. In Figure 10, temperatures are seen to decrease considerably with an increase in π‘ƒπ‘Ÿ values (for a fixed time, 𝑑=0.25) as we progress into the boundary layer regime; profiles also decay much more sharply for higher π‘ƒπ‘Ÿ values since momentum diffusion exceeds energy (heat) diffusion for π‘ƒπ‘Ÿ>1. For the case of π‘ƒπ‘Ÿ=1, the boundary layer thicknesses will be approximately of the same order of magnitude. For π‘ƒπ‘Ÿ=0.1, the profile is approximately linear for a substantial distance from the plate. Spatial velocity (𝑒) distributions, for two time values are illustrated in Figure 11, for the effect of Hartmann hydromagnetic parameter (𝑀2). This parameter represents the ratio of the hydromagnetic retarding force to the viscous hydrodynamic force in the boundary layer. The classical velocity overshoot is identified [1, 2, 4, 39] near the moving plate surface for lower values of 𝑀2 that is, 1.0 and 5.0; with 𝑀2=10.0 this overshoot is clearly suppressed owing to stronger resistance to the flow. We note that for 𝑑=1.0, the profiles are always greater in value than for 𝑑=0.25 that is, the flow is accelerated considerably with time, although velocities are strongly reduced with an increase in Hartmann number. All profiles decrease towards zero in the free stream, although this state is attained much faster for higher magnetic field values (𝑀2=10) and for shorter times. In Figure 12, the combined effects of time (𝑑) and radiation-conduction parameter (πΎπ‘Ÿ) on spatial distribution of temperature (πœƒ) through the boundary layer is shown. An increase in πΎπ‘Ÿ serves to supplement fluid thermal conductivity with radiation contribution and significantly heats the fluid-saturated regime that is, increases temperature values. Similar results were reported by Ali et al. [26], Hossain et al. [28], and very recently by Ghosh and BΓ©g [41]. A large difference is observed between the profiles computed at 𝑑=1.0 and 𝑑=0.25, indicating that thermal radiation effects are amplified at greater times, compared with smaller times where the flow is still developing. After greater times a greater quantity of thermal energy will be absorbed into the fluid regime via the imposed flux causing enhanced heating of the fluid. For example, for 𝑑=1.0, at 𝑦=2, for πΎπ‘Ÿ=5 (maximum thermal radiation effect), πœƒ reaches a value of approximately 0.65, whereas the corresponding value for 𝑑=0.25 is much lower at 0.35. Finally in Figure 13 we have plotted the spatial variation of velocity for the combined effects of radiation-conduction parameter (πΎπ‘Ÿ) and time (𝑑). Again a velocity overshoot is observed in the close vicinity of the plate; however, this overshoot is distinctly greater for the highest value of πΎπ‘Ÿ (= 2.0) and greater time values (𝑑=1.0); all profiles descend gradually to zero far from the wall. Thermal radiation therefore augments the flow that is, accelerates the flow in the porous regime. Velocities are minimized when thermal conduction swamps thermal radiation that is, for πΎπ‘Ÿ=0.1.

Table 1: Dimensionless spatial velocity distribution (𝑒) for πΊπ‘Ÿ=2.0πΎπ‘Ÿ = 1.0, 𝐾𝑝=1.0, π‘ƒπ‘Ÿ=0.7, and 𝑑=0.2 for the effect of Hartmann number (𝑀2).
Table 2: Dimensionless shear stress (πœ•π‘’/πœ•π‘¦|𝑦=0) for πΊπ‘Ÿ=2.0, 𝐾𝑝=1.0, π‘ƒπ‘Ÿ=0.7, and 𝑑=0.2 for the combined effects of radiation-conduction parameters (πΎπ‘Ÿ) and Hartmann number (𝑀2).
Figure 4: Temporal shear stress distribution (πœ•π‘’/πœ•π‘¦|𝑦=0) with πΊπ‘Ÿ=10.0, πΎπ‘Ÿ=2.0, 𝐾𝑝=1.0, and π‘ƒπ‘Ÿ=0.7 for the effect of Hartmann hydromagnetic parameter (𝑀2).
Figure 5: Temporal temperature gradient distribution (πœ•πœƒ/πœ•π‘¦|𝑦=0) with πΊπ‘Ÿ=10.0, 𝑀2=5, 𝐾𝑝=1.0, and π‘ƒπ‘Ÿ=0.7 for the effect of radiation-conduction number (πΎπ‘Ÿ).
Figure 6: Temporal shear stress distribution (πœ•π‘’/πœ•π‘¦|𝑦=0) with πΊπ‘Ÿ=10.0, 𝑀2=5, πΎπ‘Ÿ=2.0, and π‘ƒπ‘Ÿ=0.7 for the effect of Darcian drag parameter (𝐾𝑝).
Figure 7: Temporal shear stress distribution (πœ•π‘’/πœ•π‘¦|𝑦=0) with 𝐾𝑝=0.1, 𝑀2=5, πΎπ‘Ÿ=2.0, and π‘ƒπ‘Ÿ=0.7 for the effect of the Grashof number (πΊπ‘Ÿ).
Figure 8: Temporal shear stress distribution (πœ•π‘’/πœ•π‘¦|𝑦=0) with 𝐾𝑝=1, 𝑀2=5, πΎπ‘Ÿ=2.0, and πΊπ‘Ÿ = 10.0, for the effect of Prandtl number (π‘ƒπ‘Ÿ).
Figure 9: Temporal temperature gradient distribution (πœ•πœƒ/πœ•π‘¦|𝑦=0) with 𝐾𝑝=1, 𝑀2=5, πΎπ‘Ÿ=2.0, and πΊπ‘Ÿ=10.0, for the effect of Prandtl number (π‘ƒπ‘Ÿ).
Figure 10: Spatial temperature distribution (πœƒ) with 𝐾𝑝=1, 𝑀2=5, πΎπ‘Ÿ=2.0, and πΊπ‘Ÿ=10.0, for the effect of Prandtl number (π‘ƒπ‘Ÿ) at 𝑑=0.25.
Figure 11: Spatial velocity distribution (𝑒) with 𝐾𝑝=1, πΊπ‘Ÿ=10, πΎπ‘Ÿ=2.0, and π‘ƒπ‘Ÿ=0.7 for the effect of Hartmann hydromagnetic parameter at 𝑑=0.25 and 𝑑=1.
Figure 12: Spatial temperature distribution (πœƒ) with 𝐾𝑝=1, πΊπ‘Ÿ=10, 𝑀2=10, and π‘ƒπ‘Ÿ=0.7 for the effect of radiation-conduction parameter (πΎπ‘Ÿ) at 𝑑=0.25 and 𝑑=1.
Figure 13: Spatial velocity distribution (𝑒) with πΊπ‘Ÿ=10.0, 𝑀2=5, πΎπ‘Ÿ=2.0, and π‘ƒπ‘Ÿ=0.7 for the effect of radiation-conduction parameter (πΎπ‘Ÿ).

6. Conclusions

Closed form and numerical (NSM) solutions have been presented for the transient hydromagnetic natural convection boundary layer flow past a moving vertical plane adjacent to a Darcian porous regime with thermal radiation flux present. It has been shown that thermal radiation strongly increases fluid temperatures and accelerates the flow; conversely magnetic field as simulated via the Hartmann number serves to impede the flow and reduce velocity gradient (shear stress) values. The effects of both parameters are enhanced with a greater elapse of time. Darcian drag is seen to decelerate the flow, whereas increasing free convection serves to accelerate the flow owing to the assistance of thermal buoyancy forces in the regime. A rise in Prandtl number however decreases temperatures in the regime, but accelerates the flow that is, increases velocity gradient values. A velocity overshoot is observed with magnetic field effects but vanishes for very high values of the Hartmann number.


  1. G. F. Carrier and H. P. Greenspan, β€œThe time-dependent magnetohydrodynamic flow past a flat plate,” Journal of Fluid Mechanics, vol. 7, pp. 22–32, 1960. View at Google Scholar
  2. A. S. Gupta, β€œSteady and transient free convection of an electrically conducting fluid from a vertical plate in the presence of a magnetic field,” Applied Scientific Research, vol. 9, no. 1, pp. 319–333, 1960. View at Publisher Β· View at Google Scholar Β· View at Scopus
  3. R. M. Singer, β€œTransient magnetohydrodynamic flow and heat transfer,” Zeitschrift für Angewandte Mathematik und Mechanik, vol. 16, no. 4, pp. 483–494, 1965. View at Publisher Β· View at Google Scholar Β· View at Scopus
  4. I. Pop, β€œUnsteady hydromagnetic free convection flow from a vertical infinite flat plate,” Zeitschrift für Angewandte Mathematik und Mechanik, vol. 49, no. 12, pp. 756–757, 1969. View at Google Scholar
  5. C. P. Yu and H. K. Yang, β€œEffect of wall conductances on convective magnetohydrodynamic channel flow,” Applied Scientific Research, vol. 20, no. 1, pp. 16–24, 1969. View at Publisher Β· View at Google Scholar Β· View at Scopus
  6. V. V. R. Rao, β€œTransient hydromagnetic free-convection flow,” Indian Journal of Physics, vol. 45, no. 5, pp. 233–235, 1971. View at Google Scholar
  7. M. Y. Antimirov and A. A. Kolyshkin, β€œUnsteady MHD convection in a vertical circular channel,” Magnetohydrodynamics, vol. 14, no. 3, pp. 279–284, 1978. View at Google Scholar Β· View at Scopus
  8. S. Rajaram and C. P. Yu, β€œMHD channel flow in a transient magnetic field,” in ASME Symposium Forum on Unsteady Flow, 1984, Publication no. G00259.
  9. J. N. Tokis, β€œUnsteady magnetohydrodynamic free-convection flows in a rotating fluid,” Astrophysics and Space Science, vol. 119, no. 2, pp. 305–313, 1986. View at Publisher Β· View at Google Scholar Β· View at Scopus
  10. S. K. Ghosh, β€œUnsteady hydromagnetic flow in a rotating channel with oscillating pressure gradient,” Journal of the Physical Society of Japan, vol. 62, no. 11, pp. 3893–3903, 1993. View at Google Scholar Β· View at Scopus
  11. N. C. Sacheti, P. Chandran, and A. K. Singh, β€œAn exact solution for unsteady magnetohydrodynamic free convection flow with constant heat flux,” International Communications in Heat and Mass Transfer, vol. 21, no. 1, pp. 131–142, 1994. View at Google Scholar Β· View at Scopus
  12. H. A. Attia, β€œTransient MHD flow and heat transfer between two parallel plates with temperature dependent viscosity,” Mechanics Research Communications, vol. 26, no. 1, pp. 115–121, 1999. View at Google Scholar Β· View at Scopus
  13. M. A. Al-Nimr and M. K. Alkam, β€œMagnetohydrodynamics transient free convection in open-ended vertical annuli,” Journal of thermophysics and heat transfer, vol. 13, no. 2, pp. 256–265, 1999. View at Google Scholar Β· View at Scopus
  14. H. S. Takhar, A. J. Chamkha, and G. Nath, β€œUnsteady flow and heat transfer on a semi-infinite flat plate with an aligned magnetic field,” International Journal of Engineering Science, vol. 37, no. 13, pp. 1723–1736, 1999. View at Google Scholar Β· View at Scopus
  15. A. T. Eswara, S. Roy, and G. Nath, β€œUnsteady MHD forced flow due to a point sink,” Acta Mechanica, vol. 145, no. 1–4, pp. 159–172, 2000. View at Publisher Β· View at Google Scholar Β· View at Scopus
  16. A. J. Chamkha, β€œTransient hydromagnetic three-dimensional natural convection from an inclined stretching permeable surface,” Chemical Engineering Journal, vol. 76, no. 2, pp. 159–168, 2000. View at Publisher Β· View at Google Scholar
  17. B. K. Jha, β€œNatural convection in unsteady MHD couette flow,” Heat and Mass Transfer, vol. 37, no. 4-5, pp. 329–331, 2001. View at Google Scholar
  18. M. A. Seddeek, β€œEffects of radiation and variable viscosity on a MHD free convection flow past a semi-infinite flat plate with an aligned magnetic field in the case of unsteady flow,” International Journal of Heat and Mass Transfer, vol. 45, no. 4, pp. 931–935, 2001. View at Publisher Β· View at Google Scholar Β· View at Scopus
  19. M. Zakaria, β€œMagnetohydrodynamic unsteady free convection flow of a couple stress fluid with one relaxation time through a porous medium,” Applied Mathematics and Computation, vol. 146, no. 2-3, pp. 469–494, 2003. View at Publisher Β· View at Google Scholar Β· View at Scopus
  20. S. K. Ghosh and I. Pop, β€œHall effects on unsteady hydromagnetic flow in a rotating system with oscillating pressure gradient,” International Journal of Applied Mechanics and Engineering, vol. 8, no. 1, pp. 43–56, 2003. View at Google Scholar
  21. J. Z. Jordán, β€œNumerical study of an unsteady free convective magnetohydrodynamic flow of a dissipative fluid along a vertical plate subject to a constant heat flux,” International Journal of Engineering Science, vol. 44, no. 18-19, pp. 1380–1393, 2006. View at Publisher Β· View at Google Scholar
  22. O. Bég, H. S. Anwar, G. N. Takhar, and A. J. Chamkha, β€œMathematical Modeling of hydromagnetic convection from a rotating sphere with impulsive motion and buoyancy effects,” Nonlinear Analysis: Modelling and Control, vol. 11, no. 3, pp. 227–245, 2006. View at Google Scholar
  23. H. M. Duwairi, R. A. Damseh, and B. Tashtoush, β€œTransient non-boussinesq magnetohydrodynamic free convection flows over a vertical surface,” International Journal of Fluid Mechanics Research, vol. 33, no. 2, pp. 137–152, 2006. View at Publisher Β· View at Google Scholar Β· View at Scopus
  24. C. Bozkaya and M. Tezer-Sezgin, β€œBoundary element solution of unsteady magnetohydrodynamic duct flow with differential quadrature time integration scheme,” International Journal for Numerical Methods in Fluids, vol. 51, no. 5, pp. 567–584, 2006. View at Publisher Β· View at Google Scholar Β· View at Scopus
  25. R. Siegel and J. R. Howell, β€œThermal radiation heat transfer,” Hemisphere, 1993. View at Google Scholar
  26. M. M. Ali, T. S. Chen, and B. F. Armaly, β€œNatural convection radiation interaction in boundary layer flow over horizontal surfaces,” AIAA Journal, vol. 22, no. 12, pp. 1797–1803, 1984. View at Google Scholar Β· View at Scopus
  27. A. Raptis and C. V. Massalas, β€œMagnetohydrodynamic flow past a plate by the presence of radiation,” Heat and Mass Transfer, vol. 34, no. 2-3, pp. 107–109, 1998. View at Google Scholar Β· View at Scopus
  28. M. A. Hossain, D. A. S. Rees, and I. Pop, β€œFree convection-radiation interaction from an isothermal plate inclined at a small angle to the horizontal,” Acta Mechanica, vol. 127, no. 1–4, pp. 63–73, 1998. View at Google Scholar Β· View at Scopus
  29. M. A. Abd El-Naby, E. M. E. Elbarbary, and N. Y. Abdelazem, β€œFinite difference solution of radiation effects on MHD unsteady free-convection flow over vertical plate with variable surface temperature,” Journal of Applied Mathematics, vol. 2003, no. 2, pp. 65–86, 2003. View at Publisher Β· View at Google Scholar Β· View at Scopus
  30. A. Ogulu and S. Motsa, β€œRadiative heat transfer to magneto-hydrodynamic couette flow with variable wall temperature,” Physica Scripta, vol. 71, no. 4, pp. 336–339, 2005. View at Google Scholar Β· View at Scopus
  31. M. Abd-El Aziz, β€œThermal radiation effects on magnetohydrodynamic mixed convection flow of a micropolar fluid past a continuously moving semi-infinite plate for high temperature differences,” Acta Mechanica, vol. 187, no. 1–4, pp. 113–127, 2006. View at Publisher Β· View at Google Scholar Β· View at Scopus
  32. A. Ogulu and J. Prakash, β€œHeat transfer to unsteady magneto-hydrodynamic flow past an infinite moving vertical plate with variable suction,” Physica Scripta, vol. 74, no. 2, pp. 232–239, 2006. View at Publisher Β· View at Google Scholar Β· View at Scopus
  33. V. R. Prasad, N. B. Reddy, and Muthucumaraswamy, β€œTransient radiative hydromagnetic free convection flow past an impulsively started vertical plate with uniform heat and mass flux,” Theoretical Applied Mechanics, vol. 33, no. 1, pp. 31–63, 2006. View at Google Scholar
  34. J. Zueco, β€œNetwork simulation method applied to radiation and viscous dissipation effects on MHD unsteady free convection over vertical porous plate,” Applied Mathematical Modelling, vol. 31, no. 9, pp. 2019–2033, 2007. View at Publisher Β· View at Google Scholar Β· View at Scopus
  35. A. J. Chamkha, β€œUnsteady hydromagnetic natural convection in a fluid-saturated porous medium channel,” Advanced Filtration and Separation Technology, vol. 10, pp. 369–375, 1996. View at Google Scholar
  36. A. J. Chamkha, β€œUnsteady laminar hydromagnetic flow and heat transfer in porous channels with temperature-dependent properties,” International Journal of Numerical Methods for Heat and Fluid Flow, vol. 11, no. 5-6, pp. 430–448, 2001. View at Publisher Β· View at Google Scholar Β· View at Scopus
  37. O. A. Bég, H. S. Takhar, and A. K. Singh, β€œMultiparameter perturbation analysis of unsteady oscillatory magnetoconvection in porous media with heat source effects,” International Journal of Fluid Mechanics Research, vol. 32, no. 6, pp. 635–661, 2005. View at Publisher Β· View at Google Scholar Β· View at Scopus
  38. R. C. Chaudhary and A. Jain, β€œMagnetohydrodynamic transient convection flow past a vertical surface embedded in a porous medium with oscillating temperature,” Turkish Journal of Engineering and Environmental Sciences, vol. 32, no. 1, pp. 13–22, 2008. View at Google Scholar Β· View at Scopus
  39. J. A. Shercliff, A Textbook of Magnetohydrodynamics, Pergamon, Oxford, UK, 1965.
  40. V. P. Isachenko, V. A. Osipova, and A. S. Sukomel, Heat Transfer, Mir, Moscow, Russia, 1969.
  41. S. K. Ghosh and O. A. Bég, β€œTheoretical analysis of radiative effects on transient free convection heat transfer past a hot vertical surface in porous media,” Nonlinear Analysis: Modelling and Control, vol. 13, no. 4, pp. 419–432, 2008. View at Google Scholar
  42. J. Zueco, β€œUnsteady conjugate problem of a dissipative fluid in a horizontal channel with a periodic variation temperature,” Meccanica, vol. 43, no. 1, pp. 37–46, 2008. View at Publisher Β· View at Google Scholar Β· View at Scopus
  43. O. A. Bég, H. S. Takhar, J. Zueco, A. Sajid, and R. Bhargava, β€œTransient Couette flow in a rotating non-Darcian porous medium parallel plate configuration: network simulation method solutions,” Acta Mechanica, vol. 200, no. 3-4, pp. 129–144, 2008. View at Publisher Β· View at Google Scholar Β· View at Scopus
  44. O. A. Bég, J. Zueco, and H. S. Takhar, β€œUnsteady magnetohydrodynamic Hartmann-Couette flow and heat transfer in a Darcian channel with Hall current, ionslip, viscous and Joule heating effects: network numerical solutions,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 4, pp. 1082–1097, 2009. View at Publisher Β· View at Google Scholar Β· View at Scopus
  45. O. A. Bég, J. Zueco, and H. S. Takhar, β€œLaminar free convection from a continuously-moving vertical surface in thermally-stratified non-Darcian high-porosity medium: network numerical study,” International Communications in Heat and Mass Transfer, vol. 35, no. 7, pp. 810–816, 2008. View at Publisher Β· View at Google Scholar Β· View at Scopus
  46. Pspice 6.0. Irvine, California 92718. Microsim Corporation, 20 Fairbanks, 1994.
  47. L. W. Nagel, SPICE—A Computer Program to Simulate Semiconductor Circuits, chapter 4–6, UCB-ERL, M520, University of California, Berkeley, Calif, USA, 1998.