Effect of Magnetic Field on Kelvin-Helmholtz Instability in a Couple-Stress Fluid Layer Bounded Above by a Porous Layer and Below by a Rigid Surface

Krishna B. Chavaraddi^{1}^{*},
Vishwanath B. Awati^{2},
Nagaraj N. Katagi^{3},
Priya M. Gouder^{2,4}

Show more

1. Introduction

Kelvin-Helmholtz instability is one of the basic instabilities of two-fluid systems, which affects an interface. The prototypical case is that with one layer of lighter fluid overlying another of denser fluid, and the two moving horizontally in the same direction but with different velocities.

It is not uncommon for environmental fluids to be subject simultaneously to the destabilizing effect of a velocity shear and the stabilizing effect of density stratification, and, when such competition occurs, the outcome is often the so-called Kelvin-Helmholtz (KH) instability [1] . Ever since von Helmholtz [2] and Kelvin [3] developed the theory, this instability has become a standard staple of fluid mechanics, and the basic theory can be found in numerous textbooks, for example Lamb ( [4] , pp. 373-374), Turner ( [5] , pp. 93-96), Kundu ( [6] , pp. 373-381) and Scorer ( [7] , pp. 231-234) to cite a few up to the present time. Investigations into the details of the instability, including secondary instabilities [8] have been carried out perhaps with far more depth than for any other type of fluid instability.

The importance of the Kelvin-Helmholtz (KH) instability of parallel flows in laboratory, geophysical, or astrophysical systems, recognized many years ago, has generated a huge literature. The study of the Kelvin-Helmholtz instability has a long history in hydrodynamics, The basic linear stability analysis of the magnetohydrodynamvic (MHD) K-H instability was carried out long ago (Chandrasekhar [8] ). There is now also a growing literature of the nonlinear evolution of the MHD K-H instability beginning from a variety of possible initial Ñow configurations, at least in the earlier evolution stages in two dimensions. Strong magnetic fields, through their tension, are well known to stabilize the K-H instability. However, the considerable potential for much weaker fields to modify the nonlinear instability and, in particular, to reorganize the subsequent flow has only recently been emphasized.

Malik and Singh [9] investigated the nonlinear Kelvin-Helmholtz properties of (2 + 1) dimensional wave packets propagating at the interface of two superposed Ferrofluids. They considered that the fluids as moving with uniform speeds parallel to the common interface and subjected to a tangential magnetic field. They derived a nonlinear equation which governs the evolution of the amplitude of the system. The effect of a time- dependent acceleration in the presence of a tangential magnetic field on the nonlinear Kelvin-Helmholtz instability has been discussed by El-Dib [10] . El-Sayed [11] investigated the RTI problem of rotating stratified conducting fluid layer through porous medium in the presence of an inhomogenous magnetic field. This problem corresponds physically (in astrophysics) to the RTI of an equatorial section of a planetary magnetosphere or of stellar atmosphere when rotation and magnetic field are perpendicular to gravity. The KHI of two superposed viscous fluids in a uniform vertical magnetic field is discussed in the presence of effects of surface tension and permeability of porous medium by Bhatia and Sharma [12] . Following Babchin et al., [13] and Rudraiah et al., [14] , a simple theory based on Stokes and lubrication approximations is used in this study with the primary objective of using porous layer to suppress the growth rate of KHI.

In the above studies the fluid has been considered to be Newtonian. In the recent years a great deal of interest has been focused on the understanding of the couple stress effects occurring in the flow of non-Newtonian fluids through porous media. This problem appears to be, at this time, of special interest in oil reservoir engineering, where an increasing interest is being shown in the possibility of improving oil recovery efficiency from water flooding projects through mobility control with non-Newtonian displacing fluids. Consequently, it has become essential to have an adequate understanding of the couple stress effect of non-Newtonian displacing and displaced fluids in an oil displacement mechanism. Many technological processes involve the parallel flow of fluids of different viscosity, elasticity and density through porous media. Such flows exist in packed bed reactors in the chemical industry, petroleum engineering, boiling in porous media and in many other processes. Should the interface between the two fluids become unstable, a substantial increase in the resistance to the flow will result. This increase in resistance, in turn, may cause flooding in counter current packed chemical reactors and dry out in boiling porous media. In the same vein, in petroleum production engineering, such instabilities lead to emulsion formation. Hence, the knowledge of the conditions for the onset of instability will enable us to predict the limiting operation conditions of the above processes.

