Trapped solitary-wave interaction for Euler equations with low pressure region

04/02/2020 ∙ by Marcelo V. Flamarion, et al. ∙ 0

Trapped solitary-wave interaction is studied under the full Euler equations in the presence of a variable pressure distribution along the free surace. The physical domain is flattened conformally onto a strip and the computations are performed in the canonical domain. Computer simulations display solitary waves that remain trapped in a low pressure region. In terms of confinement we observe that these waves are stable for small perturbations of either their amplitudes or the pressure forcing term. Furthermore multiple solitary waves are considered within the low pressure region without escaping the low pressure region. We identify regimes in which multiple solitary waves remain trapped after several collisions. In particular we display a regime where three solitary waves are trapped and collide several times, before one escapes at a time. The remaining solitary waves stays trapped in the low pressure region.



There are no comments yet.


page 9

This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

1 Introduction

Waves generated by an external force frequently occur in nature or in laboratory experiments. These waves have been extensively studied through reduced models such as the Boussinesq and the forced Korteweg-de Vries (fKdV) (Wu [1987]; Milewski [2004]; Binder et al. [2005]; Flamarion et al. [2019]). Physically these waves can arise in the context of atmospheric flows encountering topographic obstacles or flow of water over rocks (Baines [1995]). In addition, waves can also be generated when one considers a pressure distribution moving over a free surface, as for instance, ship waves and ocean waves generated by storms (when a low pressure region moves on the surface of the ocean) (Johnson [2012]). In a laboratory experiment performed by Pratt [1984], in which two obstacles were fixed and water circulated through a flume with a variable-speed, it was observed a steady wave train between the two obstacles (see Pratt’s figure ). Here we use the term “trapped wave” to describe waves that remain trapped in a certain region, generally in between two topographic obstacles or in low pressure regions. When the trapped waves escape after sufficient time has elapsed we use the terminology temporally trapped waves.

From the fKdV point of view a great deal of attention has been paid on trapped waves over a bottom topography. For the interested reader, we mention a few recent references form which the bibliography may be useful. Grimshaw & Maleewong [2013] showed that steady trapped waves can be found in the vicinity of a holelike forcing (topography) and studied numerically and theoretically how the shape of the hole affects the formation of trapped waves. Considering a bottom topography with two bumps Lee & Whang [2015] and Lee [2018] found steady waves and used them as initial data for the fKdV equation observing that those waves remained trapped between the two bumps up to a certain period of time. They also investigated the stability of those temporally trapped waves for small perturbations in their amplitudes and in the forcing term. In the same spirit Kim & Choi [2017] studied the propagation of temporally trapped waves between two bumps submerged and showed that the trapped waves have to overcome a certain threshold of energy in order to go over the bumps. They concluded that the energy barrier does not depend on the distance between the bumps, while it increases with the height of the bumps.

To the best of our knowledge there are no articles about time-dependent trapped waves solutions for the Euler equations. In particular, we here present the first study regarding waves trapped by a pressure distribution. For steady trapped waves a study was done by Dias & Vanden-Broeck [2004], under the assumption of a two-bump topography. They found steady wave trains in between the two obstacles and showed that its shape depends on the height and width of the bumps. They also showed that as the distance between the obstacles increases the flow approaches the generalised hydraulic fall past a single bump.

In this work we consider the full Euler equations in the presence of a localised low-pressure distribution on the free surface as a proxy for a storm over a region. We found time-dependent trapped solitary waves which move back and forth in the low pressure region). We also investigate the stability of these waves for small perturbations of their amplitudes and the profile of the pressure profile with respect to its amplitude and range. Our numerical experiments indicate that single solitary waves are stable, namely remaining trapped. An additional numerical result is the existence of multiple trapped solitary waves. We have found regimes in which several solitary waves remain trapped in the low pressure region after many collisions.

The paper is organized as follows. In section we present the mathematical formulation of the non-dimensional Euler equations. In section we rewrite the Euler equations in the canonical domain, which is a uniform strip where computations are more easily performed and present the numerical method to solve the equations. The results are presented in section and the conclusion in section .

2 The Euler equations in the presence of a moving pressure

We consider a two-dimensional incompressible and irrotational flow of an inviscid fluid with constant density () in a channel of finite depth () in the presence of gravity () with a traveling pressure distribution with constant speed () on the free surface. We denote the velocity potential by and the free surface profile by . Choosing , , e as our characteristic scales for length, velocity, time and pressure, respectively, the non-dimensional governing equations are as follows:


where is the Froude number. It is convenient to consider the equations in the moving frame and by writing

we obtain


We are interested in investigating the scenario in which the wave is trapped in a low pressure region. According to our adimensionalization means atmospheric pressure, while positive pressure means a high pressure region and negative pressure means a low pressure region.

3 Conformal mapping and numerical methods

