Skip to content
BY 4.0 license Open Access Published by De Gruyter May 22, 2018

Effect of boundary roughness on nonlinear saturation of Rayleigh-Taylor instability in couple-stress fluid

  • Vishwanath B. Awati , Krishna B. Chavaraddi EMAIL logo and Priya M. Gouder
From the journal Nonlinear Engineering

Abstract

The boundary roughness effects on nonlinear saturation of Rayleigh-Taylor instability (RTI) in couple-stress fluid have been studied using numerical technique on the basis of stability of interface between two fluids of the system. The resulting fourth order ordinary nonlinear differential equation is solved using Adams-Bashforth predictor and Adams-Moulton corrector techniques numerically. The various surface roughness effects and surface tension effects on nonlinear saturation of RTI of two superposed couple-stress fluid and fluid saturated porous media are well investigated. At the interface, the surface tension acts and finally stability of the problem is discussed in detail.

1 Introduction

The growth rate of instability due to gravity associated with a dense layer overlying lighter layer, Rayleigh-Taylor instability (RTI) depends on constitutive law related to stress and strain rate. The upper fluid has a higher density compared to lower fluid, the interface is unstable to small perturbations whose amplification is well described by linear stability theory which depends on density ratio, gravity, coefficient of surface tension and viscosity (see Chandrasekhar [1]).

As the interface is destabilized by the gravity, configuration is unstable, but if there exist surface tension between two fluids, the configuration will have stabilizing effect. The abundant literature is available on linear RTI in composite two-fluid system is treated as a non-porous and fluid as Newtonian fluids. Stokes [2] has framed the theory of couple-stress for non-Newtonian fluid. The small amount of supplements in a lubricant can improve the bearing performance by increasing lubricant viscosity and producing an increase in load carrying capacity. This supplement in a lubricant decrease the coefficient of friction and also increases the temperature range in which bearing can operate. Rudraiah et al. [3] discussed the dispersion model of Gill and Sankarsubramanian [4], the effect of couple-stress is to facilitate the loss of hemoglobin, a disease known as hemolysis in blood. To comprehend the control of such situation, surface instability of RTI is studied. Sunil et al. [5] discussed the RTI of two superposed couple-stress fluids with unchanging densities in a horizontal magnetic field, they have carried out stability analysis for two highly viscous fluids of equal kinematic viscosity and equal couple-stress kinematic viscosity. Rudraiah et al. [6] analyzed the couple-stress fluid effects on the control of RTI at interface between dense fluids accelerated by lighter fluid and used the approximations to derive the growth rate of RTI. Another important application of RTI is flow past a porous layer is the design of effective porous wall insulation. Wetting of porous insulation has been noticed to persist despite many preventive devices discussed by Ognewicz and Tien [7]. This attributes to unavoidable formation of openings arising from faulty insulation and aging. The wetting phenomenon is generally believed to cause an important deterioration in the performance of wall insulation.

Although abundant literature is available on linear RTI, the work on nonlinear RTI is very scarce. The discovery of high-speed computing facility and nonlinear approach has changed significantly. If the perturbed quantities are insignificant in comparison with basic state then the product of perturbed quantities can be neglected which leads to linear stability theory. On the other hand if the perturbations are not small then their product cannot be ignored this leads to nonlinear theory. Let us give a brief review on the literature of nonlinear RTI. A mathematical model is used for both analytical and numerical approaches, of RTI [1]. McCrory et al. [8] discussed the simulation of RTI of abrasively accelerated thin-shell fusion targets and showed that nonlinear evolution exhibits spike amplitude saturation. A simple mathematical model is derived by Baker and Freeman [9] for nonlinear evolution of RTI. Babchin et al. [10] the nonlinear saturation of RTI in thin films and found that combined action of flow shear and surface tension is essence of saturation mechanism. A method of strained co-ordinates used for investigating nonlinear RTI problem given by Shivamoggi [11]. After incorporating the rectifications pointed out by Malik and Singh [12] and found a revised face for nonlinear cut-off wave number which separates the region of stability from that of instability. Mohamed and Shehawey [13] considered nonlinear electrohydrodynamic RTI in the absence of surface charges and a charge free surface separating two semi-infinite dielectric fluids inclined by a normal electric field subjected to nonlinear deformations. Allah and Yahia [14] analyzed the nonlinear RTI in the presence of magnetic field and also the heat and mass transfer using simplified formulation. Rudraiah et al. [15] discussed the various problems of nonlinear study of RTI in thin film past a porous layer.

