{"id":"582d46b4-817f-43a8-bdd1-49eae7a66936","arxiv_id":"2506.04857","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A positivity-preserving third-order active flux scheme for ideal MHD with the Godunov-Powell source term is constructed and validated.","lead":"A new numerical method for magnetized plasma flows keeps density and pressure positive while preserving third-order accuracy. It is a candidate building block for astrophysical and fusion simulation codes that need robust divergence control.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Positivity-preserving proof in §3.2 depends on an unstated geometric lemma: the box from coordinate-wise minima of shrunken vertices is asserted to lie in the admissible set, with no proof or fallback if false.","rationale":"The reader's weakest assumption identifies the same geometric step in Section 3.2 that we consider most load-bearing. A careful reading shows the PP proof for the cell average reduces to the unproved claim that the axis-aligned box [0,Lambda] from coordinate-wise minima of the shrunken vertices lies in the convex pressure-admissible set S_p. Without this, the blending coefficients theta from step (3) are not certified to keep U^{Lim,2} in G. The numerical tests provide empirical evidence but cannot cover all configurations, especially near strong shocks where the anti-diffusive flux differences are large and the scaling factors r_k can differ strongly between vertices. We do not see a mathematical contradiction, and the special structure of the vertices (axis and mixed supports) likely makes the lemma true, but the omission is a genuine gap in a central claim. Hence the conditional verdict stands; a proof or a targeted counterexample search would settle it.","tokens_in":24519,"tokens_out":34080,"duration_ms":352231,"concrete_test":"Write a standalone script that, for 10^4 random admissible states U^n and random flux-difference vectors H_L, H_R, H_D, H_U, performs the Section 3.2 construction: form the 16 vertices A_k, compute the scaling factors r_k by the first positive root of p(r A_k)=epsilon_p, build B_k, compute Lambda_I, and then (a) evaluate p at all 16 corners of the box [0,Lambda_L] x [0,Lambda_R] x [0,Lambda_D] x [0,Lambda_U], and (b) solve a linear program to test whether the whole box is contained in conv{B_k, 0}. If any corner has p < epsilon_p or the LP reports non-containment, the geometric lemma is false and the PP claim fails; if all pass, the step is empirically supported but still requires an analytical proof.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Section 3.2, after scaling each vertex of the density-admissible rectangle R_rho along the ray from the origin to obtain B_k (each in S_p), the paper sets Lambda_I = min_{k: k_I=1} B_{k,I} and asserts that this gives 'a rectangular subset inside the convex polygon with vertices B_k'. This is load-bearing: the final blending coefficients theta are taken as minima of these Lambda_I across adjacent cells, and the PP proof for the cell average update requires that every theta in [0,Lambda_L] x [0,Lambda_R] x [0,Lambda_D] x [0,Lambda_U] yields U^{Lim,2}(theta) in S_p. The coordinate-wise-minimum box is not generally contained in the convex hull of an arbitrary point set (e.g., in 2D, points (1,100) and (100,1) have m=(1,1) but (1,0) is outside conv{0,(1,100),(100,1)}). Here the special structure—axis vertices of the form (0,...,r_I,0,...) and mixed vertices—may rescue the claim, but no proof is given. If the box is not contained in S_p, the computed Lambda_I overestimate the admissible coefficients and the limited cell average can leave G even though all numerical tests pass. The round-off fallback in Remark 4.1 does not address this structural gap.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proposes a positivity-preserving active flux (AF) scheme for the two-dimensional ideal MHD equations with the Godunov-Powell source term. The cell average is updated by a conservative finite volume formula, while the nonconservative source term is discretized by a 3x3 Gauss-Lobatto quadrature based on the in-cell bi-parabolic reconstruction; the point values are evolved by the LLF flux vector splitting from [25] with a central-difference discretization of the divergence in the source term. Positivity is enforced by blending the high-order AF update with the first-order LLF scheme whose PP property was proved in [58]: a parametrized flux limiter is constructed for the cell average and a scaling limiter for the point values. A new shock sensor involving total pressure, velocity divergence and magnetic-field divergence is used to damp oscillations through additional blending. Numerical experiments (smooth accuracy tests, Orszag-Tang, rotor, blast, shock-cloud interaction, and MHD jets) demonstrate third-order convergence, positivity, oscillation control, and control of the discrete divergence error.","tokens_in":24838,"tokens_out":16195,"duration_ms":173536,"significance":"If the central claims hold, the paper is a valuable contribution: it extends the compact active flux methodology to ideal MHD with a positivity guarantee and with a Godunov-Powell source discretization that does not enlarge the stencil. The use of the externally proved first-order LLF PP property [58] as the fallback scheme is methodologically sound, and the numerical test suite is substantial, including problems with very low pressure that would crash without the limiters. The evidence for third-order accuracy and for divergence control is convincing. The main weakness is that two load-bearing PP assertions are not proved in the manuscript: the geometric containment step in Section 3.2 and the direct applicability of Lemma 3.1 to the staggered face-centered point-value LLF updates in Section 3.3. Because the paper's headline property is 'positivity-preserving', these gaps must be closed before the claims can be taken as established. No code or data is provided, so numerical reproducibility cannot be independently assessed.","major_comments":[{"comment":"In Step (2) of the parametrized flux limiter, the paper asserts that the rectangle defined by Lambda_I = min_{k_I=1} B_{k,I} is 'a rectangular subset inside the convex polygon with vertices B_k'. This containment is load-bearing: the PP proof of U^{Lim,2}_{i,j}(theta) requires every theta in the box [0,Lambda_L] x [0,Lambda_R] x [0,Lambda_D] x [0,Lambda_U] to yield an admissible state. For arbitrary point sets the coordinate-wise minimum box is not contained in the convex hull (e.g., conv{(0,0),(1,100),(100,1)} contains neither (1,0) nor (0,1)); the special ray-scaling structure may rescue the claim, but no argument is given. The round-off fallback in Remark 4.1 does not address a structural overestimate of the admissible coefficients. Please add a proof (or a precise reference) that the constructed box lies in S_p, or modify the construction so that containment is explicit.","section":"Section 3.2"},{"comment":"The statement that the three first-order LLF schemes for the point values 'are PP according to Lemma 3.1' is not justified for the face-centered updates. For U_{i+1/2,j}, the y-direction stencil uses points at y_{j-1/2}, y_j and y_{j+1/2} with spacing Delta_y/2, while the source-term divergence uses a denominator 2 Delta_y; Lemma 3.1 is stated for cell averages on a uniform mesh with spacing Delta_x, Delta_y and a central source term evaluated at distance Delta_y from the cell center. A direct application is therefore not immediate, and the point-value PP property is part of Definition 3.1. Please provide a proof (for example, by representing the update as a convex combination of Lemma 3.1 updates on suitable sub-grids) and state the corresponding CFL condition, or restrict the PP claim to the cell average with a separate numerical justification for the point values.","section":"Section 3.3"},{"comment":"The formula theta_{i,j} = min(1, (p(U^{LLF}_{i,j}) - eps_p)/(p(U^{LLF}_{i,j}) - p(U^{src}_{i,j}))) is not valid when p(U^{src}_{i,j}) >= p(U^{LLF}_{i,j}): the fraction is then negative and the coefficient leaves [0,1]. In that case concavity of the pressure gives p(U^{Lim,1}_{i,j}) >= eps_p for every theta in [0,1], so the standard fix is to set theta = 1 when the denominator is non-positive. Please add this case distinction; as written, the source-blending step of the PP construction is incomplete.","section":"Section 3.2, Step (1)"}],"minor_comments":[{"comment":"The numerical tests never state the CFL number or the time-step formula used. Since the PP property depends on the CFL condition in Lemma 3.1 and on analogous conditions for the point-value updates, please report the time-step control used in all examples.","section":"Section 4 / Section 2"},{"comment":"The discrete divergence values in equation (6) are labeled with indices (0,0), (0,1), ..., (2,2), while the quadrature points are labeled with l,m = 1,2,3. Please make the correspondence explicit to avoid confusion.","section":"Section 2.2"},{"comment":"The initial data for the Leblanc problem are listed as identical to those of the first Riemann problem in Example 4.1, with only the adiabatic index and final time changed. Please clarify whether this is intended and cite the specific Leblanc setup from [40] explicitly.","section":"Example 4.2"},{"comment":"The shock sensor strength kappa is chosen differently in nearly every test (kappa = 1, 2, 10, 50, 1000) and no sensitivity study is provided. Please comment on the sensitivity of the results to kappa and on how it should be selected in practice.","section":"Section 3.4"},{"comment":"In the formula for theta*_sigma, the denominator is written as rho(U^L_sigma) - rho(U^H_sigma); presumably U^{LLF}_sigma is intended. Please correct this typo.","section":"Section 3.3"},{"comment":"The phrase 'tapper function' should read 'taper function'.","section":"Example 4.5"}],"recommendation":"major_revision","confidential_remarks":"The main correctness risk is the unproved geometric containment lemma in Section 3.2. If the authors can supply a proof (or replace the construction with one whose containment is explicit), the paper would be a solid contribution to high-order MHD schemes. The face-centered point-value PP reduction in Section 3.3 also needs a rigorous statement. These are fixable within the scope of the manuscript, so major revision rather than rejection seems appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a real step forward for active flux methods—the first positivity-preserving AF scheme for ideal MHD with the Godunov-Powell source term—and the numerical evidence is convincing. The PP proof for the cell average has a gap, but it looks fixable.\n\nWhat's new: the source-term discretizations (Gauss-Lobatto for the cell average, central difference for point values) keep the compact stencil, and the divergence-aware shock sensor is a sensible addition. The scheme hits third order on smooth tests, positivity is demonstrated on the hard cases (Leblanc, low-pressure vortex, blast, jets), and the rotor comparison shows the source term actually controls divergence error. The first-order LLF PP property is imported from Wu–Shu, so the foundation is solid.\n\nThe soft spot is Section 3.2. After shrinking the rectangle vertices to make them pressure-admissible, the authors define Lambda_I as the coordinate-wise minimum of the shrunken vertices and assert this gives a rectangle inside the convex polygon. For arbitrary points that is false, and no proof is provided for this special structure. That is load-bearing: if the rectangle leaks out of S_p, the limited cell average is not guaranteed to be PP. The tests pass, so the claim may be true in this setting, but a referee should ask for a proof or a different limiter construction. Remark 4.1 also qualifies the strict PP claim—the 2^m 10^-8 shrinkage is an explicit fallback, minor but worth stating. No code or data is shipped, which is a modest reproducibility ding for a numerical methods paper.\n\nOverall the construction looks right and the experiments back it up. This deserves a serious referee: novel, clearly written, and the gap is fixable. I would accept after major revision.","headline":"First PP active flux scheme for ideal MHD with Godunov-Powell source term; numerics strong, but the cell-average PP proof has an unproved geometric step.","tokens_in":25346,"tokens_out":3656,"would_cite":true,"duration_ms":41365,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65M08","65M12","65M20","35L65"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that a compact, Riemann-solver-free active flux scheme can solve two-dimensional ideal magnetohrodynamics with third-order accuracy while provably keeping density and pressure positive.","keywords":["active flux method","ideal magnetohydrodynamics","positivity-preserving scheme","Godunov-Powell source term","divergence-free constraint","local Lax-Friedrichs flux vector splitting","shock sensor","third-order accuracy"],"falsifier":"A direct check would be to construct or randomly search for a set of admissible vertex states for which one vertex pressure is nonnegative only after radial shrinking, yet the coordinate-wise minimum rectangle is not contained in the pressure-admissible set; evaluating the limited scheme on that state should yield a negative-pressure cell average, contradicting the claimed guarantee. Short of that, a formal counterexample to the containment statement would settle the question, since the paper gives no proof of it.","tokens_in":24315,"feed_emoji":"🧲","tokens_out":8566,"duration_ms":72179,"temperature":0.7,"pith_summary":"This paper sets out to establish that the active flux method, a compact high-order finite-volume scheme that evolves cell averages and interface point values together, can be made safe for the two-dimensional ideal magnetohydrodynamics equations. The targeted properties are third-order accuracy, positivity of density and pressure, control of the magnetic divergence, and suppression of oscillations near strong shocks. The proposed scheme adds the Godunov-Powell source term to the conservative equations and discretizes it on the same compact stencil, so divergence errors are advected away rather than accumulated. To keep the solution admissible, the high-order update is blended with a first-order local Lax-Friedrichs scheme that is provably positivity-preserving; the cell-average blend uses a parametrized flux limiter and the point-value blend uses a scaling limiter. A reader should care because high-order MHD solvers frequently fail on near-vacuum, strongly magnetized flows, and the paper demonstrates the scheme on rotor, blast-wave, shock-cloud, and high-Mach jet problems where the un-limited scheme fails.","feed_headline":"Third-order MHD scheme keeps density and pressure positive","feed_subtitle":"By blending with a positivity-preserving baseline, the scheme survives strong shocks and near-vacuum jets.","key_machinery":"The central mechanism is the active flux semi-discretization itself: cell averages are advanced by a conservative finite-volume flux through Simpson's rule, while interface point values are advanced by local Lax-Friedrichs flux vector splitting with upwind one-sided differences; because point values are continuous across interfaces, no Riemann solver is needed. The Godunov-Powell source term $-\\left(\\nabla\\cdot B\\right)\\Psi$, with $\\Psi=(0,B_1,B_2,v\\cdot B)^\\top$, is discretized compactly: the cell-average source uses the $3\\times3$ Gauss-Lobatto quadrature on the bi-parabolic reconstruction, and the point-value divergence is a central difference built on the same stencil. Positivity is carried by two limiters: a parametrized flux limiter that blends the high-order cell-average fluxes and source with the first-order positivity-preserving LLF scheme, and a scaling limiter for point values. The shock-suppression machinery is a sensor $\\theta^s_{i+1/2,j}=\\exp\\left(-\\kappa[(\\varphi_1)(\\varphi_2)+\\varphi_3]\\right)$ built from total-pressure jumps, velocity dilatation, and a new discrete $\\nabla\\cdot B$ term, which sets local blending coefficients for both fluxes and the source term.","core_discovery":"The central claim is that the active flux scheme with the LLF flux vector splitting and the Godunov-Powell source term is third-order accurate and positivity-preserving for both cell averages and point values, while keeping the discrete divergence under control. For the cell average, the source term is discretized by a $3\\times3$ Gauss-Lobatto quadrature using the bi-parabolic reconstruction, and the flux is blended with the first-order positivity-preserving LLF flux through a parametrized flux limiter; for the point values, the flux derivative uses LLF splitting and the divergence is approximated by a central difference, with a scaling limiter enforcing admissibility. The paper further introduces a shock sensor that combines a total-pressure jump sensor, a velocity-dilatation sensor, and a new discrete magnetic-divergence term; the sensor's blending coefficient damps the high-order flux and source term near discontinuities. Numerical tests on two-dimensional problems, including smooth accuracy tests, the Orszag-Tang vortex, rotor, blast wave, shock-cloud interaction, and magnetized jets, are presented as evidence that the scheme retains third-order convergence, positivity, and stable divergence error over long times.","pith_inferences":["If the unproved rectangle-containment assertion fails, the scheme may still be numerically positivity-preserving; the paper's guarantee would weaken from a theorem to an empirical property, and a repairing argument would need to find a smaller or differently shaped admissible subset.","The same two-step blending (source first, then flux) could plausibly transfer to three-dimensional MHD, where the additional direction enlarges the parameter rectangle but the convexity argument appears to require only more bookkeeping.","The new magnetic-divergence term in the shock sensor could be reused as a standalone refinement indicator: regions where $\\varphi_3$ is large mark places where the discrete divergence has not yet been advected away, which is exactly where a divergence-cleaning or mesh-refinement step would pay off.","Because the cell-average blending is designed to keep mass conservation through unique interface fluxes, the method may be a useful building block for MHD codes that need strict conservation of mass but can tolerate the nonconservative momentum and energy updates introduced by the Godunov-Powell term."],"forward_implications":["The scheme attains third-order $\\ell^1$ convergence on smooth two-dimensional tests, including a vortex whose central pressure drops to about $5.3\\times10^{-12}$.","With the PP limitings active, the scheme runs to completion on the Leblanc shock, rotor, blast-wave, shock-cloud, and high-Mach jet problems; without them, negative pressure appears within one step or shortly after.","The Godunov-Powell source term, discretized as specified, controls divergence error: the two monitored discrete divergences grow slowly and reach a plateau in the Orszag-Tang test, and the rotor test shows distorted Mach contours when the source term is omitted from either the cell-average or the point-value update.","The shock-sensor blending is applied locally rather than globally, so high-order resolution is retained in smooth regions while oscillations near strong discontinuities are damped.","The scheme stays compact: the source-term discretizations for both cell averages and point values use only degrees of freedom in the current cell or on its edges, matching the stencil of the original active flux method."],"supporting_citations":[{"why":"Introduces the local Lax-Friedrichs flux vector splitting for the active flux point-value update and the idea of blending the high-order AF scheme with a first-order PP scheme, which this paper adapts to MHD with the Godunov-Powell source term.","marker":"[25]"},{"why":"Proves that the first-order LLF scheme with a suitable discrete Godunov-Powell source term is positivity-preserving on general meshes, the baseline whose PP property the present limiters exploit.","marker":"[58]"},{"why":"Defines the semi-discrete active flux method with bi-parabolic reconstruction and compact stencil that the MHD scheme extends.","marker":"[3]"},{"why":"Supplies the parametrized flux limiter technique for positivity preservation of high-order finite difference schemes, adapted here to the cell-average update with a source term.","marker":"[20]"},{"why":"Extends the parametrized flux limiter to the ideal MHD equations, informing the two-step limiting procedure for cell averages.","marker":"[19]"},{"why":"Shows the provably positive discontinuous Galerkin method for multidimensional ideal MHD using the Godunov-Powell source term, the closest predecessor and the source of several test problems.","marker":"[57]"},{"why":"Introduces the eight-wave formulation and shows the divergence error satisfies a transport equation, the conceptual basis for using the Godunov-Powell source term to control divergence.","marker":"[46]"},{"why":"Original parametrized maximum-principle-preserving flux limiter for high-order schemes, from which the flux-limited PP cell-average update is derived.","marker":"[60]"}],"fun_headline_variants":["Positivity-preserving Active Flux for ideal MHD","Third-order MHD AF scheme survives shocks and near-vacuum","Godunov-Powell source ensures AF MHD positivity","AF MHD scheme with shock sensor keeps density and pressure positive","Third-order MHD scheme with source term controls divergence"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"In the cell-average limiter of Section 3.2, the pressure-positivity guarantee rests on an unproved geometric assertion: after radially shrinking each vertex of the parameter rectangle to make it pressure-admissible, the rectangle formed by the coordinate-wise minima of the shrunk vertices still lies inside the pressure-admissible set; if that containment fails, the blended cell average is not guaranteed to stay positive even though the numerical tests do.","fun_headline_variants_meta":{"raw":{"variants":["Positivity-preserving Active Flux for ideal MHD","Third-order MHD AF scheme survives shocks and near-vacuum","Godunov-Powell source ensures AF MHD positivity","AF MHD scheme with shock sensor keeps density and pressure positive","Third-order MHD scheme with source term controls divergence"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000723,"raw_usage":{"total_tokens":3311,"prompt_tokens":1078,"completion_tokens":2233,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":694,"completion_tokens_details":{"reasoning_tokens":2151}},"tokens_in":694,"tokens_out":2233,"duration_ms":21638,"temperature":1.0,"reasoning_tokens":2151,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T10:32:16.197256+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct check would be to construct or randomly search for a set of admissible vertex states for which one vertex pressure is nonnegative only after radial shrinking, yet the coordinate-wise minimum rectangle is not contained in the pressure-admissible set; evaluating the limited scheme on that state should yield a negative-pressure cell average, contradicting the claimed guarantee. Short of that, a formal counterexample to the containment statement would settle the question, since the paper gives no proof of it.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the local Lax-Friedrichs flux vector splitting for the active flux point-value update and the idea of blending the high-order AF scheme with a first-order PP scheme, which this paper adapts to MHD with the Godunov-Powell source term."},{"cited_title":"Wu and C.-W","cited_arxiv_id":null,"evidence_quote":"Proves that the first-order LLF scheme with a suitable discrete Godunov-Powell source term is positivity-preserving on general meshes, the baseline whose PP property the present limiters exploit."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the parametrized flux limiter technique for positivity preservation of high-order finite difference schemes, adapted here to the cell-average update with a source term."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Extends the parametrized flux limiter to the ideal MHD equations, informing the two-step limiting procedure for cell averages."},{"cited_title":"Wu and C.-w","cited_arxiv_id":null,"evidence_quote":"Shows the provably positive discontinuous Galerkin method for multidimensional ideal MHD using the Godunov-Powell source term, the closest predecessor and the source of several test problems."},{"cited_title":"Report 94-24, ICASE, NASA Langley, 1994","cited_arxiv_id":null,"evidence_quote":"Introduces the eight-wave formulation and shows the divergence error satisfies a transport equation, the conceptual basis for using the Godunov-Powell source term to control divergence."},{"cited_title":"Xu , Parametrized maximum principle preserving flux limiters for high order schemes solving hyperbolic conservation laws: O ne-dimensional scalar problem , Math","cited_arxiv_id":null,"evidence_quote":"Original parametrized maximum-principle-preserving flux limiter for high-order schemes, from which the flux-limited PP cell-average update is derived."}],"review_version":1}