Journal topic
Nonlin. Processes Geophys., 26, 283–290, 2019
https://doi.org/10.5194/npg-26-283-2019
Nonlin. Processes Geophys., 26, 283–290, 2019
https://doi.org/10.5194/npg-26-283-2019

Research article 21 Aug 2019

Research article | 21 Aug 2019

# Explosive instability due to flow over a rippled bottom

Explosive instability due to flow over a rippled bottom
Anirban Guha1,2,* and Raunak Raj1,* Anirban Guha and Raunak Raj
• 1Environmental and Geophysical Fluids Group, Department of Mechanical Engineering, Indian Institute of Technology, Kanpur, U.P. 208016, India
• 2Institute of Coastal Research, Helmholtz-Zentrum Geesthacht, Geesthacht 21502, Germany
• *These authors contributed equally to this work.

Correspondence: Anirban Guha (anirbanguha.ubc@gmail.com)

Abstract

In this paper, we study Bragg resonance, i.e., the triad interaction between surface and/or interfacial waves with a bottom ripple, in the presence of background velocity. We show that when one of the constituent waves of the triad has negative energy, the amplitudes of all the waves grow exponentially. This is very different from classic Bragg resonance in which one wave decays to cause the growth of the other. The instabilities we observe are “explosive” and are different from normal mode shear instabilities since our velocity profiles are linearly stable. Our work may explain the existence of large-amplitude internal waves over periodic bottom ripples in the presence of tidal flow observed in oceans and estuaries.

1 Introduction

The energy exchange between two counter-propagating surface gravity waves mediated by a bottom ripple, otherwise known as the “Bragg resonance”, is a widely known phenomenon in oceanography and coastal engineering (Davies1982; Mei1985; Kirby1986). Bragg resonance strongly affects the wave spectrum in continental shelves and coastal regions (Ball1964), modifies the shore-parallel sandbars, and protects the shoreline from wave attacks . Actually, Bragg resonance is a special kind of resonant triad in which one of the constituent waves is the bottom ripple, which acts as a stationary wave . Usually, in a classical resonant triad (hereafter, “resonant” is suppressed for brevity), one wave gives energy to the other two waves via nonlinear interactions such that the individual wave amplitudes remain bounded at all times.

The present work aims at studying explosive Bragg resonance which occurs as a consequence of the presence of a velocity field. We find that contrary to typical Bragg resonances, in explosive Bragg resonance the amplitude of the waves is not limited by energy considerations. This may explain the presence of large-amplitude internal waves over periodic bottom ripples such as those in the Rotterdam waterway . Such a large-amplitude response was primarily seen during the strong flood tides and the presence of a velocity field was deemed essential for the same. Further, this mechanism of wave generation may also explain the presence of high-amplitude internal waves in continental shelves whose amplitudes are uncorrelated with tidal forcing. We also note here that in the presence of a vertically nonuniform current, a work similar to ours but for capillary–gravity waves was performed by . In the present paper, one of our arguments will be that shear is not necessary for such explosive growth and that even a uniform mean flow (flowing over a stationary bottom boundary) is sufficient.

In order to explain the explosive Bragg resonance, we first consider in Sect. 2 a single-layered flow. Here, the system consists of a surface wave propagating over a rippled bottom topography and the fluid flowing with a mean velocity profile as shown in Fig. 1a. For explosive instability, we expect this wave to grow along with its explosive counterpart while simultaneously conserving the energy of the system. After highlighting the importance of mean flow in explosive Bragg resonance and detailing various aspects of negative energy waves, we briefly examine a more realistic two-layered flow scenario (Fig. 1b) in Sect. 3. Here we explain why internal (interfacial) waves would be more susceptible to such resonances. Continuously stratified flows that can support internal gravity waves have not been studied. We expect explosive Bragg resonance to be a general phenomenon that can therefore be realized in a variety of systems.

2 Theory

A wave of wavenumber k can interact with an undulated bottom of wavenumber kb to transfer some of its energy to the wavenumbers k±kb. However, this energy transfer is maximized when the following resonance condition is satisfied:

$\begin{array}{}\text{(1a)}& & {k}_{\mathrm{i}}±{k}_{\mathrm{r}}±{k}_{\mathrm{b}}=\mathrm{0},\text{(1b)}& & {\mathit{\omega }}_{\mathrm{i}}±{\mathit{\omega }}_{r}=\mathrm{0},\end{array}$

where ω denotes frequency and the subscripts i and r, respectively, denote incident and resonant waves. For a single-layered flow (i.e., with no density variation) in the absence of any mean flow, only one such resonant set can exist at the first order of nonlinearity (Davies1982). In this case, an incident surface gravity wave of wavenumber ki resonantly interacts with the bottom of twice the wavenumber (i.e., kb=2ki) to generate exactly one surface gravity wave, which has a wavenumber kr=ki and travels in the direction opposite to the incident wave. This is the classic Bragg resonance condition for surface waves. The corresponding dispersion relation is shown in Fig. 2a. Each wave has been vectorially represented as (k,ω) (marked by arrows) in the dispersion diagram, and the corresponding coordinates (or vector tips) are marked by discs “”. In this vectorial representation, the “Bragg triad” forms a vector triangle (in fact, any resonant triad at the first order of nonlinearity forms a vector triangle in the dispersion diagram). We observe that the dispersion relation for a surface gravity wave is symmetrical about the k axis; i.e., there is no difference between the positively and negatively traveling waves apart from the direction of propagation.

