Abstract
Relying on conformal mappings we prove the logarithmic convexity of certain flow quantities associated with irrotational periodic travelling waves that propagate at the surface of water over a flat bed. These results enable us to quantify the observation that the kinetic energy and the time-period of the particle paths are larger near the surface and reduce with increasing depth.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
Regular ocean waves propagating practically without change of form at constant speed in a fixed direction are commonly observed. Their sharp crests and flat troughs, departing from sinusoidal profiles, are indicative of the nonlinear processes at play. Considerable theoretical progress in the understanding of nonlinear flow properties was achieved in the last decades by relying on an interplay between harmonic and complex analysis for irrotational two-dimensional waves at the surface of a fluid of infinite depth (see [1, 15]), setting that is considered adequate to model deep-water waves. In recent years the flat-bed case received a lot of attention (see the discussions in [3, 6, 11, 14]), and complex-analytic considerations were central to these developments. The flat-bed hypothesis is physically realistic in many circumstances: other than arctic continental shelfs with a layer of relatively shallow water that cover thousands of km\(^2\), there are even vaster deep abyssal plains in some major sea and ocean basins, featuring depth variations of less than 0.01%.
Note that water waves are fluctuations of the free surface that signal tremendous energy transport (the associated mass transport being of a much lesser degree). For this reason, a better understanding of the properties of the kinetic energy associated to these wave phenomena is of great importance (see Sect. 5 and the discussions in [5, 9] for further details). It is natural to start with the most regular wave patterns—two-dimensional periodic travelling waves in irrotational flow over a flat bed, for which complex analysis techniques prove to be very powerful in revealing aspects of the flow structure. We investigate the total kinetic energy of a particle and, in the process of doing this, we also obtain some by-products that are of independent interest. In particular, we clarify an issue relevant to the particle-path pattern. While a good theoretical understanding of the particle paths beneath such waves is already available, supported by numerical simulations and experimental results (see the discussions in [4], respectively [3, 10, 11, 16]), the uniformity in time of the particle paths, that is noticeable numerically and experimentally, has so far not been established theoretically. More precisely, each particle describes a repeated pattern. We show that the time it takes a particle to describe one of these repeating orbits is independent of the initial location of the particle, in the sense that it depends only on the corresponding streamline, thus confirming the experimental evidence (and showing the uniformity of the “elapsed time” concept introduced in [4] as the time needed for a particle initially located beneath the wave trough to regain a position beneath the wave trough). From now on we are going to refer to this specific time as to the streamline time-period. Thus, the maximal elevation of each particle is periodic in time with period equal to streamline time-period, while the horizontal position is shifted (the shift being dependent on the streamline) in the direction of wave propagation after each streamline time-period. An important step in our approach is to notice that the streamline time-period can be expressed as an integral over a horizontal segment of the modulus of a suitably defined analytic function, a feature which permits us to apply Hardy’s convexity theorem for an annulus. The link to such integrals of analytic functions helps us further in the study of the total kinetic energy of a particle over a streamline time-period.
The paper is organized as follows. In Sect. 2 we present the governing equations and a transformation to an equivalent form that better suits our purposes. Section 3 is devoted to the study of the streamline time-period of a particle and in Sect. 4 we derive some results about the total kinetic energy of a particle over a streamline time-period. Finally, some background material is provided in Sect. 5, putting the mathematical considerations made in the paper in a broader context of water-wave energy studies.
2 Preliminaries
For two-dimensional water waves it suffices to investigate the flow characteristics in a cross-section oriented towards the direction of wave propagation. We choose Cartesian coordinates (X, Y) with the X-axis pointing in the direction of wave propagation and the Y-axis oriented upwards. Let \(Y=-d\) be the flat bed, \(Y=\eta (X,t)\) be the free surface and \((U(X,Y,t),\,V(X,Y,t))\) the velocity field. Under the physically reasonable assumption of a homogeneous inviscid flow, the governing equations for a flow determined by the balance between the restoring gravity force and the inertia of the system are the equation of mass conservation
and Euler’s equations
where P(X, Y, t) is the pressure and g is the constant gravitational acceleration. The associated boundary conditions are the kinematic boundary conditions
expressing the fact that the water’s free surface and the flat bed are interfaces, and, the effect of surface tension being negligible, the dynamic boundary condition
decoupling the motion of the water from that of the air above it, with \(P_{atm}\) the constant atmospheric pressure at sea level. The absence of non-uniform currents is ensured by the irrotational character of the flow, requiring
A Stokes wave is a smooth travelling wave solution to the governing equations (1)-(6) for which there exists a period \(\lambda >0\) and a wave speed \(c>0\) such that the free surface profile \(\eta \), the fluid velocity \((U,\,V)\) and the pressure P have period \(\lambda \) in the X variable, \(\eta \) depends only on \((X-cT)\), while U, V, and P depend only on \((X-cT)\) and Y. Moreover, the wave profile is strictly monotonic between successive crests and troughs (in particular, there is a single crest and trough per period), and symmetric about the wave crest. Since we only consider smooth waves, the flow presents no stagnation points, that is,
It is known (see the discussion in [14]) that for irrotational travelling waves the breakdown of (7) can only occur at the wave crest, in which case the free surface is not a continuously differentiable curve, being real-analytic except at the crest where it is continuous with a corner containing an inner angle of 120\(^\circ \) (the setting for incipient wave breaking). In contrast to this, if (7) holds, then a global bifurcation approach ensures the existence of waves of small, moderate and large amplitude, whose symmetric profiles must be real-analytic and for which the velocity components have harmonic extensions across the free surface and the flat bed (see the discussion in [14]).
Taking advantage of the (X, t)-dependence of the form \((X-ct)\), passing to the moving frame
and writing
we can reformulate the governing equations for Stokes waves in the form
with
Without loss of generality we can assume the wave crest to be located at \(x=0\), so that \(x=\pm \lambda /2\) are trough lines and \(\eta (x)=\eta (-x)\) for all \(x \in {\mathbb R}\), with \(\eta \) increasing on \( [-\lambda /2,0]\) and decreasing on \([0,\lambda /2]\); see Fig. 1.
Structural properties of the system (9)–(15) can be used to reduce the number of unknowns without loss of information. First notice that the pressure can be eliminated by taking advantage of Bernoulli’s law: the expression
is constant throughout the fluid. A further reduction can be achieved by introducing the stream function \(\psi (x,y)\), defined up to an additive constant by
Note that the (11)–(12) ensure that \(\psi \) is harmonic throughout the fluid domain (in the moving frame)
while the kinematic boundary conditions (13)–(14) together with the non-existence of stagnation points, (16), ensure that \(\psi \) is constant and equal to its minimum value \(\psi _{min}\) on the free surface \(y=\eta (x)\), and constant equal to its maximum value \(\psi _{max}\) on the flat bed \(y=-d\). For \(y \in [-d,\eta (\lambda /2))]\) we have
and this last term turns out to vanish since (11) and the periodicity of u yield
while (14) gives
Similarly one can show that
The maximum principle for the harmonic function \((x,y) \mapsto \psi (x+\lambda ,y)-\psi (x,y)\) in the rectangular domain
now yields that \(\psi (x,y)\) is \(\lambda \)-periodic in the x-variable throughout the whole fluid domain \({{\mathcal {D}}}\). A straightforward argument using the maximum principle for the harmonic harmonic function \((x,y) \mapsto \psi (x,y)-\psi (-x,y)\) in \({{\mathcal {D}}}_0\) then shows that, due to the symmetry of the wave profile, we must have
and consequently we get
(Note that quite often the fact that u is even and v is an odd in the x-variable is assumed to be part of the symmetry of the flow pattern, while the above argumenation shows that these properties have to hold if the travelling wave profile is symmetric.) The relations (19) show that \(v(0,\cdot ) = 0\) and, since v is periodic in the x variable, we also have \(v(\pm \lambda /2,\cdot )=0\). Using these, the kinematic boundary conditions (13)–(14) together with the monotonicity of \(\eta \), we can now deduce that
by the maximum principle applied to the harmonic function v.
Let us now define the velocity potential \(\varphi (x,y)\) as the harmonic conjugate of \(\psi \), defined uniquely up to an additive constant by
From (21) we deduce that \(\varphi \) is constant and equal to its minimum value \(\varphi _{min}\) on the vertical segment \(x=\lambda /2\), and constant equal to its maximum value \(\varphi _{max}\) on \(x=-\lambda /2\). For simplicity, we shall from now on assume without loss of generality that \(\psi _{min}=\varphi _{min}=0\), so that
is the relative mass flux (that is, the mass flux relative to the uniform horizontal flow at speed c; see [2]). Since the function \((x,y) \mapsto \varphi (x+\lambda ,y)-\varphi (x,\lambda )\) is the harmonic conjugate of the identically zero harmonic map \((x,y) \mapsto \psi (x+\lambda ,y)-\psi (x,\lambda )\) throughout \({\mathcal D}\), we deduce that this map is constant. Thus
for any (x, y) in the fluid domain \({{\mathcal {D}}}\). The fact that at any depth y below the wave trough level \(\eta (\lambda /2)\) we have
means that \(-\varphi _{max}\) represents the strength of the uniform current underlying the irrotational flow, currents being defined as means of the horizontal velocity component (see [2]).
In light of (11)–(12), the function \(F=\varphi + \mathrm{i}\,\psi \), called the hodograph transform, is analytic in the interior of \({{\mathcal {D}}}\) and we have \(F'=(u-c)+\mathrm{i}(-v)\). In particular, F is a biholomorphic map between the interior of periodicity cell \(\{(x,y):\ x\in [-\lambda /2,\lambda /2] , -d \le y \le \eta (x)\}\) and the open rectangle \(\{(q,p):\ 0<q<\varphi _{max},\ 0<p < \psi _{\max }\}\), and it extends to a homeomorphism between the closures of these domains. The hodograph transform induces the following orientation-reversing conformal change of variables
which is a global diffeomorphism from the fluid domain \({{\mathcal {D}}} \) to the closure of the horizontal strip
Note that the conformal bijection \(q + \mathrm{i}p \mapsto x+ \mathrm{i}y\) from the horizontal strip to the fluid domain in the moving frame, obtained as the inverse of the change of variables (23),
has the “periodicity” properties (see Fig. 1)
where the first two relations above are mere reformulations of (22) and the periodicity of \(\psi \). The existence of a conformal mapping \(q + \mathrm{i}p \mapsto x+ \mathrm{i}y\) from a horizontal strip \({{\mathcal {R}}}_{\psi _{max}}\) to \({{\mathcal {D}}}\) that admits an extension as a homeomorphism between the closures of these domains and satisfies (25) is contingent upon the preservation of the scaling ratio \(\psi _{max}/\varphi _{max}\) (a geometric invariant). Indeed, the conformal diffeomorphism \((q+\mathrm{i}p) \mapsto \mathrm{e}^{2\pi \mathrm{i}(q+ip)/\varphi _{max}}\) from the periodicity rectangle \(\{(q,p):\ 0 \le q <\varphi _{max},\ 0 \le p \le \psi _{\max }\}\) to the annulus \(\{z \in {\mathbb C}: \ \mathrm{e}^{-2\pi \psi _{max}/\varphi _{max}} \le |z| \le 1\}\) shows that the claim is equivalent to a statement about the conformal equivalence of two such annuli, the necessary and sufficient condition being that the annuli have the same ratio of outer radius to inner radius (see [12]).
3 The streamline time-period
In this section we define the streamline time-period and establish some fundamental properties of this flow characteristic.
The streamlines (i.e. the level sets of \(\psi \)) are smooth curves, since for every \(p \in [0,\psi _{max}]\) and for every \(x \in {\mathbb R}\) there exists a unique \(f_{p}(x) \in [-d,\eta (x)]\) such that \(\psi (x,f_{p}(x))=p\) as a consequence of the implicit function theorem, with
From the symmetry (18) of \(\psi \) it follows that \(f_{p}\) is an even function of x. The periodicity of \(\psi \) in the x variable ensures that \(f_{p}\) is periodic with period \(\lambda \). Furthermore, (16) combined with (20) shows that \(f_{p}\) is increasing on \( [-\lambda /2,0]\) and decreasing on \( [0,\lambda /2]\).
Due to smoothness and the boundedness of the velocity field, the system describing the particle trajectories,
with initial data \(x(t_0)=x_0\) and \(y(t_0)=y_0\), has a unique global solution \((x(t;t_0,x_0,y_0), y(t;t_0,x_0,y_0))\) that depends smoothly on t and on the initial data \((t_0,x_0,y_0)\). It is easy to check that a solution viewed in the moving frame belongs to a level set of \(\psi \) , i.e. there is some \(p \in [0,\psi _{max}]\) such that
To indicate that a solution corresponds to the level set \(\psi =p\) we shall write \((x_{p}(t;t_0,x_0,y_0),\,y_{p}(t;t_0,x_0,y_0))\). Note that
Obviously, each \((x_0,y_0)\) belongs at time \(t_0\) to some streamline \(\psi =p\), so that \(y_0=f_{p}(x_0-ct_0)\). Thus, since \(f_{p}\) is smooth, the dependence of the solution \((x_{p}(t;t_0,x_0,y_0),\,y_{p}(t;t_0,x_0,y_0))\) to (27) on \(y_0\) is actually an implicit smooth dependence on \(x_0\) and \(t_0\). From now on we write \((x_{p}(t;t_0,x_0),\,y_{p}(t;t_0,x_0))\).
Our first result ensures that the particle trajectories viewed in the moving frame are periodic in time, with a period depending only on the streamline.
Theorem 1
Let \(p \in [0,\psi _{max}]\). Then for each \((t_0,x_0) \in {\mathbb R}^2\) the equation
has a unique solution \(T=T(p)>0\) that is independent of \((t_0,x_0)\), and the following relation holds:
Furthermore, the map \(p \mapsto \log { T}(p)\) is a convex function and the map \(p \mapsto { T}(p)\) is nonincreasing.
Proof
We first show that equation (29) has a unique solution \(T=T(p,t_0,x_0)\). To this end, let
We have
where
because there are no stagnation points in the fluid domain \({{\mathcal {D}}}\). Since \(h(0)=\lambda >0\), this implies
thus ensuring the existence of \(T=T(p,t_0,x_0)>0\) with \(h(T)=0\). The uniqueness of T follows from the strict monotonicity of h.
To see that T is independent of \((t_0,x_0)\), consider the map
Then
and, by the definition of T, we have \( \sigma (t_0)=x_0-ct_0\) and \(\sigma (T+t_0)=x_0-ct_0-\lambda \). Using these facts, we obtain the following formula for T:
where the last step follows by periodicity. It is now clear that T depends only on p. Furthermore, (30) follows easily from the periodicity of \(f_{p}\). Indeed,
Now our strategy to show the convexity of \(p \mapsto \log { T}(p)\) is to use the hodograph transform to turn the integral expression for T into an integral of the modulus of an analytic function over a horizontal segment, which, using the periodicity, can be subsequently transformed into an integral over a circle, to which we may apply Hardy’s convexity theorem.
Recall that the hodograph transform \(F=\varphi + \mathrm{i}\,\psi \) is analytic and we have \(F'=(u-c)+\mathrm{i}(-v)=\varphi _x-\mathrm{i}\varphi _y\). Inspired by this transform we perform a change of variables in the integral expression for T by setting
and
Hence
Notice that
which implies
The function \(G:=(F^{-1})'\) is analytic in the rectangle \([0,\varphi _{max}] \times [0,\psi _{max}]\) and the periodicity of \(F'\) gives
By an appropriate transformation (see Fig. 2) we now pass from the periodicity rectangle \([0,\varphi _{max}) \times [0,\psi _{max}]\) to an annulus, where we may apply Hardy’s theorem. Denote
and consider the principal branch of the logarithm,
Then (see Fig. 2)
maps \({\mathbb A}\) bijectively onto \([0,\varphi _{max}) \times [0,\psi _{max}]\), and the map \(G \circ \beta \) is analytic in the interior of \({\mathbb A}\) with the real-line segment \((-\mathrm{e}^{2\pi \frac{\psi _{max}}{\varphi _{max}}},-1)\) excised. By property (35) the jump of \(\text {Log}\) gets smoothed out by the composition with G, and hence \(G \circ \beta \) extends to a continuous map on \({\mathbb A}\). By Morera’s theorem it now follows that \(G \circ \beta \) is analytic on \({\mathbb A}\).
Let us now fix \(p \in [0,\psi _{max}]\) and set \(r:=\mathrm{e}^{\frac{2\pi }{\varphi _{max}}\,p}\). Then \(1 \le r \le \mathrm{e}^{2\pi \frac{\psi _{max}}{\varphi _{max}}}\) and we have
where the last step above follows by the change of variables \(q=\tfrac{\varphi _{max}}{2} \big (1- \tfrac{\theta }{\pi } \,\big )\). Hardy’s convexity theorem (see [13]) now ensures that the map \(r \mapsto \log M_2(G \circ \beta , r)\) is a convex function of \(\log r\) for \(r \in [1,\mathrm{e}^{2\pi \frac{\psi _{max}}{\varphi _{max}}}]\), i.e. for any \(r_1, r_2\in [1,\mathrm{e}^{2\pi \frac{\psi _{max}}{\varphi _{max}}}]\) and for every \(\alpha \in [0,1]\) we have
where \(\log r=\alpha \log r_1+(1-\alpha ) \log r_2\). Since \( \log r=\tfrac{2\pi }{\varphi _{max}}\,p\), we deduce from (36) that \(p \mapsto \log {T}(p)\) is a convex function.
The monotonicity of T now follows combining the above with the fact that v vanishes on the flat bed. Indeed, the convexity of \(p \mapsto \log {T}(p)\) implies the convexity of \(p \mapsto {T}(p)\) which, at its turn, implies that \(T'(p)\) is nondecreasing on \([0,\psi _{max}]\). Therefore, in order to show that T is nonincreasing it suffices to show that \({T}'(\psi _{max})=0\). Due to (34), we have
We claim that the evaluation of the integrand above at \(p=\psi _{max}\) vanishes. Denoting the components of \(F^{-1}\) by x and y, using the chain rule we obtain
We now evaluate the above expression at \(p=\psi _{max}\), so that all the terms containing u and v above are to be evaluated at
Since \(v=0\) on \(y=-d\) and \(u_y(\cdot ,-d)=v_x(\cdot ,-d)\) we infer that the last three of the four terms on the right side of (38) vanish. Furthermore, computing the inverse of the Jacobian matrix of F,
so that the first term on the right side of (38) vanishes as well, and we can now conclude that \({T}'(\psi _{max})=0\). Thus T is nonincreasing and, with this, the proof is complete.\(\qquad \Box \) \(\square \)
Remark 1
If \((x_0,y_0)\) is the location of a particle at time \(t=0\), then there exists a unique \(p \in [0,\psi _{max}]\) with \(y_0 = f_{p}(x_0)\), and the location of the particle at time \(t >0\) is given by the solution \((x_{p}(t;x_0),\,y_{p}(t;x_0))\) to (27) with initial data \((x_{p}(0),\,y_{p}(0))=(x_0,y_0)\). Theorem 1 ensures
so that the elevation is T(p)-periodic while the horizontal position in shifted by \(c\,T(p) - \lambda \) in the direction of wave propagation. Moreover, by Theorem 1 this shift depends monotonically on p, or, equivalently, on the minimal (or maximal) elevation of the particle, since \(\psi _y=u-c<0\) throughout the fluid domain. We infer that the particle trajectory consists of a repeating pattern, that is simply shifted horizontally after each streamline time-period T(p), and this pattern depends solely on the streamline. For the precise pattern we refer to the theoretical study [4], the numerical simulations in [3] and the experiments reported in [16]. The fact that the pattern is repeated after a suitable shift was established in [4] by considering particles that are located initially beneath the wave crest and arguing that the relation \(\psi _y=u-c<0\) ensures that any particle will occupy this position at some time. Theorem 1 via relations (39) shows the uniformity in time of the particle trajectories, a result that is to be expected by glancing at the patterns depicted in [3, 4, 16], but which was not proven hitherto. \(\square \)
4 The total kinetic energy of a particle
The main aim of this section is to investigate the total kinetic energy of an arbitrary particle over a streamline time-period. This quantity turns out to depend only on the corresponding streamline and we pursue a detailed investigation of this dependence.
Theorem 2
The total kinetic energy of a particle located initially at \((x_0,y_0)\) with \(y_0=f_{p}(x_0)\), given by
depends only on the corresponding streamline \(\psi =p\) (i.e. is independent of \(x_0\)) and the map \(p \mapsto \log {{\mathcal {K}}}(p)\) is a convex function, while the map \(p \mapsto {{\mathcal {K}}}(p)\) is nonincreasing. Furthermore, the total kinetic energy of the particle relative to the moving frame, computed over a streamline time-period, given by
is constant equal to \( \varphi _{max}/2\).
Remark 2
Theorem 2 ensures that the total kinetic energy of a particle increases with the elevation of the streamline over the flat bed (measured either below the crest or, equivalently, below the trough), i.e. the total kinetic energy of a particle above the flat bed is larger than that of a particle moving beneath it.
Proof
The proof is quite similar to the one of Theorem 1, so that we shall not provide all the details.
As already mentioned, the hodograph transform \(F=\varphi + \mathrm{i}\psi \) is analytic and we have \(F'=(u-c)+\mathrm{i}(-v)\). Using this notation we obtain
due to (28). We first make the change of variables
to get
where the last step follows by periodicity. This shows that \({{\mathcal {K}}}\) is independent of the particle, depending only on the corresponding streamline \(\psi =p\).
As before, inspired by the hodograph transform, we perform a further change of variables
Hence, by (31)–(32)–(33), we obtain
The function \(H=1+ c(F^{-1})'\) is analytic in the interior of the rectangle \([0,\varphi _{max}] \times [0,\psi _{max}]\) and the periodicity of \(F'\) gives
The last expression for \({{\mathcal {K}}}(p)\) is given by an integral on a horizontal segment of the modulus of an analytic function. The exact same argument as the one used in the proof of Theorem 1 now shows that \(p \mapsto \log {{\mathcal {K}}}(p)\) is a convex function.
The second assertion follows by taking advantage of the first one. Indeed, the convexity of \(p \mapsto \log {{\mathcal {K}}}(p)\) implies the convexity of \(p \mapsto {{\mathcal {K}}}(p)\) which, at its turn, implies that \({{\mathcal {K}}}'(p)\) is nondecreasing on \([0,\psi _{max}]\). Therefore, in order to show that \({{\mathcal {K}}}\) is nonincreasing it suffices to show that \({\mathcal K}'(\psi _{max})=0\). We have
We claim that the evaluation of the integrand above at \(p=\psi _{max}\) vanishes. From (38) we have
and since this formula is actually valid also for \(c=0\), we can now conclude that \({{\mathcal {K}}}'(\psi _{max})=0\). Thus \({{\mathcal {K}}}\) is nonincreasing.
The last assertion in our statement follows immediately by letting \(c=0\) in the above argument. More precisely, for \(c=0\), relation (40) becomes \({\mathbb K}(p,x_0)=\varphi _{max}/2\). \(\qquad \Box \)
Remark 3
Using exactly the same method as above, one can prove that the conclusion of Theorem 2 remains true for the \(L^s\)-norm of the kinetic energy, i.e. for expressions of the form
for the physical frame, respectively
for the moving frame; here \(s>0\) is a rational number. More precisely, one can show that \({{\mathcal {K}}}_s(p,x_0)\) and \({\mathbb K}_s(p,x_0)\) are independent of \(x_0\), the functions \(p \mapsto \log {{\mathcal {K}}}_s(p)\) and \(p \mapsto \log {\mathbb K}_s(p)\) are convex and the maps \(p \mapsto {{\mathcal {K}}}_s(p)\) and \(p \mapsto {\mathbb K}_s(p)\) are nonincreasing. Indeed, if \(s=m/n\) with \(m,n\in \mathbb {N},\,n\ne 0\), the analogue of (40) becomes
the analytic function to which we apply Hardy’s convexity theorem being in this case \(\frac{(F'+c)^{m}}{(F')^n}\,(F^{-1}(q,p))\). Analogously, we have
and we apply Hardy’s convexity theorem to the analytic function \(F'(F^{-1}(q,p))\). Notice that \({\mathbb K}_1(p)=\varphi _{max}/2\), as already shown in Theorem 2.\(\Box \)
Remark 4
The total kinetic energy (in the moving frame) on a streamline \(\gamma _p\), with
is given by
where s represents the arclength parameter. We can show that the map \(p \mapsto \log {{\mathcal {E}}}(p)\) is convex, while \(p \mapsto {{\mathcal {E}}}(p)\) is increasing.
The method of the proof is the same as for Theorems 1–2, so we only sketch some of the main steps. Using the derivative of the hodograph transform \(F'=(u-c) + \mathrm{i} (-v)\) together with relation (26), we infer
since \(\varphi _x=u-c<0\) by (16), (17) and (21). Performing the change of variables \(q=\tilde{q}_p(x)=\varphi (x,f_p(x))\) and using the relations (31)–(33), one finds
where the last equality follows from (42).\(\Box \)
5 Physical motivation
In order to explain the broader relevance of our considerations, let us note that the energy of a water wave comprises the kinetic energy of the orbital motion of the particles beneath the surface and the potential energy associated with the vertical rise and fall of the water surface from its still-water, undisturbed state. Ocean wave energy is becoming an important renewable source due to its vast but yet somewhat untapped potential. In particular, with wave power density averages up to 50 kW per meter of wave front near the coast, and up to 100 KW/m offshore (see [7]), ocean wave energy compares favourably with typically less energy-dense sources such as wind and solar, with the added benefit that the lower profile of the energy converters reduces their visual impact. The fact that ocean wave harvesting devices convert the kinetic energy into electricity (see [17]) or make it available directly for other purposes, like taking advantage of the associated pressure to run desalinization processes (see [8]), provides a practical motivation for a better understanding of the structure of the kinetic energy of water waves.
Change history
References
Amick, C.J., Toland, J.F.: On periodic water-waves and their convergence to solitary waves in the long-wave limit. Philos. Trans. R. Soc. A 303, 633–669 (1981)
Bühler, O.: Waves and Mean Flows. Cambridge University Press, Cambridge (2014)
Clamond, D.: Note on the velocity and related fields of steady irrotational two-dimensional surface gravity waves. Philos. Trans. R. Soc. Lond. A 370, 1572–1586 (2012)
Constantin, A., Strauss, W.: Pressure beneath a Stokes wave. Comm. Pure Appl. Math. 53, 533–557 (2010)
Constantin, O., Persson, A.-M.: A complex-analytic approach to kinetic energy properties of irrotational flows. Proc. Am. Math. Soc. 150, 2647–2653 (2022)
Craig, W., Sulem, C.: Numerical simulation of gravity waves. J. Comput. Phys. 108, 73–83 (1993)
EPRI: Mapping and assessment of the United States ocean wave energy resource (2011)
Leijon, J., Boström, C.: Freshwater production from the motion of ocean waves–a review. Desalination 435, 161–171 (2018)
Longuet-Higgins, M.S.: On the decrease of velocity with depth in an irrotational water wave. Math. Proc. Camb. Philos. Soc. 49, 552–560 (1953)
Nachbin, A., Ribeiro-Junior, R.: A boundary integral formulation for particle trajectories in Stokes waves. Discrete Contin. Dyn. Syst. 34, 3135–3153 (2014)
Nachbin, A., Ribeiro-Junior, R.: Capturing the flow beneath water waves. Philos. Trans. R. Soc. A 376, 17 (2018). (Art. 20170098)
Rudin, W.: Real and Complex Analysis. McGraw-Hill Book Co., New York (1987)
Sarason, D.: The \(H^p\) spaces of an annulus. Mem. Am. Math. Soc. 56, 78 (1965)
Strauss, W.: Steady water waves. Bull. Am. Math. Soc. 47, 671–694 (2010)
Toland, J.F.: Stokes waves. Topol. Methods Nonlinear Anal. 7, 1–48 (1996)
Umeyama, M.: Eulerian-Lagrangian analysis for particle velocities and trajectories in a pure wave motion using particle image velocimetry. Philos. Trans. R. Soc. Lond. A 370, 1687–1702 (2012)
Xie, J., Zuo, L.: Dynamics and control of ocean wave energy converters. Int. J. Dyn. Control 1, 262–276 (2013)
Acknowledgements
The author is grateful for helpful comments from the referee, which brought considerable improvement.
Funding
Open access funding provided by University of Vienna.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Constantin, O. A complex-analytic approach to kinetic energy properties of irrotational traveling water waves. Math. Z. 301, 4201–4215 (2022). https://doi.org/10.1007/s00209-022-03042-3
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00209-022-03042-3