**Abstract**

In this chapter, the magnetohydrodynamic effects on heat transfer and thermal radiation at a stagnation point flowing in a nanofluid containing different types of nanoparticles namely, copper (Cu), alumina (Al2O3) and titania (TiO2) through a porous medium have been investigated numerically. By using appropriate transformation for velocity and temperature into a set of non-linear coupled ordinary differential equations which are solved numerically. Numerical results are presented for velocity and temperature profiles for different parameters of the problem. Also, the effects of the pertinent parameters on the skin friction and the heat fluxes are obtained and discussed numerically and illustrated graphically.

**Keywords:** nanofluid, MHD, porous medium, thermal radiation, heat transfer

#### **1. Introduction**

Nanofluids have novel properties that make them potentially useful in many applications. They exhibit enhanced thermal conductivity and the convective heat transfer coefficient compared to the base fluid [1], so the nanofluids transfer heat at a higher rate than ordinary fluids (for example, water) which allows for more efficient heating or cooling while reducing energy consumption. Since nanofluid consists of very small sized solid particles, therefore in low solid concentration it is reasonable to consider nanofluid as a single-phase flow [2].

The term "nanofluid" refers to a liquid containing a suspension of submicronic solid particles (nanoparticles). This interest is generated by a variety of applications, ranging from laser-assisted drug delivery to electronic chip cooling. The term was coined by Choi [3]. The characteristic feature of nanofluids is thermal conductivity enhancement, a phenomenon observed by Masuda et al. [4]. This phenomenon suggests the possibility of using nanofluids in advanced nuclear system [5].

In recent years, numerous investigations have been conducted on the magnetohydrodynamic (MHD) flows and heat transfer because of its important applications in metallurgical industry, such as the cooling of continuous strips and filaments, drawn through a quiescent fluid and the purification of molten metals from non-metallic inclusions. It is known that the properties of the final product depend considerably on the rate of cooling during the manufacturing processes. The rate of cooling can be

controlled by drawing the strips in an electrically conducting fluid subject to a magnetic field, so that a final product of desired characteristics can be achieved [6]. One can solve this problem by using solid particles as an additive suspended into the base fluid as Choi [3] did, who first utilized and used fluids suspended by nanometer-sized solid particles. The resulting mixture referred to as a nanofluid possesses a substantially larger thermal conductivity compared to that of the traditional fluids [7]. These suspended nanoparticles can change the transport and thermal properties of the base fluid [8]. Therefore, by mixing the nanoparticles in the fluid, thermal conductivity of the fluid increases and the heat transfer capability improves.

where *J* is the electric current density, *B* ¼ *B* þ *d* is the total magnetic field, *μ<sup>f</sup>* is the magnetic permeability and *d* is the induced magnetic field. The magnetic Reynolds number of the flow is taken to be small, so that the flow induction distortion of the applied magnetic field can be neglected as in the case with most of conducting fluids. The magnetic body force *<sup>J</sup>* � *<sup>B</sup>* takes the form *<sup>σ</sup> <sup>V</sup>* � *<sup>B</sup>* � *<sup>B</sup>*,

*Heat Transfer and Thermal Radiation at a General Three-Dimensional in a Nanofluid…*

*u, Fy* ¼ �*σB*<sup>2</sup>

*∂*2 *u <sup>∂</sup>z*<sup>2</sup> � *<sup>σ</sup> <sup>B</sup>*<sup>2</sup> *ρnf*

*∂*2 *v <sup>∂</sup>z*<sup>2</sup> � *<sup>σ</sup> <sup>B</sup>*<sup>2</sup> *ρnf*

*∂*2 *T <sup>∂</sup>z*<sup>2</sup> � <sup>1</sup> ð Þ *ρc nf*

*<sup>u</sup>* ! *ue*ð Þ *<sup>x</sup> , v* ! *ve*ð Þ*<sup>y</sup> , T* ! *<sup>T</sup>*<sup>∞</sup> *as z* ! <sup>∞</sup> *:* (6)

We now make the standard boundary-layer approximation, based on a scale

*x* þ *μnf ρnf*

*x* þ *μnf ρnf*

*u* ¼ *v* ¼ *w* ¼ 0*, T* ¼ *Tw at z* ¼ 0

nanofluid and *ρnf* the density of the nanofluid, which are given by

<sup>5</sup>*=*<sup>2</sup> *, <sup>α</sup>nf* <sup>¼</sup> *knf*

Here *μnf* the viscosity of the nanofluid, *αnf* the thermal diffusivity of the

*ρCp nf*

where *qr* is the energy flux which is relative to a frame moving with the nanofluid velocity. Using the Rosseland approximation (Rashed [16]), the radiative

*qr* ¼ �*k*<sup>1</sup>

Assuming that the temperature difference within the flow is sufficiently small

*∂T*<sup>4</sup> *∂ z*

<sup>∞</sup>*<sup>T</sup>* � <sup>3</sup>*T*<sup>4</sup>

*<sup>f</sup>* þ *φ ρCp s ,*

*B* ¼ ð Þ 0*;* 0*; B* . The Lorentz force (MHD body force) has two components:

*Fx* ¼ �*σB*<sup>2</sup>

*∂u ∂x* þ *∂ v ∂y* þ *∂w*

analysis, and write the governing equations (see [6, 11]):

*V*, where *V* is velocity vector *V* ¼ ð Þ *u; v; w* and

*<sup>∂</sup><sup>z</sup>* <sup>¼</sup> <sup>0</sup>*,* (2)

*<sup>u</sup>* � *<sup>μ</sup>nf ρnf κ*

*<sup>v</sup>* � *<sup>μ</sup>nf ρnf κ*

*∂qr*

*, ρnf* ¼ ð Þ 1 � *φ ρ<sup>f</sup>* þ *φ ρs,*

<sup>∞</sup> (9)

