{"id":"fe43ed02-13b3-41c6-8432-c5edd0412395","arxiv_id":"2607.01591","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"NuDEAL framework delivers three complementary GPU-accelerated solvers (MOC/HFEM, DGMOC, DFEM-SN) for multigroup transport on unstructured meshes, with DFEM-SN showing eigenvalue errors below 50 pcm and the others matching large CPU cluster speeds on single GPUs.","lead":"The paper develops and verifies NuDEAL, a framework implementing three GPU-accelerated deterministic neutron transport solvers on unstructured meshes, tested on benchmarks and reactor models. A smart generalist might read it to understand how GPU computing can make high-fidelity whole-core reactor simulations more practical.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.3","headline":"Demonstrated problems may not reach the scale or heterogeneity of true whole-core simulations","rationale":"The reader's weakest assumption (benchmark results taken as sufficient evidence for target applications) directly matches the scale mismatch identified here. With the full manuscript now available, the concrete_test above can be performed directly on the reported problem sizes to decide whether the claim holds or needs qualification.","tokens_in":1763,"tokens_out":280,"duration_ms":22048,"concrete_test":"From the results section or tables, extract the spatial mesh size (elements or cells), number of energy groups, and total unknowns for the ABTR, Empire, and MSRE cases; if all are below ~10^6 spatial elements, the whole-core practicality claim requires an additional large-scale demonstration run.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim requires that accuracy and single-GPU performance on the reported cases (C5G7 plus ABTR, Empire, MSRE) extrapolate to practical whole-core heterogeneous advanced reactors. Empire is a microreactor and MSRE an experiment; both are orders of magnitude smaller than commercial whole-core models. If the mesh sizes, degrees of freedom, and runtimes in the results section correspond only to these smaller domains, the scalability assertion for whole-core use rests on an untested extrapolation rather than direct evidence.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The paper presents NuDEAL, a unified GPU-accelerated framework implementing three deterministic transport solvers (planar MOC/HFEM, DGMOC, and DFEM-SN) on unstructured meshes for solving the multigroup transport equation. It reports verification on the C5G7 benchmark (eigenvalue errors below 50 pcm for DFEM-SN) and applications to ABTR, Empire microreactor, and MSRE problems, with single-GPU runtimes claimed comparable to large CPU clusters, concluding that the methods enable practical whole-core simulations for heterogeneous advanced reactors.","tokens_in":1865,"tokens_out":299,"duration_ms":16715,"significance":"If the reported accuracy and performance metrics hold under independent verification, the work demonstrates concrete GPU acceleration for unstructured-mesh deterministic transport, providing a foundation for efficient high-fidelity neutronic analysis of advanced reactors and extensions to multiphysics.","major_comments":[{"comment":"Abstract and results discussion: the central claim that the solvers 'enable practical whole-core simulations for heterogeneous advanced reactors' is load-bearing but rests on extrapolation; the tested cases (C5G7 plus ABTR, Empire microreactor, MSRE) are orders of magnitude smaller than commercial whole-core models, with no reported scaling studies, larger-domain timings, or heterogeneity metrics to support the assertion.","section":"Abstract"}],"minor_comments":[],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for the constructive feedback on our manuscript. We address the single major comment below and agree that revisions are warranted to moderate the central claim.","responses":[{"response":"We agree with the referee that the tested problems are substantially smaller than commercial whole-core models and that the manuscript contains no explicit strong-scaling studies, larger-domain timings, or quantitative heterogeneity metrics to directly support the assertion. The claim was intended to reflect the observed single-GPU performance relative to large CPU clusters on the reported benchmarks together with the unstructured-mesh capability, but we recognize that this constitutes an extrapolation. We will revise the abstract, introduction, and conclusions to replace the phrasing 'enabling practical whole-core simulations' with 'providing a foundation toward practical whole-core simulations' and will add an explicit statement noting the absence of scaling studies to larger domains and the need for future work in that direction. These changes will be made without altering the reported verification and timing results.","revision_made":"yes","referee_comment":"[Abstract] Abstract and results discussion: the central claim that the solvers 'enable practical whole-core simulations for heterogeneous advanced reactors' is load-bearing but rests on extrapolation; the tested cases (C5G7 plus ABTR, Empire microreactor, MSRE) are orders of magnitude smaller than commercial whole-core models, with no reported scaling studies, larger-domain timings, or heterogeneity metrics to support the assertion."}],"tokens_in":1377,"tokens_out":311,"duration_ms":12493,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The main thing to know is that this paper delivers a practical GPU implementation of deterministic transport on unstructured meshes. It unifies MOC/HFEM, DGMOC, and DFEM-SN in one framework, adds GPU-specific optimizations such as memory alignment, compressed flux storage, and sequential sweeps, and reports concrete verification on C5G7 plus runs on ABTR, Empire, and MSRE. DFEM-SN hits eigenvalue errors below 50 pcm while the others show single-GPU times that match large CPU clusters. That is useful engineering output for reactor analysis.\n\nThe integration and the side-by-side performance data are the real additions. The methods themselves are established, but putting them together with these GPU tweaks and showing they can be swapped for accuracy versus speed is a reasonable contribution. The verification step uses standard benchmarks and reports direct error metrics, which is the right way to do it.\n\nThe soft spot is scale. Empire is a microreactor and MSRE an experiment; both are orders of magnitude smaller than commercial whole-core models. The abstract's claim that this enables practical whole-core simulations for heterogeneous advanced reactors therefore rests on extrapolation rather than direct evidence from large meshes. If the full paper contains bigger cases or explicit scaling studies, that would address the concern; otherwise the stress-test point stands.\n\nThis work is for nuclear computational physicists and engineers who need GPU transport tools on unstructured geometries. A reader looking for implementation details and benchmark comparisons would get value from the numbers and the framework description.\n\nIt deserves peer review. The verification data and performance claims are concrete enough to justify referee time, even if the scalability discussion needs tightening.","headline":"NuDEAL gives a working unified GPU framework for three transport methods with solid benchmark numbers, but the whole-core scalability claim extrapolates from smaller test cases.","tokens_in":2360,"tokens_out":410,"would_cite":false,"duration_ms":22570,"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":"GPU-accelerated deterministic solvers on unstructured meshes deliver accuracy and scalability for whole-core reactor simulations.","keywords":["deterministic neutron transport","GPU acceleration","unstructured meshes","method of characteristics","discontinuous Galerkin","finite element method","reactor core simulation","NuDEAL framework"],"falsifier":"A run of the C5G7 or MSRE problem on the same hardware showing eigenvalue error above 50 pcm for DFEM-SN or single-GPU runtimes substantially exceeding the reported CPU-cluster equivalents for MOC/HFEM or DGMOC.","tokens_in":2656,"feed_emoji":"⚛️","tokens_out":828,"duration_ms":20917,"temperature":0.7,"pith_summary":"The paper presents the NuDEAL framework containing three GPU-accelerated solvers for the multigroup transport equation on unstructured meshes: MOC/HFEM, DGMOC, and DFEM-SN. These methods are implemented with GPU optimizations such as memory alignment and sequential sweeps, then validated on the C5G7 benchmark before application to the ABTR, Empire microreactor, and MSRE problems. DFEM-SN reaches eigenvalue errors below 50 pcm while the other two methods match the speed of large CPU clusters on a single GPU. The work shows that deterministic transport on unstructured meshes can now support practical whole-core calculations for heterogeneous advanced reactors. The unified framework also sets up extensions to transient and multiphysics modeling on large GPU systems.","feed_headline":"GPU solvers achieve whole-core reactor accuracy on unstructured meshes","feed_subtitle":"NuDEAL's three methods reach eigenvalue errors below 50 pcm with single-GPU runtimes matching large CPU clusters on C5G7 and advanced reacto","key_machinery":"The NuDEAL unified framework that implements and selectively deploys MOC/HFEM, DGMOC, and DFEM-SN solvers, each optimized for GPU execution via memory alignment, compressed-flux storage, and sequential azimuthal sweeps to solve the multigroup transport equation consistently on unstructured meshes.","core_discovery":"NuDEAL unifies three complementary deterministic solvers—MOC/HFEM, DGMOC, and DFEM-SN—on unstructured meshes, all accelerated on GPUs through memory alignment, compressed-flux storage, and sequential azimuthal sweeps. Validation on C5G7 and application to ABTR, Empire, and MSRE show DFEM-SN achieving eigenvalue errors below 50 pcm while MOC/HFEM and DGMOC deliver single-GPU runtimes comparable to large CPU clusters, establishing that deterministic GPU solvers can enable practical whole-core simulations for heterogeneous advanced reactors.","pith_inferences":["The same GPU optimizations could be applied to deterministic transport in other fields such as radiative transfer or charged-particle problems.","Coupling the framework to existing multiphysics codes would allow testing whether the reported speedups persist under coupled iteration.","Users facing memory-limited problems could systematically compare the three solvers on the same mesh to quantify the accuracy-efficiency frontier.","Extension to multi-GPU or distributed GPU clusters would test whether the observed single-GPU scaling continues linearly with problem size."],"forward_implications":["DFEM-SN supplies the highest accuracy among the three solvers with eigenvalue errors below 50 pcm.","MOC/HFEM and DGMOC achieve the best computational efficiency, with single-GPU times matching large CPU clusters.","The framework permits choosing among the three solvers to trade accuracy against memory and runtime cost for a given geometry.","Whole-core simulations become practical for heterogeneous advanced reactors that require unstructured meshes.","The same code base provides a direct path to transient and multiphysics extensions on large-scale GPU hardware."],"fun_headline_variants":["NuDEAL unifies MOC/HFEM DGMOC and DFEM-SN on GPU meshes","DFEM-SN shows highest accuracy below 50 pcm in NuDEAL","MOC/HFEM and DGMOC match large CPU clusters on single GPU","NuDEAL framework tested on C5G7 and advanced reactor problems"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"The benchmark results are taken as sufficient proof that the solvers contain no significant implementation or numerical errors and that the tested problems adequately represent real advanced reactor applications.","fun_headline_variants_meta":{"raw":{"variants":["NuDEAL unifies MOC/HFEM DGMOC and DFEM-SN on GPU meshes","DFEM-SN shows highest accuracy below 50 pcm in NuDEAL","MOC/HFEM and DGMOC match large CPU clusters on single GPU","NuDEAL framework tested on C5G7 and advanced reactor problems"]},"model":"grok-4.3","cost_usd":0.009245,"raw_usage":{"total_tokens":4188,"prompt_tokens":765,"num_sources_used":0,"completion_tokens":79,"cost_in_usd_ticks":92449500,"prompt_tokens_details":{"text_tokens":765,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":3344,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":765,"tokens_out":79,"duration_ms":25783,"temperature":1.0,"reasoning_tokens":3344,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-03T02:38:48.281799+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"A run of the C5G7 or MSRE problem on the same hardware showing eigenvalue error above 50 pcm for DFEM-SN or single-GPU runtimes substantially exceeding the reported CPU-cluster equivalents for MOC/HFEM or DGMOC.","supporting_citations":[],"review_version":1}