{"id":"52b11f41-3486-4f86-9288-193784117db3","arxiv_id":"2412.05034","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For the Benney-Luke shallow-water model, Whitham modulation theory gives formulas for Mach reflection and expansion of solitons and predicts circular precursor wavefronts for moving topography.","lead":"This paper derives equations that describe how the shape and amplitude of a solitary ocean wave change when it bends or reflects, using a model that allows waves to travel in all directions. The authors predict the size of the Mach stem formed at the bend and show that waves generated by moving underwater topography spread in circles, unlike the parabolic fronts predicted by an older model.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equations (3.10b,c) use a reciprocal factor (1+a) that is inconsistent with the dispersion relation and Appendix A, and this error propagates into the Rankine–Hugoniot conditions (5.7) used for the Mach reflection predictions.","rationale":"The reader's weakest assumption concerned the unjustified decomposition of (3.10c) into two one-dimensional conservation laws. My stress-test found a more basic and more specific problem: the coefficient in (3.10b,c) appears to be the reciprocal of the coefficient that follows from the stated dispersion relation and compatibility conditions. This is not a stylistic issue; it changes the characteristic speeds, the Riemann invariants, and every downstream conservation-law statement, including the RH conditions (5.7) that determine the Mach stem and reflected-soliton amplitudes. If the factor error is real, the central analytical claim of the paper is unsupported, although the numerical results might still be correct if the code used the correct equations. Because the manuscript contains an internal inconsistency between the main text and Appendix A, the paper cannot be accepted without correction. I recommend a conditional acceptance, with the condition that the factor error be resolved and the numerical comparisons be recomputed or confirmed from the corrected equations. This is consistent with the reader's conditional verdict, but the required revision is more specific than the decomposition justification: the modulation equations themselves must be fixed and the RH predictions re-derived. Therefore I mark agreement as partial: the reader identified the decomposition as a weak point, and my concern lies one level deeper in the same chain of reasoning.","tokens_in":29443,"tokens_out":14726,"duration_ms":122230,"concrete_test":"Independently re-derive (3.10c) from (3.5)–(3.8) and confirm whether the coefficient is √(1+a)√(1+q²) or √(1+q²)/√(1+a). Then recompute the RH system (5.4)-(5.7) with the corrected B for the case a0=0.21, θ0=20° and compare the resulting (a_w, a_i, q_i) against the quoted values 0.478149, 0.0568022, and 0.797368. If the corrected system yields different values, the analytical predictions and the claimed numerical agreement must be re-examined; if the corrected system reproduces the quoted numbers, the error is typographical and can be fixed without changing the conclusions.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central analytic result for Mach reflection rests on the modulation system (3.10). Re-deriving (3.10c) from the compatibility conditions (3.5) and the soliton dispersion relation (3.8), with ω = −k√(1+a)√(1+q²) for left-going waves, gives q_t + q B_x − q_x B − B_y = 0 with B = √(1+a)√(1+q²). The published (3.10c) instead has B′ = √(1+q²)/√(1+a), i.e. a reciprocal factor (1+a). The same reciprocal error appears in (3.10b), where the correct flux is kB = k√(1+a)√(1+q²). Appendix A (A 19)–(A 20) are consistent with B, not B′, so the main text is internally inconsistent. Because the decomposition (5.5) and the classical Rankine–Hugoniot conditions (5.7) are built on (3.10c), the predicted amplitudes (a_w ≈ 0.478149, a_i ≈ 0.0568022) and the critical slope (5.10) are not supported by the printed equations. The reader flagged the decomposition of (3.10c) as unjustified; even if that decomposition were accepted, the base equation itself appears to contain a factor error, making the analytical predictions unreliable as written. The numerical agreement reported in Figure 7 may indicate the actual code used the correct B, but the manuscript must be corrected and re-verified.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a Whitham modulation theory for the two-dimensional Benney–Luke equation and uses it to analyze Mach expansion and Mach reflection of obliquely interacting line solitons, deriving critical angles and Mach-stem amplitudes that are compared with direct numerical simulations of the BL equation and with the KP equation. The paper also treats the forced BL equation for topography-generated waves and claims that the far-field precursor wavefronts are circular, in contrast to the parabolic wavefronts of the forced KP equation.","tokens_in":29814,"tokens_out":13067,"duration_ms":205664,"significance":"If the modulation system and its downstream predictions are correct, the paper would provide a useful isotropic, bidirectional counterpart to the well-developed KP modulation theory for soliton interactions, with explicit formulas for critical slopes and amplitudes and with direct BL numerics supporting the predictions. The paper is self-contained in the sense that no free parameters are fitted to the numerical solutions, the derivations are presented in detail, and the comparison with the KP limit near unit speed is a valuable check. However, the central modulation system contains a factor inconsistency that propagates into most of the analytical results, and the reduced Rankine–Hugoniot construction in Section 5 relies on an unproved decomposition. These issues make the analytical claims, as printed, unsupported by the derivation.","major_comments":[{"comment":"The fluxes in (3.10b,c) are sqrt(1+q^2)/sqrt(1+a), but the dispersion relation (3.8), omega = -k sqrt(1+a) sqrt(1+q^2) for left-going waves, together with the compatibility conditions (3.5), yields k_t - (kB)_x = 0 and q_t + q B_x - q_x B - B_y = 0 with B = sqrt(1+a) sqrt(1+q^2). This is exactly what Appendix A states in (A19)–(A20). The printed (3.10b,c) therefore contain a reciprocal factor (1+a). Because (3.10c) feeds into the characteristic speeds (4.3), the Riemann invariants (4.4), the Mach expansion formulas (4.8) and (4.20), and the Rankine–Hugoniot system (5.4)–(5.7), the analytical predictions in Sections 4 and 5 are not supported by the printed equations. Moreover, the reduced 1D system (4.2) is not consistent with either version of (3.10c): using B' gives a_y coefficient +sqrt(1+q^2)/(2(1+a)^{3/2}), whereas using B gives q_y coefficient sqrt(1+a) q/sqrt(1+q^2). The manuscript must be corrected and the numerical comparisons re-verified against the corrected system.","section":"§3, Eqs. (3.10b,c) and Appendix A (A19)–(A20)"},{"comment":"The decomposition of (3.10c) into two independent one-dimensional conservation laws (5.5) by 'alternately neglecting the x- and y-directions' is asserted without a derivation or a controlled asymptotic justification. This step is load-bearing because the classical Rankine–Hugoniot conditions (5.7) obtained from it determine the Mach stem amplitude a_w and the reflected amplitude a_i. In addition, the quantity F^(1)=1/q in (5.6) is singular at the Mach stem value q_w=0, so the first equality in (5.7) is undefined as written and no limiting procedure is stated. Even if the decomposition were accepted, the factor error in (3.10c) changes the fluxes in (5.6) and hence the algebraic system that produces a_w≈0.478149 and a_i≈0.0568022.","section":"§5.1, Eqs. (5.5)–(5.7)"},{"comment":"The circular similarity solution (6.2) is introduced with the assumptions 'k = 1 and a independent of y', but no derivation is given and the stated relation a = f^2/t - 1 does not appear to satisfy the modulation equations (4.2) for general f(t); for example, with f(t)=t the first equation of (4.2) is not satisfied identically. Since the circular (rather than parabolic) shape of the precursor wavefronts is a central claim of the paper, the far-field reduction should either be derived carefully from the corrected modulation system or explicitly presented as a numerically motivated ansatz with its range of validity stated.","section":"§6.1–6.2, Eq. (6.2)"}],"minor_comments":[{"comment":"The expansion variable eps = -omega_0 - 1 is introduced in (4.20) and (5.10) using the symbol eps, which is already used for the small parameter in the multiple-scales rescaling (3.2); please use a distinct symbol to avoid confusion.","section":"§4.3, Eq. (4.20)"},{"comment":"The captions refer to 'black dashed lines' for the linear wake-angle predictions, while the figures show black dotted lines; please unify the terminology.","section":"Figures 10 and 13"},{"comment":"The reduction from the fBL equation (2.1) to the fKP equation (2.3)–(2.4) appears to have a coefficient in the nonlinear term (3/2 versus the standard value after the rescaling (2.5)) that is not explained; please check the derivation and state the scaling conventions explicitly.","section":"§2.1, Eq. (2.3)"},{"comment":"The phase parameters kappa_1, kappa_2, kappa_3 are used to express the KP soliton amplitudes and slopes, but their normalization relative to the line soliton solution (2.9) is not defined; please add the precise definition or a reference.","section":"§5.3, Eq. (5.13)"}],"recommendation":"major_revision","confidential_remarks":"The factor inconsistency in (3.10b,c) is serious because it propagates into the paper's central analytical results in Sections 4 and 5. The reported agreement with direct BL numerics in Figures 3 and 7 cannot be assessed against the printed theory, and the same inconsistency affects the Rankine–Hugoniot system and the singular 1/q flux in (5.6). I recommend that the authors re-derive the modulation system from the dispersion relation and Appendix A, re-derive all downstream formulas, and re-run the numerical comparisons before the paper is reconsidered."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Zhan, this one is worth a serious look, but not as it stands. The paper does something that hasn't been done before: it derives Whitham modulation equations for the Benney–Luke equation and uses them to get Riemann invariants, critical angles, and a circular similarity solution for wave wakes. The numerical work is careful and the comparison to KP is thoughtful. Section 6's circular precursor prediction is a genuinely testable consequence.\n\nThe soft spot is the algebra in the modulation system. I re-derived (3.10c) from the compatibility conditions and the dispersion relation (3.8) with ω = −k√(1+a)√(1+q²). The correct flux is B = √(1+a)√(1+q²), giving q_t + qB_x − q_xB − B_y = 0. The printed (3.10c) uses B′ = √(1+q²)/√(1+a). That's a reciprocal factor. The same error appears in the classical RH conditions: in (5.6), G^(1) correctly uses B, but H^(2) uses B′. So the paper is internally inconsistent, and since (5.7) and the predicted amplitudes a_w ≈ 0.478, a_i ≈ 0.057 lean on H^(2), those numbers are not supported by the printed equations. The stress-test note says Appendix A is consistent with B; it isn't—the appendix uses yet another reciprocal, √(1+a)/√(1+q²). That makes the mess broader, not smaller. The numerical agreement in Figures 3 and 7 suggests the code used the correct equations and the wrong ones are typos, but the manuscript has to be corrected and re-verified before the claims stand.\n\nThe other soft spot is the decomposition (5.5) in Section 5.1: splitting (3.10c) into two one-dimensional conservation laws by alternately dropping x and y is asserted, not derived. The reader flagged this, and it's fair. If the decomposition is not valid, the algebraic system is underdetermined. The paper also doesn't give quantitative error measures or convergence studies, which is a minor but real omission for a numerical-theory comparison.\n\nThat said, the overall program is coherent and the numerical experiments are substantial. This deserves peer review, but it needs major revision first: correct (3.10b,c) and (5.6), re-derive the RH predictions, and justify the decomposition. I'd cite it once the corrected version is out, and I'd bring the corrected version to the reading group.","headline":"A genuine first Whitham theory for the BL equation, but the printed modulation equations contain a reciprocal factor error that undermines the Mach reflection predictions as written.","tokens_in":30286,"tokens_out":16431,"would_cite":false,"duration_ms":125917,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["76B15","76B25","35Q53","35Q51"],"pacs":["47.35.Fg","47.35.-i"],"model":"deepseek-v4-flash","headline":"Whitham modulation theory for the isotropic Benney–Luke equation yields analytic formulas for Mach reflection and Mach expansion of solitary waves, and predicts circular rather than parabolic wave wakes, all confirmed by direct numerical…","keywords":["Benney-Luke equation","Whitham modulation theory","Mach reflection","Mach expansion","solitary waves","Rankine-Hugoniot conditions","forced Benney-Luke equation","wave wakes"],"falsifier":"Directly solve the full BL equation for a reverse bent soliton with $a_0=0.21$ over a fine sweep of initial angles near the predicted critical angle from (5.10), extract the Mach stem amplitude and the triple-point velocity, and compare them to the solution of (5.4) and (5.7): a systematic deviation beyond numerical error would falsify the classical Rankine–Hugoniot part of the prediction.","tokens_in":29228,"feed_emoji":"🌊","tokens_out":9869,"duration_ms":87713,"temperature":0.7,"pith_summary":"This paper develops a Whitham modulation theory for the Benney–Luke (BL) equation, an isotropic and bi-directional model of shallow-water waves, and uses it to predict how obliquely colliding solitary waves interact. The theory produces analytical formulas for the critical angle that separates Mach reflection from regular reflection and for the amplitude of the Mach stem that forms in the resonant case; direct numerical solutions of the BL equation reproduce these predictions. Applied to the forced BL equation, the same modulation equations yield a similarity solution describing waves generated by moving topography, whose far-field wavefronts are circular arcs rather than the parabolic fronts of the forced Kadomtsev–Petviashvili equation. The results matter because they show that quantitative soliton-interaction predictions need not be tied to a unidirectional, anisotropic model, and they expose a concrete observable difference between isotropic and anisotropic descriptions of shallow-water wakes.","feed_headline":"Shallow-water wakes are circular, not parabolic, in isotropic theory","feed_subtitle":"Predictions for Mach stems and circular precursor fronts are confirmed by direct numerical simulation.","key_machinery":"The central object is the Whitham modulation system (3.10) for the BL equation, derived from the phase compatibility conditions (3.5) and the solvability condition at next-to-leading order. For soliton dynamics the system simplifies to the $2\\times 2$ hyperbolic system (4.2) whose Riemann invariants (4.4) and characteristic speeds (4.3) produce the rarefaction solutions describing Mach expansion. For discontinuous transitions, the analysis uses the modified Rankine–Hugoniot conditions (5.4) from the full conservation form (5.2) together with classical Rankine–Hugoniot conditions (5.7) from a decomposition of (3.10c) into two one-dimensional conservation laws. The wave-wake description rests on a similarity solution of the modulation equations under the assumption $k=1$ and $a$ independent of $y$, which gives the circular wavefronts in formula (6.2).","core_discovery":"The paper's central claim is that the modulation equations (3.10c,d) for the BL equation, supplemented by the modified Rankine–Hugoniot conditions (5.4) and the classical Rankine–Hugoniot conditions (5.7) obtained by decomposing (3.10c) into two one-dimensional conservation laws, determine the Mach stem and reflected soliton amplitudes and the critical angles for Mach expansion and Mach reflection. For bent solitons with initial slope $q_0<0$, rarefaction-wave solutions give a Mach stem of amplitude fixed by the integral relation (4.8) and a critical angle $\\theta_{cr}$ given by (4.13); for reverse bent solitons $q_0>0$, the algebraic system yields stem and reflected amplitudes $(a_w,q_w)$ and $(a_i,q_i)$ that agree with direct numerical simulations of the BL equation. As the soliton speed approaches unity, the BL predictions converge to those of the KP equation, while for larger speeds the two models diverge. For the forced BL equation, a slowly varying similarity solution of the same modulation equations describes far-field waves as circular arcs, matching numerical simulations for subcritical, critical, and supercritical topography speeds, with the subcritical wake angle given by formula (6.8).","pith_inferences":["The decomposition of (3.10c) into two one-dimensional conservation laws is an ansatz; a rigorous two-dimensional shock theory would be needed to confirm the classical Rankine–Hugoniot predictions beyond the parameter ranges tested numerically.","The circular-versus-parabolic far-field wavefront is an observable signature that could be tested in a shallow-water tank: a subcritical towed obstacle should show circular precursor crests, whereas the forced KP picture predicts parabolic crests.","Because internal-wave versions of the BL equation share the same structure with modified coefficients, the authors' modulation framework could transfer to stratified flows, potentially yielding the same circular-wavefront prediction for internal wave wakes.","The wake-interaction Mach stem threshold (amplitude, speed, separation) seems plausibly mappable onto the single-soliton critical-angle formula, which would turn the empirical parameter scans in section 6.3 into a quantitative prediction."],"forward_implications":["Given an incident soliton amplitude and speed, formulas (4.13), (4.20), (5.9), and (5.10) tell whether a Mach stem will form and how tall it will be, so laboratory or field observations of solitary-wave collisions can be checked against these quantitative thresholds.","At speeds close to the shallow-water wave speed, the BL and KP equations agree on critical angles and stem speeds, meaning the simpler unidirectional KP model remains reliable near resonance, while away from it the isotropic BL description is the better guide.","For moving topography at subcritical speeds, the far-field wavefront is a circle of radius roughly $t$ centered on the starting point; at critical speeds the first precursor front stays approximately circular but is pushed ahead of that circle by the moving topography.","Supercritical topography produces no upstream precursor and confines the wake to an angle $\\arcsin(1/c_b)$, a qualitative change from the subcritical regime where the wake angle spans 0 to 90 degrees as $c_b$ approaches 1.","Two co-moving topographies can generate a Mach stem where their wakes meet, with the stem appearing only for sufficiently large topography amplitude, low enough speed, and separation below about 1000; these conditions remain partly empirical in the paper."],"supporting_citations":[{"why":"Establishes the resonant-soliton mechanism for oblique collisions of shallow-water solitary waves that the paper derives analytically for the BL equation.","marker":"Miles (1977a,b)"},{"why":"Introduces Whitham modulation theory for bent solitons in the KP equation and the Mach expansion analysis that this paper adapts to the BL equation.","marker":"Ryskamp et al. (2021)"},{"why":"Derives the modified Rankine–Hugoniot conditions for KP Mach reflection that the paper extends to the non-integrable BL equation.","marker":"Ryskamp et al. (2022)"},{"why":"Introduced the BL equation, the model studied throughout the paper.","marker":"Benney & Luke (1964)"},{"why":"Shows the BL equation reduces formally to the KP equation, justifying the comparison of the two models.","marker":"Ablowitz et al. (2006)"},{"why":"Provides the Whitham modulation theory for KP and the consistency condition used in deriving the BL modulation equations (3.10).","marker":"Ablowitz et al. (2017)"},{"why":"Analyses upstream-advancing waves from moving disturbances using the forced KP equation and its modulation equations, the starting point for the wave-wake analysis.","marker":"Lee & Grimshaw (1990)"},{"why":"Derived the forced Benney–Luke equation over topography that governs the moving-topography wave problem.","marker":"Milewski (1998)"}],"fun_headline_variants":["Circular wave wakes predicted by shallow-water theory","Isotropic theory shows circular wave wakes","Mach stems and circular wakes from Benney-Luke theory","Circular precursor waves: a new shallow-water prediction","Benney-Luke model flips wake shape from parabola to circle"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The Mach-reflection predictions stand on the unproved claim that equation (3.10c) can be split into two independent one-dimensional conservation laws by alternately neglecting the $x$- and $y$-derivatives; if that split is invalid, the algebraic system for the Mach stem and reflected soliton amplitudes becomes underdetermined.","fun_headline_variants_meta":{"raw":{"variants":["Circular wave wakes predicted by shallow-water theory","Isotropic theory shows circular wave wakes","Mach stems and circular wakes from Benney-Luke theory","Circular precursor waves: a new shallow-water prediction","Benney-Luke model flips wake shape from parabola to circle"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000705,"raw_usage":{"total_tokens":3229,"prompt_tokens":1045,"completion_tokens":2184,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":661,"completion_tokens_details":{"reasoning_tokens":2108}},"tokens_in":661,"tokens_out":2184,"duration_ms":16858,"temperature":1.0,"reasoning_tokens":2108,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T20:57:40.640512+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Directly solve the full BL equation for a reverse bent soliton with $a_0=0.21$ over a fine sweep of initial angles near the predicted critical angle from (5.10), extract the Mach stem amplitude and the triple-point velocity, and compare them to the solution of (5.4) and (5.7): a systematic deviation beyond numerical error would falsify the classical Rankine–Hugoniot part of the prediction.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduced the BL equation, the model studied throughout the paper."},{"cited_title":", Maiden, M","cited_arxiv_id":null,"evidence_quote":"Introduces Whitham modulation theory for bent solitons in the KP equation and the Mach expansion analysis that this paper adapts to the BL equation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Derives the modified Rankine–Hugoniot conditions for KP Mach reflection that the paper extends to the non-integrable BL equation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows the BL equation reduces formally to the KP equation, justifying the comparison of the two models."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Whitham modulation theory for KP and the consistency condition used in deriving the BL modulation equations (3.10)."},{"cited_title":"& Grimshaw, R","cited_arxiv_id":null,"evidence_quote":"Analyses upstream-advancing waves from moving disturbances using the forced KP equation and its modulation equations, the starting point for the wave-wake analysis."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Derived the forced Benney–Luke equation over topography that governs the moving-topography wave problem."}],"review_version":1}