Effect of Exponentially Temperature-Dependent Viscosity on the Onset of Penetrative Ferro-Thermal-Convection in a Saturated Porous Layer via Internal Heating ()
1. Introduction
Ferrofluids are stable colloidal suspensions of magnetic nanoparticles in a carrier fluid such as water, hydrocarbon (mineral oil or kerosene), or fluorocarbon. The weirdness of these fluids is the combination of normal liquid behaviour with a magnetic control of their flow and properties. Presently, these fluids are in wide use in seals, bearings, magnetostatic support, jet printers, separation of nonmagnetic particles, flow control and drag reduction, dampers, actuators, sensors, transducers, and medical applications. An authoritative introduction to this fascinating subject along with their applications is provided in [1] [2] [3] .
The magnetization of ferrofluids depends on the magnetic field, temperature, and density. Hence, any variations of these quantities induce change of body force distribution in the fluid and eventually give rise to convection in ferrofluids in the presence of a gradient of magnetic field. There have been numerous studies on thermal convection in a ferrofluid layer called ferroconvection analogous to Rayleigh-Bénard convection in ordinary viscous fluids. The theory of thermal convective instability in a ferrofluid layer began with Finlayson [4] and extensively continued over the years [5] [6] [7] [8] . Nanjundappa and Shivakumara [9] have studied a variety of velocity and temperature boundary conditions on the onset of ferroconvection in an initially quiescent ferrofluid layer. Singh and Bajaj [10] have investigated a time-periodic modulation in temperatures of two horizontal rigid planes containing an initially quiescent ferrofluid layer induces time-periodic oscillations in the fluid layer at the onset of instability.
Thermal convection of ferrofluids saturating a porous medium has also attracted considerable attention in the literature owing to its importance in controlled emplacement of liquids or treatment of chemicals, and emplacement of geophysically imageable liquids into particular zones for subsequent imaging etc. Rosensweig et al. [11] have studied the stabilization of fluid penetration through a porous medium using magnetizable fluids. The stability of the magnetic fluid penetration through a porous medium in high uniform magnetic field oblique to the interface is studied. Zahn and Rosensweig [12] . The thermal convection of a ferrofluid saturating a porous medium in the presence of a vertical magnetic field is studied by Vaidyanathan et al. [13] . The laboratory-scale experimental results of the behavior of ferrofluids in porous media consisting of sands and sediments are presented by Borglin et al. [14] . Sunil et al. [15] have dealt with the theoretical investigation of the double-diffusive convection in a micropolar ferromagnetic fluid layer heated and soluted from below saturating a porous medium. Nanjundappa et al. [16] have explored a model for penetrative ferroconvection in saturated porous layer via internal heat generation.
Majority of ferrofluids are either water-based or oil-based. The viscosity of water is far more sensitive to temperature variations and oils are known to have viscosity decreasing exponentially with temperature rather than linearly. Realizing the importance, several investigators have considered exponential variation in viscosity with temperature in analyzing thermal convective instability in a horizontal fluid layer but the studies are limited to ordinary viscous fluids [17] [18] [19] [20] as well as in a layer of saturated porous medium [21] [22] [23] . To our knowledge, due attention has not been given to investigate convective instability problems involving ferrofluids despite its relevance and importance in many heat transfer applications. Shivakumara et al. [24] have investigated the onset of thermogravitational convection in a horizontal ferrofluid layer with viscosity depending exponentially on temperature.
The intent of the present study is to analyze the influence of viscosity varying exponentially with temperature on the onset of penetrative FTC in a ferrofluid saturated porous layer via internal heating in the presence of a uniform vertical magnetic field. In investigating the problem, the boundaries of the ferrofluid layer is considered to be rigid-ferromagnetic with insulated to temperature perturbations. The resulting eigenvalue problem is solved numerically by the Galerkin technique and analytically by a regular perturbation technique.
2. Mathematical Formulation
The physical configuration considered is horizontal layer of an incompressible ferrofluid of characteristic thickness d in the presence of an imposed spatially uniform magnetic field
in the vertical direction (see Figure 1).
The lower and upper boundaries are maintained at constant but different temperatures
and
, respectively. A Cartesian co-ordinate system (x, y, z) is used with the origin at the bottom and z-axis is directed vertically upward. Gravity acts in the negative z-direction,
, where
is the unit vector in the z-direction. The variation of viscosity
of the ferrofluid with temperature is assumed to be exponential, given by
(1)
where, T is the temperature,
is the reference value at the reference temperature
and
is a positive constant.
The governing equations under the Oberbeck-Boussinesq approximation are given by the following:
Mass balance:
. (2)
Linear momentum balance:
(3)
Energy balance:
(4)
Equation of state:
(5)
Maxwell’s equations in the magnetostatic limit:
,
, (6a,b)
. (6c)
Here,
the seepage velocity vector, p the pressure,
the fluid density,
the magnetic induction,
the magnetization,
the magnetic field intensity,
the magnetic permeability of vacuum, k the permeability of the porous medium,
the average temperature,
the porosity of the porous medium,
the thermal conductivity,
the specific heat at constant volume and magnetic field,
the reference density,
the thermal expansion coefficient and the subscript s represents the solid. In view of Equation (6b),
can be expressed as
(7)
where,
is the magnetic potential.
Since the magnetization depends on the magnitude of magnetic field and temperature, we have
(8)
The linearized equation of magnetic state about
and
is
(9)
where,
is the saturation magnetization,
the magnetic susceptibility,
the pyromagnetic co-efficient,
and
.
It is clear that there exists the following solution for the quiescent basic state:
(10)
where,
is the temperature gradient and the subscript b denotes the basic state. To investigate the conditions under which the quiescent solution is stable against small disturbances, we consider a perturbed state such that
(11)
where,
,
,
,
,
and
are perturbed variables and are assumed to be small. Then, we note that
(12)
Substituting Equation (9) into Equation (8) and using Equation (9) and assuming the
, we get (after dropping the primes)
(13)
(14a)
Hence on using Equation (14), we get
Again substituting Equation (11) into momentum Equation (3), linearizing, eliminating the pressure term by operating curl twice and using Equation (13) the z-component of the resulting equation can be obtained as (after dropping the primes)
(14b)
where,
and
is the horizontal Laplacian operator.
The energy balance Equation (4) after using Equation (11) and linearizing, takes the form (after dropping the primes)
(15)
where
,
and we have assumed
.
Equations 6(a, b), after substituting Equation (11) and using Equation (13), may be written as (after dropping the primes)
. (16)
Since the principle of exchange of stability is valid [5] , the normal mode expansion of the dependent variables is assumed in the form
(17)
where,
and m are wave numbers in the x and y directions, respectively.
On substituting Equation (17) into Equations (14)-(16) and non-dimesionalizing the variables by setting
(18)
where, kinematic viscosity,
, effective thermal diffusivity
and
.
Following the classical lines of linear stability theory as presented by Chandrasekhar [25] , neglecting the asterisk, the linearized and dimensionless governing equation can be written as:
(19)
(20)
. (21)
Here,
is the differential operator,
is the overall horizontal wavenumber, W is the amplitude of vertical component of velocity,
is the amplitude of temperature,
is the amplitude of magnetic potential,
is the thermal Raleigh number,
is the magnetic Rayleigh number,
is the magnetic number,
is the magnetic parameter,
is the measure of nonlinearity of fluid magnetization parameter. The typical value of
for magnetic fluids with different carrier liquids turns out to be of the order of 10−6 and hence its effect is neglected as compared to unity and
is given by
(22)
where,
is the dimensionless viscosity parameter. If the reference temperature
, then
. (23)
The boundaries are considered to be rigid-ferromagnetic and they are insulated to temperature perturbations:
at
(24)
3. Method of Solution
Equations (19)-(21) together with the corresponding boundary condition (24) constitute an eigenvalue problem with
or
as an eigenvalue. The eigenvalue problem is solved numerically using the Galerkin technique as well as analytically using a regular perturbation technique and the results so obtained are compared to know the accuracy of the methods employed.
3.1. Solution by the Galerkin Technique
The Galerkin method is used to solve this problem as explained in the book by Finlayson [26] . In this method, the test (weighted) functions are the same as the base (trial) functions. Accordingly, W,
and
are written as
(25)
where,
,
and
are unknown constants to be determined. The base functions
,
and
are generally chosen such that they satisfy the boundary condition (24) but not the differential equations. For this boundary, they are chosen as
(26)
where,
’s are the modified Chebyshev polynomials. The above trial functions satisfy the boundary condition. Multiplying Equation (19) by
, Equation (20) by
and Equation (21) by
; performing the integration by parts with respect to z between
and
and using the boundary conditions, we obtain a system of linear homogeneous algebraic equations in
,
and
. Nontrivial solution exists if and only if the characteristic determinant is equal to zero. This leads to a relation involving the characteristic equation in the form
.
The critical values of
or
are found as a function of wave number a for various values of physical parameters. It is observed that the convergence is achieved with six terms in the series expansion.
3.2. Solution by Regular Perturbation Technique
It is known that for insulated boundary conditions the onset of convection corresponds to a vanishingly small wave number (i.e. unicellular convection). The numerical calculations carried out in the previous section also corroborate this fact. Therefore, an attempt is being made to exploit this fact to obtain an analytical formula for the onset of convection using a regular perturbation technique with wave number a as a perturbation parameter. Accordingly, the variables W,
and
are expanded in powers of a2 as
(27)
Substituting Equation (27) into Equations (19)-(21) and also in the boundary conditions, and collecting the terms of zero-th order, we obtain
(28a)
(28b)
(28c)
with the boundary condition
. The solution to the zero-th order equations is:
, and
.
The first order equations are then
(29a)
(29b)
(29c)
with the boundary conditions
.
The general solution of Equation (29a) is given by
(30)
where,
,
,
,
,
,
,
with
,
,
,
.
From Equation (29b), it follows that
. (31)
Substituting for
from Equation (30) into Equation (31) and carrying out the integration leads to an expression for the critical Rayleigh number
in the form
where
with
,
,
,
,
,
,
.
As
and
, Equation (32), reduce to
(33)
These are the results for constant viscosity ferrofluids and coincide with Nanjundappa and Shivakumara [16] . When
(i.e. ordinary viscous fluid), Equation (33) reduce to the critical Rayleigh number of
, which is he known exact value documented in the literature. Equation (33) further reveal that the nonlinearity of fluid magnetization (i.e.
) has no effect on the onset of ferroconvection; a result which is observed by numerical computations carried out in the previous section. This result is similar to the one noticed in the case of constant viscosity ferrofluids [16] . Since at the onset of convection
(very large wave length), one would expect that
has no effect on the stability of the system.
4. Results and Discussion
The linear stability analysis is carried out with viscosity depending exponentially on temperature on the onset of FTC in a ferrofluid saturated porous layer. The bounding surface of the ferrofluid layer is rigid-ferromagnetic and insulated to temperature perturbations. The critical eigenvalue
or
and the corresponding wave number
are computed numerically by the Galerkin method as well as analytically by employing a regular perturbation technique for different
and B. It is noted that the critical wave number is vanishingly small and this fact is exploited to obtain an analytic expression for
using a regular perturbation technique with wave number a as a perturbation parameter. Such a study also helps in knowing the accuracy of the numerical method employed in solving the problem. The stability characteristics of the system are found to be independent of the nonlinearity of fluid magnetization parameter
. The salient features of the physical parameters on the onset of FTC are exhibited in Figures 2-7.
Figure 2 shows the variation of
as a function of viscosity parameter B for different magnetic parameter
when
. The figure clearly illustrates the strong influence of viscosity parameter B on FTC.
As a result of viscosity variation, two distinguish regions are shown,
increases negligibly small for B up to 5.17019, at which maximum value of
are reached; rapidly decreasing trends are found for
(i.e. absence of magnetic force) and decreasing slowly for
(i.e. presence of magnetic force). At maximum value of
a sublayer starts to form. It is seen that maximum
exists at
and
decreases monotonically with increasing B indicating
Figure 2. Variation of
versus B for different
when
.
Figure 3. Variation of
versus B for different
when
.
Figure 4. Locus of
versus
for different
for
.
Figure 5. Vertical velocity eigenfunction for
and
.
Figure 6. Vertical velocity eigenfunction for
and
.
Figure 7. Vertical velocity eigenfunction for
and
.
its effect is to destabilize the system. This is due to the decrease in viscosity of the ferrofluid with temperature. Moreover,
decreases quite rapidly at first then slowly and finally the curves of different
merge with increasing B. It is more so with an increase in the value of
and this is due to additive reinforcement of destabilizing magnetic force. The results for
correspond to ordinary viscous fluid and it is observed that higher heating is required to have instability in this case. Thus the combined effect of temperature dependent viscosity and magnetic forces is to reinforce together and to hasten the onset of FTC compared to their effect in isolation. The value of
at which
attains its maximum value
are tabulated in Table 1 for different values of
when
. From Table 1 it is seen that
decreases with increasing
.
Figure 3 shows variation of
with variation parameter B for various values of
when
. It is seen that
increases with increasing
and hence its effect is to delay the onset of FTC.
For a fixed thickness of the porous layer, increase in the value of
leads to decrease in the permeability of the porous medium which in turn retards the flow of ferrofluid. Therefore, higher heating and hence higher value of
is required for the onset of
in a ferrofluid saturated porous medium. For different
, Figure 3 demonstrates two distinct characteristics and which is same situation in the presence of magnetic forces with increase in B.
The value of B at which
attains its maximum value
are tabulated in Table 2 for different
when
. From Table 2 it is seen that
and
decreases with increasing B.
The locus of
and
is shown in Figure 4 for different
with
to know the simultaneous presence of buoyancy and magnetic forces on the stability of the system.
From the figures it is obvious that there is a strong coupling between
and
and the curves are slightly convex. That is, when the buoyancy force is predominant the magnetic force becomes negligible and vice-versa. From Figure 4, it is seen that an increase in
is to increase in
as well as
and thus its effect is to delay the onset of ferroconvection.
The perturbed vertical velocity eigenfunctions are presented in Figures 5-7 for different values of
, B and
, respectively.
Table 1. Values of
occurring at B for different values of
when
.
Table 2. Values of
occurring at B for different values of
when
.
As can be seen, the shape of the eigenfunction is parabolic in nature. Increasing
, decreasing B and
is to increase the vigor of the ferrofluid flow and hence their effect is to hasten the onset of ferroconvection.
5. Conclusion
The onset of ferroconvection in a ferofluid saturated porous layer with viscosity varying exponentially with temperature is studied. The viscosity parameter B exhibits a dual effect on the stability characteristics of the system. It shows stabilizing effect on the system viscosity. The viscosity parameter B exhibits a dual effect on the stability characteristics of the system. It shows stabilizing effect on the system initially but displays a reverse trend after exceeding certain value of B. The critical Rayleigh number
attains maximum value at some intermediate values of B. The effect of the increase in
and the Darcy number Da is to hasten the onset of FTC. The buoyancy and magnetic forces reinforce each other in hastening the onset of FTC. The nonlinearity of fluid magnetization parameter
has no effect on the onset of FTC. The critical eigenvalues were obtained by a regular perturbation technique and computed numerically using the Galerkin method complement with each other indicating the analytical solutions obtained are exact.
Acknowledgements
The authors acknowledge the financial support received in the form of a “Research Fund for Talented Teacher Scheme” from Vision Group of Science & Technology, Government of Karnataka, Bengaluru, India.