`ISRN ThermodynamicsVolume 2012 (2012), Article ID 214362, 9 pageshttp://dx.doi.org/10.5402/2012/214362`
Research Article

## Radiative Boundary Layer Flow in Porous Medium due to Exponentially Shrinking Permeable Sheet

1Department of Mathematics, University of Rajasthan, Jaipur 302004, India
2Department of Mathematics, Poornima University, Jaipur 303905, India

Received 24 September 2012; Accepted 10 October 2012

Academic Editors: D. Bratko and Q. Zhang

Copyright © 2012 Paresh Vyas and Nupur Srivastava. 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

This communication pertains to the study of radiative heat transfer in boundary layer flow over an exponentially shrinking permeable sheet placed at the bottom of fluid saturated porous medium. The porous medium has permeability of specified form. The fluid considered here is Newtonian, without phase change, optically dense, absorbing-emitting radiation but a nonscattering medium. The setup is subjected to suction to contain the vorticity in the boundary layer. The radiative heat flux in the energy equation is accounted by Rosseland approximation. The thermal conductivity is presumed to vary with temperature in a linear fashion. The governing partial differential equations are reduced to ordinary differential equations by similarity transformations. The resulting system of nonlinear ordinary differential equations is solved numerically by fourth-order Runge-Kutta scheme together with shooting method. The pertinent findings displayed through figures and tables are discussed.

#### 1. Introduction

Flow and heat transfer in boundary layer flow of viscous fluid due to deforming surface is pivotal in many industrial processes cutting across different realms. Crane [1] pioneered a closed form analytical solution for an incompressible fluid flow due to a linearly stretching sheet. This novel flow configuration in fluid mechanics was well received by the investigators to extract pertinent information in real analogous systems like extrusion processes, hot rolling, wire drawing and glass wire production, and so forth, to name a few. The so-called stretching surface flow problem has been studied for variety of assumptions pertaining to surface velocity (linear/power law) and surface temperature (constant/variable)/temperature flux for the simple reason that optimal quality of desired component squarely depends on the shear and heat transfer rate at the sheet. The problem has also been discussed for non-Newtonian fluid, Ferro-fluid and micropolar fluid and also by incorporating physical features like suction/injection, MHD, chemical reaction, and so forth [29]. Recently, Liao [10, 11] discovered multiple solutions for the flow over both impermeable and permeable stretching sheets. In order to complete the discussion, attempts were also made to explore flow and/or heat transfer due to rather unusual exponential stretching surface. Magyari and Keller [12] discussed boundary layer flow due to a sheet stretching exponentially and bearing a variable temperature of exponential fashion. Elbashbeshy [13] furnished numerical findings about flow and heat transfer over an exponentially stretching surface subjected to mass suction. Exponentially stretching surface in the presence of magnetic field was analysed by Al-Odat et al. [14]. Partha et al. [15] considered dissipative mixed convection flow past an exponentially stretching surface. The problem was also extended to viscoelastic fluids by Sanjayanand and Khan [16, 17]. Sajid and Hayat [18] examined radiative boundary layer flow due to exponentially stretching sheet and reported HAM series solution.

Looking at the literature of flow due to deforming boundary, one finds that flow due to shrinking boundary is still a relatively low explored area. Investigation of boundary layer flow and heat transfer past shrinking surfaces requires much attention as it finds important applications in many industries and manufacturing processes. A heat shrinking sheet provides a high quality field coating solution for welded pipe ends. A heat shrinking tape is used for wrapping junctions. Wang [19] was probably the first to report the boundary layer flow over an unsteady shrinking film. In a subsequent study, Miklavčič and Wang [20] reasoned that flow due to shrinking sheet cannot exist unless adequate suction is applied across the boundary layer to contain the vorticity. In their study they reported a pertinent finding that in shrinking flow problem one may come across multiple solutions for different suction rates. Wang [21] studied two-dimensional and axisymmetric stagnation point flow towards a shrinking sheet and concluded that solutions do not exist for larger shrinking rates. Appreciating the possibility of exploring new avenues in flow arising due to shrinking surface, some attempts have been made recently [2229]. MHD flow and heat transfer in a channel bounded by a shrinking sheet and a plate with a porous substrate was discussed by Chauhan and Agarwal [30]. Though the work of Bhattacharyya [31] throws some light on boundary layer flow and heat transfer past an exponentially shrinking sheet, however, the study of flow dynamics past an exponentially shrinking surface in porous media is still unexplored.

