Next Article in Journal
A Geometric Perspective on the Modulation of Potential Energy Release by a Lateral Potential Vorticity Gradient
Next Article in Special Issue
On the Optimal Control of Stationary Fluid–Structure Interaction Systems
Previous Article in Journal
The Propagation of Nonlinear Internal Waves under the Influence of Variable Topography and Earth’s Rotation in a Two-Layer Fluid
Previous Article in Special Issue
Acceleration Waves in Rational Extended Thermodynamics of Rarefied Monatomic Gases
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability of the Plane Bingham–Poiseuille Flow in an Inclined Channel

by
Paolo Falsaperla
,
Andrea Giacobbe
and
Giuseppe Mulone
*,†
Department of Mathematics and Computer Science, University of Catania, 95125 Catania, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Fluids 2020, 5(3), 141; https://doi.org/10.3390/fluids5030141
Submission received: 25 July 2020 / Revised: 20 August 2020 / Accepted: 25 August 2020 / Published: 27 August 2020
(This article belongs to the Special Issue Classical and Modern Topics in Fluid Dynamics and Transport Phenomena)

Abstract

:
We study the stability of laminar Bingham–Poiseuille flows in a sheet of fluid (open channel) down an incline with constant slope angle β ( 0 , π / 2 ) . This problem has geophysical applications to the evolution of landslides. In this article, we apply to this problem recent results of Falsaperla et al. for laminar Couette and Poiseuille flows of Newtonian fluids in inclined channels. The stability of the basic motion of the generalised Navier–Stokes system for a Bingham fluid in a horizontal channel against linear perturbations has been studied. In this article, we study the flows of a Bingham fluid when the channel is oblique and we prove a stabilizing effect of the Bingham parameter B. We also study the stability of the linear system with an energy method (Lyapunov functions) and prove that the streamwise perturbations are always stable, while the spanwise perturbations are energy-stable if the Reynolds number Re is less than the critical Reynolds number R c obtained solving a generalised Orr equation of a maximum variational problem.

Graphical Abstract

1. Introduction

