A Block Adaptive Mesh Refinement Algorithm for the Neutral Particle Transport Equation

Abstract We describe the development and implementation of a block-based adaptive mesh refinement (AMR) algorithm for solving the discrete ordinates neutral particle transport equation. AMR algorithms allow mesh refinement in areas of interest without requiring the extension of this refinement throughout the entire problem geometry, minimizing the number of computational cells required for calculations. The block-based AMR algorithm described here is a hybrid between traditional cell or patch-based approaches and is designed to allow an efficient parallel solution of the transport equation while still reducing the cell count. This paper discusses the data structure implementation and CPU/memory efficiency for our Block AMR method, the equations and procedures used in mapping edge fluxes between blocks of different refinement levels for both diamond and linear discontinuous spatial discretizations, effects of AMR on mesh convergence, and our approach to parallelization. Comparisons between our Block AMR method and a traditional single-level mesh are presented for a sample brachytherapy problem. The Block AMR results are shown to be significantly faster for this problem (on at least a few processors), while still returning an accurate solution.

[1]  T. Hill.  ONETRAN: a discrete ordinates finite element code for the solution of the one-dimensional multigroup transport equation , 1975 .

[2]  R. E. Alcouffe,et al.  Diffusion synthetic acceleration methods for the diamond-differenced discrete-ordinates equations , 1977 .

[3]  R. D. O'Dell,et al.  Nodal methods for discrete-ordinates transport problems in (x,y) geometry , 1981 .

[4]  P. Colella,et al.  Local adaptive mesh refinement for shock hydrodynamics , 1989 .

[5]  An adaptive weighted diamond-differencing method for three-dimensional XYZ geometry , 1993 .

[6]  E. Lewis,et al.  Computational Methods of Neutron Transport , 1993 .

[7]  R. S. Baker,et al.  DANTSYS: A diffusion accelerated neutral particle transport code system , 1995 .

[8]  J. Quirk.  A parallel adaptive grid algorithm for computational shock hydrodynamics , 1996 .

[9]  Marvin L. Adams,et al.  A Transport Synthetic Acceleration method for transport iterations , 1997 .

[10]  P. Colella,et al.  An Adaptive Mesh Refinement Algorithm for the Radiative Transport Equation , 1998 .

[11]  R. Baker,et al.  An Sn algorithm for the massively parallel CM-200 computer , 1998 .

[12]  P. Colella,et al.  A CONSERVATIVE ADAPTIVE-MESH ALGORITHM FOR UNSTEADY, COMBINED-MODE HEAT TRANSFER USING THE DISCRETE ORDINATES METHOD , 1999 .

[13]  Scott A. Turner.  Automatic Mesh Coarsening for Discrete Ordinates Codes , 1999 .

[14]  D. Rogers,et al.  Two-Dimensional Discrete Ordinates Photon Transport Calculations for Brachytherapy Dosimetry Applications , 2000 .

[15]  Raymond E. Alcouffe,et al.  PARALLEL 3-D SN PERFORMANCE FOR DANTSYSMPI ON THE CRAY T 3 D , .