Development and verification of MOC code based on Krylov subspace and efficient preconditioning techniques

Abstract Derived from the conventional characteristic sweeping operations, the method of characteristics (MOC) with matrix form has more favorable performance compared with its traditional implementation potentially. However, the advantage depends heavily upon the efficiency of the linear algebraic solver. In current study, a simple and efficient preconditioning technique is implemented in the restart version of Generalized Minimal RESidual algorithm (GMRES) to accelerate the resulted linear system. A complete code including geometry processing and algebraic solver has been developed and verified with the reference benchmark problem. Numerical results demonstrate the code can model the reference problem accurately, and the proposed preconditioning techniques based on the typical iterative method can decrease dramatically the number of iterations without introducing additional computation and storage expense.

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

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

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

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

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

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

[7]  Elizabeth R. Jessup,et al.  A Technique for Accelerating the Convergence of Restarted GMRES , 2005, SIAM J. Matrix Anal. Appl..

[8]  Satoshi Takeda,et al.  Development of efficient Krylov preconditioning techniques for multi-dimensional method of characteristics , 2012 .

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

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

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

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

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

[14]  Elmer E Lewis,et al.  Benchmark on deterministic 3-D MOX fuel assembly transport calculations without spatial homogenization , 2004 .

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