In order to deal numerically with the system (2) we proceed in the same fashion as Dyachenko et al. [1996], which roughly speaking consists in constructing a time-dependent conformal mapping

which maps a strip of width onto the fluid domain, satisfying the boundary conditions

The canonical depth is a function of , which is chosen so that the canonical and physical domains have the same length.

Let and where is the harmonic conjugate of . Their traces along is denoted by and . By , we denote the horizontal and vertical free surface coordinates at . Substituting these variables in yields when evaluated


where , , and


Fourier modes are given by

where , . More details about these computations can be seen in Ribeiro-Jr et al. [2017] and Flamarion et al. [2019]. Moreover, the adimensional mass integral quantity in the canonical system is given by the formula


In the following sections we consider the computational domain , with an uniform grid with points and step . All derivatives in are computed spectrally (Trefethen [2001]). Besides, the time evolution of the system (3) is calculated through the Runge-Kutta fourth order method with time step . Unless mentioned otherwise, we use the following parameters: , and .

4 Results

In all simulations presented in this work we consider as initial data for system (3) unforced traveling wave solutions with velocity , namely computed in the absence of pressure () and Froude number . The initial data is obtained in the same fashion as Ribeiro-Jr et al. [2017] by a Newton continuation method.

We verify our numerical method by first evolving the initial data considering and in the system (3). The method is accurate in the sense that the solution is indeed a steady wave. Furthermore, denoting the initial wave profile by and by

the absolute error we get the relative error

Now, we study time-dependent trapped solitary waves in low pressure regions. For this purpose we assume that in system (2) the pressure term on the free surface is modelled by the function

with its profile is depicted in figure 1. Unless mentioned otherwise, we assume , and .

Figure 1: Profile of the pressure and its parameters.

Figure 2 (top) shows the evolution of a trapped solitary wave in the low pressure region (darker region) as well as the dynamics of its crest (thicker black line). Initially the solitary wave crest is at . As time elapses the wave travels downstream until it reaches the high pressure region, then it reflects back and moves upstream until it reaches the high pressure region again. This motion is repeated several times. Figure 2 (bottom) displays how the amplitude of the trapped wave varies throughout time. Notice that when the wave is moving downstream its amplitude is smaller than when it is moving upstream. Besides it is remarkable that the wave behaves effectively as a travelling wave. No radiation is observed during reflection and amplitude adjustments. The same will be observed during collisions as reported below. In this simulation, the mass is conserved (3). More precisely, we have that

Figure 2: Top: Trapped solitary wave in the low pressure region (darker area). The thicker black line shows the evolution of the crest of the wave. Bottom: The amplitude of the wave as a function of time.

Some authors (Lee & Whang [2015]; Kim & Choi [2017]; Lee [2018]) have studied the stability of trapped waves for the Forced-Korteweg-de-Vries equation by disturbing its amplitude and its forcing term examining whether the waves remain trapped or not. Along the same lines we disturb the amplitude of the wave showed on figure 2.

Firstly, we consider the initial data under a perturbation with a reduction of its amplitude. The results are shown in figure 3. The wave that was previously trapped, now escapes. In our numerical simulations we found that for reductions smaller than the wave never reverses its movement, i.e, it always moves downstream.

Figure 3: Top: The figure shows the evolution of the trapped solitary wave with reduction of its amplitude. The thicker line displays the propagation of its crest. Bottom: The amplitude of the wave as a function of time.

Secondly, we increase the amplitude of the trapped wave in . This is depicted in figure 4. As we can see the perturbation is strong enough to make the wave move upstream at and leave the low pressure region at . We would like to point out that the wave remains trapped (at least until ) for perturbations smaller than . To the best of our knowledge there is no theory to explain the nontrivial behaviour of the temporally trapped solitary wave.

Figure 4: Top: The figure shows the evolution of the trapped solitary wave with accretion of its amplitude. The thicker line represents its crest propagation. Bottom: The amplitude of the wave as a function of time.

We also consider perturbations in the amplitude of the pressure . In our simulations we observed that the trapped wave is robust for reductions below . We carried out accretions of up to and the wave remained confined in the low pressure region (at least until ). For accretions above the generated waves due to the pressure may break (see Grimshaw & Maleewong [2013]) which drives this analysis to another subject that is not within the scope of our study. Figure 5 displays the propagation of the crest of some waves for different values of (top) and the heights their peaks (bottom). On one hand as increases so does the back-and-forth motion of the trapped waves. On the other hand the fluctuations of the peaks do not change much as time elapses. Note that the smallest value of delta is too weak to promote the oscillator pattern of the solitary wave.

Figure 5: Top: Propagation of the crest for different values of . Bottom: The heights of the highest peaks as function of time.

Regarding of the parameter , the initial wave remains trapped in the low pressure region up to for values of (which is approximately half the efective wavelength of the initial wave).

