Advances in Numerical Analysis

Volumeย 2011, Article IDย 304124, 17 pages

http://dx.doi.org/10.1155/2011/304124

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

^{1}Aerospace Engineering, Department of Engineering and Mathematics, Sheaf Building, Sheffield Hallam University, Sheffield S1 1WB, UK^{2}Thermal Engineering and Fluids Department, Technical University of Cartagena, Campus Muralla del Mar, 30202 Cartagena, Spain^{3}Magnetohydrodynamics, Applied Mathematics Program, Department of Mathematics, Narajole Raj College, P.O. Narajole, Midnapore 721 211, WB, India^{4}Institute 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.

#### Abstract

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 () strongly decelerates the flow and for very high strength magnetic fields , 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 at both small and large times; however this overshoot vanishes for larger strengths of the transverse magnetic field . 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, , 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: where is time. The magnetic flux density, , is defined as follows: where is the magnetic field strength and is the magnetic permeability. The generalized Ohmโs law defines the total current flow as follows: where is the velocity vector and is the electrical conductivity of the fluid. The electromagnetic retarding force, , to be incorporated into the momentum conservation equation, then takes the form: 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 The appropriate boundary conditions at the wall and in the free stream are 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, is fluid temperature, is free stream temperature, 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, the component of magnetic field in the -direction, 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, 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 may be expressed as a linear function of the temperature, a power-series expansion of about , neglecting higher order terms leads to Implementing (2.8) and (2.9) in (2.6), we arrive at the modified energy conservation equation: 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: 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, is the Darcian drag force coefficient (inverse permeability parameter), and is the Hartmann magnetohydrodynamic parameter. The transformed equations are thereby reduced to the following pair of coupled, second-order partial differential equations:

The corresponding transformed boundary conditions become
We note that the optically-thick radiative approximation is valid for relatively low values of the parameter, . The *electrically nonconducting* version (i.e., with ) 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: 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: in which .

#### 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
The electrical analogy is applied to (4.1) together with Kirchhoffโs law for the currents. To implement the boundary conditions at and , 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).

#### 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: , , , , , and , 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 . In Table 1, we observe that an increase in Hartmann number () 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., ) are usually associated with the formation of a Hartmann boundary layer [39]. In Table 2, the shear stress, , is found to be *decreased *significantly with an increase in Hartmann number () 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, , 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 , shear stress, increases (for the NSM solutions) from โ2.04602 for (thermal conduction exceeds thermal radiation) to โ2.0252 for (for which the thermal radiation and thermal conduction mode contributions are approximately the same), to โ1.9960 for (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 () are illustrated. A strong decrease is observed in shear stress from to , after which profiles, although they continue to decrease with increasing values, tend for , 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 (very strong transverse magnetic field) which becomes negative for , 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 (); 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 * 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, , which decelerates the flow and reduces the shear stress at the plate. Steady-state values are achieved faster with lower Darcian drag () than with higher Darcian drag (). 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, ) 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 () 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 , , but with . 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 , 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, ) 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 . For the case of , the boundary layer thicknesses will be approximately of the same order of magnitude. For , 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 (). 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 that is, 1.0 and 5.0; with this overshoot is clearly *suppressed* owing to stronger resistance to the flow. We note that for , the profiles are always greater in value than for 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 () 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 and , 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 , at , for (maximum thermal radiation effect), reaches a value of approximately 0.65, whereas the corresponding value for 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 (); 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 .

#### 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.

#### References

- 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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. - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - R. Siegel and J. R. Howell, โThermal radiation heat transfer,โ
*Hemisphere*, 1993. View at Google Scholar - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - J. A. Shercliff,
*A Textbook of Magnetohydrodynamics*, Pergamon, Oxford, UK, 1965. - V. P. Isachenko, V. A. Osipova, and A. S. Sukomel,
*Heat Transfer*, Mir, Moscow, Russia, 1969. - 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 - 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 - 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 - 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 - 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 - Pspice 6.0. Irvine, California 92718. Microsim Corporation, 20 Fairbanks, 1994.
- L. W. Nagel,
*SPICE—A Computer Program to Simulate Semiconductor Circuits*, chapter 4–6, UCB-ERL, M520, University of California, Berkeley, Calif, USA, 1998.