International Journal of Geophysics

Volume 2011 (2011), Article ID 207123, 17 pages

http://dx.doi.org/10.1155/2011/207123

## Nonlinear Magnetoconvection in a Sparsely Packed Porous Medium

^{1}Department of Mathematics, National Institute of Technology Warangal, Warangal 506004, India^{2}Disha Institute of Management and Technology, Satya Vihar, Vidhan Sabha-Chandrakhuri Marg, Raipur 492101, India

Received 2 July 2011; Revised 6 September 2011; Accepted 7 September 2011

Academic Editor: Steve Milan

Copyright © 2011 A. Benerji Babu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Linear and weakly nonlinear properties of magnetoconvection in a sparsely packed porous medium are investigated. We have obtained the values of Takens-Bogdanov bifurcation points and codimension two bifurcation points by plotting graphs of neutral curves corresponding to stationary and oscillatory convection for different values of physical parameters relevant to magnetoconvection in a sparsely packed porous medium near a supercritical pitchfork bifurcation. We have derived a nonlinear two-dimensional Ginzburg-Landau equation with real coefficients by using Newell-Whitehead (1969) method. The effect of the parameter values on the stability mode is investigated and shown the occurrence of secondary instabilities namely, Eckhaus and Zigzag instabilities. We have studied Nessult number contribution at the onset of stationary convection. We have also derived two nonlinear one-dimensional coupled Ginzburg-Landau-type equations with complex coefficients near the onset of oscillatory convection at a supercritical Hopf bifurcation and discussed the stability regions of standing and travelling waves.

#### 1. Introduction

Magnetoconvection in a porous medium uniformly heated from below is of considerable interest in geophysical fluid dynamics, as this phenomena may occur within the mushy layer of Earth's outer core. Earth's outer core consists of molten iron and lighter alloying element, sulphur in its molten form. This lighter alloying element present in the liquid phase is released as the new iron freezes due to supercooling onto the solid inner core. Hence we get mushy layer near the inner core boundary where the problem becomes convective instability in a porous medium [1]. The effect of geomagnetic field on the magnetoconvection instability is of interest in geophysics, particular in the study of Earth's interior where the molten liquid Iron is electrically conducting, which can become convectively unstable as a result of differential diffusion.

Magnetoconvection in an electrically conducting fluid in a nonporous medium has been studied extensively [2–8]. However, magnetoconvection in a porous medium has not received any attention inspite of its application in geophysical fluid dynamics problems. Palm et al., [9] investigated Rayleigh-Benard convection problem in a porous medium. Brand and Steinberg [10, 11] investigated convecting instabilities in binary liquid in a porous medium; However, Plam et al. [9] and Brand et al. have made use of Darcy's law ( is replaced by where is the permeability of a porous medium. for nonporous medium is infinity). They have also not considered usual convective nonlinearity. It is well known that Darcy's law breaks down in situations where in other effects like viscous shear and inertia come into play. In fact Darcy's law is applicable to densely packed porous medium. An alternative to Darcy's equation is Brinkman equation and is of the form where is the fluid viscosity and is the effective fluid viscosity. Brinkman model is valid for a sparsely packed porous medium wherein there is more window fluid to flow so that the distortion of velocity give rise to the usual shear force. Lapwood [13] was the first to suggest the inclusion of convective term in the momentum equation and study the Rayleigh-Benard convection in a sparsely packed porous medium. Recently, Tagare and Benerji [14] have investigated the problem of nonlinear convection in a sparsely packed porous medium due to thermal and compositional buoyancy.

In this paper we investigate the problem of magnetoconvection in a sparsely packed porous medium. The multiplicity of control parameters makes this system an interesting one for the study of hydrodynamic stability, bifurcation and turbulence [15]. Rudraiah [16] and Rudraiah and Vortmeyer [17] have studied both linear and steady nonlinear magnetoconvection in a sparsely packed porous medium using Brinkman model but they have taken effective viscosity same as fluid viscosity . However, experiments show that the ratio of effective viscosity to fluid viscosity takes the value ranging from 0.5 to 10.9 [18]. In Section 2, we write basic dimensionless equations in Boussinesq approximation for magnetoconvection in a sparsely packed medium by using a momentum equation with effective viscosity different from fluid viscosity. In Section 3, we study linear stability analysis. In Section 4.1, by using multiple-scale analysis of Newell and Whitehead [19], we derive two-dimensional nonlinear Ginzburg-Landau equation in complex amplitude with real coefficients near the super critical pitchfork bifurcation. In Section 4.2, we show the occurrence of secondary instabilities such as Eckhaus instability and Zigzag instability. We have also considered the effect of Nusselt number on heat transport by magnetoconvection in a sparsely packed porous medium. In Section 5, we derive two nonlinear one-dimensional coupled Ginzburg-Landau type equations with complex coefficients near the onset of oscillatory convection at a supercritical Hopf bifurcation. Following Matthews and Rucklidge [20], we have dropped slow space dependence in and obtained two coupled ordinary differential equations in and and discussed the stability regions of travelling and standing waves. By obtaining a one-dimensional Ginzburg-Landau equation in complex amplitude with complex coefficients near a supercritical Hopf bifurcation, we have shown the condition for occurrence of Benjamin-Feir-type instability [21] for travelling and standing waves. In Section 6, we write conclusions of the paper.

