- About this Journal
- Abstracting and Indexing
- Aims and Scope
- Article Processing Charges
- Articles in Press
- Author Guidelines
- Bibliographic Information
- Citations to this Journal
- Contact Information
- Editorial Board
- Editorial Workflow
- Free eTOC Alerts
- Publication Ethics
- Reviewers Acknowledgment
- Submit a Manuscript
- Subscription Information
- Table of Contents
Advances in Numerical Analysis
Volume 2011 (2011), Article ID 304124, 17 pages
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 () 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 ().
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  who considered unsteady hydromagnetic flows past a semi-infinite flat plate moving impulsively in its own plane. Gupta  considered unsteady magneto-convection under buoyancy forces. Singer  further assessed the unsteady free convection heat transfer with magnetohydrodynamic effects in a channel regime. Pop  reported on transient buoyancy-driven convective hydromagnetics from a vertical surface. Yu and Yang  investigated the influence of channel wall conductance on hydromagnetic convection. Rao  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  for a vertical pipe and Rajaram and Yu for a parallel-plate channel . Tokis  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 . Other transient MHD studies include the papers by Sacheti et al. , Attia  who included viscosity variation effects, Al-Nimr and Alkam  who considered open-ended vertical annuli, and Takhar et al.  who employed a numerical method to study flat-plate magnetohydrodynamic unsteady convection flow. Eswara et al.  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  has analyzed the unsteady MHD free three-dimensional convection over an inclined permeable surface with heat generation/absorption. Jha  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  incorporating variable viscosity effects, Zakaria  who considered a polar fluid, and Ghosh and Pop  who included Hall currents . Zueco presented network simulation solutions for the transient natural convection MHD flow with viscous heating effects. Bég et al.  studied the free convective MHD flow from a spinning sphere with impulsive motion using the Blottner difference method. Duwairi et al.  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  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 . The most popular of these simplifications remains the Rosseland diffusion approximation which has been employed by for example, Ali et al.  and later by Hossain et al. . Radiation magnetohydrodynamic convection flows are also important in astrophysical and geophysical regimes. Raptis and Massalas  considered induced magnetic field effects in their study of unsteady hydromagnetic-radiative free convection. Abd El-Naby et al.  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  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  studied the thermal radiation flux effects on unsteady MHD micropolar fluid convection. Ogulu and Prakash  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.  which included species transfer and Zueco  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  studied the transient-free convection magnetohydrodynamic boundary layer flow in a fluid-saturated porous medium channel, and later  extended this study to consider the influence of temperature-dependent properties and inertial effects on the convection regime. Bég et al.  presented perturbation solutions for the transient oscillatory hydromagnetic convection in a Darcian porous media with a heat source present. Chaudhary and Jain  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 , 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. , 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  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  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.  considered unsteady rotating Couette flow in a porous medium channel,  analyzed the magnetohydrdoynamic rotating flow in a Darcian channel with dissipation effects and Hall/ionslip currents, and  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 . 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 . 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 . 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. , Hossain et al. , and very recently by Ghosh and Bég . 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 .
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.
- 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.
- 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.
- R. M. Singer, “Transient magnetohydrodynamic flow and heat transfer,” Zeitschrift für Angewandte Mathematik und Mechanik, vol. 16, no. 4, pp. 483–494, 1965.
- 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.
- 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.
- V. V. R. Rao, “Transient hydromagnetic free-convection flow,” Indian Journal of Physics, vol. 45, no. 5, pp. 233–235, 1971.
- M. Y. Antimirov and A. A. Kolyshkin, “Unsteady MHD convection in a vertical circular channel,” Magnetohydrodynamics, vol. 14, no. 3, pp. 279–284, 1978.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- B. K. Jha, “Natural convection in unsteady MHD couette flow,” Heat and Mass Transfer, vol. 37, no. 4-5, pp. 329–331, 2001.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- R. Siegel and J. R. Howell, “Thermal radiation heat transfer,” Hemisphere, 1993.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.