All the above-said investigations regarding “shrinking surface problems” consider clear fluid region with constant thermal conductivity. The thermal conductivity is an intrinsic property of fluids which exhibits variations with temperature. Many authors have reported studies with variable thermal conductivity aspects [3234]. Khan et al. [35] studied thin film flow over a shrinking/stretching sheet taking variable viscosity and thermal conductivity into account.

The radiative thermal regime in porous medium has drawn much attention recently due to ample applications, such as, gasification of oil shale, waste heat storage in aquifers, and so forth. To be specific, in the case of gasification, large temperature gradient exists in the vicinity of the combustion regime making radiation effect dominant. However, literature is rather scanty on the radiative flow in porous media. A few relevant texts/studies have been presented [3639]. Radiative flow of Jeffery fluid with variable thermal conductivity in a porous medium was discussed by Hayat et al. [40]. Elbashbeshy and Emam [41] discussed the effects of radiation and heat transfer over an unsteady stretching surface embedded in a porous medium.

The present study envisages the flow past an exponentially shrinking sheet placed at the bottom of fluid saturated porous medium taking variable thermal conductivity and radiation into account. It is expected that the present work would help peep into analogous real world systems.

#### 2. Mathematical Model

Let us consider the steady two-dimensional boundary layer flow of optically thick viscous Newtonian fluid and associated heat transfer over a permeable sheet placed at the bottom of a fluid saturated porous medium having permeability of specific form. A Cartesian coordinate system is chosen where the -axis is taken along the sheet and -axis is normal to it. The flow is caused by the sheet shrinking in an exponential fashion. A suction is applied normal to sheet to contain the vorticity. The fluid considered here is without phase change, optically dense, absorbing-emitting radiation but a nonscattering medium. The thermal conductivity of the fluid is assumed to vary linearly with temperature. The radiation flux in the energy equation is presumed to follow Rosseland approximation. The boundary layer equations for the considered setup are with the following boundary conditions where , are the velocity components along and directions, respectively, is the permeability, is the specific heat at constant pressure, is the kinematic viscosity, is the density, and , , and are the temperature, viscosity and thermal conductivity of the fluid, respectively. Further, is the characteristic length, is the variable temperature at the sheet, is the constant reference temperature, and is the constant free stream temperature. and are the shrinking velocity of the sheet and mass transfer velocity, respectively, where is the shrinking constant and is a constant (where corresponds to mass suction).

#### 3. Similarity Transformations

We now introduce the stream function as Thus (1) is identically satisfied and the similarity transformations can be written as On using (5) and (6) we obtain the expressions for velocity components in nondimensional form as

In order to obtain the similarity solutions, it is assumed that the permeability of the porous medium takes the following form where is reference permeability.

Following Chiam [42], thermal conductivity of the fluid is assumed to vary with temperature in a linear fashion as where is the thermal conductivity variation parameter. In general for fluids such as water and air, while for fluids such as lubrication oils. The radiative heat flux in the energy equation is presumed to follow Rosseland approximation and is given by where is the Stephan-Boltzmann constant and is the mean absorption constant. It is further assumed that the temperature difference within the fluid is sufficiently small so that may be expressed as a linear function of temperature . This is done by expanding in a Taylor series about and omitting higher-order terms to yield The momentum (2) and the energy (3) thus reduce to the following nondimensional form with the boundary conditions where denote the permeability parameter, Prandtl number, and radiation parameter, respectively.

#### 4. Numerical Solution

The coupled and highly nonlinear boundary value problem (BVP) described by (12) and (13) is not amenable to analytic closed form solution, therefore we resorted to numerical solution. The BVP has been solved by fourth-order Runge-Kutta scheme together with shooting method. The essence of shooting method to solve a BVP is to first convert the BVP into system of initial value problems and make suitable guesses for the unknown quantities such that the end conditions of the solution space are satisfied.

