{"id":"10fca9bb-ea0f-4c3f-9e61-8e7cfb2e7cf0","arxiv_id":"2607.24600","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"low","formal_verification":"none","parameter_count":4,"one_line_summary":"Trixi.jl gains adaptive multiphysics coupling via boundary exchange and user coupling functions, with runtime model switching that cuts cost versus full complex-model runs.","lead":"The authors add adaptive multiphysics coupling to the DG code Trixi.jl so many domains with different equations can exchange boundary data through user-defined converters, and the active model can change during a run. That lets cheaper models cover most of the domain while a costly model tracks only where it is needed, cutting runtime in their Euler–MHD demo by about 2.5×.","discovery_kind":"new_method","skeptic_critique":{"model":"moonshotai/kimi-k3","headline":"The accuracy evidence is load-bearing but nearly unfalsifiable as constructed: every demo uses smooth, ~1%-amplitude waves, and the Euler–MHD reference test only differs from the coupled run where B≈0 — precisely where the coupling does no work.","rationale":"The reader's weakest_assumption identified the same soft spot: no coupling-consistency or stability analysis, validation limited to visual continuity and one density-difference map, hand-tuned AMS parameters. My pass confirms this is the genuinely load-bearing concern rather than a peripheral one, and sharpens it in two ways the reader only implied: first, the Euler–MHD validation is structurally near-unfalsifiable because the two models coincide wherever B is negligible, so the reported 10⁻⁶ agreement mostly certifies the AMS tuning rather than the coupling; second, the specific mathematical risk is over-determined boundary data (full-state copy without characteristic/flux treatment), which predicts order-jump reflections for shocks crossing interfaces — a regime the paper advertises (astrophysics) but never tests. This does not make the paper wrong: as an engineering proof of concept with honest scope caveats (structured meshes, proof of concept stated in §7), CONDITIONAL remains the right verdict. The missing piece is a single hard test — a discontinuity crossing the interface with a refinement study — which would either shore up the central claim or show it holds only in the smooth, weakly-coupled regime actually demonstrated. No fabrication, no circularity, and the overhead/speedup measurements, while modest in scale, are reported with enough specificity (per-step ms, AMS check cost, percentages) to be credible. Hence: agree with the reader, verdict unchanged, with the concrete shock-refinement test as the decisive check.","tokens_in":12085,"tokens_out":2123,"duration_ms":80090,"concrete_test":"Run a Brio–Wu-type MHD shock tube (or Sod shock for the Euler–polytrope pair) with the initial discontinuity placed inside the MHD/complex region, propagating across a static coupling interface into the Euler region; compare against a monolithic full-domain reference at three resolutions (e.g. 64, 128, 256 cells per subdomain). Measure L1 and L∞ error in ρ and p in a window downstream of the interface plus the amplitude of any wave reflected back from the interface. If the coupling is sound, the error should converge at the scheme's design order and the reflected amplitude should vanish under refinement; if the reflected wave scales with the shock jump and stalls at O(1) under refinement, the one-shot full-state Dirichlet exchange is the cause and the central accuracy claim fails for discontinuous solutions. As a cheap adjunct, repeat §6.1 with t_delay halved/doubled and report the max |","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim — that one-shot Dirichlet-style boundary exchange with user converters yields solutions \"close to\" a monolithic complex-model reference — rests on tests that cannot meaningfully stress the coupling. (1) The polytrope–Euler (§5.2) and 3×3 array (§5.3) tests use linear, smooth waves of amplitude ~0.01; \"no observed discontinuities\" in Figures 4 and 6 is a visual statement with no error norm, no mesh-refinement study, and no interface-reflection diagnostic. (2) The headline Euler–MHD comparison (§6.1, Fig. 11) reports max |Δρ| < 1.1×10⁻⁶ at t=8, but the initial condition is uniform ρ=1, p=1 with a smooth magnetic ring (eq. 20) confined to the MHD region. Outside the ring, MHD with B≈0 reduces identically to Euler, so by construction the two runs can only disagree where (a) weak acoustic waves cross the interface or (b) the AMS criterion mis-tracks the field — and the thresholds plus t_delay=0.2 were hand-tuned for exactly this configuration. The 10⁻⁶ figure therefore measures how well-tuned the AMS is for this one problem, not the accuracy of the coupling mechanism. (3) For hyperbolic systems, copying the full conserved state on both sides of an interface over-determines the boundary data (no characteristic-based upwinding or flux matching); with genuinely nonlinear solutions — shocks or contact discontinuities crossing the interface, which is exactly the \"highly dynamical\" astrophysics regime the abstract advertises — spurious reflections of order the jump amplitude are expected. Nothing in §§5–6 tests this regime. (4) AMS model switches also silently drop or inject magnetic energy (B set to 0 in reclaimed cells, B discarded when reverting to Euler), and no conservation or energy-budget check at switches is reported. The speedup numbers (~2.5×) are believable but measured on a low-resolution, single-threaded laptop case, so the accuracy half of the central claim is the insecure part.","agreement_with_reader":"agree"},"referee_report":{"model":"moonshotai/kimi-k3","summary":"The manuscript presents an implementation, in the open-source DG code Trixi.jl, of two capabilities: (i) coupling of n hyperbolic subsystems on different subdomains (\"mesh views\" of a parent mesh) through exchange of boundary values composed with user-defined converter functions c_ab/c_ba, and (ii) adaptive model selection (AMS), in which the subdomain assigned to each model is resized at runtime according to user criteria, with the solution remapped between models via the same converters. Three demonstrations are given: a polytropic–Euler wave coupling (§5.2), a 3×3 array of polytropic/isothermal systems (§5.3), and an adaptively coupled Euler–MHD simulation of a travelling magnetic ring (§6.1). The authors report visual continuity at interfaces, a maximum density difference of 1.1×10⁻⁶ against a full-MHD reference at t=8, and a ~3× per-step (2.5× wall-clock) speedup with coupling/AMS overhead of a few percent. The paper is a proof of concept and says so in the conclusions.","tokens_in":12569,"tokens_out":2982,"duration_ms":108805,"significance":"If the validation is strengthened, this is a useful contribution to the Trixi.jl ecosystem: runtime adaptive model selection on resizable mesh views is a genuinely practical idea for multiphysics hierarchies (Euler↔MHD being the natural astrophysics target), the minimally invasive design means non-coupled users are unaffected, and the reported overhead numbers (0.8% coupling, 4.2% AMS checking, 2.5× wall-clock saving in §6.2) are concretely documented. The converter-function abstraction (§5.1, Eqs. 4–7) is simple and general. These are real strengths. However, in its current form the accuracy evidence does not yet support the paper's central claim that the coupled solution is close to a monolithic complex-model reference, because all three tests use smooth ~1%-amplitude waves and the headline Euler–MHD comparison is constructed so that the two runs can only differ where the coupling does essentially no work. As a methods paper whose selling point is exactly the coupling mechanism, this gap is load-bearing but fixable within the manuscript's scope.","major_comments":[{"comment":"§6.1, Fig. 11: the headline accuracy figure (max |Δρ| < 1.1×10⁻⁶ at t=8 vs full MHD) does not measure what the text implies it measures. The initial condition (uniform ρ=1, p=1, magnetic ring Eq. (20) confined to the central region) means that outside the ring the full-MHD reference evolves under equations that reduce identically to Euler when B≈0. The two runs can therefore only disagree where (a) the weak (~1% amplitude) acoustic waves generated by the non-force-balanced ring cross an interface, or (b) AMS mis-tracks the field — and the extend/shrink thresholds plus t_delay = 0.2 were chosen for exactly this configuration. The 10⁻⁶ number thus primarily certifies AMS tuning on this one problem, not the accuracy of the boundary-exchange coupling. The paper needs at least one test that actually stresses the coupling: a nonlinear solution (e.g. a shock or contact wave, or at minimum a fin","section":"§6.1, Fig. 11"},{"comment":"§5.1 and §4: the coupling mechanism copies the full conserved state (composed with converters) as boundary data on both sides of the interface, evaluated at every Runge–Kutta stage. For hyperbolic systems this over-determines the boundary data: a hyperbolic IBVP at an outflow boundary admits no prescribed incoming characteristic, and imposing the full state from the neighbour is a Dirichlet condition that is well known to generate spurious reflections for non-smooth solutions. The manuscript neither discusses this nor provides any characteristic-based or flux-matching alternative (contrast the cited Roux [6] iterative value/flux matching for the heat equation). Smooth 1%-amplitude waves will not reveal this. At minimum the paper needs (i) an honest discussion of the well-posedness/reflection issue, and (ii) a diagnostic quantifying interface reflection (e.g. reflection coefficient for a","section":"§5.1, §4"},{"comment":"§4, final paragraph: the claim that exchanging boundary values at every stage of the time integrator 'preserves the convergence order of the time integration' is stated as a consequence but is never demonstrated — there is no convergence study anywhere in the paper. This is an easily checkable assertion (e.g. a manufactured solution or a smooth wave with known dispersion, with a temporal-order study at fixed spatial resolution, coupled vs uncoupled), and since it is one of the few mathematical claims in the paper it should be either verified numerically or downgraded to an expectation.","section":"§4"},{"comment":"§6, AMS remapping: when an MHD cell is converted to an Euler cell, the magnetic field (and its energy B²/2 and the GLM variable Ψ) is discarded; when converting back, B is set to zero. The manuscript does not discuss or quantify the resulting conservation error in total energy/momentum across AMS events, nor its accumulation over many switches. Since AMS is advertised for 'highly dynamical scenarios', repeated switching is the expected regime. A plot of the global conserved quantities (mass, momentum, total energy) versus time for the §6.1 run, with the jumps at AMS events identified, would either close this concern or bound it.","section":"§6, Fig. 7"}],"minor_comments":[{"comment":"Reproducibility: no code-availability statement. Please state the Trixi.jl version/branch containing the coupling features and provide the setup scripts for §§5–6 (the Trixi ecosystem convention of a reproducibility repository would be ideal).","section":"§6.2"},{"comment":"Literature context is thin for coupling software: frameworks such as preCICE, MUI, and MpCCI are standard references for partitioned multiphysics coupling and should be discussed alongside [1], including how the present boundary-copy approach compares to their coupling schemes.","section":"§1"},{"comment":"§5.2, Eqs. (13)–(15): the initial density uses exponent 1/γ_E while the polytropic side uses γ_p = 2; please double-check and comment on why this specific mixed form is the right initialization (it appears to seed a wave consistent with the Euler side, which is worth saying explicitly).","section":"§5.2"},{"comment":"§6.1: please report the actual numerical values of the AMS extend/shrink thresholds (currently only described qualitatively) and add a brief sensitivity remark — e.g. how |Δρ| in Fig. 11 changes if t_delay is halved or doubled.","section":"§6.1"},{"comment":"§6.2: the timing comparison is a single laptop run. Please state whether the numbers are single measurements or minima/medians over repeats, and note that the 2.5× factor depends on the MHD-domain fraction and will shrink as the MHD region grows (a sentence on this scaling would calibrate the reader's expectations).","section":"§6.2"},{"comment":"Figures 4, 6, 10: the colorbar ranges span ±0.01 around 1, which is appropriate, but the claim 'no discontinuities at the interface' is visual; adding a 1D cut through the interface (density or pressure profile) would make the statement checkable from the figure.","section":"Figs. 4, 6, 10"},{"comment":"Typos/phrasing: §2 'coupled at x=0 and x=−2 with x=2' is hard to parse; §4 'minimally invasive' claim is repeated several times; abstract 'arbitrarily many domains' could be softened to 'an arbitrary number of subdomains'.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The engineering contribution (mesh views, AMS, Trixi.jl integration, low overhead) appears solid and fits a methods paper, but the validation is currently calibrated so that it cannot fail: all tests are smooth and small-amplitude, and the one quantitative comparison (Fig. 11) differs from the reference only where the coupling is inactive. I do not suspect bad faith — the paper calls itself a proof of concept — but the abstract's promise of 'highly dynamical scenarios' (shocks, astrophysics) is precisely the regime the tests avoid. The revisions I request (one nonlinear interface test with error norms, a convergence check of the §4 order claim, and a conservation diagnostic for AMS switches) are all within the existing code's capability and would decide whether this is a trustworthy tool or a demonstration of tuning. If the authors cannot produce a passing nonlinear test, the honest fix is to scope the claims to smooth regimes, which would weaken but not kill the paper."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The real deliverable here is a working multiphysics layer inside Trixi.jl: mesh views over a parent mesh, user converter functions for unequal variable sets, stage-tight boundary exchange that keeps the time integrator’s order, and runtime adaptive model selection. That package is new as concrete software, not as abstract idea. Boundary multiphysics and model hierarchies already exist; what they shipped is a clean, minimally invasive design with three demos and honest timing (~3× per step, ~2.5× wall-clock vs full MHD, coupling/AMS overhead a few percent).\n\nWhat they do well is practical. Converters that map Euler↔polytrope and Euler↔MHD without shared variables are the right abstraction for hierarchical models. The 3×3 array and the moving MHD ring show the machinery actually runs. They state the scope limit (structured meshes, proof of concept) instead of overselling.\n\nThe soft spot is the accuracy half of the central claim, and the stress-test is mostly right. All demos are smooth ~1% waves. “No discontinuities” is visual; there is no interface error norm, refinement study, or reflection diagnostic. The headline |Δρ|<1.1×10⁻⁶ vs full MHD is weak evidence: outside the ring B≈0, so MHD collapses to Euler and the runs can only disagree on weak acoustics and AMS tracking—exactly what the hand-tuned thresholds and t_delay=0.2 were set for. For hyperbolic systems, copying the full state as Dirichlet data over-determines the boundary; shocks or contacts crossing the interface are the regime the abstract advertises and are untested. AMS also drops/injects magnetic energy with no conservation check. Speedup numbers look real for the laptop case they ran.\n\nThis is for people already in the Trixi/DG multiphysics lane who want a scheduling pattern for cheap/expensive models. Not a theory paper. It deserves a serious referee—methods software with shipped design and measured overhead—but referees should demand a nonlinear interface test and a clearer statement of what the coupling guarantees. I would engage if I were extending Trixi; I would not treat the accuracy claim as settled.","headline":"Useful Trixi.jl engineering for adaptive multiphysics, with believable speedups but accuracy evidence that barely stresses the coupling.","tokens_in":12492,"tokens_out":546,"would_cite":false,"duration_ms":22493,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65M60","65M99","76N15","76W05"],"pacs":[],"model":"grok-4.5","headline":"Adaptively coupling different hyperbolic physics on resizable subdomains cuts runtime while keeping interface solutions smooth.","keywords":["multiphysics coupling","hyperbolic systems","adaptive model selection","discontinuous Galerkin","boundary exchange","mesh views","CFD","simulation methods"],"falsifier":"Rerun the adaptive Euler–MHD magnetic-ring test at higher resolution and longer time; growing jumps in density or magnetic energy at the moving interfaces, or a density difference from full-MHD far above the reported 10^-6 level, would falsify the central claim.","tokens_in":12101,"feed_emoji":"⚡","tokens_out":835,"duration_ms":41383,"temperature":0.7,"pith_summary":"This paper shows how to run many hyperbolic conservation laws side by side when the equations, and even the sets of unknowns, differ. Coupling is done by exchanging boundary states through user-written converter functions, so systems need not share a single variable. Subdomains sit on lightweight mesh views of one parent mesh and can grow or shrink at runtime under user criteria, letting the expensive model follow features such as a traveling magnetic ring. In Euler–polytropic and adaptive Euler–MHD tests the interfaces stay free of visible jumps, density stays within about one part in a million of a full-MHD reference, and wall-clock time drops by roughly 2.5× because the complex model is used only where needed. A reader who simulates dynamical multiphysics cares because the method buys that flexibility at a few percent overhead instead of paying the full complex-model cost everywhere.","feed_headline":"Coupled cheap and expensive physics cut runtime 2.5×","feed_subtitle":"Moving subdomains and boundary converters track local MHD features without full-domain cost","key_machinery":"Mesh views—subdomain handles on a fixed parent mesh that can be resized without rebuilding the grid—together with user-defined coupling functions that map state vectors across interfaces that may share few or no variables.","core_discovery":"Arbitrarily many hyperbolic systems can be coupled through interface boundary exchange plus algebraic user coupling functions, and the active model on each subdomain can be switched at runtime by resizing mesh views; the resulting discrete solutions remain free of interface artifacts and stay close to a monolithic complex-model reference while cutting per-step and wall-clock cost by large factors.","pith_inferences":["Pushing adaptive model selection to per-cell granularity and combining it with mesh refinement, as the authors note as next work, would further reduce cost when expensive physics occupies only sparse pockets.","Converters between truly non-overlapping variable sets (for example kinetic distributions and fluid moments) will be a stricter test than the shared-variable Euler–MHD case shown here.","Dynamic subdomain sizes can reintroduce load imbalance across participants; without a balancing strategy the wall-clock gain may shrink on large parallel runs."],"forward_implications":["Hierarchies such as multi-component to single-component gas can be coupled with simple algebraic density and pressure converters.","Traveling features that need an expensive model only locally can be tracked by growing and shrinking the corresponding mesh view.","When converters are cheap copies or algebraic maps, coupling and adaptivity overhead stays only a few percent of runtime.","The same pattern applies to arrays of many domains and to systems whose variable counts differ strongly."],"fun_headline_variants":["Adaptive multiphysics coupling swaps models mid-run to cut cost","Boundary exchange links arbitrary hyperbolic systems on the fly","Runtime model switching on subdomains slashes DG simulation time","Hierarchy of cheap-to-complex physics coupled without interface artifacts","Mesh-view resizing enables adaptive domain physics in Trixi.jl"],"cache_read_input_tokens":128,"weakest_assumption_plain":"A single-pass exchange of converted boundary values, plus hand-chosen thresholds and a short delay before shrinking expensive domains, is enough to keep the coupled solution physically consistent and close to a full complex-model run.","fun_headline_variants_meta":{"raw":{"variants":["Adaptive multiphysics coupling swaps models mid-run to cut cost","Boundary exchange links arbitrary hyperbolic systems on the fly","Runtime model switching on subdomains slashes DG simulation time","Hierarchy of cheap-to-complex physics coupled without interface artifacts","Mesh-view resizing enables adaptive domain physics in Trixi.jl"]},"model":"grok-4.5","effort":"low","cost_usd":0.004402,"raw_usage":{"total_tokens":1228,"prompt_tokens":687,"num_sources_used":0,"completion_tokens":88,"cost_in_usd_ticks":44024000,"prompt_tokens_details":{"text_tokens":687,"audio_tokens":0,"image_tokens":0,"cached_tokens":128},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":453,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":687,"tokens_out":88,"duration_ms":6970,"temperature":1.0,"reasoning_tokens":453,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-31T10:51:31.730695+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Rerun the adaptive Euler–MHD magnetic-ring test at higher resolution and longer time; growing jumps in density or magnetic energy at the moving interfaces, or a density difference from full-MHD far above the reported 10^-6 level, would falsify the central claim.","supporting_citations":[],"review_version":1}