IJEM International Journal of Engineering Mathematics 2314-6109 2356-7007 Hindawi Publishing Corporation 10.1155/2014/296216 296216 Research Article Weak Nonlinear Double-Diffusive Magnetoconvection in a Newtonian Liquid under Temperature Modulation Bhadauria B. S. Kiran Palle Liu Yurong Department of Applied Mathematics School for Physical Sciences Babasaheb Bhimrao Ambedkar University Lucknow 226 025 India bbau.ac.in 2014 6 7 2014 2014 22 02 2014 15 06 2014 6 7 2014 2014 Copyright © 2014 B. S. Bhadauria and Palle Kiran. 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.

The present paper deals with a weak nonlinear theory of double-diffusive magnetoconvection in an electrically conducting Newtonian liquid, confined between two horizontal surfaces, under a constant vertical magnetic field, and subjected to imposed time-periodic thermal boundaries. The temperature of both walls is varied time periodic in this case. The disturbances are expanded in terms of power series of amplitude of convection, which is assumed to be small. Using nonautonomous Ginzburg-Landau equation, the Nusselt and Sherwood numbers obtained analytically and studied heat and mass transport in the system. Effect of various parameters on the heat and mass transport is discussed extensively. It is found that the effect of magnetic field is to stabilize the system. Further, it is also notified that the heat and mass transport can be controlled by suitably adjusting the external parameters of the system.

1. Introduction

Double-diffusive convection is an important fluid dynamics phenomenon that involves motions driven by two different density gradients diffusing at different rates. In double-diffusive convection, the buoyancy force is affected not only by the difference of temperatures but also by the difference of concentration of the fluid. An example of double-diffusive convection is seen in oceanography, lakes, underground water, atmospheric pollution, chemical processes, laboratory experiments, modeling of solar ponds , electrochemistry, magma chambers and sparks , Fernando and Brandt, formation of microstructure during the cooling of molten metals, fluid flows around shrouded heat-dissipation fins, migration of moisture through air contained in fibrous insulations, grain storage system, the dispersion of contaminants through water saturated soil, crystal P growth, solidification of binary mixtures, and the underground disposal of nuclear wastes. Therefore, much work has also been done on double-diffusive convection in an electrically conducted fluid layer because of its natural occurrence as mentioned above applications. Convection in planetary cores, stellar interiors, and Earth’s metallic core occurs in the presence of strong magnetic field. The study of double-diffusive magnetoconvection has recently drawn the attention of astrophysicists, geophysicists, oceanographers, engineers, and a host of others [4, 5]. The study of magnetoconvection in an electrically conducting horizontal fluid layer was motivated by astrophysical and geophysical applications; relate in some or the other way to the problems concerning the external constraints like rotation or magnetic field operative on double-diffusive systems, in particular by observation of sunspots .

Thompson  and Chandrasekhar  were the first to study the magnetoconvection in horizontal fluid layer. Lortz  was the first to study the effect of magnetic field on double-diffusive convection. His object was to clarify some of the mathematical aspects of stability criterion  but his analysis is silent about the detailed study of stability analysis. Stommel et al.  explained that the diffusion is generally a stabilizing factor in a single-component fluid. But in the case of two-component system it can act to release the potential energy in the component that is the heaviest at the top and make the system unstable. Gotoh and Yamada  studied the problem of magnetoconvection in a horizontal layer of magnetic fluid which is heated from below and cooled from above and found condition for onset of convection. Oreper and Szekely  have found that the presence of a magnetic field can suppress natural convection currents and that the strength of the magnetic field is one of the important factors in determining the quality of the crystal. Rudraiah and Shivakumara  investigated both linear and nonlinear theory of this problem in detail. They have shown that the magnetic field, under certain conditions, makes the system unstable. They have also investigated the effect of magnetic field on the physically preferred cell pattern. Rudraiah  investigated the interaction between double-diffusive convection and an externally imposed vertical magnetic field in a Boussinesq fluid. Siddheshwar and Pranesh  analyzed the role of magnetic field in the inhibition of natural convection driven by combined buoyancy and surface tension forces in a horizontal layer of an electrically conducting Boussinesq fluid with suspended particles confined between an upper free/adiabatic and a lower rigid/isothermal boundary is considered in 1 g and μg situations. Bhadauria  also studied the effect of magnetic field on thermal modulated convection in the case of porous medium. Siddheshwar et al.  performed a local nonlinear stability analysis of Rayleigh-Bénard magnetoconvection using Ginzburg-Landau equation. They showed that gravity modulation can be used to enhance or diminish the heat transport in stationary magnetoconvection.