Recently Rudraiah et al. [16] examined the non-linear study of electrohydrodynamic RTI in a composite fluid-porous layer. This problem has great impact of slip velocity at the interface between porous layer and thin film. Chavaraddi et al. [17] analyzed the hydrodynamic RTI in a composite fluid-porous layer. The main purpose of this article is to demonstrate such Rayleigh-Taylor instability in a fully nonlinear with effect of boundary roughness. We consider two-dimensional nonlinear saturation of RTI on two superposed couple-stress fluid and fluid saturated porous layer using slip condition proposed by Saffman [18] with the effect of boundary roughness.

The paper is organized as follows; the basic equations and corresponding boundary conditions are discussed in Section 2 using approximation proposed by Babchin et al. [10]. Section 3 is devoted to the study of nonlinear saturation of RTI in superposed composite couple-stress fluid and fluid saturated porous layer with the effect of boundary roughness. Finally results and discussions are drawn in Section 4.

2 Mathematical Formulation

Let us consider two superposed of fluids are compressible, immiscible and inviscid flowing down an inclined plane as shown in Figure 1. Cartesian Cartesian co-ordinate system aligned with the plane describes couple-stress fluid, x points in horizontal direction and y is vertical to it. The fluid properties of two layers are different; we use subscripts to distinguish them as ρf and ρp which denote densities in upper and lower fluids respectively.

Fig. 1 The physical configuration of the problem
Fig. 1

The physical configuration of the problem

The system involving of two semi-infinite inviscid, incompressible couple-stress fluid and fluid saturated porous media differentiated by an interface in the existence of surface roughness is considered. Since gravity is present; the effective acceleration is in the positive y-direction (upwards). Therefore, light fluid thrusts on heavy fluid (fluid saturated porous media). As far as the interface amid fluids remains uniform, that is completely horizontal and perpendicular to effective acceleration. The light fluid (couple stress fluid) has adequate pressure to hold heavy fluid against ceiling. However, small deviations are bound to occur at the interface. Whenever irregularities occur the portions of interface lie higher than average and therefore feel more pressure from light fluid is necessary to support water. Therefore, interface continues to increase in the spots. In the portions where the interface has dropped by small amounts below the average, more pressure is required for the support, and interface drops further. Even the smallest perturbation will cause the instability. If situation were reversed, heavy fluid layer was on the floor supporting the volume of light fluid then pressure gradients reverse. In the portions of interface have risen, the dense fluid will lower back to average, whereas the lowered portions will rise again. This qualitative analysis leads to conclusive statement that, if heavy fluid pushes the light fluid then the interface is stable and if light fluid pushes the heavy fluid then the interface is unstable and unstable case is RTI.

The upper fluid has density ρp, the lower one has density ρf and surface tension exist between two fluids is denoted by y.

To understand the physics of problem described here, it is simplified by following assumptions of Stokes and lubrication approximations (Babchin et al. [10]) as

  1. The clear dense liquid is homogeneous and isotropic.

  2. The thickness of film h is much lesser than thickness H of porous layer bounded above the film.

    i.e., hH

  3. The Strouhal number S is assumed to be considerably small.

  4. The surface elevation η is assumed to be small compared to film thickness h,

    i.e., ηh

  5. The fluid viscosity and thermal conductivity are assumed to be constants.

Based on the above assumptions and approximations, the basic equation leads to

q=0(2.1)
ρqt+q.q=p+μ2qλ4q(2.2)

