International Journal of Engineering Mathematics

Volume 2014 (2014), Article ID 296216, 14 pages

http://dx.doi.org/10.1155/2014/296216

Research Article

Weak Nonlinear Double-Diffusive Magnetoconvection in a Newtonian Liquid under Temperature Modulation

B. S. Bhadauria and Palle Kiran

Department of Applied Mathematics, School for Physical Sciences, Babasaheb Bhimrao Ambedkar University, Lucknow 226 025, India

Received 22 February 2014; Accepted 15 June 2014; Published 6 July 2014

Academic Editor: Yurong Liu

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.

Abstract.

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 [1], electrochemistry, magma chambers and sparks [2], Fernando and Brandt[3], 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 [6].

Thompson [7] and Chandrasekhar [8] were the first to study the magnetoconvection in horizontal fluid layer. Lortz [9] 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 [10] but his analysis is silent about the detailed study of stability analysis. Stommel et al. [11] 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 [12] 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 [13] 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 [5] 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 [14] investigated the interaction between double-diffusive convection and an externally imposed vertical magnetic field in a Boussinesq fluid. Siddheshwar and Pranesh [15] 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 [16] also studied the effect of magnetic field on thermal modulated convection in the case of porous medium. Siddheshwar et al. [17] 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 [8] and Drazin and Reid [18]. 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 [19], who was the first to consider the effect of temperature modulation on thermal instability in a horizontal fluid layer. Rosenblat and Tanaka [20] 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. [21]. Bhadauria and Bhatia [22] studied the effect of temperature modulation on thermal instability by considering rigid-rigid boundaries and different types of temperature profiles. Bhadauria [23], Malashetty and Swamy [24], and Bhadauria and Kiran [25] 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 , confined between two infinite parallel horizontal planes at and . A Cartesian frame of reference has been taken with the origin at the bottom, and the axis vertically upward direction of the fluid layer. The surfaces are maintained at a constant temperature gradient and a constant magnetic field 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

where is velocity , is a viscosity, is magnetic permeability, is magnetic viscosity, is the thermal diffusivity, is the solutal diffusivity, 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, is gravitational acceleration, while is reference density, and is the magnetic field. The externally imposed thermal boundary conditions considered in this paper are given by Venezian [19]:

where is small amplitude of modulation, 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:

The basic state is assumed to be quiescent and the quantities in this state are given by

The basic state pressure field is not required here; however the basic temperature and solute fields are governed by the following ordinary differential equations:

The solutions of (11) subjected to the boundary conditions (8)-(9) are given by

where is the study temperature field and is the oscillating part while stands for the real part. We assume finite amplitude perturbations on the basic state in the form

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:

Further we consider only two-dimensional disturbances in our study and hence the stream function and magnetic potential are introduced as and . By eliminating density and pressure terms from (14) and nondimensionalizing using the following transformations: , , , , ,  , , and , finally dropping the asterisk for simplicity, we get the nondimensionalized governing system of equations:

Here the nondimensionalizing parameters in the above equations are magnetic Prandtl number, Prandtl number, is Lewis number, is thermal Rayleigh number, is Solutal Rayleigh number, and is Chandrasekhar number. In (16), one can observe that the basic state solution influences the stability problem through the factor , which is given by

where

, , and . We assume small variations of time and rescaling it as , to study the stationary convection of the system, we write the nonlinear Equations (15)–(18) in the matrix form as given below:

The considered stress free and isothermal boundary conditions to solve the above system (21) are

where .

Figure 1: Physical configuration of the problem.

3. Finite Amplitude Equation, Heat and Mass Transport

We now introduce the following asymptotic expansions in (21):

where 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

The solutions of the lowest order of the system subject to the boundary conditions equation (22) are

where . The critical value of the Rayleigh number for the onset of magnetoconvection in the absence of temperature modulation is

which is the same as Siddheshwar et al. [17] when and we obtain classical results of Chandrasekhar [8] for without magnetic field and single component fluid layer.

At the second order, we have

The second order solutions subjected to the boundary conditions Equation (22) are obtained as follows:

The horizontally averaged Nusselt and Sherwood numbers for the stationary double-diffusive magnetoconvection (the mode considered in this problem) are given by

One must note here that is effective at and affects and through can be seen later. Therefore, substituting , , , respectively, into (29) and simplifying we obtain

At the third order, we have

where

Substituting , , and into (32), we obtain the expressions for , , , and 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

where

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 , where is the chosen initial amplitude of convection. In our calculations we may use , 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 in which the salt and heat make opposing contributions . We also consider the effect of temperature modulation to be of order ; 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:

