Loading [Contrib]/a11y/accessibility-menu.js

Application of green-naghdi equations for a nonlinear stability of a tangential-velocity discontinuity in shallow water flow.

Research Article | DOI: https://doi.org/10.31579/ 2690-0440/002

Application of green-naghdi equations for a nonlinear stability of a tangential-velocity discontinuity in shallow water flow.

  • Thi Thai Le 1

Nuremberg Campus of Technology (NCT), Nuremberg Institute of Technology Georg Simon Ohm 246b Further Strasse, 90429 Nuremberg, Germany

*Corresponding Author: Thi Thai Le, Nuremberg Campus of Technology (NCT), Nuremberg Institute of Technology Georg Simon Ohm 246b Further Strasse, 90429 Nuremberg Germany.

Citation: Thi Thai Le (2021) Application of Green-Naghdi equations for a nonlinear stability of a tangential-velocity discontinuity in shallow water flow. J. Mathematical Methods in Engineering, 2(1): Doi: 10.31579/ 2690-0440/002

Copyright: © 2021 Thi Thai Le. This is an open-access article distributed under the terms of The Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Received: 16 April 2021 | Accepted: 30 May 2021 | Published: 18 June 2021

Keywords: kelvin-helmholtz instability; tangential-velocity discontinuity; green-naghdi equations; nonlinear analyse; shallow water flow.

Abstract

It is well known that the Kelvin-Helmholtz instability (KHI) is always unstable in the absence of wave, regardless of the strength of velocity difference. But in the present of sound wave, the KHI is sup-pressed if the velocity difference is equal or greater than √8 times to the velocity of sound wave. Since a shallow water flow has an analogy with a compressible gas flow in 2D, a linear stability of the KHI in a shallow water flow was considered by Bedenzkov and Pogutse to obtain the same condition √8 of the Froude number which defines by the ratio of the velocity difference to the traveling velocity of gravity wave. In this paper, we consider the nonlinear stability of the KHI in a shallow water flow by using the Green-Naghdi equations. We obtain the critical value of the ratio of the Froude number to make the inter-face stabilized is smaller than √8 given by others. This critical value is analysed by using both of the analytical method and the Sturm’s theory.

Introduction

The discontinuity of a tangential velocity in an incompressible fluid is well known as the Kelvin Helmholtz instability. The flow is necessarily unstable, The discontinuity of a tangential velocity in an incompressible fluid is well known as the Kelvin Helmholtz instability. The flow is necessarily unstable, regardless of the velocity difference [1, 2]. However, a surprising results was shown by Landau [3] that the compressibility can suppress the KHI. The flow is stable for a large velocity difference. A mathematical analogy of sound waves in a compressible fluid and gravity waves in a shallow water flow was also mentioned (see p. 322 [3]). However, the flow of a compressible fluid is stable for only two-dimensional (2D) flow but always unstable for three-dimensional flow [4]. Cairns [5] showed that the physical mechanism for instability lies in the fact that the wave on the interface has a negative energy. The unstable region is produced by a coalescence of positive and negative energy modes. Miles [6] and Ribner [7] were the first study of over-reflexion problem related to the transmission and reflexion of a sound waves at a vortex sheet separating by two regions of constant horizontal velocity U1 and U2. This problem was extended by Fejer [8] by including the effects of hydromagnetic, then McKenzie [9] included the effects due to buoyancy. The similar stability of tangential-velocity discontinuity in a shallow water in 2D was given by Bedenzkov and Pogutse [10] since a shallow water flow has an analogy with a compressible gas flow in 2D. The horizontal length scale is assumed much greater than the vertical length scale. The analogy of stability theory of compressible fluids with that of shallow water is limited to two dimensions because the hydrostatic balance is employed in the vertical direction and that we have to consider perturbations, with wavelength     

, where H is the depth of the fluid layer, depending only on the coordinates of the horizontal plane of the liquid layer (not on the depth coordinate z) [11, 12]. The interface between two regions of fluid is stabilised√8 if the velocity difference U is equal or greater than √8 times to the velocity of gravity waves c = √gH [13, 14], with g being the gravity acceleration and H being the depth of water. However, Bedenzkov and Pogutse considered the linear stability of discontinuity surface by ignoring the higher-order of dissipative gravity waves.