El-Dib and Matoog [15] have studied the Electrorheological Kelvin-Helmholtz instability of a fluid sheet. This work deals with the gravitational stability of an electrified Maxwellian fluid sheet shearing under the influence of a vertical periodic electric field. The field produces surface charges on the interfaces of the fluid sheet. Due to the rather complicated nature of the problem a mathematical simplification is considered where the weak effects of viscoelastic fluids are taken into account. The effect of boundary roughness on Kelvin-Helmholtz instability in Couple stress fluid layer bounded above by a porous layer and below by rigid surface is studied by Chavaraddi et al., [16] . Recently, they [17] have observed the effect of surface roughness on Kelvin-Helmholtz instability in presence of magnetic field. The objective of this paper is to study the effect of magnetic field on Kelvin-Helmholtz discontinuity between two couple-stress viscous conducting fluids in a transverse magnetic field through a porous medium in the presence of the effects of surface tension at the interface.

The paper is organized as follows. The basic equations are established in Section 2 together with Maxwell’s equations. The basic equations are simplified and non-dimen- sionalized using the following Stokes and lubrication approximations in this section. The resulting dispersion relation is derived using suitable boundary and surface conditions in Section 3. The cutoff and maximum wave numbers and the corresponding maximum growth rate are also obtained in Section 3. The results are discussed in Section 4 and some important conclusions are drawn in final section of this paper.

2. Mathematical Formulation

The physical configuration is shown in Figure 1. We consider a thin target shell in the form of a thin film of unperturbed thickness h (Region 1) filled with an incompressible, viscous, poorly electrically conducting light fluid of density bounded below by a rigid surface at y = 0 and above by an incompressible, viscous poorly conducting heavy fluid of density saturating a dense porous layer of large extent compared to the shell thickness h. The co-ordinates x and y spans the horizontal and vertical directions. The interfacial is denoted by. When the interface is flat then when. The fluid velocity vector and the fluid is assumed to be non-Newtonian (couple-stress fluid), viscous electrically conducting and incompressible. The viscosity of fluid (porous medium) is given by, the porous parameter, the permeability of the porous medium and is the slip parameter at the interface. The stress gradient is related to the gravitational acceleration through the relation. The perturbed interface is along the y direction.

The basic equations for clear fluid layer (region 1) and those for porous layer (region 2) are as given below:

Region-1:

(2.2)

Maxwell’s Equations:

(2.3)

and the auxiliary equations

Figure 1. Physical configuration.

(2.4)

Region-2:

(2.5)

where the fluid velocity, the electric field, the magnetic field, the current density, the dielectric field, the magnetic induction, the electrical conductivity, k the permeability of the porous medium, p the pressure, magnetic permeability, the uniform Darcy velocity, the fluid viscosity, the couple-stress parameter and the fluid density.

The basic equations are simplified using the following Stokes and lubrication and electrohydrodynamic approximations (See Rudraiah et al. [14] ):

1) The electrical conductivity of the liquid, s, is negligibly small, i.e.,.

2) The film thickness h is much smaller than the thickness H of the dense fluid above the film. That is

3) The surface elevation h is assumed to be small compared to film thickness h. That is

4) The Strauhal number S, a measure of the local acceleration to inertial acceleration in Equation (2.2), is negligibly small.

That is

where is the characteristic velocity, the kinematic viscosity, the characteristic length and the characteristic time.

Under these approximations Equations (2.1) and (2.2) for fluid in the film, after making dimensionless using

(2.6)

become (after neglecting the asterisks for simplicity).

Region 1:

(2.7)

(2.8)

(2.9)

where is the couple-stress parameter and the Hart-

mann number.

Region 2:

(2.10)

where is the porous parameter.

3. Dispersion Relation

To find the dispersion relation, first we have to find the velocity distribution from Equation (2.8) using the following boundary and surface conditions:

(3.1)

(3.2)

where

(3.3)

(3.4)

Here is the Bond number and is the elevation of the interface.

The solution of (2.8) subject to the above conditions is

(3.5)

where

After integrating Equation (2.7) with respect to y between y = 0 and 1 and using Equation (3.5), we get

(3.6)

where

.

Then Equation (3.3), using Equations (3.6) and (3.4), becomes

. (3.7)

To investigate the growth rate, n, of the periodic perturbation of the interface, we look for the solution of Equation (3.7) in the form

(3.8)

where is the wave number and is the amplitude of perturbation of the interface.

Substituting Equation (3.8) into (3.7), we obtain the dispersion relation in the form

. (3.9)

where.

Also, Equation (3.9) can be expressed as

(3.10)

where

, ,.

Setting n = 0 in Equation (3.9), we obtain the cut-off wavenumber, in the form

(3.11)

because and are non-zero.

The maximum wavenumber, obtained from Equation (3.9)) by setting is

(3.12)