(1)the need for nonlinear stability analysis,(2)the relation of the problem to a real application,(3)the selection of all dimensionless parameters utilized in computations.

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

(1)in-phase modulation [IPM] ,(2)out-of-phase modulation [OPM] ,(3)lower boundary modulation [LBMO] .

The parameters of the system are , , , , , , , ; 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 , . There is the property of the fluid coming into picture as well as through Prandtl number . 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 , , , , , , , 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 2–7 show the individual effect of each nondimensional parameter on heat and mass transfer.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

Figure 2: versus (a) (b) (c) (d) (e) (f) , .

(a)

(b)

(c)

(d)

(e)

(f)

(g)

Figure 3: versus (a) (b) (c) (d) (e) (f) (g) .

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 4: versus (a) (b) (c) (d) (e) (f) (g) .

(a)

(b)

(c)

(d)

(e)

(f)

(g)

Figure 5: versus (a) (b) (c) (d) (e) (f) , .

(a)

(b)

(c)

(d)

(e)

(f)

(g)

Figure 6: versus (a) (b) (c) (d) (e) (f) (g) .

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 7: versus (a) (b) (c) (d) (e) (f) (g) .

(1) Figures 2(a)–7(a) show that the effect of Chandrasekhar number which is ratio of Lorentz force to viscous force is to delay the onset of convection, hence heat and mass transfer. The and 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 shows oscillatory behaviour and increment in it decreases the magnitude of both and . Hence 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 .

(2) The effect of Prandtl number 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 and Magnetic Prandtl numbers is to advance the convection and hence heat and mass transfer. Hence both and 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. [17].

(4) The effect of solutal Rayleigh number 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 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 and 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 , 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 [19], Siddheshwar et al. [17], and Bhadauria and Kiran [25].

(6) From Figures 2(g) and 5(g) we observe that has strongly stabilizing effect. .

(7) The comparison of three types of temperature modulations is given in Figures 4(h) and 7(h): .

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

(1), Figures 2(a)–7(a).(2), Figures 2(b)–7(b).(3), Figures 2(c), 3(c), and 4(c).(4), Figures 5(c), 6(c), and 7(c).(5), Figures 2(d)–7(d).(6), Figures 2(e)–7(e).(7), Figures 3(f), 4(f), 6(f), and 7(f).(8), 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.

(1)The effect of IPM is negligible on heat and mass transport in the system.(2)In the case of IPM, the effect of and is also found to be negligible on heat and mass transport.(3)In the case of IPM, the values of and increase steadily for small values of time ; however Nu and Sh become constant when is large.(4)The effect of increasing , , , is found to increase in and thus increasing heat and mass transfer for all three types of modulations.(5)The effect of increasing is to increase the value of and for the case of OPM and LBMO, hence heat and mass transfer.(6)The effect of increasing is to decrease the value of and for the case of OPM and LBMO, hence heat and mass transfer.(7)In the cases of OPM and LBMO, the natures of and remain oscillatory.(8)Initially when is small, the values of Nusselt and Sherwood numbers start with 1, corresponding to the conduction state. However as increases, and also increase, thus increasing the heat and mass transfer.(9)The values of and for LBMO are greater than those in IPM but smaller than those in OPM.(10)The effect of magnetic field is to stabilize the system.

Nomenclature

Latin Symbols

: | Amplitude of convection
: | Depth of the fluid layer
: | Acceleration due to gravity
: | Chandrasekhar number,
: | Critical wave number
: | Nusselt number
: | Sherwood number
: | Reduced pressure
: | Prandtl number,
: | Magnetic Prandtl number,
: | Lewis number,
: | Thermal Rayleigh number,
: | Solutal Rayleigh number,
: | Solute concentration
: | Critical Rayleigh number,
: | Temperature
: | Temperature difference across the fluid layer
: | Time
: | Horizontal and vertical space coordinates.

Greek Symbols

: | Coefficient of thermal expansion
: | Coefficient of solute expansion
: | Horizontal wave number
: | Perturbation parameter
: | Thermal diffusivity
: | Solutal diffusivity
: | Heat capacity ratio
: | Frequency of temperature modulation
: | Amplitude of temperature modulation
: | Dynamic coefficient of viscosity of the fluid
: | Magnetic permeability
: | Kinematic viscosity,
: | Magnetic viscosity
: | Fluid density
: | Stream function
: | Dimensionless stream function
: | Magnetic potential
: | Dimensionless magnetic potential
: | Slow time
: | Phase angle
: | Perturbed temperature.