Figure 1Schematic diagram showing waves over a rippled bottom. The profile of the bottom is ηb. (a) One-layered flow of depth H with density ρ is assumed. The free surface profile is η. The velocity profile (chosen to be linear) must create a Doppler shift between the free surface and the bottom. (b) Two-layered flow of depth $H\equiv {h}_{\mathrm{u}}+{h}_{\mathrm{l}}$ with density ρu and ρl is assumed. The free surface (interface) profile is ηu (ηl). The velocity profile is piecewise linear.

In the presence of a velocity field, however, the symmetry between the rightward and leftward traveling waves is destroyed. We find that this not only leads to a modification in the resonance conditions, but also leads to the formation of new resonant triads . We assume a single-layered flow with a “Couette”-type mean velocity profile, i.e., the mean velocity increases linearly from $\stackrel{\mathrm{‾}}{u}=\mathrm{0}$ at the bottom ($z=-H$) to $\stackrel{\mathrm{‾}}{u}=U$ at the surface (z=0); see Fig. 1a. The mean shear is denoted by $\mathrm{\Omega }\equiv \mathrm{d}\stackrel{\mathrm{‾}}{u}/\mathrm{d}z=U/H$ and is a constant in this situation. We emphasize here that, although the actual velocity profile is of some relevance, what matters most is the Doppler shift between the bottom ripples and the surface. Therefore, even a uniform current (flowing over a stationary bottom boundary) would have worked just fine. Mathematically, the dispersion relation in this case is given by

$\begin{array}{}\text{(2)}& {\mathit{\omega }}_{\mathrm{in}}^{\mathrm{2}}+\mathrm{\Omega }\mathrm{tanh}\left(kH\right){\mathit{\omega }}_{\mathrm{in}}-gk\mathrm{tanh}\left(kH\right)=\mathrm{0}.\end{array}$

Here ${\mathit{\omega }}_{\mathrm{in}}\equiv \mathit{\omega }-Uk$ is defined as the intrinsic frequency of the wave, which basically means the frequency obtained after subtracting the Doppler shift produced by the mean flow. For this paper, without loss of generality, we have considered a positive k and have allowed the frequency ω to be either positive or negative. As a consequence, ω>0 means a positively traveling wave, while ω<0 means a negatively traveling wave in the stationary reference frame. Equation (2), a quadratic equation in ωin, reveals that the product of the two roots is negative. Hence, the two branches of the dispersion curves, shown in Fig. 2b, have intrinsic frequencies of opposite signs. Note that in Fig. 2b, we have plotted ω vs. k (and not ωin vs. k) in the nondimensional form and have taken the surface velocity U>0 without loss of generality. The positive intrinsic frequency branch has been labeled 𝒮𝒢+ and the negative one 𝒮𝒢. Below we perform a detailed analysis of the dispersion curve.

## 2.1 An analysis of triads using dispersion curves

From Fig. 2b, we see that while ω for the 𝒮𝒢+ branch increases monotonically with k, the same is not true for the 𝒮𝒢 branch. For this branch, ω initially decreases with k, attains a minimum, and then starts to increase. Thus the sign of ω becomes opposite to ωin; in Fig. 2b this happens after kH≈3 for ${U}^{*}\equiv U/\sqrt{gH}=\mathrm{0.7}$. The non-monotonic behavior of ω with k allows a given wave vector (k,ω) to form multiple resonant triad sets. This is starkly different from the classic Bragg resonance, where, as mentioned earlier, only one triad condition is possible for a single-layered fluid (Fig. 2a).

The minimum frequency of the 𝒮𝒢 branch1 is labeled ωmin. The wavenumber at which the 𝒮𝒢 branch crosses the k axis (ω=0 line) is labeled kz. For every point on the dispersion curve with frequency $\mathit{\omega }<|{\mathit{\omega }}_{\mathrm{min}}|$ (shown in the shaded region in Fig. 2b), there exist three resonant triads for any given incident wavenumber. For an example, we choose a frequency ${\mathit{\omega }}_{\mathrm{0}}<|{\mathit{\omega }}_{\mathrm{min}}|$ and plot $\mathit{\omega }=±{\mathit{\omega }}_{\mathrm{0}}$ on the dispersion curve. There will be four intersections with the dispersion curve, all shown using bullets () in Fig. 2b. An interesting observation here is that three of these above-mentioned points lie on the same branch of the dispersion curve 𝒮𝒢. Any two of the four intersection points form a resonant triad with an appropriately chosen bottom ripple. If the points lie on the opposite side of the k axis, the bottom's wavenumber for resonance would be the sum of two wavenumbers involved, else it would be the difference. For a given $|{\mathit{\omega }}_{\mathrm{0}}|$ less than $|{\mathit{\omega }}_{\mathrm{min}}|$, a total of six ( 4C2) Bragg triads can be obtained. We again recall here that in the absence of mean velocity, only one triad would have been possible, as shown in Fig. 1a.

If we choose ${\mathit{\omega }}_{\mathrm{0}}>-{\mathit{\omega }}_{\mathrm{min}}$, there would be only one such triad between the 𝒮𝒢+ branch and the 𝒮𝒢 branch. In that case, since both the incident and resonant frequencies will be positive, the wavenumber of the bottom ripple would be the difference of the two wavenumbers. However, for the case when ω0<ωmin, a stably propagating wave cannot exist and there will not be any chance of resonance.

Figure 2(a) A Bragg triad in the absence of mean velocity. Wavenumber k has been nondimensionalized by 1∕H (H flow depth), and frequency ω by $\sqrt{g/H}$, i.e., ${\mathit{\omega }}^{*}\equiv \mathit{\omega }/\sqrt{g/H}$. (b) Bragg triads in the presence of “Couette-type” mean velocity with ${U}^{*}\equiv U/\left(\sqrt{gH}\right)=\mathrm{0.7}$. In the shaded region, three Bragg triads can exist for a given k. For labels, see the text.