*u,* (3)

*v,* (4)

(7)

(8)

*<sup>∂</sup> <sup>z</sup> ,* (5)

*v:*

therefore, *<sup>σ</sup> <sup>V</sup>* � *<sup>B</sup>* � *<sup>B</sup>* ¼ �*σB*<sup>2</sup>

*DOI: http://dx.doi.org/10.5772/intechopen.81575*

*u ∂u ∂x* þ *v ∂u ∂y* þ *w ∂u <sup>∂</sup><sup>z</sup>* <sup>¼</sup> *<sup>a</sup>*<sup>2</sup>

*u ∂ v ∂x* þ *v ∂ v ∂y* þ *w ∂ v <sup>∂</sup><sup>z</sup>* <sup>¼</sup> *<sup>b</sup>*<sup>2</sup>

> *u ∂T ∂x* þ *v ∂T ∂y* þ *w ∂T <sup>∂</sup><sup>z</sup>* <sup>¼</sup> *<sup>α</sup>nf*

Subject to the boundary conditions

*<sup>μ</sup>nf* <sup>¼</sup> *<sup>μ</sup><sup>f</sup>*

*ρCp* 

heat flux *qr* could be expressed by

**143**

*knf kf*

ð Þ 1 � *φ*

<sup>¼</sup> *ks* <sup>þ</sup> <sup>2</sup>*kf*

*ks* þ 2*kf*

*nf* ¼ ð Þ 1 � *φ ρCp*

� <sup>2</sup>*<sup>φ</sup> kf* � *ks*

<sup>þ</sup> *<sup>φ</sup> kf* � *ks*

where *<sup>k</sup>*<sup>1</sup> <sup>¼</sup> <sup>4</sup>*σ*<sup>∗</sup>*=*3*k*<sup>∗</sup> is the nanofluid thermal conductivity.

such that *T*<sup>4</sup> could be approached as the linear function of temperature;

*<sup>T</sup>*<sup>4</sup> ffi <sup>4</sup>*T*<sup>3</sup>

The flow and heat transfer characteristics at a stagnation point, for both twodimensional and axisymmetric, have been studied extensively in the literature [9]. These studies have been motivated by the fundamental nature of the boundary layer flows near such points by the exact applicability of similar solutions and by their relevance to the leading edge and nose regions of bodies in high speed flows. These two cases of two-dimensional and axisymmetric flows can be recognized as special cases of more general stagnation point flows. Both two-dimensional and axisymmetric flows were extended to three-dimensional by Howarth [10]. Bhattacharyya and Gupta [11] studied the flow and heat transfer in an incompressible viscous and electrically conducting fluid near a three-dimensional stagnation point of a body permeated by a uniform magnetic field. Analysis of such flows is very important in both theory and in practice. From a theoretical point of view, flows of this type are fundamental in fluid mechanics and forced convective heat transfer. On the other hand, from a practical point of view, these flows have applications in many manufacturing processes in industry such as the boundary layer along material handling conveyers, the aerodynamic extrusion of plastic sheet and petrochemical industries.

It is worth mentioning that the nanofluid model proposed by Buongiorno [12] was very recently used by Nield and Kuznetsov [13], Kuznetsov and Nield [14] and Bachok et al. [15] in their papers. Different from the above model.

The aim of the present chapter is to study the Magnetohydrodynamic effects on heat transfer and thermal radiation in an incompressible viscous fluid near the three-dimensional stagnation point of a body that is placed in a water based nanofluid through a porous medium containing different types of nanoparticles: copper (Cu), alumina (Al2O3) and titania (TiO2). Numerical results are presented for velocity and temperature profiles for different parameters of the problem.

#### **2. Mathematical analysis**

Consider a flow of an electrically conducting fluid with heat transfer flow at a stagnation point of an incompressible viscous fluid past a body that is placed in a nanofluid of ambient uniform temperature *T*<sup>∞</sup> through a porous medium in the presence of radiation has been considered, where the body surface is kept at a constant temperature *Tw*. The stagnation point is located at the origin *o* of the Cartesian coordinate system *oxyz*. It was shown by Howarth [10] that the inviscid irrational flow near *o* has the velocity components

$$u\_{\epsilon}(\mathbf{x}) = a\mathbf{x}, \quad v\_{\epsilon}(\mathbf{y}) = b\mathbf{y} \tag{1}$$

with the constants *a* > 0 and *b* ≥ 0 or *b* ≤ 0. The MHD body forces *J* � *B* the Maxwell's equations:

$$\operatorname{div} B = 0, \quad \operatorname{Curl} B = \mu\_m I \text{ and } \quad \operatorname{div} E = 0.$$

*Heat Transfer and Thermal Radiation at a General Three-Dimensional in a Nanofluid… DOI: http://dx.doi.org/10.5772/intechopen.81575*

where *J* is the electric current density, *B* ¼ *B* þ *d* is the total magnetic field, *μ<sup>f</sup>* is the magnetic permeability and *d* is the induced magnetic field. The magnetic Reynolds number of the flow is taken to be small, so that the flow induction distortion of the applied magnetic field can be neglected as in the case with most of conducting fluids. The magnetic body force *<sup>J</sup>* � *<sup>B</sup>* takes the form *<sup>σ</sup> <sup>V</sup>* � *<sup>B</sup>* � *<sup>B</sup>*, therefore, *<sup>σ</sup> <sup>V</sup>* � *<sup>B</sup>* � *<sup>B</sup>* ¼ �*σB*<sup>2</sup> *V*, where *V* is velocity vector *V* ¼ ð Þ *u; v; w* and *B* ¼ ð Þ 0*;* 0*; B* . The Lorentz force (MHD body force) has two components:

$$F\_{\mathfrak{x}} = -\sigma B^2 \mathfrak{u}, \\ F\_{\mathfrak{y}} = -\sigma B^2 \mathfrak{v}.$$

