{"id":"9862c88f-cfd1-4417-a3e3-f5ba7690e657","arxiv_id":"2508.06273","paper_version":2,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The paper constructs a third-order fully discrete Active Flux scheme for the 2D Euler equations with positivity-preserving limiting and reflecting boundary conditions.","lead":"A fully discrete Active Flux method for the two-dimensional Euler equations is presented, with new limiting that keeps pressure and density positive. The paper also treats reflecting boundaries and reports accurate results on coarse grids.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Positivity guarantee depends on the approximate evolution operator preserving an invariant region; the abstract gives no evidence this holds for the fully discrete scheme.","rationale":"The reader's weakest assumption focused on the accuracy of the bicharacteristics-based evolution operator for third-order accuracy. That is a valid risk, but the more load-bearing concern for the abstract's strongest claim is whether the positivity guarantee is actually established for the fully discrete approximate operator. The two concerns are related: both depend on properties of the approximate evolution operator, but they target different consequences (accuracy vs. positivity). Since we have only the abstract, we cannot determine whether the guarantee holds; this supports the reader's UNVERDICTED verdict. I do not see a reason to change the verdict, but I have identified a concrete technical condition that the full text must satisfy for the central claim to be credible. The proposed test—checking the invariant-region property of the discrete operator on extreme low-density states—is a standard way to validate such guarantees. The reader's weakest assumption was not exactly the same, hence 'partial' agreement.","tokens_in":525,"tokens_out":3280,"duration_ms":39530,"concrete_test":"Examine the paper's positivity proof (likely in the section on limiting) to find the exact lemma that establishes positivity of the fully discrete update. Specifically, check whether the proof relies on a property of the exact bicharacteristic evolution operator or on a property of the actual discrete approximate operator. Then run the fully discrete scheme on a two-dimensional Riemann problem with a near-vacuum state (e.g., density and pressure both ~1e-6) on a coarse grid. If the approximate operator produces negative density or pressure before limiting, and the limiter is not proven to be applied to the evolved point values before the next time step, the guarantee fails. If the proof instead contains a discrete invariant-region lemma for the approximate operator, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract's strongest claim is that the new limiting strategies guarantee positivity of pressure and density. In any fully discrete Active Flux method, the point values are updated by an approximate evolution operator derived from the method of bicharacteristics. For a positivity guarantee to be valid, either (i) that operator must, by itself, map states with positive density/pressure into states with positive density/pressure, or (ii) the limiter must provably cure any non-positive states produced by the operator, without breaking the fully discrete space-time coupling. The abstract does not state which mechanism is used, nor does it mention any condition on the approximate operator (e.g., a CFL-like invariant-region property). The bicharacteristics method in its exact form has an invariant region, but the fully discrete approximation is exactly where such properties can fail. If the proof of positivity silently assumes the exact operator, or if the limiter is applied to reconstructed states rather than to the raw evolved point values, the guarantee may not hold for the actual fully discrete algorithm. This is a load-bearing concern because the headline contribution is the positivity guarantee, not just empirical accuracy.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript (arXiv:2508.06273) proposes a fully discrete, truly multidimensional Active Flux method for the two-dimensional Euler equations. The abstract claims third-order accuracy with compact stencils in space and time, using exact or approximate evolution operators derived from the method of bicharacteristics for point-value updates. It introduces new limiting strategies that guarantee positivity of pressure and density, discusses implementation of reflecting boundary conditions, and reports accurate approximations on coarse grids. The reviewable material provided consists only of the abstract; no derivations, algorithms, or numerical results beyond the abstract are available.","tokens_in":793,"tokens_out":4263,"duration_ms":46605,"significance":"If the claims hold, this would be a meaningful contribution to Active Flux methods for multidimensional hyperbolic conservation laws. A compact, fully discrete third-order method with positivity preservation for the Euler equations and reflecting boundary conditions would be valuable. However, the significance cannot be assessed from the abstract alone; the central claims require verification in the full text. The lack of convergence tables, error norms, proof details, and numerical comparisons in the abstract prevents any substantive evaluation.","major_comments":[{"comment":"The claim of third-order accuracy is unsupported in the provided text. No convergence study, error norms, or order verification is presented. For a numerical method paper, this is a load-bearing assertion. If the full text contains the verification, the abstract should point to it; otherwise, such evidence must be added before the claim can be assessed.","section":"Abstract, sentence 3"},{"comment":"The positivity guarantee is not qualified. The fully discrete update uses an approximate evolution operator (method of bicharacteristics). For the guarantee to hold, either the operator must preserve positive density and pressure, or the limiter must provably cure non-positive states produced by the operator under the scheme's CFL condition. The abstract does not specify the mechanism or the conditions. This is load-bearing because positivity is a headline contribution. The stress-test concern that the approximate operator may not preserve invariant regions is directly relevant and must be addressed in the full text.","section":"Abstract, sentence 4"},{"comment":"The numerical evidence is described only as 'accurate approximates on coarse grids.' No test problems, error measures, computational cost, or comparisons with existing methods are given. This prevents any assessment of the strength of the numerical results and is a central piece of missing support for the paper's claims.","section":"Abstract, sentence 5"}],"minor_comments":[{"comment":"Typo: 'approximates' should be 'approximations'.","section":"Abstract, sentence 5"},{"comment":"The abstract does not state the CFL condition or the nature of the approximate evolution operator (e.g., order of approximation). Clarifying this would help the reader understand the scheme's practical regime.","section":"Abstract, sentence 2"},{"comment":"The term 'truly multidimensional' should be defined or cited, as it may carry a specific meaning in the Active Flux literature.","section":"Abstract, sentence 1"},{"comment":"No references are given; at minimum, the full text should cite prior Active Flux methods and recent work on positivity-preserving limiting for hyperbolic conservation laws.","section":"Abstract"}],"recommendation":"uncertain","confidential_remarks":"The manuscript as provided to the referee contains only the abstract. My review is necessarily limited to that material. The stress-test concern about the positivity guarantee is serious and should be communicated to the authors: the abstract must specify whether the guarantee is proven for the fully discrete approximate evolution operator or only for the exact bicharacteristics semigroup, and under what CFL-type conditions. If the full text already addresses this, the abstract should be revised to reflect it. I recommend obtaining the full manuscript before making an editorial decision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"I've only seen the abstract, so this is provisional. The paper advertises new limiting strategies that guarantee positivity of density and pressure for a fully discrete third-order Active Flux method on 2D Euler, plus reflecting boundary conditions. That is a genuine addition to a niche but active program. The method of bicharacteristics for point value updates is standard in this community, and the coarse-grid accuracy claim is plausible—Active Flux methods are designed for that. The limiter design and boundary implementation are the actual contributions, and they address a real failure mode.\n\nThe soft spot is the word 'guarantee.' For positivity to hold in the fully discrete scheme, either the approximate evolution operator itself preserves positivity, or the limiter provably restores it without breaking the compact space-time stencil. The abstract doesn't say which, and that's the load-bearing question. The stress-test concern is right: the exact bicharacteristics have an invariant region, but the fully discrete approximation is exactly where that property can fail. I'd want a precise statement of what the limiter does to the evolved point values and under what CFL-like condition the guarantee holds. Also, there's no numerical evidence in the abstract—no error norms or convergence tables—so the third-order claim is unverified here. That's normal for an abstract, but it means the referees have real work to do.\n\nOverall, this looks like a solid paper for the Active Flux / finite volume crowd. The authors are credible and the claim is plausible. I wouldn't cite it myself, but it deserves a serious referee: the positivity proof and the numerical verification should be checked carefully. If the guarantee holds for the fully discrete operator, this is a useful robustness improvement.","headline":"Solid incremental Active Flux paper; the positivity guarantee needs checking against the fully discrete operator.","tokens_in":1195,"tokens_out":2231,"would_cite":false,"duration_ms":24208,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65M08","65M12","76N15"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper presents a fully discrete, third-order Active Flux method for the 2D Euler equations in which new limiting strategies guarantee positivity of density and pressure while retaining accurate coarse-grid solutions.","keywords":["Active Flux method","finite volume method","hyperbolic conservation laws","Euler equations","positivity-preserving limiter","method of bicharacteristics","fully discrete scheme","third order accuracy"],"falsifier":"Run a standard 2D Riemann problem or double Mach reflection on a fine uniform mesh and check whether any cell ever attains negative density or pressure before the final time; a single negative cell disproves the positivity guarantee. Alternatively, compute the convergence order on a smooth isentropic vortex: an observed order substantially below three would contradict the third-order claim.","tokens_in":475,"feed_emoji":"🧮","tokens_out":4524,"duration_ms":49168,"temperature":0.7,"pith_summary":"This paper is trying to establish that the Active Flux framework, which evolves both cell averages and point values, can be made both third-order accurate and physically safe for the two-dimensional compressible Euler equations. The authors add limiting strategies that guarantee density and pressure stay positive, and they describe a way to impose reflecting boundary conditions inside the same multidimensional update. The motivation is practical: a scheme that stays accurate on coarse grids and never produces negative densities or pressures can be used for demanding gas-dynamics simulations without excessive mesh refinement or ad hoc fixes. The paper supports this with numerical experiments on standard test problems.","feed_headline":"Third-order active flux scheme keeps density and pressure positive","feed_subtitle":"Compact fully discrete solver for the 2D Euler equations stays accurate on coarse grids while blocking nonphysical states.","key_machinery":"The central object is the Active Flux reconstruction, where cell averages and point values are updated independently, and the point values are advanced using an exact or approximate evolution operator obtained from the method of bicharacteristics. This operator carries genuinely multidimensional wave information that would otherwise be lost in dimension-split finite volume methods. On top of this sit the two new limiting strategies, which control undershoots in density and pressure, and the reflecting-boundary implementation that feeds consistent one-sided data into the same update.","core_discovery":"The central claim is that the fully discrete, compact-stencil Active Flux method can be made reliable for 2D Euler flows by pairing the bicharacteristic-based point-value update with positivity-preserving limiters. The new limiters are constructed so that the numerical density and pressure cannot drop to nonphysical negative values, while the cell-average/point-value structure retains third-order accuracy. The paper also shows how to impose reflecting boundary conditions consistently within the genuinely multidimensional evolution. Numerical tests indicate that the resulting scheme yields accurate approximations even on relatively coarse grids.","pith_inferences":["The positivity-limiting strategy is formulated around conserved quantities and so could likely be carried over to other hyperbolic systems with similar convex invariant regions, such as the shallow-water or MHD equations.","The method's accuracy hinges on the fidelity of the bicharacteristic evolution operator; a systematic study of exact versus approximate operator choices would clarify the cost-accuracy trade-off for practical use.","Because the update is fully discrete, the scheme may be particularly natural for adaptive or locally refined meshes, though the paper does not test this."],"forward_implications":["Density and pressure positivity is guaranteed by construction, so simulations that would otherwise fail with nonphysical states can run to completion.","Third-order accuracy is achieved compactly in space and time without dimensional splitting, preserving genuinely multidimensional wave information.","Coarse-grid accuracy lowers the resolution needed for a given error target, reducing computational cost in practice.","Reflecting boundary conditions expand the method's applicability to wall-bounded compressible flows, such as internal aerodynamics."],"supporting_citations":[],"fun_headline_variants":["2D Euler active flux with guaranteed positive states","Active flux keeps density and pressure positive in 2D","Third-order active flux for 2D Euler with positivity limiter","Compact active flux avoids negative pressure for 2D Euler","Positivity-preserving active flux for 2D Euler flows"],"cache_read_input_tokens":2816,"weakest_assumption_plain":"Everything rests on the bicharacteristic-based evolution operator being a sufficiently accurate proxy for the true point-value evolution in the fully discrete scheme; if that operator is not accurate enough, the claimed third-order rate and the positivity guarantees in practice would degrade.","fun_headline_variants_meta":{"raw":{"variants":["2D Euler active flux with guaranteed positive states","Active flux keeps density and pressure positive in 2D","Third-order active flux for 2D Euler with positivity limiter","Compact active flux avoids negative pressure for 2D Euler","Positivity-preserving active flux for 2D Euler flows"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000248,"raw_usage":{"total_tokens":1309,"prompt_tokens":596,"completion_tokens":713,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":340,"completion_tokens_details":{"reasoning_tokens":631}},"tokens_in":340,"tokens_out":713,"duration_ms":6527,"temperature":1.0,"reasoning_tokens":631,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T22:47:43.268452+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a standard 2D Riemann problem or double Mach reflection on a fine uniform mesh and check whether any cell ever attains negative density or pressure before the final time; a single negative cell disproves the positivity guarantee. Alternatively, compute the convergence order on a smooth isentropic vortex: an observed order substantially below three would contradict the third-order claim.","supporting_citations":[],"review_version":1}