Open access peer-reviewed chapter

# Weakly Nonlinear Stability Analysis of a Nanofluid in a Horizontal Porous Layer Using a Multidomain Spectral Collocation Method

Written By

Osman A.I. Noreldin, Precious Sibanda and Sabyasachi Mondal

Submitted: May 29th, 2017 Reviewed: September 19th, 2017 Published: December 20th, 2017

DOI: 10.5772/intechopen.71066

From the Edited Volume

## Complexity in Biological and Physical Systems

Edited by Ricardo López-Ruiz

Chapter metrics overview

View Full Metrics

## Abstract

In this chapter, we present a weakly nonlinear stability analysis of the flow of a nanofluid in a porous medium with stress-free boundary conditions. Some previous studies have investigated cross-diffusion in a nanofluid layer although in most cases these studies mostly deal with linear stability analysis. It is important to study the nonlinear stability in flows subject to cross-diffusion due to the wide range of applications where such flows arise such as in hydrothermal growth, compact heat exchanges, the solidification of binary mixtures, geophysical systems, solar pond, etc. Here we consider flow between parallel plates with an applied magnetic field and zero nanoparticle flux at the boundaries. A truncated Fourier series is introduced reducing the flow equations to a Lorenz-type system of nonlinear evolution equations. The multidomain spectral method is used to solve the equations that describe the growth of the convection amplitudes. The solutions are obtained as sets of trajectories in the phase space. Some interesting spiral trajectories and their sensitivity to the Rayleigh number are given.

### Keywords

• nonlinear instability
• nanofluid flow
• porous medium
• multidomain spectral collocation method

## 1. Introduction

The enhancement of thermal conductivity of a fluid is a matter of supreme interest to engineers due to the important applications of fluids in heat transfer processes. Natural and forced convection plays an important role in heat transfer processes due to continuous molecular movements in fluid. Recent studies show that the suspension of solid nanoparticles in a fluid can substantially improve the fluid’s thermophysical properties, including thermal conductivity.

The term nanofluid describes a liquid containing a suspension of nanometer sized 1–100 nm solid particles [1]. Examples of commonly used nanoparticles include metallic particles such as Al, Cu and Ag, and oxides such as Al2O3 and CuO. The base fluid is often a common liquid such as water, ethylene, glycol, or oil. The enhancement of thermochemical properties of a fluid due to the addition of nanoparticles has been observed in experimental studies such as in [2, 3]. Researchers have investigated the influence of seven slip mechanisms, namely, inertia, Brownian diffusion, thermophoresis, diffusiophoresis, magnus effect, fluid drainage, and gravity in nanofluids. It has been shown that, in the absence of turbulence, the most significant among these mechanisms are the Brownian diffusion and thermophoresis.

The classical Rayleigh-Benard convection problem in a heated horizontal layer has been extensively studied in the literature. Among recent studies on nanofluids, Tzou [4] studied the thermal instability and natural convection in nanofluid flow using an eigenfunction expansion method. Narayana et al. [5, 6] studied convection and the stability of a Maxwell fluid in a porous medium. Yadav et al. [7] investigated thermal instability of a rotating nanofluid layer. The studies by Kuznetsov and Nield [8, 9, 10, 11] focused on thermal instability in a porous layer saturated with a nanofluid. They investigated the onset of instability in a horizontal porous layer using a model for the nanofluid that incorporated particle Brownian motion and thermophoresis. Related studies with various assumptions on the geometry and flow structure have been made by [12, 13, 14, 15]. In the last few decades, researchers have also investigated thermal instability in a horizontal nanofluid layer subject to an applied magnetic field [16, 17]. The effects of a magnetic field on convection and the onset of instability have important applications in problems such as in cooling systems, pumps, magnetohydrodynamics and generators. The experimental study by Heris et al. [18] showed that thermal efficiency could be achieved by subjecting the flow to a magnetic field. The studies by Ghasemi et al. [19] and Hamad et al. [20] focused on the flow behavior and heat transfer in an electrically conducting nanofluid under the influence of a magnetic field and subject to Brownian diffusion and thermophoresis. They used a water-based nanofluid containing different types of nanoparticles such as copper, alumina and silver in their numerical simulations. Related studies of interest include [21, 22, 23, 24]. Rana et al. [25] studied thermal convection in a Walters (Model B) fluid in a porous medium. They showed that a magnetic field may introduce oscillatory instability modes and acts to stabilize the system.

