arxiv: 2605.01751 · v1 · submitted 2026-05-03 · ⚛️ physics.flu-dyn

Recognition: unknown

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

Semyon Churilov

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 comments

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: Sketch of the flow configuration in the vertical plane. The view at source ↗

**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: Solution of equation (39): line 1 shows the left-hand side, a view at source ↗

**Figure 4.** Figure 4: The subcritical part of the phase portrait of equation (2 view at source ↗

**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: 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: 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: 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 reflectionless 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.

Desk Editor's Note private letter to a colleague

Churilov applies the Laplace cascade to generate specific reflectionless profiles for linear shallow-water waves, but only when bathymetry and current take integrable functional forms.

read the letter

The paper's main move is to take the Laplace cascade method for second-order hyperbolic equations and apply it to the linear long-wave equation that comes from shallow-water theory with variable depth h(x) and current u(x). It produces an algorithm that locates parameter values yielding traveling-wave solutions with no reflection, then works out a handful of explicit examples and lines them up against earlier results in the literature. That is the concrete advance: a systematic way to build those special profiles rather than just guessing at them. The examples are the part that actually lands, because they show what kinds of depth and flow variations satisfy the integrability chain needed for the cascade to terminate with a reflectionless solution. The comparisons give a reader a quick sense of where this sits relative to prior work on the same problem. The limitation is straightforward. The cascade only closes when the coefficient functions obey a sequence of auxiliary ODEs, so the admissible h(x) and u(x) are restricted to particular families (exponentials, powers, and similar). The abstract's claim of a general algorithm therefore holds only inside that integrable class; it does not deliver a method that starts from an arbitrary given bathymetry and current and either finds or rules out reflectionless flow. That scope issue is real but not fatal for the paper's stated purpose. The work is aimed at people who already care about exact solutions for wave equations in inhomogeneous media, whether in coastal modeling or in related hyperbolic problems. A reader who needs tools for arbitrary real-world bathymetry will find it too narrow. It is worth sending to referees because the mathematical construction is reproducible, the examples are checkable, and the comparisons are present; reviewers can then press on the exact range of applicability and ask for any missing derivation steps or numerical spot-checks.

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)

[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.
[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)

All equations in the derivation should be numbered consecutively to facilitate reference when discussing the cascade steps and auxiliary conditions.
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

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
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)

Reference graph

Works this paper leans on

21 extracted references

[1]

1986 Catastrophe Theory, 2nd edn

Arnold, V.I. 1986 Catastrophe Theory, 2nd edn. Springer

1986
[2]

2023 a Traveling internal waves in a two-layer shallow medium with variable bathymetry and current

Churilov, S.M. 2023 a Traveling internal waves in a two-layer shallow medium with variable bathymetry and current. Phys. Fluids 35, 026602

2023
[3]

2023 b Traveling internal waves in a two-layer shallow medium with variable bathymetry and current: Surface and near-bottom ﬂows

Churilov, S.M. 2023 b Traveling internal waves in a two-layer shallow medium with variable bathymetry and current: Surface and near-bottom ﬂows. Phys. Fluids 35, 076613

2023
[4]

2024 Traveling Alfv´ en waves in inhomogeneous plasm a ﬂows along mag- netic ﬁeld

Churilov, S.M. 2024 Traveling Alfv´ en waves in inhomogeneous plasm a ﬂows along mag- netic ﬁeld. Phys. Fluids 36, 087166

2024
[5]

& Stepanyants, Yu.A

Churilov, S.M. & Stepanyants, Yu.A. 2022 a Reﬂectionless wave propagation on shallow water with variable bathymetry and current, J. Fluid Mech. 931, A15 (Paper I)

2022
[6]

& Stepanyants, Yu.A

Churilov, S.M. & Stepanyants, Yu.A. 2022 b Reﬂectionless wave propagation on shallow water with variable bathymetry and current. Part 2, J. Fluid Mech. 939, A15 (Paper II)

2022
[7]

1915 Lec¸ ons sur la th´ eorie g´ en´ erale des surfaces et les appli cations g´ eom´ etriques du calcul inﬁnit´ esimal, tome 2

Darboux G. 1915 Lec¸ ons sur la th´ eorie g´ en´ erale des surfaces et les appli cations g´ eom´ etriques du calcul inﬁnit´ esimal, tome 2. Gautier-Villars

1915
[8]

& Fr¨ oman, P.O

Fr¨ oman, N. & Fr¨ oman, P.O. 1965JWKB Approximation . North-Holland
[9]

1926 Le¸ cons sur l’int´ egration des ´ equations aux d´ eriv´ ees partielles du second ordre a deux variables ind´ ependantes, tome 2

Goursat, ´E. 1926 Le¸ cons sur l’int´ egration des ´ equations aux d´ eriv´ ees partielles du second ordre a deux variables ind´ ependantes, tome 2. Libraire Scientiﬁque J. Hermann

1926
[10]

& Pelinovsky, E

Grimshaw, R., Pelinovsky, D. & Pelinovsky, E. 2010 a Homogenization of variable-speed wave equation. Wave Motion 47, 496–507

2010
[11]

& Talipova, T

Grimshaw, R., Pelinovsky, E. & Talipova, T. 2010 b Nonreﬂecting internal wave beam propagation in the deep ocean. J. Phys. Oceanogr. 40, 802–813

2010
[12]

1962 An Introduction to Phase-Integral Methods

Heading, J. 1962 An Introduction to Phase-Integral Methods . Wiley

1962
[13]

1916 Integration of diﬀerential equations with partial derivat ives of the 1st and 2nd orders

Imshenetsky, V.G. 1916 Integration of diﬀerential equations with partial derivat ives of the 1st and 2nd orders . Moscow Mathematical Society (in Russian)

1916
[14]

& Kaptsov, D.O

Kaptsov, O.V. & Kaptsov, D.O. 2023 Solutions of some wave mech anics models. Fluid Dynamics 58, 1227–1234

2023
[15]

2005 Methods for Constructing Exact Solutions of Partial Diﬀere ntial Equations

Meleshko, S.V. 2005 Methods for Constructing Exact Solutions of Partial Diﬀere ntial Equations. Springer

2005
[16]

1958 Mathematical Theory of Compressible Fluid Flow

Mises, R. 1958 Mathematical Theory of Compressible Fluid Flow . Academic. 27

1958
[17]

& Kaptsov, O

Pelinovsky, E. & Kaptsov, O. 2022 Traveling waves in shallow seas of variable depth. Symmetry 14, 1448

2022
[18]

& Didenkulova, E

Pelinovsky, E., Talipova, T., Didenkulova, I. & Didenkulova, E. 201 9 Interfacial long traveling waves in a two-layer ﬂuid with variable depth. Stud. Appl. Math. 142, 513–527
[19]

& Didenkulova, E.G

Petrukhin, N.S., Pelinovsky, E.N. & Didenkulova, E.G. 2020 Non-re ﬂective magnetohy- drodynamic waves in an inhomogeneous plasma. Radiophysics and Quantum Electronics 63, 29–39

2020
[20]

& Talipova, T.G

Petrukhin, N.S., Pelinovsky, E.N. & Talipova, T.G. 2023 Non-reﬂec ted acoustic waves in an inhomogeneous atmosphere. Radiophysics and Quantum Electronics 66, 431–440

2023
[21]

& Zaitsev, V.F

Polyanin, A.D. & Zaitsev, V.F. 2003 Handbook of Exact Solutions for Ordinary Diﬀerential Equations, 2nd edn. CRC Press. 28

2003