A variational formulation for steady surface water waves on a Beltrami flow

J. Horn; Mark D. Groves

REVIEW 2 major objections 4 minor 14 references

A variational formulation for steady surface water waves on a Beltrami flow

T0 review · 2 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Steady Beltrami water waves are equivalent to two scalar surface equations that arise from a variational principle, with an operator H(η) that generalises the Dirichlet–Neumann operator.

desk verdict Correct variational derivation, but Theorem 3.1 as printed has a spurious α-term that breaks the stated equivalence; easy fix, still worth publishing. read the letter →

arxiv 1908.09386 v1 pith:OHAF2LH2 submitted 2019-08-25 math.AP

Mark D. Groves , J. Horn This is my paper · ORCID

classification math.AP MSC 76B1535Q35

keywords BeltramiflowssteadywaterwavesvariationalprincipleDirichlet-NeumannoperatorHodge-Weyldecompositionnonlocalfree-boundaryproblem

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper claims that steady water waves riding a Beltrami flow—a three-dimensional flow whose velocity and vorticity are parallel—can be described by just two scalar functions of the horizontal coordinates: the surface elevation η and a surface potential Φ. The full free-boundary problem with curl u = αu reduces to the two Euler–Lagrange equations of a single variational functional L(η, Φ), written through a nonlocal operator H(η) that maps Φ to the normal velocity at the surface. In the irrotational limit α = 0, H(η) becomes the Dirichlet–Neumann operator and the system reduces to the classical Zakharov–Craig–Sulem formulation. A sympathetic reader would care because the result recasts a difficult three-dimensional free-boundary problem as a two-dimensional variational structure, opening the way to existence theories and numerical schemes built on the surface variables alone.

What carries the argument

The nonlocal operator H(η), defined by H(η)Φ = curl A·N|y=η where A is the unique solution of (1.15)–(1.19), is the object that carries the argument. It converts the surface potential Φ into the normal component of velocity, plays the role of the Dirichlet–Neumann operator for Beltrami flows, and gives the compact form L(η, Φ) = ∫(½ΦH(η)Φ − ∇Φ·A⋆⊥‖ + Γ(η) + ½gη² + σ((1+|∇η|²)^{1/2}−1)). Its formal self-adjointness (Lemma 3.1) makes the variational structure possible, and its analytic dependence on η, proved by flattening the domain and applying the analytic implicit-function theorem, supplies the functional-analytic foundation. The companion operator K(η)Φ = ∇Φ − α∇⊥Δ⁻¹(H(η)Φ) expresses the tangential velocity field that appears in the dynamic surface equation.

What would settle it

Numerically solve the flattened boundary-value problem (4.7)–(4.11) for a fixed surface elevation η away from zero and a value of |α| exceeding the paper's smallness threshold; if two distinct Hs solutions A emerge for the same Φ, then H(η) is multivalued, the variational principle δL = 0 is not single-valued, and the claimed reduction collapses at that parameter point.

Watch

Extended reading notes

Core claim

The central discovery is that the hydrodynamic problem (1.7)–(1.11) is equivalent to two equations for (η, Φ): H(η)Φ + u⋆·N = 0 together with a dynamic surface equation, and that these are exactly the Euler–Lagrange equations of δL(η, Φ) = 0 for the functional (1.14). Here H(η) is defined by solving a boundary-value problem for a vector potential A with curl curl A = α curl A, and setting H(η)Φ = curl A·N at the surface; the Hodge–Weyl decomposition of the tangential velocity selects the gradient component Φ as the second unknown. The variational functional combines a Woltjer-type energy for Beltrami fields with the surface-energy terms of the water-wave problem. When α = 0 the vector-potential construction returns a harmonic scalar potential, H(η) reduces to the Dirichlet–Neumann operator G(η), and the two equations become the classical steady water-wave equations in Zakharov–Craig–Sulem form.

Load-bearing premise

The whole reduction rests on the boundary-value problem for the vector potential A having a unique solution, and this is proved only for small values of the vorticity strength |α| and for surface elevations η in a neighbourhood of zero; outside that regime H(η) may not be defined and the two-surface-variable formulation may fail.

