{"id":"cde58918-98f3-4a7d-89d5-2b4974d87556","arxiv_id":"2606.28613","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":6.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":0,"one_line_summary":"NIMSTELL implements a semi-implicit MHD formulation for non-axisymmetric stellarators using nodal spectral elements and Fourier modes, verified linearly and nonlinearly on interchange and tearing instabilities against NIMROD and JOREK results.","lead":"The paper presents NIMSTELL, a semi-implicit MHD code for stellarators that combines 2D nodal spectral elements in the poloidal plane with Fourier representation in a generalized toroidal angle, generalizing the NIMROD code. Smart generalists might read it to understand advances in numerical tools for modeling unstable plasma dynamics in potential fusion devices.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.3","headline":"No significant objection identified","rationale":"Reader's weakest_assumption correctly flags the semi-implicit operator as the key unproven element, but the full text supplies targeted verifications on non-axisymmetric equilibria that directly address it; therefore the assumption does not rise to a load-bearing concern that would alter the UNVERDICTED verdict.","tokens_in":1860,"tokens_out":313,"duration_ms":18046,"concrete_test":"Re-run the W7-A tearing case (or the interchange benchmark) at the largest timestep reported in the paper but with the semi-implicit operator replaced by its axisymmetric (n=0 only) projection; if the growth rate or saturation amplitude deviates by more than the spatial convergence tolerance, the 3D-field construction is the source of the difference.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is a new semi-implicit formulation for stellarator MHD using 2D nodal spectral elements + Fourier in generalized toroidal angle, with the implicit operator drawn from the 3D ideal-MHD energy integral. The manuscript supplies linear/nonlinear verification on resonant interchange (stable-side convergence via the NIMROD stabilization) and on tearing in the W7-A rotating-ellipse case against JOREK, plus explicit statements on H1 vs H(curl) representations and preconditioning by stellarator mode families. No internal inconsistency, hidden assumption about axisymmetry, or untested regime for the semi-implicit operator is apparent from the supplied text.","agreement_with_reader":"disagree"},"referee_report":{"model":"grok-4.3","summary":"The manuscript presents a new semi-implicit formulation for nonlinear MHD computations in stellarators. It employs 2D nodal spectral elements over the poloidal plane with Fourier representation in a generalized toroidal angle, expanding geometric mappings and equilibrium fields in the same 3D basis. The semi-implicit operator is constructed from the ideal-MHD energy integral using 3D pressure and magnetic fields. The NIMSTELL implementation generalizes NIMROD, supports continuous H1 and optional H(curl) representations, and is verified linearly and nonlinearly on resonant ideal interchange (stable-side convergence) and on tearing modes in the W7-A rotating-ellipse case against JOREK. Preconditioning incorporates stellarator mode families.","tokens_in":1977,"tokens_out":570,"duration_ms":34480,"significance":"If the central claims hold, the work supplies a verified computational capability for macroscale MHD dynamics in genuinely non-axisymmetric stellarator equilibria, addressing a recognized need in fusion plasma modeling. The nodal spectral-element discretization enables systematic h- or p-refinement, the 3D energy-integral operator targets large-timestep stability, and direct comparisons to NIMROD and JOREK on standard interchange and tearing benchmarks provide concrete evidence of correctness. These elements collectively advance the toolkit for stellarator-specific nonlinear simulations.","major_comments":[{"comment":"The central claim of accuracy and stability at large timesteps for non-axisymmetric equilibria rests on the semi-implicit operator derived from the 3D ideal-MHD energy integral. The manuscript should supply quantitative data (e.g., timestep size relative to explicit CFL limits, growth-rate errors, and any numerical artifacts) for the W7-A tearing case and the interchange verification to substantiate this for genuinely 3D configurations.","section":"Verification (interchange and W7-A tearing)"},{"comment":"For the H(curl) vector-potential representation, the reported requirement of a minimum electrical resistivity to suppress numerical noise on interchange is a load-bearing limitation. The text should quantify the resistivity threshold relative to physical values and spatial resolution, and clarify whether this affects the method's applicability to ideal or low-resistivity regimes.","section":"H(curl) representation and interchange verification"}],"minor_comments":[{"comment":"Abstract contains the typographical error 'diUusive' (should read 'diffusive').","section":"Abstract"},{"comment":"The description of preconditioning by 'stellarator mode families' would benefit from an explicit statement of how the Fourier components are grouped and whether this choice is problem-dependent.","section":"Implementation and preconditioning"}],"recommendation":"minor_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for the positive evaluation and the recommendation of minor revision. We address the two major comments point by point below. Both points identify areas where the manuscript can be strengthened with additional quantitative detail and clarification; we will incorporate these changes in the revised version.","responses":[{"response":"We agree that explicit quantitative metrics would strengthen the verification claims for the 3D cases. The existing simulations contain the requested information (timestep sizes normalized to explicit CFL limits, growth-rate errors relative to NIMROD and JOREK, and notes on any observed artifacts). In the revised manuscript we will add a dedicated table or subsection presenting these data for both the resonant ideal interchange and the W7-A tearing-mode cases to directly support the large-timestep stability assertions.","revision_made":"yes","referee_comment":"[Verification (interchange and W7-A tearing)] The central claim of accuracy and stability at large timesteps for non-axisymmetric equilibria rests on the semi-implicit operator derived from the 3D ideal-MHD energy integral. The manuscript should supply quantitative data (e.g., timestep size relative to explicit CFL limits, growth-rate errors, and any numerical artifacts) for the W7-A tearing case and the interchange verification to substantiate this for genuinely 3D configurations."},{"response":"The manuscript already states that the H(curl) implementation requires a minimum resistivity to suppress noise on interchange at given resolutions. We will expand this discussion with concrete examples of the threshold (expressed, for instance, as a function of polynomial degree and element size, or in terms of magnetic Reynolds number) and will compare these values to typical physical resistivities encountered in stellarator plasmas. We will also add an explicit statement that the H1 representation is the appropriate choice for ideal or very-low-resistivity regimes, while the H(curl) option is intended for resistive MHD applications. These additions will clarify the method's scope of applicability.","revision_made":"yes","referee_comment":"[H(curl) representation and interchange verification] For the H(curl) vector-potential representation, the reported requirement of a minimum electrical resistivity to suppress numerical noise on interchange is a load-bearing limitation. The text should quantify the resistivity threshold relative to physical values and spatial resolution, and clarify whether this affects the method's applicability to ideal or low-resistivity regimes."}],"tokens_in":1598,"tokens_out":511,"duration_ms":35287,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"This paper delivers a concrete generalization of the NIMROD code to non-axisymmetric stellarator geometries. The key move is the mixed discretization—nodal spectral elements in the poloidal plane, Fourier in a generalized toroidal angle—plus carrying the equilibrium fields and mappings in the same 3D space so the semi-implicit operator can draw directly from the ideal-MHD energy integral with 3D pressure and B.\n\nIt does the basics right. Linear and nonlinear tests on resonant interchange recover the expected stabilization behavior from the parent code, and the tearing comparison against JOREK on the W7-A rotating-ellipse case gives an external benchmark. They are upfront about the two magnetic-field representations (H1 continuous versus H(curl) vector potential) and note that the latter needs a minimum resistivity to suppress noise at given resolution.\n\nThe soft spot is the leap from these tests to the claim of reliable large-timestep accuracy in genuinely 3D stellarator equilibria. The preconditioning by stellarator mode families is mentioned but not quantified, and the nonlinear evidence is still thin. If the semi-implicit operator introduces artifacts once the configuration is far from axisymmetry, that would only show up in longer or more strongly driven runs.\n\nThis is for people already running or extending NIMROD-type codes who need stellarator capability. A reader who knows the parent code and the JOREK comparisons will extract the most value. It is worth sending to peer review; the formulation is new enough and the checks are specific enough that referees can usefully push on the nonlinear robustness and scaling questions.","headline":"NIMSTELL extends NIMROD to stellarators via 2D nodal spectral elements plus Fourier toroidal with a 3D energy-integral semi-implicit operator, and the verifications on interchange and W7-A tearing look usable.","tokens_in":2467,"tokens_out":416,"would_cite":false,"duration_ms":22960,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"A semi-implicit MHD formulation with 2D nodal spectral elements and 3D ideal-MHD operator computes nonlinear stellarator dynamics at large timesteps.","keywords":["stellarator MHD","semi-implicit method","nodal spectral elements","non-axisymmetric equilibria","ideal-MHD energy integral","nonlinear computation","toroidal confinement"],"falsifier":"A run of the W7-A rotating-ellipse tearing-mode test at a timestep several times the explicit limit that produces either growing numerical noise or a nonlinear evolution visibly different from the JOREK reference solution would falsify the claim of maintained accuracy.","tokens_in":2764,"feed_emoji":"","tokens_out":819,"duration_ms":36265,"temperature":0.7,"pith_summary":"The paper introduces a computational approach for nonlinear time-dependent magnetohydrodynamics in stellarators that combines 2D nodal spectral elements in the poloidal plane with Fourier representation along a generalized toroidal angle. All geometric mappings and equilibrium fields are expanded in the same three-dimensional basis as the evolving fields to capture non-axisymmetric geometry. The semi-implicit operator is constructed directly from the ideal-MHD energy integral evaluated with the full three-dimensional pressure and magnetic field, which supports stable time steps much larger than explicit limits. The resulting NIMSTELL implementation generalizes an earlier axisymmetric code and is verified on resonant interchange modes and on tearing modes in the W7-A rotating-ellipse configuration, showing agreement with independent codes.","feed_headline":"Semi-implicit method advances stellarator MHD at large timesteps","feed_subtitle":"Nodal spectral elements in the poloidal plane and Fourier in toroidal angle use a 3D ideal-MHD operator for non-axisymmetric simulations.","key_machinery":"The semi-implicit operator derived from the ideal-MHD energy integral using 3D pressure and magnetic fields, paired with continuous H1 (or optionally H(curl)) nodal spectral elements in the poloidal plane and Fourier expansion in the generalized toroidal angle.","core_discovery":"The central claim is that a semi-implicit time advance for toroidally shaped MHD systems, discretized with 2D nodal spectral elements over the poloidal plane and Fourier modes in the generalized toroidal angle, remains accurate and stable at large timesteps when the implicit operator is taken from the ideal-MHD energy integral constructed with three-dimensional pressure and magnetic fields; the same 3D representation is used for all equilibrium and time-dependent quantities, allowing direct modeling of non-axisymmetric stellarator equilibria.","pith_inferences":["The method opens the possibility of following the evolution of MHD modes that are linearly unstable yet permit robust operation in stellarator experiments.","Because the semi-implicit operator uses the full 3D fields, the same framework could be applied to other toroidally confined devices whose equilibria depart strongly from axisymmetry.","The ability to switch between H1 and H(curl) representations within one code base allows direct comparison of how divergence control and gauge choice affect long-time nonlinear behavior."],"forward_implications":["Convergence is obtained by either increasing the number of elements (h-refinement) or the polynomial degree within each element (p-refinement).","Both the continuous H1 expansion of magnetic-field components with diffusive divergence control and the optional H(curl) vector-potential representation reproduce linear and nonlinear results from established codes on ideal interchange and tearing.","Algebraic systems arising from the implicit steps are preconditioned by including the Fourier components of each stellarator mode family.","The vector-potential option requires a minimum electrical resistivity at given spatial resolution to suppress numerical noise in interchange simulations."],"fun_headline_variants":["Semi-implicit stellarator MHD with nodal spectral elements","Stellarator MHD uses 3D ideal-MHD operator at large timesteps","Poloidal nodal elements with toroidal Fourier for stellarator MHD","Semi-implicit time advance for toroidally shaped stellarator MHD"],"cache_read_input_tokens":64,"weakest_assumption_plain":"The semi-implicit operator built from the ideal-MHD energy integral with three-dimensional fields stays accurate and stable at large timesteps for genuinely non-axisymmetric stellarator equilibria without unacceptable numerical artifacts.","fun_headline_variants_meta":{"raw":{"variants":["Semi-implicit stellarator MHD with nodal spectral elements","Stellarator MHD uses 3D ideal-MHD operator at large timesteps","Poloidal nodal elements with toroidal Fourier for stellarator MHD","Semi-implicit time advance for toroidally shaped stellarator MHD"]},"model":"grok-4.3","cost_usd":0.009881,"raw_usage":{"total_tokens":4466,"prompt_tokens":812,"num_sources_used":0,"completion_tokens":71,"cost_in_usd_ticks":98812000,"prompt_tokens_details":{"text_tokens":812,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":3583,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":812,"tokens_out":71,"duration_ms":36865,"temperature":1.0,"reasoning_tokens":3583,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-06-30T00:42:14.752803+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"A run of the W7-A rotating-ellipse tearing-mode test at a timestep several times the explicit limit that produces either growing numerical noise or a nonlinear evolution visibly different from the JOREK reference solution would falsify the claim of maintained accuracy.","supporting_citations":[],"review_version":1}