Lastly we consider multiple solitary waves as initial data for (3) in a larger low pressure region. For this purpose we set the initial data as sums of translations of . Figure 6 (left) displays two solitary waves propagating downstream at . When the left solitary wave reaches the high pressure region it reflects back upstream and collides with the other solitary wave and both reverse directions. Then this behaviour repeats on the right side of the high pressure region and we observe that after many collisions both solitary waves remain trapped in the low pressure region. Figure 6 (right) depicts three solitary waves. Initially all waves propagate downstream. Differently of the previous case after a few collisions one solitary wave escapes the low-pressure region. After another shorter series of collisions a secondary solitary wave escapes with a remaining trapped wave. To the best of our knowledge the scenarios here presented have never been contemplated in the literature, not even for the fKdV equation.

Figure 6: The figure shows the collision of solitary waves in the low pressure region. Left: Two solitary waves trapped in the low pressure area. Right: The propagation of three solitary waves. As time elapses two of the waves escapes and only one remains trapped. Parameters: and .

5 Conclusion

In this paper we have studied time-dependent trapped solitary waves solutions for the full Euler equations. Through an iterative method we computed we computed solitary waves for the Euler equations and showed that for certain regimes these waves move back and forth oscillating within the low pressure region. Our numerical investigations demonstrated that the trapped wave solutions are stable for small perturbations of either their amplitudes or the forcing term. Furthermore we discovered regimes in which multiple solitary waves can temporally be trapped in a low pressure region. The numerical method applied in this work can be extended to the case in which there is a variable topography. This is a first step in studying time-dependent dynamics for trapped waves under different types of forcing and different flow regimes. New solitary wave-interaction regimes were uncovered, with further investigations to be pursued in the future.

6 Acknowledgements

The authors are grateful to IMPA-National Institute of Pure and Applied Mathematics for the research support provided during the Summer Program of 2020 and to Prof. André Nachbin (IMPA) for his constructive comments and suggestions which improved the manuscript. M.F. is grateful to Federal University of Paraná for the visit to the Department of Mathematical Sciences. R.R.-Jr is grateful to University of Bath for the extended visit to the Department of Mathematical Sciences.


  • Baines [1995] Baines, P. 1995 Topographic effects in stratified flows. Cambridge University Press.
  • Binder et al. [2005] Binder, BJ., Dias. F., Vanden-Broeck, J.-M. 2005 Forced solitary waves and fronts past submerged obstacles, Chaos.
  • Dias & Vanden-Broeck [2004] Dias. F., Vanden-Broeck, J.-M. 2004 Trapped waves between submerged ostacles, J. Fluid Mech. 509, 93-102.
  • Dyachenko et al. [1996] Dyachenko. AL., Zakharov. VE., Kuznetsov. EA. 1996 Nonlinear dynamics of the free surface of an ideal fluid, Plasma Phys. 22, 916-928.
  • Flamarion et al. [2019] Flamarion, M. V. & Milewski, P. A. & Nachbin A. 2019 Rotational waves generated by current-topography interaction. Stud Appl Math, 142, 433-464.
  • Grimshaw & Maleewong [2013] Grimshaw, R. & Maleewong, M. 2013 Stability of steady gravity waves generated by a movingl ocalized pressure disturbance in water of finite depth. Phys Fluids, 25.
  • Johnson [2012] Johnson, R. S. 2012 Models for the formation of a critical layer in water wave propagation. Phil. Trans. R. Soc. A 370, 1638-1660.
  • Lee & Whang [2015] Lee, S., Whang, S. 2015 Trapped supercritical waves for the forced KdV equation with two bumps. Appl. Math. Model. 39, 2649-2660.
  • Lee [2018] Lee, S. 2018 Dynamics of Trapped Solitary Waves for the Forced KdV Equation. Symmetry. 10, (129).
  • Kim & Choi [2017] Kim, H., Choi, H. 2017 A study of wave trapping between two obstacles in the forced Korteweg-de Vries equation J. Eng. Math
  • Milewski [2004] Milewski, P. A. 2004 The Forced Korteweg-de Vries Equation as a Model for Waves Generated by Topography. CUBO A mathematical Journal 6 (4), 33-51.
  • Pratt [1984] Pratt, LJ. 1984 On nonlinear flow with multiple obstructions. J. Atmos. Sci. 41 , 1214-1225.
  • Ribeiro-Jr et al. [2017] Ribeiro Jr R., Milewski P.A., & Nachbin A. 2017 Flow structure beneath rotational water waves with stagnation points. J. Fluid Mech. 812, 792-814.
  • Trefethen [2001] Trefethen, L. N. 2001 Spectral Methods in MATLAB. Philadelphia: SIAM..
  • Wu [1987] Wu, TY. 1987 Generation of upstream advancing solitons by moving disturbances J Fluid Mech 184 , 75-99.