## 2.2 Negative energy waves

As shown in , the energy per unit area E of a wave (hereafter, simply referred to as “energy”) with frequency ω, amplitude a, and satisfying the dispersion relation $\mathfrak{D}\left(\mathit{\omega },k\right)=\mathrm{0}$, is given by

$\begin{array}{}\text{(3)}& E=\frac{\mathrm{1}}{\mathrm{4}}\mathit{\omega }\frac{\partial \mathfrak{D}}{\partial \mathit{\omega }}|a{|}^{\mathrm{2}}.\end{array}$

The energy of a wave in the words of is “the work done during an idealized process in which the wave is driven up by an external force applied on a surface in the fluid.” Therefore, a negative energy wave is such a wave whose introduction into the system lowers the energy of the system. For a given system, the dispersion relation $\mathfrak{D}\left(\mathit{\omega },k\right)=\mathrm{0}$ can be written in multiple forms; for every form, the factor $\partial \mathfrak{D}/\partial \mathit{\omega }$ would be different. However, only when the dispersion relation for the purpose of calculating energy is written in the form described in do we obtain the energy of the wave. In this approach an interface $z=\mathit{\eta }\left(x,t\right)$ is chosen, and the pressure p1 (p2) just above (below) this interface is written as

$\begin{array}{}\text{(4a)}& & {p}_{\mathrm{1}}\left(x,t\right)={\mathfrak{D}}_{\mathrm{1}}\left(\mathit{\omega },k\right)\mathit{\eta }\left(x,t\right),\text{(4b)}& & {p}_{\mathrm{2}}\left(x,t\right)={\mathfrak{D}}_{\mathrm{2}}\left(\mathit{\omega },k\right)\mathit{\eta }\left(x,t\right).\end{array}$

The dispersion relation is given by

$\begin{array}{}\text{(5)}& \mathfrak{D}\left(\mathit{\omega },k\right)={\mathfrak{D}}_{\mathrm{1}}\left(\mathit{\omega },k\right)-{\mathfrak{D}}_{\mathrm{2}}\left(\mathit{\omega },k\right).\end{array}$

The dispersion relation of an intermediate depth surface gravity wave in the presence of a constant mean shear current is found to be

$\begin{array}{}\text{(6)}& \mathfrak{D}\left(\mathit{\omega },k\right)=\mathit{\rho }\left[\frac{\left(\mathit{\omega }-Uk{\right)}^{\mathrm{2}}}{k\mathrm{tanh}\left(kH\right)}+\frac{\left(\mathit{\omega }-Uk\right)\mathrm{\Omega }}{k}-g\right]=\mathrm{0},\end{array}$

where ρ is the density of the fluid. Hence,

$\begin{array}{}\text{(7)}& \frac{\partial \mathfrak{D}}{\partial \mathit{\omega }}=\mathit{\rho }\frac{{\mathit{\omega }}_{\mathrm{in}}^{\mathrm{2}}+gk\mathrm{tanh}\left(kH\right)}{{\mathit{\omega }}_{\mathrm{in}}k\mathrm{tanh}\left(kH\right)}.\end{array}$

We observe that in this case, the sign of $\partial \mathfrak{D}/\partial \mathit{\omega }$ is dependent only on the sign of ωin, i.e., the intrinsic frequency of the wave. Therefore, for the 𝒮𝒢+ branch, $\partial \mathfrak{D}/\partial \mathit{\omega }$ is always positive, whereas for the 𝒮𝒢 branch, it is always negative. Further, for the 𝒮𝒢+ branch, the frequency ω is always positive. Therefore, according to Eq. (3), the energy of the 𝒮𝒢+ branch is always positive (E>0). For the 𝒮𝒢 branch, ω<0 when k<kz, which implies positive energy for k<kz. However, for k>kz, we have ω>0, but $\partial \mathfrak{D}/\partial \mathit{\omega }$ still remains negative, which implies a negative energy wave (E<0). Therefore, using the negative energy approach of Cairns, we expect that a wave on the branch 𝒮𝒢, for which k>kz, will form an explosive Bragg triad with a suitable positive energy wave.

## 2.3 An explosive Bragg resonance pair

To derive the amplitude evolution equations for Bragg resonance, we represent the constituent surface waves as (k1,ω1) and (k2,ω2). The waves are expressed in the form ${a}_{\mathrm{1}}\left(t\right)\mathrm{exp}\left[\mathrm{i}\left({k}_{\mathrm{1}}x-{\mathit{\omega }}_{\mathrm{1}}t\right)\right]+\mathrm{c}.\mathrm{c}$ and ${a}_{\mathrm{2}}\left(t\right)\mathrm{exp}\left[\mathrm{i}\left({k}_{\mathrm{2}}x-{\mathit{\omega }}_{\mathrm{2}}t\right)\right]+\mathrm{c}.\mathrm{c}$, where c.c denotes the complex conjugate. With a stationary bottom ripple with a wavenumber kb, they satisfy the resonance condition

$\begin{array}{}\text{(8a)}& & {k}_{\mathrm{1}}+{k}_{\mathrm{2}}={k}_{\mathrm{b}},\text{(8b)}& & {\mathit{\omega }}_{\mathrm{1}}+{\mathit{\omega }}_{\mathrm{2}}=\mathrm{0}.\end{array}$

The amplitude evolution equations for Bragg resonance for such a case are given by :