#### 2. Basic Equations

We consider an electrically and thermally conducting fluid saturating an infinite horizontal layer of a sparsely packed isotropic porous medium of depth with a uniform magnetic field in the vertical -direction. This layer is heated from below, the upper and lower bounding surfaces of the layer are assumed to be stress-free. Physical properties of the fluid are assumed to be constant, except for the density in the buoyancy term, so that the Boussinesq approximation is valid. The temperature difference across the stress-free boundaries is and the flow in the sparsely packed porous medium is governed by the Darcy-Lapwood-Brinkman model. The relevant basic equations are The fluid density is described by where is thermal expansion coefficient and is mean fluid density. Here is pressure, is mean fluid velocity, is temperature, is magnetic field, is porosity, is acceleration due to gravity, is permeability of porous medium, is coefficient of effective fluid viscosity, is thermal diffusivity, is magnetic permeability, and is magnetic diffusivity. Equation (3) is known as Darcy-Lapwood-Brinkman equation and is valid for . Givler and Altobelli [18] shown that the range of varies from 0.5 to 10.9. is dimensionless heat capacity and is defined as the ratio of the effective heat capacity of the porous medium to the heat capacity of the fluid. In a nonporous medium, and and (3) reduces to Navier-Stokes equation. In this paper, for sparsely packed porous medium, we consider . The conduction state is characterized by and we take the temperature perturbation as . We use the scaling Here is thermal diffusion time in a porous medium. Using (6) and (8), we can write basic dimensionless equations for magnetoconvection in a porous medium as The dimensionless parameters required for the description of the motion are Rayleigh number , thermal Prandtl number , magnetic Prandtl number , Chandrasekhar number , and Darcy number . The Curl of (10) gives where vorticity , current and The Curl of (13) in turn gives, after use of (9), Now taking the scalar product of (12), (13), and (15) with , we get, Geophysically acceptable velocities of propagating instabilities corresponding to geometric scalar variations occur only (where instabilities develop in ohmic diffusion timescale ), and 5, when the turbulent is present in the Earth's outer core. In the case of the instabilities are extremely slow depending on the thermal diffusion timescale . Using (11), (18), and (16) can be reduced to a form where here

*Boundary Conditions*

We assume that fluid is contained between and , where corresponds to boundary of solid iron core with Earth's mushy layer and corresponds to boundary of Earth's mushy layer with Earth's outer liquid core. For perfectly conducting boundary with temperature, we have
Also the normal component of the velocity would vanish on , , that is,
However, there are two more conditions to be imposed on velocity depending on the nature of the surface. In this paper we consider free-free boundary conditions, that is, on surfaces the tangential stresses vanish, which is equivalent to
where is dynamic viscosity. Since vanishes for on , , it follows that on a free surface . Hence from equation of continuity we have on for all , . In this paper we have considered only the idealized stress-free conditions on the surface and vanishing of temperature fluctuations. Thus at . and its even derivatives vanish at and .

#### 3. Linear Stability Analysis

We perform a linear stability analysis of the problem by substituting into linearized version of (19) is , and obtaining an equation We consider stress-free boundary conditions, then on , for all . Thus we can assume .

Substituting and into (27), we get where and , from relation equation (30), .

##### 3.1. Stationary Convection ()

Substituting in (28), we get here is the value of the Rayleigh number for stationary convection. The minimum value of is obtained for . where The wave number is identical to that for the single component fluid, while the threshold for the onset of stationary convection at pitchfork bifurcation is given by (34) with , where . Thus the magnetic field inhibits the onset of stationary convection.

##### 3.2. Oscillatory Convection ()