In this paper, we consider the nonlinear stability of a tangential-velocity discontinuity by using the Green-Naghdi equations (GN). The GN describe the wave propagation of fully nonlinear and weakly dispersive gravity waves on the fluid of finite depth. The GN are derived by using shallow water scaling asymptotics for a domain of small aspect ratio, however there is no restriction placed on the Froude number [15, 16, 17]. The consideration of nonlinear wave plays a vital role in predicting the shape of wave such as tsunami and tides [18]. We obtain the dispersion relation of wave frequency and other characteristics by enforcing boundary conditions at the interface of tangential-velocity discontinuity. The stability characteristic of an interface is reduced by solving the dispersion equation. Five roots of complex frequency ω are gained as a functions of the velocity discontinuity U, the traveling speed c and the order of qH with H being the depth of water and q being the wavenumber in the direction x of wave propagation. The resulting dispersion relation is calculated analytically and compared with that used the Sturm´s theory to find the number of real roots of a polynomial equation [19, 20]. Sec. 2 shows the formulation of the problem and the dispersion equation. In Sec. 3, we go into the stability characteristic by both solving the dispersion equation analytically and using the Sturm´s theory. A summary and conclusions is presented in the last section (Sec. 4).

Formulations and Dispersion Relation
We consider two-region flow which is moving in the region y > 0 and is at rest for y < 0, as shown in Figure 1. We consider only the discontinuity of tangential velocity in the direction x of wave propagation, the normal
velocity in the y direction is still continuous for entirely region. We consider small perturbation (proportional to ei(qx−ωt)) over the surface of discontinuity, suppose that
u(x,y,z,t) = U0+u˜(x,y,z,t), v(x,y,z,t) = v˜(x,y,z,t), h(x,y,t) = H+u˜(x,y,t),              (1)
Where

                                                                                          

               

By ignoring the higher order items, the irrotational Green-Naghdi equations are written as

In which, D/Dt = ∂/∂t + U0∂/∂x.

Figure 1: The top view of the interface of the tangential-velocity discontinuity in a shallow water flow. In region I (y > 0), the fluid is moving with uniform velocity U but is at rest in region II (y < 0). The surface of discontinuity in tangential velocity is horizontally perturbed to y = ζ(x,t) with an infinitesimal amplitude. By taking derivative of two last equations on x, y respectively and then substituting into the first one, we obtain easily

 

                       (4)

We seek solution in the form ei(qx−ωt) eκy, where κ = −κ1 for y > 0 and κ = κ2 for y < 0. In region I (y > 0) with U0 = U, therefore we have the relation between wavenumber κ1 and other quantities as follows

 

Similarly, in region II (y < 0) with U0 = 0, we have

The kinematic boundary condition is

               (7)

Where ζ(x,t) = aei(qx−ωt) is the position of the discontinuity surface, with a is small constant.
The normal component of the velocity are equal on both sides of the interface and equal to the movement of the surface in that direction. Thus, this gives


v˜1 = v˜2 at y = ζ(x,t) ≈ 0.                     (8)

The pressure p should be continuous at the tangential discontinuity. In the hydrostatic approximation p = −ρgh˜, the condition of pressure is reduced to that of the wave depth on both sides of the interface:

h˜1 = h˜2 at y = ζ(x,t) ≈ 0.         (9)

The horizontal displacement ζ of the interface is connected with the vertical one h˜, thus we obtain

                   (10)

Combination of (5), (6) and (10), we reduce the dispersion relation between the wave frequency ω and other quantities as follows

The first factor gives one real root which    

   does not contribute to the instability of the interface. Thus, we do not consider this root longer. The instability now is considered simply from the contribution of the second factor as

                  ( 12)

Stability Analyze of an Interface
Analytical Solution
The interface is stabilized if and only if the dispersion equation (12) has only real roots, since it is a quartic polynomial of with real coefficients of argument ω. We introduce new variable as follows 

                                                                            

                                                        (13)

The dispersion equation (12) then turns to

Equation (14) gives four analytical solutions as follow

               (15)

 

The stabilized condition reduces to

     

                                                       (16)

 

                                                                               (17)

 Therefore, the interface is stabilised if and only if    

. This critical value of Froude number to make the interface stabilised is smaller than the one √8 given by Landau [3] and the others [4, 10].
Figure 2 displays the imaginary and real part roots given in (15) for a given qH = 1. The imaginary part of roots Ω±, − tend to zero at the Froude number M = √24/5 ≈ 2.19089 < 2.822 ≈ √8 and then vanish for greater values of the Froude number. Two other roots Ω±, + are always real for any positive Froude number, i. e. the imaginary parts Im[Ω±,+] = 0. Therefore, the flow is stable if and only if the Froude number M is equal or greater than √24/5; otherwise unstable. 

Where M = U/c defines the Froude number, then we have

Froude number M


Figure 2: The imaginary and real parts of root Ω in equation (15) for a given qH = 1. Since two roots Ω±,+ are always real, their imaginary parts are zero.