We now make the standard boundary-layer approximation, based on a scale analysis, and write the governing equations (see [6, 11]):

$$
\frac{
\partial u
}{
\partial \mathbf{x}
} + \frac{
\partial v
}{
\partial y
} + \frac{
\partial w
}{
\partial \mathbf{z}
} = \mathbf{0},
\tag{2}
$$

$$
\mu \frac{\partial u}{\partial \mathbf{x}} + \nu \frac{\partial u}{\partial \mathbf{y}} + \nu \frac{\partial u}{\partial \mathbf{z}} = a^2 \mathbf{x} + \frac{\mu\_{\text{nf}}}{\rho\_{\text{nf}}} \frac{\partial^2 u}{\partial \mathbf{z}^2} - \frac{\sigma B^2}{\rho\_{\text{nf}}} u - \frac{\mu\_{\text{nf}}}{\rho\_{\text{nf}} \kappa} u,\tag{3}
$$

$$
\mu \frac{\partial v}{\partial \mathbf{x}} + v \frac{\partial v}{\partial \mathbf{y}} + w \frac{\partial v}{\partial \mathbf{z}} = b^2 \mathbf{x} + \frac{\mu\_{\text{nf}}}{\rho\_{\text{nf}}} \frac{\partial^2 v}{\partial \mathbf{z}^2} - \frac{\sigma B^2}{\rho\_{\text{nf}}} v - \frac{\mu\_{\text{nf}}}{\rho\_{\text{nf}} \kappa} v,\tag{4}
$$

$$u\frac{\partial T}{\partial \mathbf{x}} + v\frac{\partial T}{\partial \mathbf{y}} + w\frac{\partial T}{\partial \mathbf{z}} = a\_{\eta f}\frac{\partial^2 T}{\partial \mathbf{z}^2} - \frac{1}{(\rho c)\_{\eta f}}\frac{\partial q\_r}{\partial \mathbf{z}},\tag{5}$$

Subject to the boundary conditions

$$\begin{aligned} u = v = w = 0, &\ T = T\_w \quad \text{at} \quad z = \mathbf{0} \\ u \to u\_\epsilon(\mathbf{x}), &\ v \to v\_\epsilon(\mathbf{y}), &\ T \to T\_\infty \text{ as} \ z \to \infty \end{aligned} \tag{6}$$

Here *μnf* the viscosity of the nanofluid, *αnf* the thermal diffusivity of the nanofluid and *ρnf* the density of the nanofluid, which are given by

$$\begin{aligned} \mu\_{\mathfrak{n}f} &= \frac{\mu\_f}{(1-\varrho)^{5/2}}, \quad \alpha\_{\mathfrak{n}f} = \frac{k\_{\mathfrak{n}f}}{\left(\rho \, \mathbf{C}\_{\mathfrak{p}}\right)\_{\mathfrak{n}f}}, \quad \rho\_{\mathfrak{n}f} = (\mathbf{1}-\varrho)\rho\_f + \varrho\rho\_s, \\\ \left(\rho \, \mathbf{C}\_{\mathfrak{p}}\right)\_{\mathfrak{n}f} &= (\mathbf{1}-\varrho)\left(\rho \, \mathbf{C}\_{\mathfrak{p}}\right)\_f + \varrho\left(\rho \, \mathbf{C}\_{\mathfrak{p}}\right)\_s, \\\ k\_{\mathfrak{n}f} &= \frac{\left(k\_s + 2k\_f\right) - 2\varrho\left(k\_f - k\_s\right)}{\left(k\_s + 2k\_f\right) + \varrho\left(k\_f - k\_s\right)} \end{aligned} \tag{7}$$

where *qr* is the energy flux which is relative to a frame moving with the nanofluid velocity. Using the Rosseland approximation (Rashed [16]), the radiative heat flux *qr* could be expressed by

$$q\_r = -k\_1 \frac{\partial T^4}{\partial z} \tag{8}$$

where *<sup>k</sup>*<sup>1</sup> <sup>¼</sup> <sup>4</sup>*σ*<sup>∗</sup>*=*3*k*<sup>∗</sup> is the nanofluid thermal conductivity.

Assuming that the temperature difference within the flow is sufficiently small such that *T*<sup>4</sup> could be approached as the linear function of temperature;

$$T^4 \cong 4T\_{\infty}^3 T - \Im T\_{\infty}^4 \tag{9}$$

controlled by drawing the strips in an electrically conducting fluid subject to a magnetic field, so that a final product of desired characteristics can be achieved [6]. One can solve this problem by using solid particles as an additive suspended into the base fluid as Choi [3] did, who first utilized and used fluids suspended by nanometer-sized solid particles. The resulting mixture referred to as a nanofluid possesses a substantially larger thermal conductivity compared to that of the traditional fluids [7]. These suspended nanoparticles can change the transport and thermal properties of the base fluid [8]. Therefore, by mixing the nanoparticles in the fluid, thermal conductivity of

The flow and heat transfer characteristics at a stagnation point, for both twodimensional and axisymmetric, have been studied extensively in the literature [9]. These studies have been motivated by the fundamental nature of the boundary layer flows near such points by the exact applicability of similar solutions and by their relevance to the leading edge and nose regions of bodies in high speed flows. These two cases of two-dimensional and axisymmetric flows can be recognized as special cases of more general stagnation point flows. Both two-dimensional and axisymmetric flows were extended to three-dimensional by Howarth [10]. Bhattacharyya and Gupta [11] studied the flow and heat transfer in an incompressible viscous and electrically conducting fluid near a three-dimensional stagnation point of a body permeated by a uniform magnetic field. Analysis of such flows is very important in both theory and in practice. From a theoretical point of view, flows of this type are fundamental in fluid mechanics and forced convective heat transfer. On the other hand, from a practical point of view, these flows have applications in many manufacturing processes in industry such as the boundary layer along material handling conveyers, the aerody-