The classical Rayleigh-Bénard convection due to bottom heating is well known and highly explored phenomenon given by Chandrasekhar  and Drazin and Reid . Many researchers, under different physical models, have investigated thermal instability in a horizontal fluid layer with temperature modulation in the absence of double-diffusive magnetoconvection. Some of them are Venezian , who was the first to consider the effect of temperature modulation on thermal instability in a horizontal fluid layer. Rosenblat and Tanaka  studied the linear stability for a fluid in a classical geometry of Bénard by considering the temperature modulation of rigid-rigid boundaries. The first nonlinear stability problem in a horizontal fluid layer, under temperature modulation of the boundaries, was studied by Roppo et al. . Bhadauria and Bhatia  studied the effect of temperature modulation on thermal instability by considering rigid-rigid boundaries and different types of temperature profiles. Bhadauria , Malashetty and Swamy , and Bhadauria and Kiran  are the related problems.

The problem of double-diffusive magnetoconvection is called thermohaline magnetoconvection when the two diffusive mechanisms are thermal and solute. From the above literature (second paragraph) it is observed that concerning linear theory numerous studies on magnetoconvection or double-diffusive magnetoconvection have been made. But, the double-diffusive magnetoconvection is still in its infant stage in nonlinear case. Liner theory gives information about onset of convection but fails at heat and mass transport. Most of the studies are made with uniform temperature gradient across the boundaries. In general there are many practical situations in which the temperature gradient is a function of both time and space. One of the external effective mechanisms to regulate (or control) the convection by maintaining a nonuniform temperature gradient across the fluid layer. We also observe from the literature (third paragraph) that some work has been done under temperature modulation without magnetic field and double-diffusive convection which is essential in thermal and engineering science as mentioned in the literature. Therefore, from the above understanding we motivated to analyze a weak nonlinear thermohaline magnetoconvection under temperature modulation while deriving an amplitude of convection using Ginzburg-Landau model.

2. Governing Equations

