## Thermal Radiation Effect on a Porous Media under Optically Thick Approximation

**M.O Ibrahim**^{1}, **K.K. Asogwa**^{2,}, **I.J Uwanta**^{3}, **B.G Dan Shehu**^{4}

^{1}Department of Mathematics, University of Ilorin, Kwara State, Nigeria

^{2}Department of Physical and Computer Sciences, College of Natural and Applied Science McPherson University, Abeokuta, Ogun State Nigeria

^{3}Department of Mathematics, Usmanu Danfodiyo University, Sokoto, Nigeria

^{4}Sokoto Energy Research Centre, Usmanu Danfodiyo University, Sokoto, Nigeria

### Abstract

The purpose of this paper is to investigate thermal radiation effect on a porous media under optically thick approximation using Newton scheme method from Taylor series implemented numerically on MATLAB. The velocity profiles and temperature profile are studied for different physical parameters like, Porous term P, Radiation F and thermal Grashof number Ga.** **The results shows that Porous parameter increases with increasing velocity, while the trend reverses with thermal radiation and thermal Grashof number under optically thick approximation. The flow rate increases asymmetrically due to conduction.

### At a glance: Figures

**Keywords:** porous media, recurrence relation, radiation, flow rate

*American Journal of Numerical Analysis*, 2013 1 (1),
pp 15-21.

DOI: 10.12691/ajna-1-1-3

Received September 16, 2013; Revised November 05, 2013; Accepted November 25, 2013

**Copyright:**© 2013 Science and Education Publishing. All Rights Reserved.

### Cite this article:

- Ibrahim, M.O, et al. "Thermal Radiation Effect on a Porous Media under Optically Thick Approximation."
*American Journal of Numerical Analysis*1.1 (2013): 15-21.

- Ibrahim, M. , Asogwa, K. , Uwanta, I. , & Shehu, B. D. (2013). Thermal Radiation Effect on a Porous Media under Optically Thick Approximation.
*American Journal of Numerical Analysis*,*1*(1), 15-21.

- Ibrahim, M.O, K.K. Asogwa, I.J Uwanta, and B.G Dan Shehu. "Thermal Radiation Effect on a Porous Media under Optically Thick Approximation."
*American Journal of Numerical Analysis*1, no. 1 (2013): 15-21.

Import into BibTeX | Import into EndNote | Import into RefMan | Import into RefWorks |

### 1. Introduction

Fluid flow through a porous channel has been studied theoretically and experimentally by numerous authors due to its wide applications in various fields such as diffusion technology, transpiration cooling, hemodialysis processes, flow control in nuclear reactors, etc. In a pioneering work, Berman ^{[1]} studied an exact solution of the Navier-Stokes equations that describes the steady two-dimensional flow of an incompressible viscous fluid along a channel with parallel rigid porous walls, the flow being driven by uniform suction or injection at the walls. henceforth, many authors such as Terrill ^{[2]}, Makinde ^{[3]}, Robinson ^{[4]} and Brady ^{[5]} have extended and reconsidered the problem under various physical conditions. Minkowycz and Cheng ^{[6]}, Cheng and Minkowycz ^{[7]}, and Badr and Pop ^{[8]} were based on Darcy’s law which states that the volume-averaged velocity is proportion to the pressure gradient. Kaviany ^{[9]} used the line integral method to study the heat transfer from a semi-infinite flat plate embedded in a fluid-saturated porous medium. Jang and Shiang ^{[10]} studied the mixed convection along a vertical adiabatic surface embedded in a porous medium. Few studies of convective boundary-layer flows in porous media using the Darcy-Brinkman equation model are considered for the momentum equation, for example, Hsu and Cheng ^{[11]}, Rees and Vafai ^{[12]}, Nazar et al. ^{[13, 14]}, Ishak et al. ^{[15]}, and Harris et al ^{[16]}.