By making some transformations the BVP is reduced to the following system of initial value problems under the boundary conditions Here and are the initial guesses which are arbitrarily chosen and an iterative procedure is set to obtain solutions through Runge-Kutta fourth-order method. The computational procedure involved two challenges, firstly determination of , that is, maximum value of for which and at and secondly proper estimates for the unknown quantities and . It is not out of place to make remark that the “guesses” were made purely on hit and trial basis and their refinement was interpolated iteratively with the prescribed error tolerance. A grid independence study was also carried out to examine the effect of step size and the edge of the boundary layer on the solution in the quest for their optimization. The , that is, at infinity was so chosen that further changes in it do not produce changes in the values of and (constant till 10−7) vis a vis boundary conditions are satisfied. A step size of was found to be satisfactory for a convergence criterion of 10−7 in all cases.

#### 5. Results and Discussion

The numerical computations carried out for different sets of values of the parameters entering into the problem have been depicted through graphs and tables. The analysis, besides the effects of parameters on flow and heat transfer, also aims to compute the range of mass suction for parameter-values required for laminar flow on the shrinking sheet. Here we are inclined to recall the works of Miklavčič and Wang [20] and Fang and Zhang [24] who discussed flow due to linear shrinking sheet in clear Newtonian fluid. They reported that to facilitate the steady flow due to shrinking sheet mass suction parameter of strength greater or equal to 2 is required. Recently, Bhattacharyya [31] computed values of suction velocity for the flow due to exponentially shrinking sheet in clear fluid regime () and reported that laminar flow does not exist if the suction velocity is less than 2.266684. Here it is emphasized that due to exponential shrinking of the sheet, the vorticity grows rapidly in the boundary layer and tends to go beyond the boundary layer, hence an adequate mass suction is required for existence of the laminar flow. Further, here it is worth noting that in the present analysis the sheet is placed at the bottom of porous medium which offers resistance to fluid traversal inside the porous matrix and contributes in some way in fluid stability, hence one may expect rather low suction rate as compared to the case of clear fluid situation. In present investigation, we did rigorous computations to extract the minimum suction rate required for ensuring the laminar flow in porous medium due to exponentially shrinking sheet. We found that in case of the porous medium one requires lesser suction values ( when ) to contain the vorticity in the boundary layer depending upon the porous material chosen (see Table 1). The table shows that similarity solution exists for the exponential shrinking sheet placed at the bottom of the fluid saturated porous medium if when ; when . Further it is emphasised that when , then the present problem as far as momentum equation is concerned reduces to that of Bhattacharyya [31]. We see that when , then it is found that . This is in perfect agreement with Bhattacharyya [31]. Further, Miklavčič and Wang [20] observed that in flow problems involving shrinking surfaces, one may encounter nonunique solutions. In the present analysis we get dual solution regime. These solutions are closely dependent on mass suction velocity at the surface.

Table 1: Variation of minimum suction velocity required for flow phenomena with the permeability parameter .

In order to solve the system of initial value problems we required proper estimates for the unknown quantities and . A large number of guesses for were tested for different values of permeability parameter . These guesses were chosen on purely hit and trial basis such that the end conditions got satisfied. An error tolerance of magnitude 10−7 was prescribed for the computations. Tables 2, 3, 4, and 5 display rigorously computed values of for dual solution regime for momentum equation and the range of mass suction for various set of values of permeability parameter.

Table 2: Estimates for guesses of for both the solutions of the velocity profile when K = 10.
Table 3: Estimates for guesses of for both the solutions of the velocity profile when K = 20.
Table 4: Estimates for guesses of for both the solutions of the velocity profile when K = 100.
Table 5: Estimates for guesses of for both the solutions of the velocity profile when K = 1000.

The profiles for the skin friction coefficient and rate of heat transfer versus suction velocity parameter drawn for various values of permeability parameter are shown in Figures 1 and 2, respectively. We find dual solution regimes for both and for the different values of suction velocity parameter. Further Figure 1 depicts that the skin friction coefficient registers increase with an increase in suction velocity for the first solution whereas it decays in the case of second solution for increasing values of . The figure also shows that with increasing values of permeability parameter , for the first solution decays whereas for the second solution it increases. From Figure 2 we see that increases for the first solution with increase in whereas for the second solution it decreases. Moreover with the increasing values of the permeability parameter , increases for the second solution but insignificant changes occur in case of first solution.

Figure 1: Profiles for skin friction coefficient () for variation in permeability parameter .
Figure 2: Profiles for temperature gradient at the sheet () for variation in permeability parameter .

Figure 3 shows the variation in with respect to the suction velocity parameter . For the first solution decreases numerically within the boundary layer for increasing values of suction parameter and becomes zero at the edge of the momentum boundary layer for every value . The profiles of second solution exhibit a similar pattern.