We consider an electrically conducting fluid layer of depth d , confined between two infinite parallel horizontal planes at z = 0 and z = d . A Cartesian frame of reference has been taken with the origin at the bottom, and the z axis vertically upward direction of the fluid layer. The surfaces are maintained at a constant temperature gradient Δ T / d and a constant magnetic field H b k is applied across the fluid layer (as shown in Figure 1). Under the Boussinesq approximation, the dimensional governing equations for the study of double-diffusive magnetoconvection in an electrically conducting fluid layer are defined as (1) · q = 0 , (2) · H = 0 , (3) q t + ( q · ) q = 1 ρ 0 p + ρ ρ 0 g - μ ρ 0 2 q + μ m ρ 0 H · H , (4) γ T t + ( q · ) T = κ T 2 T , (5) S t + ( q · ) S = κ S 2 S , (6) H t + ( q · ) H - ( H · ) q = ν m 2 H , (7) ρ = ρ 0 [ 1 - α ( T - T 0 ) + β ( S - S 0 ) ] , where q is velocity ( u , v , w ) , μ is a viscosity, μ m is magnetic permeability, ν m is magnetic viscosity, κ T is the thermal diffusivity, κ S is the solutal diffusivity, T is temperature, α is thermal expansion coefficient, β is solute expansion coefficient, and γ is the ration of heat capacity. For simplicity γ is taken to be unity in this paper, ρ is density, g = ( 0,0 , - g ) is gravitational acceleration, while ρ 0 is reference density, and H is the magnetic field. The externally imposed thermal boundary conditions considered in this paper are given by Venezian : (8) T = T 0 + Δ T 2 [ 1 + ϵ 2 δ cos ( ω t ) ] at z = 0 = T 0 - Δ T 2 [ 1 - ϵ 2 δ cos ( ω t + θ ) ] at z = d , where δ is small amplitude of modulation, Δ T is the temperature difference across the fluid layer, ω is the modulation frequency, and θ is the phase difference. Since we are not considering cross diffusion terms, the walls of the fluid layer are assumed to be maintained at constant solute concentration as defined below: (9) S = S 0 + Δ S at z = 0 , = S 0 at z = d . The basic state is assumed to be quiescent and the quantities in this state are given by (10) q b = 0 , ρ = ρ b ( z , t ) , p = p b ( z , t ) , T = T b ( z , t ) , S = S b ( z ) . The basic state pressure field is not required here; however the basic temperature and solute fields are governed by the following ordinary differential equations: (11) T b t = κ T 2 T b z 2 , 2 S b z 2 = 0 . The solutions of (11) subjected to the boundary conditions (8)-(9) are given by (12) T b ( z , t ) = T s ( z ) + ϵ 2 δ Re [ T 1 ( z , t ) ] , S b = S 0 + Δ S ( 1 - z d ) , where T s ( z ) is the study temperature field and T 1 ( z , t ) is the oscillating part while Re stands for the real part. We assume finite amplitude perturbations on the basic state in the form (13) q = q b + q , ρ = ρ b + ρ , S = S b + S , p = p b + p , T = T b + T , H = H b + H , where the primes denote the quantities at the perturbations. Substituting (13) into (1)–(6) and using the basic state results, we obtain the following equations: (14) · q = 0 , · H = 0 , q t + ( q · ) q - μ m ρ 0 ( H · ) H = 1 ρ 0 p - ρ ρ 0 g k + ν 2 q + μ m ρ 0 H b H z , T t + ( q · ) T + w T b z = κ T 2 T , S t + ( q · ) S + w S b z = κ S 2 S , H t + ( q · ) H - ( H · ) q - H b w z = ν m 2 H , ρ = - ρ 0 ( α T - β S ) . Further we consider only two-dimensional disturbances in our study and hence the stream function ψ and magnetic potential Φ are introduced as ( u , w ) = ( ψ / z , - ψ / x ) and ( H x , H z ) = ( Φ / z , - Φ / x ) . By eliminating density and pressure terms from (14) and nondimensionalizing using the following transformations: ( x , y , z ) = d ( x * , y * , z * ) , ψ = κ T ψ * , t = ( d 2 / κ T ) t * , q = ( κ T / d ) q * , T = Δ T    T * ,   S = Δ S    S * , Φ = d H b Φ * , and ω = ( κ T / d 2 ) ω * , finally dropping the asterisk for simplicity, we get the nondimensionalized governing system of equations: (15) - 4 ψ + R a T T x - Ras S x - Q Pm z ( 2 Φ ) = - 1 Pr t ( 2 ψ ) + 1 Pr ( ψ , 2 ψ ) ( x , z ) - Q Pm ( Φ , 2 Φ ) ( x , z ) , (16) - ψ x T b z - 2 T = - T t + ϵ 2 δ f 2 ( z , t ) ψ x + ( ψ , T ) ( x , z ) , (17) - ψ x S b z - 1 Le 2 S = - S t + ( ψ , S ) ( x , z ) , (18) - ψ z - Pm 2 Φ = - Φ t + ( ψ , Φ ) ( x , z ) . Here the nondimensionalizing parameters in the above equations are Pm = ν m / κ T magnetic Prandtl number, Pr = ν / κ T Prandtl number, Le = κ T / κ S is Lewis number, R a T = α g Δ T d 3 / ν κ T is thermal Rayleigh number, Ras = β g Δ S d 3 / ν κ T is Solutal Rayleigh number, and Q = μ m H b 2 d 2 / ρ 0 ν ν m is Chandrasekhar number. In (16), one can observe that the basic state solution influences the stability problem through the factor T b / z , which is given by (19) T b z = - 1 + ϵ 2 δ [ f 2 ( z , t ) ] , where (20) f 2 ( z , t ) = Re [ f ( z ) e - i ω t ] , f ( z ) = [ A ( λ ) e λ z + A ( - λ ) e - λ z ] , A ( λ ) = ( λ / 2 ) ( ( e - i θ - e - λ ) / ( e λ - e - λ ) ) , and λ = ( 1 - i ) ω / 2 . We assume small variations of time and rescaling it as τ = ϵ 2 t , to study the stationary convection of the system, we write the nonlinear Equations (15)–(18) in the matrix form as given below: (21) [ - 4 R a T x - Ras x - Q Pm z 2 x - 2 0 0 x 0 - 1 Le 2 0 - z 0 0 - Pm 2 ] [ ψ T S Φ ] = [ - 1 Pr t ( 2 ψ ) + 1 Pr ( ψ , 2 ψ ) ( x , z ) - Q Pm ( Φ , 2 Φ ) ( x , z ) - T t + ϵ 2 δ f 2 ( z , t ) ψ x + ( ψ , T ) ( x , z ) - S t + ( ψ , S ) ( x , z ) - Φ t + ( ψ , Φ ) ( x , z ) ] . The considered stress free and isothermal boundary conditions to solve the above system (21) are (22) ψ = 0 = 2 ψ , Φ = D Φ = 0 at z = 0 , ψ = 0 = 2 ψ , Φ = D Φ = 0 at z = 1 , where D = / z .

Physical configuration of the problem.

3. Finite Amplitude Equation, Heat and Mass Transport