The stability of laminar flows in an inclined channel is important in many geophysical and industrial applications. Recently, Falsaperla et al. [1,2] described the stability of laminar flows of the steady solutions of a sheet of a Newtonian fluid (open channel) both in the horizontal case [1], and down an incline with constant slope angle β [2]. In a reference frame O ,   i ,   j ,   k with i parallel to the incline and k orthogonal to the plane of the channel, let z be the coordinate of the axis orthogonal to the incline, the basic steady motion is given by the velocity field U ( z ) i and the pressure p ( z ) . For a laminar Newtonian fluid in a horizontal layer, U ( z ) is linear for Couette flows and parabolic for Poiseuille flows, while, for an inclined channel, U ( z ) is a parabolic function which vanishes at the bottom of the channel and whose derivative with respect to z vanish at the top.
Choosing an appropriately weighted L 2 -energy equivalent to the classical energy norm, we prove in [1] that the plane Couette and Poiseuille flows are nonlinearly stable with respect to streamwise perturbations for any Reynolds number (see also [3]). We also prove in [4] that the plane Couette and Poiseuille flows are nonlinearly stable (and linearly stable in the energy norm) if the Reynolds number is less then the critical value
Re O r r / sin φ ,
when the perturbation is a tilted perturbation in the x direction, with x a coordinate associated to a new reference frame with i = i cos φ + j sin φ , j = j cos φ i sin φ , and k = k . The number Re O r r is the Orr [5] critical Reynolds number for spanwise perturbations which, for the Couette flow, is Re O r r = 177.22 and for the Poiseuille flow is Re O r r = 175.31 . The results in [1] improve those obtained by Joseph [6], who found for streamwise perturbations a critical nonlinear value of 82.6 in the plane Couette case, and those obtained by Joseph and Carmi [7], who found the value 99.1 for plane Poiseuille flow for streamwise perturbations. The results we obtained in [1] are in a good agreement with the experiments and the numerical simulations. In a recent paper [2], we studied the stability of laminar Poiseuille flows for a Newtonian fluid in an oblique channel. That work is a preliminary investigation to model debris flows down an incline, and this work is its natural application to fluids which model those of realistic mudflows.
De Blasio, in his monograph, “Introduction to the Physics of Landslides” [8] introduces, in Chapter 4, non-Newtonian fluids, mudflows, and debris flows. In particular, he considers the behaviour of a mudflow described by Bingham rheology. He gives a first example of non-Newtonian fluid flow pattern considering a laminar stationary flow of a mud layer of thickness D along an inclined plane of infinite length [8] (pp. 100–103). He obtains the basic motion, but he does not investigate its stability. Basic flows of Bingham type in inclined channels arise also in other settings, for instance lubrication flows [9,10].
Bingham fluids are typically applied to drilling muds. They are subjected to a yield stress in addition to a plastic viscosity. The rheological properties are well-presented in the introduction of the paper by Frigaard et al. [11]. The stability (both in the linear and nonlinear case) of the horizontal laminar Poiseuille flows of a Bingham fluid has been investigated during the last 15 years by C. Nouar, I.A. Frigaard and co-workers [11,12,13,14,15,16]. In [12] the authors write: “Parallel duct flows of slurries and suspensions are relatively common in the petroleum and mining industries, where prediction of the flowregime can be an important design parameter in hydraulic systems. These complex fluids are commonly visco-plastic, meaning that they are characterised rheologically by having a yield stress, below which no deformation takes place”. Among other results, in [11,12,13,14,15,16], the plug (unyielded) zone and the yielded zone are introduced, and it is observed that the basic motion is the Bingham–Poiseuille flow, with basic velocity field U ( z ) parabolic in the yielded zone and constant in the plug zone. In [14] the authors determine the yield surface and write the perturbation equations to the basic motion. They also linearise around the basic motion, and analyse the long-time behaviour of the disturbance which is transformed in an eigenvalue problem. They finally solve the eigenvalue problem with the appropriate boundary conditions using a Chebyshev collocation method, and investigate the dependence of the Reynods number on the Bingham parameter. In particular, they obtain the linear stability of the basic motion and prove that the Bingham number has a stabilizing effect.
The main aim of this work is to study the stability of the basic motion of the Bingham fluid down an inclined open channel as given in De Blasio, [8]. In particular, we study the stability of the linearised system of perturbations, and generalise to Bingham flows the stability results obtained for Newtonian fluids in [2]. Precisely, we prove that the basic motion is linearly stable for any Reynolds number and we prove (as in [11,12,13,14,15,16]) a stabilizing effect of the Bingham parameter. We also consider the linear stability studying the evolution in time of the energy of the perturbations, which amounts of investigating the stability in the energy norm with the Lyapunov method introduced by Reynolds [4] and Orr [5]. We prove that the streamwise perturbations are always stable while the spanwise perturbations are energy-stable if the Reynolds number Re is less than a critical Reynolds number R c obtained by solving a generalised Orr equation of a maximum variational problem.
We also observe that Allouche et al. [17] focussed on Newtonian and generalised Newtonian (in particular, Carreau) fluid film flows down an incline and compared the thresholds of oblique waves to those of two-dimensional waves for a given slope in order to reach the dominant instability in a given flow configuration. In particular they consider the shear-thinning and the shear-thickening cases and performed a temporal linear stability study on the problem. However, they do not study the stability of Bingham fluids.
In this article, we consider the stability in the linear context (both numerically and in energy). In a future paper we plan to investigate the nonlinear stability of Bingham flows and the stability of flows that obey to the Herschel–Bulkley rheology, which is even more appropriate in applications to landslides. We also plan to investigate, as in [1], the stability with respect to oblique perturbations and its relation with the experiments.
The layout of the paper is the following. In Section 2, following [8], we introduce the governing equations for a laminar open flow of a Bingham fluid down an inclined channel. In Section 3, following [11,12,13,14,15,16], we write the linearised equations for a perturbation to the basic motion U ( z ) with the appropriate boundary conditions, and we solve the generalised Sommerfeld equations for the eigenvalue problem with the Chebyshev collocation method. In Section 4 we study the stability of the linearised perturbations with the Lyapunov method by using an energy norm, we find sufficient conditions for stability in energy for the streamwise and the spanwise perturbations and, for the latter, we write the Euler–Lagrange equations for the maximum problem and we solve the related generalised Orr equation. In Section 5, we draw some conclusions.

2. Laminar Open Flow of a Bingham Fluid Down an Inclined Channel