Sreekanth *et al*. ^{[17]} examined the effect of permeability variation on free convective flow past a vertical porous wall in a porous medium when the permeability varies with respect to time. Subsequently, Singh *et al*. ^{[18]} studied hydro magnetic free convective and mass transfer flow of a viscous stratified fluid considering the variation in permeability with direction. The problem of unsteady oscillatory laminar free convection flow of an electrically conducting fluid through a porous medium along a porous hot vertical plate with time dependent suction in the presence of heat source or sink was analyzed by Kumar *et al. *^{[19]} while Acharya *et al*. ^{[20]} discussed the magnetic field effects on the free convention and mass transfer flow through a porous medium with constant suction and constant heat flux. An experimental and analytic investigation was carried out by Fand et al. ^{[21]} to study the free convection heat transfer from a horizontal cylinder embedded in porous medium consisting of randomly packed glass spheres and the medium is saturated by water or silicon oil. Murthy and Singh ^{[22]} have modified the viscous dissipation effect on the flow of an incompressible fluid in a saturated porous medium. They applied the Forchheimer-Darcy model for momentum equation. Neild ^{[23]} suggested a modified formula to model viscous dissipation in non-Darcy porous medium. He concluded that modelling of viscous dissipation is related with local drag modelling. Recently Taiwo and Ogunlaran ^{[24]} studied Numerical solution of fourth order linear ordinary differential equations by cubic spline collocation tau method.

Makinde ^{[25]} studied Hermite-Pad´e approximation approach to nonlinear wall driven steady flow in a tube, Recently Ibrahim *et al** *^{[27]} presented and analysed Radiative Effect on MHD Fluid Flow in a Vertical Channel under Optically Thick Approximation.

This research work is an improved work of Ibrahim ^{[28]} to include porous term. Hence it is proposed to study thermal radiation effect on porous media under optically thick approximation.

### 2. Formulation of the Problem

**Figure**

**1**

**.**Geometry of the problem

The geometry of the problem shows thermal radiation effect on a porous media under optically thick approximation formulated, analysed and solved numerically. The - axis is taken along the plate in the vertically upward direction and also the -axis is taken normal to the plate. That the flow is fully developed, the velocity and temperature fields are symmetrical about the central line of the channel in a porous media. The temperature of the walls is the same and is maintained at a constant temperature. The viscosity, the thermal conductivity and specific heat are independent of temperature and the essential influence of the variation in density is included in the body force term. Steady flow equations are momemtum and energy equation.

(1) |

(2) |

where u and are the velocity of the fluid, T is the fluid temperature, β is the thermal expansion of the fluid is the fluid density, is the specific heat capacity, v is the viscosity of the fluid, k is the thermal conductivity and the temperature of the wall , permeability of the porous medium parameter and the kinematic viscosity.

In this research work the mathematical formulation has Porous term which is not included in the work of Ibrahim (1997) where the thermal diffusivity and is the radiative, using the Rosseland differential approximation.

(3) |

The boundary conditions are;

(4) |

On introducing the following non-dimensional quantities

(5) |

Substituting the non-dimensional quantities of equation (5) into (1) to (2), leads to

(6) |

(7) |

Where

Equation (6) and (7) leads to

(8) |

(9) |

### 3. Solution to the Problem

To solve equations (8), subjected to the boundary conditions of (9), the solutions are obtained for temperature and velocity flow. using the Newton scheme method from Taylor series to solve the fourth order non linear problem.

(10) |

Thus from (8) we have

(11) |

(12) |

(13) |

(14) |

(15) |

Substituting equations (11) to (15) into equation (10) we obtain

(16) |

Expanding equations (16) where

(17) |

(18) |

Collecting the terms involving

On the L.H.S and the terms involving neglecting

On the R.H.S gives

(19) |

Where is our initial solution and is the assumed solution.

Let

Considering N = 7

(20) |

Putting (20) into (19)

(21) |

Simplifying by collecting terms in we obtain.

(22) |

We now collate equations (22) at point , where

In the problem under consideration, N=7, a =-1 and b=+1.

The boundary conditions are

Thus, we obtain five equations from (22).

Let be our initial solution

We now substitute for into (22) and then solve the equation simultaneously to obtain the unknowns using Gaussian elimination method with partial pivoting using MATLAB software to obtain the unknown constants

But

(23) |

Determination of the Temperature profile, to obtain the solution for the temperature

in the interval since the various values of and are known from the MatLab Program

Determination of the Velocity profile

The velocity distribution of the flow can be obtained from equation (6)

(24) |

(25) |

We have our assumed solution to be

Substituting the expression into (6) gives

(26) |

Evaluating equation (6) in the interval will enable us to determine the effect of different radiation parameter.

Determination of the Non-dimensional flow rate

The non-dimensional flow rate through the channel per unit width is given by

(27) |

Using the equation (6) and evaluating the integral will enable us to determine the flow rate.

Determination of Heat Transfer coefficient

Thus heat transfer coefficient due to thermal conduction is given by

(28) |

### 4. Results and Discussion