We now introduce the following asymptotic expansions in (21): (23) R a T = R 0 c + ϵ 2 R 2 + ϵ 4 R 4 + , ψ = ϵ ψ 1 + ϵ 2 ψ 2 + ϵ 3 ψ 3 + , T = ϵ T 1 + ϵ 2 T 2 + ϵ 3 T 3 + , S = ϵ S 1 + ϵ 2 S 2 + ϵ 3 S 3 + , Φ = ϵ Φ 1 + ϵ 2 Φ 2 + ϵ 3 Φ 3 + , where R 0 c is the critical value of the Rayleigh number at which the onset of convection takes place in the absence of temperature modulation. Now we solve the system for different orders of ϵ .

At the lowest order, we have (24) [ - 4 R 0 c x - Ras x - Q Pm z ( 2 ) x - 2 0 0 x 0 - 1 Le 2 0 - z 0 0 - Pm 2 ] [ ψ 1 T 1 S 1 Φ 1 ] = [ 0 0 0 0 ] . The solutions of the lowest order of the system subject to the boundary conditions equation (22) are (25) ψ 1 = A ( τ ) sin ( k c x ) sin ( π z ) , T 1 = - k c δ 2 A ( τ ) cos ( k c x ) sin ( π z ) , S 1 = - k c δ 2 Le A ( τ ) cos ( k c x ) sin ( π z ) , Φ 1 = π Pm δ 2 A ( τ ) sin ( k c x ) cos ( π z ) , where δ 2 = k c 2 + π 2 . The critical value of the Rayleigh number for the onset of magnetoconvection in the absence of temperature modulation is (26) R 0 c = δ 2 ( δ 4 + Q π 2 ) + Ras k c 2 Le k c 2 , which is the same as Siddheshwar et al.  when Ras = 0 and we obtain classical results of Chandrasekhar  for without magnetic field and single component fluid layer.

At the second order, we have (27) [ - 4 R 0 c x - Ras x - Q Pm z 2 x - 2 0 0 x 0 - 1 Le 2 0 - z 0 0 - Pm 2 ] [ ψ 2 T 2 S 2 Φ 2 ] = [ R 21 R 22 R 23 R 24 ] , R 21 = 0 , R 22 = ψ 1 x T 1 z - ψ 1 z T 1 x , R 23 = ψ 1 x S 1 z - ψ 1 z S 1 x , R 24 = ψ 1 x Φ 1 z - ψ 1 z Φ 1 x . The second order solutions subjected to the boundary conditions Equation (22) are obtained as follows: (28) ψ 2 = 0 , T 2 = - k c 2 8 π δ 2 A 2 ( τ ) sin ( 2 π z ) , S 2 = - k c 2 L e 2 8 π δ 2 A 2 ( τ ) sin ( 2 π z ) , Φ 2 = - π 2 8 k c P m 2 δ 2 A 2 ( τ ) sin ( 2 k c x ) . The horizontally averaged Nusselt Nu and Sherwood Sh numbers for the stationary double-diffusive magnetoconvection (the mode considered in this problem) are given by (29) Nu ( τ ) = 1 + [ ( k c / 2 π ) 0 2 π / k c ( T 2 / z ) d x ] z = 0 [ ( k c / 2 π ) 0 2 π / k c ( T b / z ) d x ] z = 0 , Sh ( τ ) = 1 + [ ( k c / 2 π ) 0 2 π / k c ( S 2 / z ) d x ] z = 0 [ ( k c / 2 π ) 0 2 π / k c ( d S b / d z ) d x ] z = 0 . One must note here that f 2 ( z , τ ) is effective at O ( ϵ 2 ) and affects Nu ( τ ) and Sh ( τ ) through A ( τ ) can be seen later. Therefore, substituting T 2 , T b , S 2 , respectively, into (29) and simplifying we obtain (30) Nu ( τ ) = 1 + k c 2 4 δ 2 A 2 ( τ ) , Sh ( τ ) = 1 + k c 2 L e 2 4 δ 2 A 2 ( τ ) . At the third order, we have (31) [ - 4 R 0 c x - Ras x - Q Pm z 2 x - 2 0 0 x 0 - 1 Le 2 0 - z 0 0 - Pm 2 ] [ ψ 3 T 3 S 3 Φ 3 ] = [ R 31 R 32 R 33 R 34 ] , where (32) R 31 = - 1 Pr τ ( 2 ψ 1 ) - R 2 T 1 x - Ras S 1 x - Q Pm ( Φ 1 z x ( 2 Φ 2 ) - Φ 2 x z ( 2 Φ 1 ) ) , R 32 = - T 1 τ + δ f 2 ( z , τ ) ψ 1 x + ψ 1 x T 2 z , R 33 = - S 1 τ + ψ 1 x S 2 z , R 34 = - Φ 1 τ - ψ 1 z Φ 2 x . Substituting ψ 1 , T 1 , and T 2 into (32), we obtain the expressions for R 31 , R 32 , R 33 , and R 34 easily. For the existence of third order solution of the system we apply the solvability condition which leads to arrive at the nonautonomous Ginzburg-Landau equation for stationary mode of convection with time-periodic coefficients in the form (33) A 1 A ( τ ) - A 2 A ( τ ) + A 3 A ( τ ) 3 = 0 , where (34) A 1 = [ δ 2 Pr + R 0 c k c 2 δ 4 - Ras k c 2 L e 2 δ 4 - Q π 2 Pm δ 2 ] , A 2 = [ R 2 k c 2 δ 2 - 2 R 0 c k c 2 δ 2 δ I 1 ] , A 3 = [ R 0 c k c 4 8 δ 4 + Q π 4 k c 2 2 P m 2 δ 4 - Q π 4 4 P m 2 δ 2 - Ras k c 4 L e 3 8 δ 4 ] , I 1 = 0 1 f 2 ( z , τ ) sin 2 ( π z ) d z . The Ginzburg-Landau equations given in (33) are Bernoulli equation and obtaining its analytical solution is difficult, due to its nonautonomous nature. So that, it has been solved numerically using the in-built function NDSolve of Mathematica 8, subjected to the initial condition A ( 0 ) = a 0 , where a 0 is the chosen initial amplitude of convection. In our calculations we may use R 2 = R 0 c , to keep the parameters to the minimum.