Consider the stationary laminar motion of a mud layer of thickness D along an inclined plane of infinite length with constant slope angle β , 0 < β < π / 2 [8] (p. 100) in a reference frame O ,   i ,   j ,   k with associated coordinates x , y , z . The motion takes place in a channel D 1 = R x , y 2 × [ 0 , D ] z of thickness D. The x-axis is taken along the slope direction while the z-axis is perpendicular to the slope. The y-axis is also parallel to the slope and orthogonal to the slope direction x. The channel extends indefinitely in the x , y directions and has a finite depth D in the z direction (see Figure 1).
We consider the flow of an incompressible Bingham fluid with yield stress τ 0 and plastic viscosity (Bingham viscosity) μ B . The wall placed in z = 0 is rigid (the velocity of the flow vanishes) while that placed at z = D is stress-free (the shear stress is zero). The governing equations are (cf. [11,12,13,14,15,16]):
U t + U · U = g 1 ρ p ¯ + 1 ρ · τ ( U ) , · U = 0 ,
where U is the velocity field, ρ is the constant density, p ¯ is the pressure, τ is the deviatoric extra-stress tensor, ∇ is the gradient operator, and g is the gravity, g = g sin β i g cos β k . The velocity vector is of the form U = U i + V j + W k , where U, V, W are the velocity components, and i , j , k are unit vectors in the streamwise (x) and spanwise (y) directions, and normal to the wall (z), respectively. As in [11,12,13,14,15,16], we assume the constitutive relations between the stress and the strain rate for a Bingham fluid
τ = η γ ˙ τ τ 0 γ ˙ = 0 τ < τ 0 ,
where γ ˙ and τ are the strain rate and deviatoric stress tensors respectively. The yield stress is denoted by τ 0 and η denotes the effective viscosity, defined by
η = μ B + τ 0 γ ˙ 1 ,
where μ B is the Bingham viscosity (see De Blasio [8]), also called limiting viscosity (cf. [11]), τ 0 is the yield stress, γ ˙ and τ are respectively the second invariants of the strain rate and the deviatoric stress tensors. They are given by γ ˙ = γ ˙ i j γ ˙ i j / 2 , τ = τ i j τ i j / 2 , γ ˙ i j = U i , j + U j , i , and the Einstein convention on summation of repeated indices has been assumed.
Following De Blasio [8], we assume that “the thickness of the mud flow does not change along the x coordinate, hence, U and p ¯ are independent of x”. So, we consider laminar shear flows U = U ( z ) i , p ¯ = p ¯ ( z ) which are solutions to the stationary equations
U · U = g 1 ρ p ¯ + 1 ρ · τ ( U ) , · U = 0 .
Equation (2) becomes
g sin β + 1 ρ τ x z z = 0 1 ρ p ¯ z + g cos β = 0 .
Solving these equations, we have
τ x z = ρ g sin β z + A ¯ , p ¯ = ρ g cos β z + B ¯ ,
where A ¯ and B ¯ are integration constants that can be determined by imposing the chosen boundary conditions, that, in our case, are vanishing shear stress at the upper surface of the mud, i.e., z U ( D ) = 0 , and pressure equals to the atmospheric pressure p a t m . Applying Bingham Equation (1) and denoting τ = τ x z , we have
μ B U z = τ τ 0 = ρ g sin β ( D z ) τ 0 , if τ τ 0 , U z = 0 , if τ τ 0 .
Following De Blasio [8] (p. 101), the first of two, Equation (3), can now be integrated from the base up to the level where the right-hand side becomes zero. This defines a layer known as the shear layer (the yielded zone) whose thickness D S can be found from Equation (3) imposing τ = τ 0 . We find
z = D S = D τ 0 ρ g sin β .
From Equation (3) and the definition of D S we have
μ B U z = τ τ 0 = ρ g sin β ( D z ) τ 0 , if 0 z D S U z = 0 , if z D S .
We note that for z D S the velocity U ( z ) increases with height. For z D S , as De Blasio observes, the weight of the overlying material is not sufficient to shear the material, as it is obvious from the second of Equation (4). As a consequence, the material moves rigidly, U ( z ) = c 0 ( c 0 a constant) like a plug. The thickness of this plug layer (unyielded zone) is
D P = D D S = τ 0 ρ g sin β .
Integrating (4) and adopting the boundary conditions
U ( 0 ) = 0 , U ( D S ) z = 0 .
We have the basic Bingham-Poiesuille flow, p ¯ = p a t m + ρ g cos β ( D z ) , and
U ( z ) = ρ g sin β μ B z D S z 2 if 0 z D S , ρ g sin β D S 2 2 μ B if z D S ,
These equations are non-dimensionalized using a length scale D, a velocity scale U 0 = U ( D S ) = ( ρ g sin β D S 2 ) / ( 2 μ B ) , and a pressure–stress scale ρ U 0 2 . Introducing the non-dimensional Reynolds and Bingham numbers
Re = ρ U 0 D μ B , B = τ 0 D μ B U 0 ,
the non-dimensional form of (3) are
τ = 1 Re U z + B , if τ B Re U z = 0 , if τ B Re ,
while the non-dimensional form of U ( z ) is
U ( z ) = 1 ( z z S ) 2 z S 2 , if 0 z z S , 1 , if z S z 1 ,
where z S = D S / D is a positive number that represents the depth of the yield surface. Such depth can be computed directly, using the velocity U 0 and the Bingham number B, and solving the equation B z S 2 = 2 ( 1 z S ) , using the momentum conservation equation in the plug region, where we have τ 0 = | τ w | , with τ w the wall shear stress | τ w |   =   τ ( z = 1 ) in the yielded region. In non-dimensional form, we have
B Re ( 1 z S ) = 1 Re U z z = 0 + B Re = B Re 1 + 2 B z S B = 2 ( 1 z S ) z S 2 z S = 2 1 + 1 + 2 B .
We observe that in the limit when B tends to 0, the fluid becomes a Newtonian fluid, and z S tends to 1. In such a limit case there is no plug region.

3. Linear Stability Analysis

