pith. sign in

arxiv: 2605.01751 · v2 · pith:CXSEE3V7new · submitted 2026-05-03 · ⚛️ physics.flu-dyn

Traveling surface wave propagation on shallow water with variable bathymetry and current

Semyon Churilov This is my paper

Pith reviewed 2026-05-09 16:44 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords reflectionless wavesshallow water wavesvariable bathymetryLaplace cascade methodsurface wave propagationinhomogeneous flowschannel with variable cross-sectionlong linear waves
0
0 comments X

The pith

The Laplace cascade method yields a general algorithm to locate parameters for reflectionless long surface waves in shallow channels with variable bathymetry and current.

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

The paper develops a systematic algorithm based on the Laplace cascade method for integrating second-order hyperbolic equations to identify conditions in inhomogeneous media that permit waves to propagate with virtually no reflection. This is applied specifically to long linear surface waves traveling in a channel whose cross-section varies with position. A reader would care because such arrangements would allow wave energy to travel long distances efficiently, which is relevant to natural coastal processes and to engineering efforts that seek to limit wave damage. The work illustrates the algorithm with several concrete solutions and places them against earlier results on the same problem.

Core claim

By applying the Laplace cascade method to the governing wave equation, the paper constructs a general procedure that finds the spatial profiles of bathymetry and current for which the inhomogeneous shallow-water system supports exact traveling-wave solutions without reflection or scattering. The procedure is then specialized to the case of long linear surface waves in a channel of variable cross-section, producing explicit families of reflectionless flows that are compared with previously known examples.

What carries the argument

The Laplace cascade method applied to the second-order hyperbolic wave equation, which systematically generates the variable coefficients (bathymetry and current) that admit reflectionless traveling solutions.

If this is right

  • Explicit profiles of depth and current exist that support long-distance transmission of wave energy with negligible scattering.
  • The same algorithm applies to other inhomogeneous hyperbolic systems that admit reflectionless solutions.
  • Coastal channels can in principle be shaped so that incoming waves do not reflect and therefore do not amplify local loads on structures.
  • The representative solutions recover or extend known special cases from earlier literature on the same problem.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Laboratory flume experiments with controlled bottom topography and mean flow could directly test whether the predicted profiles indeed produce near-zero reflection.
  • The method could be adapted to acoustic or electromagnetic waveguides with varying cross-sections to engineer reflectionless propagation in those domains.
  • If the linear solutions remain approximately reflectionless when weak nonlinearity is restored, the approach would offer a starting point for designing low-reflection coastal defenses.

Load-bearing premise

The physical setup of variable bathymetry and current must produce a wave equation whose coefficients can be arranged so that the Laplace cascade method finds exact reflectionless traveling-wave solutions.

What would settle it

Take one of the explicit bathymetry-current profiles produced by the algorithm, solve the linear wave equation numerically for an incident wave packet, and measure whether the reflected amplitude is numerically zero across a range of wave frequencies.

Figures

Figures reproduced from arXiv: 2605.01751 by Semyon Churilov.