4. Results and Discussion

External regulation of thermal instability is important to study the double-diffusive convection in a fluid layer. The objective of this paper is to consider two such candidates, namely, vertical magnetic field and temperature modulation for either enhancing or inhibiting convective heat and mass transport as is required by a real application. The present paper deals with double-diffusive magnetoconvection under temperature modulation by using Ginzburg-Landau equation. It is necessary to consider a nonlinear theory to analyze heat and mass transfer which is not possible by the linear theory. We consider the direct mode ( κ S / κ T < 1 otherwise Hopf mode ) in which the salt and heat make opposing contributions ( κ T κ S ) . We also consider the effect of temperature modulation to be of order O ( ϵ 2 ) ; this leads to small amplitude of modulation. Such an assumption will help us in obtaining the amplitude equation of convection in a rather simple and elegant manner and is much easier to obtain than in the case of the Lorenz model. We give the following features of the problem before our results:

the need for nonlinear stability analysis,

the relation of the problem to a real application,

the selection of all dimensionless parameters utilized in computations.

We consider the following three types of temperature modulation on the boundaries of the system:

in-phase modulation [IPM] ( θ = 0 ) ,

out-of-phase modulation [OPM] ( θ = π ) ,

lower boundary modulation [LBMO] ( θ = - i ) .

The parameters of the system are Q , Pr , Le , Pm , Ras , θ , δ , ω ; these parameters influence the convective heat and mass transfer. The first five parameters are related to the fluid layer, and the last three parameters concern the external mechanisms of controlling convection. The effect of temperature modulation is represented by amplitude δ , which lies around 0.3 due to the assumption. The effect of electrical conductivity and magnetic field comes through Pm , Q . There is the property of the fluid coming into picture as well as through Prandtl number Pr . Further, the modulation of the boundary temperature is assumed to be of low frequency. At low range of frequencies the effect of frequency on onset of convection as well as on heat and mass transport is minimal. This assumption is required in order to ensure that the system does not pick up oscillatory convective mode at onset due to modulation in a situation that is conductive otherwise to stationary mode. It is important at this stage to consider the effect of Q , Pr , Le , Pm , Ras , θ , δ , ω on the onset of convection. The heat and mass transfer of the problem are quantified by the Nusselt and Sherwood numbers which are given in (30). Figures 27 show the individual effect of each nondimensional parameter on heat and mass transfer.

Nu versus τ (a) Q (b) Pr (c) Le (d) Pm (e) Ras (f) δ , ω .

Nu versus τ (a) Q (b) Pr (c) Le (d) Pm (e) Ras (f) δ (g) ω .

Nu versus τ (a) Q (b) Pr (c) Le (d) Pm (e) Ras (f) δ (g) ω .

Nu versus τ (a) Q (b) Pr (c) Le (d) Pm (e) Ras (f) δ , ω .

Sh versus τ (a) Q (b) Pr (c) Le (d) Pm (e) Ras (f) δ (g) ω .