In this chapter, we give a weakly nonlinear stability analysis of a nanofluid layer with an applied magnetic field, stress free boundary conditions and under the assumption of zero nanoparticle flux at the boundary. The studies by Kuznetsov and Nield [9] and Nield and Kuznetsov [10, 11] investigated cross-diffusion in a nanofluid layer. However, these studies mostly presented a linear stability analysis. It is important to study the nonlinear regime for a nanofluid flow subject to cross-diffusion due to the wide range of applications where such flows may arise. Typical examples may be found in hydrothermal growth, compact heat exchanges, solidification of binary mixtures, geophysical systems, and so on. Hence, with this in mind, we studied the finite amplitude convection in a nanofluid flows subject to cross-diffusion. By introducing a truncated Fourier series, a Lorenz-type system of seven nonlinear differential equations is obtained. The recent multidomain spectral method is used to solve the nonlinear equations. This method is accurate and very easy to implement compared to older methods such as finite difference methods. An analysis of heat and mass transfer for different parameters such as the Prandtl number, the Dufour and thermophoresis is presented.

## 2. Mathematical formulation

Consider viscous incompressible MHD nanofluid flow in an infinitely extended horizontal porous layer, confined between two boundaries at z = 0 and z = h . The layer is heated from below and cooled from above, see Figure 1. A Cartesian frame of reference is chosen in which the z-axis is vertically upward. The boundaries are perfectly conducting. The temperature at the lower and upper walls is T c and T h , respectively with T h > T c . The Oberbeck-Boussinesq approximation and the Darcy law are assumed to be applicable. The continuity equation, momentum equation, energy equation, concentration equation and volumetric fraction nanoparticle equation, which describe the above configuration in dimensionless form, are given as

V = 0 , E1
Da Pr V t = P + Da 2 V V + QV e ̂ z Rm e ̂ z + RaT e ̂ z + RsC e ̂ z Rnϕ e ̂ z , E2
T t + V T = 2 T + N B Les ϕ T + N A N B Les T T + Du 2 C , E3
C t + V C = 1 Le 2 C + Sr 2 T , E4
1 σ ϕ t + 1 ε V ϕ = 1 Les 2 ϕ + N A Les 2 T , E5

subject to the boundary conditions

V = 0 , T = 1 , C = 1 ϕ z + N A T z = 0 at z = 0 , E6
V = 0 , T = 0 , C = 0 ϕ z + N A T z = 0 at z = 1 , E7

where V is the fluid velocity, T is the temperature, C is the solute concentration and ϕ is the volumetric fraction of nanoparticles. The dimensionless parameters are the Darcy number (modified by the viscosity ratio) Da , Prandtl number Pr , Hartmann-Darcy number Q , thermal Rayleigh-Darcy number Ra , nanoparticle Rayleigh number Rn and the basic density Rayleigh number Rm . The parameter N A is a modified diffusivity ratio, Le is the Lewis number, Rs is solutal Rayleigh number, N B is a modified nanoparticle density increment and Du is a modified Dufour parameter. The parameter Les is the thermo-nanofluid Lewis number, ν is the kinematic viscosity and Sr is a modified Soret parameter. These parameters have the form

Da = μ ˜ K μ h 2 , Pr = μ ρ f α m , Q = δ B 0 2 K μ , Ra = ρ f βKhg T h T c μα m , Les = α m D B , E8
Rn = ρ p ρ f ϕ 1 ϕ 0 gKh μα m , Rm = ρ p ϕ 0 + 1 ϕ 0 ρ f gKh μα m , Le = α m D S , E9
N B = ε ρc p ρc f ϕ 1 ϕ 0 , N A = D T T h T c D B T c ϕ 1 ϕ 0 , Rs = ρ f βKhg C h C c μα m , E10
Du = σ D TC C h C c α m T h T c , Sr = σ D CT T h T c α m C h C c , E11

