{"id":"d677504c-d544-4e0c-a49e-1e645f312347","arxiv_id":"1908.09626","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":4,"one_line_summary":"A new boundary-condition-flexible spectral formulation for pipe-flow linear stability matches Meseguer-Trefethen accuracy and reproduces the inviscid lift-up optimal disturbance patterns.","lead":"The paper builds a new spectral formulation for computing the stability of flow through a pipe, replacing the radial coordinate with its square so the unknown functions stay smooth at the pipe center. The method matches the best existing calculations for the slowest-decaying disturbances and keeps the same equations when boundary conditions change.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Undisclosed eigenvalue filter in the Appendix B code discards all modes with decay rate exceeding α (or 1), so the claimed 'spectrum' is a filtered subset and the central accuracy claim may not cover the complete spectrum.","rationale":"The reader's verdict of CONDITIONAL is reasonable, and its three requested items (boundary-condition demonstration, filter documentation, and condition-number sentence correction) are real. Among these, the undisclosed eigenvalue filter is the most load-bearing because it directly affects the central claim that the system 'produces a spectrum that is as accurate as Meseguer & Trefethen'. The reader listed documentation of the mode-selection threshold as one of the conditional items, so there is partial agreement. However, the reader's designated weakest assumption was the Priymak-Miyazaki analyticity ansatz, which I do not consider a serious correctness risk: physical pipe-flow perturbations are smooth in Cartesian coordinates, so after factoring r^|n|-1 the remaining fields must be analytic in y=r^2; the regularity conditions therefore select the physical branches and should not delete any physical modes. The condition-number sentence is a typo ('closer' should be 'less close') and the boundary-condition flexibility, while not numerically demonstrated, is structurally consistent because the wall conditions appear only as row replacements and do not enter the derivation of the operators or the centreline regularity conditions. The filter, by contrast, is part of the actual code that produces every reported result; it is arbitrary, undocumented, and may remove physical modes or hide the insufficiency of the regularity conditions. Because the accuracy of the least-decaying modes is still convincingly validated, I do not recommend changing the verdict, but the paper should disclose and justify the filter or remove it and demonstrate that the full spectrum is still accurate.","tokens_in":25608,"tokens_out":15755,"duration_ms":156371,"concrete_test":"Remove the filter in the Appendix B code by commenting out the lines `ok=((imag(e)>-1*al)&(al>=1))|((imag(e)>-1)&(al<1)); e=e(ok); q=q(:,ok);` and recompute for α=1, n=1, Re=3000 and for α=10, n=7, Re=2000, with N varying from 47 to 200. Determine whether eigenvalues with -imag(omega) > α converge with increasing N to fixed, physical values, and compare any converged discarded modes against an independent computation (e.g., Meseguer-Trefethen or a filtered-out Burridge-Drazin run). If converged physical modes are being discarded, the spectrum claim is overstated; if the discarded eigenvalues all diverge or are clearly spurious, then the threshold should be documented as a convergence criterion in the text rather than left as a hidden filter.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim in the abstract is that the system 'produces a spectrum that is as accurate as Meseguer & Trefethen'. The full MATLAB code in Appendix B, however, contains a filter that is never described in the body of the paper: after solving the generalized eigenvalue problem, the line `ok=((imag(e)>-1*al)&(al>=1))|((imag(e)>-1)&(al<1)); e=e(ok); q=q(:,ok)` retains only eigenvalues with -imag(omega) ≤ α for α ≥ 1, and -imag(omega) ≤ 1 for α < 1. Every spectrum displayed in Figs. 1 and 2 is therefore a filtered subset. If the discarded high-decay eigenvalues are physical, the method does not produce a complete spectrum. If they are spurious, then the regularity conditions (22)-(27) do not by themselves remove all non-analytic branches; the filter is performing essential work. In either case, the unqualified claim that the system produces an accurate spectrum is overstated. The validations in Table 1 and the claim that all 41 eigenvalues from [21] match are made on the filtered subset only, so they do not test the discarded region. This is a load-bearing gap because the advertised capability is the spectrum itself, and the code as provided does not deliver it without an arbitrary, undocumented threshold.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript proposes a spectral collocation formulation for the linear stability of Hagen-Poiseuille pipe flow. The authors use the Priymak-Miyazaki ansatz to factor out the centerline power-law behavior of the perturbation velocities, introduce analytic variables φ and Ω in the radial variable y=r², and derive the governing equations (15)-(20) together with the centerline regularity conditions (22)-(27). The claimed advantages are spectral accuracy comparable to Meseguer & Trefethen and the ability to impose different wall boundary conditions directly on the unknowns without reformulating the operators. Validation is provided through agreement with eigenvalues from Meseguer & Trefethen, Schmid & Henningson, and Priymak & Miyazaki in Table 1, and with the exact characteristic roots of the α=n=0 Stokes modes in Table 2. The paper also derives an inviscid algebraic-growth problem in §4 and shows the associated optimal vortex and streak patterns.","tokens_in":25842,"tokens_out":7777,"duration_ms":85387,"significance":"The core technical contribution is plausible, and the numerical validation is strong where it applies. Table 1 shows agreement to 8-14 significant digits against three independent codes, and Table 2 is a clean internal consistency check against exact characteristic relations to about 14 digits. The inviscid optimal-growth derivation in §4 is parameter-free and reproduces the known lift-up/streak phenomenology without fitted constants. The inclusion of the complete MATLAB code in Appendix B is also a practical strength. However, the central claim that the system produces a spectrum as accurate as Meseguer & Trefethen is not supported for the full spectrum, because the code silently discards all eigenvalues whose decay rate exceeds α (or 1). The significance of the paper therefore depends on either demonstrating that the discarded modes are spurious or, if they are physical, restricting all claims to the least-decaying part of the spectrum.","major_comments":[{"comment":"Appendix B contains an undocumented filter that removes eigenvalues before the reported spectra are produced. In both the n≠0 block and the n==0 block, after solving the generalized eigenproblem, the code executes `ok=((imag(e)>-1*al)&(al>=1))|((imag(e)>-1)&(al<1)); e=e(ok); q=q(:,ok)`. This discards every mode with decay rate -Im(ω) > α for α≥1, and -Im(ω) > 1 for α<1; no part of §3.1 or §3.2 describes this operation. The filter is not merely removing numerical noise: for the α=n=0 case, the exact characteristic relations (35)-(36) have infinitely many physical roots with -Im(ω) growing without bound, so the code demonstrably deletes physical Stokes modes. Consequently, the abstract's unqualified claim that the system produces a spectrum as accurate as Meseguer & Trefethen is not supported; the precise statement is that it computes a decay-rate-limited subset. The authors should either remove the filter and compare the complete spectra with reference codes, show that the discarded modes are spurious, or explicitly and consistently restate the claim for the least-decaying part of the spectrum.","section":"Appendix B"},{"comment":"The statement that all 41 eigenvalues listed in [21] matched at similar accuracy is not documented. Table 1 lists only the least-decaying eigenvalue for each parameter set, so the reader cannot see which 41 modes were compared, what threshold was used, or whether any of them lay in the region removed by the Appendix B filter. Because this comparison is the main evidence for the headline accuracy claim, the authors should tabulate or plot the complete comparison, including the decay rates of the compared modes.","section":"§3.2"},{"comment":"The advertised advantage over Meseguer & Trefethen, namely that boundary conditions can be modified without changing the formulation, is asserted but never demonstrated. All computed results use the no-slip conditions (21); no example with a different wall condition, such as slip, transpiration, or a time-dependent condition, is given. Since this flexibility is a central contribution of the paper, at least one demonstration with a non-no-slip boundary condition, ideally with a comparison against an independent method for that condition, is needed before the claim can be evaluated.","section":"§1 and §5"}],"minor_comments":[{"comment":"In the display following Eq. (48), the upper limit of the integration appears as y, although y is the integration variable; this should presumably be 1.","section":"§4, Eq. (49)"},{"comment":"The paragraph comparing condition numbers refers to Eqs. (42)-(42); the second equation number should be (43).","section":"§3.2"},{"comment":"There is a typo in the sentence 'such requirement is satisfied by the the analytic Jν'; the duplicated 'the' should be removed.","section":"§2.2"},{"comment":"The option `opts.tol=1e-21` asks for a tolerance far below double-precision machine epsilon; this cannot improve the computed eigenvalues and should either be removed or explained in a comment.","section":"Appendix B"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The two things to know: the paper actually delivers a new spectral formulation for pipe flow stability that matches Meseguer-Trefethen accuracy with a simpler, BC-flexible structure, and it prints the full MATLAB code. The catch is that the code contains a filter that silently drops eigenvalues, so the 'spectrum' shown is a subset, and the text never tells you.\n\nWhat's genuinely good: the 2-tuple analytic-in-y formulation is new, and the regularity argument based on pole multiplicity is a real contribution. The validation is strong — Table 1 agrees with three independent codes at 8–14 digits, and Table 2 nails the alpha=n=0 characteristic roots to ~14 digits. The inviscid optimal-growth section is a clean, parameter-free derivation that reproduces the lift-up/streak picture without viscosity. That is solid work.\n\nThe soft spots are real but not fatal. The filter is the big one: in Appendix B, after solving the generalized eigenproblem, only eigenvalues with -imag(omega) <= alpha (or <=1 for alpha<1) are kept. The body text does not mention this anywhere. As a result, the claim 'the system produces a spectrum' is overstated, and Figures 1 and 2 show only the filtered subset. The validation tables are not affected, but the reader of the paper has no way to know whether the discarded region is physical or spurious. That is a gap that needs to be fixed, either by documenting and justifying the threshold or by showing the discarded modes are spurious.\n\nThe boundary-condition flexibility is asserted but never demonstrated: no example with a non-no-slip wall condition appears, so the main advertised advantage remains unproven. Also, the condition-number sentence in Section 3.2 contradicts its own numbers: cond ~10^18 for the present system vs ~10^21 for Burridge-Drazin means the present system is less, not 'much closer', to singularity. Minor wording issue, but confusing.\n\nWho is this for: anyone working on linear or nonmodal stability of pipe flow, and anyone building spectral methods with arbitrary boundary conditions. It deserves a serious referee — the core method is sound and the code is transparent enough that the issues are discoverable. I'd suggest asking the authors to document the mode-selection filter, demonstrate the BC flexibility with one worked example, and fix the condition-number sentence.","headline":"A genuinely accurate and new spectral formulation for pipe flow, but the code's undocumented eigenvalue filter overstates the 'spectrum' claim.","tokens_in":26441,"tokens_out":3774,"would_cite":true,"duration_ms":33784,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65M70","76E05","65L15"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper derives a two-variable spectral system for pipe-flow stability that fixes the centreline regularity conditions once, so wall boundary conditions can be swapped without reformulating the operators, while matching reference…","keywords":["pipe flow stability","linear stability","spectral collocation","centreline regularity","analytic centreline ansatz","algebraic growth","optimal disturbances","circular pipe flow"],"falsifier":"Solve the same linear stability problem with an independent discretisation that does not impose the analyticity-in-y ansatz — for example, a high-order finite-difference scheme in the original radial coordinate r on a very fine grid, allowing non-analytic centreline behaviour — and compare the full spectra. If any eigenvalue appears there that the present regularity conditions exclude, especially at high azimuthal wavenumber |n| or high axial wavenumber α, the ansatz has discarded a physical mode and the central claim fails; if none appears, the regularity conditions are confirmed to cut exactly the spurious branches.","tokens_in":25242,"feed_emoji":"🌊","tokens_out":11386,"duration_ms":103952,"temperature":0.7,"pith_summary":"The paper sets out a spectral method for the linear stability of flow through a circular pipe that keeps the accuracy of the best existing approach while making one previously rigid part of the calculation flexible: the wall boundary condition. The state is written in two unknowns, a factored radial velocity and a factored radial vorticity, chosen so that both are analytic in the squared radius y=r². The regular singularity at the pipe centre then yields exactly the right number of regularity conditions to complement the wall conditions and make the eigenvalue problem well posed, and those conditions are derived once and for all. For axisymmetric, axially constant disturbances the computed eigenvalues agree with the analytical characteristic relations to about 14 digits. If the claim holds, the same equations can be reused with slip, permeable, compliant, or time-dependent walls without re-deriving the stability operators.","feed_headline":"Pipe-flow spectra at top accuracy, with wall conditions free to change","feed_subtitle":"A two-variable formulation computes spectra to ~14 digits and lets you swap wall conditions without re-deriving operators.","key_machinery":"The load-bearing object is the pair (φ, Ω) of factored radial velocity and radial vorticity, together with the coordinate change y=r². Under the standard power-law centreline ansatz, φ and Ω are analytic in y near the centre, so ordinary Chebyshev polynomials of all orders can represent them. The pipe centre is a regular singularity — a point where the coefficient of the highest derivative vanishes, allowing the equations to admit non-analytic branches — and its pole multiplicity is such that the regularity limits of the equations, plus one differentiated equation, produce exactly three centreline conditions for each sixth-order case. That exact count is what lets the three wall conditions be imposed directly on the unknowns, so different wall laws can be substituted without changing the operators.","core_discovery":"The central claim is that the linearised disturbance equations for pipe flow can be reduced to a sixth-order system in two unknowns — or a fourth-order system for the axisymmetric, axially constant case — with the wall boundary conditions imposed directly on the unknowns rather than built into the basis functions. The unknowns are chosen so that, after factoring out the known power-law behaviour at the centreline, they are analytic functions of y=r². The regular singularity at y=0 has a pole multiplicity that supplies exactly the missing number of centreline conditions: requiring all y-multiplied operator terms to vanish there gives two regularity conditions, and one further condition obtained by differentiating the governing equation ensures fourth-order differentiability of the leading unknown. The resulting spectra match the established high-accuracy benchmark, and in the special case of axisymmetric axially constant modes the eigenvalues match the roots of the analytical characteristic equations to about 14 digits. In the inviscid limit for axially constant modes, the same working variables give an explicit solution with linear-in-time growth of streamwise velocity, reproducing the counter-rotating vortex and streak patterns familiar from plane shear flows.","pith_inferences":["If the centreline analyticity assumption transfers, the same regularity-counting strategy should work for other radially inhomogeneous base flows (annular pipes, swirling pipe flow, heated or stratified pipes), because the pole-counting argument is tied to the structure of the singularity, not to the no-slip law.","The boundary-condition flexibility invites a concrete test: implement the same operators with a slip-length or suction/blowing condition and compare against an independent solver; agreement would confirm that the claimed flexibility is real, and disagreement would localise the failure in the regularity conditions.","The explicit inviscid solution implies a quantitative prediction for optimal streaks — their wall distance should shrink as the azimuthal wavenumber n grows, since v'=r^l φ0 peaks nearer the wall — which could be checked against nonlinear DNS at moderate Reynolds numbers.","A practical consequence the authors do not spell out: users must apply the stretching map for large α or Re because φ and Ω may develop boundary-layer behaviour at the centre; the transition point could be detected by monitoring the condition number of the discretised operator."],"forward_implications":["The same system of equations can be applied to a range of wall conditions (slip, permeable, compliant, time-dependent) by changing only the boundary rows, not the formulation.","Reference-level accuracy is reached with modest discretization — a 93-by-93 system at N=48 in the standard test cases — and the y-stretching map controls round-off for large axial wavenumber and Reynolds number.","In the axisymmetric, axially constant case, the computed eigenvalues match the analytical characteristic relations to about 14 digits, providing a clean validation benchmark.","The inviscid axially constant analysis gives an explicit solution φ(t,y)=φ0(y), Ω(t,y)=Ω0(y)−2in U_y t φ0(y), showing algebraic growth and optimal vortex/streak patterns like those in plane shear flows.","The pole-counting argument guarantees the conditions are neither under- nor over-specified, so the method does not introduce spurious spectral branches through incorrect regularity treatment."],"supporting_citations":[{"why":"supplies the centreline power-law ansatz and the technique of factoring the leading radial behaviour before discretisation","marker":"[19]"},{"why":"defines the reference spectral benchmark whose accuracy the present method matches","marker":"[21]"},{"why":"provides the high-Reynolds benchmark spectra used for visual validation","marker":"[22]"},{"why":"introduces the two-variable (φ, Ω) formulation whose boundary conditions the present method generalises and with which it is compared","marker":"[14]"},{"why":"gives the analytical characteristic relations for axisymmetric axially constant modes used as the double-precision validation benchmark","marker":"[15]"},{"why":"supplies the earlier viscous optimal-growth computation whose near-wall streak patterns the inviscid result reproduces","marker":"[7]"},{"why":"provides the plane-shear-flow algebraic-growth solution whose pipe analogue is derived in Section 4","marker":"[24]"}],"fun_headline_variants":["14-digit pipe spectra, boundary conditions on demand","Two unknowns, six orders, precise pipe stability with flexible BCs","Analytic variables unlock easy boundary tweaks in pipe stability","Pipe flow spectra: double precision, no operator re-derivation","Boundary-flexible pipe stability with analytic unknowns"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation assumes every physical disturbance is analytic in r² at the pipe centre once the known power-law factor is removed; if a genuine mode has fractional-power or logarithmic centreline behaviour, the regularity conditions imposed here would delete it from the spectrum.","fun_headline_variants_meta":{"raw":{"variants":["14-digit pipe spectra, boundary conditions on demand","Two unknowns, six orders, precise pipe stability with flexible BCs","Analytic variables unlock easy boundary tweaks in pipe stability","Pipe flow spectra: double precision, no operator re-derivation","Boundary-flexible pipe stability with analytic unknowns"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000221,"raw_usage":{"total_tokens":1434,"prompt_tokens":912,"completion_tokens":522,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":528,"completion_tokens_details":{"reasoning_tokens":442}},"tokens_in":528,"tokens_out":522,"duration_ms":6620,"temperature":1.0,"reasoning_tokens":442,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:48:28.177385+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the same linear stability problem with an independent discretisation that does not impose the analyticity-in-y ansatz — for example, a high-order finite-difference scheme in the original radial coordinate r on a very fine grid, allowing non-analytic centreline behaviour — and compare the full spectra. If any eigenvalue appears there that the present regularity conditions exclude, especially at high azimuthal wavenumber |n| or high axial wavenumber α, the ansatz has discarded a physical mode and the central claim fails; if none appears, the regularity conditions are confirmed to cut exactly the spurious branches.","supporting_citations":[{"cited_title":"Priymak, T","cited_arxiv_id":null,"evidence_quote":"supplies the centreline power-law ansatz and the technique of factoring the leading radial behaviour before discretisation"},{"cited_title":"Meseguer, L","cited_arxiv_id":null,"evidence_quote":"defines the reference spectral benchmark whose accuracy the present method matches"},{"cited_title":"Meseguer, L","cited_arxiv_id":null,"evidence_quote":"provides the high-Reynolds benchmark spectra used for visual validation"},{"cited_title":"Burridge, P","cited_arxiv_id":null,"evidence_quote":"introduces the two-variable (φ, Ω) formulation whose boundary conditions the present method generalises and with which it is compared"},{"cited_title":"Salwen, C","cited_arxiv_id":null,"evidence_quote":"gives the analytical characteristic relations for axisymmetric axially constant modes used as the double-precision validation benchmark"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the earlier viscous optimal-growth computation whose near-wall streak patterns the inviscid result reproduces"},{"cited_title":"Ellingsen, E","cited_arxiv_id":null,"evidence_quote":"provides the plane-shear-flow algebraic-growth solution whose pipe analogue is derived in Section 4"}],"review_version":1}