Sh versus τ (a) Q (b) Pr (c) Le (d) Pm (e) Ras (f) δ (g) ω .

(1) Figures 2(a)7(a) show that the effect of Chandrasekhar number Q which is ratio of Lorentz force to viscous force is to delay the onset of convection, hence heat and mass transfer. The Nu and Sh start with one and for small values of time τ increase and become constant for large values of time τ in the case of [IPM] given in Figures 2(a) and 5(a). In the case of [OPM, LBMO] the effect of Q shows oscillatory behaviour and increment in it decreases the magnitude of both Nu and Sh . Hence Q has stabilizing effect in all the three types of modulations given in Figures 3(a), 4(a), 6(a), and 7(a), so that heat and mass transfer decrease with Q .

(2) The effect of Prandtl number Pr is to advance the convection and hence heat and mass transfer, which has similar behaviour of (1), given in Figures 2(b)7(b).

(3) The effect of Lewis Le and Magnetic Prandtl Pm numbers is to advance the convection and hence heat and mass transfer. Hence both Le and Pm have destabilizing effect of the system given in Figures 2(c)7(c) and 2(d)7(d) and have similar behaviour of (1); these results are earlier obtained by Siddheshwar et al. .

(4) The effect of solutal Rayleigh number Ras is to increase Nu and Sh so that heat and mass transfer. Hence it has destabilizing effect in all the three types of modulations which is given by the Figures 2(e)7(e) and has similar behaviour of (1). Though the presence of a stabilizing gradient of solute will prevent the onset of convection, the strong finite-amplitude motions, which exist for large Rayleigh numbers, tend to mix the solute and redistribute it so that the interior layers of the fluid are more neutrally stratified. As a consequence, the inhibiting effect of the solute gradient is greatly reduced and hence fluid will convect more and more heat and mass when Ras is increased.

(5) In the case of [IPM] we observe no effect of amplitude δ and frequency ω of modulation which is given by the Figures 2(f) and 5(f). But in the case of [OPM-LBMO] the increment in δ leads to increment in magnitude of Nu and Sh hence heat and mass transfer given in Figures 3(f), 4(f), 6(f), and 7(f); the increment in ω shortens the wavelength and decreases in magnitude of Nu , Sh and hence heat and mass transfer given in Figures 3(g), 4(g), 6(g), and 7(g) which are the results obtained by Venezian , Siddheshwar et al. , and Bhadauria and Kiran .

(6) From Figures 2(g) and 5(g) we observe that Q has strongly stabilizing effect. [ Nu / Sh ] Q = 0 > [ Nu / Sh ] Q 0 .

(7) The comparison of three types of temperature modulations is given in Figures 4(h) and 7(h): [ Nu / Sh ] IPM < [ Nu / Sh ] LBMO < [ Nu / Sh ] OPM .

(8) The results of this work can be summarized as follows from Figures 27.

[ Nu / Sh ] Q = 25 < [ Nu / Sh ] Q = 15 < [ Nu / Sh ] Q = 10 , Figures 2(a)7(a).

[ Nu / Sh ] Pr = 0.5 < [ Nu / Sh ] Pr = 1.0 < [ Nu / Sh ] Pr = 1.5 , Figures 2(b)7(b).

[ Nu ] Le = 1.2 < [ Nu ] Le = 4.2 < [ Nu ] Le = 6.2 , Figures 2(c), 3(c), and 4(c).

[ Sh ] Le = 1.2 < [ Sh ] Le = 1.4 < [ Sh ] Le = 1.6 , Figures 5(c), 6(c), and 7(c).

[ Nu / Sh ] Pm = 1.2 < [ Nu / Sh ] Pm = 1.4 < [ Nu / Sh ] Pm = 1.6 , Figures 2(d)7(d).

[ Nu / Sh ] Ras = 20 < [ Nu / Sh ] Ras = 50 < [ Nu / Sh ] Ras = 100 , Figures 2(e)7(e).

[ Nu / Sh ] δ = 0.1 < [ Nu / Sh ] δ = 0.3 < [ Nu / Sh ] δ = 0.5 , Figures 3(f), 4(f), 6(f), and 7(f).

[ Nu / Sh ] ω = 100 < [ Nu / Sh ] ω = 50 < [ Nu / Sh ] ω = 30 < [ Nu / Sh ] ω = 5 , Figures 3(g), 4(g), 6(g), and 7(g).

5. Conclusions

The effect of temperature modulation on weak nonlinear double-diffusive magnetoconvection has been analyzed by using the nonautonomous Ginzburg-Landau equation. The following conclusions are drawn from previous analysis.