It is worth mentioning that the nanofluid model proposed by Buongiorno [12] was very recently used by Nield and Kuznetsov [13], Kuznetsov and Nield [14] and

The aim of the present chapter is to study the Magnetohydrodynamic effects on

Consider a flow of an electrically conducting fluid with heat transfer flow at a stagnation point of an incompressible viscous fluid past a body that is placed in a nanofluid of ambient uniform temperature *T*<sup>∞</sup> through a porous medium in the presence of radiation has been considered, where the body surface is kept at a constant temperature *Tw*. The stagnation point is located at the origin *o* of the Cartesian coordinate system *oxyz*. It was shown by Howarth [10] that the inviscid

div*B* ¼ 0*, CurlB* ¼ *μmJ and* div*E* ¼ 0

*ue*ð Þ¼ *x a x, ve*ð Þ¼ *y b y* (1)

heat transfer and thermal radiation in an incompressible viscous fluid near the three-dimensional stagnation point of a body that is placed in a water based nanofluid through a porous medium containing different types of nanoparticles: copper (Cu), alumina (Al2O3) and titania (TiO2). Numerical results are presented for velocity and temperature profiles for different parameters of the problem.

the fluid increases and the heat transfer capability improves.

*Nanofluid Flow in Porous Media*

namic extrusion of plastic sheet and petrochemical industries.

**2. Mathematical analysis**

**142**

irrational flow near *o* has the velocity components

with the constants *a* > 0 and *b* ≥ 0 or *b* ≤ 0.

The MHD body forces *J* � *B* the Maxwell's equations:

Bachok et al. [15] in their papers. Different from the above model.

Further, we look for a solution of Eqs. (2)–(5) of the form

$$\begin{aligned} u &= a \text{x} f'(\eta), \quad v = b \text{y} g'(\eta), \quad w = -\sqrt{a\nu\_f}(f + c\text{g}), \\ \theta(\eta) &= \frac{(T - T\_{\infty})}{(T\_w - T\_{\infty})}, \quad \eta = \sqrt{a/\nu\_f} \text{z} \end{aligned} \tag{10}$$

Using Eqs. (10), (15), and (16), we obtain

ð Þ 1 � *φ*

*kf θ=* ð Þ 0

*==*ð Þ <sup>0</sup> , *<sup>g</sup>==*ð Þ <sup>0</sup> and �*θ<sup>=</sup>*

**Figures 2, 3** and **5a–c**, display the velocity *f*

*Effect of nanoparticle volume fractions on (a) the velocity profile f*

*temperature profile θ η*ð Þ *at* c *= 0.5 and* �*0.5 for copper-water nanofluid.*

<sup>5</sup>*=*<sup>2</sup> *f*

*==*ð Þ <sup>0</sup> *,* Re1*=*<sup>2</sup>

*Heat Transfer and Thermal Radiation at a General Three-Dimensional in a Nanofluid…*

*<sup>w</sup>=aν<sup>f</sup>* is the local Reynolds number.

obtained using the modified fourth order Runge-Kutta method along with

Numerical solutions to the nonlinear ordinary differential Eqs. (11)–(14) were

Nachtsheim-Swigert shooting technique [17]. In order to gain physical insight, the velocity and the temperature profiles have been discussed by assigning numerical values to the parameter, encountered in the problem which the numerical results are tabulated and displayed with the graphical illustrations **Figure 1**. We find the

ters, namely the nanoparticle volume fraction *φ* and the ratio of the gradient of velocities in the i- and *x*-directions *c*, where 0 < *c* < 1 displays the nodal stagnation point of attachment and �1 < *c* < 0 displays the saddle points of attachment. Three types of nanoparticles were considered, namely, copper (Cu), alumina (Al2O3) and titania (TiO2). Following Oztop and Abu-Nada [18], the value of the Prandtl number *Pr* is taken as 6.2 (for water) and the volume fraction of nanoparticles is from 0 to 0.2 0ð Þ ≤*φ*≤0*:*2 in which *φ* corresponds to the regular Newtonian fluid. The numerical results are summarized in **Figures 2**–**6**. It is worth mentioning that we

have used data related to the thermophysical properties of the fluid and nanoparticles as listed in **Table 1** [18] to compute each case of the nanofluid.

profiles under the different parameters of the problem as magnetic field parameter, porous medium parameter and nanoparticle parameter. It is observed that the velocity profile decreases, while the temperature profile increases with the increase in each of magnetic field parameter, porous medium parameter and nanoparticle parameter. The effects of the nanoparticle volume fractions on the velocity and the temperature profiles are shown in **Figure 4a**–**c**, respectively. It is observed that the

*=*

*=*

*<sup>x</sup> Cfy* <sup>¼</sup> <sup>1</sup>

ð Þ 1 � *φ*

<sup>5</sup>*=*<sup>2</sup> *<sup>g</sup>==*ð Þ <sup>0</sup>

ð Þ 0 for some values of the governing parame-

ð Þ*<sup>η</sup> , g<sup>=</sup>*ð Þ*<sup>η</sup>* and the temperature *θ η*ð Þ

ð Þ*<sup>η</sup> , (b) the velocity profile <sup>g</sup><sup>=</sup>*ð Þ*<sup>η</sup> and (c) the*

(17)

*<sup>x</sup> Cfx* <sup>¼</sup> <sup>1</sup>

*DOI: http://dx.doi.org/10.5772/intechopen.81575*

Re�1*=*<sup>2</sup> *<sup>x</sup> Nux* ¼ � *knf*

Re1*=*<sup>2</sup>

where Re*<sup>x</sup>* <sup>¼</sup> *<sup>u</sup>*<sup>4</sup>

missing slopes *f*

**Figure 1.**

**145**

