{"id":"fed970e2-d5a4-4973-a610-d00998dc75cd","arxiv_id":"2607.14403","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Upwind embedded-boundary SBP operators with interior accuracy up to order 9 and boundary order 4 are error-minimized and verified on a 3D wave equation with spherical excision, converging near the designed rates (~5.5 total L2).","lead":"A new family of upwind finite-difference operators lets wave equations run on simple Cartesian grids with holes cut out at arbitrary positions, reaching up to 9th-order interior and 4th-order boundary accuracy. The paper demonstrates design-order convergence on a 3D spherical-excision wave test, promising cheaper, simpler black-hole excision codes in numerical relativity.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"3D robustness/convergence claim rests on exact-solution injection in the boundary layer (Sec. IV.A); the boundary closure is not actually exercised, so stability for genuine IBVPs remains untested.","rationale":"The paper's central claim, as summarized in the abstract and conclusion, is that the newly derived embedded-boundary upwind SBP operators are robust and accurate for arbitrarily shaped domains, specifically in view of black-hole excision. The operator-existence part is supported by the construction procedure (ansatz, accuracy conditions, error minimization) and the promised supplementary coefficient tables. However, the 3D numerical evidence, which is the only evidence for the use-case claim, uses a manufactured-solution setup with exact-solution injection in the layer of active points inside the excision boundary and exact characteristic BCs at the boundary. This means the numerical scheme never has to handle the feedback between the embedded boundary and the interior evolution on its own. The acknowledged lack of a 2D/3D stability proof (Section III.B) then becomes consequential: without a genuine IBVP test, the paper does not demonstrate that the operators can be used for the intended applications. This is the most load-bearing weakness because it affects the interpretation of the headline convergence order and the robustness claim. Other concerns, such as the duplicated text in Section III.A.2 or the absence of a direct dispersion-relation comparison to prior embedded operators, are manuscript-quality issues or secondary claims that do not undermine the operator derivation itself. The reader's weakest_assumption identified the same issue: 3D stability is unproven and the test is solution-injected. I agree, and I recommend keeping the verdict CONDITIONAL: the operator-existence and formal-accuracy claims can stand, but the paper must either provide the promised coefficients (to verify the derivation) and add a more realistic stability/convergence test that does not rely on exact injection in the boundary layer. The proposed concrete test would directly settle whether the boundary closure is stable with characteristic data derived from the interior, rather than from the exact solution.","tokens_in":18036,"tokens_out":6213,"duration_ms":71052,"concrete_test":"Repeat the D±9−4 cubed-sphere excision run at the same five resolutions and t=500, but replace the exact evolution of the active interior-boundary points with one-sided extrapolation from the interior numerical solution, and inject incoming characteristics computed from the instantaneous interior solution rather than the exact reference. Add a small random perturbation to initial data. Monitor E_L2 and a discrete energy over t in [0,5000]. If any run shows persistent error growth or the convergence order drops from ~5.5 toward b+1 or lower, the robustness claim fails; if stability is retained and order ~5.5 persists, the concern is resolved.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The only 3D demonstration (Section IV.B) is not a genuine initial-boundary-value test. In Section IV.A, active grid points just inside the excision sphere are evolved using exact time derivatives of the reference solution, and both the inner and outer boundary conditions are injected as exact incoming characteristics. Thus the embedded-boundary closure is never required to be stable on its own: any unstable mode that would result from the coupling between Cartesian stencils and the arbitrarily placed boundary is suppressed because the layer of points feeding the boundary stencils is prescribed analytically. The claimed 'robustness and accuracy' for arbitrarily shaped domains and black-hole excision is therefore supported only by a manufactured-solution test with a filled-in ghost zone. This does not refute the operator-existence claim, but it removes the central evidence for the paper's headline use case. The paper itself concedes in Section III.B that 2D/3D stability with embedded boundaries cannot be proven and that more challenging geometries may require stronger dissipation. The convergence rates (b+1.5, b+1) may reflect the exterior stencils and exact data, not the actual behavior of the boundary closure with characteristic data derived from the interior solution.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper constructs new upwind summation-by-parts (SBP) finite-difference operators for embedded boundaries on Cartesian grids. The construction in Section III builds a parameterized ansatz for (Q+,H), imposes polynomial accuracy conditions (Eqs. 18-19), minimizes a boundary-truncation-error function f (Eq. 21), and verifies positive-definiteness of H and negative semi-definiteness of S. This yields minimal-width operators D±_{p-b} for interior orders p=2..9 with boundary order up to 4, plus dispersion-relation-preserving (DRP) variants. The operators are then tested in Section IV on the 3D curvilinear scalar wave equation in a multiblock cubed-sphere grid with an inner spherical excision, with convergence rates reported for the total L2, inner-block L2, and L∞ norms. The paper claims these are the highest-order embedded-boundary finite-differencing operators derived to date and argues that they are suitable for black-hole and worldtube excision on Cartesian grids.","tokens_in":18241,"tokens_out":8976,"duration_ms":90632,"significance":"If the operator construction and the accompanying claims are accepted, this is a useful and nontrivial contribution: it extends embedded-boundary SBP from the previous order range to interior order 9 and boundary order 4, introduces upwind and DRP variants, and provides explicit operators with positive norms and dissipative S. The supplementary materials with polynomial coefficients and a Mathematica notebook are a strength, as is the fact that the closure parameters are determined by an internal error-minimization criterion rather than fitted to convergence data. However, the numerical demonstration in Section IV does not, as designed, exercise the embedded-boundary closure as a genuine initial-boundary-value problem, because active points just inside the excision sphere are evolved with exact time derivatives of the reference solution. This weakens the paper's central robustness claim for realistic excision applications, although the operator-existence and accuracy claims appear defensible.","major_comments":[{"comment":"The text calls Eqs. (18)-(19) the '2(b+1) accuracy conditions', but the conditions are displayed only for q∈[1,b], which is 2b conditions. More importantly, starting at q=1 omits the zero-th order consistency condition (Q± + (1/2)B)X^0 = 0 (with X^{-1} understood as zero), i.e. that the derivative of a constant vanishes. Without this condition the boundary row-sum constraints are not enforced, and the stated uniqueness and free-parameter counts in Section III.A and Table I do not follow from the equations as written. The authors should correct the stated q-range (probably q=0,...,b) and confirm that the actual derivation and the supplementary materials impose the q=0 condition.","section":"III.A, Eqs. (18)-(19)"},{"comment":"The 3D convergence/robustness test does not exercise the embedded-boundary closure as a genuine IBVP. Section IV.A states that points just inside the excision boundary are active but 'evolved via the exact time derivatives obtained by differentiating the reference solution', and both the inner and outer boundary data are injected as exact incoming characteristics. Thus the stencils that include these boundary-layer points are fed with analytically prescribed values rather than values generated by the PDE, so any unstable boundary mode is suppressed by construction. The convergence rates in Figures 4-15 validate the interior stencils, multiblock interfaces, and boundary interpolation against exact data, but they do not demonstrate stability of the embedded-boundary operator when coupled with characteristic data derived from the interior solution. Please add a test in which all active poin","section":"IV.A and IV.B"},{"comment":"The paper itself concedes in Section III.B that 2D/3D stability with embedded boundaries 'cannot be proven' and relies on 'the hope that with appropriate numerical dissipation, the scheme can be stabilized.' Section IV.A also states that the test geometry avoids trapped points, slender geometries, and narrow domains. Given these limitations, the conclusion in Section V that the operators 'can be used to evolve hyperbolic problems in first order form with arbitrarily shaped boundaries' is broader than the evidence supports. Either provide additional stability evidence for a representative non-benign configuration (for example an inclined or tilted boundary, or a domain requiring the derivative-extrapolation treatment for trapped points), or narrow the stated claims to the tested spherical-excision setting with injected inner solution.","section":"III.B and IV.A"}],"minor_comments":[{"comment":"Typo: 'dissapative' should be 'dissipative' in the phrase 'asymmetric dissapative scheme'.","section":"III.B"},{"comment":"Typo: 'nontheless' should be 'nonetheless'.","section":"IV.C"},{"comment":"The distinction between the minimal-width operators and the DRP operators is hard to follow in plain text because the two families appear with nearly identical notation. Please make the notational difference explicit in all table captions, figure captions, and the text.","section":"III.A.2 / Table I"},{"comment":"The sentence describing active points inside the excision boundary should be more prominent: it is the key reason the subsequent convergence test is not a genuine IBVP test for the embedded boundary closure. At minimum, the authors should state what is injected and why this does not test feedback from the interior solution.","section":"IV.A"}],"recommendation":"major_revision","confidential_remarks":"The operator construction appears sound and the supplementary materials are a strength. The main obstacle is evidential: the 3D test is a manufactured-solution test with a filled-in boundary layer, so the robustness claim for excision applications is not yet supported. This is repairable with additional experiments or a substantial tightening of the claims; I do not see grounds for rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis is a genuine step forward: the first upwind and DRP embedded-boundary SBP operators, with interior orders 2–9 and boundary closures up to order 4. The construction pipeline is sound—ansatz, accuracy conditions, error-function minimization—and SBP holds by construction. The operators have negative semi-definite S and positive H on the stated α ranges, and the 3D convergence tests match the designed orders (b+1.5, b+1) on a cubed-sphere grid with an inner excision sphere. That is a solid, useful result.\n\nThe main soft spot is the test design. In Section IV.A, points just inside the excision boundary are evolved using exact time derivatives of the reference solution, and both inner and outer boundaries are injected as exact incoming characteristics. That means the layer feeding the boundary stencils is prescribed analytically. Any unstable mode that would arise from coupling the Cartesian stencils to an arbitrarily placed boundary is suppressed because the boundary data is exact, not computed from the interior evolution. The measured convergence rates therefore reflect the interior stencils and interpolation accuracy, but the stability of the embedded boundary closure in a self-consistent evolution is not demonstrated. The paper honestly states that 2D/3D stability with embedded boundaries cannot be proven and that dissipation is hoped to suffice; that is fine, but the numerical evidence for that hope is thinner than the abstract suggests.\n\nThere are smaller issues. Section III.A.2 contains duplicated text from III.A.1: the D±4−2/5−2 and D±6−3/7−3 descriptions are copied nearly verbatim with only the labels changed. No baseline comparison to the centered embedded operators of [13] is given, which makes the claimed spectral improvement harder to evaluate. And “highest order ... derived to date” is plausible but a little strong without a broader literature sweep.\n\nWho gets value: people working on high-order finite differences with embedded boundaries—black hole excision, worldtube excision, aeroacoustics. The operator families are new, the derivation is careful, and the explicit acknowledgement of limitations is a credit. But the central robustness claim for realistic excision is not fully substantiated by the present numerical test.\n\nRecommendation: send to peer review. It merits serious referee time, but a good referee should push for a genuine IBVP test without exact-solution injection, or at least a much more prominent scope limitation. Also require the promised supplementary materials and fix the copy-paste text.","headline":"New upwind/DRP embedded-boundary SBP operators are a real advance, but the 3D test injects exact solution near the excision boundary, so stability for genuine IBVPs remains unproven.","tokens_in":18844,"tokens_out":3730,"would_cite":true,"duration_ms":41372,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65M06","65M12","83-08"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that new upwind SBP operators embed arbitrary boundaries at up to 9th interior order and demonstrate ~5.5-order convergence in 3D wave tests.","keywords":["upwind summation by parts","embedded boundary","finite difference","high order accuracy","excision boundary","numerical relativity","dispersion relation preserving","diagonal norm"],"falsifier":"Take the D±_{9-4} operator and evolve the scalar wave equation with an excision sphere that is shallow relative to a grid line (so trapped points appear), or with a moving sphere whose velocity is comparable to the wave speed, or with characteristic data injected from a curved metric rather than a plane wave; if the L2 error grows without bound or the convergence rate drops below design as resolution increases, the stability claim fails for that regime.","tokens_in":17812,"feed_emoji":"🌀","tokens_out":5487,"duration_ms":48735,"temperature":0.7,"pith_summary":"This paper tries to establish that upwind summation-by-parts operators can be built for embedded boundaries—surfaces that cut through a Cartesian grid rather than aligning with grid lines—at interior accuracy orders p=2 through 9, with boundary closure accuracy up to order 4. The construction uses a parameterized ansatz for the derivative-pair (Q+, H), enforces the 2(b+1) accuracy conditions, and sets leftover free parameters by minimizing boundary error. If correct, these are the highest-order embedded-boundary finite-difference methods derived to date, and the 3D curvilinear wave-equation tests with a spherical excision show total L2 convergence near order 5.5. The payoff is that black-hole excision and worldtube coupling in numerical relativity could run on plain Cartesian grids without boundary-conforming coordinate transforms.","feed_headline":"Upwind embedded-boundary operators reach 9th order","feed_subtitle":"New SBP family keeps Cartesian grids and excises spheres, showing ~5.5-order convergence in 3D wave tests.","key_machinery":"The carrying object is the embedded-boundary upwind SBP pair (Q±, H): two derivative operators sharing a diagonal norm H and a boundary operator B = -e_l e_l^T + e_r e_r^T built from p-th-order accurate interpolation vectors. With Q+ = -(Q-)T, SBP holds automatically, and the scheme's dissipation is the matrix S = H(D+ - D-)/2, which is negative semi-definite and independent of the boundary position alpha. The construction enforces the accuracy conditions (Q± + 1/2 B) X^q - H X^{q-1} = 0 for q=1..b, then minimizes the total boundary error f over the free parameters to pin down a unique operator.","core_discovery":"The central discovery is a family of diagonal-norm upwind SBP operators for embedded boundaries, labeled D±_{p-b} and DRP variants, for interior orders 2..9 and boundary orders b up to 4. The operators satisfy SBP by construction via the identity (Q± + 1/2 B) X^q = H X^{q-1} for q=1..b, and the remaining free parameters are fixed by minimizing the total boundary error f = sum over q=b+1..2b-1 of the L2 norms of the error vectors. The paper shows each operator has a positive-definite norm on a length-one interval of the boundary position parameter alpha, has S negative semi-definite (so it supplies numerical dissipation), and converges in the 3D excision test at the designed b+1.5/b+1 rates,","pith_inferences":["The stability evidence in 3D rests on a single benign geometry: a convex sphere with no trapped points, no slender gaps, and exact time derivatives injected at near-boundary points. A moving or shallow-angle excision boundary could excite embedded-boundary modes that the fixed dissipation may not damp, so the robustness claim is narrower than the operator existence claim.","Since the dissipation matrix is alpha-independent, implementers could tabulate operator coefficients once per alpha range and interpolate, avoiding repeated symbolic derivation in production codes—an optimization the paper does not pursue.","The DRP operators sacrifice accuracy and time-step for better phase-velocity fidelity; for long-distance gravitational-wave propagation they might outperform the maximal-order operators, but this needs a dedicated study with distant observers.","The reported convergence order ~5.5 for the 9th-order interior operator is limited by the 4th-order boundary closure; pushing boundary order b to 5 or higher would likely raise the asymptotic rate, but the paper does not show such operators."],"forward_implications":["Black-hole excision in numerical relativity could use a stationary Cartesian grid, eliminating the need for dual-frame coordinate evolution when the excision boundary moves.","Worldtube excision—coupling an analytical interior to a numerical exterior through characteristic boundary conditions—becomes a high-order option on Cartesian grids.","The upwind dissipation S is built into the operators and independent of alpha, so no separate artificial-dissipation operator is needed for stability in the tested regime.","Odd-order operators have better spectral radii than their even-order counterparts at the same interior stencil width, meaning larger stable time steps and equal memory loads, making them attractive for memory-bound codes.","Setting alpha = 0 recovers traditional boundary-conforming SBP operators, so the new operators generalize standard codes rather than replacing them."],"fun_headline_variants":["Upwind embedded-boundary operators reach 9th order","New SBP family improves wave equations on Cartesian grids","Upwind SBP schemes tame embedded boundaries","9th-order accuracy for complex domains on simple grids","Embedded boundaries get a high-order upwind boost"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The scheme's stability in 2D/3D with embedded boundaries is assumed rather than proven—energy stability is shown only in 1D, and the 3D robustness tests use a friendly sphere geometry with exact-solution injection at boundary points, avoiding the trapped-point and narrow-domain issues that can arise for moving or shallow boundaries.","fun_headline_variants_meta":{"raw":{"variants":["Upwind embedded-boundary operators reach 9th order","New SBP family improves wave equations on Cartesian grids","Upwind SBP schemes tame embedded boundaries","9th-order accuracy for complex domains on simple grids","Embedded boundaries get a high-order upwind boost"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000556,"raw_usage":{"total_tokens":2489,"prompt_tokens":753,"completion_tokens":1736,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":497,"completion_tokens_details":{"reasoning_tokens":1661}},"tokens_in":497,"tokens_out":1736,"duration_ms":12300,"temperature":1.0,"reasoning_tokens":1661,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T02:11:35.906591+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the D±_{9-4} operator and evolve the scalar wave equation with an excision sphere that is shallow relative to a grid line (so trapped points appear), or with a moving sphere whose velocity is comparable to the wave speed, or with characteristic data injected from a curved metric rather than a plane wave; if the L2 error grows without bound or the convergence rate drops below design as resolution increases, the stability claim fails for that regime.","supporting_citations":[],"review_version":1}