The effect of IPM is negligible on heat and mass transport in the system.

In the case of IPM, the effect of δ and ω is also found to be negligible on heat and mass transport.

In the case of IPM, the values of Nu and Sh increase steadily for small values of time τ ; however Nu and Sh become constant when τ is large.

The effect of increasing Pr , Le , Pm , Ras is found to increase in Nu and Sh thus increasing heat and mass transfer for all three types of modulations.

The effect of increasing δ is to increase the value of Nu and Sh for the case of OPM and LBMO, hence heat and mass transfer.

The effect of increasing ω is to decrease the value of Nu and Sh for the case of OPM and LBMO, hence heat and mass transfer.

In the cases of OPM and LBMO, the natures of Nu and Sh remain oscillatory.

Initially when τ is small, the values of Nusselt and Sherwood numbers start with 1, corresponding to the conduction state. However as τ increases, Nu and Sh also increase, thus increasing the heat and mass transfer.

The values of Nu and Sh for LBMO are greater than those in IPM but smaller than those in OPM.

The effect of magnetic field is to stabilize the system.

Nomenclature Latin Symbols A :

Amplitude of convection

d :

Depth of the fluid layer

g :

Acceleration due to gravity

Q :

Chandrasekhar number, Q = μ m H b 2 d 2 / ρ 0 ν ν m

k c :

Critical wave number

Nu ( τ ) :

Nusselt number

Sh ( τ ) :

Sherwood number

p :

Reduced pressure

Pr :

Prandtl number, Pr = ν / κ T

Pm :

Magnetic Prandtl number, Pm = ν m / κ T

Le :

Lewis number, Le = κ T / κ S

R a T :

Thermal Rayleigh number, R a T = α g Δ T d 3 / ν κ T

Ras :

Solutal Rayleigh number, Ras = β g Δ S d 3 / ν κ T

S :

Solute concentration

R 0 c :

Critical Rayleigh number, R 0 c = ( δ 2 ( δ 4 + Q π 2 ) + Ras k c 2 Le ) / k c 2

T :

Temperature

Δ T :

Temperature difference across the fluid layer

t :

Time

( x , z ) :

Horizontal and vertical space coordinates.

Greek Symbols α :

Coefficient of thermal expansion

β :

Coefficient of solute expansion

δ 2 :

Horizontal wave number δ 2 = k c 2 + π 2

ϵ :

Perturbation parameter

κ T :

Thermal diffusivity

κ S :

Solutal diffusivity

γ :

Heat capacity ratio γ = ( ρ c ) m / ( ρ c ) f

ω :

Frequency of temperature modulation

δ :

Amplitude of temperature modulation

μ :

Dynamic coefficient of viscosity of the fluid

μ m :

Magnetic permeability

ν :

Kinematic viscosity, ν = μ / ρ 0

ν m :

Magnetic viscosity

ρ :

Fluid density

ψ :

Stream function

ψ * :

Dimensionless stream function

Φ :

Magnetic potential

Φ * :

Dimensionless magnetic potential

τ :

Slow time τ = ϵ 2 t

θ :

Phase angle

T :

Perturbed temperature.

Other Symbols 2 :

( 2 / x 2 ) + ( 2 / y 2 ) + ( 2 / z 2 ) .

Subscripts b :

Basic state

c :

Critical

0 :

Reference value.

Superscripts :

Perturbed quantity

* :

Dimensionless quantity.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was done during the lien sanctioned to the author by Banaras Hindu University, Varanasi, to work as professor of Mathematics at Department of Applied Mathematics, School for Physical Sciences, Babasaheb Bhimrao Ambedkar Central University, Lucknow, India. The author B.S. Bhadauria gratefully acknowledges Banaras Hindu University, Varanasi, for the same. Further, the author Palle Kiran gratefully acknowledges the financial assistance from Babasaheb Bhimrao Ambedkar Central University as a research fellowship.