In order to perform the linear stability analysis, as in [11,12,13,14,15,16], we introduce ( ϵ u , ϵ p ) , infinitesimal perturbations to the basic flow ( U , p ¯ ) , with ϵ < < 1 and | u | and p bounded. The linearised perturbation equation are
ϵ u t + ϵ [ u · U + U · u ] = ϵ p + · [ τ ( U + ϵ u ) τ ( U ) ] , · u = 0 .
wherever the yield stress is exceeded, the effective viscosity of the perturbed flow is expanded about the basic flow (see [14]). Denoting with μ the non dimensional effective viscosity, we have
μ ( U ( z ) + ϵ u ( x , y , z , t ) ) = μ ( U ( z ) + ϵ u , ϵ v , ϵ w ) = 1 Re 1 + B A / 2 ,
where
A = γ ˙ i j ( U + ϵ u ( x , y , z , t ) ) γ ˙ i j ( U + ϵ u ( x , y , z , t ) ) .
A simple calculation shows that
A = 4 ϵ u x 2 + 2 ϵ u y + ϵ v x 2 + 2 ( U ( z ) + ϵ u ) z + ϵ w x ) 2 + 4 ( ϵ v y ) 2 + 2 ( ϵ v z + ϵ w y ) 2 + 4 ( ϵ w z ) 2 ,
where A depends on the number ϵ , ( u , v , w ) are the components of u , and the subscript indicate partial differentiation with respect to the variables. Expanding μ ( U ( z ) + ϵ u ( x , y , z , t ) ) in powers of ϵ , we have
μ ( U ( z ) + ϵ u ( x , y , z , t ) ) = μ ( U ( z ) ) + ϵ μ ( U ( z ) + ϵ u ( x , y , z , t ) ) | ϵ = 0 ϵ + O ( ϵ 2 ) .
From (6) it follows that
μ ϵ | ϵ = 0 = B 4 Re 1 ( U ) 2 | U | A ϵ | ϵ = 0 ,
where U is the derivative of U with respect to z. Finally, at the first order in ϵ we have:
μ ( U ( z ) + ϵ u ( x , y , z , t ) ) = 1 Re 1 + B | U | B | U | U ( u z + w x ) ϵ .
We recall that
τ ( U + ϵ u ) τ ( U ) = μ ( U + ϵ u ) ( γ ˙ ( U ) + ϵ γ ˙ ( u ) ) μ ( U ) γ ˙ ( U ) .
Computing γ ˙ ( U ) and by observing that U 0 , at the first order in ϵ , we have
τ ( U + ϵ u ) τ ( U ) = ϵ 1 Re 1 + B U γ ˙ ( u ) ( u z + w x ) B Re ( U ) 2 γ ˙ ( U ) .
From this relation, we easily obtain
1 ϵ τ ( U + ϵ u ) τ ( U ) i j = β i j ,
with
β i j = Re 1 1 + B U γ ˙ i j ( u ) if ( i , j ) ( 1 , 3 ) or ( 3 , 1 ) Re 1 γ ˙ i j ( u ) if ( i , j ) = ( 1 , 3 ) or ( 3 , 1 ) .
Collecting (5), (8) and (9), we have
u t + u · U + U · u = p + ( β i j ) i e j · u = 0
where we have used the convention 1 = x , 2 = y , 3 = z and e 1 = i , e 2 = j , e 3 = k . The perturbation Equation (10) in the open domain D S = R 2 × ( 0 , z S ) , written with respect to the components ( u , v , w ) of u , are
u t + w U + U u x = p x + 1 Re Δ u + B Re U ( u x x + u y y w x z ) , v t + U v x = p y + 1 Re Δ v + B Re 1 U ( v x x + v y y + v z z ) + 1 U ( v z + w y ) , w t + U w x = p z + 1 Re Δ w + B Re 1 U ( w y y + w z z + v y z + w z z ) + 2 1 U w z ,
where Δ is the Laplacian. We assume that the perturbations u , p are periodic in x and y of periods 2 π / a , 2 π / b , with wave number a and b in the directions x and y. We will later consider the limit cases of streamwise perturbations, i.e., perturbations which do not depend on x and are periodic only on y, and spanwise perturbations, i.e., perturbations which do not depend on y and are periodic only on x.
Equation (11) must be completed with a set of boundary conditions. Since, for a yield stress fluid U = 2 ( z z S ) / z S 2 and 1 / U is singular as z z S , the effective viscosity goes to + at the yield surface. It follows that some care is required in deriving the boundary conditions at z = z S as explained for example in [11,14]. As in [11], we assume that the yield surface z = z S is perturbed to z ^ = z S + ϵ h ( x , y , t ) , with h ( x , y , t ) periodic in x and y of periods 2 π / a , 2 π / b . The boundary conditions can be obtained as in [11,13,14,15,16] and turn out to be
u ( x , y , 0 , t ) = v ( x , y , 0 , t ) = w ( x , y , 0 , t ) = 0 , u ( x , y , z S , t ) = v ( x , y , z S , t ) = w ( x , y , z S , t ) = 0 , v ( x , y , z s , t ) z = w ( x , y , z s , t ) z = 0 , u ( x , y , z s , t ) z = h U ( z s ) .
The boundary conditions at z = 0 are rigid, as classically found in Newtonian fluids. The Dirichlet boundary conditions u = 0 at the yield surface can be deduced from the fact that the unyielded plug zone is constrained to move as a rigid body according to the Bingham model [13]. The Neumann conditions come from linearisation of the condition γ ˙ i j ( U + ϵ u ) = 0 at the perturbed yield surface, onto the unperturbed yield surface position [13]. In particular, following [13], we note that the conditions z v ( x , y , z s , t ) = z w ( x , y , z s , t ) = 0 , are necessary due to the singular behaviour of 1 / U on z = z S in the non-Newtonian part of (11). These boundary conditions ensure that (11) are well-defined as z z S . We thus have that v and w are asymptotic at z S to ( z z S ) α , with α > 1 . This fact will be useful in the Lyapunov investigation, in particular when using the boundary conditions in the periodicity cell Ω = [ 0 , 2 π / a ] × [ 0 , 2 π / b ] × [ 0 , z S ] .
In order to study the linear stability with the spectral method, we assume that the perturbations have the form
( u , v , w , p , h ) = [ u ( z ) , v ( z ) , w ( z ) , p ( z ) , h ] e i ( a x + b y ) + C t ,
where C is a complex number whose real part gives the decay rate of the perturbation.
We observe that, in [13], Frigaard and Nouar prove a Squire-like (see [18]) theorem: “if the equivalent eigenvalue bounds for a Newtonian fluid yield a stability criterion, then the same stability criterion is valid for the Bingham fluid flow, but with reduced wavenumbers and Reynolds numbers”. Assuming this, we consider in what follows only streamwise or only spanwise perturbations.
In the case of streamwise perturbations a = 0 and b 0 , we hence have
C u + w U ( z ) = Re 1 ( D 2 b 2 ) u Re 1 B U b 2 u C v = i b p + Re 1 ( D 2 b 2 ) v + Re 1 B 1 U ( D 2 b 2 ) v + 1 U ( v + i b w ) C w = p + Re 1 ( D 2 b 2 ) w + Re 1 B 1 U ( b 2 w + 2 w + i b v ) + 2 1 U w i b v + w = 0 .
An easy calculation proves that the streamwise perturbations are always stable (the real part of C is always less than 0). For the spanwise perturbations a 0 and b = 0 , and we have
C u + i a U u + w U ( z ) = i a p + Re 1 ( D 2 a 2 ) u Re 1 B U ( a 2 u + i a w ) C v + i a U v = Re 1 ( D 2 a 2 ) v + Re 1 B 1 U ( D 2 a 2 ) v + 1 U v C w + i a U w = p + Re 1 ( D 2 a 2 ) w + 2 Re 1 B 1 U w + 1 U w i a u + w = 0 .
The variable v appears only in the second equation of system (12), and it can proven that v cannot destabilize the basic motion. The remaining three equations depend only on u , w and on the pressure p. Using the continuity equation, one can further eliminate the variables u , p and derive a system of one equation in the variable w. This can be obtained by subtracting i a times the z-derivative of the first equation of (12) and adding the third equation of (12) multiplied by a 2 (this amounts to writing the equation of the third component of the doube-curl of u ). In this way one obtains the generalised Orr–Sommerfeld equation
( C + i a U ) ( w a 2 w ) i a w U ( z ) = Re 1 ( w 2 a 2 w + a 4 w ) 4 a 2 B Re 1 w U + 1 U w ,
with boundary conditions
w ( x , y , 0 ) = w z ( x , y , 0 ) = 0 , w ( x , y , z S ) = w z ( x , y , z S ) = 0 .
Numerical calculations show that the linear perturbations are always stable, i.e., the real part of C is negative for any Reynolds number. The Chebyshev collocation method has been used up to 100 polynomials for w. This transforms the differential equation in an algebraic generalized eigenvalue system which also takes into account the given boundary conditions. Figure 2 shows that the maximum real part of the spectrum of C is negative and decreasing with the Bingham number. Here, we have considered the significative ranges a [ 0 , 2.6 ] , Re [ 1000 , 11000 ] . We also obtain the expected fact that as Re tends to 0, the real part of C tends to .