where ρ f , ρ p , μ ˜ , β 1 , β 2 , κ m , δ , ε and K are the fluid density, nanoparticle density, effective viscosity of porous medium, thermal volumetric expansion coefficient of the fluid, solutal volumetric expansion coefficient, the thermal conductivity of porous medium, the electrical conductivity, the porosity, and permeability of porous medium, respectively. The gravitational acceleration is denoted by g and D B is the Brownian diffusion coefficient, D T is the thermophoresis diffusion coefficient, D S is the solutal diffusion coefficient, D TC is the Dufour parameter and D CT is the Soret parameter. The heat capacity of the fluid is ρc f , ρc p is the effective heat capacity of the nanoparticle, ρc m is the effective heat capacity of the porous medium and B 0 is the uniform magnetic field strength.

The basic state is the time independent solution of Eqs. (1)(5). Solving these equations with boundary conditions, we obtain

T b = 1 z , C b = 1 z , ϕ b = ϕ 0 + N A z . E12

## 3. Weakly nonlinear stability analysis

In this section, we restrict the analysis to the case of two-dimensional disturbances. We define the stream function Ψ by the equations

u = Ψ z , w = Ψ x .

Eqs. (1)(5) may now be simplified by introducing the truncated Fourier series

Ψ = A 11 sin αx sin πz , T = B 11 cos αx sin πz + B 02 sin 2 πz , E13
C = C 11 cos αx sin πz + C 02 sin 2 πz , ϕ = N A D 11 cos αx sin πz + D 02 sin 2 πz , E14

where A 11 , B 11 , B 02 , C 11 , C 02 , D 11 and D 02 are amplitudes that depend on time. This leads to the Lorenz-type system of nonlinear ordinary differential equations

Y ̇ 1 = Pr Da BY 1 N Y 2 + Y 4 + N A Rn R Y 6 E15
Y ̇ 2 = RY 1 Y 2 DuY 4 Y 1 Y 3 E16
Y ̇ 3 = 1 2 Y 1 Y 2 G Y 3 + DuY 5 E17
Y ̇ 4 = RY 1 1 Les Y 4 SrY 2 Y 1 Y 5 E18
Y ̇ 5 = 1 2 Y 1 Y 2 G 1 Les Y 5 + SrY 3 E19
Y ̇ 6 = σR ε Y 1 σ Le Y 6 Y 2 N A σ ε Y 1 Y 7 E20
Y ̇ 7 = N A σ 2 ε Y 1 Y 7 G Le Y 7 Y 3 E21

subject to Y n 0 = Y n 0 for n = 1 , 2 , , 7 . The following variables have been introduced in the equations above:

Y 1 = απ γ A 11 , Y 2 = π RB 11 , Y 3 = π RB 02 , Y 4 = π RC 11 , Y 5 = π RC 20 , Y 6 = π RD 11 , Y 7 = π RD 20 , t = γt , R = α 2 γ 3 Ra , G = 4 π 2 γ and N = Rs Ra ,
B = Da γ 2 + γ α 2 Q γ 2 .

Eqs. (15)(21) give an approximate description of the full dimensional nonlinear system. An analytical solution of the system of nonlinear ordinary differential Eqs. (15)(21) is not possible for the general time variable t . However, it is possible to discuss the stability of the nonlinear system of equations. The system of equations is uniformly bounded in time and dissipative in the phase space. We can easily show that

i = 1 7 Y ̇ i Y i = DaB Pr + 1 + G + Les 1 + GLes 1 + σ Le 1 + GLe 1 . E22

This is always true if B 0 . As has been shown in previous studies, the trajectories may be attracted to a fixed point, limit cycle or other attractor. For a set of initial points in the phase space occupying a region V 0 at time t = 0 , after a time t > 0 , the end point of the corresponding trajectories fills a volume

V t = V 0 exp DaB Pr + 1 + G + Les 1 + GLes 1 + σ Le 1 + GLe 1 t . E23

Eq. (23) shows that the volume decays exponentially with time. Further, it can be noted that the system of Eqs. (15)(21) are invariant under the transformation

S Y 1 Y 2 Y 3 Y 4 Y 5 Y 6 Y 7 S Y 1 Y 2 Y 3 Y 4 Y 5 Y 6 Y 7 . E24

