{"id":"2bcfdf03-d03e-4040-b2aa-1dd0791bbdf1","arxiv_id":"2607.10176","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Wall-modeled WPTS recovers flat-plate transition skin-friction and log-layer velocity on a ~0.2 M-cell mesh better than GKS by adding tunable particle non-equilibrium transport.","lead":"The authors write down the full kinetic equations for wave-particle turbulence simulation (WPTS) and couple it to a standard wall model. On a coarse mesh the method predicts flat-plate transition skin friction and mean velocity closer to DNS than a pure gas-kinetic scheme, by letting stochastic particles carry unresolved non-equilibrium transport.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.5","headline":"The claimed superiority of WM-WPTS rests on a single hand-tuned free coefficient that is not independently constrained.","rationale":"The Reader correctly isolates the hand-tuned mixing-length coefficient as the weakest link supporting the central numerical claim. The derivation of the wave-particle equations is clean and the wall-model coupling is a legitimate first step, yet the only quantitative evidence that WPTS outperforms GKS on under-resolved transitional flow is obtained after that single free parameter is adjusted to match the DNS transition location. Because the paper presents no independent constraint on C_m^{2} and only one flow configuration, the superiority claim remains conditional on successful tuning. My concrete test simply asks whether the same coefficient works without re-tuning; a negative answer would confirm that the improvement is not yet predictive. No stronger internal inconsistency is present, so the Reader’s CONDITIONAL verdict is unchanged.","tokens_in":10068,"tokens_out":509,"duration_ms":6766,"concrete_test":"Re-run the M2 calculation with C_m^{2} fixed at the value obtained from a separate, non-transitional calibration (e.g., a fully-developed channel or the same plate at a different Re) and report the resulting c_f curve and U(y) at x=9.9. If the transition location again requires re-tuning of C_m^{2} by more than ~20 %, the claim of robust multi-scale superiority over GKS is not supported.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The strongest claim is that WM-WPTS on coarse mesh M2 produces a skin-friction curve and mean-velocity profile that match DNS and are markedly better than WM-GKS under identical settings. That improvement appears only after the free coefficient C_m^{2} in the mixing-length closure (Eq. 21) is set to 0.032; when C_m^{2} = 0 the method reverts exactly to the failing GKS baseline (Fig. 2 and surrounding text). No a-priori estimate, grid-convergence study, or independent calibration of C_m^{2} is supplied, nor is any other transitional or wall-bounded case shown. Consequently the numerical success is not yet a demonstration that the wave-particle decomposition itself supplies predictive multi-scale modelling; it remains a demonstration that a suitably tuned particle-generation rate can restore transition location on one mesh for one flow.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The paper derives the complete kinetic model equations of wave-particle turbulence simulation (WPTS) from an additive decomposition of the distribution function into wave (grid-resolved) and particle (stochastic, subgrid) components, interprets the generation/annihilation and forcing terms, and shows that the combined equation recovers a BGK-type relaxation that conserves mass, momentum and energy. It then couples WPTS to an equilibrium wall model and applies the resulting WM-WPTS method to a forced flat-plate transition at Ma = 0.7, Re = 5e4. On a coarse mesh (~0.23 M cells) the method, with a mixing-length closure for the turbulent time scale, produces a skin-friction curve and a mean-velocity profile that match DNS data more closely than wall-modeled GKS under identical settings, thereby claiming multi-scale predictive capability for transitional wall-bounded flows.","tokens_in":10402,"tokens_out":1144,"duration_ms":20698,"significance":"If the wave-particle decomposition and its non-equilibrium particle transport genuinely supply a grid-adaptive multi-scale model that automatically reduces to a kinetic NS solver in laminar regions, the framework would offer a useful alternative to conventional LES/RANS hybrids for transitional flows. The explicit derivation of the model equations (Section 2) and the first demonstration of wall-model coupling are concrete contributions that strengthen the theoretical foundation of earlier WPTS work. The numerical improvement over GKS on a deliberately under-resolved mesh is of practical interest, provided the free coefficients can be constrained independently of the target data.","major_comments":[{"comment":"Section 4.2 and Eq. (21): the turbulent time scale that controls particle generation, annihilation and free-transport weight is closed by a mixing-length formula whose free coefficient C_m^{2} is hand-tuned to 0.032 so that the skin-friction rise coincides with the DNS target (Fig. 2). When C_m^{2} = 0 the method reverts exactly to the failing WM-GKS baseline. No a-priori estimate, grid-convergence study, independent calibration on a different flow, or sensitivity analysis is supplied. Consequently the claimed superiority of WM-WPTS rests on a single adjustable parameter fitted to the very quantity being predicted, which undercuts the assertion of predictive multi-scale modelling.","section":"§4.2, Eq. (21), Fig. 2"},{"comment":"The entire numerical validation consists of one forced-transition case at a single Mach and Reynolds number, compared on two meshes against one DNS reference. No additional transitional, fully turbulent, or separated wall-bounded flow is shown, nor is any demonstration that the same C_m^{2} (or a fixed procedure for choosing it) works elsewhere. This single-case, single-parameter success is insufficient to support the abstract’s claim of “considerable promise or transitional flow simulations.”","section":"§4, Abstract"},{"comment":"Particle sampling (Eqs. 12–13) assumes the production of turbulent kinetic energy to be C_0 e^{-Δt/τ_n} \rho e with C_0 = 1 and e = ½U^{2}. This choice is ad-hoc, dimensionally inconsistent with a true TKE production term, and is never varied or justified against measured production budgets. Because the particle fraction (and therefore the non-equilibrium transport) depends directly on this assumption, its influence on the reported skin-friction and mean-velocity results should be quantified.","section":"§3, Eqs. (12)–(13)"}],"minor_comments":[{"comment":"Figure captions and text repeatedly use “steamwise” for “streamwise” (Figs. 1b, 2b and surrounding prose).","section":"Figs. 1–2"},{"comment":"In §4.2 the label “WM-WPTSC 2 m=0” is a typesetting artifact; it should read “WM-WPTS with C_m^{2} = 0 (i.e., WM-GKS)”.","section":"§4.2"},{"comment":"The high-density modification α(ρ) and the damping length y_ref = 0.15 appear without reference or sensitivity test; a brief justification or citation would help reproducibility.","section":"§4.2, Eq. (21)"},{"comment":"Table 1 lists wall units for the DNS mesh R1 but does not state the friction velocity used to non-dimensionalize the coarser meshes M1/M2; consistency of the reported y+ values should be clarified.","section":"Table 1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a natural incremental step from the authors’ recent WPTS papers; the principal novelty is the explicit equation derivation and the wall-model coupling. The tuning issue is serious enough that the numerical claims should not be accepted without additional independent calibration or multi-case evidence, but the theoretical sections are sound and the paper is within scope for a fluids journal."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The new material here is the explicit wave-particle kinetic equations (3–6) and the first wall-modelled WPTS run on a forced flat-plate transition. Both are useful. The equations are written carefully: they reduce exactly to GKS when the particle weight vanishes, the force terms cancel so that the combined distribution still satisfies the BGK conservation properties, and the free-transport flux is constructed so that surviving particles are not double-counted. That part is solid and was missing from the earlier WPTS papers.\n\nOn the coarse mesh (~0.23 M cells) the method with C_m^{2} = 0.032 recovers a skin-friction rise and a log-layer mean profile that sit close to the DNS reference and are clearly better than pure wall-modelled GKS under identical settings. The wall-model coupling itself is standard and does what it is supposed to do. So the numerical claim in the abstract is true for this one case.\n\nThe soft spot is exactly the one the stress-test flags. Particle generation and lifetime are controlled by a mixing-length eddy viscosity whose free coefficient C_m^{2} is dialled until the transition location matches DNS. Set it to zero and you recover the delayed GKS result. No a-priori estimate, no second transitional case, no comparison with any established wall-modelled LES, and no code or data release. The free-parameter list is longer than one would like (C_0, y_ref, density limiter, wall-model activation thresholds). That does not make the equations wrong; it does mean the present evidence is still a demonstration that a suitably tuned particle rate can restore transition, not yet a demonstration of predictive multi-scale modelling.\n\nWho should read it: people already working with kinetic schemes or looking for hybrid wave-particle ideas for transitional wall flows. It is not yet a paper that reorganises the broader LES/RANS conversation. I would send it to referees; the derivation is clean enough and the numerical improvement real enough to deserve a proper look, provided the authors are pressed on calibration and additional cases. I would not cite it myself until the coefficient is constrained independently or more flows appear.","headline":"Clean first derivation of the WPTS equations plus a wall-modelled flat-plate demo that works only after the particle-generation coefficient is hand-tuned to the DNS transition location.","tokens_in":10967,"tokens_out":548,"would_cite":false,"duration_ms":6932,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"Wave-particle turbulence simulation with a wall model predicts flat-plate transition on coarse grids by letting stochastic particles carry unresolved non-equilibrium transport.","keywords":["wave-particle turbulence simulation","wave-particle decomposition","flat-plate transition","wall model","gas-kinetic scheme","non-equilibrium transport","subgrid-scale modeling","coarse-grid simulation"],"falsifier":"Rerun the identical coarse-mesh flat-plate case with the same wall model but a different free coefficient (or none) and check whether the skin-friction rise still locks onto the direct-simulation location and whether the mean-velocity profile still recovers a clean log layer.","tokens_in":10930,"feed_emoji":"🌊","tokens_out":642,"duration_ms":11290,"temperature":0.7,"pith_summary":"This paper derives the full set of model equations for wave-particle turbulence simulation (WPTS) by splitting the fluid distribution into a wave part that evolves the large, grid-resolved motions and a particle part that evolves the unresolved motions through free streaming, forcing, and annihilation. It then couples WPTS to a simple equilibrium wall model so that near-wall resolution can be relaxed. On a coarse mesh of roughly two hundred thousand cells the method recovers the skin-friction rise and the mean velocity log-layer of a forced flat-plate transition that previously required a multi-million-cell direct simulation, while a pure gas-kinetic scheme on the same mesh delays transition badly. The practical payoff is a multi-scale kinetic framework that automatically injects more particles only where the grid can no longer resolve the turbulence, offering a route to transitional and wall-bounded flows without extreme mesh refinement.","feed_headline":"Particles recover flat-plate transition on a coarse mesh","feed_subtitle":"Wall-modeled wave-particle scheme matches DNS skin friction where pure gas-kinetic fails","key_machinery":"Wave-particle decomposition of the distribution function: the wave component supplies the multi-scale interface flux for resolved motion while stochastic particles, generated and annihilated according to a local turbulent time scale, supply the non-equilibrium free-transport contribution that models sub-grid scales.","core_discovery":"The complete kinetic model equations of WPTS follow directly from a wave-particle decomposition of the distribution function; once these equations are closed by a mixing-length turbulent time scale and coupled to an equilibrium wall model, the resulting scheme on a coarse mesh yields skin-friction and mean-velocity profiles for flat-plate transition that match direct-numerical-simulation data far better than the pure gas-kinetic scheme under identical grid and wall-model settings.","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["Wave-particle equations yield coarse-mesh flat-plate transition","Wall-modeled WPTS recovers DNS skin friction where GKS fails","Particle component models SGS non-equilibrium on coarse grids","WPTS wall model matches DNS velocity for flat-plate transition","Wave-particle split enables accurate coarse-grid transition"],"cache_read_input_tokens":128,"weakest_assumption_plain":"The free coefficient that multiplies the mixing-length formula for the turbulent time scale must be hand-tuned so that particles appear at the right rate; without that tuning the method reverts to ordinary gas-kinetic behaviour and the transition is delayed.","fun_headline_variants_meta":{"raw":{"variants":["Wave-particle equations yield coarse-mesh flat-plate transition","Wall-modeled WPTS recovers DNS skin friction where GKS fails","Particle component models SGS non-equilibrium on coarse grids","WPTS wall model matches DNS velocity for flat-plate transition","Wave-particle split enables accurate coarse-grid transition"]},"model":"grok-4.5","effort":"low","cost_usd":0.003714,"raw_usage":{"total_tokens":1189,"prompt_tokens":759,"num_sources_used":0,"completion_tokens":87,"cost_in_usd_ticks":37140000,"prompt_tokens_details":{"text_tokens":759,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":343,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":759,"tokens_out":87,"duration_ms":4396,"temperature":1.0,"reasoning_tokens":343,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-14T13:43:02.652619+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Rerun the identical coarse-mesh flat-plate case with the same wall model but a different free coefficient (or none) and check whether the skin-friction rise still locks onto the direct-simulation location and whether the mean-velocity profile still recovers a clean log layer.","supporting_citations":[],"review_version":1}