4. Linear-Energy Stability of Perturbations

In this section we use the linear Lyapunov method, also typically called energy-stability method. In the periodicity cell Ω , we introduce the scalar product ( , ) and the L 2 ( Ω ) norm · . We have the energy identities:
d d t u 2 2 + ( w U , u ) + ( U u x , u ) = ( p x , u ) + 1 Re ( Δ u , u ) + B Re 1 U ( u x x + u y y w x z ) , u d d t v 2 2 + ( U v x , v ) = ( p y , v ) + 1 Re ( Δ v , v ) + B Re 1 U Δ v , v + 1 U ( v z + w y ) , v d d t w 2 2 + ( U w x , w ) = ( p z , w ) + 1 Re ( Δ w , w ) + B Re 1 U ( w y y + 2 w z z + v y z ) + 2 U w z , w .
These identities can be simplified by applying integration by parts and using the boundary conditions. Since the term ( U u x , u ) integrates to zero
( U u x , u ) = U , x u 2 2 = lim z ¯ z S U u 2 2 z = 0 z ¯ = 0 .
We have the simplified identities
d d t u 2 2 + ( w U , u ) = ( p x , u ) 1 Re u 2 + B Re 1 U ( u x x + u y y w x z ) , u
d d t v 2 2 = ( p y , v ) 1 Re v 2 + B Re 1 U Δ v , v + 1 U ( v z + w y ) , v
d d t w 2 2 = ( p z , w ) 1 Re w 2 + B Re 1 U ( w y y + 2 w z z + v y z ) + 2 U w z , w .

4.1. Stability of Streamwise Perturbations

