{"id":"c36537c9-9769-42df-8f2f-5d838160e0a1","arxiv_id":"1908.01770","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A joint position-velocity Fokker-Planck framework for advective-diffusive transport is introduced, with closure-based velocity-PDF equations whose stationary solutions match the Eulerian velocity distribution in Poiseuille flow.","lead":"This paper derives a Fokker-Planck equation for the joint probability density of particle position and advective velocity in steady heterogeneous flow, then reduces it to a velocity-only evolution equation. The reduction requires closures, and the paper tests an equilibrium closure on simple shear flows, reproducing known results.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The proposed evolution equation is exact but unclosed; its practical use depends on the equilibrium closure (Eq. 14), which is validated only by a self-consistency check, not against the true time-dependent dynamics in heterogeneous media.","rationale":"The paper's strongest contribution is the exact marginal equation (10), which I verified by direct integration of the joint Langevin dynamics; the negative sign on DΔu and the irrelevance of the cross-diffusion term after spatial averaging are consistent. The joint equation (6) has a matrix-orientation error in the cross term (should be G^T, not G), but this does not affect the marginal equation, so it is not the critical issue. The load-bearing fragility is the equilibrium closure (Section 3.2). Eq. (10) is unclosed because the conditional expectations are functionals of the concentration field. Replacing them by their equilibrium values (14) is accurate only at long times or for nearly uniform plumes. The paper's Fig. 1 is a self-consistency check: it shows the closed equation relaxes to the correct stationary PDF, which any closure with the correct equilibrium would do. It does not test the time-dependent path. For the intended application to heterogeneous porous media, the closure has no evidence. Thus the paper supports the exact derivation but not yet the practical closed model. This warrants a conditional acceptance pending direct validation in a heterogeneous flow.","tokens_in":136,"tokens_out":16558,"duration_ms":366910,"concrete_test":"Simulate the Langevin equation (1) with D>0 in a 2D periodic, incompressible velocity field with non-trivial iso-velocity geometry, e.g., u=(sin y, sin x) on [0,2π]^2, starting from a localized Gaussian blob of particles. Estimate F(v,t) by particle histograms at several times. Compute the equilibrium conditional coefficients using Eq. (14) (level-set integrals of G and GG^T) and solve the closed Eq. (10) with the same initial F. Compare the two evolutions; if the L2 relative error exceeds, say, 10% at any time before the distribution approaches equilibrium, the equilibrium closure is not adequate for localized plumes. This directly tests the load-bearing assumption.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim—that Eq. (10) provides a rigorous evolution equation for the Lagrangian velocity PDF—is technically correct, but the equation is unclosed: the conditional expectations G^{x|v}, ∆u^{x|v}, and GG^T^{x|v} depend on the evolving concentration field c(x,t) (Eq. 11). The paper's equilibrium closure (Eq. 14) replaces these with long-time, uniform-concentration averages. This is valid only when c is nearly uniform on each iso-velocity surface, a condition met at asymptotic times, not for the localized plumes that motivate the method. The only numerical test (Fig. 1) solves the closed equation and checks that its stationary limit matches the Eulerian PDF; this is a necessary consistency check but cannot distinguish a correct closure from any other that shares the same equilibrium. No comparison to direct particle simulation or to the exact unclosed dynamics is provided. In heterogeneous media with non-invertible flow maps (multiple x for a given v), the closure error is likely to be significant at finite times. Thus the paper does not yet establish that the proposed closed model predicts the time-dependent Lagrangian velocity PDF in the intended applications.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a probability-density-function (PDF) framework for particle transport in steady heterogeneous velocity fields. Starting from the Itō SDE dX = u(X) dt + sqrt(2D) dW, the authors apply Itō's lemma to U = u(X), obtain a joint position-velocity Fokker-Planck equation, and spatially marginalize it into an equation for the Lagrangian velocity PDF F(v,t) with conditional-expectation coefficients. They discuss closures: an equilibrium closure replacing conditional expectations by uniform-concentration values, a perturbation expansion around equilibrium, and an IEM/Gaussian closure. For two-dimensional shear flows they derive reduced equations, and for Hagen-Poiseuille flow they obtain the closed equation (24) whose stationary solution is the Eulerian velocity PDF, with a numerical illustration of convergence.","tokens_in":14282,"tokens_out":23311,"duration_ms":227251,"significance":"The paper offers a useful and well-motivated route from a particle-resolved SDE to velocity-space transport equations for Lagrangian velocity statistics in porous media, a problem of current interest for CTRW and anomalous-dispersion modeling. The Poiseuille reduction (Eqs. 22-25) is a clean exact benchmark, and the conditional-expectation and IEM closure framework is a sensible adaptation of standard PDF methods. The derivation is, however, an application of textbook Itō calculus rather than a fundamentally new technique, and the paper currently contains internal algebraic inconsistencies in the central equations (see major comments). No code or machine-checked proofs are provided; the numerical test is illustrative. If the sign and cross-derivative issues are repaired and the equilibrium closure is tested in a genuinely heterogeneous flow, the framework would be a solid contribution.","major_comments":[{"comment":"The sign of the DΔu term is incorrect. Marginalizing the conservative form (7), whose velocity-space divergence is ∇v·[(-Gv - DΔu + DGG^T∇v)f], gives ∂F/∂t = -∇v·[(G^{x|v} v + DΔu^{x|v})F] plus the diffusion term, not -∇v·[(G^{x|v} v - DΔu^{x|v})F]. The printed sign is consequential: inserting the Poiseuille data into Eq. (10) with Gv=0 and Δu = -2u0/L² yields a drift -2Du0/L² ∂F/∂u, opposite to the +2Du0/L² ∂F/∂u in Eq. (24), and the stationary solution is then not Eq. (25). Equation (22), which the authors actually solve, uses the correct sign, so the paper is internally inconsistent with its own central equation.","section":"Section 3, Eqs. (9)-(10), (15)-(17)"},{"comment":"The three displayed forms of the joint Fokker-Planck equation are not equivalent as written. With G = ∇u as defined, the Itō SDE (5) has a noise matrix whose off-diagonal block is G^T, so the mixed derivative in the standard Fokker-Planck equation is 2D ∇x·∇v:(G^T f), not 2D ∇x∇v:[G f] as in Eq. (6). In addition, Eq. (8), claimed to be an equivalent non-conservative form, has no mixed derivative at all; expansion of the correct Eq. (6) retains a mixed second derivative and an additional drift containing ∇x·G. The mixed term vanishes under the spatial marginalization used in Section 3 for periodic or infinite domains, so the marginal equation can be repaired, but the joint-PDF equation and the equivalence statement need correction.","section":"Section 2, Eqs. (6)-(8)"},{"comment":"The numerical test does not validate the equilibrium closure that is the paper's main closure tool. For the Hagen-Poiseuille profile, the map y ↦ u(y) is invertible, so the conditional density f(y|u;t) is δ(y-y(u)) for all times; consequently M_{σ'}(u;t)=σ' and M_{σ²}(u;t)=σ²(y(u)) are time-independent even without the equilibrium closure, and Eq. (22) is exact. Fig. 1 therefore only demonstrates convergence of the exact reduced equation to its stationary solution. It gives no evidence about the accuracy of closure (14) in heterogeneous, non-invertible flow fields, where the conditional expectations in Eq. (10) depend on the evolving concentration. The authors should either add a direct comparison against particle tracking (or the exact unclosed dynamics) in a non-invertible flow, or explicitly limit the claims of predictive closure to the shear-flow cases.","section":"Section 3.2 and Section 5, Fig. 1"}],"minor_comments":[{"comment":"The phrase 'in in deterministic' should read 'in deterministic'.","section":"Abstract"},{"comment":"For consistency with SDE (5), the x-space advection term should be u(x)f, not v f; the equality v f = u(x)f holds only on the support of the joint delta distribution.","section":"Eq. (6)"},{"comment":"Clarify that the transverse velocity component v in the shear-flow dynamics is taken to be zero; otherwise the statement that 'the term σv ∂f/∂u can be disregarded' is not justified.","section":"Section 4, Eqs. (20)-(21)"},{"comment":"The sentence 'Since f*(x|v)=0' is confusing because f* is not a probability density; rephrase as an explicit assumption that the conditional-density fluctuation vanishes.","section":"Section 3.3"},{"comment":"State the no-flux boundary conditions in velocity space and give the numerical time-stepping details for the Chebfun solution; the paper cites Chebfun but provides no discretization or convergence information.","section":"Fig. 1 and Eq. (24)"},{"comment":"The claim 'For the first time, in this work, we propose a rigorous approach...' should be softened, since velocity-PDF/Fokker-Planck derivations based on Itō calculus and conditional expectations have a long history in turbulence and stochastic coarse-graining.","section":"Section 5"},{"comment":"Typographical errors: 'Ito' for 'Itō' throughout, 'Smoluchowki' in Section 5, 'would noto be fully closed' in Section 3.4, and inconsistent notation for dimensionless times in Fig. 1.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The sign error in Eq. (10) is serious but clearly local and fixable; the paper's own correct Poiseuille derivation indicates the intended form. I recommend major revision rather than rejection. The authors should also temper the 'first time' novelty claim; the derivation is textbook Itō calculus applied to a specific problem. If the authors can add one validation in a heterogeneous non-invertible flow, the paper would be considerably stronger."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear [Colleague],\n\nPunchline: this paper is a mostly clean derivation of a marginal equation for Lagrangian velocity PDFs in porous media, but the only test of its practical closure is a self-consistency check, not a comparison against real dynamics. It deserves serious review, not a desk reject.\n\nWhat is new: the authors start from the Ito SDE for (X, u(X)), derive a joint position-velocity Fokker-Planck equation, then marginalize over space to get Eq. (10), an evolution equation for F(v,t) with conditional-expectation coefficients. That framing is new in the porous-media literature. The Poiseuille example in Sec. 4.2 is a nice sanity check: under the equilibrium closure the stationary solution reproduces the known Eulerian velocity PDF. The paper is also honest with itself—it names the closure problem and flags future work.\n\nWhere it is soft. First and most important: the stress-test concern is correct. Eq. (10) is unclosed; the practical model is the equilibrium closure of Sec. 3.2, which is valid only when concentration is nearly uniform on iso-velocity surfaces. The only numerical test, Fig. 1, solves the closed equation and shows relaxation to equilibrium; that is necessary but not sufficient. It cannot discriminate between a correct closure and any other closure with the same stationary solution. For localized plumes in heterogeneous media, there is no evidence the finite-time dynamics are captured. Second, there is an apparent error in the cross-diffusion term of the joint FP equation (Eq. 6). It does not affect the marginal equation, since it integrates to a boundary term, but Eqs. (6)–(8) are presented as equivalent and at least one is miswritten. Third, the claim of doing this 'for the first time' is a bit strong: the math is standard stochastic calculus and conditional-expectation closures are common in turbulence PDF methods. The novelty is the application, not the tool.\n\nBottom line: a serious theoretical contribution for people building velocity-evolution and CTRW models in porous media, but it needs a corrected Eq. (6) and at least one test of the equilibrium closure against direct particle simulation in a simple heterogeneous medium before it can be called validated. I'd bring it to a reading group and would cite the derivation.\n\nBest,\n[Your name]","headline":"A worthwhile theoretical derivation of a velocity-PDF evolution equation for porous media, but the equilibrium closure is tested only by self-consistency and the joint FP has a likely typo; deserves revision rather than rejection.","tokens_in":14766,"tokens_out":6179,"would_cite":true,"duration_ms":64038,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper derives, for the first time, a rigorous evolution equation for the Lagrangian velocity PDF in porous media, using stochastic calculus rather than empirical velocity models.","keywords":["porous media","Lagrangian velocity PDF","Fokker-Planck equation","stochastic differential equations","hydrodynamic dispersion","continuous time random walk","conditional expectations","Poiseuille flow"],"falsifier":"Run pore-scale particle-tracking simulations in a heterogeneous 2D or 3D porous medium starting from a spatially localized injection, compute the marginal Lagrangian velocity PDF over time, and compare it with the solution of the equilibrium-closed equation (24 or its generalisation); if the predicted relaxation to the Eulerian velocity PDF is systematically slower or faster than the simulated one before the concentration becomes uniform, the closure is falsified.","tokens_in":13775,"feed_emoji":"💧","tokens_out":6542,"duration_ms":64535,"temperature":0.7,"pith_summary":"The authors aim to replace empirical models of how particle velocities change as solutes drift through heterogeneous porous media with a derivation from stochastic calculus. They extend the usual position-only description to a joint position–velocity probability density, write the corresponding Itô/Fokker-Planck equation, and then average over space to obtain an evolution equation for the marginal Lagrangian velocity PDF. In the equilibrium closure—valid when the solute concentration has spread uniformly—the equation closes and, for Poiseuille flow, its stationary solution is $F(v)\\propto(1-u/u_0)^{-1/2}$, exactly the Eulerian velocity PDF. A sympathetic reader would care because this supplies a principled basis for continuous-time-random-walk velocity models and for closures in anomalous dispersion.","feed_headline":"Porous-flow velocity distributions derived from stochastic calculus","feed_subtitle":"A Fokker-Planck derivation reproduces the Eulerian velocity PDF in Poiseuille flow, replacing empirical velocity models.","key_machinery":"The carrying object is the joint position–velocity PDF $f(x,v;t)$, which restores trajectory information lost in the concentration-only advection-diffusion description. Its evolution is generated by the Itô expansion $dU=\\nabla u\\,dX+\\tfrac12 dX^T\\nabla^2 u\\,dX$ of the Lagrangian velocity, producing the coupled SDE system (5) and the Fokker-Planck equation (7) with velocity-space drift $(Gv + D\\Delta u)$ and diffusion $D GG^T$. Marginalising over space moves the closure problem into the velocity-space terms, which are expressed through conditional expectations $G^{x|v}$, $\\Delta u^{x|v}$, $GG^T{}^{x|v}$; Bayes's rule turns these into concentration-weighted integrals over iso-velocity surfaces. The equilibrium closure replaces the concentration by its uniform long-time value, making the coefficients functions of $v$ alone; in shear flows the marginal equation reduces to an explicit diffusion equation in velocity space whose stationary solution is the Eulerian velocity PDF.","core_discovery":"The central claim is that the marginal Lagrangian velocity PDF $F(v,t)$ satisfies a Fokker-Planck-type equation (Eq. 10) whose drift and diffusion coefficients are conditional expectations of the velocity gradient and Laplacian given velocity. Because the position PDF is just the concentration field and the velocity is a deterministic function of position, these conditional expectations can be written as integrals over the iso-velocity surfaces weighted by concentration (Bayes formula, Eq. 11). Under the equilibrium closure, the coefficients become the time-independent Eulerian conditional averages of Eq. (14), yielding the closed equation (24) in two-dimensional Hagen-Poiseuille flow. The paper shows the stationary solution of this closed equation is $F\\propto(1-u/u_0)^{-1/2}$, which matches the Eulerian velocity PDF obtained by uniform sampling of the parabolic profile. This establishes, on the paper's own terms, a first-principles route from the advection-diffusion equation to Lagrangian velocity statistics, without assuming a velocity Markov model.","pith_inferences":["I would expect the unclosed time-dependent conditional expectations to matter most in the early life of a localized plume; a direct particle-tracking test in a three-dimensional pore geometry, comparing measured $F(v,t)$ with the equilibrium-closed solution, would show how fast the closure becomes valid.","The same marginalisation could be applied to filtered (locally averaged) equations, not just global averages; that would produce a coupled system where concentration heterogeneities feed back into velocity statistics—an untested extension the paper leaves open.","Because the stationary solution is the Eulerian velocity PDF, the paper implies that any discrepancy between Lagrangian and Eulerian velocity statistics must be a transient, diffusion-driven relaxation phenomenon; this could be checked experimentally with particle tracking in bead packs.","The analogy with turbulent mixing suggests standard mixing closures (interaction-by-exchange-with-the-mean type) could be ported to porous-media velocity PDFs, yielding relaxation-rate parameters that are measurable from acceleration statistics."],"forward_implications":["The evolution of Lagrangian velocity statistics in steady heterogeneous flows can in principle be computed from the flow field and its gradients, not fitted to data.","The derived equation gives CTRW dispersion models a rigorous velocity-evolution ingredient instead of an ad hoc one.","In Poiseuille-type flows the equilibrium closure is exact in the long-time limit, so the closed equation reproduces the known Eulerian velocity distribution $F(v)\\propto(1-u/u_0)^{-1/2}$.","Because the spatial advection term needs no closure, the approach shifts the modelling difficulty into velocity-space conditional averages, where explicit closures such as the IEM approximation can be applied.","The formulation extends naturally to inertial particles through a drag-force SDE with Stokes number, opening a route to second-order velocity models."],"supporting_citations":[{"why":"supplies the Itô calculus and Fokker-Planck construction used to derive the joint position-velocity equation.","marker":"[27]"},{"why":"gives the standard stochastic-processes reference for the Fokker-Planck equation and Ito calculus.","marker":"[18]"},{"why":"provides the Poiseuille-flow Eulerian velocity distribution that the equilibrium stationary solution is designed to reproduce.","marker":"[47]"},{"why":"motivates the problem by modelling Lagrangian velocity evolution within a continuous-time random walk framework.","marker":"[11]"},{"why":"offers an empirical stochastic model for Lagrangian velocity series that the derived equation aims to ground.","marker":"[36]"},{"why":"introduces the PDF-method and interaction-by-exchange-with-the-mean closure conventions the paper adapts.","marker":"[43]"}],"fun_headline_variants":["Exact velocity PDF without velocity Markov assumption","Fokker-Planck derivation reproduces Eulerian velocity PDF","From advection-diffusion to Lagrangian velocity statistics","New closure yields exact dispersion velocity distribution","Stochastic model for velocity PDF, no empirical input"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The equilibrium closure assumes the solute concentration has already spread uniformly across the medium, so that time-dependent conditional averages can be replaced by their long-time, concentration-independent values; for a localized plume in a heterogeneous medium this is not yet justified, and the paper offers no direct simulation of that regime.","fun_headline_variants_meta":{"raw":{"variants":["Exact velocity PDF without velocity Markov assumption","Fokker-Planck derivation reproduces Eulerian velocity PDF","From advection-diffusion to Lagrangian velocity statistics","New closure yields exact dispersion velocity distribution","Stochastic model for velocity PDF, no empirical input"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000216,"raw_usage":{"total_tokens":1465,"prompt_tokens":1008,"completion_tokens":457,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":624,"completion_tokens_details":{"reasoning_tokens":383}},"tokens_in":624,"tokens_out":457,"duration_ms":5173,"temperature":1.0,"reasoning_tokens":383,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:14:31.931929+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run pore-scale particle-tracking simulations in a heterogeneous 2D or 3D porous medium starting from a spatially localized injection, compute the marginal Lagrangian velocity PDF over time, and compare it with the solution of the equilibrium-closed equation (24 or its generalisation); if the predicted relaxation to the Eulerian velocity PDF is systematically slower or faster than the simulated one before the concentration becomes uniform, the closure is falsified.","supporting_citations":[{"cited_title":"Stochastic Methods: A Handbook for the Natural and Social Sciences, 2009","cited_arxiv_id":null,"evidence_quote":"gives the standard stochastic-processes reference for the Fokker-Planck equation and Ito calculus."},{"cited_title":"Pope and S.B","cited_arxiv_id":null,"evidence_quote":"provides the Poiseuille-flow Eulerian velocity distribution that the equilibrium stationary solution is designed to reproduce."},{"cited_title":"Continuous time random walks for the evolution of la- grangian velocities","cited_arxiv_id":null,"evidence_quote":"motivates the problem by modelling Lagrangian velocity evolution within a continuous-time random walk framework."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"introduces the PDF-method and interaction-by-exchange-with-the-mean closure conventions the paper adapts."}],"review_version":1}