Exact dispersion relation for linear surface waves on arbitrary vertical shear

This paper derives an exact, implicit dispersion relation for linear surface waves on an arbitrary vertical shear profile by employing a Green's function framework for the Rayleigh equation, providing a general solution that reduces to known analytical and asymptotic limits.

The Gist

Imagine you are standing on a pier, watching waves roll toward the shore. Usually, we think of waves as simple rhythmic motions on the surface. But in the real ocean, there is a hidden complexity: beneath that surface, the water isn't just sitting still. It’s moving in currents, and those currents often move at different speeds at different depths—like a multi-layered cake of moving water.

This paper, written by Kjell S. Heinrich and Simen Å. Ellingsen, provides a new mathematical "master key" to understand how these deep, invisible currents change the way surface waves behave.

Here is the breakdown of their discovery using everyday analogies.

1. The Problem: The "Invisible Hand" of the Current

Think of a wave as a runner sprinting across a field. If the field is flat and still, the runner’s speed is easy to predict. But what if the field is actually a giant, moving treadmill? And what if the treadmill isn't just one speed, but has different layers—the top layer moves fast, the middle layer moves medium, and the bottom layer moves slow?

The runner (the wave) is constantly being pushed or pulled by these different layers of movement. Because the wave isn't just a flat line—it has depth—it "feels" all those different speeds at once. For decades, scientists have struggled to write a single, perfect equation that accounts for every possible way a current can change speed from the surface to the bottom. They usually had to make "guesses" or simplifications to make the math work.

2. The Solution: The "Path-Ordered" Recipe

The authors didn't want to make guesses. They wanted an exact answer. To do this, they borrowed a trick from Quantum Physics.

In quantum physics, particles don't just move in straight lines; they interact with everything in their path. The authors treated the "curvature" of the current (the way the speed changes) like a series of obstacles or "scattering events."

The Analogy: The Musical Journey
Imagine you are a musician playing a melody (the wave) while walking through a series of rooms. In each room, there is a different wind blowing (the current shear).

• In the first room, the wind pushes you left.
• In the second, it pushes you right.
• In the third, it swirls.

To know exactly how your melody sounds when you exit the building, you can't just average the wind speeds. You have to account for the order in which you encountered the rooms. If you hit the "left wind" room first and then the "right wind" room, the result is different than if you had gone in reverse.

The authors created a mathematical tool called a "path-ordered exponential." This is essentially a highly sophisticated "travel log" that tracks exactly how the wave is modified by every single layer of the current, in the exact order it encounters them, from the bottom of the ocean up to the surface.

3. Why This Matters: The "Effective Depth" Trick

One of the coolest things they discovered is that a complex, swirling current can be thought of as a "phantom depth."

Imagine you are swimming in a pool that is 10 feet deep. Suddenly, a strong current starts moving at different speeds. Even though the floor of the pool hasn't moved, the wave starts acting as if the pool is suddenly 12 feet deep or only 8 feet deep. The current "tricks" the wave into thinking the ocean floor is at a different level. This "effective depth" gives oceanographers a much more intuitive way to visualize what is happening.

4. The Big Picture: Why should we care?

Why spend all this time on math? Because waves are the "messengers" of the ocean.

• Predicting Rogue Waves: If we want to know if a massive, ship-breaking "rogue wave" is forming, we need to know exactly how the currents are shaping the waves.
• Protecting Ships and Coasts: Understanding how waves hit a ship or a sea wall depends on knowing their true speed and power.
• Pollution Tracking: If we want to know where plastic waste or an oil spill will go, we have to understand the relationship between the surface waves and the currents below.

In short: The authors have moved us from using "blurry snapshots" of wave behavior to a "high-definition video." They have provided a mathematical framework that is exact, handles any kind of current, and works perfectly, whether the water is shallow or deep.

Technical Summary

Technical Summary: Exact Dispersion Relation for Linear Surface Waves on Arbitrary Vertical Shear

Authors: Kjell S. Heinrich and Simen Å. Ellingsen (NTNU)

1. Problem Statement

