The discrete ordinates method is well-suited to implementation with block-structured adaptive mesh refinement (AMR). AMR meshes group points into logically-rectangular patches, and provide the benefits of adaptivity without sacrificing the efficiency and geometric regularity of regular grids. In particular, these meshes preserve the directional ordering of points required for explicit ordinate transport sweeps. A number of algorithmic issues must be addressed to make such a method practical. These include sequencing of ordinates and grids for parallel execution, simultaneous solution of the transport equation on multiple levels of the grid hierarchy, implicit coupling to the fluid energy, and conservation of energy in a time-dependent context where different grid levels advance with different timesteps. The author discusses these and other issues and present example calculations in two and three spatial dimensions.

Advanced numerical simulations that use adaptive mesh refinement (AMR) methods have now become routine in engineering and science. Originally developed for computational fluid dynamics applications these methods have propagated to fields as diverse as astrophysics, climate modeling, combustion, biophysics and many others. The underlying physical models and equations used in these disciplines are rather different, yet algorithmic and implementation issues facing practitioners are often remarkably similar. Unfortunately, there has been little effort to review the advances and outstanding issues of adaptive mesh refinement methods across such a variety of fields. This book attempts to bridge this gap. The book presents a collection of papers by experts in the field of AMR who analyze past advances in the field and evaluate the current state of adaptive mesh refinement methods in scientific computing.

Block-structured adaptive mesh refinement (AMR) uses a mesh structure built up out of locally-uniform rectangular grids. In the BoxLib parallel framework used by the Raptor code, each processor operates on one or more of these grids at each refinement level. The decomposition of the mesh into grids and the distribution of these grids among processors may change every few timesteps as a calculation proceeds. Finer grids use smaller timesteps than coarser grids, requiring additional work to keep the system synchronized and ensure conservation between different refinement levels. In a paper for NECDC 2002 I presented preliminary results on implementation of parallel transport sweeps on the AMR mesh, conjugate gradient acceleration, accuracy of the AMR solution, and scalar speedup of the AMR algorithm compared to a uniform fully-refined mesh. This paper continues with a more in-depth examination of the parallel scaling properties of the scheme, both in single-level and multi-level calculations. Both sweeping and setup costs are considered. The algorithm scales with acceptable performance to several hundred processors. Trends suggest, however, that this is the limit for efficient calculations with traditional transport sweeps, and that modifications to the sweep algorithm will be increasingly needed as job sizes in the thousands of processors become common.

The discrete ordinates form of the radiative transport equation (RTE) is spatially discretized and solved using an adaptive mesh refinement (AMR) algorithm. This technique permits the local grid refinement to minimize spatial discretization error of the RTE. An error estimator is applied to define regions for local grid refinement; overlapping refined grids are recursively placed in these regions; and the RTE is then solved over the entire domain. The procedure continues until the spatial discretization error has been reduced to a sufficient level. The following aspects of the algorithm are discussed: error estimation, grid generation, communication between refined levels, and solution sequencing. This initial formulation employs the step scheme, and is valid for absorbing and isotopically scattering media in two-dimensional enclosures. The utility of the algorithm is tested by comparing the convergence characteristics and accuracy to those of the standard single-grid algorithm for several benchmark cases. The AMR algorithm provides a reduction in memory requirements and maintains the convergence characteristics of the standard single-grid algorithm; however, the cases illustrate that efficiency gains of the AMR algorithm will not be fully realized until three-dimensional geometries are considered.

Over recent years, Adaptive Mesh Refinement (AMR) algorithms which dynamically match the local resolution of the computational grid to the numerical solution being sought have emerged as powerful tools for solving problems that contain disparate length and time scales. In particular, several workers have demonstrated the effectiveness of employing an adaptive, block-structured hierarchical grid system for simulations of complex shock wave phenomena. Unfortunately, from the parallel algorithm developer's viewpoint, this class of scheme is quite involved; these schemes cannot be distilled down to a small kernel upon which various parallelizing strategies may be tested. However, because of their block-structured nature such schemes are inherently parallel, so all is not lost. In this paper we describe the method by which Quirk's AMR algorithm has been parallelized. This method is built upon just a few simple message passing routines and so it may be implemented across a broad class of MIMD machines. Moreover, the method of parallelization is such that the original serial code is left virtually intact, and so we are left with just a single product to support. The importance of this fact should not be underestimated given the size and complexity of the original algorithm. Quirk, James J. and Hanebutte, Ulf R. Unspecified Center NAS1-19480; RTOP 505-90-52-01...