The problem of thermal radiation effect on a porous media under optically thick approximation formulated, analysed and solved numerically. In order to point out the effects of physical parameters namely; thermal Grashof number Ga, radiation parameter F, Porous term P. on the flow patterns, the computation of the flow fields are carried out.The values of velocity and temperature are obtained for the physical parameters as mention. The velocity profiles has been studied and presented in Figure 2 to Figure 4 The effect of velocity for different values of Porous parameter (P = 2, 5, 10, 20) is presented in Figure 2. The trend shows that the velocity increases with increasing Magnetic parameter. The effect of velocity for different values of radiation (F = 2, 5, 10) is also presented in Figure 3. It is then observed that the velocity decreases with increasing values of radiation. The effect of velocity for different values of thermal Grashof number (Ga= 1, 2, 3) is also presented in Figure 4. It is then observed that the velocity decreases with increasing values of thermal Grashof number

**Figure**

**2**

**.**

**Velocity profiles for different values of P**

**Figure**

**3**

**.**

**Velocity profiles for different values of F**

**Figure**

**4**

**.**velocity profiles for different values of Ga

**Figure**

**5**

**.**

**Trend for temperature profiles against y**

**Figure**

**6**

**.**Trend of heat transfer coefficient due to thermal conduction h

In Figure 5 It is observed that temperature rise with increasing values of while Figure 6 shows that trend of heat transfer coefficient due to thermal conduction h increases with the rate of change of at less than 5 cooling sets in.

### 5. Summary and Conclusion

Thermal radiation effect on a porous media under optically thick approximation has been studied. .In order to point out the effects of physical parameters namely; thermal Grashof number Ga, radiation parameter F, Porous parameter P. are presented graphically. It is observed that the velocity increases with increasing Porous parameter, while radiation and thermal Grashof number decreases with increasing values of radiation and thermal Grashof number respectively.

It also observed that temperature rise with increasing values of and trend of heat transfer coefficient due to thermal conduction h increases with the rate of change of at less than 5 cooling sets in.

### References

[1] | Berman, A. S.: “Laminar flow in channels with porous walls,” J. Appl. Phys., 24, 1232-1235. 1953. | ||

In article | CrossRef | ||

[2] | Terrill R. M.: “Laminar flow in a uniformly porous channel with large injection,” Aeronaut. Q. 1965. | ||

In article | |||

[3] | Robinson W. A.: “The existence of multiple solutions for the laminar flow in a uniformly porous channel with suction at both walls,” J. Eng. Math. 10, 23 1976. | ||

In article | CrossRef | ||

[4] | Makinde, O.D.: “Extending the utility of perturbation series in problems of laminarflow in a porous pipe and diverging channel,”.J. Austral. Math. Soc. Ser B 40, 1-11 1998. | ||

In article | |||

[5] | Brady, J. F.: “Flow development in a porous channel and tube” Phys. Fluids, 27, 1061-1067, 1984. | ||

In article | CrossRef | ||

[6] | Minkowycz W. J. and Cheng P.: “Free convection about a vertical cylinder embedded in a porous medium,” International Journal of Heat and Mass Transfer, vol. 19, no. 7, pp. 805-813, 1976. | ||

In article | CrossRef | ||

[7] | Cheng P. and Minkowycz, W. J.: “Free convection about a vertical flat plate embedded in a porous medium with application to heat transfer from a dike,” Journal of Geophysical Research, vol. 82, no. 14, pp. 2040-2044, 1977. | ||

In article | CrossRef | ||

[8] | Badr H. M. and Pop, I.: “Combined convection from an isothermal horizontal rod buried in a porous medium,” International Journal of Heat and Mass Transfer, vol. 31, no. 12, pp. 2527-2541, 1988. | ||

In article | CrossRef | ||

[9] | Kaviany M.: “boundary-layer treatment of forced convection heat transfer from a semi-infinite flat plate embedded in porous media,” Journal of Heat Transfer, vol. 109, no. 2, pp. 345-349, 1987. | ||

In article | CrossRef | ||

[10] | Jang J.Y. and Shiang C.T.: “The mixed convection plume along a vertical adiabatic surface embedded in a non-Darcian porous medium,” International Journal of Heat and Mass Transfer, vol. 40, no. 7, pp. 1693-1699, 1997. | ||

In article | CrossRef | ||

[11] | C. T. Hsu and P. Cheng: “The Brinkman model for natural convection about a semi-infinite vertical flat plate in a porous medium,” International Journal of Heat and Mass Transfer, vol. 28, no. 3, pp. 683-697, 1985. | ||

In article | CrossRef | ||

[12] | Rees D. A. and Vafai K.: “Darcy-Brinkman free convection from a heated horizontal surface,” Numerical Heat Transfer, vol. 35, no. 2, pp. 191-204, 1999. | ||