Figure 1
Figure 1. Figure 1: Sketch of the flow configuration in the vertical plane. The view at source ↗
Figure 2
Figure 2. Figure 2: Solutions to equation (25) for (a) c0 = 0.7 and (b) c0 = 2; (c) solution to equation (27) for c0 = 2. Subcritical branches are marked by 1, supercritical ones – by 2. Here, the flow velocity vanishes at ξ = ξ0 as well. The left-hand side of (27) has a minimum at a = c0, so at the critical point ξ = ξc = ξ0 − 1 3 + 1 c 2 0 + c 4 0 − 1 c 3 0 arctan c0 > ξ0 there is a fold-type singularity, and the solution i… view at source ↗
Figure 3
Figure 3. Figure 3: Solution of equation (39): line 1 shows the left-hand side, a view at source ↗
Figure 4
Figure 4. Figure 4: The subcritical part of the phase portrait of equation (2 view at source ↗
Figure 5
Figure 5. Figure 5: The plots of E(τ ) (black lines 1a, b, c) and 1/E(τ ) (blue lines 2a, b, c) for (a) τ∗ = 0.4 and (b) τ∗ = 1.6. Negative branches are shown by dashes. figure 2), inherent to flows of rank 1, arises because upon passing this point, the dependence of the coordinate x on the parameter τ ceases to be monotonic (see the first formula (47) and equation (42)), dτ = −Cds = C c dx c 2 − U2 ≡ c dξ c 2 − U2 ≡ c 3dξ c … view at source ↗
Figure 6
Figure 6. Figure 6: Currents without critical points: (a) – flows of class Bˆ+, c0 = 0.5, (b) – flows of class Cˆ+, c0 = 2; τ∗ = 1.6. 17 view at source ↗
Figure 7
Figure 7. Figure 7: Currents with two critical points: (a) – flows of class Bˆ+, (b) – flows of class Cˆ+; τ∗ = 0.4, c0 = 1. 4.1.2 Global flows The structure of the solution (47) dictates a somewhat different approach to the existence of global flows than that adopted in §§ 3.2.2 and 3.2.3. As already noted, the domain of definition of each flow described by such a solution lies in one of three intervals of τ and is bounded b… view at source ↗
Figure 8
Figure 8. Figure 8: Solutions to equations (a) (56) and (b) (62). Curves 1 correspond to subcritical branches, and curves 2 – to supercritical ones. the left-hand side of which has a minimum at the critical point a = c0 (ξa = ξa0 + 8 7 c 7/2 0 ). The plot a(ξa) is shown in figure 8 (a). To consider the global flows, it is convenient to introduce the functions G(ξa) = " a(ξa) c(ξa) #1/4 and g(ξa) = 1 G(ξa) and represent the se… view at source ↗
read the original abstract

Energy transmission over long distances by waves is a key mechanism for many natural processes. This possibility arises when an inhomogeneous medium is arranged in such a manner that it enables a certain type of wave to propagate with virtually no reflection or scattering. By application of the Laplace cascade method for integrating second-order hyperbolic equations, a general algorithm for finding the parameters of inhomogeneous non-reflecting flows is proposed. The algorithm is applied to the problem of long linear surface waves propagation in a channel with variable cross-section. The general analysis of the problem is illustrated by a few representative solutions and compared with the results of previous studies. The results obtained may be of interest to mitigate the possible impact of waves on ships, marine engineering constructions, and human coastal activities.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper claims that application of the Laplace cascade method for integrating second-order hyperbolic equations yields a general algorithm for determining the parameters of inhomogeneous reflectionless flows. This algorithm is applied to the linear long-wave equation for surface waves propagating in shallow water with variable bathymetry and current (variable channel cross-section), with the analysis illustrated by representative solutions that are compared to results from prior studies.

Significance. If the algorithm correctly identifies reflectionless configurations, the work would provide a systematic mathematical procedure for constructing special classes of bathymetry and current profiles that permit traveling waves without reflection. This could be useful for understanding energy transmission in inhomogeneous media and for applications in coastal engineering aimed at reducing wave impacts on structures and vessels.

major comments (2)
  1. [Abstract] Abstract: The description of the method as a 'general algorithm for finding the parameters of inhomogeneous reflectionless flows' applied to the wave problem overstates the scope. The Laplace cascade requires the potential term in the canonical hyperbolic form to satisfy a sequence of auxiliary ODEs, so the procedure constructs only special integrable profiles (e.g., exponential or power-law families) rather than solving the inverse problem for arbitrary given h(x) and u(x). The paper must explicitly state the restricted class of bathymetry/current profiles for which the algorithm succeeds.
  2. [Section presenting the general analysis and algorithm] Section presenting the general analysis and algorithm: No explicit derivation steps are supplied showing how the linear long-wave equation is reduced to canonical form, how the cascade integrations are performed, or how the resulting solutions are verified to produce zero reflection over the full parameter range. Without these steps, error analysis, or direct checks against the original wave equation, the central claim that the constructed flows are reflectionless cannot be assessed.
