Abstract:In this research, an innovative numerical method is presented to solve the fractional space-time multigroup neutron diffusion model. The method employs the fractional-order Bernstein polynomials for temporal calculations and higher order finite difference schemes for spatial discretization. Two major drawbacks should be addressed to solve such fractional diffusion model. The first is that the spatial discretization produces a large system of ordinary differential equations with a fundamental square matrix of huge size. To address this problem, the used numerical method should be simple and fast to overcome such difficulty. Therefore, a straightforward recurrence relation for the Bernstein coefficient estimations is introduced. The second is that the fractional operator consumes more computational times due to its persistent memory, especially when using smaller values of the anomalous diffusion exponents. This problem is overcome by employing the technique of automatic time step size control to reduce the number of iterations and CPU times. Also, without causing significant change in the system complexity, a higher order central finite difference scheme with accuracy order equals O(Δ4) is employed for spatial discretization of the core interior mesh points. To further improve the stationary state results, a new scheme with accuracy order equals O(Δ3) is introduced using seven point stencils at the core boundaries. It is proved that using these higher order schemes at both the interior and boundary mesh points significantly improves the accuracy of the steady state calculations. The method is applied to homogenous and heterogonous benchmark reactors with varying anomalous subdiffusion exponents. Numerical simulation of assembly power densities and temperature distributions for LRA BWR reactor including adiabatic heating and Doppler feedback is presented. It is proved that for heterogeneous reactors, the effects of the fractional orders diminish throughout the transient for step, ramp perturbations and also with the presence of thermal-hydraulic temperature feedback.

Hermite methods for multi-energy group neutron kinetics model via quasi-static approach and piecewise functions

Simultaneous Approach for Simulation of A High-Temperature Gas-Cooled Reactor

Neutron transport methods for multiphysics heterogeneous reactor core simulation in Griffin

A Simulation Method for Nuclear Power Plant with High-Temperature Gas-Cooled Reactor

Eigenvalue Solvers for Modeling Nuclear Reactors on Leadership Class Machines

Reduced-order methods for neutron transport kinetics problem based on proper orthogonal decomposition and dynamic mode decomposition

Neutron Generation Based Method for Monte Carlo Three-Dimensional Reactor Time-Dependent Simulation

Dynamic Monte Carlo Method for Nuclear Reactor Kinetics Calculations

Fractional Bernstein polynomial method for solving time-fractional neutron diffusion systems

Numerical solutions for nonlinear mixed integro-differential equation and the neutron transport equation

An efficient solution of the multi-term multi-harmonic electron Boltzmann equation for use in global models

Application of Improved Quasi-Static Method to High Temperature Gas Cooled Reactor Simulator

Monte-Carlo Methods for the Neutron Transport Equation

Analysis of the 3D Neutron Kinetics for MHTGR

Method of Characteristics for 3D, Full-Core Neutron Transport on Unstructured Mesh

Development of a neutron space–time kinetics solver with improved quasi-static method based on OpenFOAM

A novel discontinuous-Galerkin deterministic neutronics model for Fusion applications: development and benchmarking

A Consistent One-Dimensional Multigroup Diffusion Model for Molten Salt Reactor Neutronics Calculations

Application of Monte Carlo Method to Solve the Neutron Kinetics Equation for a Subcritical Assembly

Multiphysics for nuclear energy applications using a cohesive computational framework

Adoption of a fast dynamic model to simulate transient operations of MNSR