Editorial extensions

If this is right

The three-dimensional steady Beltrami water-wave problem is fully captured by two scalar surface equations, so any solution of the variational principle automatically satisfies the kinematic and dynamic boundary conditions.
The generalised Dirichlet–Neumann operator H(η) is formally self-adjoint and depends analytically on η, so perturbation and bifurcation methods for the surface variables can be applied in the small-α, small-η regime.
In the irrotational limit the formulation reproduces the classical Zakharov–Craig–Sulem equations, making the Beltrami case a genuine extension rather than an unrelated model.
The variational principle merges the classical fixed-domain variational principle for Beltrami fields with surface-energy terms, making energy-based existence methods available for steady waves riding Beltrami flows.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

If the well-posedness of the vector-potential boundary-value problem can be established beyond the small-α, small-η regime, the two-surface-variable formulation would provide a variational or Hamiltonian setting for large-amplitude steady Beltrami waves, paralleling the role of the Zakharov–Craig–Sulem equations in irrotational theory.
The operator H(η) is a natural target for numerical simulation: replacing the full three-dimensional elliptic solve with evaluations of H(η) would lower the computational dimension, provided the defining boundary-value problem can be solved efficiently for general η.
The same Hodge–Weyl selection of the gradient component Φ might adapt to other helical or force-free flows, such as magnetohydrostatic equilibria, where a reduction to surface variables could be sought.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Core idea is right: the paper derives a two-scalar variational formulation for steady Beltrami water waves, reducing the free-boundary problem to two equations for (η,Φ). The Section 2 derivation is careful, and the α=0 limit reproduces the Zakharov-Craig-Sulem system. That part stands up.

But the printed Theorem 3.1 contains a spurious α-term in the dynamic equation. Substituting (2.5) into (2.6) and expressing everything in terms of K and H cancels the term -α H(H+K·∇η)/(1+|∇η|^2). The correct equation has no α-term. As printed, the claimed equivalence between the two-scalar system and (1.7)-(1.11) fails. The stress-test is right. The good news is that the variational derivation in Section 2 is unaffected, and the fix is deleting one term from a display. This is an algebraic slip, not a methodological failure.

The functional-analytic part is solid: Lemma 4.1 and Theorem 4.1 establish existence, uniqueness, and analyticity of H(η) for small |α| and small η, an honest limitation the paper states. Proposition 4.7 has a minor typo in its statement, but the Green's function argument is believable. Citations are clean; the distinction from Lokharu-Wahlen is accurately drawn.

The paper deserves peer review. The variational principle is a useful structural contribution, and the error is trivial to repair. I would ask the authors to correct Theorem 3.1, re-check the display, and then accept. The smallness restriction should stay front and center, but it is not a flaw.

Referee Report

2 major / 4 minor

Summary. The paper proposes a two-scalar reduction for steady three-dimensional water waves on a Beltrami flow. The velocity perturbation is written as v = curl A, and a nonlocal operator H(η), defined through a boundary-value problem for A, maps the surface potential Φ to the normal velocity at the free surface. The authors show formally that the free-boundary problem is equivalent to two equations for (η, Φ): a kinematic-type equation H(η)Φ + u⋆·N = 0 and a scalar dynamic equation, and that these are the Euler-Lagrange equations of a variational functional L. For α = 0 the operator is claimed to reduce to the Dirichlet-Neumann operator and the system to the Zakharov-Craig-Sulem formulation. Section 4 supplies a Hodge-Weyl decomposition, well-posedness of the defining boundary-value problem for small |α|, and analyticity of H(η) for small η.

Significance. If the central equivalence were correct, this would be a valuable extension of the Zakharov-Craig-Sulem formulation to rotational Beltrami flows, with a clean variational structure for a genuinely three-dimensional free-boundary problem. The paper contains explicit formal computations, an independent α = 0 benchmark, and a rigorous local well-posedness and analyticity analysis, which are genuine strengths. However, the displayed two-scalar system in Theorem 3.1 contains a spurious α-term, and the quoted classical irrotational limit has a sign inconsistency. These issues are local and readily fixable, but until corrected the main equivalence statement is not reliable as printed.