minor comments (2)
  1. All equations in the derivation should be numbered consecutively to facilitate reference when discussing the cascade steps and auxiliary conditions.
  2. The comparison with previous studies would benefit from a table summarizing the specific bathymetry/current profiles tested and the quantitative measures of reflection (or lack thereof) obtained.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review of our manuscript. The comments have helped us identify areas where the presentation can be clarified and strengthened. We address each major comment below and have revised the manuscript to incorporate the suggested improvements.

read point-by-point responses

  1. Referee: [Abstract] Abstract: The description of the method as a 'general algorithm for finding the parameters of inhomogeneous reflectionless flows' applied to the wave problem overstates the scope. The Laplace cascade requires the potential term in the canonical hyperbolic form to satisfy a sequence of auxiliary ODEs, so the procedure constructs only special integrable profiles (e.g., exponential or power-law families) rather than solving the inverse problem for arbitrary given h(x) and u(x). The paper must explicitly state the restricted class of bathymetry/current profiles for which the algorithm succeeds.

    Authors: We agree that the original abstract could be read as implying broader applicability than the method provides. The Laplace cascade method yields reflectionless solutions only for those profiles where the potential term satisfies the required sequence of auxiliary ODEs, resulting in specific integrable families. In the revised manuscript we have updated the abstract to state explicitly that the algorithm identifies parameters of reflectionless flows for special classes of bathymetry and current profiles (such as exponential and power-law families) rather than for arbitrary prescribed h(x) and u(x). This revision removes any overstatement of scope while preserving the description of the algorithm's utility within its applicable domain. revision: yes

  2. Referee: [Section presenting the general analysis and algorithm] Section presenting the general analysis and algorithm: No explicit derivation steps are supplied showing how the linear long-wave equation is reduced to canonical form, how the cascade integrations are performed, or how the resulting solutions are verified to produce zero reflection over the full parameter range. Without these steps, error analysis, or direct checks against the original wave equation, the central claim that the constructed flows are reflectionless cannot be assessed.

    Authors: We acknowledge that the original manuscript presented the reduction and verification steps in a condensed manner. To address this concern we have expanded the relevant section with a new subsection that provides the explicit sequence of transformations: (i) reduction of the linear long-wave equation to canonical hyperbolic form, (ii) the successive integrations performed via the Laplace cascade, and (iii) direct substitution of the resulting solutions back into the original wave equation to confirm that the reflection coefficient vanishes identically over the admissible parameter range. We have also added a brief consistency check against known analytic solutions from the literature and a short discussion of the error bounds inherent to the long-wave approximation. These additions allow the reader to follow and verify the reflectionless property without ambiguity. revision: yes

Circularity Check

0 steps flagged

No circularity detected; derivation relies on external Laplace cascade method applied to derived wave equation.

full rationale

The paper proposes an algorithm based on the established Laplace cascade method for hyperbolic PDEs to identify special bathymetry/current profiles permitting reflectionless long-wave propagation. This is an application of a known external mathematical technique to the shallow-water wave equation (derived from standard linearization), with representative solutions compared to prior literature. No step reduces a claimed prediction or result to a fitted input or self-citation by construction; the central claim is the systematic location of integrable cases rather than a tautological redefinition of inputs. The method's restriction to special functional forms is a scope limitation, not circularity.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work rests on standard assumptions of linear shallow-water wave theory and the applicability of the Laplace cascade method to the resulting hyperbolic system; no new free parameters or invented entities are introduced in the abstract.

axioms (2)
  • domain assumption Long linear surface waves on shallow water obey a second-order hyperbolic PDE system.
    Standard modeling assumption for the regime stated in the title and abstract.
  • domain assumption Inhomogeneous media can be arranged to support reflectionless traveling waves.
    Core premise required for the algorithm to produce useful solutions.

pith-pipeline@v0.9.0 · 5415 in / 1266 out tokens · 33509 ms · 2026-05-09T16:44:55.234353+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.