$\begin{array}{}\text{(9a)}& & \frac{\mathrm{d}{a}_{\mathrm{1}}}{\mathrm{d}t}={\mathit{\beta }}_{\mathrm{1}}{a}_{\mathrm{b}}{\stackrel{\mathrm{‾}}{a}}_{\mathrm{2}},\text{(9b)}& & \frac{\mathrm{d}{a}_{\mathrm{2}}}{\mathrm{d}t}={\mathit{\beta }}_{\mathrm{2}}{a}_{\mathrm{b}}{\stackrel{\mathrm{‾}}{a}}_{\mathrm{1}},\end{array}$

where a1 and a2 are the complex amplitudes of the waves involved and ab is the complex amplitude of the bottom ripple with overbars denoting the complex conjugates. The coefficients β1 and β2 in Eqs. (9a) and (9b) can be derived following Fredholm's alternative, the procedure of which has been elaborated in . In this particular case of uniform shear, the coefficients are as follows:

$\begin{array}{}\text{(10a)}& & {\mathit{\beta }}_{\mathrm{1}}=\mathrm{i}\frac{{k}_{\mathrm{1}}{\left({\mathit{\omega }}_{\mathrm{1}}-U{k}_{\mathrm{1}}\right)}^{\mathrm{2}}\left({\mathit{\omega }}_{\mathrm{2}}-U{k}_{\mathrm{2}}\right)}{\mathrm{cosh}{k}_{\mathrm{1}}H\mathrm{sinh}{k}_{\mathrm{2}}H\left[g{k}_{\mathrm{1}}\mathrm{tanh}{k}_{\mathrm{1}}H+{\left({\mathit{\omega }}_{\mathrm{1}}-U{k}_{\mathrm{1}}\right)}^{\mathrm{2}}\right]},\text{(10b)}& & {\mathit{\beta }}_{\mathrm{2}}=\mathrm{i}\frac{{k}_{\mathrm{2}}{\left({\mathit{\omega }}_{\mathrm{2}}-U{k}_{\mathrm{2}}\right)}^{\mathrm{2}}\left({\mathit{\omega }}_{\mathrm{1}}-U{k}_{\mathrm{1}}\right)}{\mathrm{cosh}{k}_{\mathrm{2}}H\mathrm{sinh}{k}_{\mathrm{1}}H\left[g{k}_{\mathrm{2}}\mathrm{tanh}{k}_{\mathrm{2}}H+{\left({\mathit{\omega }}_{\mathrm{2}}-U{k}_{\mathrm{2}}\right)}^{\mathrm{2}}\right]}.\end{array}$

The coefficients β1 and β2 are purely imaginary. Furthermore, it can be easily seen that the signs of β1 and β2, respectively, depend on the signs of ω1Uk1 and ω2Uk2 only, i.e., only on the intrinsic frequencies of the respective waves. Thus the sign of the product β1β2 depends only on the sign of the product of the intrinsic frequencies of the waves forming the resonant pair. As noted before, the intrinsic frequency of the 𝒮𝒢+ branch is positive for all k, while that of the 𝒮𝒢 branch is negative for all k.

We can also express Eqs. (9a) and (9b) as

$\begin{array}{}\text{(11a)}& & \frac{{d}^{\mathrm{2}}{a}_{\mathrm{1}}}{\mathrm{d}{t}^{\mathrm{2}}}=-{\mathit{\beta }}_{\mathrm{1}}{\mathit{\beta }}_{\mathrm{2}}|{a}_{\mathrm{b}}{|}^{\mathrm{2}}{a}_{\mathrm{1}},\text{(11b)}& & \frac{{d}^{\mathrm{2}}{a}_{\mathrm{2}}}{\mathrm{d}{t}^{\mathrm{2}}}=-{\mathit{\beta }}_{\mathrm{2}}{\mathit{\beta }}_{\mathrm{1}}|{a}_{\mathrm{b}}{|}^{\mathrm{2}}{a}_{\mathrm{2}}.\end{array}$

Hence, for explosive growth,

$\begin{array}{}\text{(12)}& -{\mathit{\beta }}_{\mathrm{1}}{\mathit{\beta }}_{\mathrm{2}}>\mathrm{0}⇒\left({\mathit{\omega }}_{\mathrm{1}}-U{k}_{\mathrm{1}}\right)\left({\mathit{\omega }}_{\mathrm{2}}-U{k}_{\mathrm{2}}\right)>\mathrm{0}.\end{array}$

This basically means that the intrinsic frequencies have to be of the same sign, or in other words, the two points must lie on the same branch of the dispersion curve in Fig. 2b. Furthermore, Eqs. (8a) and (8b) assume that the actual frequencies of the resonant pair must be of opposite signs. Therefore we deduce that for Bragg resonance to give rise to explosive growth, if the actual frequencies of the waves are of opposite signs, then the intrinsic frequencies must have the same sign. This is precisely the condition satisfied by all the Bragg triads formed by the waves with $\mathrm{0}<\mathit{\omega }<|{\mathit{\omega }}_{\mathrm{min}}|$ on the 𝒮𝒢 branch, with the waves on the same branch with ${\mathit{\omega }}_{\mathrm{min}}<\mathit{\omega }<\mathrm{0}$.

## 2.4 Explosive Bragg resonance from the negative energy perspective

To understand how the waves in Sect. 2.3 constitute a positive energy–negative energy pair, we write the coefficients β1 and β2 in terms of $\partial \mathfrak{D}/\partial \mathit{\omega }$. Using Eq. (7) along with Eqs. (10a) and (10b), we obtain