We obtain the possible stationary points of the nonlinear system of equations by setting Y ̇ i = 0 for i = 1 , 2 , , 7 . One of these stationary points is Y i = 0 and by linearizing about this point, we obtain the Jacobian matrix

A = PrB Da Pr Da 0 PrN Da 0 PrN A Rn DaR 0 R 1 0 Du 0 0 0 0 0 G 0 GDu 0 0 R Sr 0 Les 1 0 0 0 0 0 GSr 0 GLes 1 0 0 σR ε σ Le 1 0 0 0 σ Le 1 0 0 0 GLe 1 0 0 0 GLe 1 . E25

The eigenvalues of the above matrix depend on the various parameters. For the specific parameters R = 10 3 , Da = 20 , Pr = 10 , N = 25 , Du = 0.2 , Sr = 3 , Les = 10 , Le = 5 , σ = 0.05 , G = 3   and ε = 0.04 , the characteristic polynomial is

P λ = λ 7 + 21.41 λ 6 + 651.5 λ 5 + 5391.7 λ 4 + 12772.232 λ 3 370.962 λ 2 2996.712 λ + 545.8

with eigenvalues

λ 1 = 0.2955056985 ,   λ 2 = 0.2402382976 ,   λ 3 = 0.6139990637 ,   λ 4 = 4.886000936 , λ 5 = 5 ,   λ 6 = 5.7228719981 21.9033954659 i ,   λ 7 = 5.7228719981 + 21.9033954659 i .

This stationary point is a saddle point. Nonetheless, because the eigenvalues depend on various parameters, we cannot make general conclusions as to the stability of the system. We note, however, that if we denote the trace of the matrix A by T and the determinant d , then

T = DaB Pr + 1 + G + Les 1 + GLes 1 + σ Le 1 + GLe 1 , E26

and

d = σ Pr G 3 DaLeLes DuNRSr Le DuNRSr ε + BDuSr Le + DuRSr Le + NR LeLes + NR Lesε B LeLes R LeLes . E27

The trace is always negative, but the sign of determinant depends on the parameter values. If d < 0 then

1 N εDuSrRLes + εB NRLe DuSrLes + ε + Le NLe < B + R ε , E28

## 4. Method of solution

To study the influence of various physical parameters on the average Nusselt and Sherwood numbers, we solved the nonlinear system of Eqs. (15)(21) numerically using the multidomain spectral collocation method. This is a novel technique for solving nonlinear initial value problems and parabolic equations with large time domains. It has been suggested in the literature that the method gives better accuracy compared to other methods such as finite difference and Runge-Kutta methods [26]. To apply the multidomain spectral collocation to the nonlinear system of equations, we first divide the interval 0 T into subintervals Ω i = t i 1 t i for i = 1 , 2 , , p . The transformation

t = t i t i 1 2 τ + t i + t i 1 2 E29

is used to transform each subinterval Ω i into the interval 1 1 . The system of Eqs. (15)(21) can be written in the form

dY 1 i dt = Pr Da BY 1 i N Y 2 i + Y 4 + N A Rn Ra Y 6 i , E30
dY 2 i dt = RY 1 i Y 2 i DuY 4 i Y 1 i Y 3 i , E31
dY 3 i dt = 1 2 Y 1 i Y 2 i G Y 3 i + DuY 5 i , E32
dY 4 i dt = RY 1 i 1 Les Y 4 i SrY 2 i Y 1 i Y 5 i , E33
dY 5 i dt = 1 2 Y 1 i Y 2 i G 1 Les Y 5 i + SrY 3 i , E34
dY 6 i dt = σR ε Y 1 i σ Le Y 6 i Y 2 i N A σ ε Y 1 i Y 7 i , E35
dY 7 i dt = N A σ 2 ε Y 1 i Y 7 i G Le Y 7 i Y 3 i , E36

subject to

Y n i t i 1 = Y n i 1 t i 1 for n = 1 , 2 , , 7 . E37

The first step in using the multidomain spectral collocation method (MDSCM) concerns the quasilinearization of Eqs. (30)(36) leading to a system of equations in the form

n = 1 7 a j n r i Y n , r + 1 i dY j , r + 1 i dt = R jr i , E38