For the oscillatory convection and from (28), will be complex. But the physical meaning of requires it to be real. The condition that is real implies that imaginary part of (28) is zero, that is, where and are given by (30) and (31). For oscillatory convection since , for oscillatory convection . For , (35) implies that is a double zero corresponding to Takens-Bogdanov bifurcation point. For oscillatory convection, we havewhere . A necessary condition for is However, this is not sufficient condition and one must have in addition At Takens-Bogdanov bifurcation point , , and is a double zero at where The Takens-Bogdanov bifurcation point occurs when the neutral curves for Hopf and pitchfork bifurcation meet and only a single wave number is present, namely, . If then for all the first instability to set in is an oscillatory convection. The asymptotic values of and for large Chandrasekhar number are From the monotonic dependence of and on , we may conclude that for , there exists a such that for the onset of first instability will be stationary convection at pitchfork bifurcation while for it will be oscillatory convection at Hopf bifurcation. and for , we have above condition (41) gives codimension-two bifurcation point. However, there is no simple formula to give at the codimension-two bifurcation point by assuming as an independent variable, such kind of interesting result is not available in Chandrasekhar [2]. In Figures 1 and 2, each solid line stands for stationary convection (pitchfork bifurcation) and dotted line stands for oscillatory convection (Hopf bifurcation). In Figures 1 and 2, we have showed the effect of several physical parameters, like , , , , , , and on the onset of both stationary convection and oscillatory convection when a physical parameter increases for the remaining fixed parameters, the onset of instabilities increases, that is, the onset of stationary convection and oscillatory convection inhibit when a parameter increases with the remaining fixed parameters.

#### 4. Onset of Stationary Convection at Supercritical Pitchfork Bifurcation

##### 4.1. Derivation of Two-Dimensional Nonlinear Ginzburg-Landau Equation Using Newell-Whitehead [19] Method

In this section the evolution of a general pattern is developed by means of a multiple scale analysis used by Newell and Whitehead [19]. A small amplitude convection cell is imposed on the basic flow. If this amplitude is of the size then the interaction of the cell with itself forces a second harmonic and mean state correction of size and then in turn drives an correction to the fundamental component of the imposed roll. A solvability criteria for this correction yields the one-dimensional nonlinear Ginzburg-Landau equation of the complex valued amplitude of the imposed disturbance with real coefficients. To simplify the problem we assume the formulation of cylindrical rolls with axis parallel to -axis, so that -dependence disappears from (19). The -dependence is contained entirely in the and functions, which ensures that stress-free boundary conditions are satisfied. We use the expansion parameter as For the values of close to threshold value that is, , the structure of the slow length scales will be insensitive to , but a slow modulation in space and time is possible by making use of the band of the unstable solutions and linear growth rate is likely to saturate due to nonlinear effects. This behavior can be analyzed by writing solutions of (9)–(12) in power series as where with the first approximation is given by the eigenvector of the linearized problem: where , here C.C. stands for complex conjugate, is the critical mode for the linear problem at and . The complex amplitude depends on the slow variables , , , and to be scaled by introducing multiple scales and these formally separate the fast and slow dependent variables in . It should be noted that the difference in scaling in the two directions reflects the inherent symmetry breaking of instability which was chosen here with wave vector in -direction. The differential operators can be expressed as with the assumption (46), the operators (20) and (21) are transformed into a set of linear in homogeneous equations. The solvability conditions for the latter yields the amplitude equation using (44) in the linear operator (20) can be written as where Similarly nonlinear term is given by substituting (47), (51), and (43) into (19), we get by equating the coefficients of ; Equation (48) gives the critical Rayleigh number for the onset of stationary convection In (53), , and hence . From equation of continuity we find that . The relevant equations for and are form (56) and (44), we get Equation (12) gives relevant equation for as From (58) and (44), we get Similarly we have . The solvability criterion of (54) gives the amplitude equation which can be written as where Equation (60) is two-dimensional, nonlinear time-dependent Ginzburg-Landau equation describing the effect of magneticfield in a sparsely packed porous medium near the onset of stationary convection at supercritical pitchfork bifurcation. Here is always positive for and for any but if then is positive only if . Thus for supercritical pitchfork bifurcation is always positive. For , decreases as increases and becomes zero at . and are always positive. is positive only if The pitchfork bifurcation is supercritical if and subcritical if . At , we get tricritical bifurcation point [22] (see Figure 3). Dropping the time dependence from (60), we get since , the solution of (63) is given by where

##### 4.2. Long Wavelength Instabilities (Secondary Instabilities)

The secondary Instabilities arising in nonequilibrium systems do not exhibit strict symmetries but may show spatially slow deformations of the cellular structures. Further, there are secondary instabilities like Eckhaus and Zigzag instabilities, such phenomena are studied using evolution equations for amplitudes which are slowly varying in time as well as in space. These envelope equations can be derived by the method of Newell and Whitehead [19]. The two-dimensional Ginzburg-Landau equation (60), can be written in fast variables , and , as In order to study the properties of a structure with a given phase winding number , we substitute into (66) and we obtain The steady-state uniform solution of (68) is Let be an infinitesimal perturbation from a uniform steady-state solution given by (69). Now substituting into (68) and equating real and imaginary parts, we obtain We analyze (71) by using normal modes of the form Putting (72) in (71) we get, Here , . On solving (73) we get, whose roots are real. Here is defined as Solution is clearly negative, thus the corresponding mode is stable and if is positive then rolls can be unstable. Symmetry considerations help us to restrict the study of to a domain .

