Track truncation method for improving the computational efficiency of 2D matrix MOC

Abstract The method of characteristics (MOC) with matrix form has more favorable performance compared with its standard application. However, it is observed that the advantage decreases when the physical problem size becomes larger. The spatial domain decomposition method (DDM) is a better remedy to reduce computation burden for large scale problem, but increasing the outer iterations due to the decoupling between flat source regions. To overcome this drawback, the straightforward multi-domain coupling technique is studied and utilized in the current study based on a track truncation method (TTM). The resulting linear system is solved using preconditioned generalized minimal residual (PGMRES) method, a variant of Krylov subspace method. Numerical results for benchmark problems demonstrate the proposed acceleration technique is efficient in terms of memory usage and computation time without compromising accuracy for both hexagonal and rectangular geometries. For the large scale problem, it shows an ideal linear relationship between the memory/time required and the physical problem size rather than quadratic growth as before.

[1]  Han Gyu Joo,et al.  Methods and performance of a three-dimensional whole-core transport code DeCART , 2004 .

[2]  Ser Gi Hong,et al.  CRX: A code for rectangular and hexagonal lattices based on the method of characteristics , 1998 .

[3]  Thomas J. Downar,et al.  An optimally diffusive Coarse Mesh Finite Difference method to accelerate neutron transport calculations , 2016 .

[4]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[5]  C. T. McDaniel,et al.  Optimal polar angles and weights for the characteristics method , 1995 .

[6]  L. Cao,et al.  A 2D/1D coupling neutron transport method based on the matrix MOC and NEM methods , 2013 .

[7]  Simone Santandrea,et al.  Treatment of Boundary Conditions in Trajectory-Based Deterministic Transport Methods , 2002 .

[8]  Liangzhi Cao,et al.  An acceleration technique for 2D MOC based on Krylov subspace and domain decomposition methods , 2011 .

[9]  Jim E. Morel,et al.  Angular Multigrid Preconditioner for Krylov-Based Solution Techniques Applied to the Sn Equations with Highly Forward-Peaked Scattering , 2012 .

[10]  Han Gyu Joo,et al.  Whole Core Transport Calculation Employing Hexagonal Modular Ray Tracing and CMFD Formulation , 2008 .

[11]  Kord Smith,et al.  FULL-CORE, 2-D, LWR CORE CALCULATIONS WITH CASMO-4E , 2002 .

[12]  S. Kosaka,et al.  Transport Theory Calculation for a Heterogeneous Multi-Assembly Problem by Characteristics Method with Direct Neutron Path Linking Technique , 2000 .

[13]  T. Downar,et al.  Theory and analysis of accuracy for the method of characteristics direction probabilities with boundary averaging , 2015 .

[14]  Brendan Kochunas,et al.  Overview of development and design of MPACT: Michigan parallel characteristics transport code , 2013 .

[15]  Kan Wang,et al.  Matrix Method of Characteristics Based on Modular Ray Tracing , 2013 .

[16]  Minjun Peng,et al.  Development and verification of MOC code based on Krylov subspace and efficient preconditioning techniques , 2017 .

[17]  Chuntao Tang,et al.  Development and verification of an MOC code employing assembly modular ray tracing and efficient acceleration techniques , 2009 .

[18]  Ser Gi Hong,et al.  Method of characteristic direction probabilities for heterogeneous lattice calculation , 1999 .

[19]  Brendan Kochunas,et al.  Model and analysis of performance for the method of characteristics direction probabilities with boundary averaging , 2015 .

[20]  Benoit Forget,et al.  The OpenMOC method of characteristics neutral particle transport code , 2014 .

[21]  Farzad Rahnema,et al.  Simplified two and three dimensional HTTR benchmark problems , 2011 .

[22]  Generalized Coarse-Mesh Rebalance Method for Acceleration of Neutron Transport Calculations , 2005 .