subject to

Y n , r + 1 i t i 1 = Y n , r + 1 i 1 t i 1 for n = 1 , 2 , , 7 . E39

where a j n r i and R jr i for j = 1 , 2 , , 7 are given in the Appendix. Having linearized the equations, the second step is to integrate Eqs. (30)(36). To this end, we use the Gauss-Lobatto nodes

τ j i = cos πj Nc , for j = 0 , 1 , , Nc . E40

We approximate the derivatives of the unknown functions Y n , r + 1 i t at the collocation points by

dY n , r + 1 i dt τ j i = k = 0 Nc D jk Y n , r + 1 i τ j j = DU n , r + 1 i j , E41

where D = 2 D / t i t i 1 , D is the Chebyshev differentiation matrix and

U n , r + 1 i = Y n , r + 1 i τ 0 i Y n , r + 1 i τ Nc i T ,

is a vector of the unknown functions at the collocation points. Substituting Eq. (41) into Eqs. (38) and reducing the result into matrix form, we obtain

AU n , r + 1 i = R n i , U n , r + 1 i τ N c i 1 = U n i τ N c i 1 , n = 1 , 2 , , 7 . E42

where the matrices A = A ij and R n i are given in the Appendix.

## 5. Heat and mass transfer

The study of heat and mass transfer in a horizontal nanofluid layer heated from below and cooled from above has important engineering applications. We define the rate of heat transfer by the average Nusselt number Nu t where

Nu t = 1 + α 2 π 0 2 π α T z dx α 2 π 0 2 π α T b z dx z = 0 + Du 1 + α 2 π 0 2 π α C z dx α 2 π 0 2 π α C b z dx z = 0 . E43

Substituting Eqs. (12) and (13) into Eq. (43), we obtain

Nu t = 1 + 2 R Y 3 + Du 1 + 2 R Y 5 . E44

Similarly, the rate of mass transfer stated in terms of the average Sherwood number is

Sh t = 1 + 2 R Y 5 + Sr 1 + 2 R Y 3 E45

## 6. Results and discussion

We have studied the weakly nonlinear instability of nanofluid flow in a horizontal layer with stress free boundary conditions. For numerical simulations, the parameter values were chosen from the literature on nanofluid flow such as [4, 7]. In the literature, the critical Rayleigh number is found when the Darcy number is very large. In this study, we investigated the critical Rayleigh number for low Darcy numbers.

The method of solution described in Section 4 was used to solve Eqs. (15)(21). All computations are carried out up to a value of maximum time t max = 1 , and solutions are obtained using initial conditions selected in the neighborhood of stationary points. Periodic solution sets were obtained for the system of nonlinear equations. We determined the rate of heat and mass transfer as functions of time for different parameter values. The results are shown in Figures 24. Figure 2 shows the effect of the Dufour and Soret parameters on the Nusselt and Sherwood numbers with time t . Figure 2(a) shows how the heat transfer coefficient changes with both the Dufour parameter and time. The heat transfer coefficient increases with the Dufour parameter but eventually settles to a steady value with time. In Figure 2(b), the Soret parameter is similarly shown to enhance the mass transfer coefficient. We investigated the effect of the Prandtl and Lewis numbers (see Figures 3 and 4). An increase in the Lewis number enhances both heat and mass transfer in a nanofluid layer heated from below. However, Figure 3 shows that increasing the Prandtl number reduces the amplitude of oscillatory heat and mass transfer. The Prandtl number can lead to both positive and negative contributions to the Nusselt and Sherwood numbers. It is interesting to note that our investigation shows that the magnetic field parameter has very little effect on the heat and mass transfer for this type of flow.

Figures 511 show the effect of the Rayleigh number on the trajectories projected onto the Y i Y j phase planes. The solution sets provide a visual representation of the system’s behavior with every phase point on the phase space representing the physical state of the system. The convective solution sets for different values of R have been presented with the trajectories projected onto the Y i Y j phase planes. These trajectories spiral toward the fixed point for Rayleigh numbers from 10 2 to 10 4 . The solution sets give spiral phase portraits as R increases and for the high Rayleigh numbers, the trajectories spiral many times before they reach a fixed point.