To prove that the linear streamwise perturbations are always stable, we begin observing that, for streamwise perturbations, x 0 . The solenoidality of the kinetic perturbation hence reads v y + w z = 0 . From (14) and (15), we have
1 2 d d t v 2 + w 2 = 1 Re ( v 2 + w 2 ) + B Re 1 U Δ v , v + 1 U ( v z + w y ) , v + B Re 1 U ( w y y + w z z ) + 2 U w z , w .
From (16), we have
1 U Δ v , v = 1 U ( v y y + v z z ) , v = = 1 U v y , v y + lim z ¯ z S 1 U v z v z = 0 z ¯ 1 U v z , v 1 U v z , v z = Ω v y 2 + v z 2 U d Ω Ω 1 U z v 2 2 d Ω = Ω v y 2 + v z 2 U d Ω lim z ¯ z S 1 U v 2 2 z = 0 z ¯ + Ω 1 U v 2 2 d Ω .
We observe that as z ¯ z S we have
1 U v ( z ¯ z S ) α 1 , and 1 U v 2 ( z ¯ z S ) 2 ( α 1 ) .
Hence, the boundary terms vanish. Proceeding in the same way for the other terms, we easily obtain:
1 2 d d t v 2 + w 2 = 1 Re ( v 2 + w 2 ) B Re Ω 1 U ( v y 2 + w y 2 + v z 2 + w z 2 ) d Ω + Ω 1 U w 2 d Ω .
Since
1 U = z s 2 2 ( z z S ) > 0 and 1 U = z S 2 ( z z S ) 3 > 0 in 0 < z < z S ,
Poincaré inequality allows us to prove that v and w decay exponentially to 0 as t + . By using (13), it easy to see that also u goes to 0 as t + . Therefore, the energy norm of the kinetic perturbation tends to 0 for any Reynolds and Bingham number. We conclude that the streamwise perturbations are always linearly stabilizing.

4.2. Stability of Spanwise Perturbations

To investigate the spanwise perturbations we begin observing that, in this case, y 0 . It follows that the solenoidality of the kinetic perturbation reads u x + w z = 0 . From (14), integrating by parts, and using the boundary conditions as before, we have:
d d t v 2 2 = Re 1 v 2 + Re 1 B 1 U Δ v , v + 1 U v z , v = Re 1 v 2 Re 1 B Ω 1 U ( v x 2 + v z 2 ) d Ω .
Applying the Poincaré inequality and observing that 1 / U = z s 2 / ( 2 ( z z S ) ) > 0 for z < z S , we have that v goes exponentially to 0 as t + . The energy equations of u and w become:
d d t u 2 2 = ( w U , u ) ( p x , u ) 1 Re u 2 2 B Re Ω 1 U u x 2 d Ω , d d t w 2 2 = ( p z , w ) 1 Re w 2 2 B Re Ω 1 U w z 2 d Ω .
Adding the two equations together we have
d d t u 2 2 + w 2 2 = ( w U , u ) 1 Re ( u 2 + w 2 ) 2 B Re Ω 1 U ( u x 2 + w z 2 ) d Ω .
From this equation, we infer that the critical Reynolds number for spanwise perturbations is greater than the critical Reynolds number for Newtonian fluid (cf. [2]). Indeed, the term multiplied by the Bingham number is always non-positive, hence that term is always stabilizing.
The linear energy stability can be studied as in [1,2] using the classical energy
V ( t ) = 1 2 [ u 2 + v 2 + w 2 ] .
In order to obtain sufficient stability condition, we use the Orr–Reynolds energy identity:
V ˙ = ( w U , u ) 1 Re u 2 2 B Re Ω 1 U ( u x 2 + w z 2 ) d Ω B Re Ω 1 U ( v x 2 + v z 2 ) d Ω = I 1 Re D ,
where V ˙ is the orbital time derivative (i.e., the Lagrangian derivative computed along the solutions of (11), where now y 0 ), and
I = ( w U , u ) and D = u 2 + 2 B Ω 1 U ( u x 2 + w z 2 ) d Ω + B Ω 1 U ( v x 2 + v z 2 ) d Ω .
We have
V ˙ = I 1 Re D = I D 1 Re D 1 R ¯ 1 Re D ,
where
1 R ¯ = max S I D
and S is the space of the kinematically admissible fields
S = { u , v , w H 1 ( Ω ) , u = v = w = 0 on z = 0 and z = z S , u x + v y + w z = 0 , D ( u , v , w ) > 0 } .
In the definition of S , the space H 1 ( Ω ) is the Sobolev space of all functions which are in L 2 ( Ω ) together with their first generalized derivatives. The Euler–Lagrange equation of this maximum problem is given by
R ¯ ( U ( z ) w i + U ( z ) u k ) 2 Δ u 4 B 1 U u x x i + 1 U w z z k + 1 U w z k = λ ,
where λ is a Lagrange multiplier.
Since, for spanwise perturbations, v 0 and y 0 , by taking the third component of the double-curl of (17) and by using the solenoidality condition u x + w z = 0 , we obtain the generalised Orr equation [5]
R ¯ 2 U ( z ) w x + 2 U ( z ) w x z + Δ Δ w + 4 B 1 U w x x z z + 1 U w x x z = 0 ,
with boundary conditions w = w z on z = 0 and z = z s . We solve this boundary value problem with the Chebyshev collocation method, and we obtain the critical energy-linear Orr Reynolds number, R c .
In Figure 3 we have plotted the critical Orr-Reynolds number as a function of the Bingham number. The graphic shows the stabilizing effect of yield stress (Bingham). We note that in the limiting case as B 0 we obtain the critical value for the Newtonian fluid (note that here the depth of the layer is D while in [2] it is 2 D , therefore in the limiting case we obtain a critical value which is twice the size of the corresponding critical value of the Orr–Reynolds number obtained with the same boundary conditions).