$\begin{array}{}\text{(13a)}& & {\mathit{\beta }}_{\mathrm{1}}=i\mathit{\lambda }{\left(\frac{\partial \mathfrak{D}}{\partial \mathit{\omega }}\right)}_{{\mathit{\omega }}_{\mathrm{1}},{k}_{\mathrm{1}}}^{-\mathrm{1}},\text{(13b)}& & {\mathit{\beta }}_{\mathrm{2}}=i\mathit{\lambda }{\left(\frac{\partial \mathfrak{D}}{\partial \mathit{\omega }}\right)}_{{\mathit{\omega }}_{\mathrm{2}},{k}_{\mathrm{2}}}^{-\mathrm{1}},\end{array}$

where $\mathit{\lambda }=\left({\mathit{\omega }}_{\mathrm{1}}-U{k}_{\mathrm{1}}\right)\left({\mathit{\omega }}_{\mathrm{2}}-U{k}_{\mathrm{2}}\right)\left[\mathrm{sinh}\left({k}_{\mathrm{1}}H\right)\mathrm{sinh}{k}_{\mathrm{2}}H{\right]}^{-\mathrm{1}}$. Therefore Eqs. (9a) and (9b) can be written as

$\begin{array}{}\text{(14a)}& & {\left(\frac{\partial \mathfrak{D}}{\partial \mathit{\omega }}\right)}_{{\mathit{\omega }}_{\mathrm{1}},{k}_{\mathrm{1}}}\frac{\mathrm{d}{a}_{\mathrm{1}}}{\mathrm{d}t}=i\mathit{\lambda }{a}_{\mathrm{b}}{\overline{a}}_{\mathrm{2}},\text{(14b)}& & {\left(\frac{\partial \mathfrak{D}}{\partial \mathit{\omega }}\right)}_{{\mathit{\omega }}_{\mathrm{2}},{k}_{\mathrm{2}}}\frac{\mathrm{d}{a}_{\mathrm{2}}}{\mathrm{d}t}=i\mathit{\lambda }{a}_{\mathrm{b}}{\overline{a}}_{\mathrm{1}}.\end{array}$

Although we have written the expressions in this form for a special case of a single-layered flow, such a form is very general. In fact, for the case of triad interactions involving three waves, a similar set of equations was obtained by . Using Eqs. (13a) and (13b) along with the condition for explosive instability, i.e., Eq. (12), we get

$\begin{array}{}\text{(15)}& {\left(\frac{\partial \mathfrak{D}}{\partial \mathit{\omega }}\right)}_{{\mathit{\omega }}_{\mathrm{1}},{k}_{\mathrm{1}}}{\left(\frac{\partial \mathfrak{D}}{\partial \mathit{\omega }}\right)}_{{\mathit{\omega }}_{\mathrm{2}},{k}_{\mathrm{2}}}>\mathrm{0}.\end{array}$

This basically means that energy coefficients $\partial \mathfrak{D}/\partial \mathit{\omega }$ of the two waves are of the same sign. However, we know from Eqs. (8a) and (8b) that the frequencies are of opposite signs. Hence, the energy of the two waves, given by Eq. (3), must be of opposite signs for an explosive instability to occur.

## 2.5 Another explosive Bragg resonance pair

Here we consider the case where the two frequencies are of the same sign. The resonance condition in this case is given by

$\begin{array}{}\text{(16a)}& & {k}_{\mathrm{2}}-{k}_{\mathrm{1}}={k}_{\mathrm{b}},\text{(16b)}& & {\mathit{\omega }}_{\mathrm{2}}-{\mathit{\omega }}_{\mathrm{1}}=\mathrm{0}.\end{array}$

The amplitude evolution equation is found to be

$\begin{array}{}\text{(17a)}& & \frac{\mathrm{d}{a}_{\mathrm{1}}}{\mathrm{d}t}={\mathit{\beta }}_{\mathrm{1}}{\stackrel{\mathrm{‾}}{a}}_{\mathrm{b}}{a}_{\mathrm{2}},\text{(17b)}& & \frac{\mathrm{d}{a}_{\mathrm{2}}}{\mathrm{d}t}={\mathit{\beta }}_{\mathrm{2}}{a}_{\mathrm{b}}{a}_{\mathrm{1}},\end{array}$

where β1 and β2 remain the same as that in Eqs. (11a) and (11b). The above equations can be written as

$\begin{array}{}\text{(18a)}& & \frac{{d}^{\mathrm{2}}{a}_{\mathrm{1}}}{\mathrm{d}{t}^{\mathrm{2}}}={\mathit{\beta }}_{\mathrm{1}}{\mathit{\beta }}_{\mathrm{2}}|{a}_{\mathrm{b}}{|}^{\mathrm{2}}{a}_{\mathrm{1}},\text{(18b)}& & \frac{{d}^{\mathrm{2}}{a}_{\mathrm{2}}}{\mathrm{d}{t}^{\mathrm{2}}}={\mathit{\beta }}_{\mathrm{2}}{\mathit{\beta }}_{\mathrm{1}}|{a}_{\mathrm{b}}{|}^{\mathrm{2}}{a}_{\mathrm{2}},\end{array}$

which implies that the condition for explosive growth is

$\begin{array}{}\text{(19)}& {\mathit{\beta }}_{\mathrm{1}}{\mathit{\beta }}_{\mathrm{2}}>\mathrm{0}⇒\left({\mathit{\omega }}_{\mathrm{1}}-U{k}_{\mathrm{1}}\right)\left({\mathit{\omega }}_{\mathrm{2}}-U{k}_{\mathrm{2}}\right)<\mathrm{0}.\end{array}$

