{"id":"4ddc395e-d66a-427f-8f66-ede97d03b3e8","arxiv_id":"1908.09386","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"Steady surface water waves on a Beltrami flow are equivalent to two scalar equations for the surface elevation and a surface potential, derived from a variational principle via a nonlocal operator that generalizes the Dirichlet-Neumann operator.","lead":"This paper rewrites the equations for steady water waves that ride a Beltrami flow as two equations for two scalar surface functions, generalizing the classical Zakharov-Craig-Sulem formulation. The new variational principle provides a compact mathematical framework that could support Hamiltonian and numerical treatments of three-dimensional waves with constant vorticity.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 3.1's dynamic equation contains a spurious term; the printed two-scalar system is not equivalent to (1.7)-(1.11).","rationale":"The reader's weakest_assumption concerned the small-|α|, small-η well-posedness of the operator H(η). That is a legitimate limitation and the paper states it. However, a more immediate correctness risk appears in the main theorem: the dynamic Euler-Lagrange equation as displayed in Theorem 3.1 does not follow from the equations derived in Section 2. A direct computation from (2.5)-(2.6) eliminates the nonlocal α-term and yields a dynamic equation without the printed −α H(H+K·∇η)/(1+|∇η|^2) correction. This is not a matter of scope or smallness but of internal algebraic consistency. The variational principle itself may still be salvageable, since Section 2 derives correct equations, but the claimed compact two-equation reduction, which is the paper's headline result, is misstated. The appropriate remedy is to correct or remove the spurious term and re-verify the subsequent conclusions. This warrants a conditional rather than unconditional acceptance, pending that correction.","tokens_in":21713,"tokens_out":31092,"duration_ms":276469,"concrete_test":"Independently recompute the substitution from (2.5)-(2.6) into the variables H(η)Φ and K(η)Φ, using v_‖ = K, v·N = H, v₂ = (H+K·∇η)/(1+|∇η|^2), and u⋆·N = −H. The correct reduction should produce 1/2|K|^2 − (H+K·∇η)^2/(2(1+|∇η|^2)) + K·u⋆_h + gη − σ(...) = 0, with no term proportional to α. A symbolic-algebra verification of this elimination, or a numerical evaluation on one smooth test pair (η,Φ) satisfying (2.5), will settle whether the α-term in Theorem 3.1 is spurious.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's central reduction claims that (1.7)-(1.11) is equivalent to the two equations H(η)Φ + u⋆·N = 0 and the dynamic equation displayed in Theorem 3.1. That dynamic equation contains the term −α H(H+K·∇η)/(1+|∇η|^2), which, according to a direct substitution from the paper's own Euler-Lagrange equations (2.5)-(2.6), should not survive. Indeed, after imposing (2.5), the terms (curl A)₂(−∇·A^⊥_‖ + ∇η·u⋆_‖) and α∇^⊥Δ⁻¹(∇·A^⊥_‖ + ∇·A⋆^⊥_‖)·u⋆_‖ both vanish, leaving 1/2|v|^2 + v_h·u⋆_h + gη − σ(...) = 0. Writing v_‖ = K, v·N = H, v₂ = (H+K·∇η)/(1+|∇η|^2), v_h = K − v₂∇η, and using u⋆·N = −H (from (2.5)) gives 1/2|v|^2 + v_h·u⋆_h = 1/2|K|^2 − (H+K·∇η)^2/(2(1+|∇η|^2)) + K·u⋆_h, with no α-term. The displayed −α H(H+K·∇η)/(1+|∇η|^2) is therefore inconsistent with (2.5)-(2.6) and the identities stated in the proof of Theorem 3.1. As printed, the two-scalar system is not equivalent to the original free-boundary problem, even though the variational derivation in Section 2 may be correct.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","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 η.","tokens_in":22106,"tokens_out":16754,"duration_ms":144282,"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":[{"comment":"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.","section":"Theorem 3.1 and §2, Eq. (2.6)"},{"comment":"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).","section":"Introduction, Zakharov-Craig-Sulem limit"}],"minor_comments":[{"comment":"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.","section":"§4(c), Proposition 4.7"},{"comment":"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.","section":"Abstract and §4"},{"comment":"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.","section":"§4(a)"},{"comment":"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.","section":"Theorem 3.1"}],"recommendation":"major_revision","confidential_remarks":"The spurious α in Theorem 3.1 looks like a typographical or rewriting error rather than a flaw in the Section 2 variational derivation; nevertheless it affects the main theorem as printed. The sign inconsistency in the α = 0 limit should also be resolved. Both issues are fixable within the manuscript's scope, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","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.\n\nBut 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.\n\nThe 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.\n\nThe 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.","headline":"Correct variational derivation, but Theorem 3.1 as printed has a spurious α-term that breaks the stated equivalence; easy fix, still worth publishing.","tokens_in":22600,"tokens_out":6440,"would_cite":true,"duration_ms":47932,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["76B15","35Q35"],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["Beltrami flows","steady water waves","variational principle","Dirichlet-Neumann operator","Hodge-Weyl decomposition","nonlocal operator","free-boundary problem"],"falsifier":"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.","tokens_in":21538,"feed_emoji":"🌊","tokens_out":7096,"duration_ms":64400,"temperature":0.7,"pith_summary":"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.","feed_headline":"Beltrami water waves collapse to two scalar surface equations","feed_subtitle":"A nonlocal operator H(η) turns the 3-D free-boundary problem into a pair of equations; α=0 recovers the classical wave theory.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the classical surface-variable formulation that the Beltrami result must reduce to when α = 0.","marker":"[1]"},{"why":"Supplies the Dirichlet–Neumann-operator form of the irrotational equations that H(η) generalises.","marker":"[2]"},{"why":"Supplies the variational principle for irrotational water waves whose Beltrami analogue is constructed here.","marker":"[3]"},{"why":"Provides the Hamiltonian formulation with vorticity and the characterisation of the space Xη used in the well-posedness proof.","marker":"[5]"},{"why":"Provides the classical variational principle for force-free fields underlying the bulk term of the functional L.","marker":"[6]"},{"why":"Gives an alternative vector-potential variational principle for Beltrami water waves that the present formulation refines and compares with.","marker":"[9]"},{"why":"Supplies the elliptic boundary-value estimates used for the flattened problem in Section 4.","marker":"[11]"},{"why":"Supplies the analytic implicit-function theorem used to prove analyticity of H(η).","marker":"[14]"}],"fun_headline_variants":["Beltrami surface waves reduce to a pair of scalar equations","Variational principle yields two equations for Beltrami water waves","Nonlocal operator H(η) yields two-field water-wave equations","Two scalar unknowns describe steady waves on Beltrami flow","Variational reduction of Beltrami water waves to two surface equations"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Beltrami surface waves reduce to a pair of scalar equations","Variational principle yields two equations for Beltrami water waves","Nonlocal operator H(η) yields two-field water-wave equations","Two scalar unknowns describe steady waves on Beltrami flow","Variational reduction of Beltrami water waves to two surface equations"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001121,"raw_usage":{"total_tokens":4657,"prompt_tokens":929,"completion_tokens":3728,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":545,"completion_tokens_details":{"reasoning_tokens":3644}},"tokens_in":545,"tokens_out":3728,"duration_ms":28436,"temperature":1.0,"reasoning_tokens":3644,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:12:51.484130+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"1968 Stability of periodic waves of finite amplitude on the surface of a deep fluid","cited_arxiv_id":null,"evidence_quote":"Supplies the classical surface-variable formulation that the Beltrami result must reduce to when α = 0."},{"cited_title":"1993 Numerical simulation of gravity waves","cited_arxiv_id":null,"evidence_quote":"Supplies the Dirichlet–Neumann-operator form of the irrotational equations that H(η) generalises."},{"cited_title":"1967 A variational principle for a fluid with a free surface","cited_arxiv_id":null,"evidence_quote":"Supplies the variational principle for irrotational water waves whose Beltrami analogue is constructed here."},{"cited_title":"2015 Well-posedness and shallow-water stability for a new Hamiltonian formulation of the water waves equations with vorticity","cited_arxiv_id":null,"evidence_quote":"Provides the Hamiltonian formulation with vorticity and the characterisation of the space Xη used in the well-posedness proof."},{"cited_title":"1958 A theorem on force-free magnetic fields","cited_arxiv_id":null,"evidence_quote":"Provides the classical variational principle for force-free fields underlying the bulk term of the functional L."},{"cited_title":"2019 A variational principle for three-dimensional water waves over Beltrami flows","cited_arxiv_id":null,"evidence_quote":"Gives an alternative vector-potential variational principle for Beltrami water waves that the present formulation refines and compares with."},{"cited_title":"2013 The Water Waves Problem: Mathematical Analysis and Asymptotics","cited_arxiv_id":null,"evidence_quote":"Supplies the elliptic boundary-value estimates used for the flattened problem in Section 4."},{"cited_title":"2003 Analytic Theory of Global Bifurcation","cited_arxiv_id":null,"evidence_quote":"Supplies the analytic implicit-function theorem used to prove analyticity of H(η)."}],"review_version":1}