Figures 58 show the phase portraits projected onto the Y i Y j - plane correspond to a simple spiral for R = 100 . As R is increased to 104, the complexity of the trajectories increases leading to certain chaotic forms. Figures 811 show the trajectories in the three-dimensional phase space. Here, we observe similar solution sets as in the two-dimensional phase portraits.

Figures 12 and 13 show the streamline, isotherm and isoconcentration contours in the nanofluid flow for different values of the Darcy number and buoyancy ratio. Figure 12 displays the streamlines for various values of the buoyancy ratio term. Two different eddies are observed. The clockwise and anticlockwise flows are shown via negative and positive stream function values, respectively. The anticlockwise rotating flow occupies the largest area of the nanofluid layer.

For low buoyancy ratio parameters, the flow structure is significantly influenced by the buoyancy within the whole enclosure. Increasing the buoyancy ratio causes the boundary layer thickness to become thinner. Also, a high buoyancy ratio changes the flow structure, and this impacts significantly on the concentration field, which builds up a vertical stratification in the enclosure. It is interesting to note that for N = 25 , the effect of the solutal buoyancy force is in the opposite direction of the thermal buoyancy force. The isothermal and isoconcentration profiles are situated toward the left wall, while for N = 1 , the thermal and solutal buoyancy forces are equal. For N = 25 , the effect of solutal buoyancy force is in the same direction as the thermal buoyancy force. In such cases, the isothermal and isoconcentration contours are mostly toward the right wall.

We observe that when N = 25 , the stream function values in the central eddies increase because the thickness of the boundary layer increases with the buoyancy ratio. The streamlines and the flow behavior are affected by the change in the buoyancy ratio, but the flow pattern remains unaltered. As N decreases from 1 to −25, the streamlines become very dense to the left side of nanofluid layer while when N increases from 1 to 25, the streamlines are less so. The buoyancy forces that drive the nanofluid motion are mainly due to the temperature gradient.

Three different types of eddies are observed for the isoconcentration contours when N = 25 . Of these, two have a clockwise rotation and one is anticlockwise. It is seen that the small eddy at the right bottom edge is diminished as N decreases from 1 to 25 . Here, the concentration boundary layer decreases due to increasing N values, hence the buoyancy ratio has a significant influence on the concentration gradient. As the buoyancy ratio N increases from 1 to 25 the isoconcentrations become very dense at the bottom of nanofluid layer.

The effect of the Darcy number on the nanofluid flow in the porous medium is shown in detail in Figure 13. The streamline patterns are similar to those depicted in Figure 12. However, as Da increases from 0.05 to 0.07, the rotation of the streamlines changes. Similarly, the isotherm patterns change with increasing Darcy numbers. The value of the center eddies increases with increasing Da. Increasing Da has the effect of increasing the effective fluid viscosity and reducing the thermal and solutal boundary layers.

## 7. Conclusion

We have investigated the onset of thermal instability in a horizontal porous layer of infinite extent in a cross-diffusive nanofluid flow. The focus of the study has been on stress free boundary conditions with zero nanoparticle flux at the wall. A multidomain spectral collocation method was used to solve the system of nonlinear evolution equations. As the Rayleigh number increases to 104, the trajectories spiral many times before reaching a fixed point. The nanofluid convection regime is complex for Rayleigh numbers higher than R = 10 4 , and the flow pattern presents difficulties in interpreting correctly.

Additionally, a change in system parameters, such as an increase in the flow Lewis number, improves the rate of heat and mass transfer in the nanofluid saturated porous media. The Dufour parameter has the effect of increasing heat transfer, while increasing the Soret parameter increases the rate of mass transfer.

## Acknowledgments

The authors are grateful to the University of KwaZulu-Natal and the Claude Leon Foundation, South Africa for financial support.

The terms a j n r i and R jr i for j = 1 , 2 , , 7 in Eq. (38) are given by