where q⃗ = (u, v) is the fluid velocity, λ is the couple-stress parameter, p is the pressure, px is the pressure gradient, μ is the fluid viscosity and ρ is the fluid density of the fluid.

The relevant boundary conditions becomes

  1. Roughness condition:

    β1uy=uaty=0.(2.3)
  2. Slip condition:

    uy=αk at y=h.Saffman[17](2.4)
  3. Kinematic condition:

    v=ηt+uηxaty=h.(2.5)
  4. Dynamic condition:

    p=ηγ2ηx2aty=h.(2.6)
  5. Couple-stress conditions:

    2uy2=0aty=0andh.(2.7)

Let us non-dimensionalize the above equations by using

x=xh,y=yh,u=uδh2/μ,v=vδh2/μ,p=pδh,β1=β1h,σ=hk,ηηh.(2.8)

The above assumptions and approximations (i.e.Stokes and lubrication approximations), that the heavy fluid in porous layer is almost stationary because of creeping flow approximation and substituting Eq. (2.8) into Eqs. (2.1) – (2.2), we get (after neglecting the asterisks for simplicity)

0=ux+vy(2.9)
0=px+2uy2M024uy4(2.10)
0=py(2.11)

where M0=λ/μh2 is the couple-stress parameter.

3 Dispersion Relation

To justify the use Stokes and lubrication approximations and derive the closed form of dispersion relation from appropriate governing equation (2.10) subjected to boundary and surface conditions.

  1. Boundary roughness condition

    β1uy=uaty=0.(3.1)
  2. Saffmann slip condition

    uy=βuaty=1.(3.2)

    where β1 is the roughness parameter, σp=hk the porous parameter and β = αpσp.

  3. Couple-stress conditions:

    2uy2=0aty=0and1.(3.3)
  4. Kinematic condition:

    v=ηt+uηxaty=1.(3.4)
  5. Dynamic condition:

    p=η1B2ηx2aty=1.(3.5)

    where B = δh2/y is the Bond number and η = η(x, y, t) is the elevation of interface.

    The solution of Eq. (2.10) subject to the boundary conditions (3.13.3)

    u={C1cosh(y/M0)+C2sinh(y/M0)+12y2+M02+Ay+B}P(3.6)

    where P=px,a1=1M0 sinh(1/M0) + β cosh(1/M0), a2 = 1M0cosh(1/M0)+βsinh(1/M0),C1=M02,C2 = M02[cosh(1/M0)1]sinh(1/M0)

    A={122+β(1+2M02)+M02a1+M02[1cosh(1/M0)]sinh(1/M0)a2β1M0[1cosh(1/M0)]sinh(1/M0)M02β}/1+β(1β1)B=M02β1M0cosh(1/M0)1sinh(1/M0)Aβ1.

    Integrating Eq. (2.9) with respect to y between y = 0 and 1 and using Eq. (3.6), we get

    v(1)=v1=01uxdy=2px2N.(3.7)

    where N=2M03{1cosh(1/M0)}sinh(1/M0)+16(1+M02)+A2+B.

    By using Eqs. (3.5) and (3.7) then Eq. (3.4) becomes

    ηt+N1ηx+1B3ηx3ηx=N2ηx2+1B4ηx4.(3.8)

    where N1=C1Cosh(1/M0)+C2sinh(1/M0)+12+M02+A+B. Eq. (3.8) is a fourth order nonlinear ordinary differential equation and now let us analyzes the interface evolution by this equation. The process described here is quite different from a process wherein the film is bounded by a fluid with moving boundaries instead of porous layer discussed by Babchin et al. [10]. Eq. (3.8) is not amenable to solve by using analytical method and hence it can be solved numerically using 4th order central differences in space and time with periodic boundary conditions as explained below. For time- integration of Eq. (3.8), we use Adams-Bashforth predictor and Adams-Moulton corrector methods of 4th order [19]. For the spatial derivatives are described by the following central difference formulae of fourth order accuracy

    ηx112Δxηi28ηi1+8ηi+1ηi+2.(3.9)
    2ηu2112Δx2ηi2+16i130ηi+16ηi+1ηi+2.(3.10)
    3ηx318Δx3ηi38ηi2+13ηi113i+1+8ηi+2ηi+3.(3.11)
    4ηx416Δx4ηi3+12ηi239ηi1+56ηi39ηi+1+12ηi+2ηi+3.(3.12)

