-□ □-
Розглядаеться задача про поширення внутрштх хвиль для идеальной нестисливог ридини. Гiдродинамiчна система складаеться з трьох шарiв сктченног товщини, як не 3Mi-шуеться та обмежен зверху твердою кришкою, а знизу твердим дном. На поверхнях роздту ридких середовищ дiе сила поверхневого натягу. Постановка задачi була здшс-нена в безрозмiрному виглядi. В якостi малого параметру використовуеться коефщент нелiнiйностi, який дорiв-нюе видношенню характерног амплтуди до характерног довжини хвилi.
Розв'язки лшшног задачi шукаються виглядi прогресив-них хвиль. На основi цих розв'язтв отримано дисперсшне спiввiдношення як умову розв'язуванностi системи лтш-них диференщальних рiвнянь. Виявлено инування двох характерних мод (дшсних коретв дисперсшного рiвнян-ня). Проаналiзовано графши коретв дисперсшного рiв-няння в залежностi вид рiзних фiзичних та геометричних параметрiв системи. Встановлено, що товщини шарiв не впливають на дисперсю хвиль, тодi як змта поверхневого натягу та вгдношення густин значно впливають на умови проходження хвиль. Розглянуто хвильовi пакети в лтшнш постановщ, що е суперпозищею гармончних хвиль близьког довжини. Виявлено, що амплтуда обвгдног хви-льового пакету на нижнш поверхш контакту залишаеть-ся синусогдальною, тодi як на верхнш змшюеться за бшьш складним законом.
Задача поширення внутрштх хвиль вздовж поверхш трьох ргдких шарiв можуть моделювати сильно страти-фжований термоклин в океат. Вивчення впливу поверхневого натягу також може бути використане при розроб)щ нових технологш, пов'язаних з використанням трьохргд-ких шарiв, що не змшуються
Ключовi слова: идеальна нестислива ргдина, внутрш-н хвилi, тришарова гiдродинамiчна система, дисперсшне рiвняння
-□ □-
UDC 532.59
|DOI: 10.15587/1729-4061.2018.1326911
ANALYSIS OF CONDITIONS FOR THE PROPAGATION OF INTERNAL WAVES IN A THREE-LAYER FINITE-DEPTH LIQUID
Yu. Hurtovyi
PhD, Associate Professor* Е-mail: hurtovyy@gmail.com V. Naradovyi PhD, Senior Lecturer* Е-mail: naradvova1986@gmail.com V. Bohdanov Postgraduate student* Е-mail: vovka14@bigmir.net * Department of Applied Mathematics, Statistics and Economics Volodymyr Vynnychenko Central Ukrainian State Pedagogical University Shevchenka str., 1, Kropyvnytskyi, Ukraine, 25006
1. Introduction
Fluids with discrete or continuous density stratification are encountered in numerous applied problems. Perturbation of such an inhomogeneous fluid causes propagation of internal waves with interesting characteristic properties. As a rule, mathematical formulas of such problems consist of a system of nonlinear differential partial derivative equations. Nonlinearity of such systems brings about significant mathematical difficulties in the problem solution. Therefore, some small parameter that permits rejection of the equation members having insignificant effect and the problem statement linearization itself is singled out in the problem statement. The so-called dispersion relation which binds the frequency and the wave number is the condition for solvability of such a system. In a general case, the dispersion relation has several real roots which are called modes. In linear problems, solution is considered as a superposition of various modes. Thus, the dispersion relation indicates conditions of waves propagation for a concrete hydrodynamic system. The conditions of internal wave propagation are quite sensitive in relation to physical and geometric parameters of the system.
Relevance of this study direction was determined by the necessity of developing a theoretical apparatus for the problems of detecting internal waves in various media. The problems of internal wave propagation along the surfaces of contact of three liquid layers can model strongly depth-stratified thermocline regions in the ocean. In shallow seas, vertical density stratification has a three-layered structure with a prominent seasonal pycnocline at a depth of about several tens of meters and a major pycnocline at a greater depth. For example, the Black Sea and the Baltic Sea have a more or less permanent three-layer structure formed by flow of fresh water on the surface and penetration of the saltiest water to the bottom layers. A noticeable three-layer density stratification is also found in the South China Sea. However, many important issues remain unexplored. In a three-layer fluid, specific classes of nonlinear waves, so called envelope solitons, can spread. They are still poorly studied both analytically and numerically. All this shows urgency of the problem of studying the internal waves in a three-layer fluid. The knowledge of the effect of surface tension can also be applied in the development of new technologies associated with the use of three non-mixing liquid layers.
©
2. Literature review and problem statement
A new type of stability of wave solutions related to the models of Korteweg-de-Vries type was studied in [1] under condition of a low dispersion of waves. Existence of exponentially increasing solutions of the linearized problem with the use of asymptotic perturbation methods was established. The criteria of spectral instability of single waves in the linear case were found. Besides, nonlinear stability and spectral instability of solutions for a specific set of parameters were analyzed. However, this paper was devoted to the study of mathematical aspects of wave propagation with no relation to concrete models.
A study of propagation of a single internal wave down a sloping bottom in a three-layer hydrodynamic system based on the Korteweg-de-Vries equation with variable coefficients was presented in [2]. Precise and numerical solutions of the problem were considered and compared. Two model configurations for which the system energy change was numerically analyzed were considered. Existence of the phenomenon of fading waves was established. This model is based on the Korteweg-de-Vries equation which enables study of classical solitons but does not take into account the potential appearance of envelope solitons described by the nonlinear Schrodinger equation.
Spread of internal waves over obstacles was studied in [3]. It was established for narrow obstacles that the energy of the wave flow increases in an upstream flow. The situation is just opposite for wide ridges. It was shown that the pressure field makes a significant contribution to the total energy flow. A linear theory was developed that is valid for general stratification and the surface flow that accurately predicts the wave field. This paper has described internal waves in the presence of a flow and obstacles but did not take into account formation of internal waves under the influence of gravitation and capillary forces.
Overturning regular waves were studied in [4] using a numerical model of fluid dynamics based on the Navi-er-Stokes equations. The constructed model well predicts vortices and overlapping of waves but does not describe the phenomena of propagation and interaction of internal waves.
Dynamics of collision of two internal waves in a two-layer fluid in a model based on the Navier-Stokes equation was investigated in [5]. It was shown that collision of two waves of small and medium amplitude leads to a phase shift and generation of a dispersive wave. It was established that the maximum amplitude of a wave formed after collision can exceed the cumulative amplitude of the collided waves. Wave interaction can also be studied by other methods. One of such methods is the multi-scale method which leads to evolutionary envelope equations.
Work [6] was devoted to the study of single Rosby waves in a stratified fluid. A new model was constructed with the use of the generalized Boussinesq equation. Based on this equation, a generalized Euler-Lagrange equation was obtained with the use of the Agraval methods and reduction of perturbations. Exact solutions were found for the latter. To analyze the dissipation phenomenon, approximate solutions were constructed in the model with the use of a new iterative method. The waves that have a considerable length as regards thickness of the layers were considered in this paper but the presence of waves with lengths commensurate with the geometric parameters of the system was not considered.
Exact solutions were constructed in [7] in the form of a progressive wave for the generalized connected nonlinear Korteweg-de-Vries equations based on an extended algebraic method. The results obtained describe generation and evolution of such waves as well as their interaction and stability. Korteweg-de-Vries equations of higher orders were considered in this paper. Solution of these equations gives soliton waves. Conditions for propagation of such waves were not taken into account
Solutions in the form of singular waves of the fifth-order Corteweg-de-Vries equation were studied in [8]. In this equation, in addition to the quadratic nonlinearity and dispersion of the third order, there is a cubic nonlinearity and a linear dispersion of the fifth order as well as two nonlinear dispersion members. An exact solution was obtained and the dependence of its amplitude and velocity on the parameters of the equation was studied. It was shown that the resulting solution can be a soliton. It was indicated that the presence of nonlinear dispersion members may materially affect the single wave existence, profile and resistance to small perturbations. The results obtained are applicable to the study of internal and surface waves in liquids as well as waves in other media. A single-layer model with a free surface and propagation of internal waves at the interface of two media were considered in this paper. However, the presence of two or more contact surfaces required construction and study of new models.
Interaction of small internal wave packets with a large internal wave was studied in [9]. The main result was that the indicated interaction with a single wave leads to a complete destruction of the packet which consists of waves having small length in comparison with the length of single waves. In the future, a study of interaction of wave packets in hydrodynamic systems of a more complex structure would be of interest.
Capillary-gravitational and interphase waves in a two-layer fluid were studied in [10] with the use of the Stokes wave theory. The study was conducted by the method of perturbations of the third order. It was established that solutions of the third order depend essentially on the surface tension, density and thickness of the layers. In a linear approximation, the constructed model coincides with classical results. Similar studies are also important for three-layer liquids in which a significant dependence of wave amplitudes on various system parameters were also revealed.
Propagation of waves in a two-layer liquid of finite depth with a free surface was studied in [11]. Internal and surface waves were considered with the use of numerical simulation by the finite element method. The ratio of amplitudes of two waves for various values of density and thickness of liquids was analyzed. It should be noted that the analysis of the amplitude ratio is facilitated when the solutions are obtained in an analytical form.
Internal weakly nonlinear waves of finite amplitude in a three-layer liquid of a finite depth were studied in [12]. A pair of Korteweg-de-Vries equations with coefficients depending on the liquid parameters was obtained. As a result of introduction of an additional small parameter, two Gardner equations for two contacting surfaces were obtained. Soliton solutions were constructed and studied for the obtained equations depending on various parameters of the system. Cases are considered when Gardner equations can be written without nonlinear terms. This model does not take into account propagation of wave packets with their inertia described by nonlinear Schr dinger equations.
Work [13] is devoted to the study of single waves in a three-layer fluid with a constant density of each layer. A completely nonlinear numerical scheme based on integral equations was presented. Based on the constructed scheme, various types of waves were studied. The phenomenon of wave overturning was established when the upper layer thickness is much larger than that of the lower layers. This study was devoted to the numerical analysis of wave propagation but analytical solutions describing wave motions in a three-layer liquid are more suitable for revealing new effects.
Linear waves and nonlinear wave interactions in a three-layer fluid were studied in [14] with the use of analytical and numerical methods. It was established that single waves cannot be described mathematically as solitons because of insignificant energy loss. However, since these losses are small and the waves in such a system are disper-sionally connected, the single waves manifest themselves as solitons. In this case, an issue of studying another type of wave motion in three-layer systems appear, in particular, propagation of waves whose lengths are commensurate with the geometric parameters of the system.
A model of a two-layer incompressible ideal fluid was studied in [15]. Solution of the original nonlinear problem was constructed in the form of power series for a small parameter. For the coefficients of these series, there are corresponding linear (homogeneous and inhomogeneous) problems solution of which is represented in the form of series of previously constructed eigenfunctions of a nonlinear boundary value problem. Another approach with the use of power series for a small parameter is a multi-scale method which enables obtaining of information about wave motions for various time and spatial scales.
Propagation of internal separated waves in a two-layer liquid medium bounded above and below with solid surfaces without taking into account surface tension was studied in [16]. Solution was found in the form of generalized power series for a parameter depending on the magnitude inverse to the Froude number. It is interesting to study influence of surface tension on conditions of wave propagation in a three-layer system.
The attention deserving theoretical studies in [17] where the model equations derived from the Euler equations describing evolution of internal gravitational waves in a two-layer liquid that are non-viscous and non-mixed have been solved without imposing any limitation on the amplitude. Also, areas of suitability of asymptotic approximations for strongly-nonlinear wave packets of internal waves in the layer - layer system were indicated in this paper A similar study of application of the multi-scale method is relevant for three-layer systems as well.
The problems of propagation and stability of waves in a two-layer hydrodynamic system were considered in [18, 19]. The issue of modulation stability were considered and the system parameters at which wave packets will be modulation-unstable were established. Evolution of wave packets is described in these papers by means of a nonlinear Schroding-er equation using the multi-scale method. Similar studies are also important for a three-layer liquid.
A model of propagation and interaction of waves in a two-layered hydrodynamic system with a free surface was constructed in [20-22]. In particular, conditions for the propagation of waves at various density ratios in the hydro-dynamic system were revealed in the first approximation. Dependencies of potential frequencies on the lower layer
thickness and the wave number were also analyzed. Interaction of internal and surface waves was studied. Influence of presence of surface tension on the contact surfaces was revealed. Study of the dispersion modes makes it possible to determine structure of the wave motions on the contact surfaces which is also relevant for a three-layer system of a finite depth.
A problem statement for an ideal incompressible three-layer liquid with a solid bottom and a solid cover was considered in this study. At the same time, the characteristic scales of physical quantities were introduced which enabled the mathematical equations to be presented in a dimension-less form and the small parameter (coefficient of nonlinear-ity) be highlighted. Solution of a linear system leads to a fairness of the relation which is a dispersion relation. The dependence of the modes of the dispersion relation on physical and geometric parameters of the system was analyzed and the wave motion on the contact surfaces was studied.
3. The aim and objectives of the study
This study objective was to analyze propagation and interaction of waves along the contact surfaces in a three-layer hydrodynamic "a layer with solid bottom - a layer - a layer with a cover" system. This will enable a more detailed assessment of the wave processes and interaction in three-layer systems.
To achieve the objective, the following tasks were formulated:
- to perform mathematical statement of the problem and linearize it;
- to find solutions of the linearized problem and derive the dispersion relation;
- to construct dispersion diagrams and analyze the roots of the dispersion relation depending on the physical and geometric parameters of the system being studied;
- to construct wave packets on the contact surfaces and study their shape for various modes of the dispersion relation.
4. The statement and solution of the problem of propagation of waves in a three-layer hydrodynamic system
4. 1. The mathematical problem statement and the method of study
The problem of propagation of three-dimensional wave packets of finite amplitude on the surface of the liquid layer was studied.
Q1 = {(x,z): |x| <~, -h < z<0} with density pi, of the medium liquid layer
Q2 = {(x,z): |x| <~, 0<z<h,} with density p2, and the upper liquid layer
Q3 = {(x,z): |x| <~, h2 <z<h2 + h3}
with density p3. Velocities of the liquids, uj, in Qj are expressed in terms of gradient of the potential qj, j=1, 2, 3.
The layers Q and Q2 are separated by the contact surface z=ni(x, t) and the layers Q2 and Q3 are separated by the contact surface z=h2+n2(x, t). When solving, the forces of surface tension on the contact surfaces are taken into consideration. Gravitational force is directed perpendicular to the interface in the negative z-direction. Liquids are considered to be ideal and incompressible and the wave motions under study are vortex-free and potential (Fig. 1):
boundary conditions (conditions of impenetrability)
dî3 = 0 at 2 = h2 + h3, dz 2 3
= 0 at z = -h.
dz
Introduce dimensionless quantities in (1) using the following formulas:
x = Lx', z = Hz',
^1,2 = ani^ t =
4gH
t ',
Fig. 1. The problem statement
gaL ,
îi,2,3
4gH
The mathematical problem statement in a plane case has the following form: Laplace equation
d>, d V -
A9j = + = 0, uj = V9j in Qj, j = 1,2,3,
1 dx2 dz2 i^i j
kinematic conditions on contact surfaces
3n1 j_ dn1 j
dt dz dx dx
dt dz dx dx dynamic conditions on contact surfaces
at z = n,( x, t ), j = 1,2, at z = h2 + n2(x,t), j = 2,3,
(1)
T1,2 L P1 gT1 ,2, P1,2,3 p1p1,2,3 ,
where stroked variables are dimensional, and H, L, a are characteristic thickness of the layer, length, and wave amplitude, respectively, p1 is the density of the lower layer and g is the acceleration of gravity. In a dimensionless form, the statement (1) will take the following form (the strokes are omitted below):
; _ i av
= 0
dx2 p dz2 in Q., j = 1,2,3,
dn -1 dî.=-adni dîz.
dt P dz dx dx
Pi ^ "P2 ^ + g (Pi P2 )ni + 0.5pi
-0.5p2
dîi dx
dîi dz
\ - T [i+i3n J21
) Ux J
fdîi 12+fdîi 1 ^ dx ) I dz J
-3/2
2
dx2
= 0
at
z = ni( x, t ),
dî2 dî3
P2^ _P3^ + g n2(P2 -P3)-
+h2(p2 -P3) +0.5P
-0.5P
dîi dx
-T
dîa
dx
i + |dn2
dx
-3/2
dîa
dz
dx2
2
= 0
dî2
dz
2
at z = ani( x, t ), j = i,2,
dn. -1 dîz.=-adni dî.
dt P dz dx dx
at
z = h2 + an2( x, t ), j = 2,3,
dî dî2
Pi^dt^ "P2 +(Pi-P2)ni"
+0.5Pia -0.5P2a
dîi J2+i fdîi
dx ) p^dz
dî2 J2+i fdî!
dx ) pt dz
2
2
-T
i+1«^
3/2
dx2
= 0
(2)
at
at
z = h2 + n2( x, t ),
z = ani( x, t ),
+ ail2(p2-Pa) +
dt 3 dt +h2(p2 -pa) + 0.5p2a
dîi 12+1 fdîi dx J p t 92
2
-0.5paa2
-T2a
12 +1 fd% dx J p t 92
-, N2H-3/2
dn
2
1+1«^
dx2
= 0
at
2 = h2 + an2( x, t ),
boundary conditions (conditions of impenetrability)
^^ = 0 at 2 = h2 + ha, d2
^ = 0 at 2 = -h, d2
In the statement (2), there are two parameters: a = —
H2 L
(the coefficient of nonlinearity), and P = (the coefficient
of dispersion). Nonlinear problem (2) will be studied using the method of multi-scale developments. Let us assume for the future that the coefficient of nonlinearity a<< 1, and P = 1. Then the unknown potentials of velocity and deviations of the contact surfaces can be represented as:
n (x, t ) = £a"-1n„ (x0, xi, x 2, t0, ti, t2) + 0(aa), i =1,2,
n=1
a
( x, 2, t ) = (x0, x1, x 2, t0, t1, t2) + 0(aa),
j = 1,2,3,
(3)
where xk=akx and tk=akt (k=0, 1, 2) are scale variables.
Substitute (3) into (2) and equate expressions with the same degrees a to obtain linear approximations of the problem under study.
The problem of the first approximation (at a0)
j, xo + =0 in " j , j = 1,2,3
n11,t0 + j = 0 at 2 = 0, j = 1,2,
n21,t0 +j = 0 at 2 = h-2, j = 2,3,
P1^11t0 -P2^21,t0 +(p1 -p2)^11 - 71^11x0 X0 =0 at 2 = 0,
P2^21,t0 -Pa^31,t0 +(P2 -P3)n21 -72n21,x0x0 =0 at 2 = K
9312 = 0 at 2 = h2 + h3,
= 0 at 2 = -h1. (4)
The problem of the second approximation (at a1)
jx0x0 +jz = -2^j1,x0x, in , j = 1,2,3
n12,t0 j 2,2 = -n11A -n11,x0 9;1,x0 +n119 j1,22 at 2 = 0 j = 1,2>
n22,t0 -9j2,2 = -n21,t1 -n21,x0 9;1,x0 +n219;1,22 at 2 = h2, j = 2,3,
p1%2,t0 -p2^22,t0 +(p1 -p2)n12 -71n12,x0x0 =
= -p1(911t1 +n11911t0 2 ) +p2(921A +
+n11921t02 ) - °.5p1(9l21,x0 +921,2 )" +0.5p2(92 +^21,2 ) +2T1n
11, x0 x1
at 2 = 0,
p2^22,t0 -pa932t0 +(p2 -pa)n22 -T2n22,x0x0 = = -p2(921A +n21921t02 ) +p3(931,t1 + +n21931t02 ) - °.5p2(921,x0 +9212 ) + +0.5p3(921,x0 +9^1,2 ) +2T2n21,x0x1
at 2 = h2,
9322 = 0 at 2 = h2 + h3, 9122 = 0 at 2 = -h1.
The problem of the third approximation (at a2)
9j3,x0x0 +9j322 = -9j1,x1 x1 - 29j1,x0x1 - 29j1x0x2
in Q., j = 1,2,3,
j
n13t0 —9j3,2 = n11,t2 —n12,t1 —n11,x0 9j1x1 —
-n11,x19j1,x0 -n12,x0 9j1,x0 -n11,x0 9;2,x0 +
+ n119 j2,22 +n129j1,22 + 0.5n219 j1,222 ^nHx,, 9 j1,x0 2
at 2 = 0, j = 1,2,
n23,t0 -9j3,2 =-n21t2 -n22A -n21,x0 9j1,x1 -n21,x19j1,x0 -n22,x0 9j1,x0 -n21,x0 9j2,x0 + +n219j2,22 + n229j1,22 + °.5n219j122 - n21n21,x0 9j1,x02
at 2 = h2, j = 2,3,
p1913t0 -p2923t0 +(p1 -p2)n13 - T1n13,x0x0 =
= -p1(911fe +912,t1 +n12911,t0 2 +
+^1911,^2 + °.5ni22911,t0 22 ) +
+p2(921,t2 + 922,t1 +n12921,t0 2 +n11921,t12 +
+°.5ni22921,t0 22) -p1(n11912t0 2 + 911,x0 911,x1 +
+911,x0 912,x1 +n11911,x0 911,x0 2 +
+911,2911,2 + ^1911*91122) + p2(n11922,t0 2 +
+921x0 921,x1 +921,x0 922x +
+n11921,x0 921x02 + 921,2921,2 + n11921,2921,22 ) +
+T1n11,x1 x1 + 2T1n12,x0x1 + 2T1n12,x0x2 -
-1.5T1(n11,x0)2 n11,x0 x0
(5)
(6)
2
at z = 0,
P2Î23,t0 -P3Î33,t0 +(P2 -P3)n23 - x0 =
= -P2(Î2U2 +Î22t +n22Î2it0z + +n2iÎ2i,ii + 0.5n^2Î2i,i0zz ) +
+P3(Î3U2 +Î32ti +n22Î3i,ioz +n2iÎ3i,iiz + ^n^U.zz )" -P2(n2iî22,i0z +Î2i,x0Î2i,xi +Î2i,x0Î22,xi +
+n2iÎ2i,x0Î2W + Î2i,zÎ2i,z + n2iÎ2i,zÎ2i,J +
c = sinh2 (kh2 )((i - P2 )k + Tik3 )((P2 - P3 )k + T2k3 ).
According to (8), there are two pairs of independent solutions to the problem (3): - for the o>i mode:
n(i = A exp(i( kx -w?t));
n? =
sinh(kh2 ) ((i - P2 )k + Tik3 - w?2 coth(khi ) - P2 w?2 coth(kh, ))
+P3 (n2iî32,t0z + î3i,x0 î3i,x1 + î3i,x0 î32,x1 + n2i93i,*0 93i,x0z + xA exp(i(kx -w?t ));
2
w? P2
+Î31,zÎ31,z +n2lÎ31,zÎ31,zz) + ^21,*,*, +
+2T2n22,x0x1 + 2T2n22,x0x2 - ^(^l^ )*n21,x0x0 at z = h2,
î33z = 0 at z = h2 + hg, î13z = 0 at z = -h?.
The problems (4) to (6) are linear with respect to the functions that are summands in the decompositions (3). Since all terms in the decompositions (3), except for the first, enter with coefficients a in the first degree and above, therefore the solutions of the first linear approximation (4) give the main contribution to solution of the problem (2).
4. 2. Solutions to the problem of the first approximation and the dispersion relation
In a solution of the first approximation (4), a dispersion relation is obtained which simultaneously is a condition of solvability of the problem (4) in the form:
sinh2(kh2)
(1 -P2)k + T?k3 --w2 (coth(kh? ) + P2 coth(kh2 )) (P2 -P3)k + T2k3 + +w2 (coth(kh, ) - P3 coth(kh ))
= 0.
(7)
2a
2a
î(1) = -
î21 =
iw
k sinh(kh?)
cosh(k(h? + z))A exp(i(kx - w?t));
(9)
iw
-cosh(k(h2 - z))-
k sinh( kh2 )
i ((1 - P2 )k + T?k3 - w?2 coth(khj ) - P2w?2 coth(kh, )) kP2w?
cosh(kz)
x A exp(i(kx - w?t ));
(1) i sinh(kh2 )((1 - P2 )k + T?k3 - w?2 coth(kh? ) - P2w?2 coth(kh, )) î31 kP2w? sinh(kh3 )
xcosh(k(h2 + h3 - z))Aexp(i(kx- w?t));
- for the ro2 mode:
n(2) =
w22P2
sin h(kh2)((1 -P2 )k + T?k3 - w22 coth(khj ) - P2w22 coth(kh2 )) x B exp(i(kx -w2t ));
n(2?) = B exp(i( kx -w?t));
îS?)
iw23P2
k sinh2(kh?)((1 - P2 )k + T?k3 - w22 coth(kh? ) - P2w22 coth(kh2 )) xB cosh(k(h? + z))exp(i(kx - w2t)); (10)
î22i)=-
iw
k sinh(kh2)
The dispersion relation (7) can be written as a biquadratic equation relative to ro which has two x pairs of roots:
w2 P2
2 -b Wb2 - 4ac 2 -b-Vb2 - 4ac w42 =- and w2 =---, (8)
sin h(kh2)((1 -P2 )k + T?k3 - w22 coth(khj ) - P2w22 coth(kh2 )) xcosh(k(h2 - z)) + cosh(kz) x B exp(i(kx - w2t ))
where a, b, c are the coefficients at ro4, ro2 and a free term in the corresponding biquadratic equation
a = p2 -sinh2(kh2)(coth(kh?) + P2 coth(kh2)) x x (coth(kh2 ) - P3 coth(kh )),
î32i) =
k sinh^^)
cosh(k(h2 + h3 - z))Aexp(i(kx-w?t)).
b = sinh2(kh2)
((1 - P2 )k + T?k3 )(coth(kh2 ) - P3 coth(kh3 )) -- ((P2 - P3 )k + T2k3 ) (coth(khj ) + P2 coth(kh2 ))
Solutions (9) and (10) are independent while n2? in (9) is a response wave on the upper contact surface to a wave riff, which has amplitude A and propagates with a frequency ro1 on the lower contact surface. In (10), r|(f is a response wave on the bottom contact surface to a wave n(2i\ which has amplitude B and propagates with a frequency ro2 on the upper surface of the contact.
Analysis of the roots of the dispersion relation will be made hereinafter depending on various geometric and physical parameters of the system under study. The issue of wave motion on the contact surfaces is also considered.
5. Analysis of the dispersion relation
5. 1. Dispersion diagrams for various physical and geometric parameters of the system
Let us study dependence of the dispersion relation (8) modes, ro1 and ro2, on the change of the ratio of layer densities. Fix p3 = 0.7 and change p2 in the range from 0.7 to 1 in steps of 0.05 at the following values of the system parameters:
value from 0.7 to 0.9 leads to a gradual decrease in the ro2 value (Fig. 3, b).
Next, let us analyze the effect of the surface tension on the change of the ro1 and ro2 modes. Fix the value of T2 = 0 and change the value of T1 in a range from 0 to 0.2 in 0.05 steps (Fig. 4). The change of the surface tension, T1, on the bottom contact surface nearly does not affect the change of the ro1 mode (Fig. 4, a). At the same time, taking into account surface tension on the lower contact surface significantly affects variation of the ro2 mode for the wave numbers k > 0.5, corresponding to the gravitational-capillary and capillary waves (Fig. 4, b). A similar pattern was also observed with consideration of the surface tension on the upper surface of the contact, T2 (Fig. 5).
hi = h = h3 = 1, T = T2 = 0.
Cases p2 = 0.7 and p2 = 1 will not be considered because they are marginal and will be considered below.
a b
Fig. 2. Dependence of the dispersion relation modes on the wave number for various values P2 = {0.75,0.8,0.85,0.9,0.95}: №1 (a); (b)
If density of the middle layer, p2, is the mean value of the densities of the lower and upper layers, then the roi graph is higher than the graphs for other P2 values (Fig. 2, a) and the ro2 graph is lower than other graphs (Fig. 2, b) (for wave numbers k > 2). If the value of p2 becomes closer to p1 or p3, then the ro1 graph descends and the ro2 graph ascends.
0.6 o.s
0.4 CD 0.3 0.2 0.1 0
3 4
if
3 4 k
a b
Fig. 3. Dependence of the dispersion relation modes on the wave number for various values: ro1 (a); ro2 (b)
Fix the value of p2 = 0.9 and change p3 in a range from 0.7 to 0.9 in 0.05 steps (Fig. 3). At the fixed value of p2, the change of p3 from 0.7 to 0.8 has little effect on the value (Fig. 3, a). Further approximation of the p3 value to the p2 value leads to the ro1 zeroing. This will be described when considering the marginal cases. Change of the p3
-0-0.05-0.1 -0.15
-0-0.05-0.1 -0.15
a b
Fig. 4. Dependence of the dispersion relation modes on the wave number for various values: ro1 (a); ro2 (b)
b
Fig. 5. Dependence of the dispersion relation modes on the wave number for various values: coi (a); co2 (b)
a b
Fig. 6. Dependence of the dispersion relation modes, ro1 and ro2, on the wave number and: thickness of the lower layer, h1 (a); thickness of the upper layer, h3 (b)
Analyzing the effect of geometric parameters of the system, i. e. thickness of the lower layer, h1 and thickness of the upper layer, h 3, the study has shown that the change of these parameters has no significant effect on
a
the change of the and ro2 modes for fixed wave numbers (Fig. 6).
5. 2. Marginal cases of the system under study
Let us consider the marginal cases in which the studied three-layer "a layer with a solid bottom - a layer - a layer with a solid cover" system degenerates into a two-layer "a layer with a solid bottom - a layer with a solid cover" system. The first marginal case arises under conditions P2 = P3 and T2 = 0. The second marginal case arises under the conditions P? = P2 and T = 0. In both cases, the first pair of roots, ro^ of the dispersion relation becomes equal to zero and the second pair of roots, w 2, degenerates into a pair of roots of the dispersion relation of the "a layer with a solid bottom - a layer with a solid cover" system.
(1 -P)k + Tk3
coth(kh? ) + P coth(kh2 )
[18].
5. 3. Analysis of wave motion on the contact surfaces
According to (9) and (10), the wave motion on the lower and upper contact surfaces can be represented in the first approximation as:
nil = n(?i+n(i), n2i = n2?+n22i),
(ii)
that is, as a linear combination of waves for various modes. Let us consider propagation of wave packets on the lower and upper contact surfaces which are a linear combination of waves of the same amplitude and a close length:
n? = nii(ki) + nii(k2> + nii(k3X n2 =n21(kl) +n21(k2> +n21(k3>.
(12)
Fig. 7 shows appearance of wave packets (12) on the contact surfaces for the following system parameters:
h = h2 = h3 = 1, T1 = T2 = 0, p2 = 0.95,
p3 = 0.92, A = 0.05, B = 0.025.
Fig. 7. Wave packets on the contact surfaces at
k1 = 0.009, k2 = 0.01, k3 = 0.011
Fig. 8 shows the case of propagation of the wave packet on the lower and upper contact surfaces which is a linear combination of waves with frequency ro1.
Fig. 9 shows the case of propagation of a wave packet which is a linear combination of waves with a ro2 frequency on the lower and upper contact surfaces.
n i
Fig. 8. Propagation of wave packets that correspond to the mode in a three-layer system
Fig. 9. Propagation of wave packets in a three-layer system corresponding to the mode
In the case of waves of a similar length corresponding to the mode passing on the lower and upper contact surfaces (Fig. 8), the wave packets formed on both contact surfaces have a clear shape, i. e. amplitude of the envelope of such wave packets is sinusoidal. A similar situation was observed in passing of a group of waves of similar length corresponding to the ro2 mode (Fig. 9). In the presence of groups of waves of a similar length, for both and ro2 modes on the upper and lower contact surfaces, the wave packets formed on the lower contact surface have a clearer shape than on the upper contact surface (Fig. 7). That is, amplitude of the envelope on the lower contact surface remains sinusoidal whereas it changes according to a more complicated law on the upper contact surface.
6. Discussion of the results obtained in the study of a three-layer hydrodynamic system
The study of a three-layer system with solid cover and bottom was carried out for a model in which the wave amplitude is much smaller than the wave length (the coefficient of
nonlinearity a = -0- is much smaller than one) and the wave
lengths are close to the thickness of the middle layer (the dis-
persion coefficient P = ~jT is close to one). It is interesting to
study other cases, for example, when P is also close to zero. As a result of application of the multi-scale method to the third order, it was possible to obtain the first three approximations of the problem under consideration. Let us derive dispersion relation and solutions of the first approximation constructed on the basis of its roots.
w2 =
During the study of the dispersion relation roots, it was found that the ro1 and ro2 modes differed in their sensitivity to changes of the system physical parameters. In particular, the ro1 mode was more sensitive to the change of the ratio of the p2 densities than the ro2 mode while the ro2 mode was more sensitive to the change of p3, T1 and T2 parameters. The study has also shown that both modes were insensitive to the change of geometric parameters of the system, namely thicknesses of the lower, h1, and upper, h3, layers.
It was also found that w2 > w1 for the same wave numbers. A question arises about the contribution of waves with ro1 and ro2 frequencies in their motion on the upper and lower contact surfaces, i. e. at which system parameters the waves with above frequencies make a significant contribution to the wave motion on each contact surface. To answer this question, it is necessary to study the amplitude ratio of the waves on the lower and upper contact surfaces. It is the subject of further studies.
In the course of study of the dispersion relation, two marginal cases have been established in which the three-layer system under study degenerated into a two-layer system. In this case, the ro1 mode became equal to zero and the dispersion relation degenerated to the corresponding dispersion relation for the "a layer with a solid bottom - a layer with a solid cover" system [18].
The structure of wave motion on the contact surfaces during passage of groups of waves of similar length was analyzed. It was found that the amplitude of the envelope of the wave packets on the lower contact surface was sinusoidal in the first approximation and it changed on the upper surface according to a more complicated law. A detailed study of evolution of the wave-packet envelopes requires an analysis of higher approximations (5) and (6) which should result in the study of nonlinear Schrodinger equations of the second order.
It was planned to continue the study of a three-layer hydrodynamic system for obtaining approximations of the second and third order using the multi-scale method.
7. Conclusions
1. The mathematical statement of the problem of propagation of internal wave packets in a three-layer hydrody-namic "a layer with a solid bottom - a layer - a layer with a free surface" system was fulfilled. A transition to dimension-less quantities was made during which two parameters were found: nonlinearity coefficient and dispersion coefficient. On the basis of this statement, a linear approximation of the problem under study was obtained.
2. Solutions of the linear problem in the form of progressive waves were obtained. A dispersion relation was derived based on these solutions. Existence of two characteristic modes (the dispersion relation roots) was established.
3. The dispersion relation roots of the on the diagrams (w,k) were analyzed depending on various physical and geometric parameters of the system. It was found that the change of geometric parameters of the system (thickness of the lower and the upper layers) did not affect the dispersion relation modes. The change of the ratio of densities affected both modes but w1 to a greater extent. The second mode, w2, was more sensitive to the change of parameters p3, T1 and T2.
4. The shape of wave packets on the contact surfaces for both modes which are roots of the dispersion relation was investigated. Wave packets in a linear formulation which are a superposition of harmonic waves of a similar length were considered. It was found that the amplitude of the wave packet envelope on the lower contact surface remained sinusoidal whereas it changed according to a more complicated law on the upper contact surface.
References
1. Pava J. A. Stability properties of solitary waves for fractional KdV and BBM equations // Nonlinearity. 2018. Vol. 31, Issue 3. P. 920-956. doi: 10.1088/1361-6544/aa99a2
2. Yuan C., Grimshaw R., Johnson E. The evolution of second mode internal solitary waves over variable topography // Journal of Fluid Mechanics. 2017. Vol. 836. P. 238-259. doi: 10.1017/jfm.2017.812
3. Lamb K. G., Dunphy M. Internal wave generation by tidal flow over a two-dimensional ridge: energy flux asymmetries induced by a steady surface trapped current // Journal of Fluid Mechanics. 2017. Vol. 836. P. 192-221. doi: 10.1017/jfm.2017.800
4. Water wave overwash of a step / Skene D. M., Bennetts L. G., Wright M., Meylan M. H., Maki K. J. // Journal of Fluid Mechanics. 2018. Vol. 839. P. 293-312. doi: 10.1017/jfm.2017.857
5. Frontal collision of internal solitary waves of first mode / Terletska K., Jung K. T., Maderich V., Kim K. O. // Wave Motion. 2018. Vol. 77. P. 229-242. doi: 10.1016/j.wavemoti.2017.12.006
6. Lu C., Fu C., Yang H. Time-fractional generalized Boussinesq equation for Rossby solitary waves with dissipation effect in stratified fluid and conservation laws as well as exact solutions // Applied Mathematics and Computation. 2018. Vol. 327. P. 104-116. doi: 10.1016/j.amc.2018.01.018
7. Seadawy A. Stability Analysis of Traveling Wave Solutions for Generalized Coupled Nonlinear KdV Equations // Applied Mathematics & Information Sciences. 2016. Vol. 10, Issue 1. P. 209-214. doi: 10.18576/amis/100120
8. Khusnutdinova K. R., Stepanyants Y. A., Tranter M. R. Soliton solutions to the fifth-order Korteweg-de Vries equation and their applications to surface and internal water waves // Physics of Fluids. 2018. Vol. 30, Issue 2. P. 022104. doi: 10.1063/ 1.5009965
9. Xu C., Stastna M. On the interaction of short linear internal waves with internal solitary waves // Nonlinear Processes in Geophysics. 2018. Vol. 25, Issue 1. P. 1-17. doi: 10.5194/npg-25-1-2018
10. Third-order stokes wave solutions of the free surface capillary-gravity wave and the interfacial internal wave / Meng R., Cui J., Chen X., Zhang B., Zhang H. // China Ocean Engineering. 2017. Vol. 31, Issue 6. P. 781-787. doi: 10.1007/s13344-017-0089-z
11. Min E.-H., Koo W. Hydrodynamic characteristics of internal waves induced by a heaving body in a two-layer fluid // Ocean Engineering. 2017. Vol. 145. P. 290-303. doi: 10.1016/j.oceaneng.2017.09.017
12. Propagation regimes of interfacial solitary waves in a three-layer fluid / Kurkina O. E., Kurkin A. A., Rouvinskaya E. A., Soomere T. // Nonlinear Processes in Geophysics. 2015. Vol. 22, Issue 2. P. 117-132. doi: 10.5194/npg-22-117-2015
13. Rusas P.-O., Grue J. Solitary waves and conjugate flows in a three-layer fluid // European Journal of Mechanics - B/Fluids. 2002. Vol. 21, Issue 2. P. 185-206. doi: 10.1016/s0997-7546(01)01163-3
14. Weidman P. D., Nitsche M., Howard L. Linear Waves and Nonlinear Wave Interactions in a Bounded Three-Layer Fluid System // Studies in Applied Mathematics. 2011. Vol. 128, Issue 4. P. 385-406. doi: 10.1111/j.1467-9590.2011.00540.x
15. Peregudin S. I. Vnutrennie i poverhnostnye volny v sloisto-neodnorodnoy zhidkosti // Mater. Mezhdunar. konf. "Differenc. uravneniya i ih pril.". Saransk, 1995. P. 269-276.
16. Bontozoglou V. Weakly nonlinear Kelvin-Helmholtz waves between fluids of finite depth // International Journal of Multiphase Flow. 1991. Vol. 17, Issue 4. P. 509-518. doi: 10.1016/0301-9322(91)90046-6
17. Choi W., Camassa R. Fully nonlinear internal waves in a two-fluid system // Journal of Fluid Mechanics. 1999. Vol. 396. P. 1-36. doi: 10.1017/s0022112099005820
18. Selezov I. T., Avramenko O. V., Gurtoviy Yu. V. Osobennosti raspostraneniya volnovyh paketov v dvuhsloynoy zhidkosti konechnoy glubiny // Prykladna hidromekhanika. 2005. Vol. 7, Issue 1. P. 80-89.
19. Selezov I. T., Avramenko O. V., Gurtoviy Yu. V. Ustoychivost' volnovyh paketov v dvuhsloynoy gidrodinamicheskoy sisteme // Prykladna hidromekhanika. 2006. Vol. 8, Issue 4. P. 60-65.
20. Nonlinear interaction of internal and surface gravity waves in a two-layer fluid with free surface / Selezov I. T., Avramenko O. V., Gurtovyi Y. V., Naradovyi V. V. // Journal of Mathematical Sciences. 2010. Vol. 168, Issue 4. P. 590-602. doi: 10.1007/s10958-010-0010-2
21. Avramenko O. V., Naradovyi V. V., Selezov I. T. Conditions of Wave Propagation in a Two-Layer Liquid with Free Surface // Journal of Mathematical Sciences. 2015. Vol. 212, Issue 2. P. 131-141. doi: 10.1007/s10958-015-2654-4
22. Avramenko O. V., Naradoviy V. V. Analysis of propagation of weakly nonlinear waves in a two-layer fluid with free surface // Eastern-European Journal of Enterprise Technologies. 2015. Vol. 4, Issue 7 (76). P. 39-44. doi: 10.15587/1729-4061.2015.48282