a 1 1 r i = γ 3 B , a 1 2 r i = γ 3 , a 1 4 r i = γ 3 N , a 1 6 r i = γ 3 N A Rn Ra , a 1 3 r i = a 1 5 r i = a 1 7 r i = 0 , a 2 1 r i = R Y 3 , r i , a 2 2 r i = 1 , a 2 3 r i = Y 1 , r i , a 2 4 r i = Du , a 2 5 r i = a 2 6 r i = a 2 7 r i = 0 a 3 1 r i = απ 2 u 2 , r i , a 3 2 r i = 1 2 Y 1 , r i , a 3 3 r i = G , a 3 5 r i = GDua 3 4 r i = a 3 6 r i = a 3 7 r i = 0 , a 4 1 r i = R Y 5 , r i , a 4 2 r i = Sr , a 5 3 r i = GSr a 4 4 r i = 1 Les , a 4 5 r i = Y 1 , r i , a 4 3 r i = a 4 6 r i = a 4 7 r i = 0 , a 5 1 r i = 1 2 Y 4 , r i , a 5 4 r i = 1 2 Y 1 , r i , a 5 5 r i = G Les , a 5 2 r i = a 5 6 r i = a 5 7 r i = 0 , a 6 1 r i = γ 1 γ 2 Y 7 , r i , a 6 2 r i = σ Le , a 6 6 r i = σ Le , a 6 7 r i = γ 2 Y 1 , r i , a 6 3 r i = a 6 4 r i = a 6 5 r i = 0 , a 7 1 r i = γ 2 2 Y 6 , r i , a 7 3 r i = G Le , a 7 6 r i = γ 2 2 Y 1 , r i , a 7 7 r i = G Le , a 7 2 r i = a 7 4 r i = a 7 5 r i = 0 , R 1 r i = 0 , R 2 r i = Y 1 , r i Y 3 , r i R 3 r i = 1 2 Y 1 , r i Y 2 , r i , R 4 r i = Y 1 , r i Y 5 , r i , R 5 r i = 1 2 Y 1 , r i Y 4 , r i , R 6 r i = σ N A ε Y 1 , r i Y 7 , r i , R 7 r i = σ N A 2 ε Y 1 , r i Y 6 , r i ,

where γ 1 = σR ε , γ 2 = N A σ ε and γ 3 = Pr Da .

The matrices A ij in Eq. (42) are given by

A nn = diag a n , n ) r i D , A 12 = diag a 1 2 r i , A 16 = diag a 1 6 r i , A 13 = A 14 = A 15 = A 17 = O , A 21 = diag a 2 1 r i , A 23 = diag a 2 3 r i , A 24 = diag a 2 4 r i , A 25 = A 26 = A 27 = O , A 31 = diag a 3 1 r i , A 32 = diag a 3 2 r i , A 35 = diag a 3 5 r i , A 34 = A 36 = A 37 = O , A 41 = diag a 4 1 r i , A 42 = diag a 4 2 r i , A 45 = diag a 4 5 r i , A 43 = A 46 = A 47 = O , A 51 = diag a 5 1 r i , A 53 = diag a 5 3 r i , A 54 = diag a 5 4 r i , A 52 = A 56 = A 57 = O , A 61 = diag a 6 1 r i , A 62 = diag a 6 2 r i , A 67 = diag a 6 7 r i , A 63 = A 64 = A 65 = O , A 71 = diag a 7 1 r i , A 73 = diag a 7 3 r i , A 76 = diag a 7 6 r i , A 72 = A 74 = A 75 = O ,

where O is an N + 1 × N + 1 matrix of zeros and diag is an N + 1 × N + 1 diagonal matrix.

## References