Here η (i − 2) stands for value of η at the position x − 2Δ x. The integer i indicates the ith grid point. The initial condition used in numerical integration is a sine-wave with wave number and is of the form

ηx,0=η0sinx0<xπ2.(3.13)

Here the amplitude η0 is assumed to be small. In our numerical computation η0 = 10− 4, the non-dimensional form and periodic boundary conditions have been applied in x-direction.

4 Results and Discussion

The boundary roughness effects on nonlinear saturation of Rayleigh-Taylor instability (RTI) in a couple-stress fluid have analyzed based on stability of the interface of system between two fluids. The governing fourth order nonlinear ordinary differential Eq. (3.8) is solved numerically on the interval 0<xπ2 with periodic boundary conditions using Adams-Bashforth predictor and Adams-Moulton corrector method for space and time. This method involves the approximation of spatial derivatives with central finite difference scheme for time step and consequences are depicted in Figs. 27. The process of study is to explain the distortion of wave profile i.e. η(x, t). It can also be represented in the following words ( in our interface interpretation of equation); every point of interface moves in horizontal x-direction with velocity which is proportional to that point and elevation η . Figures 24 predicts the spatial structure of growth rate of interface in terms of nNLnNL,Max=1ηηt1ηmaxηmaxt at an early stage (initial line) before instability occurs for wave number = 0.75 and other parameters defined earlier. Also, Figures 24 shows six waves are contained in interval for the instance of small wave numbers and peaks reflect the position of wave modes. Figure 5 and 6 predicts that whether the full numerical solutions which deviate from a simple harmonic behaviour of η, namely η = η0 sin(ℓ x) for = 0.75, where only a portion of the wave is shown.

Fig. 2 Shape of interface growth rate shape of interface growth rate for a roughness parameter β1 = 0, 3.3 × 10−3, 3.3 × 10−2 with fixed values of M0 = 0.3, B = 0.02 and β = 2. (from top to bottom)
Fig. 2

Shape of interface growth rate shape of interface growth rate for a roughness parameter β1 = 0, 3.3 × 10−3, 3.3 × 10−2 with fixed values of M0 = 0.3, B = 0.02 and β = 2. (from top to bottom)

Fig. 3 Different values of M0 = 0.2, 0.25, 0.3, and 0.35 with fixed values of B = 0.02, β = 2, ℓ = 0.75 and β1 = 3.3 × 10− 6.
Fig. 3

Different values of M0 = 0.2, 0.25, 0.3, and 0.35 with fixed values of B = 0.02, β = 2, = 0.75 and β1 = 3.3 × 10− 6.

Fig. 4 Different values of B = 0.01, 0.02, 0.03 and 0.04 with fixed values of β = 2, M0 = 0.3, ℓ= 0.75 and β1 = 3.3 × 10–6.
Fig. 4

Different values of B = 0.01, 0.02, 0.03 and 0.04 with fixed values of β = 2, M0 = 0.3, = 0.75 and β1 = 3.3 × 10–6.

Fig. 5 Comparison between full numerical wave number of ℓ̂ = 0.75 and three roughness. Parameters β1 = 0, 3.3 × 10– 3, 3.3 × 10– 2 with fixed values of M0 = 0.3, B = 0.02 and β = 2. (from top to bottom)
Fig. 5

Comparison between full numerical wave number of ℓ̂ = 0.75 and three roughness. Parameters β1 = 0, 3.3 × 10– 3, 3.3 × 10– 2 with fixed values of M0 = 0.3, B = 0.02 and β = 2. (from top to bottom)

Fig. 6 Comparison between full numerical wave number of ℓ̂ = 0.75 and the couple stress. Parameter M0 = 0.2, 0.25, 0.3, 0.35 with fixed values of β = 2, B = 0.02 and β1 = 3.3 × 10– 6. (from top to bottom)
Fig. 6