5. Conclusions

In this article, we study the linear stability of laminar Bingham–Poiseuille flows for an inclined sheet of fluid (open channel) with constant slope angle β . The basic flow is parabolic in the shear region (the fluid there behaves like a Newtonian fluid) and is constant in the plug region. In the horizontal case, this problem has been investigated by many authors. In this work we follow the approach of Frigaard, Nouar and co-workers and we perform a novel investigation of the linear stability of the laminar solution which takes into account the inclination of the channel and a novel investigation of the critical value for energy stability, which possibly may suggest an experimental instability threshold.
The time evolution of linear perturbations of the basic solution is investigated both with spectral methods and with energy Lyapunov methods and assuming that the more destabilizing perturbations are two-dimensional (Squire-like Theorem). By using a Chebyshev collocation method, we prove that both the streamwise and the spanwise perturbations are stable for any Reynolds and Bingham numbers. We also show that the Bingham number has a stabilizing effect. This stability result, in the inclined channel, is similar to that of the Newtonian case [2] for the classical Couette motion in the horizontal layer.
The experiments show turbulence when the Reynolds numbers are sufficiently large (cf. also [19]), this turbulence cannot be accounted for by the spectral methods. The energy methods, instead, indicate the possible insurgence of instability for large enough Reynolds numbers (Couette paradox). In [1] we have already approached this question by introducing the energy method for a particular set of tilted perturbations, with tilt angles that appear in the experiments, and we have found a good agreement with the experimental data. We also apply to this problem the same ideas, and we obtain that the streamwise perturbations are always energy stable, while the spanwise perturbations are energy stable if the Reynolds number is less than a critical number which depends on the Bingham number. This suggests that there could be oblique perturbations that may physically destabilize the basic motion at the onset of instability.
The energy methods are also very important in the investigation of the nonlinear equations, since the nonlinear terms could be crucial in developing instabilities. For this reason we think that it is important, in a future work, to study the nonlinear stability problem and look for the most physical destabilizing perturbations and their possible relation with landslides.
Finally, we note that if we had introduced a Reynolds number independent of the inclination β , we would have obtained stability results with a different Reynolds number, inversely proportional to sin β . This means that small inclinations, as expected, are stabilizing.

Author Contributions

Conceptualization, P.F., A.G. and G.M.; methodology, P.F., A.G. and G.M.; software, P.F. and A.G.; formal analysis, P.F., A.G. and G.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the national project PRIN of the Italian Ministry for University and Research, grant number 2017YBKNCE, by the University of Catania, grant number PTRDMI-53722122113, and by PON SCN 00451 CLARA—CLoud plAtform and smart underground imaging for natural Risk Assessment, Smart Cities and Communities and Social Innovation.

Acknowledgments

