Zhigang Li , Junjie Pan , Shanfang Huang , Yingwei Wu , Xiafeng Zhou , Wei Lu , Bangyang Xia , Wei Zeng , Shenglong Qiang
{"title":"Cell-wise static coupling calculation of a new neutronics/thermal-hydraulics coupling benchmark for plate-type PWR core","authors":"Zhigang Li , Junjie Pan , Shanfang Huang , Yingwei Wu , Xiafeng Zhou , Wei Lu , Bangyang Xia , Wei Zeng , Shenglong Qiang","doi":"10.1016/j.anucene.2025.111895","DOIUrl":"10.1016/j.anucene.2025.111895","url":null,"abstract":"<div><div>To support the verification of cell-wise neutronics/thermal hydraulics (N-TH) coupling codes for Pressurized Water Reactor (PWR), the Nuclear Power Institute of China (NPIC) collaborated with multiple research institutions in the nuclear energy field to design and develop COPHP, a high-parameter PWR full core cell-wise N-TH coupling benchmark. This paper elaborates on Part II of the benchmark—cell-wise static N-TH coupling—including case design, modeling, results, and deviation analysis. Key conclusions are as follows: 1) Consistency was observed among the results from four participating codes (RMC/CORTH-R, OpenMC, CORCA-SPn/CORTH-R, and CORCA-SPn/TH1D) across 32 cases covering two core configurations. 2) Compared to RMC/CORTH-R: OpenMC exhibited a maximum <em>k<sub>eff</sub></em> deviation of −29.6 pcm and a maximum cell-wise power weighting error (PWE) of 0.185 %; CORCA-SPn/CORTH-R showed a maximum <em>k<sub>eff</sub></em> deviation of −474.6 pcm and a maximum PWE of 3.179 %; CORCA-SPn/TH1D yielded a maximum <em>k<sub>eff</sub></em> deviation of 1310.4 pcm and a maximum PWE of 5.245 %. 3) Two-phase flow at the core outlet in some cases of COPHP-B significantly affected the N-TH coupling results of CORCA-SPn/TH1D. 4) This benchmark suite provides comprehensive results for high-fidelity N-TH coupling in high-parameter PWRs, catering to diverse validation needs.</div></div>","PeriodicalId":8006,"journal":{"name":"Annals of Nuclear Energy","volume":"226 ","pages":"Article 111895"},"PeriodicalIF":2.3,"publicationDate":"2025-09-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145109119","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Md. Rahat Khan , Abid Hossain Khan , Ying Yin Tsui , Md. Zahurul Islam
{"title":"Automatic generation control of a nuclear-renewable hybrid power system using optimal PID controller","authors":"Md. Rahat Khan , Abid Hossain Khan , Ying Yin Tsui , Md. Zahurul Islam","doi":"10.1016/j.anucene.2025.111874","DOIUrl":"10.1016/j.anucene.2025.111874","url":null,"abstract":"<div><div>This study investigates the Automatic Generation Control (AGC) of a nuclear-renewable hybrid power system using optimized PID controllers combined with Superconducting Magnetic Energy Storage (SMES). The system is modeled in MATLAB Simulink using state-of-the-art transfer functions, with the pressurized water reactor (PWR) simulated using a point kinetics model to accurately capture dynamic behavior. Six load variation scenarios involving changes of 100–200 MW in the NPP and 5–20 MW in the renewables are analyzed to assess system stability. Optimization algorithms such as Particle Swarm Optimization (PSO), Genetic Algorithm (GA), and Grey Wolf Optimizer (GWO) are employed to fine-tune the PID gains. The results show improved frequency control: the NPP stays within 49.99–50.14 Hz, an improvement over the previous range of 49.97–50.48 Hz, while the renewables stabilize at 49.99–50.05 Hz, a significant improvement from 46.98–72.30 Hz. These findings confirm improved grid stability and control during sudden load changes.</div></div>","PeriodicalId":8006,"journal":{"name":"Annals of Nuclear Energy","volume":"226 ","pages":"Article 111874"},"PeriodicalIF":2.3,"publicationDate":"2025-09-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145109121","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Analysis of the continuous operational characteristics of the control rod eddy current retarder","authors":"Hongyu Wu, He Yan, Yujie Dong, Xingzhong Diao","doi":"10.1016/j.anucene.2025.111888","DOIUrl":"10.1016/j.anucene.2025.111888","url":null,"abstract":"<div><div>The control rod eddy current retarder (CRECR) serves as the primary damping device during the control rod dropping process in the high temperature gas-cooled reactor, converting the gravitational potential energy into eddy current losses inside the conductor disk of CRECR. In order to analyze the torque performance of CRECR during continuous operation, a multi-node coupling method for electromagnetic and thermal simulation was proposed, and an electromagnetic-thermal coupling simulation model was established. Subsequently, the model was experimentally verified and the influence of speed, operating time and compensation claddings were analyzed. The torque deviation remains consistently below 7 % within 300 s. Additionally, continuous operation at speeds exceeding 300 rpm resulted in non-ignorable changes in component temperature and torque, necessitating careful consideration in engineering tests. Furthermore, a torque rebound phenomenon was observed during continuous operation, indicating that using compensation claddings can effectively mitigate the impact of thermal fade and improve the torque stability.</div></div>","PeriodicalId":8006,"journal":{"name":"Annals of Nuclear Energy","volume":"226 ","pages":"Article 111888"},"PeriodicalIF":2.3,"publicationDate":"2025-09-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145109122","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A unified surrogate model for enhanced photon shielding through an all-natural-element multi-group flux dataset","authors":"Junyi Chen, Chenghao Cao, Shaoning Shen, Tianyuan Guo, Jingang Liang","doi":"10.1016/j.anucene.2025.111887","DOIUrl":"10.1016/j.anucene.2025.111887","url":null,"abstract":"<div><div>Efficient gamma shielding analysis is vital for nuclear safety. However, traditional Point Kernel approaches relying on outdated ANSI databases face significant limitations in providing comprehensive and fine-grained shielding analysis. This study addresses this by creating a detailed multi-group photon flux dataset for 92 nuclides via Monte Carlo simulations. The dataset’s complexity renders traditional modeling techniques ineffective. We introduce a novel generative-reconstruction surrogate model, combining a conditional Generative Adversarial Network (cGAN) and a fine-tuned UNet, both enhanced with self-attention mechanisms. This model predicts complex multi-group photon shielding parameter fields. Verification shows the model accurately predicts parameter fields, with 95% of samples achieving an average relative deviation below 20%. Predicted relative flux, converted to buildup factors, aligns well with Monte Carlo truth and ANSI values, confirming reliability and improved conservatism. This approach offers an efficient, accurate alternative for photon shielding calculations, proposing a new approach for data and computation.</div></div>","PeriodicalId":8006,"journal":{"name":"Annals of Nuclear Energy","volume":"226 ","pages":"Article 111887"},"PeriodicalIF":2.3,"publicationDate":"2025-09-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145099862","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Yiwei Wu , Qufei Song , Yuyang Shen , Hui Guo , Yao Xiao , Hanyang Gu
{"title":"Development of a MC/MOC two-step method for sodium fast reactor analysis with transport-corrected flux-moment multigroup cross-section generation method","authors":"Yiwei Wu , Qufei Song , Yuyang Shen , Hui Guo , Yao Xiao , Hanyang Gu","doi":"10.1016/j.anucene.2025.111853","DOIUrl":"10.1016/j.anucene.2025.111853","url":null,"abstract":"<div><div>New fast reactor designs could feature complex geometries and axial heterogeneity. The MC/MOC two-step method has the advantage of the high adaptability of Monte Carlo (MC) and Method of Characteristics (MOC) to geometry, offering a precise solution for such reactors. However, limited studies exist on applying the MOC method to 3D sodium fast reactors. MC-generated multi-group cross-sections, which ignore the anisotropy with respect to neutron direction, introduce bias in transport calculations. This paper develops a two-step MC/MOC approach, integrating the direct 3D MOC method and the Transport-corrected flux-moment homogenization technique (TC-MHT) for 3D reactor problems. The scheme is verified on the sodium-cooled fast reactor MET-1000, the results confirm the feasibility of the method for both homogeneous and heterogeneous core geometries. The TC-MHT method significantly reduces calculation bias, with eigenvalue bias below 200 pcm, heterogeneous pin power bias under 3.1%, and homogeneous assembly power bias under 1.4%.</div></div>","PeriodicalId":8006,"journal":{"name":"Annals of Nuclear Energy","volume":"226 ","pages":"Article 111853"},"PeriodicalIF":2.3,"publicationDate":"2025-09-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145099858","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Honghang Chi, Yongshi Li, Yuchen Xie, Yahui Wang, Yu Ma
{"title":"General intrusive neutron transport reduced-order model applicable to various deterministic methods","authors":"Honghang Chi, Yongshi Li, Yuchen Xie, Yahui Wang, Yu Ma","doi":"10.1016/j.anucene.2025.111854","DOIUrl":"10.1016/j.anucene.2025.111854","url":null,"abstract":"<div><div>The intrusive reduced-order model (IROM) has been developed for detailed neutron transport solving to reduce its expansive computational costs. These techniques take neutron angular fluxes with consistent spatial angular discretization format as snapshots, which limits the generality of existing computing software, especially commercial software. In this work, a general IROM (GIROM) for neutron transport equation (NTE) is proposed. The neutron scalar fluxes from any deterministic method can be directly used as snapshots, which effectively improves the universality of neutron transport IROM. Based on the angular integrated NTE, the corresponding neutron convection term can be calculated by neutron scalar flux. Both the neutron scalar flux and convection term are taken as snapshots to eliminate the dependency on snapshot discretize formats. Numerical results show that the GIROM is compatible with the neutron scalar flux of various deterministic neutron transport solvers, and the computational efficiency can be effectively improved.</div></div>","PeriodicalId":8006,"journal":{"name":"Annals of Nuclear Energy","volume":"226 ","pages":"Article 111854"},"PeriodicalIF":2.3,"publicationDate":"2025-09-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145099857","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Feiyang Li , Youwei Zeng , Pengcheng Zhao , Zijing Liu , Wei Li
{"title":"Reliability analysis of a lead-bismuth cooled passive system based on AL-I surrogate model","authors":"Feiyang Li , Youwei Zeng , Pengcheng Zhao , Zijing Liu , Wei Li","doi":"10.1016/j.anucene.2025.111875","DOIUrl":"10.1016/j.anucene.2025.111875","url":null,"abstract":"<div><div>Passive residual heat removal systems ensure the safe operation of lead–bismuth fast reactors. However, the resistance of such systems is similar to natural driving forces, while small fluctuations in the surrounding environment and material parameters can cause system failure; thus, analyzing the reliability of passive residual heat removal systems is important for lead–bismuth cooling. This study utilizes the passive system in the lead–bismuth eutectic loop of the TALL-3D experimental facility and proposes a reliability analysis based on the active learning-integration (AL-I) surrogate model. The AL-I surrogate model is constructed first, and single-failure and multiple-failure region validations are performed to ensure accuracy and robustness of the model. Subsequently, the sensitivity and reliability of the TALL-3D non-energetic system is determined. The active learning ensemble surrogate model only needs 99 low-cost numerical calculations to obtain a reliable result with a failure rate of 0.0650%. This model not only significantly reduces the computational resources and time costs, but also allows high-precision failure probability assessments. Therefore, this study shows that the AL-I surrogate model is advantageous for lead–bismuth cooled non-energetic waste heat discharge systems and offers solid technical support for engineering such systems.</div></div>","PeriodicalId":8006,"journal":{"name":"Annals of Nuclear Energy","volume":"226 ","pages":"Article 111875"},"PeriodicalIF":2.3,"publicationDate":"2025-09-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145099859","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Zhaohao Wang , Jinkai Zhang , Tianliang Hu , Di Yun , Duoyu Jiang , Da Li , Lixin Chen , Wenbo Liu
{"title":"FCMI analysis of multiphysics coupling of fuel rods at the end of cycle of Xi’an pulsed reactor","authors":"Zhaohao Wang , Jinkai Zhang , Tianliang Hu , Di Yun , Duoyu Jiang , Da Li , Lixin Chen , Wenbo Liu","doi":"10.1016/j.anucene.2025.111867","DOIUrl":"10.1016/j.anucene.2025.111867","url":null,"abstract":"<div><div>The Xi’an Pulse Reactor (XAPR), being a small-scale experimental reactor, operates within a very complex environment. Consequently, accurate estimation of neutron-thermo-mechanical coupling of the fuel is a crucial step to prevent operation, particularly pulse operation, outside fuel thermo-mechanical safety margin. To obtain precise temperature and stress distributions within the XAPR, a loosely coupled approach based on the finite element method was employed. In the present work, neutron transport equations were calculated using OpenMC, and the Multiphysics Object-Oriented Simulation Environment (MOOSE), an open-source multi-physics coupling platform, was employed for heat transfer and mechanics calculations. The physical modules have also been verified by comparing calculation results against corresponding experimental data, confirming the reliability of the performed calculations. In addition, we conducted 3-dimensional (3-D) calculations of the bending effects of fuel rods due to Pellet-Clad Mechanical Interaction (PCMI) under different operating conditions and discussed the influence of different pulse state initiation time on the bending of the fuel.</div></div>","PeriodicalId":8006,"journal":{"name":"Annals of Nuclear Energy","volume":"226 ","pages":"Article 111867"},"PeriodicalIF":2.3,"publicationDate":"2025-09-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145099861","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Sixian Zhu , Jien Ma , Zheng Wang , Keming Bi , Xin Yan , Youtong Fang
{"title":"Modeling and analysis of electromagnetic, fluid and energy characteristics in flat linear induction pump","authors":"Sixian Zhu , Jien Ma , Zheng Wang , Keming Bi , Xin Yan , Youtong Fang","doi":"10.1016/j.anucene.2025.111841","DOIUrl":"10.1016/j.anucene.2025.111841","url":null,"abstract":"<div><div>Conductive liquid metals and their magnetohydrodynamic (MHD) systems are crucial in nuclear energy. Flat linear induction pumps (FLIPs) enhance safety and reliability working as drive components, and the magnetic-fluid coupling behavior within them deserves further study. This study develops a 1D analytical model for rapid performance estimation, incorporating end effects and half-filled slots. In finite element analysis (FEA), 2D/3D decoupled models and 3D fully-coupled MHD models are presented. Experimental results from a small-scale FLIP prototype are compared with the models. Results show that flow structure minimally affects the magnetic field, while electromagnetic forces significantly influence the liquid behavior. At low flow rates, reverse flow occurs locally. These findings advance understanding of FLIP’s coupled physics and energy conversion, aiding optimization for engineering applications.</div></div>","PeriodicalId":8006,"journal":{"name":"Annals of Nuclear Energy","volume":"226 ","pages":"Article 111841"},"PeriodicalIF":2.3,"publicationDate":"2025-09-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145099860","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"An implicit numerical method for the three-dimensional two-fluid model with large phase changes in pressurized water reactors","authors":"Na Young Song, Han Young Yoon","doi":"10.1016/j.anucene.2025.111873","DOIUrl":"10.1016/j.anucene.2025.111873","url":null,"abstract":"<div><div>An implicit numerical method is proposed for solving the three-dimensional two-fluid model with a large phase, which is widely used in two-phase flow analysis of Pressurized Water Reactors (PWRs). The phase change process is implicitly modeled by coupling it with the pressure equation derived from the momentum and mass conservation equations. Simultaneously, the convection and diffusion terms of the governing equations are treated implicitly. To reduce the size of the system matrix, the computation procedure is divided into two separate steps: the “phase-link” step and the “space-link” step. This implicit scheme has been implemented in CUPID, a three-dimensional two-phase flow analysis code developed specifically for PWR applications. The present method was verified using a conceptual problem that simulates two-phase flow blowdown and refill phenomena occurring during a Loss of Coolant Accident (LOCA). The calculations demonstrated stability even with a large Courant–Friedrichs–Lewy (CFL) number. Subsequently, the implicit scheme was validated against the FLECHT-SEASET reflood experiment. The robustness and accuracy of the implicit scheme were confirmed by comparing the simulation results with experimental data across different CFL numbers. The differences between the results obtained at various CFL numbers were negligible.</div></div>","PeriodicalId":8006,"journal":{"name":"Annals of Nuclear Energy","volume":"226 ","pages":"Article 111873"},"PeriodicalIF":2.3,"publicationDate":"2025-09-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145099850","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}