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.
(a)
(b)
(c)
(d)
(a)
(b)
(c)
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 .
(a)
(b)
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 , , here , , , and are functions of time since and are functions of . Thus and are positive functions. Substituting the definitions of and and , , into (107) and (108) we get, Equations (109) and (111) not contain phase term, so we take these two equations for the future discussions. We have (109) and (111) as since and are positive functions. Put Now we discuss the stability of equilibrium points of (114). We get four equilibrium points like (conduction state), ( amplitude of left travelling waves, here we get , and we get one condition from i.e., ), ( amplitude of right travelling waves, here we get , and we get one condition from i.e., ), and for and we get (this gives condition for standing waves. At standing waves we have , so ). For the pair of (104) and (105), we do not get (modulated waves). Now the Jacobian of and is given by
If real parts of all eigenvalues of the Jacobian are negative at an equilibrium point, then that point is a stable equilibrium (Lyapunov's theorem or principle of linearized stability). Some valuable conditions for travelling waves and standing waves are travelling waves are stable if and . Standing waves are stable if and (i) if , then and (ii) if , then .
The stability regions of travelling waves and standing waves are summarized in Figure 6. Here is total amplitude and defined as . We do not distinguish between left travelling waves and right travelling waves. For rest state (steady state) , for travelling waves , for standing waves . Travelling waves are supercritical if and standing waves are supercritical if . Figure 6(a) is drawn for stable travelling wave conditions and Figure 6(b) is drawn for stable standing wave conditions in -plane. The symbols and in Figures 6(a) and 6(b) indicate that both roots of Jacobian are negative and at least one root is positive between two roots. In Figures 6(a) and 6(b), travelling wave solution and standing wave solution bifurcate simultaneously from the steady-state solution ( at this bifurcation point).
(a)
(b)
(c)
In these Figures 6(a) and 6(b), steady-state solution is stable for and unstable . These figures show that for both travelling waves and standing waves are supercritical. When travelling waves and standing waves bifurcate supercritically then at most one solution among travelling waves and standing waves will be stable. Thus, for (Figure 6(a)) travelling waves are stable and (Figure 6(b)) standing waves are stable. In more detail we reproduce results of the stability analysis of equilibrium solutions in Figure 6(c), which is plotted in -plane. From this figure we can observe that travelling waves are subcritical for and standing waves are subcritical for . In Figure 7, We study the stability regions of travelling waves and standing waves at the onset of Hopf bifurcation. The stability regions of standing waves and travelling waves increases when increases for fixed parameters. For a fixed if we get initially travelling waves at the onset of oscillatory convection then they are replaced by standing waves as increases.
(a)
(b)
(c)
(d)
5.2. Long Wavelength Instabilities for the Onset of Travelling Wave Convection (Benjamin-Feir Instability)
For right travelling wave and , for left travelling wave and . Thus for travelling waves we get a single amplitude equation from (102) and (103), given as For standing waves and we get a single amplitude equation from (102) and (103), given as Equation (117) possesses a family of planar wave solutions and solutions containing phase singular points, which describes weakly nonlinear wave phenomena [25]. We study the Benjamin-Feir instability of travelling waves from complex Ginzburg- Landau equation (116) can be written as where . The phase winding solutions are obtained by substituting into (118), and equating real and imaginary parts we get Here is constant and . We consider a modulated solution in the form: . Substituting the modulated into (118) which gives It is possible to conduct a general investigation of the linear stability of , but this is very difficult task, and therefore our primary concern here is to treat the stability of the uniformly oscillating solution . Inserting into (120) and equating real and imaginary parts we get Consider and in the growth rate of disturbances. Using solutions of , and into (121) and (122) we get, solving (123) and (124), we get There will be an instability only when a root of (125) is possible, that is, when travelling waves or standing waves are stable. The instability of waves against long wavelength longitudinal modes is often called the Benjamin-Feir instability. Thus we get Benjamin-Feir instability for travelling waves when . Similarly by considering (118) instead of (117) and proceeding in the same way we get Benjamin-Feir instability for standing waves when .
6. Conclusions
In this paper we have considered both linear and weakly nonlinear analysis of magnetoconvection in a sparsely packed porous medium in Earth's outer core by using free-free (stress-free) boundary conditions. Even though free-free boundary conditions cannot be achieved in laboratory, one can use it in geophysical fluid dynamic applications to Earth's outer core since they allow simple trigonometric eigenfunctions. Our goal is to identify the region of parameter values, for which rolls emerge at the onset of convection.
Following Chandrasekhar [2], we have described the stationary convection and oscillatory convection as curves and versus wave numbers. The critical wave numbers for stationary convection and oscillatory convection are . For the problem of magnetoconvection in a sparsely packed porous medium, we get Takens-Bogdanov bifurcation point and codimension-two bifurcation point. In the case of linear theory both marginal and overstable motions are discussed. In Figures 1 and 2, is shown that the effect of Chandrasekhar number and porous parameter is to make the system more stable. By drawing stability boundaries in the Rayleigh number plane it is shown that the effect of magnetic field and porous parameter is to decrease the region of stabilities. In the nonlinear equation (60), gives the Takens-Bogdanov bifurcation point at and when , (60) is not valid. The pitchfork bifurcation is supercritical if subcritical if . and we get tricritical point if . We have obtained from (60), long wave length instabilities, namely, Eckhaus and Zigzag instabilities. From (60) which is valid only for , we have calculated Nusselt number and studied heat transport by convection. We have also derived two one-dimensional nonlinear coupled Ginzburg-Landau type equations, namely, (98) at the onset of oscillatory convection at supercritical Hopf bifurcation. We have computed stability regions of SW and TW at both Hopf bifurcation. The conditions for SW and TW are and or , respectively. TW exist if and they are supercritical if . SW exist if and SW are supercritical if . When both SW and TW are supercritical then at most one equilibrium solution is stable. At Takens-Bogdanov bifurcation point we get both TW and SW. By deriving one-dimensional Ginzburg-Landau equations with complex coefficients, namely, (116) and (117), we have shown the existence of Benjamin-Feir-type of instability for both TW and SW. Near the Takens-Bogdanov bifurcation point the conducting state becomes unstable against both stationary and oscillatory mode, that is, the real parts of two eigenvalues pass through zero simultaneously. This violates the assumption made for deriving amplitude equations (60) and (98). Instead a new equation, which is second order in time, has to be used near the Takens-Bogdanov bifurcation point.