major comments (2)

[Theorem 3.1 and §2, Eq. (2.6)] The dynamic equation displayed in Theorem 3.1 is not equivalent to the Euler-Lagrange system (2.5)–(2.6). Let S = ∇·A⊥_‖ + ∇·A⋆⊥_‖. The first equation of Theorem 3.1 is exactly S = 0. Under S = 0 the first two terms of (2.6) vanish, so (2.6) reduces to 1/2|v|² + v_h·u⋆_h + gη − σ(...) = 0. Using v·N = H, v‖ = K, v_2 = (H + K·∇η)/(1 + |∇η|²), and u⋆_h·∇η = −u⋆·N = H, one obtains 1/2|K|² − (H+K·∇η)²/(2(1+|∇η|²)) + K·u⋆_h − H(H+K·∇η)/(1+|∇η|²) + gη − σ(...) = 0. The factor α multiplying the last displayed term in Theorem 3.1 is spurious; no α survives after imposing (2.5). Since the α = 0 limit cannot detect this error, the two-scalar system as printed is not equivalent to (1.7)–(1.11) for α ≠ 0.
[Introduction, Zakharov-Craig-Sulem limit] The classical irrotational system quoted in the introduction is inconsistent with the α = 0 limit of the new formulation. From (2.5) with α = 0 and u⋆ = c one obtains HΦ + u⋆·N = G(η)Φ − c·∇η = 0, i.e., G(η)Φ = c·∇η, whereas the introduction states G(η)ξ + c·∇η = 0. This sign difference propagates into the dynamic equation and means that the claimed reduction to the Zakharov-Craig-Sulem equations, as written there, is not a correct benchmark. Either the sign convention in L0 or the displayed classical equations should be corrected so that the α = 0 check is internally consistent with (1.10)–(1.11).

minor comments (4)

[§4(c), Proposition 4.7] The proof of Proposition 4.7 relies on 'straightforward calculations' for the Green's-matrix estimates, but those estimates are a substantial part of the functional-analytic argument; a short derivation of the key estimates (A 1)–(A 5) would improve verifiability.
[Abstract and §4] The abstract and introduction state the two-scalar reduction without qualification, but Theorem 4.1 proves existence and uniqueness of the defining boundary-value problem only in a neighbourhood of η = 0 and for |α| < π/(2h). The claims should be phrased as local in amplitude and vorticity.
[§4(a)] In the uniqueness statement for the Hodge-Weyl decomposition, 'unique functions in ˙H¹(R)' should presumably read 'unique functions in ˙H¹(R²)'; please correct this typo.
[Theorem 3.1] After correcting the spurious α-term, it would be helpful to display the simplified dynamic equation, since the printed formula is likely to be copied by readers implementing the two-scalar system.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the variational reduction is self-contained and benchmarked against the classical irrotational limit.

full rationale

The paper's derivation chain builds the Beltrami water-wave problem from the stated boundary-value problem for the vector potential A, defines the nonlocal operator H(eta) by that same boundary-value problem, and then verifies that the Euler-Lagrange equations of L(eta, Phi) are algebraically equivalent to the original kinematic and dynamic boundary conditions. The identification H(eta)Phi = div A^perp_| = curl A . N is a definition followed by a calculation, not a fitted input or a prediction forced by construction. The irrotational limit alpha = 0 is checked against the independent, known Zakharov-Craig-Sulem formulation and the Dirichlet-Neumann operator, providing an external benchmark. The existence and uniqueness theory for A is proved in the paper using standard tools (Lax-Milgram, implicit function theorem, Green's function estimates), with no load-bearing appeal to an unverified self-citation. No parameter is fitted, no uniqueness theorem is imported from the authors' prior work, and no known result is merely renamed as a new one. Accordingly, the central claim is self-contained rather than circular; any possible algebraic issue in the displayed dynamic equation would be a correctness concern, not a circularity concern.