Other Symbols

:
.

Subscripts

: | Basic state
: | Critical
: | 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.

References

A. Akbarzadeh and P. Manins, “Convective layers generated by side walls in solar ponds,” Solar Energy, vol. 41, no. 6, pp. 521–529, 1988.

H. E. Huppert and R. S. J. Sparks, “Double-diffusive convection due to crystallization in magmas,” Annual Review of Earth and Planetary Sciences, vol. 12, pp. 11–37, 1984.

H. J. S. Fernando and A. Brandt, “Recent advances in double diffusive convection,” Applied Mechanics Reviews, vol. 47, pp. c1–c7, 1994.

J. S. Turner, Buoyancy Effects in Fluids, University Of Cambridge, 1979.

N. Rudraiah and I. S. Shivakumara, “Double-diffusive convection with an imposed magnetic field,” International Journal of Heat and Mass Transfer, vol. 27, no. 10, pp. 1825–1836, 1984.

J. H. Thomas and N. O. Weiss, “The theory of sunspots,” in Sunspots: Theory and Observations, p. 428, Kluwer Academic Publishers, Dordrecht, The Netherlands, 1992.

W. B. Thompson, “Thermal convection in a magnetic field,” Philosophical Magazine, vol. 42, pp. 1417–1432, 1951.

S. Chandrasekhar, Hydrodynamic and Hydromagnetic Stability, Oxford University Press, London, UK, 1961.

D. Lortz, “A stability criterion for steady finite amplitude convection with an external magnetic field,” Journal of Fluid Mechanics, vol. 23, pp. 113–128, 1965.

W. V. R. Malkus and G. Veronis, “Finite amplitude cellular convection,” Journal of Fluid Mechanics, vol. 4, pp. 225–260, 1958.

H. Stommel, A. B. Arons, and D. Blanchard, “An oceanographical curiosity: the perpetual salt fountain,” Deep Sea Research, vol. 3, no. 2, pp. 152–153, 1956.

K. Gotoh and M. Yamada, “Thermal convection of a horizontal layer of magnetic fluids,” Journal of the Physical Society of Japan, vol. 51, no. 9, pp. 3042–3048, 1982.

G. M. Oreper and J. Szekely, “The effect of an externally imposed magnetic field on buoyancy driven flow in a rectangular cavity,” Journal of Crystal Growth, vol. 64, no. 3, pp. 505–515, 1983.

N. Rudraiah, “Double-diffusive magnetoconvection,” Pramana, vol. 27, no. 1-2, pp. 233–266, 1986.

P. G. Siddheshwar and S. Pranesh, “Magnetoconvection in fluids with suspended particles under 1g and μg,” Aerospace Science and Technology, vol. 6, no. 2, pp. 105–114, 2002.

B. S. Bhadauria, “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, vol. 130, no. 5, Article ID 052601, 2008.

P. G. Siddheshwar, B. S. Bhadauria, P. Mishra, and A. K. Srivastava, “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, vol. 47, no. 5, pp. 418–425, 2012.

P. G. Drazin and W. H. Reid, Hydrodynamic Stability, Cambridge University Press, Cambridge, UK, 2004.

G. Venezian, “Effect of modulation on the onset of thermal convection,” Journal of Fluid Mechanics, vol. 35, no. 2, pp. 243–254, 1969.

S. Rosenblat and G. A. Tanaka, “Modulation of thermal convection instability,” Physics of Fluids, vol. 14, no. 7, pp. 1319–1322, 1971.

M. N. Roppo, S. H. Davis, and S. Rosenblat, “Bénard convection with time-periodic heating,” The Physics of Fluids, vol. 27, no. 4, pp. 796–803, 1984.

B. S. Bhadauria and P. K. Bhatia, “Time-periodic heating of Rayleigh-Bénard convection,” Physica Scripta, vol. 66, no. 1, pp. 59–65, 2002.

B. S. Bhadauria, “Time-periodic heating of Rayleigh-Benard convection in a vertical magnetic field,” Physica Scripta, vol. 73, no. 3, pp. 296–302, 2006.

M. S. Malashetty and M. Swamy, “Effect of thermal modulation on the onset of convection in a rotating fluid layer,” International Journal of Heat and Mass Transfer, vol. 51, no. 11-12, pp. 2814–2823, 2008.

B. S. Bhadauria and P. Kiran, “Heat transport in an anisotropic porous medium saturated with variable viscosity liquid under temperature modulation,” Transport in Porous Media, vol. 100, no. 2, pp. 279–295, 2013.