Comparison between full numerical wave number of ℓ̂ = 0.75 and the couple stress. Parameter M0 = 0.2, 0.25, 0.3, 0.35 with fixed values of β = 2, B = 0.02 and β1 = 3.3 × 10– 6. (from top to bottom)

Fig. 7 Comparison between full numerical wave number of ℓ̂ = 0.75 and the Bond number B = 0.01, 0.02, 0.03, 0.04with fixed values of M0 = 0.3, β1 = 3.3 × 10– 6and β = 2. (from top to bottom)
Fig. 7

Comparison between full numerical wave number of ℓ̂ = 0.75 and the Bond number B = 0.01, 0.02, 0.03, 0.04with fixed values of M0 = 0.3, β1 = 3.3 × 10– 6and β = 2. (from top to bottom)

Figures 5(a)5(c) depicts that spatial structures of interface for t>0 and Figures 2(a)2(c) describes interface at t = 0. Initial profile is symmetric in the interval 0 ≤ x ≤ 16π . But for t > 0, we notice that symmetry is lost, since the maximum moves farther from one of zeros and closer to other one. However, this process of steepening of forward faces of interface profile does not result in breakup of interface because of the effects of surface tension, roughness parameter and couple-stress parameter. Initially these parameters may have negligible effects but as time progresses these parameters play a vital role in stabilization even though negligible effect in reducing the asymmetry of system by roughness parameter. In Figures 5(a)5(c) only a portion of the wave is shown as it move from top to bottom and clearly notice that symmetry can be obtained with effect of roughness parameter, hence reduces the growth rate of RTI at interface considerably for t > 0.

In Figures 6(a)6(d) it is clear that interface profile is symmetric in the interval 0 ≤ x ≤ 16π for initial time t = 0 and as it moves from Figures 5(a)5(c) (top to bottom), it is observed that symmetry can be recovered slowly for increasing small roughness parameter β1. Therefore, it has negligible effect on reducing asymmetry of the system because of resistance offered by surface roughness in that process a part of kinetic energy is converted into potential energy.

Also, Figures 4(a)4(d) shows the interface is symmetric in the interval 0 ≤ x ≤ 16π for t = 0. For t > 0, the initial symmetry can be lost as shown in Fig. 6(a) and increasing Bond number the symmetry can be recovered and results are depicted in Figs. 7(b)7(d). Hence for suitable choice of Bond number the symmetry can maintained at the interface and hence reduce the growth rate of RTI at interface. It may also be noted that full numerical solution for all the cases is not possible due to limitation of numerical scheme.

Acknowledgement

The authors (VBA and PMG) wish to thank respectively the VC / Principal of their institutions for their encouragement and support in doing research. Also the author (KBC) wish to thank VGST, Department of Science & Technology, Government of Karnataka under Project grant no. VGST/SMYSR (2014-15)/GRD-433/2015-16, Dated: 28.04.2015 for their support and encouragement.

References

[1] S. Chandrasekhar, Hydrodynamic and Hydro Magnetic Stability, Dover Publication, New York, (1961).Search in Google Scholar

[2] V.K. Stokes(1966), Couple stresses in Fluids, Phys. Fluids. 9, 1709.10.1007/978-3-642-82351-0_4Search in Google Scholar

[3] N. Rudraiah, S. R. Kasiviswanath and P.N. Kaloni(1991), Generalised dispersion in a synovial fluid of human joints, Biorehology 28, 207-209.10.3233/BIR-1991-283-412Search in Google Scholar

[4] Gill and Sankarsubramanian(1970), Exact analysis of unsteady convective diffusion, Proc. Roy. Soc. A316, London, 347.Search in Google Scholar

[5] Sunil, R.C. Sharma and R.S. Chandel(2002), On superposed couple-stress fluids in porous medium in hydro magnetics, Z. Naturforsch 57a, 955-960.10.1515/zna-2002-1208Search in Google Scholar