1. 1. Choi SU, Eastman JA. Enhancing Thermal Conductivity of Fluids with Nanoparticles, Technical Report. IL, United States: Argonne National Lab; 1995
2. 2. Buongiorno J. Convective transport in nanofluids. Journal of Heat Transfer. 2006;128(3):240-250
3. 3. Masuda H, Ebata A, Teramae K. Alteration of thermal conductivity and viscosity of liquid by dispersing ultra-fine particles. Netsu Bussei. 1993;4:227-233
4. 4. Tzou D. Thermal instability of nanofluids in natural convection. International Journal of Heat and Mass Transfer. 2008;51(11):2967-2979
5. 5. Narayana M, Sibanda P, Motsa S, Lakshmi-Narayana P. Linear and nonlinear stability analysis of binary maxwell fluid convection in a porous medium. Heat and Mass Transfer. 2012;48(5):863-874
6. 6. Narayana M, Sibanda P, Siddheshwar PG, Jayalatha G. Linear and nonlinear stability analysis of binary viscoelastic fluid convection. Applied Mathematical Modelling. 2013;37(16):8162-8178
7. 7. Yadav D, Bhargava R, Agrawal G. Numerical solution of a thermal instability problem in a rotating nanofluid layer. International Journal of Heat and Mass Transfer. 2013;63:313-322
8. 8. Kuznetsov A, Nield D. Effect of local thermal non-equilibrium on the onset of convection in a porous medium layer saturated by a nanofluid. Transport in Porous Media. 2010;83(2):425-436
9. 9. Kuznetsov A, Nield D. Thermal instability in a porous medium layer saturated by a nanofluid: Brinkman model. Transport in Porous Media. 2010;81(3):409-422
10. 10. Nield D, Kuznetsov A. The onset of double-diffusive convection in a nanofluid layer. International Journal of Heat and Fluid Flow. 2011;32(4):771-776
11. 11. Nield D, Kuznetsov A. Thermal instability in a porous medium layer saturated by a nanofluid: A revised model. International Journal of Heat and Mass Transfer. 2014;68:211-214
12. 12. Bhadauria B. Double diffusive convection in a porous medium with modulated temperature on the boundaries. Transport in Porous Media. 2007;70(2):191-211
13. 13. Bhadauria B, Siddheshwar P, Kumar J, Suthar OP. Weakly nonlinear stability analysis of temperature/gravity-modulated stationary rayleigh–bénard convection in a rotating porous medium. Transport in Porous Media. 2012;92(3):633-647
14. 14. Gresho P, Sani R. The effects of gravity modulation on the stability of a heated fluid layer. Journal of Fluid Mechanics. 1970;40(4):783-806
15. 15. Kiran P, Bhadauria B. Chaotic convection in a porous medium under temperature modulation. Transport in Porous Media. 2015;107(3):745-763
16. 16. Chandrasekhar S. Hydrodynamic and Hydromagnetic Stability. Oxford University Press; 1961
17. 17. Thompson W. Cxliii. Thermal convection in a magnetic field. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science. 1951;42(335):1417-1432
18. 18. Heris SZ, Salehi H, Noie S. The effect of magnetic field and nanofluid on thermal performance of two-phase closed thermosyphon (tpct). International Journal of Physical Sciences. 2012;7(4):534-543
19. 19. Ghasemi B, Aminossadati S, Raisi A. Magnetic field effect on natural convection in a nanofluid-filled square enclosure. International Journal of Thermal Sciences. 2011;50(9):1748-1756
20. 20. Hamad M, Pop I, Ismail AM. Magnetic field effects on free convection flow of a nanofluid past a vertical semi-infinite flat plate. Nonlinear Analysis: Real World Applications. 2011;12(3):1338-1346
21. 21. Bergman MI, Fearn DR. Chimneys on the earth’s inner-outer core boundary? Geophysical Research Letters. 1994;21(6):477-480
22. 22. Furmanski P, Banaszek J. Modelling of the mushy zone permeability for solidification of binary alloys. Materials Science Forum. Trans Tech Publ. 2006;508:411-418
23. 23. Gupta U, Ahuja J, Wanchoo R. Rayleigh-bénard convection of nanofluids with magnetic field and permeability effects. Procedia Engineering. 2015;127:325-332
24. 24. Imomnazarov KK. Modified darcy laws for conducting porous media. Mathematical and Computer Modelling. 2004;40(1–2):5-10
25. 25. Rana G, Kango S, Chadha K. Magneto-thermal convection in Walters’ (model b’) elastico-viscous fluid saturated by a Darcy-Brinkman porous medium. Engineering Mechanics. 2014;21:425-435
26. 26. Motsa S, Dlamini P, Khumalo M. A new multistage spectral relaxation method for solving chaotic initial value systems. Nonlinear Dynamics. 2013;72(1–2):265-283

Written By

Osman A.I. Noreldin, Precious Sibanda and Sabyasachi Mondal

Submitted: May 29th, 2017 Reviewed: September 19th, 2017 Published: December 20th, 2017