Figure 3: Dual velocity () profiles for variation in suction velocity .

The variation in for different values of permeability parameter , when the suction parameter is set as is shown in Figure 4. We see that for the first solution increases numerically with the increasing values of whereas decreases in the case of second solution for the increasing values of . The figure also reveals that effect of permeability is more significant on the second solution as compared to the first solution.

Figure 4: Dual velocity () profiles for variation in permeability parameter .

Variations in shear stress in the boundary layer region for the suction parameter values are shown in Figure 5. For the first solution, shear stress in the vicinity of the wall increases with the increasing values of suction parameter but at some critical distance from the wall, decays with the increasing values of . However, opposite trends in are observed in the case of second solution for increasing values of .

Figure 5: Shear stress () profiles for variation in suction velocity .

Figure 6 depicts the effect of permeability parameter on shear stress in the boundary layer. We observe that for the first solution, decreases near the wall for the increasing values of . For the second solution, increases with increasing near the wall. However, variations in are reversed for variation in in both the solutions at a certain distance from the wall. The figure also reveals that effect of on is more pronounced for the second solution as compared to that of the first one.

Figure 6: Shear stress () profiles for variation in permeability parameter .

Figure 7 exhibits temperature profiles for various suction values. It is clear from the figure that the temperature for the first solution decays with the increasing values of suction parameter whereas in the second solution case we find that increases with the increasing values of . We also observe that for the same set of suction values, the thermal boundary layer is thicker for the second solution compared to that of first solution.

Figure 7: Temperature () profiles for variation in suction velocity .

Figure 8 displays variation in temperature for varying values of permeability parameter . We see that the effect of permeability on is insignificant for the first solution whereas permeability has pronounced effect on temperature for second solution.

Figure 8: Dual temperature () profiles for variation in permeability parameter .

Figure 9 shows the variation of temperature for different values of the Prandtl number Pr. We see that for both solutions, decays with the increasing Pr values. However it is interesting to note that thermal boundary layer for the first solution is thicker as compared to that of second solution for all Pr values under consideration.

Figure 9: Dual temperature () profiles for variation in Prandtl number Pr.

Prandtl number physically characterizes relative importance of momentum diffusion to thermal diffusion in the flow field. Low density fluids such as liquid metals (e.g., mercury, ) enjoy low Prandtl number values, that is, less than unity. In thermal regime with low Prandtl number values heat diffuses at faster rate. On the contrary, oil having low thermal conductivity (high Pr values) gives rise to thinner boundary layer.

The effect of the radiation parameter on is shown in Figure 10. The radiation parameter is reciprocal of Stark number (also known as Stephan number) which measures the relative importance of thermal radiation transfer to the conduction heat transfer. Thus larger values of mean dominance of thermal radiation over conduction. Consequently, larger values of indicate that larger amount of radiative heat energy is poured into the flow field, causing rise in . We conclude from the figure that increase in radiation parameter results in an increase in fluid temperature.

Figure 10: Dual temperature () profiles for variation in radiation parameter .

Figure 11 displays the effect of thermal conductivity perturbation parameter on temperature. It is observed that increases with an increase in parameter .

Figure 11: Dual temperature () profiles for variation in the parameter .

#### Acknowledgment

The authors are thankful to the learned referees for offering their valuable suggestions to improve the paper.

#### References

