Simple and Efficient Parallelization Method for MOC Calculation

A parallelization method that does not entail much data communication and allowseasy implementation was developed for the method of characteristics (MOC). In the parallelization method, azimuthal angles were grouped to compute, on the same processor, theangular flux before and after the reflection on the outer boundary. The method can be easily applied to existing MOC codes without the highly technical knowledge on efficientdata communication in parallel computing. It was implemented into the GALAXY code and numerical comparisons were performed. As a result, computation speedup by a factor of 9 to 10 and good parallel efficiencies of 70 to 80% were achieved by using twelve processors. The speedup benefits the practical reactor core design using MOC codes. Other thanthe parallel efficiency, the proposed method allows easy implementation because no drastic change in the conventional algorithm for existing sequential process codes is needed. This is carried out with the absence of data communication between processors during the inner iteration of MOC with the proposed method.

[1]  Ser Gi Hong,et al.  Acceleration and Parallelization of the Method of Characteristics for Lattice and Whole-Core Heterogeneous Calculations , 2000 .

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

[3]  Anthony J. Baratta,et al.  KRAM, A lattice physics code for modeling the detailed depletion of gadolinia isotopes in BWR lattice designs , 1990 .

[4]  K. D. Lathrop,et al.  TWOTRAN-II: an interfaced, exportable version of the TWOTRAN code for two- dimensional transport , 1973 .

[5]  D. Knott,et al.  Validation of the CASMO-4 transport solution , 1993 .

[6]  日本原子力研究所,et al.  MVP/GMVP II : General purpose Monte Carlo codes for neutron and photon transport calculations based on continuous energy and multigroup methods , 2005 .

[7]  Akio Yamamoto,et al.  Neutron Transport Models of AEGIS: An Advanced Next-Generation Neutronics Design System , 2007 .

[8]  K. S. Smith Nodal method storage reduction by nonlinear iteration , 1983 .

[9]  A. Leonard,et al.  The method of characteristics in general geometry , 1995 .

[10]  A. Yamamoto,et al.  Derivation of Optimum Polar Angle Quadrature Set for the Method of Characteristics Based on Approximation Error for the Bickley Function , 2007 .

[11]  Toshikazu Takeda,et al.  Transport Calculations of MOX and UO2 Pin Cells by the Method of Characteristics , 1998 .

[12]  Robert B. Ross,et al.  Using MPI-2: Advanced Features of the Message Passing Interface , 2003, CLUSTER.

[13]  Akio Yamamoto,et al.  Development of the new pin-by-pin core calculation method with embedded heterogeneous assembly calculation , 2006 .

[14]  Mohamed Dahmani,et al.  Parallel Solver Based on the Three-Dimensional Characteristics Method: Design and Performance Analysis , 2005 .

[15]  Yasunori Kitamura,et al.  Convergence Improvement of Coarse Mesh Rebalance Method for Neutron Transport Calculations , 2004 .

[16]  Jasmina L. Vujic,et al.  ANEMONA — a neutron transport code for general geometry reactor assemblies based on the method of characteristics and R-function solid modeler , 2001 .

[17]  K. Yamaji Development of a new lattice physics code GALAXY for flexible geometry representation in next generation core analysis system , 2007 .