## On Hamiltonian systems integrable in elliptic functions that describe waves over underwater banks and ridges

Yu.V.Brezhnev, A.V.Tsvetkova

Abstract . We discuss the 4-dimensional Hamiltonian systems that describe waves over underwater banks and ridges. The systems are exactly integrable in terms of elliptic functions and of solutions to nontrivial transcendental equations involving the elliptic integrals (Weierstrass' ζ -function).

## Contents

## 1 Exact analytic solutions

## 2 Visualizing the fronts. An algorithm

## 3 Proof of theorem 1 in a nutshell

2

4

5

## Introduction

In this paper, we discuss the problem of integration of the following Hamiltonian system with two degrees of freedom:

$$\begin{array} {$$

where C x ( ) is a positive smooth function.

We show that such systems are integrable in elliptic functions for functions D x ( ) of a certain form; we shall write them out further below. The problem of integration is interesting not only as an independent problem, but also from the standpoint of applications. In particular, such Hamiltonian systems arise in the study of rapidly varying solutions of the two-dimensional wave equation. For example, in the case of long linear gravity surface waves in an infinite basin with varying bottom given by D x ( ) &gt; 0, the function C x ( ) = √ gD x ( ) ( g is the free-fall acceleration), and the corresponding equation has the form [1, 2]

$$\frac { \partial ^ { 2 } u } { \partial t ^ { 2 } } - \$$

The initial (or some other additional) conditions for system (1) are generated by the conditions for equation (2). We consider the initial conditions for system (1) corresponding to the Cauchy problem for (2) with localized initial conditions

$$u | _ { t = t _ { o } } = V \left ( \frac { x -$$

∣ where V ( y ) is a smooth function decreasing rapidly enough as | y | → ∞ , the parameter l characterizes the source size, and x ◦ is a point in whose neighborhood the initial perturbation is localized. Such conditions for (2) give the following Cauchy problem [3, 4, 5] for system (1):

$$p | _ { t _ { o } } = \mathbf n ( \psi ) = \$$

The Hamiltonian system we consider makes it possible to describe the water waves.

By X ψ,t , P ( ) ( ψ, t ) we denote the solutions of problem (1), (4). At each moment of time t , the ends of these trajectories determine the following smooth closed curves

$$\Gamma _ { t } = \{ x = X ( \psi, t ), \$$

in the 4-D phase space R 4 xp , which are called wave fronts in the phase space . According to [3, 4, 5], at each t , the asymptotic solution of problem (2), (3) is determined by Γ t and turns out to be localized in a neighborhood of the curves γ t = { x = X ψ,t , ψ ( ) ∈ S 1 } -the projections of Γ t on R 4 xp -which are called fronts in the configuration (physical) space . In contrast to Γ , the t curves γ t may be non-smooth and have the self-intersection points. Under t = 0, the set γ t is the point x ◦ .

As mentioned above, we here we will show that, in the cases where the bottom is shaped as an underwater bank, i.e., the depth of the water layer is determined by the function

$$D ( \varrho, \varphi ) = D ( \varrho ) \equiv \frac { b + \varrho ^ { 2 } } { a + \varrho ^ { 2 } }, \quad a > b > 0 \text{ are constants, } \ \varrho \text{ is the polar radius, } \quad ( 5 )$$

and in the case where the underwater ridge-shape is determined by the function

$$D ( x _ { 1 }, x _ { 2 } ) = \frac { x _ { 1 } ^ { 2 } + b } { x _ { 1 } ^ { 2 } + a }, \quad a > b > 0 \text{ are constants} \,.$$

The above-written Hamiltonian system is integrable in elliptic functions. The exact analytic solutions obtained in these cases are given in sect. 1. These formulas can be used to construct an asymptotic solution of problem (2), (3) and asymptotic solutions of stationary equations with localized right-hand sides, for example, the Helmholtz equation [6]

$$- \omega ^ { 2 } v - \varepsilon ^ { 2 } \langle \nabla, C ^ { 2 } ( x ) \nabla \rangle v = \frac { 1 } { \varepsilon } V \left ( \frac { x - x ^ { 0 } } { \varepsilon } \right ), \quad \omega = \text{const}, \ \ x \in \mathbb { R } ^ { 2 }.$$

The corresponding formulas will be given in the expanded version of the work. In sect. 2, we describe an algorithm for constructing fronts, which makes it possible to visualize them and facilitates dealing with the obtained expressions. In sect. 3, we outline the proofs.

## 1. Exact analytic solutions

Since the function determining the bottom shape is symmetric in both cases under study, we can without loss of generality assume that x ◦ = -ξ 0 , where ξ &gt; 0.

( ) Let us consider the situation where the function describing the depth of the water layer has the form (5), i.e., determines the underwater bank. We assume that the free-fall acceleration g = 1. Then, in the polar ( /rho1, ϕ )-coordinates { x 1 = /rho1 cos ϕ, x 2 = /rho1 sin ϕ } , the Hamiltonian becomes

$$\mathcal { H } ( \varrho, \varphi ; u, v ) = \sqrt { u ^ { 2 } + \frac { v ^ { 2 } } { \varrho ^ { 2 } } } \cdot \sqrt { \frac { \varrho ^ { 2 } + b } { \varrho ^ { 2 } + a } },$$

where u := p /rho1 , v := p ϕ are the momenta corresponding to the system ( /rho1, ϕ . ) The initial conditions (4) in these coordinates become

$$u | _ { t _ { o } } = - \cos \psi, \quad v | _ { t _ { o } } = - \xi \sin \psi, \quad \varrho | _ { t _ { o } } = \xi, \quad \varphi | _ { t _ { o } } = \pi.$$

Theorem 1. The solution of the Hamiltonian system with Hamiltonian (7) and initial conditions (8) is given by the expressions

$$\text{ations (\delta ) is given by the expressions} \\ \left \{ \begin{array} {$$

$$\wp ( \mathfrak { p } _ { \circ } ) = \xi ^ { 2 } - \delta + a, \quad \wp ( \varkappa ) = \frac$$

The expressions for the constants δ, α , and β in terms of the parameters of the problem have the form where

$$\L _ { \L } o r m & \delta & = \frac { 1 } { 3 } \left ( \mathfrak { h } + 2 \, a - b \right$$

$$h ^ { 2 } = \mathfrak { h } = \xi ^ { 2 } \, \frac { \xi ^ { 2 } + a } { \xi ^ { 2 } + b$$

Here, as always in the sequel, the notation ℘ z ( ) = ℘ z ( ; α, β ) , ζ ( z ) = ζ ( z ; α, β ) , σ z ( ) = σ z ( ; α, β ) is used for the Weierstrass elliptic functions [7] .

Remark 1. Note that p ◦ , p , and κ need not be real. They lie on the edges of the parallelogram (0 , ω, ω + ω , ω ′ ′ ), where ω is the pure real period of the Weierstrass ℘ -function and ω is the pure imaginary one. In this case, p lies on the same edge of the parallelogram as p ◦ . We also note that equations (10) have infinitely many roots. For definiteness, we assume that the roots are taken from the first positive branch. In sect. 2, we describe an algorithm for constructing fronts, which illustrates these formulas.

In the case where the function determining the depth of the water layer has the form (6) and describes the underwater ridge, the solutions are given by the following theorem.

Theorem 2. The solution of the Hamiltonian system (1) with Hamiltonian

$$\mathcal { H } ( x, y ; u, v ) & = \sqrt { p _ { x _ { 1 } } ^ { 2 } + p _ { x _ { 2 }$$

and initial conditions (4) has the form

$$\begin{cases} \quad x _ { 1 } = - \sqrt { \varphi ( \mathfrak { p } ; \alpha, \beta ) + \delta - a }, \quad x _ { 2 } = \sqrt { 1 - \mathfrak { h } ^ { - 1 } t } + \sqrt { \mathfrak { h } - 1 } \left ( b - a \right ) \cdot ( \mathfrak { p } - \mathfrak { p } _ { 0 } ), \\ p _ { x _ { 1 } } = \frac { 1 } { 2 } \, \frac { \gamma } { \sqrt { \mathfrak { h } - 1 } } \frac { \varphi ^ { \prime } \left ( \mathfrak { p } ( t ) \right ) } { \varphi ( \mathfrak { p } ( t ) ) + \delta - a + b } \, - 2 \sqrt { \varphi ( \mathfrak { p } ( t ) ) + \delta - a }, \quad p _ { x _ { 2 } } = \gamma. \end{cases} \\ \text{The variables $\mathfrak { p }$ and $\mathfrak { p } _ { 0 }$ are solutions of the transcendental equations} \\ \cdot \end{cases}$$

$$\frac { 1 } { h } \, t = \delta \cdot ( \mathfrak { p } - \mathfrak { p } _ { \circ } ) - \zeta ( \mathfrak { p } ; \alpha, \beta ) + \zeta ( \mathfrak { p } _ { \circ } ; \alpha, \beta ), \quad \mathfrak { p } _ { \circ } = \wp ^ { - 1 } ( \xi ^ { 2 } - \delta + a ; \alpha, \beta ). \quad ( 1 4 )$$

The expressions for the constants δ, α , and β in terms of the parameters of the problem have the form where

$$\text{for} m & & \delta = \frac { 1 } { 3 } \left ( ( b - a ) \, \mathfrak { h } + 3 \, a - 2 \, b \right ), \quad \alpha = \frac { 4 } { 3 } \left ( ( b - a ) ^ { 2 } \mathfrak { h } ^ { 2 } - ( b - a ) \, b \, \mathfrak { h } + b ^ { 2 } \right ), \\ & & \beta = \frac { 4 } { 2 7 } \left ( ( b - a ) \, \mathfrak { h } + b \right ) \left ( ( b - a ) \, \mathfrak { h } - 2 \, b \right ) \left ( 2 \, ( b - a ) \, \mathfrak { h } - b \right ), \\ \text{where} & & \varepsilon ^ { 2 } \, \pm \, h$$

$$h ^ { 2 } = \mathfrak { h } = \frac { \xi ^ { 2 } + b } { ( \xi ^ { 2 } + b ) - ( \xi ^ { 2 } + a ) \sin ^ { 2 } \psi }, \quad \gamma = \sin \psi.$$

## 2. Visualizing the fronts. An algorithm

Here we describe an algorithm for constructing, in particular, the fronts γ t by using the software package Wolfram Mathematica . This algorithm makes it easy to visualize the formulas obtained in the preceding section and illustrate the application of these formulas. On the other hand, the problem of visualizing fronts is also interesting from the standpoint of applications, because, as mentioned above, the asymptotic solution of problem (2), (3) is localized in a neighborhood of the front [3, 4, 5].

Figure 1 . The graph of y = ℘ x ( ) -ξ 2 + δ -a under ψ = 9 10 π

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00938v1-figure-317ec98e-7e58-4eb1-8a2b-d1ca027c793b.jpg)

## Algorithm

Step -1: We fix the following parameters of the problem: the constants a = 100, b = 1 determining the bottom shape and the constant ξ = 0 5 characterizing the source position. .

Step -2: For each fixed ψ ∈ [0 , π ], there exist infinitely many roots of the equation

$$\wp ( p _ { \circ } ) = \xi ^ { 2 } - \delta + a.$$

As p ◦ 1 ( ψ ) and p ◦ 2 ( ψ ) we take, respectively, the first and second positive roots (see Fig. 1). In the figure, the red point corresponds to p ◦ 1 and the blue point, to p ◦ 2 . We determine κ ( ψ ) by the formula (10), where we again take the first positive root.

Figure 2 . The graph of the front at t = 100

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00938v1-figure-83597cce-a3ff-4f54-a5bf-95331e88d76b.jpg)

Remark 2. As already mentioned, the roots of (17) need not be real. We must seek the complex roots in one of the forms ip ◦ , p ◦ + ω ′ , and ip ◦ + ω , where ω is the pure real period of the Weierstrass ℘ and ω ′ is its pure imaginary period, p ◦ is real; as p ◦ 1 and p ◦ 2 we must again take the first two positive roots.

Step -3: We fix a moment of time t . Let p i ( ψ, t ) ( i = 1 2) be the first positive root of the , equation t = δ · [ p -p ◦ i ( ψ )) -ζ ( p ) + ζ ( p ◦ i ( ψ )]. If p ◦ i is not real, then we seek p i in the same manner as p ◦ i .

Step -4: For each fixed t , we plot the curves in polar coordinates /rho1, ϕ , using formulas (9). In the case ψ ∈ [ 0 , π 2 ] , we take p ◦ 1 and p 1 as parameters, and in the case ψ ∈ [ π 2 , π ] , we take p ◦ 2 and p 2 .

## 3. Proof of theorem 1 in a nutshell

Step -5: By symmetry, we reflect the obtained graph about the horizontal axis. Figure 2 shows the front at t = 100; the dashed curve in this figure corresponds to ψ ∈ [ 0 , π 2 ] and the solid curve, to ψ ∈ [ π 2 , π ] .

Solving the Hamiltonian system (1) with Hamiltonian (7) and initial conditions (8), we get

$$v = \gamma \ \left ( = \text{const} \right ), \quad \left ( u ^ { 2 } + \frac { \gamma ^ { 2 } } { \varrho ^ { 2 } } \right ) \frac { \varrho ^ { 2 } + b } { \varrho ^ { 2 } + a } = \frac { \gamma ^ { 2 } } { h ^ { 2 } } \ \left ( = \text{const} \right ).$$

Making use of the integral H ( /rho1, ϕ u, v ; ) = γ/h , we obtain the following autonomous dynamics for the variable /rho1

$$\varrho \, \dot { \varrho } & = \frac { \sqrt { ( \varrho ^ { 2 } + b ) \left [ \varrho ^ { 2 } ( \varrho ^ { 2 } + a ) - ( \varrho ^ { 2 } + b ) h ^ { 2 } \right ] } } { ( \varrho ^ { 2 } + a ) } \,. \\ \mathfrak { r } \,. & \,. & \,. & \,. & \,. & \,. & \,. & \,.$$

Making the change z = /rho1 2 + a and applying the shift s = z -δ to reduce the integral to the canonical Weierstrass form we derive an expression for /rho1 in system (9). Using the obtained result and the integral H = γ/h , we obtain the formula for u . Substituting the obtained expression into the Hamiltonian system,

$$\int \frac { ( s + \delta ) \, d s } { \int ^ { \frac { ( s + \delta ) \, d s } { \sqrt { 4 \, s ^ { 3 } - \alpha \, s - \beta } } } } = t, \\ \intertext { system (9). } \text{ using the obtained res}.. \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \end{$$

one obtains the equation

$$\dot { \varphi } = h \left ( 1 + \frac { b } { \wp ( \mathfrak { p } ) + \delta - a } \right ) \frac { 1 } { \wp ( \mathfrak { p } ) + \delta } \,. \\. \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot$$

Its integration leads to the logarithmic elliptic integral that can be easily evaluated [7]. Simplifying the obtained expression, we arrive at the solutions (9).

## Acknowledgments

The authors wish to express their gratitude to S. Yu. Dobrokhotov for the statement of the problem and to M. Babich and M. Pavlov for useful discussions.

## References

- [1] J. J. Stoker. Waves on Water. The Mathematical Theory with Applications . New York-London, Interscience (1957).
- [2] E. N. Pelinovskii. Hydrodynamics of Tsunami Waves. Nizhnii Novgorod (1996) (in Russian).
- [3] S. Yu. Dobrokhotov, S.Ya.Sekerzh-Zenkovich, B. Tirozzi &amp; B. Volkov. Explicit asymptotics for tsunami waves in framework of the piston model. Russ. J. Earth Sciences (2006) 8 (ES403), 1-12.
- [4] S. Yu. Dobrokhotov, A.I. Shafarevich &amp; B.Tirozzi. Localized wave and vortical solutions to linear hyperbolic systems and their application to linear shallow water equations. Russ. J. Math. Phys. (2008) 15 (2), 192-221.
- [5] S. Yu. Dobrokhotov &amp; V.E.Nazaikinskii. Asymptotics of localized wave and vortex solutions of a linearized system of shallow water equations. In: Actual Problems of Mechanics (2015), 98-139. Nauka, Moscow (in Russian).
- [6] A.Yu.Anikin, S. Yu.Dobrokhotov, V. E.Nazaikinskii &amp; M. Rouleux. The Maslov canonical operator on a pair of Lagrangian manifolds and asymptotic solutions of stationary with localized right-hand side. Doklady Math. (2017) 96 (1), 406-410.
- [7] N.I. Akhiezer. Elements of the Theory of Elliptic Functions. Transl. Math. Monogr. 79 . AMS Providence, R. I. (1990).

Yu.V.Brezhnev

Tomsk State University

Tomsk, 634050 Russia

E-mail :

brezhnev@mail.ru

A.V.Tsvetkova Ishlinsky Institute for Problems in Mechanics RAS Moscow, 119526 Russia Moscow Institute of Physics and Technology Dolgoprudnyi, Moscow Region, 141701 Russia E-mail : annatsvetkova25@gmail.com
Yu. Brezhnev A. Tsvetkova
On Hamiltonian systems integrable in elliptic functions that describe waves over underwater banks and ridges
We discuss the 4-dimensional Hamiltonian systems that describe waves over
underwater banks and ridges. The systems are exactly integrable in terms of
elliptic functions and of solutions to nontrivial transcendental equations
involving the elliptic integrals (Weierstrass' $\zeta$-function).