{"id":"5fb443a8-3d34-4fec-9027-15f57d1e14e6","arxiv_id":"2608.05398","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A 3D multiphysics simulation of the ARC breeding blanket predicts about 243 mg tritium inventory and a 30-minute steady-state build-up time, with transport dominated by turbulence-enhanced mixing.","lead":"This paper couples three open-source simulation codes to model tritium movement through the molten salt blanket of an ARC-class fusion reactor, predicting a total inventory near 243 mg that builds up in about 30 minutes. It demonstrates a reusable workflow for testing how design choices affect tritium handling and safety.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 72-sector symmetry extrapolation is untested; if toroidal flow or source asymmetries exist, the 243 mg, 1.0 mg/s, and 30 min values do not scale linearly.","rationale":"The reader's CONDITIONAL verdict is appropriate. My stress-test agrees that the sector-to-full extrapolation is the most load-bearing concern. The other candidate, the numerical stabilisation scheme, is partially mitigated by the code-to-code agreement (FESTIM 243 mg versus OpenFOAM 240 mg), which uses a different discretisation and no artificial diffusion, so the central value is not a mere artifact of delta. The abstract's statement that inventories are 'governed primarily' by the stabilisation scheme is overstated relative to the paper's own sensitivity results and should be corrected, but it does not invalidate the 243 mg estimate. The sector symmetry assumption, by contrast, has no internal check and affects every system-level number. The paper's explicit lower-bound framing and its disclosure of neglected physics are creditworthy. Because the concern is already reflected in the reader's verdict, no change is needed.","tokens_in":16814,"tokens_out":8869,"duration_ms":83606,"concrete_test":"Run the same coupled OpenFOAM-FESTIM transport on a 10-degree sector (two adjacent 5-degree wedges) with the same mesh and boundary-condition strategy, and compare the per-degree steady-state inventory and outlet flux with the 5-degree case; if the per-degree values differ by more than about 5%, the linear 72-times extrapolation is invalid and the system-level numbers need revision. If feasible, also tally OpenMC tritium production in a full 360-degree model and compare per-wedge integrals to confirm toroidal source uniformity.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3 states that all system-level results are obtained by simulating a 5-degree toroidal sector and multiplying by 72, 'assuming perfect geometric and flow symmetry.' The CFD model imposes symmetry boundary conditions on the sector's toroidal faces, and the OpenMC source is treated as toroidally uniform. This is the load-bearing step between the simulated wedge and every headline number. No test is provided for the three embedded claims: (i) the 72 inlet/outlet pairs are hydraulically independent, (ii) the neutron source and heat deposition are toroidally uniform, and (iii) the geometry is exactly periodic with no manifolds or ports breaking symmetry. If any sector differs in flow rate, temperature, or tritium generation, the true inventory and outlet flux will not equal 72 times the sector value. The paper acknowledges the assumption but does not bound its error; it is not covered by the code-to-code comparison, which validates only the numerical solution on the same sector.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a fully open-source multiphysics workflow coupling OpenMC (neutron transport), OpenFOAM (thermal hydraulics), and FESTIM (tritium transport) to model tritium behavior in a 5-degree toroidal sector of the ARC liquid immersion blanket. From this sector, system-level results are obtained by scaling to 72 identical sectors: a steady-state tritium inventory of approximately 243 mg, an outlet mass flux of approximately 1.0 mg/s, an outlet-averaged concentration of approximately 0.90 mg/m^3, and a build-up time of approximately 30 minutes. The inventory is explicitly labeled a lower bound because of the zero inlet concentration and the neglect of tritium in solid structures. A code-to-code comparison with an OpenFOAM passive scalar solver agrees to within about 1% on inventory and outlet flux, and a sensitivity study varies the turbulent Schmidt number and the numerical stabilization parameter. The paper is candid about the modeling limitations, including the absence of solid trapping, permeation, and chemical effects.","tokens_in":17068,"tokens_out":5886,"duration_ms":55687,"significance":"If the central results are reliable, the paper provides a useful component-scale demonstration of a transparent, extensible, and fully open-source toolchain for tritium transport modeling in ARC-class blankets. The main strengths are the release of the two converter packages (foam2dolfinx and openmc2dolfinx), the careful code-to-code verification, the explicit lower-bound caveat, and the identification of flow-stagnation regions as tritium accumulation hotspots. The prediction of a ~243 mg inventory and ~1.0 mg/s throughput is a concrete, falsifiable result that can inform early design iteration. However, the headline numbers rest on the assumption of perfect 72-fold toroidal symmetry, which is acknowledged but not quantified, and the sensitivity discussion somewhat overstates the role of the numerical stabilization parameter relative to the turbulent Schmidt number.","major_comments":[{"comment":"All system-level values (243 mg inventory, 1.0 mg/s outlet flux, 30 min build-up) are obtained by simulating one 5-degree sector and multiplying by 72, under the stated assumption of perfect geometric and flow symmetry. This assumption is load-bearing for every headline number, yet no test or error bound is provided. The symmetry boundary conditions on the toroidal faces of the sector and the toroidally uniform treatment of the OpenMC source enforce the assumption rather than verify it. A concrete test, such as a two-sector simulation with perturbed inlet conditions, a coarse full-torus run, or a sensitivity study in which the sector's inlet velocity or tritium source is varied within a plausible range, would bound the extrapolation error. Without such a test, the reported system-level values remain conditional on an unquantified geometric and hydraulic symmetry.","section":"Section 3 and Section 2.2"},{"comment":"The abstract claims that predicted inventories are 'governed primarily by the numerical stabilisation scheme,' but the data in Figure 9 show that varying the stabilization parameter δ from 0.01 to 10 changes the inventory by -6.2% to +0.7%, while varying the turbulent Schmidt number Sct from 0.3 to 1.3 changes the inventory by -5% to +3%. These ranges overlap and are of comparable magnitude; the conclusion that the stabilization scheme is the dominant factor is therefore an overstatement. The text should be reworded to say that both parameters have modest effects, with δ having a somewhat larger influence in one direction, and the abstract should be aligned with this more balanced conclusion.","section":"Abstract and Section 3.2 (Figure 9)"},{"comment":"No mesh convergence or grid sensitivity study is reported for either the OpenFOAM CFD run or the FESTIM transport solve. The 1% agreement between FESTIM and OpenFOAM demonstrates that two discretizations of the same equation on the same underlying mesh agree, but it does not establish that the mesh adequately resolves the velocity and concentration fields. A refinement study, or at least a report of the number of cells and element sizes in the computational meshes, is needed to support the claimed absolute accuracy of the 243 mg inventory and 1.0 mg/s flux.","section":"Sections 2.2, 2.3, and 3.2"}],"minor_comments":[{"comment":"Equation (3) defines the artificial diffusion as Dart = δ h ||u||, which is an isotropic artificial diffusion term, but the text describes it as 'analogous to streamline upwind Petrov–Galerkin (SUPG) or artificial diffusion methods.' SUPG is a consistent streamline-upwind scheme, not isotropic diffusion, so the paper should clarify which stabilization is actually implemented, or use a more precise description.","section":"Section 2.3, Eq. (3)"},{"comment":"The inlet turbulent kinetic energy and specific dissipation rate are specified as kin = 0.12 m^2/s^2 and ωin = 0.6 1/s, said to correspond to 20% turbulence intensity, but no turbulent length scale or hydraulic diameter is given; without this, the ωin value cannot be independently reproduced.","section":"Section 2.2, Table 3"},{"comment":"The color scale and label placement in Figure 9 make it difficult to associate the OpenFOAM comparison points with specific Sct values; a clearer legend or annotations would improve readability and help the reader verify the claimed agreement.","section":"Section 3.2 and Figure 9"},{"comment":"The paper states that the reported inventory is a 'conservative minimum'; since 'conservative' is ambiguous in a safety context, consider using 'lower bound' consistently to avoid implying that a lower inventory is necessarily conservative from a regulatory perspective.","section":"Section 3.2"},{"comment":"The paper compares the FESTIM outlet flux with the neutronic production rate and reports a 1.5% imbalance, but the absolute value of the production rate is not stated; reporting it would put the imbalance in context.","section":"Section 3.2"}],"recommendation":"major_revision","confidential_remarks":"The paper fits the scope of physics.comp-ph and makes a useful open-source contribution. The main technical weakness is the unverified 72-sector symmetry extrapolation that underlies all system-level numbers; this is fixable with additional analysis or a clear uncertainty bound. The sensitivity claim in the abstract also needs to be toned down. The code-to-code comparison and the candid treatment of limitations are strengths. I recommend major revision rather than rejection because the central workflow and the code-to-code verification are sound, and the required additional evidence is within the scope of a revised manuscript."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick verdict: worth reading closely, and worth refereeing. The genuinely new piece is the integrated OpenMC–OpenFOAM–FESTIM workflow, with two small converter packages (foam2dolfinx, openmc2dolfinx) that pass fields between solvers without interpolation loss. The headline result is the first 3D component-scale tritium transport simulation of the ARC liquid immersion blanket, giving 243 mg inventory, ~1.0 mg/s outlet flux, and ~30 min build-up time. The code-to-code verification against an OpenFOAM passive-scalar solver is credible: inventory 243 vs 240 mg, flux 1.00 vs 1.01 mg/s, and similar transient rise. The authors are also candid about the idealized inlet condition and lower-bound framing, which is more transparent than most modelling papers.\n\nSoft spots, in order of real importance:\n\n- The abstract overstates the sensitivity result. It says inventories are 'governed primarily by the numerical stabilisation scheme', but the δ sweep changes inventory by only a few percent, and the text itself calls the results 'relatively robust' to δ. That wording should be fixed.\n\n- The 72-sector symmetry extrapolation is load-bearing for every system-level number, yet it is untested. The CFD imposes symmetry on the sector faces, the neutron source is toroidally uniform, and the paper assumes all 72 sectors are hydraulically independent. The code-to-code comparison validates only the sector. An honest error bar on 243 mg would require at least a few perturbed-sector runs or a quantified argument about toroidal uniformity. The authors acknowledge the assumption but don't bound it.\n\n- 'Fully open-source' overclaims: the converter packages are released, but the full case (CAD, OpenMC tallies, OpenFOAM case, FESTIM script, commit hash) is not shipped. Reproducibility is partial until that appears.\n\nNone of these are fatal. The physical story—turbulence-enhanced diffusion dominates, stagnation regions accumulate, residence time matters more than local source strength—is well supported by the sector-level results and the verification. The 243 mg number should be read as order-one-hundred milligrams, not a precise prediction.\n\nThe paper is for fusion blanket designers and anyone assembling partitioned multiphysics workflows. It deserves peer review, and the revision is straightforward: correct the abstract, add a symmetry-uncertainty discussion, and release the case files. I'd recommend conditional acceptance.","headline":"A solid, transparent multiphysics modelling paper that delivers verified sector-level tritium predictions but overreaches in the abstract and on the 72-sector extrapolation.","tokens_in":17542,"tokens_out":3882,"would_cite":true,"duration_ms":33046,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A coupled multiphysics simulation of the ARC fusion blanket predicts a steady-state tritium inventory of about 243 mg in the liquid breeder, an outlet flux of about 1.0 mg/s, and a 30-minute build-up time.","keywords":["ARC","breeding blanket","tritium","multiphysics","FESTIM","molten salt","FLiBe","turbulent transport"],"falsifier":"Run the same workflow on a second sector with a perturbed neutron source or different inlet/outlet plenum boundary conditions and compare the per-sector inventory and outlet flux; a material difference would invalidate the 72-fold scaling that yields 243 mg and 1.0 mg/s.","tokens_in":16586,"feed_emoji":"⚛️","tokens_out":13050,"duration_ms":95259,"temperature":0.7,"pith_summary":"This paper constructs a fully open-source, component-scale model of tritium transport in the ARC fusion reactor's liquid FLiBe breeding blanket by coupling neutron transport, thermal hydraulics, and hydrogen-isotope transport solvers. Under steady-state conditions the model predicts a total tritium inventory of about 243 mg in the liquid breeder, an outlet tritium flux of about 1.0 mg/s, and a build-up time of about 30 minutes. Tritium is bred and carried by the molten salt, so knowing where it accumulates and how fast it leaves the blanket is central to fuel self-sufficiency and to sizing extraction systems. The paper also finds that turbulence-enhanced diffusion dominates molecular diffusion, so flow stagnation zones become tritium hot spots and well-mixed turbulent regions stay lean.","feed_headline":"ARC blanket holds 243 mg of tritium and releases 1.0 mg/s","feed_subtitle":"Open-source multiphysics model pins down tritium inventory, flux, and build-up time for the ARC blanket.","key_machinery":"The load-bearing mechanism is the advection–diffusion equation $\\partial c_m/\\partial t = \\nabla\\cdot(D_{\\mathrm{eff}}\\nabla c_m) + S - \\nabla\\cdot(\\mathbf{u}c_m)$ for the mobile tritium concentration $c_m$, with an effective diffusivity $D_{\\mathrm{eff}} = D + D_{\\mathrm{turb}} + D_{\\mathrm{art}}$. Here $D$ is the Fickian molecular diffusivity with Arrhenius temperature dependence, $D_{\\mathrm{turb}} = \\nu_t/Sc_t$ is the turbulence-enhanced diffusivity built from the CFD kinematic turbulent viscosity $\\nu_t$ and turbulent Schmidt number $Sc_t$, and $D_{\\mathrm{art}} = \\delta h\\|\\mathbf{u}\\|$ is an artificial diffusion added for numerical stability in the continuous Galerkin solve. The tritium source $S$ is taken from a neutron-transport tally and the velocity and turbulence fields from a CFD simulation of the same 5-degree sector, transferred into the finite-element solver through dedicated conversion tools. Because the effective diffusivity exceeds the molecular value by several orders of magnitude, this equation is what converts the computed flow structure into the predicted tritium distribution.","core_discovery":"On the paper's own terms, the central finding is a set of steady-state values for an ARC-class blanket sector: with zero tritium concentration at the inlet, representing a perfectly efficient extraction system, the liquid breeder holds roughly 243 mg of tritium, the outlet tritium flux is roughly 1.0 mg/s, and steady-state throughput is reached in about 30 minutes, with peak local concentrations around $2.5\\ \\mathrm{mg\\,m^{-3}}$ and an outlet-averaged concentration near $0.90\\ \\mathrm{mg\\,m^{-3}}$. These system-level numbers come from simulating a single 5-degree toroidal sector and scaling to the full 72-sector blanket. The paper argues that tritium concentration is governed by the balance of local production and residence time: the fast coolant-channel jet stays tritium-lean despite high generation rates, while recirculation and stagnation zones accumulate tritium. A code-to-code comparison with an independent finite-volume passive scalar solver yields outlet fluxes of 1.00 versus 1.01 mg/s and inventories of 243 versus 240 mg, which the authors present as mutual confirmation of the finite-element result.","pith_inferences":["If the assumed 5-degree sector symmetry fails in a real blanket (different neutron flux, plenum geometries, or wall temperatures across sectors), the 72-fold scaling that produces 243 mg and 1.0 mg/s will not hold, and a multi-sector simulation would reveal the error.","The artificial-diffusion stabilisation is not exactly conservative, with the finite-element outlet flux about 1.5% below the production rate at steady state, so replacing it with a discontinuous Galerkin scheme could shift the reported numbers slightly.","Because the model omits tritium trapping and permeation in solid structures, and the paper cites ~20% inventory increases from trapping in comparable blankets plus potentially larger neutron-damage effects, the full system inventory including structures could be several times the 243 mg liquid-only value.","A direct experimental check would be to measure the outlet tritium concentration of an ARC-like FLiBe loop under the modelled inlet conditions; agreement with the predicted roughly $0.90\\ \\mathrm{mg\\,m^{-3}}$ would support the turbulence and stabilisation treatment, while strong disagreement would point to missing physics such as trapping, chemistry, or a nonzero inlet concentration."],"forward_implications":["The 243 mg inventory is a lower bound: a real extraction system would leave a finite inlet concentration, raising the equilibrium concentration, inventory, and outlet concentration.","The predicted 30-minute build-up time is comparable to expected plasma pulse lengths, so in pulsed operation the blanket may need several pulses before reaching steady-state tritium throughput.","Since turbulence-enhanced diffusion dominates, flow stagnation zones are the main tritium accumulation risk and are the natural targets for geometry or flow optimisation.","The modest sensitivity to the turbulent Schmidt number across 0.3–1.3 means the bulk outlet behaviour is robust to this physical uncertainty, while the numerical stabilisation parameter has the larger effect and local mesh refinement is preferable to artificial diffusion.","The close code-to-code agreement between the finite-element and finite-volume solutions suggests the reported inventory, outlet flux, and build-up time are not artefacts of one discretisation scheme."],"supporting_citations":[{"why":"Defines the ARC immersion blanket concept and its operating conditions, providing the geometry and plasma parameters used here.","marker":"[1]"},{"why":"The earlier axisymmetric CFD and tritium transport estimate for the ARC blanket that this work extends and compares against.","marker":"[3]"},{"why":"The FERMI project supplied the coupled neutronics and thermal-hydraulics workflow, the OpenFOAM blanket model, and the heat-deposition mapping.","marker":"[16]"},{"why":"FESTIM is the finite-element hydrogen transport solver that this paper extends to solve the tritium advection–diffusion problem.","marker":"[18]"},{"why":"OpenMC performs the neutron transport simulation that produces the spatial tritium generation and heat deposition source terms.","marker":"[19]"},{"why":"OpenFOAM provides the CFD solver and turbulence closure used to generate velocity, temperature, and turbulent viscosity fields.","marker":"[20]"},{"why":"The hydrogen transport properties database supplies the molecular diffusivity parameters for tritium in FLiBe used in the model.","marker":"[35]"},{"why":"The independent OpenFOAM passive scalar transport solution used for the code-to-code verification of inventory and outlet flux.","marker":"[38]"},{"why":"The system-level tritium fuel cycle analysis that gives the blanket residence-time and build-up context for the roughly 30-minute result.","marker":"[41]"}],"fun_headline_variants":["Open-source model: ARC blanket holds 243 mg tritium, steady in 30 min","FESTIM-coupled simulation: ARC tritium 243 mg, steady in 30 min","Multiphysics model: ARC tritium 243 mg, buildup in 30 min","Open-source coupling predicts ARC tritium at 243 mg, steady in 30 min"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The full-blanket numbers are a linear extrapolation from one 5-degree sector to 72 identical sectors, so the entire prediction rests on perfect geometric and flow symmetry among sectors.","fun_headline_variants_meta":{"raw":{"variants":["Open-source model: ARC blanket holds 243 mg tritium, steady in 30 min","FESTIM-coupled simulation: ARC tritium 243 mg, steady in 30 min","Multiphysics model: ARC tritium 243 mg, buildup in 30 min","Open-source coupling predicts ARC tritium at 243 mg, steady in 30 min"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000934,"raw_usage":{"total_tokens":4024,"prompt_tokens":1001,"completion_tokens":3023,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":617,"completion_tokens_details":{"reasoning_tokens":2926}},"tokens_in":617,"tokens_out":3023,"duration_ms":19758,"temperature":1.0,"reasoning_tokens":2926,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T13:52:16.203436+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same workflow on a second sector with a perturbed neutron source or different inlet/outlet plenum boundary conditions and compare the per-sector inventory and outlet flux; a material difference would invalidate the 72-fold scaling that yields 243 mg and 1.0 mg/s.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"FESTIM is the finite-element hydrogen transport solver that this paper extends to solve the tritium advection–diffusion problem."},{"cited_title":"URLwww.openfoam.org","cited_arxiv_id":null,"evidence_quote":"OpenFOAM provides the CFD solver and turbulence closure used to generate velocity, temperature, and turbulent viscosity fields."},{"cited_title":"Sircar, J","cited_arxiv_id":null,"evidence_quote":"The independent OpenFOAM passive scalar transport solution used for the code-to-code verification of inventory and outlet flux."}],"review_version":1}