Hence, if the actual frequencies of the waves are of the same sign, the intrinsic frequencies must be of opposite signs for explosive instability to occur. As we have mentioned earlier, having opposite signs of intrinsic frequency in this case means that the waves must be on different branches, i.e., one on 𝒮𝒢+ and the other on 𝒮𝒢, but have the same sign of actual frequency. As an example, this will correspond to the Bragg triad formed by the two black dots in the upper-half plane (i.e., ω>0) in Fig. 2b. Using Eqs. (13a), (13b) and (19) along with Eqs. (16a) and (16b), we will indeed find that the energy of the two waves must be of opposite signs for the occurrence of explosive instability.

Figure 3Growth of amplitudes of a wave on the 𝒮𝒢 branch which resonates a wave on the same branch. The incident wave wavenumber is k1H=0.06, while the resonant wave has the wavenumber krH=3.13. Other relevant parameters are kbH=3.19, ${\mathit{\omega }}_{\mathrm{1}}/\sqrt{g/H}=-\mathrm{0.0213}$, ${\mathit{\omega }}_{\mathrm{2}}/\sqrt{g/H}=\mathrm{0.0213}$, $U/\sqrt{gH}=\mathrm{0.70}$, N=2048, M=2, ${a}_{\mathrm{1}s}/H=\mathrm{2.5}×{\mathrm{10}}^{-\mathrm{5}}$, ${a}_{\mathrm{b}}/H=\mathrm{0.025}$, and ${T}_{\mathrm{1}}/\mathrm{\Delta }T=\mathrm{1024}$. Here, “N” denotes the number of Fourier modes used in numerical simulation, whereas “M” denotes the order of nonlinearity taken into account. (See for details.)

## 2.6 Generalized analysis

We have already shown the existence of Bragg resonant triads that give rise to explosive growth for a single-layered flow with shear present. We have shown for this case that for explosive instability, the dispersion curve must cross the ω=0 axis for explosive instability (assuming the reference frame attached to the bottom). Furthermore, if the intrinsic frequency of a wave is ωin, then for it to change sign, a velocity U in the opposite direction must be present, implying that for some value of k, $|Uk|>|{\mathit{\omega }}_{\mathrm{in}}|$. Given the fact that the frequency of a gravity wave varies as ${\mathit{\omega }}_{\mathrm{in}}\sim {k}^{\mathrm{1}/\mathrm{2}}$ and that for a vorticity wave, ωink0, for some value of k, $|Uk|$ will exceed $|{\mathit{\omega }}_{\mathrm{in}}|$. Therefore, mathematically, an explosive triad will always exist for any value of U. However, as can be seen from Eq. (10a), the coefficients β1 and β2 are proportional to (coshk1Hsinhk2H)−1 and (sinhk1Hcoshk2H)−1, respectively; hence, for higher values of k, they rapidly tend towards 0. Thus, for very low velocities, explosive triad conditions may be satisfied for a very large value of k, but the growth rate tends to zero in such cases. Physically, higher values of kH refer to the deep water case where waves are far away from the bottom and hence do not feel the bottom's effect, even though the resonance condition is satisfied, technically. Further, in a general setting where different interfaces are moving at different velocities, it may not be possible to form a polynomial equation in terms of intrinsic frequency. In any case, the energy of waves can always be found out and, hence, the necessary condition for the existence of an explosive Bragg pair is that a branch of the dispersion curve with negative energy must exist.

Finally, we point out that for the case of explosive triads between three waves (discussed in ), the existence of negative energy did not necessarily mean explosive interaction. The explosive instability occurs when out of the three waves involved, the wave of greatest frequency has energy of opposite sign from the other. However, in this case of Bragg resonance, only two waves are involved, with the third wave being the bottom ripple with zero frequency and the other two waves having the same magnitude of frequency. Therefore, if the bottom is at rest, Bragg resonance involving waves with opposite energy will necessarily mean an explosive growth.

Figure 4(a) Dispersion curves for a uniform current (flowing over a stationary bottom boundary). ${U}_{\mathrm{u}}^{*}={U}_{\mathrm{l}}^{*}={U}_{\mathrm{b}}^{*}=\mathrm{0.1}$, ${h}_{\mathrm{u}}/{h}_{\mathrm{l}}=\mathrm{1}/\mathrm{3}$, and ${\mathit{\rho }}_{\mathrm{u}}/{\mathit{\rho }}_{\mathrm{l}}=\mathrm{0.95}$. (b) Dispersion curves for shear in the bottom layer and single density fluid. ${U}_{\mathrm{u}}^{*}={U}_{\mathrm{l}}^{*}=\mathrm{0.1}$, ${U}_{\mathrm{b}}^{*}=\mathrm{0}$, ${h}_{\mathrm{u}}/{h}_{\mathrm{l}}=\mathrm{1}/\mathrm{3}$, and ρu=ρl. Dashed curves have negative energy. U has been non-dimensionalized with $\sqrt{gH}$, where $H={h}_{\mathrm{u}}+{h}_{\mathrm{l}}$.

We have also numerically simulated a case of explosive instability using a higher-order spectral (HOS) code , which was extended by to incorporate shear into it. Both the incident wave and the resonant waves are on 𝒮𝒢. The incident wave with a wavenumber k1H=0.06 interacted with the bottom ripples of wavenumber kbH=3.19 to resonate a wave with wavenumber k2H=3.13. Other relevant parameters are mentioned in the caption of the figure. It can be seen that within 20 time periods, the amplitude of the resonant wave has become over 10 times that of the initial amplitude of the incident wave.

