Abstract:Coupled Monte Carlo (MC) and thermal-hydraulic analysis is valuable as a design or reference tool but can be slow, especially when implemented in a Picard iteration. Previous work has developed a novel prediction block to achieve convergence with fewer MC simulations. The prediction works in two stages: (1) a surrogate-like model predicts macroscopic cross sections on the fly and (2) a reduced-order neutronic model predicts the flux response to the updated cross sections. The main challenge with the prediction block is that the reduced-order neutronic model cannot reproduce the spatial flux distribution with high fidelity. This paper investigates the well-established Jacobian-free Newton Krylov (JFNK) method to preserve equivalency between a homogeneous (nodal diffusion) solution and a high-fidelity transport (MC) solution. Instead of performing multiple computationally consuming MC simulations, the nonlinear iterative approach iterates on correction parameters, e.g., assembly discontinuity factors (ADFs) or super homogenization (SPH) factors, using unexpensive nodal solutions. The JFNK approach does not require additional overhead from the MC solver to generate flux tallies. Further, the approach iterates on diffusion solutions produced directly from a desired code, thus ensuring that the parameters are compatible with that code. The approach is applied to correct a nodal diffusion solution for a realistic three-dimensional pressurized water reactor core. The results obtained in the paper show that the method is very successful in reproducing the heterogeneous solution (up to 2.5% difference in assembly flux for SPH and 0.3% for ADFs) without needing to modify the source code of the nodal diffusion solver. In addition, the results show that ADFs yield the best agreement and are also stable (i.e., weakly varying) when thermal-hydraulic fields are perturbed.

A Second Moment Method for k-Eigenvalue Acceleration with Continuous Diffusion and Discontinuous Transport Discretizations

A Variant of the Second-Moment Method for k-Eigenvalue Calculations

Consistent Second Moment Methods with Scalable Linear Solvers for Radiation Transport

High-Order Finite Element Second Moment Methods for Linear Transport

DG-IMEX Method for a Two-Moment Model for Radiation Transport in the $\mathcal{O}(v/c)$ Limit

A Discontinuous Galerkin Finite Element Method Based Axial SN for the 2D/1D Transport Method

Practical considerations for the adoption of Anderson acceleration in nonlinear diffusion accelerated transport

A novel fast iterative moment method for near-continuum flows

Local time stepping methods and discontinuous Galerkin methods applied to diffusion advection reaction equations

A new local and explicit kinetic method for linear and non-linear convection-diffusion problems with finite kinetic speeds: II. Multi-dimensional case

Efficient finite element method for 2D singularly perturbed parabolic convection diffusion problems with discontinuous source term

Convergence Study of Traditional 2D/1D Coupling Method for K-Eigenvalue Neutron Transport Problems with Fourier Analysis

Acceleration of Source Iteration using the Dynamic Mode Decomposition

Realizability-Preserving DG-IMEX Method for the Two-Moment Model of Fermion Transport

Anderson acceleration. Convergence analysis and applications to equilibrium chemistry

Scalable explicit implementation of anisotropic diffusion with Runge-Kutta-Legendre super-time-stepping

Multi-Level Solution Strategies for Implicit Time Discontinuous Galerkin Discretizations

The Development and Application of the Moment Method in the Gas Kinetic Theory

Poisson quadrature method of moments for 2D kinetic equations with velocity of constant magnitude

On the Use of the Jacobian-Free Newton Krylov Method to Generate One-Group Discontinuity and Super Homogenization Factors for Full-Core Neutron Diffusion Simulations

A Hybrid Discontinuous Galerkin Scheme for Multi-scale Kinetic Equations