[6] N. Rudraiah and G. Chandrashekara(2010), Effect of couple-stress on the growth rate of Rayleigh-Taylor instability at the interface in a finite thickness couple-stress fluid, Journal of Applied Fluid Mechanics, 3(1), 83-89.10.36884/jafm.3.01.11882Search in Google Scholar

[7] Y. Ognewicz and C.L. Tien(1981), Int. J. Heat and Mass Transfer,421-429, 1981.10.1016/0017-9310(81)90049-1Search in Google Scholar

[8] R.L. McCrory, L. Montierth, R.L. Morse and C.P. Verdon(1981), Nonlinear evolution of ablation-driven Rayleigh-Taylor instability, Phys. Rev. Lett.,46, 336-339.10.1103/PhysRevLett.46.336Search in Google Scholar

[9] L. Baker and J.R. Freeman(1981), Heuristic model of the nonlinear Rayleigh-Taylor instability, J. Applied Phys, 52(2), 655-663.10.1063/1.328793Search in Google Scholar

[10] A.J. Babchin, A.L. Frenkel, B.G. Levich and G.I. Shivashinsky(1983), Nonlinear saturation of Rayleigh-Taylor Instability in thin films, Phys. Fluids, 26, 3159-3161.10.1063/1.864083Search in Google Scholar

[11] B.K. Shivamoggi(1983), Comments on nonlinear Rayleigh-Taylor instability, Astrophys. Space Sci, 97(1), 221-221.10.1007/BF00684626Search in Google Scholar

[12] S.K. Malik and M. Singh(1983), Non-linear Rayleigh-Taylor instability, Astrophys. Space Science, 90(1), 45-50.10.1007/BF00651545Search in Google Scholar

[13] A.E.M.A. Mohamed and E.S.F. El-Shehawey(1983), Nonlinear electro hydrodynamic Rayleigh-Taylor instability-II: A perpendicular field producing surface charge, Phys. Fluids, 26, 1724-1730.10.1063/1.864371Search in Google Scholar

[14] M.H.O. Allah and A.A. Yahia(1991), Nonlinear Rayleigh-Taylor instability in the presence of magnetic field and mass and heat transfer, Astrophys and Space Science, 181, 183-188.10.1007/BF00639087Search in Google Scholar

[15] N. Rudraiah, C. Wagner, G. S. Evans and R. Friedrich (1998), Nonlinear study of Rayleigh-Taylor instability in thin films past a porous layer, Indian J. Pure Appl. Math., 29(4), 417-429.Search in Google Scholar

[16] N. Rudraiah, K.B. Chavaraddi, I.S. Shivakumara (2011), Nonlinear study of electrohydro- dynamic Rayleigh–Taylor instability in a composite fluid–porous layer, International Journal of Non-Linear Mechanics, 46, 1009-1016.10.1016/j.ijnonlinmec.2011.04.017Search in Google Scholar

[17] K. B. Chavaraddi, I.S. Shivakumara, N. Rudraiah(2014), Nonlinear Study of Hydrodynamic Rayleigh-Taylor Instability in a Composite Fluid-Porous Layer, Gen. Math. Notes, Vol. 25, No. 1, November 2014, pp. 75-90.Search in Google Scholar

[18] P.G. Saffman(1971), On the Boundary Condition at the Surface of a Porous Medium, Studies in Appl. Math., 50, pp. 93-101.10.1002/sapm197150293Search in Google Scholar

[19] S.C. Chapra and R.P. Canale, Numerical Methods for Engineers (2ndEdition), McGraw Hill Pub. Corp., New York,1998.Search in Google Scholar

Received: 2017-10-09
Accepted: 2018-03-17
Published Online: 2018-05-22
Published in Print: 2019-01-28

© 2019 V.B. Awati et al., published by De Gruyter.

This work is licensed under the Creative Commons Attribution 4.0 Public License.

Downloaded on 29.3.2024 from https://www.degruyter.com/document/doi/10.1515/nleng-2018-0031/html
Scroll to top button