We also thank the group GNFM of INdAM for financial support.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Falsaperla, P.; Giacobbe, A.; Mulone, G. Nonlinear stability results for plane Couette and Poiseuille flows. Phys. Rev. E 2019, 100, 013113. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Falsaperla, P.; Giacobbe, A.; Mulone, G. Stability of laminar flows in an inclined open channel. Ricerche Mat. 2020. [Google Scholar] [CrossRef]
  3. Moffatt, K. Fixed points of turbulent dynamical systems and suppression of nonlinearity. In Whither Turbulence; Lumley, J.L., Ed.; Lecture Notes in Physics; Springer: Berlin/Heidelberg, Germany, 1990; Volume 357, pp. 250–257. [Google Scholar]
  4. Reynolds, O. An experimental investigation of the circumstances which determine whether the motion of water shall be direct or sinuous, and of the law of resistance in parallel channels. Proc. R. Soc. Lond. 1883, 35, 84–99. [Google Scholar]
  5. Orr, W.M.F. The stability or instability of the steady motions of a perfect liquid and of a viscous liquid. Proc. R. Irish Acad. A 1907, 27, 69–138. [Google Scholar]
  6. Joseph, D.D. Eigenvalue bounds for the Orr-Sommerfeld equation. J. Fluid Mech. 1966, 33, 617–621. [Google Scholar] [CrossRef]
  7. Joseph, D.D.; Carmi, S. Stability of Poiseuille flow in pipes, annuli and channels. Quart. Appl. Math. 1969, 26, 575–591. [Google Scholar] [CrossRef] [Green Version]
  8. De Blasio, F.V. Introduction to the physics of landslides. In Lecture Notes on the Dynamics of Mass Wasting; Springer: Berlin/Heidelberg, Germany, 2001. [Google Scholar] [CrossRef]
  9. Liu, K.F.; Mei, C.C. Slow spreading of a sheet of Bingham fluid on an inclined plane. J. Fluid Mech. 1989, 207, 505–529. [Google Scholar] [CrossRef]
  10. Fusi, L.; Farina, A. Flow of a Bingham fluid in a non inclined channel. J. Non-Newton. Fluid Mech. 2016, 238, 24–32. [Google Scholar] [CrossRef]
  11. Frigaard, I.A.; Howinson, S.D.; Sobey, I.J. On the stability of Poiseuille flow of a Bingham fluid. J. Fluid Mech. 1994, 263, 133–150. [Google Scholar] [CrossRef]
  12. Nouar, C.; Frigaard, I. Nonlinear stability of Poiseuille flow of a Bingham fluid: Theoretical results and comparison with phenomenological criteria. J. Non-Newton. Fluid Mech. 2001, 100, 127–149. [Google Scholar] [CrossRef]
  13. Frigaard, I.; Nouar, C. On three-dimensional linear stability of Poiseuille flow of Bingham fluids. Phys. Fluids 2003, 15, 2843–2851. [Google Scholar] [CrossRef]
  14. Nouar, C.; Kabouya, N.; Dusek, J.; Mamou, M. Modal and non-modal linear stability of the plane Bingham–Poiseuille flow. J. Fluid Mech. 2007, 577, 211–239. [Google Scholar] [CrossRef]
  15. Métivier, C.; Frigaard, I.A.; Nouar, C. Nonlinear stability of the Bingham Rayleigh–Bénard Poiseuille flow. J. Non-Newton. Fluid Mech. 2009, 158, 127–131. [Google Scholar] [CrossRef]
  16. Nouar, C.; Bottaro, A. Stability of the flow of a Bingham fluid in a channel: Eigenvalue sensitivity, minimal defects and scaling laws of transition. J. Fluid Mech. 2010, 642, 349–372. [Google Scholar] [CrossRef] [Green Version]
  17. Allouche, M.H.; Millet, S.; Botton, V.; Henry, D.; Ben Hadid, H.; Rousset, F. Stability of a flow down an incline with respect to two-dimensional and three-dimensional disturbances for Newtonian and non-Newtonian fluids. Phys. Rev. E 2015, 92, 063010. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Squire, H.B. On the stability of three-dimensional disturbances of viscous flow between parallel walls. Proc. R. Soc. A 1933, 142, 621–628. [Google Scholar]
  19. Güzel, B.; Burghelea, T.; Frigaard, I.A.; Martinez, D.M. Observation of laminar–turbulent transition of a yield stress fluid in Hagen–Poiseuille flow. J. Fluid Mech. 2009, 627, 97–128. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Inclined layer of an angle β . The light grey layer is the shear layer which is moving with a parabolic law. The wall corresponding to z = 0 is the rigid terrain. The superior layer, the plug region D S z D , is rigidly moving at a constant velocity U ( D S ) . D is the depth of the layer, D S is the depth of the shear layer and D D S is the depth of the plug region. g is the gravity.
Figure 1. Inclined layer of an angle β . The light grey layer is the shear layer which is moving with a parabolic law. The wall corresponding to z = 0 is the rigid terrain. The superior layer, the plug region D S z D , is rigidly moving at a constant velocity U ( D S ) . D is the depth of the layer, D S is the depth of the shear layer and D D S is the depth of the plug region. g is the gravity.
Fluids 05 00141 g001
Figure 2. Real part of the time decay coefficient, Re ( C ) , when a runs from 0 to 2.6 , and the Reynolds number Re is in the interval [ 1000 , 11000 ] . The meshed plane corresponds to Re ( C ) = 0 , the surfaces corresponds to the set of points ( a , Re , Re ( C ) ) , for sample Bingham values B = 0 (purple), B = 2 (green) and, B = 5 (cyan).
Figure 2. Real part of the time decay coefficient, Re ( C ) , when a runs from 0 to 2.6 , and the Reynolds number Re is in the interval [ 1000 , 11000 ] . The meshed plane corresponds to Re ( C ) = 0 , the surfaces corresponds to the set of points ( a , Re , Re ( C ) ) , for sample Bingham values B = 0 (purple), B = 2 (green) and, B = 5 (cyan).
Fluids 05 00141 g002
Figure 3. Critical Orr–Reynolds numbers versus Bingham number. The plot shows the stabilizing effect of the Bingham number.
Figure 3. Critical Orr–Reynolds numbers versus Bingham number. The plot shows the stabilizing effect of the Bingham number.
Fluids 05 00141 g003

Share and Cite

MDPI and ACS Style

Falsaperla, P.; Giacobbe, A.; Mulone, G. Stability of the Plane Bingham–Poiseuille Flow in an Inclined Channel. Fluids 2020, 5, 141. https://doi.org/10.3390/fluids5030141

AMA Style

Falsaperla P, Giacobbe A, Mulone G. Stability of the Plane Bingham–Poiseuille Flow in an Inclined Channel. Fluids. 2020; 5(3):141. https://doi.org/10.3390/fluids5030141

Chicago/Turabian Style

Falsaperla, Paolo, Andrea Giacobbe, and Giuseppe Mulone. 2020. "Stability of the Plane Bingham–Poiseuille Flow in an Inclined Channel" Fluids 5, no. 3: 141. https://doi.org/10.3390/fluids5030141

APA Style

Falsaperla, P., Giacobbe, A., & Mulone, G. (2020). Stability of the Plane Bingham–Poiseuille Flow in an Inclined Channel. Fluids, 5(3), 141. https://doi.org/10.3390/fluids5030141

Article Metrics

Back to TopTop