In this paper,we propose a numerical calculation model of the multigroup neutron diffusion equation in 3D hexagonal geometry using the nodal Green's function method and verified it.We obtained one-dimensional tran...In this paper,we propose a numerical calculation model of the multigroup neutron diffusion equation in 3D hexagonal geometry using the nodal Green's function method and verified it.We obtained one-dimensional transverse integrated equations using the transverse integration procedure over 3D hexagonal geometry and denoted the solutions as a nodal Green's functions under the Neumann boundary condition.By applying a quadratic polynomial expansion of the transverse-averaged quantities,we derived the net neutron current coupling equation,equation for the expansion coefficients of the transverse-averaged neutron flux,and formulas for the coefficient matrix of these equations.We formulated the closed system of equations in correspondence with the boundary conditions.The proposed model was tested by comparing it with the benchmark for the VVER-440 reactor,and the numerical results were in good agreement with the reference solutions.展开更多
Higher-order modes of the neutron diffusion/transport equation can be used to study the temporal behavior of nuclear reactors and can be applied in modal analysis, transient analysis, and online monitoring of the reac...Higher-order modes of the neutron diffusion/transport equation can be used to study the temporal behavior of nuclear reactors and can be applied in modal analysis, transient analysis, and online monitoring of the reactor core. Both the deterministic method and the Monte Carlo(MC) method can be used to solve the higher-order modes. However, MC method, compared to the deterministic method, faces challenges in terms of computational efficiency and α mode calculation stability, whereas the deterministic method encounters issues arising from homogenization-related geometric and energy spectra adaptation.Based on the higher-order mode diffusion calculation code HARMONY, we developed a new higher-order mode calculation code, HARMONY2.0, which retains the functionality of computing λ and α higher-order modes from HARMONY1.0, but enhances the ability to treat complex geometries and arbitrary energy spectra using the MC-deterministic hybrid two-step strategy. In HARMONY2.0, the mesh homogenized multigroup constants were obtained using OpenMC in the first step,and higher-order modes were then calculated with the mesh homogenized core diffusion model using the implicitly restarted Arnoldi method(IRAM), which was also adopted in the HARMONY1.0 code. In addition, to improve the calculation efficiency, particularly in large higher-order modes, event-driven parallelization/domain decomposition methods are embedded in the HARMONY2.0 code to accelerate the inner iteration of λ∕α mode using OpenMP. Furthermore, the higher-order modes of complex geometric models, such as Hoogenboom and ATR reactors for λ mode and the MUSE-4 experiment facility for the prompt α mode, were computed using diffusion theory.展开更多
In this electronic article we use the one-dimensional multigroup neutron diffusion equation to reconstruct the neutron flux in a slab reactor from the nuclear parameters of the reactor, boundary and symmetry condition...In this electronic article we use the one-dimensional multigroup neutron diffusion equation to reconstruct the neutron flux in a slab reactor from the nuclear parameters of the reactor, boundary and symmetry condition, initial flux and?keff. The diffusion equation was solved analytically for one single homogeneous fuel region and for two regions considering fuel and reflector. To validate the method proposed, the results obtained in this article were compared using reference methods found in the literature.展开更多
Machine learning-based modeling of reactor physics problems has attracted increasing interest in recent years.Despite some progress in one-dimensional problems,there is still a paucity of benchmark studies that are ea...Machine learning-based modeling of reactor physics problems has attracted increasing interest in recent years.Despite some progress in one-dimensional problems,there is still a paucity of benchmark studies that are easy to solve using traditional numerical methods albeit still challenging using neural networks for a wide range of practical problems.We present two networks,namely the Generalized Inverse Power Method Neural Network(GIPMNN)and Physics-Constrained GIPMNN(PC-GIPIMNN)to solve K-eigenvalue problems in neutron diffusion theory.GIPMNN follows the main idea of the inverse power method and determines the lowest eigenvalue using an iterative method.The PC-GIPMNN additionally enforces conservative interface conditions for the neutron flux.Meanwhile,Deep Ritz Method(DRM)directly solves the smallest eigenvalue by minimizing the eigenvalue in Rayleigh quotient form.A comprehensive study was conducted using GIPMNN,PC-GIPMNN,and DRM to solve problems of complex spatial geometry with variant material domains from the fleld of nuclear reactor physics.The methods were compared with the standard flnite element method.The applicability and accuracy of the methods are reported and indicate that PC-GIPMNN outperforms GIPMNN and DRM.展开更多
In this study, based differential equations methods are used to solve equations because these methods are dependent on boundary value data more than other mathematical equations. We have calculated neutron flux, criti...In this study, based differential equations methods are used to solve equations because these methods are dependent on boundary value data more than other mathematical equations. We have calculated neutron flux, criticality and geometrical eigenvalue by using multi-group method and solving the neutron diffusion equation for finite and infinite cylindrical and spherical reactors in this study. For the calculation of the total neutron flux cross sections, we need the neutron diffusion equation. Thus, we have established the relationship between neuron flow and cross-section of neuron depending on neutron energy. Critical calculations have been made by comparing the results with MNCP (montecarlo n-partical) simulation methods. For necessary computer calculations, the programme, Wolfram-Matematica-7 has been used.展开更多
In this work,two approaches,based on the certified Reduced Basis method,have been developed for simulating the movement of nuclear reactor control rods,in time-dependent non-coercive settings featuring a 3D geometrica...In this work,two approaches,based on the certified Reduced Basis method,have been developed for simulating the movement of nuclear reactor control rods,in time-dependent non-coercive settings featuring a 3D geometrical framework.In particular,in a first approach,a piece-wise affine transformation based on subdomains division has been implemented for modelling the movement of one control rod.In the second approach,a“staircase”strategy has been adopted for simulating themovement of all the three rods featured by the nuclear reactor chosen as case study.The neutron kinetics has been modelled according to the so-called multi-group neutron diffusion,which,in the present case,is a set of ten coupled parametrized parabolic equations(two energy groups for the neutron flux,and eight for the precursors).Both the reduced order models,developed according to the two approaches,provided a very good accuracy comparedwith high-fidelity results,assumed as“truth”solutions.At the same time,the computational speed-up in the Online phase,with respect to the fine“truth”finite element discretization,achievable by both the proposed approaches is at least of three orders of magnitude,allowing a real-time simulation of the rod movement and control.展开更多
文摘In this paper,we propose a numerical calculation model of the multigroup neutron diffusion equation in 3D hexagonal geometry using the nodal Green's function method and verified it.We obtained one-dimensional transverse integrated equations using the transverse integration procedure over 3D hexagonal geometry and denoted the solutions as a nodal Green's functions under the Neumann boundary condition.By applying a quadratic polynomial expansion of the transverse-averaged quantities,we derived the net neutron current coupling equation,equation for the expansion coefficients of the transverse-averaged neutron flux,and formulas for the coefficient matrix of these equations.We formulated the closed system of equations in correspondence with the boundary conditions.The proposed model was tested by comparing it with the benchmark for the VVER-440 reactor,and the numerical results were in good agreement with the reference solutions.
基金supported by the National Natural Science Foundation of China(No.U2267207)Science and Technology on Reactor System Design Technology Laboratory(No.KFKT-05-FWHTWU-2023004).
文摘Higher-order modes of the neutron diffusion/transport equation can be used to study the temporal behavior of nuclear reactors and can be applied in modal analysis, transient analysis, and online monitoring of the reactor core. Both the deterministic method and the Monte Carlo(MC) method can be used to solve the higher-order modes. However, MC method, compared to the deterministic method, faces challenges in terms of computational efficiency and α mode calculation stability, whereas the deterministic method encounters issues arising from homogenization-related geometric and energy spectra adaptation.Based on the higher-order mode diffusion calculation code HARMONY, we developed a new higher-order mode calculation code, HARMONY2.0, which retains the functionality of computing λ and α higher-order modes from HARMONY1.0, but enhances the ability to treat complex geometries and arbitrary energy spectra using the MC-deterministic hybrid two-step strategy. In HARMONY2.0, the mesh homogenized multigroup constants were obtained using OpenMC in the first step,and higher-order modes were then calculated with the mesh homogenized core diffusion model using the implicitly restarted Arnoldi method(IRAM), which was also adopted in the HARMONY1.0 code. In addition, to improve the calculation efficiency, particularly in large higher-order modes, event-driven parallelization/domain decomposition methods are embedded in the HARMONY2.0 code to accelerate the inner iteration of λ∕α mode using OpenMP. Furthermore, the higher-order modes of complex geometric models, such as Hoogenboom and ATR reactors for λ mode and the MUSE-4 experiment facility for the prompt α mode, were computed using diffusion theory.
文摘In this electronic article we use the one-dimensional multigroup neutron diffusion equation to reconstruct the neutron flux in a slab reactor from the nuclear parameters of the reactor, boundary and symmetry condition, initial flux and?keff. The diffusion equation was solved analytically for one single homogeneous fuel region and for two regions considering fuel and reflector. To validate the method proposed, the results obtained in this article were compared using reference methods found in the literature.
基金partially supported by the National Natural Science Foundation of China(No.11971020)Natural Science Foundation of Shanghai(No.23ZR1429300)Innovation Funds of CNNC(Lingchuang Fund)。
文摘Machine learning-based modeling of reactor physics problems has attracted increasing interest in recent years.Despite some progress in one-dimensional problems,there is still a paucity of benchmark studies that are easy to solve using traditional numerical methods albeit still challenging using neural networks for a wide range of practical problems.We present two networks,namely the Generalized Inverse Power Method Neural Network(GIPMNN)and Physics-Constrained GIPMNN(PC-GIPIMNN)to solve K-eigenvalue problems in neutron diffusion theory.GIPMNN follows the main idea of the inverse power method and determines the lowest eigenvalue using an iterative method.The PC-GIPMNN additionally enforces conservative interface conditions for the neutron flux.Meanwhile,Deep Ritz Method(DRM)directly solves the smallest eigenvalue by minimizing the eigenvalue in Rayleigh quotient form.A comprehensive study was conducted using GIPMNN,PC-GIPMNN,and DRM to solve problems of complex spatial geometry with variant material domains from the fleld of nuclear reactor physics.The methods were compared with the standard flnite element method.The applicability and accuracy of the methods are reported and indicate that PC-GIPMNN outperforms GIPMNN and DRM.
文摘In this study, based differential equations methods are used to solve equations because these methods are dependent on boundary value data more than other mathematical equations. We have calculated neutron flux, criticality and geometrical eigenvalue by using multi-group method and solving the neutron diffusion equation for finite and infinite cylindrical and spherical reactors in this study. For the calculation of the total neutron flux cross sections, we need the neutron diffusion equation. Thus, we have established the relationship between neuron flow and cross-section of neuron depending on neutron energy. Critical calculations have been made by comparing the results with MNCP (montecarlo n-partical) simulation methods. For necessary computer calculations, the programme, Wolfram-Matematica-7 has been used.
基金We acknowledge CINECA and Regione Lombardia LISA computational initiative,for the availability of high performance computing resources and support.G.Rozza acknowledges INDAM-GNCS national activity group and NOFYSAS program of SISSA.
文摘In this work,two approaches,based on the certified Reduced Basis method,have been developed for simulating the movement of nuclear reactor control rods,in time-dependent non-coercive settings featuring a 3D geometrical framework.In particular,in a first approach,a piece-wise affine transformation based on subdomains division has been implemented for modelling the movement of one control rod.In the second approach,a“staircase”strategy has been adopted for simulating themovement of all the three rods featured by the nuclear reactor chosen as case study.The neutron kinetics has been modelled according to the so-called multi-group neutron diffusion,which,in the present case,is a set of ten coupled parametrized parabolic equations(two energy groups for the neutron flux,and eight for the precursors).Both the reduced order models,developed according to the two approaches,provided a very good accuracy comparedwith high-fidelity results,assumed as“truth”solutions.At the same time,the computational speed-up in the Online phase,with respect to the fine“truth”finite element discretization,achievable by both the proposed approaches is at least of three orders of magnitude,allowing a real-time simulation of the rod movement and control.