3 Explosive resonance in a two-layered flow

In this section, we consider a two-layered flow with a density ρu and mean shear Ωu in the upper layer and a density ρl and mean shear Ωl below it. By a two-layered flow we do not necessarily mean that there has to be a density difference between the two layers. We simply mean that there has to be either a shear jump or a density jump (or both) between the two layers, which may lead to a perturbation vorticity generation at the interface. We have already discussed the theory of explosive Bragg resonance; this section mainly shows that in a two-layered setting, achieving the explosive Bragg resonance needs significantly less velocity. In a two-layered density stratification without any velocity field, there exist four different modes of wave propagation – two surface modes and two interfacial modes – propagating symmetrically in both directions. Whereas, for a one-layered setting, the energy transfer was limited to the surface only, in this case, waves on the pycnocline may also participate in the energy exchange. The intrinsic frequency of a gravity wave on a pycnocline is significantly less than that of the surface because the density contrast at the air–water interface is far more than the density contrast at the pycnocline. Therefore, the condition for explosive instability, i.e., $|Uk|>|{\mathit{\omega }}_{\mathrm{in}}|$, may be satisfied at significantly lower values of U and k for the interfacial mode. As can be seen from Fig. 4a, a uniform current with a small Froude number of 0.1 may also lead to the negative energy branch (ℐ𝒢) and, consequently, any Bragg triad involving a negative energy branch will be prone to explosive growth. One such Bragg triad has been shown using discs “” in the figure.

Figure 5(a) Incident interfacial mode wave packet with k1H=2.1 at t=0. (b) The incident mode grows while simultaneously resonating with a surface mode with k2H=0.1. Other relevant parameters are ${k}_{\mathrm{b}}H=\mathrm{2.0},{U}_{\mathrm{u}}^{*}={U}_{\mathrm{l}}^{*}=\mathrm{0.1565},{U}_{\mathrm{b}}^{*}=\mathrm{0},{h}_{\mathrm{u}}/{h}_{\mathrm{l}}=\mathrm{1}/\mathrm{3},{\mathit{\rho }}_{\mathrm{u}}/{\mathit{\rho }}_{\mathrm{l}}=\mathrm{0.95}$, ${\mathit{\omega }}_{\mathrm{1}}^{*}={\mathit{\omega }}_{\mathrm{2}}^{*}=\mathrm{0.1095}$, N=1024, M=2, ai=0.01hu, ab=0.1H, and ${T}_{\mathrm{1}}/\mathrm{\Delta }T=\mathrm{512}$. U has been non-dimensionalized with $\sqrt{gH}$, where $H={h}_{\mathrm{u}}+{h}_{\mathrm{l}}$.

Figure 6Growth of amplitudes of a wave on the ℐ𝒢+ branch which resonates a wave on the 𝒮𝒢+ branch. Parameters are the same as in Fig. 5. The subscript “i” indicates that the amplitudes are measured at the pycnocline.

Further, in cases where the density ratio limits towards one (R→1), there cannot exist a gravity wave at the interface. However, if there is a shear jump, then the interface may support a vorticity wave. Thus, there still may be energy transfer from the surface to the vorticity interface through Bragg resonance. The other way such an energy transfer can happen is through a wave–triad interaction, which was the theme of the paper by . For wave–triad interaction involving a vorticity wave, no such explosive instability was reported by , but for Bragg resonance, we find that the resonant interactions involving the vorticity wave are prone to explosive instability if shear is in the bottom layer. For a piecewise linear velocity profile with ${U}_{\mathrm{u}}={U}_{\mathrm{l}}=U>\mathrm{0},{U}_{\mathrm{b}}=\mathrm{0}$ for which ${\mathrm{\Omega }}_{\mathrm{u}}=\mathrm{0},{\mathrm{\Omega }}_{\mathrm{l}}=U/{h}_{\mathrm{l}}$, the dispersion relation is given by

$\begin{array}{}\text{(20)}& {\mathit{\omega }}_{\mathrm{in}}^{\mathrm{3}}+\frac{{\mathrm{\Omega }}_{\mathrm{l}}{T}_{\mathrm{l}}}{\left(\mathrm{1}+{T}_{\mathrm{u}}{T}_{\mathrm{l}}\right)}{\mathit{\omega }}_{\mathrm{in}}^{\mathrm{2}}-\frac{gk\left({T}_{\mathrm{u}}+{T}_{\mathrm{l}}\right)}{\left(\mathrm{1}+{T}_{\mathrm{u}}{T}_{\mathrm{l}}\right)}{\mathit{\omega }}_{\mathrm{in}}-\frac{gk{\mathrm{\Omega }}_{\mathrm{l}}{T}_{\mathrm{u}}{T}_{\mathrm{l}}}{\left(\mathrm{1}+{T}_{\mathrm{u}}{T}_{\mathrm{l}}\right)}=\mathrm{0},\end{array}$

where Tu=tanhkhu and Tl=tanhkhl.

It can be easily seen from the above Eq. (20), which is a cubic equation in ωin for Ωl>0, that the product of the three roots is positive. Further, we know that one gravity wave branch has a positive intrinsic frequency, while the other has negative intrinsic frequency. This basically means that the vorticity wave branch must have a negative intrinsic frequency. However, for any value of a positive velocity U, the frequency of vorticity waves becomes positive. The dispersion relation showing the negative energy branch (ℐ𝒱) has been plotted in Fig. 4b.

