CC BY
13
RUSSIAN JOURNAL OF EARTH SCIENCES, VOL. 20, ES4006, doi:10.2205/2020ES000734, 2020
SPECIAL ISSUE
Dedicated to Efim Pelinovsky on the occasion of his 75th Birthday
July 2020 • Moscow, Russia
Long internal ring waves in a two-layer fluid with an upper-layer current
Karima Khusnutdinova1
Received 10 June 2020; accepted 14 June 2020; published 12 July 2020.
We consider a two-layer fluid with a depth-dependent upper-layer current (e.g. a river inflow, an exchange flow in a strait, or a wind-generated current). In the rigid-lid approximation, we find the necessary singular solution of the nonlinear first-order ordinary differential equation responsible for the adjustment of the speed of the long interfacial ring wave in different directions in terms of the hypergeometric function. This allows us to obtain an analytical description of the wavefronts and vertical structure of the ring waves for a large family of the current profiles and to illustrate their dependence on the density jump and the type and the strength of the current. In the limiting case of a constant upper-layer current we obtain a 2D ring waves' analogue of the long-wave instability criterion for plane interfacial waves. On physical level, the presence of instability for a sufficiently strong current manifests itself already in the stable regime in the squeezing of the wavefront of the interfacial ring wave in the direction of the current. We show that similar phenomenon can also take place for other, depth-dependent currents in the family. KEYWORDS: Stratified shear flows; internal waves; long-wave instability.
Citation: Khusnutdinova, Karima (2020), Long internal ring waves in a two-layer fluid with an upper-layer current, Russ. J. Earth. Set., 20, ES4006, doi:10.2205/2020ES000734.
1. Intr°ducti°n water waves in Cartesian [Ablowitz and Segur, 1979],
cylindrical [Johnson, 1980], and elliptic-cylindrical Long-wave models of the Korteweg-de Vries (KdV) [Khusnutdinova et al., 2013] geometries and intertype have found numerous useful applications in nal waves, with a possible background shear flow, the studies of the oceanic nonlinear surface and in- in Cartesian geometry (see [Grimshaw, 2015] and ternal waves (see, for example, [Grimshaw, 1998; references therein), as well as cylindrical Korteweg-Helfrieh and Melville, 2006; Grimshaw et al., 2007; de Vries (cKdV)-type models for surface waves Grimshaw et al., 2010] and references therein). without shear flow [Miles, 1978; Johnson, 1980] Two-dimensional generalizations of these models and on a shear flow [Johnson, 1990] and inter-developed in the context of fluids include versions nal waves without shear flow [Lipovskii, 1985] and of the Kadomtsev-Petviashvili (KP) equation for on a shear flow [Khusnutdinova and Zhang, 2016a;
Khusnutdinova and Zhang, 2016b]. The latter mod-■W?!'", oTLMltKhlm!icalnS!ien1c!f;„Loughbor- els are relevant to the description of the nearly annular internal waves generated in straits, river-sea
i
ough University, Loughborough, Great Britain
Copyright 2020 by the Geophysical Center RAS. http://rjes.wdcb.ru/doi/2020ES000734-res.html
interaction zones as well as waves scattered by localized topographic features (see [Vlasenko et al.,
ES4006
1 of 13
2009; Nash and Moum, 2005; Vlasenko et al., 2013; Khusnutdinova and Zhang, 2016a; and references therein). Some relevant initial-value problems have been considered in [Weidman and Zakhem, 1978; Ramirez et al., 2002; McMilan and Sutherland, 2010; Khusnutdinova and Zhang, 2016b; Grimshaw, 2019].
Several recent studies have predicted various effects of shear flows on two-dimensional surface ring and ship waves (see [Johnson, 1990; Ellingsen, 2014a; Ellingsen, 2014b; Khusnutdinova and Zhang, 2016a; Li and Ellingsen, 2019] and references therein) and some of these predictions have been recently confirmed in laboratory experiments [Smeltzer et al., 2019]. The study of the effects of a piecewise-constant shear flow on long surface and interfacial ring waves in a two-layer fluid has shown a striking difference in the shapes of the wavefronts of these waves: while the surface ring waves were elongated in the direction of the current, the interfacial ring waves were squeezed in that direction [Khusnutdinova and Zhang, 2016a]. Overall, the study of the effects of the shear flow on internal ring waves is in its infancy. It the present paper we aim to build a rather large family of current profiles amenable to theoretical analysis which could be used to approximate some currents present in natural settings, and to elucidate their effects on the long internal ring waves.
subject to the free surface and rigid bottom boundary conditions:
w = ht + uhx + vhy at z = h(x,y,t), p = pa at 2 = h(x,y,t), w = 0 at 2 = 0.
Here, u, v, w are the velocity components in x, y, z directions respectively, p is the pressure, p is the density, g is the gravitational acceleration, z = h(x, y, t) is the free surface height (z = 0 at the bottom), and pa is the constant atmospheric pressure at the surface. We assume that in the basic state uo = uq(z), vo = wo = 0, poz = -pog, h = ho. Here uo(z) is a horizontal shear flow in the x-direction, and po = po(z) is a stable background density stratification. The vertical particle displacement ( is used as an additional dependent variable, which is defined by the equation
(t + u(x + v(y + w(z = w,
subject to the surface boundary condition
( = h — ho at 2 = h(x,y,t),
where ho is the unperturbed fluid depth.
The problem is considered using the following non-dimensional set of variables:
2. Modal Equations for Ring Waves
In this section we overview the derivation of the far-field set of modal equations for long ring waves in a stratified fluid over a parallel depth-dependent shear flow derived and studied in [Khus-nutdinova and Zhang, 2016a; Khusnutdinova and Zhang, 2016b]. A ring wave propagates in an in-viscid incompressible fluid, described by the set of Euler equations:
p(ut + uux + vuy + wuz) + px = 0, p(vt + UVX + VVy + wvz) + py = 0, p(wt + uwx + vwy + wwz) + pz + pg = 0, Pt + upx + vpy + wpz = 0, ux + vy + wz = 0,
x ^ Xx, y ^ Xy, z ^ ho z, t ^ — t,
„ * , * , h0C*
U ^ C U, V ^ C V, w ^ —— w,
X
(po,p) ^ p*(po,p), h ^ ho + ar], r h0
P ^ Pa + P*Po(s)g ds + p*ghop, Jz
where A is the wave length, a is the wave amplitude, c* = yjgh0 is the long-wave speed of surface waves, p* is the dimensional reference density of the fluid, while p0(z) is the non-dimensional function describing stratification in the basic state, and •q = rj(x, y, t) is the non-dimensional free surface perturbation. Non-dimensionalisation leads to the appearance of two small parameters in the problem, the amplitude parameter e = a/h0 and the wavelength parameter ô = h0/X. The maximal balance condition ô2 = e has been imposed in [Khusnutdinova and Zhang, 2016a].
The problem is then solved in the moving cylindrical coordinate frame (moving at a constant speed c: a natural choice is the speed of the shear flow at the bottom, as follows from the derivation). We consider deviations from the basic state and use the same notations u and v for the projections on the new coordinate axis, scaling the appropriate variables by the amplitude parameter e:
X ^ ct + r cos в, у ^ r sin в, z ^ z, и ^ uq(z) + e(u cos в — v sin в),
v ^ e(u sin в + v cos 6), w ^ ew, p ^ ep, p ^ pQ + ep.
t ^ t,
[Khusnutdinova and Zhang, 2016a]
u\ = —Аф-uoz cos в — v1 = Афиоz sin в —
kF
к2 + к'2 k'F л^
Аф,
Афх
к2 + к
W1 = Aç Рф,
Po
Pi
к2 + к'2 Pi = —poz Аф. щ = Аф at
AF 2фг
z = 1,
where the function (f) = <fi(z, d) satisfies the following set of modal equations:
The modal equations are obtained by looking for a solution of the problem in the form of asymptotic multiple-scales expansions of the form ( = (j + e(2+..., and similar expansions for other variables, where
PoF 2 k2+k'-
Фг) — p0 г Ф = 0,
F 2
к2+к'
2 фг — ф = 0 at z = 1, ф = 0 at z = 0, and F = —s + («о — c)(k cos в — к' sin (
(1) (2) (3)
z
Ci = А(£,В.,в)ф(г,в),
with the appropriate set of fast and slow variables:
i = rk(ß) — st, R = егк(в), в = в,
The speed of the moving coordinate frame с is fixed to be equal to the speed of the shear flow at the bottom, с = uo(0). Then, F = — s = 0 at z = 0, and the condition Fф = 0 at z = 0 which appears as a result of the derivation implies (3).
The amplitude function R, в) is then found by considering the equations at O(e). It satisfies a cylindrical Korteweg - de Vries (cKdV)-type equation [Khusnutdinova and Zhang, 2016a]
where we define to be the wave speed in the absence of a shear flow (with k(ff) = 1). When a shear flow is present the function k(d) is responsible for the adjustment of the wave speed in a particular direction, and is to be determined. The formal range of asymptotic validity of the model is defined by the conditions £ ~ R ~ 0(1). To leading order, the wavefront at any fixed moment of time is described by the equation k( ) = constant, and for the sake of definiteness we consider outward propagating ring waves, requiring that the function k = k(ff) > 0.
To leading order, assuming that perturbations of the basic state are caused only by the propagating wave, the motion is described by the solution
A Ae
piAR + P2AAÇ + рзА^ + pi — + p5 = 0,
where the coefficients pi,i = 1, 5 are given in terms of solutions of the modal equations (1) - (3) by the following formulae:
p1 = 2s / роРф2 dz,
P2 = —3 I роР2ф1 dz,
рз = — (к2 + к'2)
poF2ф2 dz,
1
1
о
1
о
^ = - l { № - 3k'2)F2 - 4k'(k2 + k'2){uo - c) sin OF
- sin2 0(uo - c)2(k2 + k'2)2) + (k'F + (k2 + k'2)(uo - c) sin 0) j dz,
2k (1
= -k2 + k,2 J PoF^2z[k'F + (uo - c)(k2 + k'2) sin6] dz.
In this paper we are concerned with the analysis of the modal equations for a two-layer fluid with the upper-layer current in the rigid-lid approximation, which will allow us to describe, to leading order, the wavefronts and vertical structure of the long interfacial ring waves. Our primary goal is to analyse the sensitivity of the shape of the wavefront to the variability of the background shear flow in the bulk of the fluid layer, but we will also illustrate the 3D modal functions and discuss the onset of the long-wave instability in the limiting case of a piecewise-constant current. In what follows, the free surface condition (2) is eventually replaced with the rigid-lid approximation
0 = 0 at 2 = 1.
01 =
A(k2 + k'2)
O-11 I)
h = — p1d
Pi
1 - (k2 + k'2) f %
FZ\
(4)
(5)
where F1 = -s + U(z)(k cos d - k' sin d) and A is an arbitrary constant. The derivative of this solution is discontinuous at z = d and must satisfy the jump condition:
F (d)2 k2 + k '2
[po4>] - [po]<P(d) = 0,
(6)
where F(d) = —s, yielding the following nonlinear first-order ordinary differential equation for the function k(d):
3. Two-layer Fluid with an Upper-Layer Current
We consider a two-layer fluid (see Figure 1. The density in the upper layer is p\, and the density in the lower payer is p2- The flow in the lower layer is assumed to have constant speed, and then without any loss of generality we can assume that this speed is equal to zero. The upper-layer flow is described by the function U(z). Thus,
Uo =
U(z), if d < z < 1, 0, if 0 < z<d.
d - oc j: D
t) (k2 + «2)2 +
(7)
s2 I " dz \ , 2 . i2\
(1 + 7 if) ^ + 2) - 7 = 0
We will assume the continuity of uo at 2 = d, i.e. U(d) = 0. This is a generalisation of the case studied in [Johnson, 1990; 1997] for surface waves in a homogeneous fluid, it is a possible model for a river inflow, an exchange flow in a strait, or a windgenerated current, for example.
The continuous solution of the modal equations (1)-(3) is given (in respective layers) by
Figure 1. Schematic of the problem formulation.
z
d
If there is no current, i.e. U(z) = 0, then F1 = —s,k = 1,k' = 0, and the equation (7) reduces to an algebraic equation for the speed :
s4 - s2 - Pi-2 d(1 - d)=0,
P2
implying
s 2 61,2
1 1 + 4 ^ d(ï-d)
2
(8)
giving us the speed of the surface and interfacial ring waves in the absence of any current.
In the rigid-lid approximation, the solution (4), (5) is replaced with
$1 = -
$2 = -
A(k2 + k'2) f1 dz_ Pi Jz F2, A(k2 + k'2)z f1 dz
1 d
F12
pid + p2S2 f* ^
k2 + k'2 =
(P2 -Pi)dJ1 f When U(z) = 0, this equation gives
s2 = ( P2 - Pi)d(1 - d)
U(z) = j(z - d)a,
Figure 2. Upper-layer currents U(z) = 7(z — d) (blue, solid), U(z) = 7(z — d)2 (red, dotted) and U(z) = 7yfz—~d (red, dot-dashed) for d = 0.7 and U (1) = 0.0045.
(9) a = 2,7 = 0.05 when a = 2 (all currents have the (10)
same strength U(z) = 0.0045 on the surface 2 = 1). The equation (11) takes the form
where A is an arbitrary constant, while the jump condition (6) yields
k2 + k!2 =
,2 _ Pid + P2S21[K(9)]
(P2 - Pi)dl[K(&)]■
(14)
where
(11)
(12)
n K (*)]=/ f = -2Fi(.
1 - d 1 1
2^ 2, -, 1 + -,
a
a
K( )
(1 - )
Pid + P2 (1 — d)
the speed of the interfacial ring wave in the absence of the background current, and in the rigid-lid approximation. The equation (11) and the speed (12) can be formally obtained from (7) and (8) as an approximation when p2 — pi ^ p 1, p2.
4. Wavefronts and Vertical Structure
Let us now consider the family of the upper-layer current profiles described by the function
and K(d) = k(d) cos 9 — k!(d) sinö. Here, the function I[K(0)] is given in terms of the hypergeometric function 2F\.
The general solution of the equation (14) has the form
k(ff) = a cos d + b(a) sin d,
where
and
b2 = -a2 +
2 Pid + P2S2 1(a)
(P2 - Pi)d 1(a)
(15)
(16)
(13)
where 7 and a are some positive constants. In particular, the upper-layer current is shown in Figure 2 for a = 1, a = 1 and a = 2 for d = 0.7 and 7 = 0.015 when a = 1, 7 = 0.00821584 when
-i
1(a) =
d
Jd [ s -7a(z - d)a\2 1 -d 2Fi( 2,1, 1 + 1 ^ (1 -df)
a
a
d
The far-field wavefronts of the ring waves at a fixed moment of time are described by the curves
H(r, d, t) = rk(ff) — st = constant,
where k(ff) is the singular solution of the equation (11) (i.e. the envelope of the general solution found by requiring ^ = 0). In what follows all wavefronts are plotted for rk(0) = 50 ( R ~ 1,e ~ 0.02).
In this general setting we assume that the current is sufficiently weak, so that there exists a part of the wavefront which is able to propagate in the upstream direction, and therefore d e [—ft, ft], but we will also discuss other possible regimes for the limiting case of a constant upper-layer current. It is sufficient to define the solution for d e [0, t] because of the symmetry of the problem. For the family of currents described by (13), the singular solution can be found explicitly in parametric form
some stronger currents [Khusnutdinova and Zhang, 2016a].
It is worth noting that in many natural cases the hypergeometric function featured in the solution reduces to elementary functions. In particular, for a = 1 (i.e. U(z) = 7(z — d)), we have
2F1 (2,1, 2,^ (1 -d))
1
1 - ^(1 -dy
for a = 1 (i.e. U(z) = 7\/2 - d),
2F1 (2, 2, 3,^ VT-d) =
2
1-d
/1-71
+
ln(1 -
^(1 - ^^r-d l2a2
and for a = 2 (i.e. U(z) = 7(z - d)2),
k(a) = a cos 9(a) + b(a) sin d(a),
=
I p1d + P2S2 1(a) (p2 - Pi)d 1(a)
- a2
(17)
(18)
2
2 ■ 2 -T (1 -*2) 1
2(1 - ^(1 -d)2)
+
2 ViT (1 -d)
arctanL I — (1 - d)].
0( a) =
arctan
arctan
2b(p2 - P1)d I2(a)
2a(p2 - p1 )d I2(a) + (1 - p2s21(a))I'(a)
_2b(p2 - P1)d I2(a)_
2a(p2 - p1 )d I2(a) + (1 - p2s21(a))I'(a)
if a e [ao,amax] + K if a e [an
e ^l]^
K (19) M (0 e [-,k]);
where
I'(a) =
1- d
aas
2
L[1 - ^(1 -d)»]2
- 2F1I 2
2
1, 1 + 1 (1 -dy-
a a
Here, the parameter a takes values in the interval [ amin ,amax] which is found by requiring b2 > 0 in (16). We consider the outward propagating ring waves, and therefore require k(ff) > 0 for all d. Then, the interval must contain a = 0 since a should take both positive and negative values to allow k(ff) to be positive at both 9 = 0 and d = t. The value ao e [amin, amax] is found from the condition
db(a) 1
In the first case, for U(z) = 7(z — d), the singular solution can be rewritten in the form k = k(d) as follows
m = \ 1 +
1
PUS
_2( P2 - P1)_
PUS
2(P2 - Pi)
cos ,
da
tan
0.
where de [--,-].
(20)
When 7 = 0, 1(a) = > 0 yielding a e [—1,1] and ao = 0. By continuity, the real solution will continue to exist at least for a sufficiently small , while the flow might become unstable for
This solution has been used to test the Mathe-matica code for the general solution (17) - (19). The code was then used to plot other figures of the wavefronts of the ring waves shown below.
1
x
1
1
2
Following Johnson, 1990, we note that the local wave speed in the normal direction to the wave-front, is given by
H
\VH\ Vk2 + k'2'
Therefore, to avoid the appearance of critical layers we require that
° =U (z)cos(d + a),
Vk2 + k'2
where a is the angle between the radial direction and the normal to the wavefront, and cos a = ^k2+k,2. This condition is equivalent to
Fx = — s + U(z)(k cos в — k' sin 9) = 0
for d < z < 1 since F2 = — s = 0. Following [Khusnutdinova and Zhang, 2016a], we can obtain a simple sufficient condition for the absence of critical layers. Indeed,
Fie = —l(z — d)a(k + k'') sin в,
where we assume that 7 > 0, a > 0,d < z < 1 and k + k" > 0 on the selected singular solution (outward propagating wave). Then, Fx has a maximum at = 0, and we require that
Fi < Fi\e=o = —s + ф — d)ak(0) < 0,
implying
7(2 - d)a <
k(0)
which we replace with a stronger condition 7(2 - d)a < 7(1 - d)a < s < "
k(0)
Thus, in order to avoid the appearance of critical layers, in what follows we impose this constraints on 7 for the examples of the shear flow.
In Figure 3 we show the wavefronts of the interfacial ring waves on the linear current
U( ) = ( — d).
The plots on the left are obtained for p1 = 1, p2 = 1.0001, while the plots on the right are for a greater density jump, p1 = 1, p2 = 1.001, with d = 0.7 in all plots. The parameter 7 takes two values 71 = 0 (blue, solid) and 72 = 0.015 (red, dashed). For the
Figure 3. Wavefronts of interfacial ring waves for U(z) = 7(z — d); p1 = 1,p2 = 1.0001 (top) and p2 = 1.001 (bottom); d = 0.7;7 = 0 (blue, solid), 7 = 0.015 (red, dashed).
same strength of the shear flow, the wavefronts appear to be convected by the flow and slightly elongated in the direction of the flow, which is more noticeable in the plots on the top, i.e. for the smaller density jump (and therefore slower interfacial waves; indeed, ^ = 0.00458251 and s2 = 0.0144892).
In Figure 4 we show the wavefronts of the interfacial ring waves on the current
U (z) = 7yfz—~d. Here again the plots on the top are obtained for
(a)
Figure 4. Wavefronts of interfacial ring waves for U(z) = 7\/z — d; pi = 1, p2 = 1.0001 (top) and p2 = 1.001 (bottom); d = 0.7;7 = 0 (blue, solid), 7 = 0.00821584 (red, dashed).
Pi = 1,p2 = 1.0001, while the plots on the bottom are for p1 = 1, p2 = 1.001, with d = 0.7 in all plots. The parameter 7 takes the values 71 =0 (blue, solid) and 72 = 0.00821584 (red, dashed). Thus, the first value is the same as before, while the second value is decreased in order to have the same strength of the current U(z) = 0.0045 on the surface z = 1. The wavefronts appear to be mainly convected in the downstream direction which is again much more pronounced in the plots on the top, i.e. for the smaller density jump.
In Figure 5 we show the wavefronts of the inter-
Figure 5. Wavefronts of interfacial ring waves for U(z) = 7(z - d)2; pi = 1,p2 = 1.0001 (top) and p2 = 1.001 (bottom); d = 0.7;7 = 0 (blue, solid), 7 = 0.05 (red, dashed).
facial ring waves on the current
U(z) =7(z -d)2.
The plots on the top are obtained for pi = 1,p2 = 1.0001, while the plots on the bottom are for p1 = 1,p2 = 1.001, with d = 0.7 in all plots. The parameter 7 takes the values 7i = 0 (blue, solid) and 72 = 0.05 (red, dashed). Thus, the first value is the same as before, while the second value is increased in order to have the same strength of the current U(z) = 0.0045 on the sur-
face = 1. The wavefronts appear to be mainly elongated in the downstream direction which is again much more pronounced in the plots on the top, i.e. for the smaller density jump.
Figure 6. Wavefronts of interfacial ring waves for Pi = 1,p2 = 1.0001, d = 0.7; 7 = 0 (blue, solid), 7 = 0.05 for a = 2 (red, dotted), 7 = 0.015 for a = 1 (red, dashed), and 7 = 0.00821584 for a = 2 (red, dot-dashed). For all three currents U(1) = 0.0045.
In Figure 6 we compare the wavefronts of the interfacial ring waves on the currents
U (z) = 7 (z — d)
where a = -, 1, 2
2' '
for the same set of parameters. Here, 1 = 1, 2 = 1.0001, d = 0.7 and 7 = 0 (blue, solid), 7 = 0.05 for a = 2 (red, dotted), 7 = 0.015 for a = 1 (red, dashed) and 7 = 0.00821584 for a = 1 (red, dot-dashed). All currents have the same strength U( ) = 0.0045 on the surface = 1. It appears that the current with a = 2 convects the ring further downstream than the currents with a = 1 and a = 2, while the current with a = 2 has a stronger effect on the shape of the wavefront elongating it in the direction of the shear flow.
It is also instructive to analyze the 3D vertical structure of the internal wave field by illustrating the dependence of the modal functions in the upper and lower layers on and . The plots in Figure 7 and Figure 8 show the upper- and lower-layer modal functions 01 and 02, respectively, for the linear current U(z) = 7(z — d). The modal functions
Figure 7. Upper-layer modal function 01 of interfacial ring waves for d = 0 (top), d = | (middle) and d = n (bottom). Here, p1 = 1,p2 = 1.0001, d = 0.7; U(z) = 7(z — d) with 7 = 0 (blue, solid) and 7 = 0.015 (red, dashed).
are given by the formulae (9), (10), and they have been normalized to be equal to 1 at d = ^ (orthogonal direction to the current, where the velocity field in the fluid is least affected by the current) and z = d (i.e. on the interface). Here, Pi = 1,P2 = 1.0001, d = 0.7 and 7 = 0 (blue, solid) or 7 = 0.015 (red, dashed). The upper-layer
7 = 0.015, the greatest changes are in the downstream direction, compared to the case when there is no background shear flow, but there is also considerable variation in the upstream direction. The variation in the direction orthogonal to the current is less significant in the upper layer, and negligible in the lower layer. The upper-layer vertical structure in the downstream direction develops a sharp gradient near the interface.
Finally, we note that for the same surface strength U(1) = 0.0045 there exist currents in this family which appear to have wavefronts squeezed in the direction of the shear flow, in contrast to the behaviour illustrated in previous plots. Indeed, this happens, for example, for a = 1/4 and 7 = 0.0060804 (where we continue to choose a in the form of positive and negative powers of 2). The vertical profile of the current, and the corresponding deformation of the wavefront of the ring wave compared to the case when there is no shear flow are shown in Figure 9. This behaviour is similar to the squeezing which was previously reported for internal ring waves propagating in a two-layer fluid with the piecewise-constant shear flow [Khusnutdinova and Zhang, 2016a], where it was linked to the presence of the long-wave instability for plane waves tangent to the ring at 9 = 0 and 9 = ж for a sufficiently strong current [Ovsyannikov, 1979; Ovsyannikov, 1985].
Indeed, the upper-layer current U(z) = 7(z-d)a tends to U(z) =7 as a ^ 0. In this limit, the equation (11) takes the form
т,2 ¡/2 p1 d[—s+7(kcos 9 —к' sin Ö)]2+p2(1 — d) s2
К +гь =-
Figure 8. Lower-layer modal function <fi2 of interfacial ring waves for 9 = 0 (top), 9 = | (middle) and 9 = n (bottom). Here, p1 = 1,p2 = 1.0001, d = 0.7; U(z) = 7(2 — d) with 7 = 0 (blue, solid) and 7 = 0.015 (red, dashed).
modal function $i(z, 9) is shown in Figure 7 for three fixed values of 9: 9 = 0 (downstream direction), 9 = 2 (orthogonal direction to the current) and 9 = k (upstream direction). The lower-layer modal function $2(z, 9) is shown in Figure 8 for the same values of 9. It is evident that the vertical structure strongly depends on the direction. For
( P2 — Pi )d(1 — d)
(21)
and coincides with the equation obtained in [Khusnutdinova and Zhang, 2016a] for the case of a two-layer fluid with the constant current U1 in the upper layer, and another constant current U2 in the lower layer if we let U1 — U2 = 7. Here, s2 is given by (11), and the general solution has the form (15), where
a2+b'2 =cq2a2 — 2aqsa+1, a =
P1
( P2 — P1)(1 — d)
> 0.
The singular solution of (21) has been found in [Khusnutdinova and Zhang, 2016a] under the assumption that the current is sufficiently weak. We
a and b is a parabola), then 1
m =
2\[lks cos 9
(cos2 d + as2 sin2 0), (23)
where d e (-—, — ) .
V 2 2/
Here, k(d) ^ x> as d ^ ±^ leading to the presence of a stationary point on the wavefront at the origin. Finally, if
1 2 2 = 1
a < 7 < Tc = a(1 — as2)
(hyperbolic regime, the locus of parameters a and is a hyperbola), then the singular solution has two branches (corresponding to the right and left parts of the wavefront):
K, (0) = ^
/i \/1—a(1—as2)j2 /n
cos g -x (24)
X a/cos2 d - (a- 1) sin2 d
where e [- arctan
1
'aj2 - 1
arctan
1
'aj2 - 1
Figure 9. Upper-layer current U(z) = 7(z — d)1/4 for d = 0.7 and U(1) = 0.0045 (top), and wave-fronts of interfacial ring waves (bottom) for p\ = 1, p2 = 1.0001, d = 0.7; 7 = 0 (blue, solid) and 7 = 0.0060804 (red, dashed).
It is now instructive to compare the singular solutions (20) and (22) - (24), for U(z) = 7(2 -d) and U(z) = 7, respectively. The solution (20) is real-valued for all values of parameters of the problem, while (22) - (24) is real-valued only for
Y
< =
1
( P2 -P1MP1d + P2(1 -d)\
a(1 - as2)
P1P2
now identify three regimes depending on the strength of the current, and provide detailed analysis up to the onset of the long-wave instability.
If 72 < \ (elliptic regime, the locus of parameters a and is a circle) then the singular solution can be written in the form
m =
/
1 - a(1 - às2)72 cos2 d + (1 - a~f2) sin2 d
1 - a^2 sin20\ a^s
1 - aj2 J 1 - aj2 where de [-—,—}.
cos , (22)
If 72 = i (parabolic regime, the locus of parameters
This condition coincides with the rigid-lid version of the criterion of long-wave instability for co- and counter-propagating plane waves on that shear flow [Ovsyannikov, 1979; Ovsyannikov, 1985] (see also [Bontozoglou, 1991; Boonkasame and Milewski, 2014; Lannes and Ming, 2015] and references therein).
The deformation of the wavefront of a ring wave for the increasing strength of the current is shown in Figure 10. As the strength of the current 7 approaches the critical value 7C, the wavefront appears to be more and more squeezed, until the ring collapses at the critical value. In the elliptic regime a part of the wavefront propagates in the upstream direction, in the parabolic regime the wavefront has a stationary point at the origin, while in the hyperbolic regime all parts of the wavefront propagate in
x
Figure 10. Wavefront of an interfacial ring wave for U(z) = 7, and pi = 1,p2 = 1.0001, d = 0.7; 7 = 0.0045 < = 0.00547723 (red, dashed), 7 =
(red, dotted), <7 = 0.0099 < = 0.00999965 (red, dot-dashed).
the downstream direction, and for each value of 9 in the interval
11 9 e I — arctan —, , arctan —, I
l /Z o ZT ' rz o
a/ <5:72 — 1
a/ <5:72 — 1
the singular solution has two branches.
Since the family of currents approaches the piece-wise-constant current as a ^ 0, we anticipate that the squeezing of the wavefront of the ring wave will persist for the values of a smaller than a certain threshold value. This observation invites studies of the long-wave instability of the plane and ring waves propagating on the background of such currents, but this is beyond the scope of our present paper.
5. Conclusion
In this paper we considered long internal ring waves in a two-layered fluid with a rather general depth-dependent parallel upper-layer current in the rigid-lid approximation. Our main aim was to obtain an analytical solution of the nonlinear firstorder ordinary differential equation responsible for the adjustment of the speed of the long interfacial ring waves in different directions, which was key to the subsequent analysis of the wave field. For a large family of the upper-layer current profiles
described by the function
U(z) = 7(2 - d)a,
where d is the position of the interface, while 7 and a are some positive constant parameters, an explicit analytical solution was obtained in terms of the hypergeometric function. In many natural cases the solution reduces to elementary functions, and we considered examples with a = 2,1, 2 and 1 (in the latter case the corresponding formula in terms of elementary functions is rather long, therefore we did not show it).
The constructed solution has allowed us to illustrate the effects of such shear flows on the wave-fronts and vertical structure of interfacial ring waves with an emphasis on the effects of the density jump and the type and the strength of the current. For the same strength of the shear flow, all flows had greater effect on internal waves for smaller values of the density jump (i.e. slower internal waves). For the currents with a = 1,2, 2 with the same strength on the surface, the current with a = 2 convected the ring further downstream than both other currents, while the current with a = 2 had greater effect on the shape of the wave-front than both other currents, with all currents elongating the ring in the direction of the current. The a = 2 current and other currents with a < 1 could be more representative of the river inflows and exchange flows in straits, while the a = 2 and other currents with a > 1 could be closer to the wind-generated currents. While all currents had the same strength on the surface, the variation of their behaviour in the bulk of the layer had a profoundly different effect on the propagation of the ring waves. The vertical structure of the internal wave field is strongly three-dimensional with the greatest changes due to the current in both layers being in the downstream direction.
We also showed that for the same surface strength, the current with a = | leads to the squeezing of the wavefront of the ring wave in the direction of the current, similarly to the previously reported behaviour for the ring waves on a piecewise-constant current, which we now revisited and described in detail up to the onset of the long-wave instability. We conjecture that this behaviour in the stable regime is generally indicative of the presence of the long-wave instability for sufficiently strong currents in our family with a less than some threshold
value (possibly, a < 1) which we hope to address in our future studies.
Acknowledgments. I would like to use this opportunity to thank Efim Pelinovsky for years of support, advice and friendship.
References
Ablowitz, M. J., H. Segur (1979), On the evolution of packets of water waves, J. Fluid Mech., 92, 691715. Crossref Bontozoglou, V. (1991), Weakly-nonlinear Kelvin-Helmholtz waves between fluids of finite depth, Intl. J. Multiphase Flow, 17, 509-518. Crossref Boonkasame, A., P. A. Milewski (2014), The stability of large-amplitude shallow interfacial non-Boussinesq flows, Stud. Appl. Math., 133, 182-213. Crossref
Ellingsen, S. A. (2014a), Ship waves in the presence of uniform vorticity, J. Fluid Mech., 14-2, R2. Cross-ref
Ellingsen, S. A. (2014b), Initial surface disturbance on a shear current: the Cauchy-Poisson problem with a twist, Phys. Fluids, 26, 082104. Crossref Li, Y., S. A. Ellingsen (2019), A framework
for modelling linear surface waves on shear currents in slowly varying waves, J. Geophys. Res.: Oceans, 124, 2527-2545. Crossref Grimshaw, R. H. J., L. A. Ostrovsky, V.I. Shrira, et al. (1998), Long nonlinear surface and internal gravity waves in a rotating ocean, Surv. Geophys., 19, 289-338. Crossref Grimshaw, R., E. Pelinovsky, T. Talipova (2007), Modelling internal solitary waves in the coastal ocean, Surv. Geophys., 28, 273-298. Crossref Grimshaw, R., E. Pelinovsky, T. Talipova, et al. (2010) Internal solitary waves: propagation, deformation and disintegration, Nonlin. Proc. Geophys., 11, 633649. Crossref Grimshaw, R. (2015), Effect of a background shear current on models for nonlinear long internal waves, Fund. Prikl. Gidrofiz., 8, 20-23. Grimshaw, R. (2019), Initial conditions for the cylindrical Korteweg-de Vries equation, Stud. Appl. Math., 143, 176-191. Crossref Helfrich, K. R., W. K. Melville (2006), Long nonlinear internal waves, Annu. Rev. Fluid Mech., 38, 395-425. Crossref Johnson, R. S. (1980), Water wave and Korteweg - de Vries equations, J. Fluid Mech., 91, 701-719. Crossref
Johnson, R. S. (1990), Ring waves on the surface of shear flows: a linear and nonlinear theory, J. Fluid Mech., 215, 145==160. Crossref
Johnson, R. S. (1997), A Modern Introduction
to the Mathematical Theory of Water Waves, 445 pp. Cambridge University Press, Cambridge. Khusnutdinova, K. R., C. Klein, V. B. Matveev, et al. (2013), On the integrable elliptic cylindrical Kadomtsev-Petviashvili equation, Chaos, 23, 013126. Crossref Khusnutdinova, K. R., X. Zhang (2016a), Long ring waves in a stratified fluid over a shear flow, J. Fluid Mech., 794, 17-44. Crossref Khusnutdinova, K. R., X. Zhang (2016b), Nonlinear ring waves in a two-layer fluid, Physica D, 333, 208221. Crossref Lannes, D., M. Ming (2015), The Kelvin-
Helmholtz instabilities in two-fluids shallow water models, Hamiltonian Partial Differential Equations and Applications, Fields Institute Communications 75, p. 46, Springer, New York. Crossref Lipovskii, V. D. (1985), On the nonlinear internal wave theory in fluid of finite depth, Izv. Akad. Nauk SSSR, Ser. Fiz, 21, 864-871. McMilan, J. M., B. R. Sutherland (2010), The life-cycle of axisymmetric internal solitary waves, Nonlin. Proc. Geophys., 17, 443-453. Crossref Miles, J. W. (1978), An axisymmetric Boussinesq
wave, J. Fluid Mech., 84, 181-191. Crossref Nash, J. D., J. N. Moum (2005), River plums as a source of large amplitude internal waves in the coastal ocean, Nature, 437, 400-403. Crossref Ovsyannikov, L. V. (1979), Two-layer'shallow water' model, J. Appl. Math. Tech. Phys., 20, 127-135. Crossref
Ovsyannikov, L. V., et al. (1985), Nonlin-
ear Problems in the Theory of Surface and Internal Waves, Nauka, Novosibirsk. Ramirez, C., D. Renouard, Yu. A. Stepanyants (2002), Propagation of cylindrical waves in a
rotating fluid, Fluid Dynam. Res., 30, 169-196. Crossref
Smeltzer, B. K., E. Esoy, S. A. Ellingsen (2019), Observation of surface wave patterns modified by subsurface shear currents, J. Fluid Mech., 873, 508-530. Crossref
Vlasenko, V., J. C. Sanchez Garrido, N. Staschuk, et al. (2009), Three-dimensional evolution of large-amplitude internal waves in the Strait of Gibraltar, J. Phys. Oceanogr., 39, 2230-2246. Crossref Vlasenko, V., N. Staschuk, et al. (2013), Generation of baroclinic tides over an isolated underwater bank, J. Geophys. Res., 118, 4395-4408. Crossref Weidman, P. D., R. Zakhem, et al. (1988), Cylindrical solitary waves, J. Fluid Mech., 191, 557-573. Crossref
Corresponding author:
K. R. Khusnutdinova, Department of Mathematical Sciences, Loughborough University, Loughborough, Great Britain ([email protected])