because and are different from zero.

The corresponding maximum growth rate, n_{m}, is

(3.13)

Similarly, using, we obtain

(3.14)

and hence

. (3.15)

The growth rate given by Equation (3.9) is computed numerically for different values of parameters and the results are presented graphically in Figures 2-5.

Figure 2. Growth rate, n versus the wavenumber, for different values of Hartmann number M when and.

Figure 3. Growth rate, n versus the wavenumber, for different values of Couple-stress parameter when and.

Figure 4. Growth rate, n versus the wavenumber, for different values of Bond number B, when and.

Figure 5. Growth rate, n versus the wavenumber, for different values of porous parameter when and.

4. Results and Discussion

In this study we have shown the effect of physical parameters involved in the problem on effect of magnetic field on surface instability of KH type in a couple- stress fluid layers bounded above by a porous layer and below by a rigid boundary. Numerical calculations were performed to determine the growth rate at different wavenumbers for various fluid properties like couple stress parameter M_{0}, Hartmann number M, Bond number B and porous parameter. We have plotted the dimensionless growth rate of the perturbation against the dimensionless wavenumber for some of the cases only. In the linear stage, all perturbed values grow exponentially in agreement with the dispersion relation Equation (3.9). At this stage the interface between the layers acquires a sinusoidal shape of small amplitude.

We have investigated the role of the magnetic field on the two-layer channel flow problem, demonstrated that either destabilization or stabilization can be obtained and presented growth rates in situations where the magnetic field is stabilizing over a broad range of wavenumbers for increasing in Hartmann number M in Figure 3 where and. The increasing the Hartmann ratio results in slightly increasing the critical wavenumber and decreasing the maximum growth rate. It thus has a stabilizing effect for the selected values of input parameters due to the increased in Hartmann ratio (Lorentz force to viscous force).

Also, when fix all the input parameters we find that the higher the couple-stress parameter the more stable the interface is. In Figure 2, we have plotted the growth rate against the wavenumber in the case where and for different values of the couple-stress parameter M_{0}. Increasing the couple-stress ratio results in slightly increasing the critical wavenumber and decreasing the maximum growth rate this is because of the action of the body couples on the system. Thus, it has a stabilizing effect for the selected values of input parameters due to the increased in the couple-stress parameter.

In addition, we have investigated the effect of the surface tension of the fluid on the instability of the interface. In our sample calculations, we have taken, , and and varied the Bond number B. For this input parameters, the critical wavenumber and maximum growth rate decreased as the ratio of the Bond number B decreased from 0.4 to 0.1 as observed in Figure 4. The Bond number is reciprocal of surface tension and thus showing that an increase in surface tension decreases the growth rate and hence make the interface more stable.

However, in order to understand the effect of the porous properties on the instability, we now fix values of other parameters, , and and vary the ratios of the porous parameters. Figure 5 displays the results of our calculations, showing that increasing the ratio of porous parameters s_{p} from 20 to 100 (and thus increasing the Darcy resistance compared to the viscous force) increases the critical wavelength and decreases the maximum growth rate, thus having a stabilizing effect by this parameter. We conclude that an increase in also stabilizes the KHI due to the resistance offered by the solid particles of the porous layer to the fluid.

5. Conclusions

Kelvin-Helmholtz instability is one of the basic mechanisms, which influence the two-phase flow. When relative velocity is larger than the critical velocity, the instability occurs. With instability, flow regime is changed, and also interface surface is significant enlarged, which influence the heat and mass transfer. With linearised Navier-Stokes equations we can analytically predict onset of instability and wavelength for inviscid flow.

We have studied the linear stability of a two-fluid flow in a channel where the fluids are assumed to be Newtonian with different fluid properties (Hartmann number, couple-stress ratio, Surface tension and porous parameter) and subjected to magnetic field normal to their interface. For this purpose, we have derived and then linearized the equations of motion where the interaction between the hydrodynamic and couple- stress problems occurs through the stress balance at the fluid interface. The growth rate of the perturbation was then computed by using the normal mode method and its variation studied as a function of the dimensionless parameter Hartmann number M, couple-stress parameter M_{0}, as well as Bond number B and porous parameter. While two layer flows in channels of small dimensions are rather stable, the instability of the fluid-porous interface is highly desirable in certain cases, particularly for chemical industry, in petroleum production engineering applications where the mixing of reagents are crucial steps in the process. However, in systems of larger scale, the instability of the fluid-porous interface in a channel is often an undesired physical phenomenon. In such situations, controlling the flow requires the stabilization of the interface. In searching for a method capable of either stabilizing a potentially unstable interface or destabilizing a potentially stable one, we have investigated the role of the magnetic field on the two-layer channel flow problem, demonstrated that either destabilization or stabilization can be obtained and presented growth rates in situations where the magnetic field is stabilizing over a broad range of wavenumbers for increasing in Hartmann number M as same behavior observed by varying the couple-stress parameter. But in the case of variation in Bond number is to increase in surface tension decreases the growth rate and hence make the interface more stable. Also we conclude that the increase in the porous parameter is to decrease the growth rate showing thereby the stabilizing effect on the interface.