In order to better visualize explosive resonance, we have also simulated an explosive case in a two-layered flow. The incident wave is a relatively short interfacial mode with a wavenumber k1H=2.1, whereas the resonant wave is a surface mode with the wavenumber k2H=0.1. The bottom ripple has a wavenumber kbH=2.0, the velocity at the bottom is zero and the velocity at both the surface and the interface is $U/\sqrt{gH}=\mathrm{0.1565}$. Figure 5a shows the wave packet at t=0 and Fig. 5b is the wave packet at the final time. The resonant wave with a large wavelength can also be observed.

Further, we have also plotted the variation of Fourier transform of the interface with time in Fig. 6. Increase in both the amplitudes pertaining to k=k1 and k=k2 shows “explosive” growth.

4 Conclusions

Bragg resonance has been traditionally understood in the absence of any background velocity field, and in such scenarios, the amplitude of one wave decays to cause a growth in the amplitude of the other wave. However, in the presence of a velocity field, we show that it is possible to have an exponential growth in the amplitudes of the waves. Notably, even though the presence of a velocity field is necessary, this exponential growth is not a consequence of linear instability because the velocity field chosen is not linearly unstable to perturbations. Further, this simultaneous growth happens without any violation of the law of conservation of energy similar to the explosive instability arising due to wave–triad interaction .

We have explored the possibilities of explosive triads (i.e., all the waves involved in the system grow while keeping the energy conserved) where one of the involved waves is the bottom ripple. Although we have shown it for certain velocity profiles, the fundamental reason is the Doppler shift of the waves with respect to the bottom. Unlike wave–triad interaction, where involvement of a negative energy wave may lead to explosive growth, in the case of Bragg resonance, the presence of such a wave will lead to explosive growth. For a single-layered flow, the velocity required for the formation of an explosive triad is moderately high; however, when a pycnocline is present, explosive instabilities can occur even for small velocities. This is because of low intrinsic frequency of the interfacial gravity wave due to which even a low velocity can Doppler shift the intrinsic frequency to change the sign of the observed frequency.

Data availability
Data availability.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements
Acknowledgements.

Anirban Guha thanks the Alexander von Humboldt Foundation for funding support.

Financial support
Financial support.

The article processing charges for this open-access publication were covered by a Research Centre of the Helmholtz Association.

Review statement
Review statement.

This paper was edited by Roger Grimshaw and reviewed by two anonymous referees.

References

Alam, M. R., Liu, Y., and Yue, D. K. P.: Bragg resonance of waves in a two-layer fluid propagating over bottom ripples, Part I, Perturbation analysis, J. Fluid Mech., 624, 191–224, 2009a. a

Alam, M. R., Liu, Y., and Yue, D. K. P.: Bragg resonance of waves in a two-layer fluid propagating over bottom ripples, Part II, Numerical simulation, J. Fluid Mech., 624, 225–253, 2009b. a

Alford, M. H., Mickett, J. B., Zhang, S., MacCready, P., Zhao, Z., and Newton, J.: Internal Waves on the Washington Continental Shelf, Oceanography, 25, 66–79, https://doi.org/10.5670/oceanog.2012.43, 2012. a

Ball, F. K.: Energy transfer between external and internal gravity waves, J. Fluid Mech., 19, 465–478, 1964. a

Cairns, R. A.: The role of negative energy waves in some instabilities of parallel flows, J. Fluid Mech., 92, 1–14, 1979. a, b, c, d

Craik, A. D. D. and Adam, J. A.: “Explosive” resonant wave interactions in a three-layer fluid flow, J. Fluid Mech., 92, 15–33, 1979. a, b, c, d

Davies, A. G.: The reflection of wave energy by undulations on the seabed, Dynam. Atmos. Ocean, 6, 207–232, 1982. a, b

Drivas, T. and Wunsch, S.: Triad resonance between gravity and vorticity waves in vertical shear, Ocean Model., 103, 87–97, 2016. a, b

Elgar, S., Raubenheimer, B., and Herbers, T. H. C.: Bragg reflection of ocean waves from sandbars, Geophys. Res. Lett., 30, 16–1, https://doi.org/10.1029/2002GL016351, 2003. a

Grimshaw, R.: Wave action and wave-mean flow interaction, with application to stratified shear flows, Annu. Rev. Fluid Mech., 16, 11–44, 1984. a

Heathershaw, A. D. and Davies, A. G.: Resonant wave reflection by transverse bedforms and its relation to beaches and offshore bars, Mar. Geol., 62, 321–338, 1985. a

Kirby, J. T.: A general wave equation for waves over rippled beds, J. Fluid Mech., 162, 171–186, 1986. a

Mei, C. C.: Resonant reflection of surface water waves by periodic sandbars, J. Fluid Mech., 152, 315–335, 1985. a

Pietrzak, J. D., Kranenburg, C., and Abraham, G.: Resonant internal waves in fluid flow, Nature, 344, 844, https://doi.org/10.1038/344844a0, 1990. a

Raj, R. and Guha, A.: On Bragg resonances and wave triad interactions in two-layered shear flows, J. Fluid Mech., 867, 482–515, 2019. a, b, c, d, e

Voronovich, A. G., Lobanov, E. D., and Rybak, S. A.: On the stability of gravitational-capillary waves in the presence of a vertically nonuniform current, Izv. Atmos. Ocean Phys., 16, 220–222, 1980. a

Note that in this case we have taken Ω>0; had we taken Ω<0, we would have got a maximum in the 𝒮𝒢+ curve rather than getting a minimum in the 𝒮𝒢 curve. Basically, the dispersion curve in that case would be a mirror image of Fig. 2 about the line ω=0.