{"id":"323f9ac5-ee29-4283-9012-4352a4050063","arxiv_id":"2608.07345","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A projection-prolongation mesh-based filter lets explicit RKDG methods on spherical-polar meshes run with the larger time steps of an auxiliary merged mesh, with a 1D equivalence proof and 2D/3D demonstrations.","lead":"This paper introduces a mesh-based filter for discontinuous Galerkin simulations in spherical coordinates that lets simulators take much larger time steps near the origin and poles. The authors prove the filter exactly matches a coarser-mesh calculation in one dimension and show large speedups in 2D and 3D tests.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 1D equivalence theorem is sound, but the production 2D/3D filter uses a metric-weighted projection and a slope-limiter-before-filter ordering that the authors concede falls outside Theorem 1; stability inheritance for the tested method is therefore unproven and should remain conditional.","rationale":"I read the paper in good faith and found the 1D theoretical contribution internally consistent: Lemma 1 is a short moment argument, and Theorem 1 telescopes the fine-cell residuals correctly because the numerical flux is single-valued and the merged basis function is continuous across internal fine-cell interfaces. The numerical results in Tables 1 and 2 independently support the theorem, since the filtered errors match the coarser-mesh errors exactly as expected. My concern is not with the theorem but with how far the paper's language extends it. The manuscript itself flags the decisive limitation in Section 4: the slope limiter is applied on the fine mesh before filtering, an ordering that the authors explicitly say falls outside Theorem 1 and for which TVD stability is not guaranteed. In addition, the 2D/3D projection uses metric-weighted inner products and a moment-matched merged metric, so even the unlimited algebraic equivalence is not proved for the production algorithm. These departures matter because the paper's broadest claim is that the filtered method inherits the stability and accuracy of the merged-mesh RKDG discretization; the reported 2D/3D tests are encouraging but do not establish that claim for general problems with stronger shocks or different flow orientations. I therefore maintain the reader's CONDITIONAL verdict rather than raising it to reject: the 1D result is solid and the 2D/3D evidence is promising, but the stability-inheritance assertion for the practical method remains unproven and should be stated conditionally or supported by an explicit equivalence or stability analysis for the implemented algorithm.","tokens_in":28252,"tokens_out":12888,"duration_ms":133191,"concrete_test":"Run a 2D/3D equivalence check with the production projection from Equations (40) and (48), k=1, no limiter: for a small periodic problem or a single merged cell, initialize the fine solution as the prolongation of a merged-mesh polynomial, take one SSPRK3 step with the merged-mesh time step on both the filtered fine mesh and the merged mesh, and compare the projected fine update with the merged-mesh update to machine precision. Then repeat the same check with the Section 4 ordering, applying the slope limiter on the fine mesh before filtering. If the first comparison fails, the metric-weighted projection breaks the equivalence claimed by Theorem 1; if only the second fails, the limiter ordering is the operative gap. Either outcome would show that the stability-inheritance claim must be restricted or the production implementation changed to match the theorem's hypotheses.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central theoretical result is Theorem 1 in Section 2.3.4, proved for the unweighted L2 projection-prolongation in Equation (16) with no fine-mesh limiter and under exact quadrature. The abstract's claim that the filtered method inherits accuracy and stability is stated as a consequence of this equivalence. The 2D/3D implementation that produces the headline speedups departs from these hypotheses in two concrete ways. First, the projection is metric-weighted: Equations (38)-(40) replace (16) with sqrt(gamma) on the fine side and a moment-matched sqrt(hat-gamma) on the merged side, and the Euler equations include geometric source terms in (28); no theorem shows that the residual-sum argument of Theorem 1 survives this replacement. Second, Section 4 says the slope limiter is applied on T before the mesh-based filter and explicitly states this ordering 'falls outside the conditions for the equivalence established in Theorem 1; consequently, stability properties of the RKDG scheme on the merged mesh, such as the TVD property, are not guaranteed to carry over.' Thus the substantial time-step gains in Figures 10 and 12 and the good agreement in Figures 8, 11, and 12 are empirical observations for the selected tests, not consequences of the proved equivalence. The 1D result itself appears correct and is independently supported by the convergence checks in Tables 1 and 2; the gap is the unproven bridge from that theorem to the practical algorithm.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a mesh-based filtering procedure for explicit Runge-Kutta discontinuous Galerkin (RKDG) methods on spherical-polar meshes. The method constructs an auxiliary merged mesh whose elements combine fine cells so that extreme angular-to-radial aspect ratios are removed, then evolves the solution on the original fine mesh using the larger time step allowed by the merged mesh, applying a projection-prolongation filter after each Runge-Kutta stage. In one dimension, Theorem 1 proves that, under exact quadrature and without limiters, one Forward Euler step on the fine mesh followed by projection equals one Forward Euler step on the merged mesh; Remark 3 extends this stage-by-stage to explicit Runge-Kutta methods. Numerical tests for linear transport and Burgers' equation confirm the expected convergence orders and show that the unfiltered fine-mesh update with the large time step blows up. In two and three dimensions, the authors introduce angular merging criteria (Eqs. 37 and 45) and a metric-weighted projection using sqrt(gamma) (Eqs. 38-40 and 46-48), and they demonstrate significantly larger stable time steps on a two-dimensional Riemann problem, a shock-bubble interaction, and a three-dimensional Sedov-Taylor blast wave. Section 4 explicitly discloses that the slope limiter is applied on the fine mesh before the filter, which the authors state falls outside the conditions of Theorem 1.","tokens_in":28547,"tokens_out":6660,"duration_ms":63998,"significance":"The one-dimensional equivalence theorem is a genuine and clearly presented contribution: it shows that projection-prolongation with a single-valued numerical flux reproduces exactly the merged-mesh DG update, so the filtered method inherits stability and accuracy in that idealized setting. The numerical validation is clean, including the blow-up of the unfiltered large-time-step update, and the paper is commendably explicit about the scope of its theorem. However, the practical two- and three-dimensional algorithm departs from the theorem in two load-bearing ways: the metric-weighted projection is not covered by the 1D analysis, and the production code applies the slope limiter on the fine mesh before filtering, an ordering the authors concede is outside the equivalence. Consequently, the abstract's and Section 5's claims that the filtered method 'inherits the accuracy and stability properties' are stronger than what the paper establishes for the tested 2D/3D method. The 2D/3D speedups in Figures 10 and 12 are valuable empirical results, but the mathematical bridge from Theorem 1 to those results remains incomplete.","major_comments":[{"comment":"The production implementation applies the slope limiter on T before the mesh-based filter, and the text states that this ordering 'falls outside the conditions for the equivalence established in Theorem 1; consequently, stability properties of the RKDG scheme on the merged mesh, such as the TVD property, are not guaranteed to carry over.' Because every 2D/3D result in Figures 8-12 is produced with this limiter ordering, the stability and robustness claims for the practical method are empirical only. Please either prove a stability statement for this ordering, modify the implementation so that the limiter is applied on the merged mesh before prolongation to T, or explicitly restrict the inheritance claim to the cases covered by Theorem 1 and Remark 5.","section":"§4, first paragraph"},{"comment":"The metric-weighted projection is not covered by the one-dimensional analysis, and there is a concrete obstruction to the natural extension. Lemma 1 relies on the fact that the fine projection P annihilates u-Pu against every polynomial of degree at most k on a fine cell. In Eq. (39), sqrt(hat-gamma) is a degree-k polynomial on the merged cell, so for a merged test function phi of degree k, the integrand phi sqrt(hat-gamma) restricted to a fine cell has degree up to 2k. Consequently, bP[P[u]] = bP[u] is not guaranteed, and the induction base in Remark 2 fails unless additional conditions are imposed. For the same reason, the residual-sum argument in Theorem 1 does not automatically extend to the flux terms in Eq. (40), which involve sqrt(gamma) and sqrt(hat-gamma) on the two sides. The 2D/3D equivalence should either be proved under explicit hypotheses on the metric representation, or the claims should be weakened to empirical observations.","section":"§3.1.2, Eqs. (38)-(40)"},{"comment":"The momentum equation in the Euler system contains a geometric source term (1/2) Pi^{ik} partial_j gamma_{ik}, and Remark 4 asserts that source terms can be included 'provided it is evaluated with the same quadrature and representation on the fine and merged meshes.' In the implementation, the fine-mesh and merged-mesh metric representations differ (sqrt(gamma) versus sqrt(hat-gamma) defined by moment matching), so the source-term evaluation on the fine mesh is not identical to the merged-mesh source-term evaluation. The equivalence for the momentum equation is therefore not established, and the conservation and accuracy statements for the 3D Euler tests should be presented as numerical evidence rather than as consequences of Theorem 1.","section":"§3.2.2 and Eq. (28b)"}],"minor_comments":[{"comment":"The left-panel labels appear to list 'I5' twice in the fine-mesh sequence, which is likely a typo for I6.","section":"Figure 1 caption"},{"comment":"The notation x_{p-1/2} and x_{p+1/2} for the endpoints of a merged cell is used without definition; defining these endpoints would improve readability.","section":"§2.3.4, proof of Theorem 1"},{"comment":"The notation 'sin hat-theta' is ambiguous because hat-theta denotes a merged cell, not a specific angle; clarify how sin theta is evaluated on a merged cell, for example by using the cell midpoint or a quadrature rule.","section":"§3.2.1, Eq. (45)"},{"comment":"The sentence 'For the first order method, we include results when the fine-mesh solution is advanced with cDelta-t without applying the mesh-based filter' is followed by blow-up values in Table 1; stating explicitly that these values are included to demonstrate instability would help the reader.","section":"§2.4.1"},{"comment":"The abstract contains the typo 'ans-stage' instead of 'an s-stage'; this should be corrected in the final version.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The one-dimensional theorem is sound and the numerical experiments are convincing, but the paper's broad claims about inheriting stability in 2D/3D are not yet supported by analysis. The metric-weighted projection gap is more than a presentation issue: the natural extension of Lemma 1 fails without additional hypotheses, so the paper needs either a new theorem covering the weighted case or a clear restriction of the theoretical claims. The empirical results are valuable even if the theoretical scope is narrowed, so I would encourage a revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing to know about this paper is that the core theoretical contribution holds up, and it is genuinely new: a projection-prolongation filter that lets an RKDG scheme on a fine mesh take time steps dictated by a coarser merged mesh, with a proof in 1D that the filtered fine-mesh update is exactly the merged-mesh DG update. Theorem 1 is correct as far as I can tell; the cancellation of internal fluxes works because the numerical flux is single-valued and the test function is continuous across the merged cell, and Lemma 1 sets up the induction. The 1D linear transport and Burgers tests match the theorem: filtered solutions with N cells match unfiltered solutions with N/2 cells, and the no-filter blow-up with the large step is a nice sanity check.\n\nCredit also for the 2D/3D merging criteria, which sensibly merge along theta and phi so that proper cell lengths stay comparable to Delta r, and for the metric-weighted projection, which they show conserves mass and energy to ~1e-13 in the shock-bubble test. The reported speedups are striking: ~18x fewer steps on the 2D Riemann at Ntheta=64, ~13x on the 256x256 shock-bubble, and ~100x on the 64x8x8 Sedov.\n\nThe soft spots are real but they are in the bridge, not the 1D theorem. First, the production 2D/3D filter uses a metric-weighted L2 projection (Eqs. 38-40) and a slope-limiter-before-filter ordering; the authors explicitly state in Section 4 that this ordering 'falls outside the conditions for the equivalence established in Theorem 1,' so the TVD property is not inherited. The abstract's sentence about inheriting accuracy and stability is too strong as written: it is proved only in 1D under exact quadrature, and the 2D/3D claims are empirical. Second, the 400x speedup for 64x16x16 is extrapolated, not computed; code and data are not provided, so reproducibility currently rests on the description. Third, the 1D analysis assumes exact quadrature, while the implementation uses Gauss quadrature; the gap is likely small but it is not discussed.\n\nWho should read this: anyone working on DG in spherical-polar coordinates or on geometric CFL alleviation in astrophysics. The 1D theorem and filter construction are solid contributions. The multi-D claims need either a proof or a clear empirical validation with relaxed claims. A serious referee should be assigned; the gaps are addressable and the paper would benefit from revision.","headline":"Sound 1D equivalence theorem for a DG cell-merging filter; the 2D/3D stability claims are empirical and should be labeled as such.","tokens_in":29091,"tokens_out":1971,"would_cite":true,"duration_ms":17548,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65M60","65M12","35L65"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper proves that a projection–prolongation filter applied after every Runge–Kutta stage makes a fine-mesh discontinuous Galerkin update exactly equivalent to evolving on a coarser merged mesh, letting simulations take the merged…","keywords":["discontinuous Galerkin methods","Runge–Kutta time integration","spherical-polar coordinates","CFL time-step restriction","mesh merging","projection-prolongation filter","Euler equations","cell anisotropy"],"falsifier":"Take the one-dimensional Burgers equation with quadratic DG and SSPRK3, and at each stage apply a slope limiter on the fine mesh before the projection–prolongation step; compare the resulting nodal values with a run that applies the same limiter on the merged mesh before prolongation. Any difference shows the practical limiter ordering used in two and three dimensions breaks the Theorem 1 equivalence, pinpointing the mechanism behind the reported shock overshoot.","tokens_in":28034,"feed_emoji":"⚡","tokens_out":9963,"duration_ms":77585,"temperature":0.7,"pith_summary":"The paper proposes a post-processing filter for explicit Runge–Kutta discontinuous Galerkin (RKDG) methods on spherical-polar meshes, where converging coordinate lines create extremely small cells near the origin and poles and force tiny stable time steps. The filter projects the fine-mesh solution onto an auxiliary \"merged mesh\" with larger cells, evolves on the fine mesh using the merged mesh's larger time step, and prolongs back after each RK stage. In one dimension the paper proves that this filtered update is exactly equivalent to evolving the RKDG discretization directly on the merged mesh, so the filtered method inherits the merged mesh's accuracy and stability. In two and three dimensions the same idea is implemented with metric-weighted projections, and numerical tests on the Euler equations show accurate solutions with time steps independent of angular resolution and up to roughly 100 times larger (with a projected factor of about 400 in three dimensions).","feed_headline":"Filtering each Runge–Kutta stage matches merged-mesh evolution","feed_subtitle":"Removes the tiny-cell CFL bottleneck at spherical-polar origin and poles, with 3D speedups projected near 400x.","key_machinery":"The load-bearing object is the projection–prolongation filter: an $L^2$ projection of the fine-mesh DG solution onto the merged mesh, followed by restriction of the merged polynomial back onto the fine cells. In nodal form it is assembled from cell merging weights $c^{(m)}_{i,\\alpha} = \\hat\\ell^{(p)}_i(x^{(m)}_\\alpha)$, which express each merged-cell basis function as a linear combination of fine-cell basis functions, so the update on any fine cell is a weighted sum of nodal values from all fine cells in the merged cell. The proof of Theorem 1 uses the fact that the merged basis functions are valid test functions on the fine cells and that, summed over a merged cell, the single-valued numerical fluxes at internal interfaces cancel, leaving only the merged-cell boundary terms. In two and three dimensions the projection is weighted by the metric volume factor $\\sqrt{\\gamma}$, with a merged representation $\\sqrt{\\hat\\gamma}$ matched by moments, to preserve mass and energy.","core_discovery":"The central discovery is Theorem 1: if at time $t^n$ the fine-mesh DG solution is the prolongation of a merged-mesh solution, then one Forward Euler step on the fine mesh with the merged-mesh time step, followed by projection onto the merged mesh, produces exactly the same solution as one Forward Euler step on the merged mesh. Because the merged basis functions are valid test functions on the fine cells, and because single-valued numerical fluxes at internal fine-cell interfaces cancel when summed over a merged cell, the projected fine-mesh residual equals the merged-mesh residual. The argument extends stage by stage to any explicit Runge–Kutta method, provided the projection–prolongation is applied at each inner stage. Consequently the filtered fine-mesh method is equivalent to the RKDG method on the merged mesh and inherits its stability and accuracy properties; in particular, evolving with the merged-mesh time step is stable even though that time step exceeds the fine-mesh CFL limit.","pith_inferences":["The equivalence argument is coordinate-agnostic, since it only uses element-wise polynomial spaces and single-valued fluxes, so the same filter should transfer to other structured curvilinear meshes with geometric anisotropy, such as cylindrical grids, provided a metric-weighted projection is used.","A testable extension: apply the slope limiter on the merged mesh before prolongation, as Remark 5 suggests, and compare the resulting 1D shock profiles with the fine-mesh-limited ordering used in the paper's 2D/3D tests; this would show whether the observed shock overshoot is an artifact of the limiter ordering.","Because the filtered time step becomes set by radial resolution rather than angular resolution, the merging criterion can be calibrated problem-by-problem, and the relaxed-criterion results in the paper indicate how to measure that accuracy-versus-speed tradeoff.","An implicit–explicit extension is natural: run explicit stages with the merged-mesh time step while performing implicit solves on the merged mesh, which would reduce cost for stiff source terms such as neutrino kinetics in core-collapse supernova simulations."],"forward_implications":["The filtered RKDG method inherits the merged mesh's accuracy and stability in one dimension: with $N$ fine cells it reproduces the error of an $N/2$-cell merged computation while using a time step twice as large for pairwise merging.","Angular refinement no longer reduces the stable time step: in the two-dimensional Riemann test the number of time steps stays roughly constant as $N_\\theta$ increases from 16 to 64, and the filtered $128\\times64$ run uses about $3.6\\times10^3$ steps versus $6.3\\times10^4$ for the unfiltered mesh.","Conservation is preserved: relative changes in total mass and energy remain at $O(10^{-13})$ in the shock–bubble test, matching the unfiltered RKDG method.","In three dimensions the benefit is larger: filtered time steps are essentially independent of angular resolution, about 25 times larger on the $64\\times4\\times4$ mesh and 100 times on the $64\\times8\\times8$ mesh, with a projected factor of about 400 on the $64\\times16\\times16$ mesh.","The merging criterion can be relaxed to trade time-step size for angular resolution: the relaxed $128\\times64$ mesh uses more steps than the standard merged mesh but has error closer to the unfiltered solution."],"supporting_citations":[{"why":"Supplies the RKDG framework, the TVD stability property, and the CFL table used for the one-dimensional time steps.","marker":"B. Cockburn & C.-W. Shu 2001"},{"why":"Foundational RKDG discretization that the filter is inserted into.","marker":"B. Cockburn & C.-W. Shu 1998"},{"why":"Provides the SSPRK3 time integrator used in the numerical tests.","marker":"C.-W. Shu & S. Osher 1988"},{"why":"Finite-volume cell-merging approach that motivates the auxiliary merged-mesh idea.","marker":"R. Asaithambi & K. Mahesh 2017"},{"why":"Finite-volume mesh-coarsening strategy in spherical-polar coordinates that the filter generalizes to DG.","marker":"B. Müller 2015"},{"why":"Cell-merging technique in spherical-polar finite volume that inspired the merged-mesh construction.","marker":"K. Nakamura et al. 2019"},{"why":"Fourier-filtering alternative; its discussion of conservation and monotonicity in prolongation frames the DG filter design.","marker":"B. Müller et al. 2019"},{"why":"Provides the curvilinear-coordinate DG formulation and polynomial metric representation used by the two- and three-dimensional implementation.","marker":"D. Pochik et al. 2021"},{"why":"Describes the structured-mesh DG code base into which the filter is integrated for the numerical experiments.","marker":"E. Endeve et al. 2019"}],"fun_headline_variants":["Stage-wise filter equals merged mesh for polar RKDG","Polar RKDG: filter each stage, take bigger time steps","Mesh filter matches merged-mesh evolution in RKDG","Avoid CFL bottleneck: filtered stages mimic merged mesh","Spherical-polar RKDG: filter stages to lift timestep limits"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The practical two- and three-dimensional implementation applies a slope limiter on the fine mesh before the filter, an ordering the authors explicitly note falls outside the conditions of Theorem 1, so the stability observed in those tests depends on the unproven assumption that the one-dimensional equivalence extends to multiple dimensions and survives that limiting order.","fun_headline_variants_meta":{"raw":{"variants":["Stage-wise filter equals merged mesh for polar RKDG","Polar RKDG: filter each stage, take bigger time steps","Mesh filter matches merged-mesh evolution in RKDG","Avoid CFL bottleneck: filtered stages mimic merged mesh","Spherical-polar RKDG: filter stages to lift timestep limits"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000226,"raw_usage":{"total_tokens":1482,"prompt_tokens":974,"completion_tokens":508,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":590,"completion_tokens_details":{"reasoning_tokens":421}},"tokens_in":590,"tokens_out":508,"duration_ms":4858,"temperature":1.0,"reasoning_tokens":421,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T05:53:52.867674+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the one-dimensional Burgers equation with quadratic DG and SSPRK3, and at each stage apply a slope limiter on the fine mesh before the projection–prolongation step; compare the resulting nodal values with a run that applies the same limiter on the merged mesh before prolongation. Any difference shows the practical limiter ordering used in two and three dimensions breaks the Theorem 1 equivalence, pinpointing the mechanism behind the reported shock overshoot.","supporting_citations":[{"cited_title":"J., et al","cited_arxiv_id":null,"evidence_quote":"Describes the structured-mesh DG code base into which the filter is integrated for the numerical experiments."}],"review_version":1}