Acknowledgements

This work is supported by the VGST, Department of Science & Technology, Government of Karnataka under grant no. VGST/SMYSR(2014-15)/GRD-433/2015-16, Dated : 28.04.2015 and the authors (NNK, VBA and PMG) wishes to thank respectively the Director/VC/Principal of their institutions for their encouragement and support in doing research.

References

[1] De Silva, I.P.D., Fernando, H.J.S., Eaton, F. and Hebert, D. (1996) Evolution of Kelvin-Helmholtz Billows in Nature and Laboratory. Earth and Planetary Science Letters, 143, 217-231.

http://dx.doi.org/10.1016/0012-821X(96)00129-X

[2] Von Helmholtz, H.L.F. (1868) Uber discontinuierliche Flüssigkeitsbewegungen. Monatsberichte der k?nigl. Akad. der Wissenschaften zu Berlin, 215-228.

[3] Kelvin, L. (1871) Hydrokinetic Solutions and Observations. Philadelphia Magazine, 10, 155-168.

[4] Lamb, H. (1945) Hydrodynamics. 6th Edition, Dover Publications, New York, 738 p.

[5] Turner, J.S. (1973) Buoyancy Effects in Fluids. Cambridge University Press, Cambridge, 368 p.

http://dx.doi.org/10.1017/CBO9780511608827

[6] Kundu, P.K. (1990) Fluid Mechanics. Academic Press, Cambridge, 638 p.

[7] Scorer, R.S. (1997) Dynamics of Meteorology and Climate. John Wiley & Sons, Hoboken, 686 p.

[8] Chandrasekhar, S. (1961) Hydrodynamic and Hydromagnetic Stability. Oxford University Press, New York.

[9] Malik, S.K. and Singh, M. (1988) Nonlinear Focusing and the Kelvin-Helmholtz Instability in Ferrofluid/Nonmagnetic Fluid Systems. Physics of Fluids, 31, 1069.

http://dx.doi.org/10.1063/1.867018

[10] El-Dib, Y.O. (1996) Nonlinear Stability of Kelvin-Helmholtz Waves in Magnetic Fluids Stressed by a Time-Dependent Acceleration and a Tangential Magnetic Field. Journal of Plasma Physics, 55, 219.

http://dx.doi.org/10.1017/S0022377800018808

[11] El-Sayed, M.F. (2002) Effect of Variable Magnetic Field on the Stability of a Stratified Rotating Fluid Layer in Porous Medium. Czechoslovak Journal of Physics, 50, 607.

http://dx.doi.org/10.1023/A:1022854217365

[12] Bhatia, P.K. and Sharma, A. (2003) KHI of Two Viscous Superposed Conducting Fluids. Proceedings of the National Academy of Sciences, India, 73, 497.

[13] Babchin, A.J., Frenkel, A.L., Levich, B.G. and Shivashinsky, G.I. (1983) Nonlinear Saturation of Rayleigh-Taylor Instability in Thin Films. Physics of Fluids, 26, 3159.

http://dx.doi.org/10.1063/1.864083

[14] Rudraiah, N., Mathad, R.D. and Betigeri, H. (1997) The RTI of Viscous Fluid Layer with Viscosity Stratification. Current Science, 72, 391.

[15] El-Dib, Y.O. and Matoog, R.T. (2005) Electrorheological Kelvin-Helmholtz Instability of a Fluid Sheet. Journal of Colloid and Interface Science, 289, 223-241.

http://dx.doi.org/10.1016/j.jcis.2005.03.054

[16] Chavaraddi, K.B., Awati, V.B. and Gouder, P.M. (2015) Effect of Surface Roughness on Kelvin-Helmholtz Instability in Presence of Magnetic Field. International Journal of Engineering Sciences & Research Technology (IJESRT), 4, 525-534.

[17] Chavaraddi, K.B., Katagi, N.N., Awati, V.B. and Gouder, P.M. (2014) Effect of Boundary Roughness on Kelvin-Helmholtz Instability in Couple Stress Fluid Layer Bounded Above by a Porous Layer and Below by Rigid Surface. International Journal of Computational Engineering Research, 4, 35-43.