**4. Results and discussion**

.

where *c* = *b*/*a* is the ratio of the gradient of velocities in the *y*- and *x*-directions, and primes denote differentiation with respect to *z*. Substituting Eqs. (8)–(10) in Eqs. (2)–(5), we obtain the following nonlinear ordinary differential equations:

$$\left(\frac{1}{\left(1-\rho\right)^{5/2}\left(1-\rho\left(1-\rho\_{\mathrm{s}}/\rho\_{\mathrm{f}}\right)\right)}\right)\left[f^{\prime\prime\prime}-\mathcal{M}(1-\rho)^{5/2}f^{\prime}-\mathcal{S}f^{\prime}\right]+(\mathcal{f}+\mathsf{c}\boldsymbol{\mathfrak{g}})f^{\prime\prime\prime}-f^{\prime2}+\mathbbm{1}=\mathsf{0},\tag{11}$$

$$\left(\frac{1}{(\mathbf{1}-\boldsymbol{\varrho})^{5/2}\Big(\mathbf{1}-\boldsymbol{\varrho}\Big(\mathbf{1}-\boldsymbol{\rho}\_{i}/\rho\_{f}\Big)\Big)}\right)\Big[\mathbf{g}^{\prime//} - \mathbf{M}(\mathbf{1}-\boldsymbol{\varrho})^{5/2}\mathbf{g}^{\prime} - \mathbf{S}\mathbf{g}^{\prime}\Big] + (\mathbf{f}+\boldsymbol{\varepsilon}\mathbf{g})\mathbf{g}^{\prime//} + c\Big(\mathbf{1}-\mathbf{g}^{\prime2}\Big) = \mathbf{0},\tag{12}$$

$$\left(\frac{\left(\left(k\_{\text{uf}}/k\_f\right)+Q\right)}{P\_r\left(1-\varrho\left(1-\left(\left(\rho\mathbf{C}\_p\right)\_s/\left(\rho\mathbf{C}\_p\right)\_f\right)\right)\right)}\right)\theta^{/\left(}+\left(f+\epsilon g\right)\theta^{\prime}=\mathbf{0},\tag{13}$$

With the appropriate boundary conditions:

$$\begin{aligned} f'(\mathbf{0}) &= \mathbf{0}, \mathbf{g}'(\mathbf{0}) = \mathbf{0}, f(\mathbf{0}) = -c \, \mathbf{g}(\mathbf{0}), \theta(\mathbf{0}) = \mathbf{1}, \\ f'(\infty) &= \mathbf{1}, \mathbf{g}'(\infty) = \mathbf{1}, \theta(\infty) = \mathbf{0}. \end{aligned} \tag{14}$$

where the prime denotes a partial differentiation with respect to *<sup>η</sup>*, *<sup>M</sup>* <sup>¼</sup> *<sup>σ</sup> <sup>B</sup>*<sup>2</sup> *=aρ<sup>f</sup>* is the Magnetic field parameter, *S* ¼ *ν<sup>f</sup> =aκ* is the Porous medium parameter and *<sup>Q</sup>* <sup>¼</sup> <sup>4</sup>*k*1*T*<sup>3</sup> <sup>∞</sup> is the Thermal radiation parameter.

There is no loss of generality in the requirement that j j *a* ≥j j *b* with *a*>0. Clearly *b* ¼ 0 corresponds to the plane stagnation point flow case, while *b* ¼ *a* is the axisymmetric case. The case 0 < *c* < 1 displays the nodal stagnation point of attachment (*b*>0), while �1 < *c* < 0 displays the saddle points of attachment (*b* < 0) [5].

#### **3. Skin-friction coefficient and Nusselt number**

The parameters of engineering interest for the present problem are the local skin-friction coefficients *Cfx* and *Cfy*, along the *x*- and *y*-directions, respectively, and the local Nusselt number *Nux*, which are defined as

$$\mathbf{C}\_{\mathbf{f}\mathbf{x}} = \frac{\mathbf{r}\_{w\mathbf{x}}}{a\mathbf{x}\rho\_f u\_w^2}, \quad \mathbf{C}\_{\mathbf{f}\mathbf{y}} = \frac{\mathbf{r}\_{w\mathbf{y}}}{b\mathbf{y}\rho\_f u\_w^2}, \quad \mathbf{N}u\_\mathbf{x} = \frac{u\_w^2 q\_w}{a k\_f (T\_w - T\_\infty)}\tag{15}$$

where *τwx* and *τwy* are the surface shear stresses along the *x*- and *y*-directions, respectively, and *qw* is the surface heat flux, which are given by

$$\pi\_{w\mathbf{x}} = \mu\_{\eta f} \left( \frac{\partial u}{\partial \mathbf{z}} \right)\_{\mathbf{z}=0}, \quad \pi\_{w\mathbf{y}} = \mu\_{\eta f} \left( \frac{\partial v}{\partial \mathbf{z}} \right)\_{\mathbf{z}=0}, \quad q\_w = -k\_{\eta f} \left( \frac{\partial T}{\partial \mathbf{z}} \right)\_{\mathbf{z}=0} \tag{16}$$

*Heat Transfer and Thermal Radiation at a General Three-Dimensional in a Nanofluid… DOI: http://dx.doi.org/10.5772/intechopen.81575*

Using Eqs. (10), (15), and (16), we obtain

$$\begin{aligned} \text{Re}\_{\mathbf{x}}^{1/2} \mathbf{C}\_{\mathbf{f}\mathbf{x}} &= \frac{\mathbf{1}}{\left(\mathbf{1} - \boldsymbol{\varrho}\right)^{5/2}} f^{\prime \prime}(\mathbf{0}), \quad \text{Re}\_{\mathbf{x}}^{1/2} \mathbf{C}\_{\mathbf{f}\mathbf{y}} = \frac{\mathbf{1}}{\left(\mathbf{1} - \boldsymbol{\varrho}\right)^{5/2}} \mathbf{g}^{\prime \prime}(\mathbf{0}) \\\ \text{Re}\_{\mathbf{x}}^{-1/2} \mathbf{N} u\_{\mathbf{x}} &= -\frac{k\_{\eta \mathbf{f}}}{k\_{\mathbf{f}}} \boldsymbol{\theta}^{\prime}(\mathbf{0}) \end{aligned} \tag{17}$$

where Re*<sup>x</sup>* <sup>¼</sup> *<sup>u</sup>*<sup>4</sup> *<sup>w</sup>=aν<sup>f</sup>* is the local Reynolds number.

## **4. Results and discussion**

.

Further, we look for a solution of Eqs. (2)–(5) of the form

ð Þ*<sup>η</sup> , v* <sup>¼</sup> *byg<sup>=</sup>*

*, η* ¼

*===* � *<sup>M</sup>*ð Þ <sup>1</sup> � *<sup>φ</sup>*

<sup>A</sup> *<sup>g</sup>===* � *<sup>M</sup>*ð Þ <sup>1</sup> � *<sup>φ</sup>*

ð Þ*η , w* ¼ � ffiffiffiffiffiffiffi

*z*

ffiffiffiffiffiffiffiffiffi *a=ν<sup>f</sup>* q

> 5*=*2 *f*

5*=*2

h i

*<sup>=</sup>* ð Þ¼ <sup>0</sup> <sup>0</sup>*, g<sup>=</sup>*ð Þ¼ <sup>0</sup> <sup>0</sup>*, f* ð Þ¼� <sup>0</sup> *c g*ð Þ <sup>0</sup> *, <sup>θ</sup>* ð Þ¼ <sup>0</sup> <sup>1</sup>*,*

where the prime denotes a partial differentiation with respect to *<sup>η</sup>*, *<sup>M</sup>* <sup>¼</sup> *<sup>σ</sup> <sup>B</sup>*<sup>2</sup>

There is no loss of generality in the requirement that j j *a* ≥j j *b* with *a*>0. Clearly

The parameters of engineering interest for the present problem are the local skin-friction coefficients *Cfx* and *Cfy*, along the *x*- and *y*-directions, respectively, and

> *byρ<sup>f</sup> u*<sup>2</sup> *w*

where *τwx* and *τwy* are the surface shear stresses along the *x*- and *y*-directions,

*∂ v ∂ z* � �

*z*¼0

*, Cfy* <sup>¼</sup> *<sup>τ</sup>wy*

respectively, and *qw* is the surface heat flux, which are given by

*, τwy* ¼ *μnf*

is the Magnetic field parameter, *S* ¼ *ν<sup>f</sup> =aκ* is the Porous medium parameter and

*b* ¼ 0 corresponds to the plane stagnation point flow case, while *b* ¼ *a* is the axisymmetric case. The case 0 < *c* < 1 displays the nodal stagnation point of attachment (*b*>0), while �1 < *c* < 0 displays the saddle points of attachment (*b* < 0) [5].

*<sup>g</sup><sup>=</sup>* � *S g<sup>=</sup>*

1

*<sup>=</sup>* ð Þ¼ <sup>∞</sup> <sup>1</sup>*, g<sup>=</sup>*ð Þ¼ <sup>∞</sup> <sup>1</sup>*, <sup>θ</sup>* ð Þ¼ <sup>∞</sup> <sup>0</sup>*:* (14)

*, Nux* <sup>¼</sup> *<sup>u</sup>*<sup>2</sup>

*, qw* ¼ �*knf*

*wqw akf*ð Þ *Tw* � *T*<sup>∞</sup>

> *∂T ∂ z* � �

*z*¼0

*<sup>=</sup>* � *S f <sup>=</sup>* h i

where *c* = *b*/*a* is the ratio of the gradient of velocities in the *y*- and *x*-directions, and primes denote differentiation with respect to *z*. Substituting Eqs. (8)–(10) in Eqs. (2)–(5), we obtain the following nonlinear ordinary differential equations:

*aν<sup>f</sup>* <sup>p</sup> ð Þ *<sup>f</sup>* <sup>þ</sup> *c g ,*

þ ð Þ *f* þ *c g f*

*==* � *<sup>f</sup>*

<sup>þ</sup> ð Þ *<sup>f</sup>* <sup>þ</sup> *c g <sup>g</sup>==* <sup>þ</sup> *<sup>c</sup>* <sup>1</sup> � *<sup>g</sup><sup>=</sup>*<sup>2</sup>

<sup>A</sup>*θ==* <sup>þ</sup> ð Þ *<sup>f</sup>* <sup>þ</sup> *c g <sup>θ</sup><sup>=</sup>* <sup>¼</sup> <sup>0</sup>*,* (13)

(10)

*<sup>=</sup>*<sup>2</sup> <sup>þ</sup> <sup>1</sup> <sup>¼</sup> <sup>0</sup>*,*

(11)

(12)

¼ 0*,*

*=aρ<sup>f</sup>*

(15)

(16)

� �

*<sup>u</sup>* <sup>¼</sup> *axf <sup>=</sup>*

1

*Nanofluid Flow in Porous Media*

1

<sup>5</sup>*=*<sup>2</sup> <sup>1</sup> � *<sup>φ</sup>* <sup>1</sup> � *<sup>ρ</sup>s=ρ<sup>f</sup>* � � � �

> 0 @

> > *f*

*f*

<sup>5</sup>*=*<sup>2</sup> <sup>1</sup> � *<sup>φ</sup>* <sup>1</sup> � *<sup>ρ</sup>s=ρ<sup>f</sup>* � � � �

ð Þ 1 � *φ*

ð Þ 1 � *φ*

*<sup>Q</sup>* <sup>¼</sup> <sup>4</sup>*k*1*T*<sup>3</sup>