###### 4.2.1. Longitudinal Perturbations and Eckhaus Instability

Inserting into (75), we get since the roots are real and their sum always negative, the pattern is stable as long as both roots are negative, that is, their product is positive. The cell pattern becomes unstable when the product is negative, that is, when for this requires , this condition defines the domain of Eckhaus instability. The above condition implies that the most unstable wave vector tends to zero, when .

###### 4.2.2. Transverse Perturbations and Zigzag Instability

Let us consider into (75), we get where . The two eigenmodes are uncoupled and we have , for one of them. The other is amplified when This implies that , the above condition defines the domain of the Zigzag Instability. When from below the wave vector of the instability also tends to zero while the growth rate varies as . We have studied the effect of magnetic field on long wavelength instabilities. We have observed that Eckhaus instability and Zigzag instability regions increases when increases (see Figure 4).

##### 4.3. Heat Transport by Convection

The maximum of steady amplitude is denoted by which is given as Equation (81) is obtained from (64) with . We use to calculate Nusselt number . To discuss the heat transfer near the neutral region, we express it through the Nusselt number is defined as , which is the ratio of the heat transported across any layer to the heat which would be transported by conduction alone. Here is the rate of heat transfer per unit area and is defined as In (82), angular brackets correspond to a horizontal average. The Nusselt number can be calculated in terms of amplitude and is given as From (83), we get conduction for and convection for . Since the amplitude equation is valid for , which is possible for (supercritical pitchfork bifurcation), thus we get for . We get convection for and conduction for . In stationary convection increases implies that heat conducted by steady mode increases. In the problem of double diffusive convection in porous medium with magnetic field, depends on , and . We have computed for different values of , for some fixed values of other parameters and observed that increases as decreases (see Figures 5(a) and 5(b)). This implies that magnetic field inhibits the heat transport. The parameters , , , , , and show the same result as shows on .

#### 5. Oscillatory Convection at the Supercritical Hopf Bifurcation

The existence of a threshold (critical value of Rayleigh number for the onset of oscillatory convection ) and a cellular structure (critical wave number ) are main characteristics of the oscillatory convection. In this section we treat the region near the onset of oscillatory convection. Here the axis of cylindrical rolls is taken as -axis, so that -dependence disappears from equation . The -dependence contained entirely in and functions which ensure that the free-free boundary conditions are satisfied. The purpose of this section is to derive coupled one-dimensional nonlinear time-dependent Ginzburg-Landau type equations near the onset of oscillatory convection at supercritical Hopf bifurcation. We introduce as We assume that is a solution to linearized equation , which satisfies free-free boundary conditions. Here denotes the amplitude of left travelling wave of the roll and denotes the amplitude of right travelling wave of the roll, which depends on slow space and time variables [23] and assume that , . The differential operators can be expressed as The solution of basic equations can be sought as power series in , where with the first approximation given by eigenvector of the linearized problem: where , , and , here and are complex conjugate of and .

We expand the linear operator and nonlinear term as the following power series substituting (87) and (88) into , for each order of , we get Here where Equation (91) is linear problem. We get critical Rayleigh number for the onset of oscillatory convection by using the zeroth-order solution in (91). At and gives where is the group velocity and is real. Hence from (92), we get . From equation of continuity we find that . Substituting the zeroth-order and first-order approximation into (56) and (58) we get, where , , , and and , and are complex conjugate of and , respectively.

Equation (93) is solvable when , one requires that its right-hand side be orthogonal to , which is ensured that if the coefficients of in are equal to zero. This implies that where Here . It should be noted that and are of order and and are of order . If in , and then these expressions match with the coefficients , , , and of Ginzburg-Landau equation at the onset of stationary convection. From (96), we get and , where . Equations (98) can be written as

Let , where are periods of , respectively. Expansion (88) remains asymptotic for times only if an appropriate solvability condition holds. This condition obtained integrating (100) over and (101) over , we get

##### 5.1. Travelling Wave and Standing Wave Convection

To study the stability regions of travelling waves and standing waves, Coullet et al. [24]. we proceed as follows.

On dropping slow variable from (102) and (103), we get a pair of first ODE's Put Then (104) and (105) take the following form Consider and (we can write a complex number in the amplitude and phase form), where , and