By using the Sturm’s Theorem
Next, we go to find the condition of the dispersion equation (12) to have all real roots by using the Sturm’s Theorem (see appendix B) since the flow is stable if and only if (12) has only real roots. The Sturm’s sequences are constructed as following:

,

                  p2(ω) = −remainder(p0,p1),      (18) p3(ω) = −remainder(p1,p2),
                    p4(ω) = −remainder(

,

 

Here remainder (pj, pj+1) is the remainder of the devision of polynomial pj to pj+1 for j = 1 ,..., 3. All roots are real if and only if all highest coefficients of pi (ω), i = 0,...,4 must have same signs. After simple calculation, we find the highest coefficients of p0 (ω), p1 (ω), p2 (ω), p3 (ω) are always positive. Thus, we just need p4 (ω) positive, that is

                        (19)

In which M is the Froude number defined in the section 3.1. This result coincides with that given in (17). In other words, we found that the critical value of the Froude number M to make the interface stabilised is √8(1 +

)    which is smaller than that one √8 given by others.

Conclusion
We have considered the nonlinear stability of the Kelvin-Helmholtz instability in a shallow water flow by using the Green-Naghdi equations. The stability characteristic of the interface obtains by analysing the dispersion relation of wave frequency and others. The analytical solution is compared with that used the Sturm’s theory. We obtain a coincidence of two these methods and show that the flow is stable with an amount of velocity discontinuity smaller than that given by previous researches [3, 4, 10].

The stability condition of flow depends on both the strength of velocity difference and the depth of considered region. For a given velocity difference, the flow is more stable in deeper regions and is less stable in shallower regions. This work was done at Nuremberg Campus of Technology of Technische Hochschule Nu¨rnberg Georg Simon Ohm. The author would like to thank to Prof. Frank Ebinger and Department for supporting this research.

A The irrotational Green-Naghdi equations
The Green Naghdi equations (GN) is the fully nonlinear shallow-water
waves whose amplitude is not necessarily small and represents a higher-order correction to the classical shallow-water equations. In shallow-water approximation, the velocity field of three dimensional long wavelength propagation can approximate over the depth as follows [17]

                        (20)

Then we can derive the irrotational Green-Naghdi equations [17]

           

        (21)

Where ∇ = (∂/∂x, ∂/∂y) is horizontal gradient; γ = (∇·u¯)2−∇·u¯t−u¯·∇[∇·u¯] is the vertical acceleration at the free surface and can be derived as

                                   (22)

B Sturm’s theorem
Here we introduce the Sturm’s theorem [19, 20] which expresses the number of distinct real roots of a univariate polynomial p located in an interval in terms of the number of changes of signs of the values of the Sturm’s sequence at the bounds of the interval. Applied to the interval of all the real numbers, it gives the total number of real roots of p. The Sturm sequence is a finite sequence of polynomials, applying Euclid’s algorithm to p and its derivative:

p2(x)=−remainder(p0,p1),
p3(x)=−remainder(p1,p2),
: :
pm=−remainder(pm−2,pm−1),

0=−remainder (pm−1, pm).

where remainder(pj,pj+1) is the remainder of the polynomial long division of pj by pj+1, and where m is the minimal number of polynomial divisions (never greater than deg(p)) needed to obtain a zero remainder.
Let σ(ξ) denote the number of sign changes (ignoring zeroes) in the sequence p0(ξ),p1(ξ),p2(ξ),··· ,pm(ξ).
Sturm’s theorem then states that for two real numbers a < b, the number of distinct real roots of p in the interval [a,b] is σ(a) − σ(b).
Choosing a = −∞, b = ∞, then the total number of real roots of a polynomial is equal to σ (−∞) –σ (∞). That is to say, all roots of a polynomial of degree m are real, if and only if σ (−∞) – σ (∞) = m.

As 0 ≤ σ (−∞) ≤ m, 0 ≤ σ (∞) ≤ m, so σ (−∞) = m,σ(∞) = 0.
Since the sign of a polynomial is decided by the term with highest degree as x → ±∞. Thus all highest coefficients of the polynomials in the Sturm sequence must have the same signs (all positive or all negative).
In conclusion, we have the following Lemma: Suppose p(x) is a univariate polynomial with real coefficients, m is the highest degree of p(x). p0 (ξ), p1 (ξ), p2 (ξ),··· ,pm(ξ) is the Sturm sequence of p(x). If all highest coefficients of pi(x),i = 0,1,··· ,m have the same signs (all positive or all negative), then all roots of the polynomial p(x) are real.

References

a