In article | CrossRef | ||

[13] | Nazar R., Amin N., Filip D., and Pop I..: “The Brinkman model for the mixed convection boundary layer flow past a horizontal circular cylinder in a porous medium,” International Journal of Heat and Mass Transfer, vol. 46, no. 17, pp. 3167-3178, 2003. | ||

In article | CrossRef | ||

[14] | Nazar R., Amin N., Pop and I.: “Unsteady mixed convection boundary layer flow near the stagnation point on a vertical surface in a porous medium,” International Journal of Heat andMass Transfer, vol. 47, no. 12-13, pp. 2681-2688, 2004. | ||

In article | |||

[15] | Ishak A., Nazar R., and Pop I.: “Steady and unsteady boundary layers due to a stretching vertical sheet in a porousmedium using Darcy-Brinkman equation mode,” International Journal of AppliedMechanics and Engineering, vol. 11, pp. 623-637, 2006. | ||

In article | |||

[16] | Harris S. D., Ingham D. B., and Pop I.: “Mixed convection boundary-layer flow near the stagnation point on a vertical surface in a porous medium: brinkman model with slip,” Transport in PorousMedia, vol. 77, no. 2, pp. 267-285, 2009. | ||

In article | CrossRef | ||

[17] | Sreekanth S., Venkataramana S. and Ramakrishna S.: “Effect of the permeability variation on free convective flow past a vertical porous wall in porous medium when the permeability varies with time,” Acta Ciencia Indica, 22M, Pp 267-275 1996. | ||

In article | |||

[18] | Singh A. K. and Singh N. P.: “Hydromagnetic free convective and mass transfer flow of a viscous stratified fluid considering the variation in permeability with direction,” Indian Journal of Pure and Applied Mathematics, 34, Pp 429-442 2003. | ||

In article | |||

[19] | Kumar A., Chand B. and Kaushik: “Unsteady oscillatory laminar free convection flow of an electrically conducting fluid through a porous medium along a porous hot vertical plate with time dependent suction,” Journal of Academy of Mathematics, 24 (1), Pp 339-354 2002. | ||

In article | |||

[20] | Acharya M., Dash C. G. and Singh L. P.: “Magnetic field effects on the free convection and mass transfer flow through a porous medium with constant suction and constant heat flux,” Indian journal of Pure and Applied Mathematics, 31 (1), Pp 1-18 2000. | ||

In article | |||

[21] | Fand, R. M., Steinberger, T. E. and Cheng, P.: “Natural convection heat transfer from a horizontal cylinder embedded in a porous medium,” Int. J. Heat Mass Transfer, Vol. 29, pp. 119 1986. | ||

In article | CrossRef | ||

[22] | Murthy, P. V. S. N. and Singh, P. “Effect of viscous dissipation on a non- Darcy natural convection regime,” Int. J. Heat Mass Transfer, Vol. 40, pp. 1251-1260 1997. | ||

In article | CrossRef | ||

[23] | Nield, D. A.: “Resolution of a paradox involving viscous dissipation and non-linear drag in a porous medium,” Transport in porous media, Vol. 41, pp. 349-357. 2000. | ||

In article | CrossRef | ||

[24] | Taiwo O.A. and Ogunlaran O.M.: ‘‘Numerical solution of fourth order linear ordinary differential equations by cubic spline collocation tau method,’’ Journal of Mathematics and Statistics 4 (4): 264-268, 2008. | ||

In article | CrossRef | ||

[25] | Makinde O. D.: “Hermite–Pad´e approximation approach to nonlinear wall driven steady flow in a tube,” Nonlinear Studies Vol. 13, No. 3, pp. 175-184, 2006 | ||

In article | |||

[26] | Chinyoka T. and Makinde O.D.: “Analysis of entropy generation rate in an unsteady porous channel flow with Navier Slip and convective cooling,” entropy 15, 2081-2099; 2013. | ||

In article | |||

[27] | M.O Ibrahim, K.K.Asogwa, I.J Uwanta, and B.G Dan Shehu, “Radiative Effect on MHD Fluid Flow in a Vertical Channel under Optically Thick Approximation.” American Journal of Modeling and Optimization 1, no. 3 (2013): 36-42. | ||

In article | |||

[28] | Ibrahim M.O: ‘‘Radiation effect on laminar convective flow of a radiating gas in a vertical channel under optically thick limit approximation,’’ M.Sc thesis , University of Ilorin, Nigeria. 1997. | ||

In article | |||