{"id":"6da1e327-0393-455b-80c4-a1301472c02c","arxiv_id":"2411.09314","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"The paper derives relaxation-rate conditions that eliminate orientation-dependent advection errors in lattice Boltzmann models.","lead":"This paper analyzes the lattice Boltzmann method and derives parameter conditions that remove directional errors when moving fluid patterns are advected. The result gives practitioners explicit tuning rules for more accurate simulations of advection-diffusion and simple fluid flows.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (47) as printed states sigma6 = sigma8 = (1/12)sigma4, but this does not cancel the D2Q13 angular advection term; substitution into Eq. (44) leaves 10055(sigma4^2 - 1)/157080 V sin(4theta).","rationale":"The reader's conditional verdict is justified, and the strongest specific worry is sharper than the nonlinear-extrapolation caveat. The full text states the D2Q13 condition in a form that fails against the paper's own Eq. (44). Since this condition is one of the two headline results in the abstract, the manuscript cannot be accepted as is; it needs a corrected equation or a demonstration that a different condition was used. This is a local, checkable algebraic issue rather than a judgment call about how far linear analysis extends. I do not see evidence of a manufactured result: the reciprocal condition appears to be what the calculation requires, so the paper is probably a typographical/convention error away from its intended claim, and the D2Q9 half (Eq. 35) is internally consistent. However, the printed inconsistency plus the acknowledged linear-only analysis and missing derivations keep the appropriate verdict at CONDITIONAL. If the corrected condition passes the symbolic re-derivation and a tuned D2Q13 simulation is shown, the remaining concern would be the linear-to-nonlinear extension, which is the reader's original caveat and not contradicted by this report.","tokens_in":16745,"tokens_out":8597,"duration_ms":78538,"concrete_test":"Independently substitute both candidate forms of Eq. (47) into Eqs. (44) and (46). With the printed form sigma6 = sigma8 = (1/12)sigma4, the sin(4theta) coefficient in Eq. (44) is 10055(sigma4^2 - 1)/157080 and the f2(theta) coefficient in Eq. (46) is (31/102)(sigma4^2 - 1); with the reciprocal form sigma6 = sigma8 = 1/(12 sigma4), both vanish and Eq. (48) follows. Use a symbolic re-derivation of the D2Q13 linearized dispersion relation from Sections 4-5 to determine which form is correct, then print the correct formula and rerun the D2Q13 plane-wave comparison in Section 6 with that condition.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing claim for D2Q13 is the tuning condition (47). Read literally, it says sigma6 = sigma8 = (1/12)sigma4. Substitution into the transverse-wave phase velocity (44) gives v_phi = [10055(sigma4^2 - 1)/157080] V sin(4theta), so the sin(4theta) anomaly survives for generic sigma4. The same substitution in (46) leaves a nonzero f2(theta) term, (31/102)(sigma4^2 - 1) V f2(theta), plus a further nonzero residual. The condition that actually cancels the angular terms is sigma6 = sigma8 = 1/(12 sigma4): this reciprocal form kills the sin(4theta) term in (44), cancels the f2(theta) part of (46), and reproduces (48). If the reciprocal form is what the authors used, the displayed formula needs correction; if the literal form is meant, the central D2Q13 isotropy result is false. Either way, the announced 'situations accurate to third order' rest on an equation that, as printed, is internally inconsistent with the paper's own dispersion analysis. The D2Q9 condition (35) does check out: with sigma3 = sigma4 and 12 sigma4 sigma6 = 1, the coefficients g1, g2, g3 and h1, h3 vanish in (30)-(31). The nonlinear-validity caveat stated in the conclusion remains a separate restriction.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes discrete, orientation-dependent advection errors in lattice Boltzmann schemes, focusing on odd-order space derivatives. Using the so-called Berlin algorithm for Taylor expansion of the LBE in moment space, it derives equivalent equations up to third order in space derivatives for advection-diffusion (D2Q5, D3Q15, D3Q19) and athermal fluids (D2Q9, D2Q13), and proposes relaxation-rate conditions, expressed in Hénon parameters, that remove the angular dependence of the advection phase velocity. For D2Q9 athermal fluids the proposed conditions are Eq. (35), σ3 = σ4 and σ4σ6 = 1/12; for D2Q13 the paper proposes Eq. (47) and then states the resulting phase velocity Eq. (48). The paper also presents plane-wave dispersion comparisons and Gaussian vortex simulations. I verified that the D2Q9 condition (35) does make the third-order matrix isotropic as claimed, but the D2Q13 condition as printed is internally inconsistent; the correct condition is the reciprocal form σ6 = σ8 = 1/(12σ4), as explained below.","tokens_in":17054,"tokens_out":8110,"duration_ms":72571,"significance":"If the D2Q13 condition is corrected, the paper is a useful contribution: it gives explicit, analytic Hénon-parameter conditions that remove orientation-dependent third-order advection errors, rather than fitting parameters to simulations, and it supports the linear analysis with plane-wave simulations agreeing to within 0.15 percent. The D2Q9 result (35)-(37) is verified analytically and is a genuine strength. The main weaknesses are that the D2Q13 central condition is misprinted as written, the third-order matrices and final formulas are asserted without derivation, and the claimed third-order accuracy is established only in the linearized regime; the finite-amplitude simulations are qualitative consistency checks. These issues are fixable, but the printed central claim for D2Q13 is currently false.","major_comments":[{"comment":"Equation (47) as printed, σ6 = σ8 = (1/12)σ4, does not cancel the angular terms in the D2Q13 dispersion relations. Substituting this condition into Eq. (44) leaves vφ = [10055(σ4^2 - 1)/157080] V sin(4θ), and substituting it into Eq. (46) leaves a nonzero f2(θ) term plus a further nonzero velocity-dependent term; Eq. (48) does not follow. The condition that cancels the sin(4θ) term in Eq. (44), removes the f2(θ) term in Eq. (46), and reproduces Eq. (48) is σ6 = σ8 = 1/(12σ4). This is obtained directly by solving Eq. (44) with σ6 = σ8 = σ, which gives 120660 σ4 σ = 10055, i.e., σ4 σ = 1/12. As printed, the central D2Q13 isotropy claim is false, so Eq. (47) must be corrected and the plane-wave table in Section 6 and Figure 4 re-examined under the intended condition.","section":"Section 5, 'Athermal fluid simulated with D2Q13', Eq. (47)"},{"comment":"The third-order matrix N3 for D2Q9 and the D2Q13 phase-velocity formulas are stated as results of the Berlin algorithm without derivation; the text says only that 'higher orders Ml are cumbersome and not given here.' Given the inconsistency in Eq. (47), the reader cannot check whether the stated formulas are correct without redoing the entire expansion. Please include the explicit N3 matrices for both D2Q9 and D2Q13, or provide a supplementary derivation or reproducible code that generates these formulas.","section":"Section 5, Eqs. (30)-(31) and (44)-(46)"},{"comment":"The abstract states that the paper proposes situations 'accurate to third order in space derivatives,' but Sections 4 and 5 analyze only linearized plane waves about a uniform velocity, as the conclusion itself acknowledges. The Gaussian vortex simulations in Figures 3 and 4 are finite-amplitude and are interpreted qualitatively, without an error measure tied to the third-order analysis. The third-order accuracy claim should be explicitly qualified as holding in the linearized regime, and the simulations should be described as consistency checks rather than validations of third-order accuracy for nonlinear flows.","section":"Abstract and Section 6"}],"minor_comments":[{"comment":"There are numerous typographical errors, e.g., 'generationg' on page 2, 'extention', 'usefull', 'ancellation', and 'basedof' in the appendices; the manuscript would benefit from a careful proofreading pass.","section":"Introduction and throughout"},{"comment":"The domain sizes '1012' and '3632' should likely read '101×101' and '363×363'; please clarify the notation.","section":"Section 6, Figures 2 and 4"},{"comment":"The table for D2Q13 plane-wave advection does not list the relaxation parameters (σ4, σ6, σ8, c1, q) used in the simulation. Please provide these values so the reader can verify which isotropy condition was actually applied.","section":"Section 6, plane-wave table"},{"comment":"The caption of Figure 1 does not define the plotted quantity h nor identify which curve is solid and which is dashed; please make the definition and legend explicit.","section":"Figure 1 and Eq. (49)"},{"comment":"'relaxation rates5' should be 'relaxation rate s5', and the table caption should state which moments correspond to the non-conserved relaxation rates s5, s6, s11, s14, s16, and s17.","section":"Appendix 2, Eq. (66) and Table 7"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern is confirmed: Eq. (47) as printed is internally inconsistent with the paper's own dispersion formulas. The reciprocal condition σ6 = σ8 = 1/(12σ4) is the one that makes Eqs. (44), (46), and (48) mutually consistent, and the plane-wave table in Section 6 appears to have been computed with that reciprocal form. This is therefore fixable, but the manuscript as submitted cannot be accepted with the central D2Q13 claim stated as it is. The absence of derivations for N3 and the D2Q13 formulas contributed to this error and should be addressed in revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Punchline: this is a real extension of the Lallemand–Dubois–Luo rotational-invariance program, redirecting it to odd-order advection defects, and the D2Q9 part is clean. But the D2Q13 headline condition (47) is wrong as printed: with sigma6 = sigma8 = (1/12)sigma4, Eqs. (44) and (46) retain (10055/157080)(sigma4^2 - 1)V sin(4theta) and (31/102)(sigma4^2 - 1)V f2(theta) plus other residuals. The condition that actually cancels the angular terms is sigma6 = sigma8 = 1/(12sigma4). I checked the arithmetic; the stress-test note lands. Either the display has a typesetting ambiguity or the central D2Q13 isotropy claim is false.\n\nCredit where it is due. The paper gives explicit third-order anomalous-advection coefficients for D2Q9, D2Q13, D3Q15 and D3Q19 rather than fitting them to data; the D2Q9 condition (35) checks out and yields the claimed isotropic N3. The plane-wave table for D2Q13 shows simulation agreeing with theory to within 0.15%, which is meaningful evidence for the linear analysis. The Gaussian-dot and vortex pictures illustrate the orientation defect clearly. The authors also honestly state in the conclusion that the analysis is linearized, so the nonlinear validity of the tuning conditions is flagged as open rather than hidden.\n\nSoft spots, in proportion. The third-order matrices N3 and the final isotropy formulas are asserted, not derived; the manuscript refers to earlier papers. That is a self-containedness problem for a paper whose whole point is a new tuning recipe. No code or data are shipped, so the reported numbers are not independently checkable except by reimplementation. The most serious issue is Eq. (47), because it is load-bearing. In addition, the D2Q13 phase-velocity expression (48) does not follow from (46) even with the reciprocal substitution unless c1 = 0, so there may be a second typo in the constant term. The D3 appendices have the same asserted-formula style.\n\nWho should read it: LBM practitioners who want to remove grid-orientation artifacts in advection-dominated simulations will get value from the D2Q9 recipe and from the general framework. The paper deserves a serious referee: the framework is sound, the linear validation is real, and the issues look correctable rather than fatal. A referee should require the derivation of N3, corrected displays, and a note on the reciprocal versus literal reading of (47) before acceptance.","headline":"Useful extension of the authors' LBM advection-error program, but the D2Q13 tuning condition is inconsistent as printed and needs correction before the paper can be trusted.","tokens_in":17618,"tokens_out":4909,"would_cite":false,"duration_ms":45581,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65M75","76M28"],"pacs":["02.70.Ns","05.20.Dd","47.11.+j"],"model":"deepseek-v4-flash","headline":"The paper shows that the directional bias of lattice Boltzmann advection is a third-order error that can be removed by tuning relaxation rates, making linear advection in D2Q9 and D2Q13 rotationally invariant to third order in space…","keywords":["Lattice Boltzmann method","Taylor expansion method","Hénon parameters","anomalous advection","rotational invariance","D2Q9 model","D2Q13 model","advection-diffusion"],"falsifier":"Simulate a Gaussian vortex in D2Q9 with the tuned rates $\\sigma_3=\\sigma_4$ and $\\sigma_4\\sigma_6=1/12$, at a finite Mach number such as 0.1, along several lattice directions, and compare the vorticity after many time steps; residual angle-dependent distortion that grows with the initial amplitude would show that the linear tuning does not control nonlinear advection errors. Equivalently, measure the phase velocity of a finite-amplitude plane wave as a function of angle and amplitude and check whether the $\\sin 4\\theta$ term stays zero outside the linear regime.","tokens_in":16533,"feed_emoji":"🌀","tokens_out":8249,"duration_ms":75573,"temperature":0.7,"pith_summary":"Lattice Boltzmann simulations of moving patterns can distort them depending on the angle between the flow and the lattice grid. This paper attributes that distortion to odd-order space derivatives in the equivalent partial differential equation of the scheme, and shows that the distortion can be tuned away by choosing relaxation rates appropriately. For D2Q9, setting $\\sigma_3=\\sigma_4$ and $\\sigma_4\\sigma_6=1/12$ makes the third-order advection terms independent of angle; for D2Q13, the analogous condition is $\\sigma_6=\\sigma_8=\\sigma_4/12$. With an additional special choice the residual phase-velocity error vanishes, so linear advection becomes accurate to third order in space derivatives. Simulations of an advected Gaussian dot and vortex confirm the rotational symmetry predicted by the tuned parameters.","feed_headline":"Lattice advection can be made rotationally accurate to third order","feed_subtitle":"Specific relaxation rates in D2Q9/D2Q13 remove angle-dependent advection error, confirmed by vortex simulations.","key_machinery":"The argument is carried by the Taylor-expansion method for lattice Boltzmann schemes, which iteratively eliminates time derivatives and produces an equivalent partial differential equation whose successive space-derivative terms expose the discrete errors. Working in moment space with orthogonal polynomial moments, the authors express relaxation rates through parameters $\\sigma_i=1/s_i-1/2$, and use a plane-wave dispersion analysis to isolate the third-order advection correction. The identity that matters is the isotropy condition for the third-order matrix $N_3$: once the $\\theta$-dependent coefficients vanish, the model's advection is rotationally invariant to third order.","core_discovery":"The central claim is that anomalous advection in lattice Boltzmann schemes is a third-order discrete effect, expressible as an orientation-dependent correction to the phase velocity of plane waves, and that this correction can be eliminated algebraically rather than by enlarging the velocity set. In the linearized analysis around a uniform velocity, the correction takes the form $A(\\theta) k^2$, where $\\theta$ is the angle between the wave vector and the grid. The paper derives the exact dependence of $A$ on the relaxation parameters for D2Q9 and D2Q13, and shows that the conditions $\\sigma_3=\\sigma_4$ with $\\sigma_4\\sigma_6=1/12$ (D2Q9) or $\\sigma_6=\\sigma_8=\\sigma_4/12$ (D2Q13) make $A$ independent of $\\theta$. With $\\sigma_4=1/\\sqrt{12}$ in D2Q13 the remaining phase velocity vanishes, giving a model whose linear advection is rotationally invariant and accurate to third order in space derivatives; the same method yields analogous cancellation conditions for D3Q15 and D3Q19 advection-diffusion models.","pith_inferences":["The same tuning conditions should carry over to nonlinear flows only if the neglected nonlinear odd-order terms are small compared with the linear ones; repeating the vortex simulations at higher Mach number would test this directly.","Because the anomalous-advection correction scales as $k^2$ of the wave vector, the rotational distortion is largest for small-scale features, so the tuning matters most in under-resolved or high-wave-number simulations.","The Appendix result that no single choice of rates simultaneously removes anomalous advection and hyper-diffusivity in the 3D TRT diffusion models suggests a practical trade-off: a user must decide whether to prioritize isotropic advection or fourth-order dissipation.","One could apply the same moment-space expansion to other lattice Boltzmann variants, such as D3Q27 or thermal models, to derive their own third-order advection isotropy conditions."],"forward_implications":["In D2Q9, the choices $\\sigma_3=\\sigma_4$ and $\\sigma_4\\sigma_6=1/12$ eliminate the angle dependence of the third-order advection matrix, so a linearly advected Gaussian vortex keeps its rotational symmetry.","In D2Q13, $\\sigma_6=\\sigma_8=\\sigma_4/12$ removes the angular dependence of the additional phase velocity, and $\\sigma_4=1/\\sqrt{12}$ removes the residual phase velocity, yielding third-order accuracy.","For three-dimensional advection-diffusion, D3Q15 and D3Q19 admit two parameter families that suppress anomalous advection, and in the two-relaxation-times case both models share the condition $\\sigma_1=1/\\sqrt{12}$.","The analysis provides explicit formulas for the anomalous-advection factor $g(k)=k\\cdot V(1+h k^2)$, allowing numerical predictions of the advection error for arbitrary wave-vector orientation, as verified in the paper's plane-wave simulations."],"supporting_citations":[{"why":"Supplies the simplified Taylor-expansion algorithm used to iterate the equivalent equations to third order.","marker":"[1]"},{"why":"Establishes the equivalent partial differential equation method for a Boltzmann scheme on which the paper's derivation rests.","marker":"[2]"},{"why":"Provides the higher-order lattice Boltzmann scheme framework and the quartic condition for cancelling hyper-viscosity, used at fourth order.","marker":"[4]"},{"why":"Defines the shifted relaxation parameters $\\sigma_i=1/s_i-1/2$ used to express all the isotropy conditions.","marker":"[6]"},{"why":"Introduces the linear-relaxation moment collision model in which non-conserved moments relax with rates $s_i$.","marker":"[7]"},{"why":"The earlier dispersion, dissipation, isotropy, Galilean invariance, and stability study that this work extends.","marker":"[9]"},{"why":"Gives the exact Gaussian solution for advection-diffusion and vortex decay used as the rotational-invariance test in simulations.","marker":"[10]"}],"fun_headline_variants":["Lattice Boltzmann advection made rotationally invariant to third order","Exact cancellation of angle-dependent advection in LBM","D2Q9 and D2Q13 achieve rotationally invariant advection","Third-order advection error removed by relaxation tuning","LBM advection corrected to third order without extra velocities"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire parameter prescription is derived from a linearized analysis around a uniform velocity, assuming slowly varying smooth fields, and the paper's own conclusion states that only linearized situations were treated; the extension to actual nonlinear vortex flows rests on the unverified assumption that nonlinear odd-order advection errors either vanish or obey the same cancellation conditions.","fun_headline_variants_meta":{"raw":{"variants":["Lattice Boltzmann advection made rotationally invariant to third order","Exact cancellation of angle-dependent advection in LBM","D2Q9 and D2Q13 achieve rotationally invariant advection","Third-order advection error removed by relaxation tuning","LBM advection corrected to third order without extra velocities"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000912,"raw_usage":{"total_tokens":3909,"prompt_tokens":924,"completion_tokens":2985,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":540,"completion_tokens_details":{"reasoning_tokens":2911}},"tokens_in":540,"tokens_out":2985,"duration_ms":20057,"temperature":1.0,"reasoning_tokens":2911,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T20:46:47.874728+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate a Gaussian vortex in D2Q9 with the tuned rates $\\sigma_3=\\sigma_4$ and $\\sigma_4\\sigma_6=1/12$, at a finite Mach number such as 0.1, along several lattice directions, and compare the vorticity after many time steps; residual angle-dependent distortion that grows with the initial amplitude would show that the linear tuning does not control nonlinear advection errors. Equivalently, measure the phase velocity of a finite-amplitude plane wave as a function of angle and amplitude and check whether the $\\sin 4\\theta$ term stays zero outside the linear regime.","supporting_citations":[{"cited_title":"Augier, F","cited_arxiv_id":null,"evidence_quote":"Supplies the simplified Taylor-expansion algorithm used to iterate the equivalent equations to third order."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the equivalent partial differential equation method for a Boltzmann scheme on which the paper's derivation rests."},{"cited_title":"Dubois, P","cited_arxiv_id":null,"evidence_quote":"Provides the higher-order lattice Boltzmann scheme framework and the quartic condition for cancelling hyper-viscosity, used at fourth order."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the shifted relaxation parameters $\\sigma_i=1/s_i-1/2$ used to express all the isotropy conditions."},{"cited_title":"d'Humi\\`eres","cited_arxiv_id":null,"evidence_quote":"Introduces the linear-relaxation moment collision model in which non-conserved moments relax with rates $s_i$."},{"cited_title":"Lallemand, L-S","cited_arxiv_id":null,"evidence_quote":"The earlier dispersion, dissipation, isotropy, Galilean invariance, and stability study that this work extends."},{"cited_title":"Landau, E.M","cited_arxiv_id":null,"evidence_quote":"Gives the exact Gaussian solution for advection-diffusion and vortex decay used as the rotational-invariance test in simulations."}],"review_version":1}