0 @

0 @ *θ η*ð Þ¼ ð Þ *<sup>T</sup>* � *<sup>T</sup>*<sup>∞</sup> ð Þ *Tw* � *T*<sup>∞</sup>

> 1 A *f*

1

*knf =kf* � � <sup>þ</sup> *<sup>Q</sup>* � �

<sup>∞</sup> is the Thermal radiation parameter.

**3. Skin-friction coefficient and Nusselt number**

the local Nusselt number *Nux*, which are defined as

*Cfx* <sup>¼</sup> *<sup>τ</sup>wx axρ<sup>f</sup> u*<sup>2</sup> *w*

> *∂u ∂ z* � �

*z*¼0

*τwx* ¼ *μnf*

**144**

� � *s = ρCp* � � *f*

� � � � � �

*Pr* 1 � *φ* 1 � *ρCp*

With the appropriate boundary conditions:

Numerical solutions to the nonlinear ordinary differential Eqs. (11)–(14) were obtained using the modified fourth order Runge-Kutta method along with Nachtsheim-Swigert shooting technique [17]. In order to gain physical insight, the velocity and the temperature profiles have been discussed by assigning numerical values to the parameter, encountered in the problem which the numerical results are tabulated and displayed with the graphical illustrations **Figure 1**. We find the missing slopes *f ==*ð Þ <sup>0</sup> , *<sup>g</sup>==*ð Þ <sup>0</sup> and �*θ<sup>=</sup>* ð Þ 0 for some values of the governing parameters, namely the nanoparticle volume fraction *φ* and the ratio of the gradient of velocities in the i- and *x*-directions *c*, where 0 < *c* < 1 displays the nodal stagnation point of attachment and �1 < *c* < 0 displays the saddle points of attachment. Three types of nanoparticles were considered, namely, copper (Cu), alumina (Al2O3) and titania (TiO2). Following Oztop and Abu-Nada [18], the value of the Prandtl number *Pr* is taken as 6.2 (for water) and the volume fraction of nanoparticles is from 0 to 0.2 0ð Þ ≤*φ*≤0*:*2 in which *φ* corresponds to the regular Newtonian fluid. The numerical results are summarized in **Figures 2**–**6**. It is worth mentioning that we have used data related to the thermophysical properties of the fluid and nanoparticles as listed in **Table 1** [18] to compute each case of the nanofluid.

**Figures 2, 3** and **5a–c**, display the velocity *f =* ð Þ*<sup>η</sup> , g<sup>=</sup>*ð Þ*<sup>η</sup>* and the temperature *θ η*ð Þ profiles under the different parameters of the problem as magnetic field parameter, porous medium parameter and nanoparticle parameter. It is observed that the velocity profile decreases, while the temperature profile increases with the increase in each of magnetic field parameter, porous medium parameter and nanoparticle parameter.

The effects of the nanoparticle volume fractions on the velocity and the temperature profiles are shown in **Figure 4a**–**c**, respectively. It is observed that the

#### **Figure 1.**

*Effect of nanoparticle volume fractions on (a) the velocity profile f =* ð Þ*<sup>η</sup> , (b) the velocity profile <sup>g</sup><sup>=</sup>*ð Þ*<sup>η</sup> and (c) the temperature profile θ η*ð Þ *at* c *= 0.5 and* �*0.5 for copper-water nanofluid.*

**Figure 2.**

*Effect of magnetic field parameter on (a) the velocity profile f =* ð Þ*<sup>η</sup> , (b) the velocity profile <sup>g</sup><sup>=</sup>*ð Þ*<sup>η</sup> and (c) the temperature profile θ η*ð Þ*.*

**Figure 5.**

*profile θ η*ð Þ*.*

**Figure 6.**

**Table 1.**

**147**

*Effect of nanoparticle on (a) the velocity profile f*

*DOI: http://dx.doi.org/10.5772/intechopen.81575*

*Effect of the thermal radiation parameter on the temperature profile θ η*ð Þ*.*

*Thermophysical properties of fluid and nanoparticles [18].*

*=*

*Heat Transfer and Thermal Radiation at a General Three-Dimensional in a Nanofluid…*

ð Þ*<sup>η</sup> , (b) the velocity profile <sup>g</sup><sup>=</sup>*ð Þ*<sup>η</sup> and (c) the temperature*

#### **Figure 3.**

*Effect of porous medium parameter on (a) the velocity profile f =* ð Þ*<sup>η</sup> , (b) the velocity profile <sup>g</sup><sup>=</sup>*ð Þ*<sup>η</sup> and (c) the temperature profile θ η*ð Þ*.*

#### **Figure 4.**

*Effect of nanoparticle volume fractions on (a) the velocity profile f =* ð Þ*<sup>η</sup> , (b) the velocity profile <sup>g</sup><sup>=</sup>*ð Þ*<sup>η</sup> and (c) the temperature profile θ η*ð Þ*.*

velocity and the temperature profiles increase with the increase of the nanoparticle volume fractions. The effects of thermal radiation parameter on the temperature profile are shown in **Figure 6**; also, we found that the temperature profiles increase with the increase of thermal radiation parameters. i.e. as expected, since the effect of thermal radiation is to decrease the rate of energy transport to the fluid.

*Heat Transfer and Thermal Radiation at a General Three-Dimensional in a Nanofluid… DOI: http://dx.doi.org/10.5772/intechopen.81575*

**Figure 5.**

*Effect of nanoparticle on (a) the velocity profile f =* ð Þ*<sup>η</sup> , (b) the velocity profile <sup>g</sup><sup>=</sup>*ð Þ*<sup>η</sup> and (c) the temperature profile θ η*ð Þ*.*


**Table 1.**

*Thermophysical properties of fluid and nanoparticles [18].*