1. L. J. Crane, “Flow past a stretching plate,” Zeitschrift für Angewandte Mathematik und Physik ZAMP, vol. 21, no. 4, pp. 645–647, 1970.
2. S. J. Liao, “On the analytic solution of magnetohydrodynamic flows of non-Newtonian fluids over a stretching sheet,” Journal of Fluid Mechanics, no. 488, pp. 189–212, 2003.
3. W. A. Khan and I. Pop, “Boundary-layer flow of a nanofluid past a stretching sheet,” International Journal of Heat and Mass Transfer, vol. 53, no. 11-12, pp. 2477–2483, 2010.
4. A. Ishak, R. Nazar, and I. Pop, “Mixed convection stagnation point flow of a micropolar fluid towards a stretching sheet,” Meccanica, vol. 43, no. 4, pp. 411–418, 2008.
5. K. R. Rajgopal, T. Y. Na, and A. S. Gupta, “Flow of a visco-elastic fluid over a stretching sheet,” RheolActa, vol. 23, pp. 213–215, 1984.
6. P. H. Veena, S. Abel, K. Rajagopal, and V. K. Pravin, “Heat transfer in a visco-elastic fluid past a stretching sheet with viscous dissipation and internal heat generation,” Zeitschrift fur Angewandte Mathematik und Physik, vol. 57, no. 3, pp. 447–463, 2006.
7. K. Ahmad, R. Nazar, A. Ishak, and I. Pop, “Unsteady three-dimensional boundary layer flow due to a stretching surface in a micropolar fluid,” International Journal for Numerical Methods in Fluids, vol. 28, pp. 118–122, 2011.
8. H. I. Andersson and O. A. Valnes, “Flow of a heated ferrofluid over a stretching sheet in the presence of a magnetic dipole,” Acta Mechanica, vol. 128, no. 1-2, pp. 39–47, 1998.
9. D. S. Chauhan and R. Agrawal, “MHD flow through a porous medium adjacent to a stretching sheet: numerical and an approximate solution,” European Physical Journal Plus, vol. 126, no. 5, pp. 1–14, 2011.
10. S. Liao, “A new branch of solutions of boundary-layer flows over an impermeable stretched plate,” International Journal of Heat and Mass Transfer, vol. 48, no. 12, pp. 2529–2539, 2005.
11. S. J. Liao, “A new branch of solutions of boundary-layer flows over a permeable stretching plate,” International Journal of Non-Linear Mechanics, vol. 42, no. 6, pp. 819–830, 2007.
12. E. Magyari and B. Keller, “Heat and mass transfer in the boundary layers on an exponentially stretching continuous surface,” Journal of Physics D, vol. 32, no. 5, pp. 577–585, 1999.
13. E. M. A. Elbashbeshy, “Heat transfer over an exponentially stretching continuous surface with suction,” Archives of Mechanics, vol. 53, no. 6, pp. 643–651, 2001.
14. R. A. M. Q. Al-Odat, T. A. Damesh, and T. A. Al-Azab, “Thermal boundary layer on an exponentially stretching continuous surface in the presence of magnetic field effect,” International Journal of Applied Mechanics and Engineering, vol. 11, pp. 289–299, 2006.
15. M. K. Partha, P. V. S. N. Murthy, and G. P. Rajasekhar, “Effect of viscous dissipation on the mixed convection heat transfer from an exponentially stretching surface,” Heat and Mass Transfer, vol. 41, no. 4, pp. 360–366, 2005.
16. E. Sanjayanand and S. K. Khan, “On heat and mass transfer in a viscoelastic boundary layer flow over an exponentially stretching sheet,” International Journal of Thermal Sciences, vol. 45, no. 8, pp. 819–828, 2006.
17. S. K. Khan and E. Sanjayanand, “Viscoelastic boundary layer flow and heat transfer over an exponential stretching sheet,” International Journal of Heat and Mass Transfer, vol. 48, no. 8, pp. 1534–1542, 2005.
18. M. Sajid and T. Hayat, “Influence of thermal radiation on the boundary layer flow due to an exponentially stretching sheet,” International Communications in Heat and Mass Transfer, vol. 35, no. 3, pp. 347–356, 2008.
19. C. Y. Wang, “Liquid film on an unsteady stretching sheet, Quart,” Applied Mathematics, vol. 48, pp. 601–610, 1990.
20. M. Miklavčič and C. Y. Wang, “Viscous flow due to a shrinking sheet,” Quarterly of Applied Mathematics, vol. 64, no. 2, pp. 283–290, 2006.
21. C. Y. Wang, “Stagnation flow towards a shrinking sheet,” International Journal of Non-Linear Mechanics, vol. 43, no. 5, pp. 377–382, 2008.
22. T. Hayat, Z. Abbas, and N. Ali, “MHD flow and mass transfer of a upper-convected Maxwell fluid past a porous shrinking sheet with chemical reaction species,” Physics Letters, Section A, vol. 372, no. 26, pp. 4698–4704, 2008.
23. T. Fang, “Boundary layer flow over a shrinking sheet with power-law velocity,” International Journal of Heat and Mass Transfer, vol. 51, no. 25-26, pp. 5838–5843, 2008.
24. T. Fang and J. Zhang, “Closed-form exact solutions of MHD viscous flow over a shrinking sheet,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 7, pp. 2853–2857, 2009.
25. T. G. Fang, J. Zhang, and S. S. Yao, “Viscous flow over an unsteady shrinking sheet with mass transfer,” Chinese Physics Letters, vol. 26, no. 1, Article ID 014703, 2009.
26. N. F. M. Noor, S. Awang Kechil, and I. Hashim, “Simple non-perturbative solution for MHD viscous flow due to a shrinking sheet,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 2, pp. 144–148, 2010.
27. T. Hayat, Z. Abbas, and M. Sajid, “On the analytic solution of magnetohydrodynamic flow of a second grade fluid over a shrinking sheet,” Journal of Applied Mechanics, Transactions ASME, vol. 74, no. 6, pp. 1165–1171, 2007.
28. T. Fang, W. Liang, and C. F. F. Lee, “A new solution branch for the Blasius equation-A shrinking sheet problem,” Computers and Mathematics with Applications, vol. 56, no. 12, pp. 3088–3095, 2008.
29. N. F. Mohd and I. Hashim, “MHD flow and heat transfer adjacent to a permeable shrinking sheet embedded in a porous medium,” Sains Malaysiana, vol. 38, no. 4, pp. 559–565, 2009.
30. D. S. Chauhan and R. Agrawal, “MHD flow and heat transfer in a channel bounded by a shrinking sheet and a plate with a porous substrate,” Journal of Engineering Physics and Thermophysics, vol. 84, no. 5, pp. 1034–1046, 2011.
31. K. Bhattacharyya, “Boundary layer flow and heat transfer over an exponentially shrinking sheet,” Chinese Physics Letters, vol. 28, no. 7, Article ID 074701, 2011.
32. B. S. Dandapat, B. Santra, and K. Vajravelu, “The effects of variable fluid properties and thermocapillarity on the flow of a thin film on an unsteady stretching sheet,” International Journal of Heat and Mass Transfer, vol. 50, no. 5-6, pp. 991–996, 2007.
33. P. Vyas and A. Rai, “Radiative flow with variable thermal conductivity over a noniIsothermal Stretching sheet in a porous medium, Int,” Journal of Contemporary Mathematical Sciences, vol. 5, pp. 2685–2698, 2010.
34. U. Sarma and G. C. Hazarika, “Effects of variable viscosity and thermal conductivity on heat and mass transfer flow along a vertical plate in the presence of a magnetic field,” Latin-American Journal of Physics Education, vol. 5, pp. 100–106, 2011.
35. Y. Khan, Q. Wu, N. Faraz, and A. Yildirim, “The effects of variable viscosity and thermal conductivity on a thin film flow over a shrinking/stretching sheet,” Computers and Mathematics with Applications, vol. 61, no. 11, pp. 3391–3399, 2011.
36. P. Vyas and N. Srivastava, “Radiative MHD flow over a non-isothermal stretching sheet in a porous medium,” Applied Mathematical Sciences, vol. 4, no. 49–52, pp. 2475–2484, 2010.
37. P. Vyas and A. Ranjan, “Dissipative MHD boundary-layer flow in a porous medium over a sheet stretching nonlinearly in the presence of radiation,” Applied Mathematical Sciences, vol. 4, no. 61–64, pp. 3133–3142, 2010.
38. D. S. Chauhan and V. Kumar, “Radiation effects on unsteady flow through a porous medium channel with velocity and temperature slip boundary conditions,” Applied Mathematical Sciences, vol. 6, no. 33–36, pp. 1759–1769, 2012.
39. D. S. Chauhan and R. Agarwal, “MHD coupled flow and heat transfer across a porous layer due to an oscillating plate with radiation,” Afrika Mathematica. In press.
40. T. Hayat, S. A. Shehzad, M. Qasim, and A. Alsaedi, “Radiative flow with variable thermal conductivity in porous medium,” Zeitschrift für Naturforschung, vol. 67, no. 3-4, pp. 153–159, 2012.
41. E. M. A. Elbashbeshy and T. G. Emam, “Effects of thermal radiation and heat transfer over an unsteady stretching surface embedded in a porous medium in the presence of heat source or sink,” Thermal Science, vol. 15, no. 2, pp. 477–485, 2011.
42. T. C. Chiam, “Heat transfer in a fluid with variable thermal conductivity over a linearly stretching sheet,” Acta Mechanica, vol. 129, no. 1-2, pp. 63–72, 1998.