Assumptions & free parameters 0 free parameters · 3 assumptions · 0 invented entities

No free parameters are fitted: g, sigma, c, and alpha are physical inputs, and H(eta) is defined, not fitted. The axioms are the physical model, the geometric and decay assumptions on the free surface, and standard analytic tools. No new physical entities are introduced.

assumptions (3)

domain assumption The velocity field is a strong Beltrami flow, curl u = alpha u with fixed constant alpha, and the fluid is inviscid, incompressible, with unit density.
The entire problem (1.1)-(1.5) is posed under this model, and the variational formulation is developed within this class of flows.
domain assumption The free surface is a graph y = eta(x,z) with eta > -h, and the perturbation v is evanescent as |(x,z)| tends to infinity.
The Hodge-Weyl decomposition and the nonlocal operator H(eta) are defined on R^2 using this decay and geometric setting, as described in Sections 1 and 4(a).
standard math Standard analytic tools: Lax-Milgram lemma, analytic implicit function theorem, Sobolev embedding theorems, and the Hodge-Weyl decomposition for L^2 vector fields.
These unproved background results are used in Lemmata 4.1 and 4.2 and Theorem 4.1 to establish existence, uniqueness, and analytic dependence.

how reviews work

0 comments

Cite this review

Pith. "Pith review of A variational formulation for steady surface water waves on a Beltrami flow." pith.science (2026). https://pith.science/paper/OHAF2LH2

@misc{pith2026190809386,
  author       = {Pith},
  title        = {Pith review of: A variational formulation for steady surface water waves on a Beltrami flow},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/OHAF2LH2}},
  note         = {Machine review of arXiv:1908.09386}
}

abstract