velocity and the temperature profiles increase with the increase of the nanoparticle volume fractions. The effects of thermal radiation parameter on the temperature profile are shown in **Figure 6**; also, we found that the temperature profiles increase with the increase of thermal radiation parameters. i.e. as expected, since the effect

*=*

of thermal radiation is to decrease the rate of energy transport to the fluid.

**Figure 2.**

**Figure 3.**

**Figure 4.**

**146**

*temperature profile θ η*ð Þ*.*

*temperature profile θ η*ð Þ*.*

*temperature profile θ η*ð Þ*.*

*Nanofluid Flow in Porous Media*

*Effect of magnetic field parameter on (a) the velocity profile f*

*Effect of porous medium parameter on (a) the velocity profile f*

*Effect of nanoparticle volume fractions on (a) the velocity profile f*

*=*

*=*

ð Þ*<sup>η</sup> , (b) the velocity profile <sup>g</sup><sup>=</sup>*ð Þ*<sup>η</sup> and (c) the*

ð Þ*<sup>η</sup> , (b) the velocity profile <sup>g</sup><sup>=</sup>*ð Þ*<sup>η</sup> and (c) the*

ð Þ*<sup>η</sup> , (b) the velocity profile <sup>g</sup><sup>=</sup>*ð Þ*<sup>η</sup> and (c) the*


2. The velocity and the temperature profiles increase with the increase of the

*Heat Transfer and Thermal Radiation at a General Three-Dimensional in a Nanofluid…*

3. For the effects of thermal radiation parameter on the temperature profile, we found that the temperature profiles increase with the increase of thermal

4.For an increase in the magnetic field parameter *M*, we observe that the skinfriction coefficient along *x*- and *y*-directions and the local Nusselt number decrease when *c* = 0.5, but at *c* = 0.5, the skin-friction coefficient along *x*direction and the local Nusselt number decrease, while the skin-friction

5. About the effects of the porous medium parameter *S*, we observe that the skinfriction coefficient and the local Nusselt number decrease when *c* = 0.5 and

6.The local Nusselt number decreases when *c* = 0.5 and *c* = 0.5, with the

*c* the ratio of the gradient of velocities in the *y*- and *x*-directions

*qr* the energy flux which is relative to a frame moving with the

nanoparticle volume fractions.

*DOI: http://dx.doi.org/10.5772/intechopen.81575*

coefficient along *y*-direction increases.

increase of the thermal radiation parameter *Q*.

*Cfx*, *Cfy* skin friction in the *y*- and *x*-directions

*k*<sup>l</sup> the nanofluid thermal conductivity *k*\* mean absorption coefficient

nanofluid velocity

*αnf* the thermal diffusivity of the nanofluid

*T*<sup>∞</sup> ambient uniform temperature *S* porous medium parameter Re*<sup>x</sup>* the local Reynolds number *u* velocity in the *x*-direction *v* velocity in the *y*-direction *w* velocity in the *z*-direction *x* horizontal distance *y* vertical distance *z* normal distance

*knf* the thermal conductivity of the nanofluid *kf* the thermal conductivities of the fluid

*ks* the thermal conductivities of the solid fractions

*d* the induced magnetic field

*M* magnetic field parameter *Q* thermal radiation parameter

*qw* the surface heat flux *T* temperature distribution *T*<sup>w</sup> constant temperature

*η* similarity variable

*Nu*<sup>x</sup> Nusselt number

radiation parameters.

*c* = 0.5.

**Nomenclature**

Greek symbols

**149**

*a* constant *B* magnetic field *b* constant

#### **Table 2.**

*Numerical of the values skin-friction coefficient (*Re<sup>1</sup>*=*<sup>2</sup> *<sup>x</sup> Cfx) and (*Re1*<sup>=</sup>*<sup>2</sup> *<sup>x</sup> Cfy) along x- and y-directions and the local Nusselt number (*Re�1*=*<sup>2</sup> *<sup>x</sup> Nux) with <sup>M</sup>*, *<sup>S</sup>*, *<sup>φ</sup> andQ , when Pr* <sup>¼</sup> <sup>6</sup>*:*2*.*

From **Table 2**, the numerical values of the skin-friction and the local Nusselt number are given in **Table 2**. For an increase in the magnetic field parameter *M*, we observe that the skin-friction coefficient along *x*- and *y*-directions and the local Nusselt number decrease when *c* = 0.5, but at *c* = �0.5, the skin-friction coefficient along *x*-direction and the local Nusselt number decrease, while the skin-friction coefficient along *y*-direction increases. While, with an increase in the porous medium parameter *S*, we observe that the skin-friction coefficient and the local Nusselt number decrease when *c* = �0.5 and *c* = �0.5.

With an increase in nanoparticle volume fractions *φ*, we observe that the skinfriction coefficient increases and the local Nusselt number increases when *c* = 0.5, but at *c* = �0.5, the skin-friction coefficient along *y*-direction and the local Nusselt number decrease, while the skin-friction coefficient along *x*-direction increases. While, with an increase in the thermal radiation parameter *Q*, we observe that the local Nusselt number decreases when *c* = �0.5 and *c* = �0.5.

#### **5. Conclusions**

In the present chapter, the magnetohydrodynamic effects on heat transfer and thermal radiation at a stagnation point flowing in a nanofluid containing different types of nanoparticles namely, copper (Cu), alumina (Al2O3) and titania (TiO2) through a porous medium have been investigated numerically. By using appropriate transformation for velocity and temperature into a set of non-linear coupled ordinary differential equations which are solved numerically, the governing equations were:

1. With the increasing values of the magnetic field parameter, porous medium parameter and nanoparticle parameter, it is observed that the velocity profile decreases, while the temperature profile increases.

*Heat Transfer and Thermal Radiation at a General Three-Dimensional in a Nanofluid… DOI: http://dx.doi.org/10.5772/intechopen.81575*