The paper addresses a fundamental problem in fluid mechanics and oceanography: determining the dispersion relation for linear surface gravity waves propagating on a horizontal mean current with an arbitrary vertical shear profile $U(z)$.

While classical solutions exist for specific cases (e.g., constant vorticity or deep-water limits), most real-world oceanographic scenarios involve complex, non-linear vertical velocity profiles. Previous approaches have relied either on:

• Approximations: Such as weak-shear or weak-curvature expansions, which fail when shear is strong or the profile is highly irregular.
• Numerical Eigenvalue Solvers: Such as the Direct Integration Method (DIM), which are computationally intensive and do not provide a closed-form analytical understanding of the physics.

The goal of this work is to derive a single, exact, implicit equation that relates the wave frequency ($\omega$), wavenumber ($k$), and the velocity profile $U(z)$ without requiring the calculation of an unknown vertical eigenfunction.

2. Methodology

The authors employ a sophisticated mathematical framework to bypass the traditional eigenvalue problem:

• Green’s Function Framework: Instead of solving for the wave velocity field directly, the problem is recast using a Green’s function $\hat{G}(z)$ for the Rayleigh equation (the inviscid Orr-Sommerfeld equation). They model the surface disturbance as coming from a submerged point source.
• The "Interaction Picture" Approach: Inspired by quantum mechanics, the authors split the differential operator into a "solvable" part (representing the shear-free case) and a "difficult" part (representing the curvature/shear effects).
• Path-Ordered Exponentials: To handle the non-commuting nature of the matrices arising from the vertical evolution of the Green's function, the authors derive the solution using a path-ordered exponential. This allows them to account for how curvature at different depths "interacts" with itself as the wave propagates vertically.
• Wronskian Matching: They construct two fundamental solutions—one satisfying the bottom boundary condition and one satisfying the free-surface condition—and match them at the surface using a Wronskian. The zeros of this Wronskian define the dispersion relation.

3. Key Contributions

• Exact Implicit Dispersion Relation: The primary contribution is the derivation of Equation (25):
  $$\sigma^2 = \left( \frac{g + \Upsilon k^2}{k} - \sigma U'_\gamma(0) \right) F(k, \sigma)$$
  where $F(k, \sigma)$ is a function of the ratio $r = Y_1(0)/Y_2(0)$, derived from the path-ordered integral of the velocity profile.
• Dyson Series Representation: The authors show that the solution can be expressed as a Dyson series (similar to multiple scattering in particle physics). This provides a systematic way to view the dispersion relation as a series of "scatterings" caused by the curvature of the current at various depths.
• Unified Framework: The solution is "all-encompassing." It is shown to reduce to:
  • The classical no-shear solution.
  • The constant vorticity (linear shear) solution.
  • The Shrira (1993) deep-water limit.
  • Known weak-shear and weak-curvature approximations.

4. Results and Validation

The authors validate their exact analytical formula through two primary methods:

• Analytical Limits: They successfully recover the exact solutions for constant vorticity and the specific class of exponential/trigonometric profiles studied by Peregrine (1976).
• Numerical Benchmarking: They implemented their formula numerically and compared it against the Direct Integration Method (DIM) using complex, sixth-order polynomial velocity profiles (mimicking wind-induced currents). The results showed identical agreement between their implicit formula and the DIM, with errors limited only by machine precision.

5. Significance

This paper provides a significant theoretical advancement for several reasons:

1. Mathematical Elegance: It transforms a complex eigenvalue problem into a single-equation implicit problem, eliminating the need to solve for the vertical structure of the wave.
2. Computational Efficiency: The path-ordered product approach (discretized as a product of propagators) offers a highly efficient and accurate way to evaluate dispersion in complex flows.
3. Physical Insight: The "multiple scattering" analogy and the "effective depth" interpretation ($F = \tanh k(h + \tilde{h})$) provide new intuitive ways to understand how current curvature modifies wave propagation.
4. Broad Application: The method is applicable to a wide range of geophysical problems, including wave forecasting, rogue wave formation, and the study of wave-current interactions in coastal and deep-ocean engineering.