This paper considers steady surface waves `riding' a Beltrami flow (a three-dimensional flow with parallel velocity and vorticity fields). It is demonstrated that the hydrodynamic problem can be formulated as two equations for two scalar functions of the horizontal spatial coordinates, namely the elevation $\eta$ of the free surface and the potential $\Phi$ defining the gradient part (in the sense of the Hodge-Weyl decomposition) of the horizontal component of the tangential fluid velocity there. These equations are written in terms of a nonlocal operator $H(\eta)$ mapping $\Phi$ to the normal fluid velocity at the free surface, and are shown to arise from a variational principle. In the irrotational limit the equations reduce to the Zakharov-Craig-Sulem formulation of the classical three-dimensional steady water-wave problem, while $H(\eta)$ reduces to the familiar Dirichlet-Neumann operator.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

14 extracted references · 14 canonical work pages

[1]

1968 Stability of periodic waves of ﬁnite amplitude on the surface of a deep ﬂuid

Zakharov VE. 1968 Stability of periodic waves of ﬁnite amplitude on the surface of a deep ﬂuid. Zh. Prikl. Mekh. Tekh. Fiz.9, 86–94. English translationJ. Appl. Mech. Tech. Phys.9, 190–194

work page 1968
[2]

1993 Numerical simulation of gravity waves

Craig W, Sulem C. 1993 Numerical simulation of gravity waves. J. Comp. Phys. 108, 73–83

work page 1993
[3]

1967 A variational principle for a ﬂuid with a free surface

Luke JC. 1967 A variational principle for a ﬂuid with a free surface. J. Fluid Mech. 27, 395–397

work page 1967
[4]

2015 A kinematic conservation law in free surface ﬂow

Gavrilyuk S, Kalisch H, Khorsand Z. 2015 A kinematic conservation law in free surface ﬂow. Nonlinearity 28, 1805–1821

work page 2015
[5]

2015 Well-posedness and shallow-water stability for a new Hamiltonian formulation of the water waves equations with vorticity

Castro A, Lannes D. 2015 Well-posedness and shallow-water stability for a new Hamiltonian formulation of the water waves equations with vorticity. Indiana U. Math. J. 64, 1169–1270

work page 2015
[6]

1958 A theorem on force-free magnetic ﬁelds

Woltjer L. 1958 A theorem on force-free magnetic ﬁelds. Proc. Natl. Acad. Sci. USA 44, 489–491

work page 1958
[7]

1991 On Woltjer’s variational principle for force-free ﬁelds.J

Laurence P , Avellaneda M. 1991 On Woltjer’s variational principle for force-free ﬁelds.J. Math. Phys. 32, 1240–1253

work page 1991
[8]

1984 Impulse, ﬂow-force and variational principles

Benjamin TB. 1984 Impulse, ﬂow-force and variational principles. IMA J. Appl. Math. 32, 3–68

work page 1984

Show all 14 references

[9]

2019 A variational principle for three-dimensional water waves over Beltrami ﬂows

Lokharu E, Wahlén E. 2019 A variational principle for three-dimensional water waves over Beltrami ﬂows. Nonlinear Analysis 184, 193–209

2019
[10]

1988 Mathematical Analysis and Numerical Methods for Science and Technology, Volume 3: Spectral Theory and Applications

Dautray R, Lions JL. 1988 Mathematical Analysis and Numerical Methods for Science and Technology, Volume 3: Spectral Theory and Applications. New York: Springer-Verlag

1988
[11]

2013 The Water Waves Problem: Mathematical Analysis and Asymptotics

Lannes D. 2013 The Water Waves Problem: Mathematical Analysis and Asymptotics . Number 188 in Mathematical Surveys and Monographs. Providence, R.I.: American Mathematical Society

2013
[12]

1975 Sobolev Spaces

Adams RA. 1975 Sobolev Spaces. New York: Academic Press

1975
[13]

1997 Lectures on Nonlinear Hyperbolic Differential Equations

Hörmander L. 1997 Lectures on Nonlinear Hyperbolic Differential Equations . Heidelberg: Springer-Verlag

1997
[14]

2003 Analytic Theory of Global Bifurcation

Buffoni B, Toland JF. 2003 Analytic Theory of Global Bifurcation . Princeton, N. J.: Princeton University Press

2003

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.

[1] [1]

1968 Stability of periodic waves of ﬁnite amplitude on the surface of a deep ﬂuid

Zakharov VE. 1968 Stability of periodic waves of ﬁnite amplitude on the surface of a deep ﬂuid. Zh. Prikl. Mekh. Tekh. Fiz.9, 86–94. English translationJ. Appl. Mech. Tech. Phys.9, 190–194

work page 1968

[2] [2]

1993 Numerical simulation of gravity waves

Craig W, Sulem C. 1993 Numerical simulation of gravity waves. J. Comp. Phys. 108, 73–83

work page 1993

[3] [3]

1967 A variational principle for a ﬂuid with a free surface

Luke JC. 1967 A variational principle for a ﬂuid with a free surface. J. Fluid Mech. 27, 395–397

work page 1967

[4] [4]

2015 A kinematic conservation law in free surface ﬂow

Gavrilyuk S, Kalisch H, Khorsand Z. 2015 A kinematic conservation law in free surface ﬂow. Nonlinearity 28, 1805–1821

work page 2015

[5] [5]

2015 Well-posedness and shallow-water stability for a new Hamiltonian formulation of the water waves equations with vorticity

Castro A, Lannes D. 2015 Well-posedness and shallow-water stability for a new Hamiltonian formulation of the water waves equations with vorticity. Indiana U. Math. J. 64, 1169–1270

work page 2015

[6] [6]

1958 A theorem on force-free magnetic ﬁelds

Woltjer L. 1958 A theorem on force-free magnetic ﬁelds. Proc. Natl. Acad. Sci. USA 44, 489–491

work page 1958

[7] [7]

1991 On Woltjer’s variational principle for force-free ﬁelds.J

Laurence P , Avellaneda M. 1991 On Woltjer’s variational principle for force-free ﬁelds.J. Math. Phys. 32, 1240–1253

work page 1991

[8] [8]

1984 Impulse, ﬂow-force and variational principles

Benjamin TB. 1984 Impulse, ﬂow-force and variational principles. IMA J. Appl. Math. 32, 3–68

work page 1984