Akbarzadeh A. Manins P. Convective layers generated by side walls in solar ponds Solar Energy 1988 41 6 521 529 10.1016/0038-092X(88)90055-2 2-s2.0-38249030098 Huppert H. E. Sparks R. S. J. Double-diffusive convection due to crystallization in magmas Annual Review of Earth and Planetary Sciences 1984 12 11 37 2-s2.0-0021638860 Fernando H. J. S. Brandt A. Recent advances in double diffusive convection Applied Mechanics Reviews 1994 47 c1 c7 Turner J. S. Buoyancy Effects in Fluids 1979 University Of Cambridge Rudraiah N. Shivakumara I. S. Double-diffusive convection with an imposed magnetic field International Journal of Heat and Mass Transfer 1984 27 10 1825 1836 10.1016/0017-9310(84)90164-9 ZBLl0595.76043 2-s2.0-0021510508 Thomas J. H. Weiss N. O. The theory of sunspots Sunspots: Theory and Observations 1992 Dordrecht, The Netherlands Kluwer Academic Publishers 428 Thompson W. B. Thermal convection in a magnetic field Philosophical Magazine 1951 42 1417 1432 MR0044996 Chandrasekhar S. Hydrodynamic and Hydromagnetic Stability 1961 London, UK Oxford University Press MR0128226 Lortz D. A stability criterion for steady finite amplitude convection with an external magnetic field Journal of Fluid Mechanics 1965 23 113 128 10.1017/S002211206500126X MR0183227 ZBLl0137.44804 Malkus W. V. R. Veronis G. Finite amplitude cellular convection Journal of Fluid Mechanics 1958 4 225 260 10.1017/S0022112058000410 MR0135012 Stommel H. Arons A. B. Blanchard D. An oceanographical curiosity: the perpetual salt fountain Deep Sea Research 1956 3 2 152 153 10.1016/0146-6313(56)90095-8 2-s2.0-49749181781 Gotoh K. Yamada M. Thermal convection of a horizontal layer of magnetic fluids Journal of the Physical Society of Japan 1982 51 9 3042 3048 10.1143/JPSJ.51.3042 2-s2.0-0020175513 Oreper G. M. Szekely J. The effect of an externally imposed magnetic field on buoyancy driven flow in a rectangular cavity Journal of Crystal Growth 1983 64 3 505 515 10.1016/0022-0248(83)90335-4 2-s2.0-0001350499 Rudraiah N. Double-diffusive magnetoconvection Pramana 1986 27 1-2 233 266 10.1007/BF02846340 2-s2.0-51249170765 Siddheshwar P. G. Pranesh S. Magnetoconvection in fluids with suspended particles under 1g and μg Aerospace Science and Technology 2002 6 2 105 114 10.1016/S1270-9638(01)01144-0 ZBLl1006.76545 2-s2.0-0036253593 Bhadauria B. S. Combined effect of temperature modulation and magnetic field on the onset of convection in an electrically conducting-fluid-saturated porous medium Journal of Heat Transfer 2008 130 5 052601 10.1115/1.2885871 2-s2.0-44349132809 Siddheshwar P. G. Bhadauria B. S. Mishra P. Srivastava A. K. Study of heat transport by stationary magneto-convection in a Newtonian liquid under temperature or gravity modulation using Ginzburg-Landau model International Journal of Non-Linear Mechanics 2012 47 5 418 425 10.1016/j.ijnonlinmec.2011.06.006 2-s2.0-84861184165 Drazin P. G. Reid W. H. Hydrodynamic Stability 2004 Cambridge, UK Cambridge University Press 10.1017/CBO9780511616938 MR2098531 Venezian G. Effect of modulation on the onset of thermal convection Journal of Fluid Mechanics 1969 35 2 243 254 10.1017/S0022112069001091 2-s2.0-0014661886 Rosenblat S. Tanaka G. A. Modulation of thermal convection instability Physics of Fluids 1971 14 7 1319 1322 10.1063/1.1693608 ZBLl0224.76038 2-s2.0-0001737835 Roppo M. N. Davis S. H. Rosenblat S. Bénard convection with time-periodic heating The Physics of Fluids 1984 27 4 796 803 10.1063/1.864707 MR741611 2-s2.0-0021406360 Bhadauria B. S. Bhatia P. K. Time-periodic heating of Rayleigh-Bénard convection Physica Scripta 2002 66 1 59 65 10.1238/Physica.Regular.066a00059 2-s2.0-0036644886 Bhadauria B. S. Time-periodic heating of Rayleigh-Benard convection in a vertical magnetic field Physica Scripta 2006 73 3 296 302 10.1088/0031-8949/73/3/010 2-s2.0-33745487906 Malashetty M. S. Swamy M. Effect of thermal modulation on the onset of convection in a rotating fluid layer International Journal of Heat and Mass Transfer 2008 51 11-12 2814 2823 10.1016/j.ijheatmasstransfer.2007.09.031 ZBLl1144.80334 2-s2.0-41949107886 Bhadauria B. S. Kiran P. Heat transport in an anisotropic porous medium saturated with variable viscosity liquid under temperature modulation Transport in Porous Media 2013 100 2 279 295 10.1007/s11242-013-0216-0 MR3121813 2-s2.0-84882319361