{"id":"e10b354f-be21-492b-8259-8c9abac44c23","arxiv_id":"2508.01890","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A fractional Navier-Stokes model with a Laplacian of order 1/3 is introduced and tested numerically, but remains preliminary and unvalidated.","lead":"This paper proposes a fractional version of the Navier-Stokes equations to model turbulent flows with memory effects. It tests this idea on simplified model equations and a 3D simulation, but leaves calibration and boundary conditions as open problems.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claim that fractional Laplacian order s=1/3 reproduces inertial-range scaling is asserted, not established: the accessible validation uses 1D/Burgers problems, and the 3D experiment is reported only for TKE decay, not for the energy spectrum.","rationale":"The reader's weakest-assumption analysis is the same one I would emphasize: the fractional order 1/3 is load-bearing but not derived or independently validated. The abstract's own limitation list confirms that calibration is open, which strengthens rather than weakens the concern. I do not treat disagreement with existing turbulence-modeling consensus as a problem in itself; the issue is internal to the evidence chain. The paper is honestly labeled preliminary, and the 3D pseudo-spectral experiment is a real step, but a TKE decay curve cannot certify an inertial-range scaling law. Since the full text is corrupted, an equation-level check is not possible; based on the accessible abstract, there is no direct error that would justify rejection, but the central claim is unverified. The reader's UNVERDICTED verdict is therefore appropriate, and my concern does not move it.","tokens_in":10980,"tokens_out":5255,"duration_ms":66540,"concrete_test":"Run the 3D pseudo-spectral fNSE solver at two resolutions (e.g., 128^3 and 256^3) with s=1/3 and matched Reynolds number, and compute the angle-averaged energy spectrum E(k) and compensated spectrum E(k)k^{5/3} from forced or decaying snapshots. Repeat with s=0.2 and s=0.5. If an inertial-range k^-5/3 plateau appears only near s=1/3 and is robust to resolution, the claim gains support; if no plateau appears or the slope is insensitive to s, the choice s=1/3 is not validated. Also report the TKE decay exponent and compare it with a standard NSE DNS at the same resolution and Reynolds number.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that s=1/3 be tied to the 3D inertial range, either by a derivation from the NSE dynamics or by direct spectral validation. The abstract provides neither. It reports validation on 1D advection-diffusion and Burgers equations; Burgers is a poor surrogate for the NSE inertial range because its shock-dominated spectrum is ~k^-2, not k^-5/3. The 3D pseudo-spectral run is described only for turbulent kinetic energy decay, a global quantity that is largely insensitive to the spectral slope. The manuscript itself lists 'Calibrating fractional parameters' as an unresolved challenge, so the value 1/3 is not independently justified. The time-fractional-order parameter is similarly uncalibrated. This does not make the claim false, but it means the key identification between s=1/3 and inertial-range scaling is currently an unsupported assumption rather than a demonstrated result.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a fractional generalization of the Navier-Stokes equations (fNSE) in which the stress-strain relation is replaced by a nonlocal fractional Laplacian and a Caputo time-fractional derivative. The central assertion is that a fractional Laplacian of order 1/3 together with a time-fractional derivative captures non-Markovian energy transfer and the inertial-range scaling of turbulence. Validation is attempted through numerical solutions of 1D advection-diffusion and Burgers equations, a 1D heat equation with Caputo derivative, and a 3D pseudo-spectral simulation of incompressible NSE used to examine turbulent kinetic energy decay. The paper candidly lists open problems, including boundary-condition treatment, hybridization with LES/RANS, and calibration of fractional parameters.","tokens_in":99,"tokens_out":2741,"duration_ms":48874,"significance":"If the fNSE with α=1/3 and a time-fractional derivative were shown to reproduce the 3D inertial-range spectrum and non-Markovian cascade dynamics, it would provide a novel and potentially useful turbulence-closure direction. The paper is honest about its preliminary nature and explicitly flags calibration and boundary-condition issues as unresolved, which is a strength: the limitations are not hidden. However, the current evidence is far from establishing the central claim. The validation uses 1D surrogates, the 3D result is reported only through a global decay curve, and the fractional-order parameters are not derived or calibrated. The contribution is therefore best regarded as a plausible ansatz needing substantial further work rather than a demonstrated framework.","major_comments":[{"comment":"The claim that 'the fractional Laplacian of order 1/3 and time fractional derivative capture non Markovian energy transfer' is asserted without a derivation or an independent calibration. The manuscript itself lists 'Calibrating fractional parameters' as an unresolved challenge, which means α=1/3 and the time-fractional order are free parameters rather than predictions of the theory. To make the claim load-bearing, the authors must either derive α=1/3 from the dynamics of the NSE or show by systematic calibration against DNS or experiments that this value is selected, while other values are not.","section":"Abstract"},{"comment":"The numerical demonstrations on the 1D advection-diffusion equation and the Burgers equation do not test the 3D inertial-range claim. Burgers turbulence is shock-dominated with a k^{-2} energy spectrum, not the k^{-5/3} Kolmogorov spectrum; matching Burgers behavior therefore provides essentially no evidence that the fractional Laplacian of order 1/3 reproduces the 3D inertial range. The spectral behavior of the fNSE must be checked directly in 3D or, at minimum, on a model that shares the 3D scaling (e.g., a shell model), and the authors should state explicitly what the 1D results can and cannot establish.","section":"1D validation (Burgers/advection-diffusion)"},{"comment":"The 3D pseudo-spectral run is described only in terms of turbulent kinetic energy decay, which is a global quantity that is largely insensitive to the spectral slope. No energy spectrum E(k) from the 3D run is shown, and no quantitative comparison against DNS, theory, or the standard k^{-5/3} scaling is reported. Consequently the central identification between α=1/3 and inertial-range scaling currently rests on an unsupported assumption. The authors should present the 3D energy spectrum and, where possible, compare it with DNS or a well-established closure at matched Reynolds number.","section":"3D pseudo-spectral simulation (TKE decay)"},{"comment":"If the value α=1/3 was selected because it makes the fractional Laplacian's Fourier symbol k^{2α} produce a spectral slope consistent with the inertial range, then later demonstrations that the fNSE exhibits that scaling would partly rediscover the input rather than independently confirm it. The logical status of α=1/3 must be clarified: is it an ansatz to be tested, a derived result, or a fitted value? This distinction is important for interpreting all subsequent numerical experiments and should be stated explicitly in the derivation section.","section":"Choice of α=1/3 and potential circularity"}],"minor_comments":[{"comment":"The version of the manuscript provided to me is severely corrupted by a character-encoding problem; most of the body text appears as mojibake, making it impossible to verify equations, section numbering, or the detailed presentation. The authors should ensure that a correctly encoded, readable PDF is deposited.","section":"Entire manuscript text"},{"comment":"The 3D pseudo-spectral run lacks essential details that a reader needs to assess reliability: mesh resolution, Reynolds number, initial conditions, numerical method, and whether the TKE decay is compared with any reference solution. Please add these parameters and, ideally, error bars or convergence checks.","section":"Numerical experiments"},{"comment":"The title and abstract use categorical language ('capture non Markovian energy transfer') for results that the paper itself describes as preliminary and parameter-uncalibrated. Qualifying these statements as conjectures or working hypotheses would better match the evidence presented.","section":"Abstract and title"},{"comment":"The reference list appears garbled and incomplete in the provided text; please verify that all citations are readable and correctly formatted.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is an early-stage exploration and the scientific claims are not yet supported by the evidence shown. The main scientific gap is the unsupported choice of α=1/3 coupled with the absence of a 3D spectral validation; these are fixable within the manuscript's scope if the authors add a derivation or a direct spectral test and clarify the status of the fractional parameters. The severe text corruption in the supplied version also needs to be fixed before the paper can be properly evaluated."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a speculative proposal, not a demonstrated result. The authors propose a fractional Navier-Stokes model with a fractional Laplacian of order 1/3 plus a Caputo time-fractional term, and claim this combination captures inertial-range scaling and non-Markovian energy transfer. That specific combination is new as far as I can tell, and the abstract is candid about what is unresolved.\n\nWhat's good: they actually do numerics - 1D advection-diffusion and Burgers problems, and a 3D pseudo-spectral run for TKE decay. They also explicitly list calibration, boundary conditions, and LES/RANS hybridization as open challenges. That kind of honesty is underrated.\n\nThe soft spots are significant. The central claim - that s=1/3 reproduces the inertial range - is asserted, not derived or validated. The 1D/Burgers tests are poor surrogates for the 3D inertial range; Burgers gives a k^-2 spectrum, not k^-5/3. The 3D run is reported only for TKE decay, a global quantity that won't discriminate between spectral slopes. There's no direct energy spectrum comparison, no error bars, no DNS or experimental reference. The fractional orders are uncalibrated, and the paper says calibration is future work. So the key identification between 1/3 and the inertial range is currently an assumption, not a finding. I don't see a derivation from NSE dynamics that would single out 1/3; it looks like it's chosen to match the desired scaling, which makes the later demonstrations somewhat circular, though not in a parameter-fitting sense.\n\nThe reader's take is largely right. I'd add that the paper's own abstract lists the missing pieces, so nobody is being deceived. It's a research program sketch, not a completed work.\n\nWho's this for? Someone working in fractional models of turbulence might find it a useful pointer, but only as a starting point. I wouldn't cite it as evidence for anything yet.\n\nRecommendation: A serious editor should probably not send this to peer review in its current form - the central claim lacks support. But it's not a cranky paper; if the authors come back with an energy-spectrum comparison and a derivation of the fractional order, it could be worth a look. I'd treat it as a preprint to watch, not a paper to referee.","headline":"A speculative fractional NSE proposal whose central claim about inertial-range scaling is asserted, not yet demonstrated; the paper is honest but preliminary.","tokens_in":11703,"tokens_out":2841,"would_cite":false,"duration_ms":31057,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["47.27.Ak"],"model":"deepseek-v4-flash","headline":"A fractional generalization of Navier-Stokes, with an order-1/3 fractional Laplacian and a Caputo time derivative, is claimed to capture non-Markovian energy transfer and inertial-range scaling in turbulence.","keywords":["fractional Navier-Stokes equation","turbulence modeling","fractional Laplacian","Caputo derivative","non-Markovian energy transfer","inertial range scaling","turbulent kinetic energy decay","pseudo-spectral method"],"falsifier":"A direct test would compare the fNSE's energy spectrum against the Kolmogorov $k^{{-5/3}}$ inertial-range law in a resolved 3D simulation, or check whether the turbulent kinetic energy decay exponent matches direct numerical simulation; if the order-1/3 choice produces a measurably different spectral slope, the central claim is falsified.","tokens_in":10792,"feed_emoji":"🌀","tokens_out":7126,"duration_ms":69574,"temperature":0.7,"pith_summary":"This paper proposes a fractional generalization of the Navier-Stokes equations in which the stress-strain relation is replaced by a fractional constitutive law encoding nonlocal spatial interactions and memory. The authors' central numerical claim is that a fractional Laplacian of order 1/3 together with a time-fractional derivative captures non-Markovian energy transfer and inertial-range scaling in turbulence. They support the formulation with numerical experiments on the 1D advection-diffusion equation, the Burgers equation, the transient heat equation with a Caputo derivative, and a pseudo-spectral 3D incompressible Navier-Stokes simulation of turbulent kinetic energy decay. The paper is explicitly preliminary and lists boundary conditions, LES/RANS hybridization, and parameter calibration as unresolved challenges.","feed_headline":"Fractional Navier-Stokes equation captures turbulence memory","feed_subtitle":"Order-1/3 fractional Laplacian plus a time-memory term reproduces inertial-range energy transfer in numerical tests.","key_machinery":"The load-bearing object is the fractional Navier-Stokes equation (fNSE), obtained by replacing the classical stress-strain relation with a fractional-order constitutive relation. Its two defining operators are the fractional Laplacian of order 1/3, which supplies nonlocal spatial interactions, and the Caputo time-fractional derivative, a time derivative that remembers past states through a power-law kernel and thereby supplies memory (non-Markovian) effects. The argument is carried numerically: a pseudo-spectral method integrates the incompressible fNSE in a 3D periodic box to demonstrate turbulent kinetic energy decay, while the 1D advection-diffusion, Burgers, and heat equations are used to test the spatial and temporal fractional terms separately.","core_discovery":"The central claim is that the Navier-Stokes equations, when generalized through a fractional constitutive relation, yield a fractional Navier-Stokes equation (fNSE) that behaves as a turbulence model with memory. With the fractional Laplacian set to order 1/3 and a Caputo-type time-fractional derivative included, the model is stated to capture non-Markovian energy transfer and the scaling of the inertial range. The numerical evidence is twofold: simplified one-dimensional problems (advection-diffusion, Burgers, transient heat with Caputo derivative) validate the spatial and temporal fractional terms individually, and a pseudo-spectral solution of the incompressible fNSE in a 3D periodic domain demonstrates turbulent kinetic energy decay behavior. The authors present this as a preliminary analysis rather than a closed theory, and they explicitly identify boundary-condition enforcement, hybridizing with LES/RANS, bridging to Lagrangian averaged models, and calibrating the fractional parameters as open problems.","pith_inferences":["Beyond the paper: if the order 1/3 is tied to the -5/3 power law by dimensional analysis, one could derive the fractional order from Kolmogorov scaling rather than treat it as a fitted parameter; the paper does not present such a derivation.","Beyond the paper: the memory term may act as an implicit subgrid-scale model, so the fNSE could be compared directly with large-eddy simulation in the same numerical setup instead of only being hybridized with it in future work.","Beyond the paper: a decisive test would be to run the fNSE in forced, statistically stationary turbulence and check whether the energy flux stays constant across scales, something the paper's decay-only setup does not examine."],"forward_implications":["If the fNSE reproduces inertial-range scaling and non-Markovian transfer, it can serve as a single-equation alternative to eddy-viscosity closures that model unresolved stress with local damping.","The same fractional formulation is portable to other transport equations, since the authors validate it on 1D advection-diffusion, Burgers, and heat equations before applying it to 3D turbulence.","The 3D pseudo-spectral run shows that the fNSE is numerically integrable in a periodic box, giving a concrete baseline for future fractional turbulence studies.","The paper's own list of open problems—boundary conditions, LES/RANS hybridization, coupling to Navier-Stokes-alpha, and parameter calibration—defines the work needed before this approach can predict engineering flows."],"supporting_citations":[],"fun_headline_variants":["Fractional Navier-Stokes: turbulence with memory","Non-Markovian Navier-Stokes models turbulence memory","Fractional Laplacian and time-derivative: memory in turbulence","Turbulence energy transfer with fractional Navier-Stokes memory","Fractional NSE: a memory-aware turbulence framework"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The claim that order 1/3 is the right fractional Laplacian order to reproduce inertial-range scaling is assumed, not derived, and the paper lists calibrating fractional parameters as future work; if a different order fits the data as well or better, the model's predictive power is undercut.","fun_headline_variants_meta":{"raw":{"variants":["Fractional Navier-Stokes: turbulence with memory","Non-Markovian Navier-Stokes models turbulence memory","Fractional Laplacian and time-derivative: memory in turbulence","Turbulence energy transfer with fractional Navier-Stokes memory","Fractional NSE: a memory-aware turbulence framework"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000507,"raw_usage":{"total_tokens":2487,"prompt_tokens":977,"completion_tokens":1510,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":593,"completion_tokens_details":{"reasoning_tokens":1427}},"tokens_in":593,"tokens_out":1510,"duration_ms":16650,"temperature":1.0,"reasoning_tokens":1427,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T05:18:28.617419+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct test would compare the fNSE's energy spectrum against the Kolmogorov $k^{{-5/3}}$ inertial-range law in a resolved 3D simulation, or check whether the turbulent kinetic energy decay exponent matches direct numerical simulation; if the order-1/3 choice produces a measurably different spectral slope, the central claim is falsified.","supporting_citations":[],"review_version":1}