Thesis (Ph.D)--Mechanical Engineering, Georgia Institute of Technology, 2010. / Committee Chair: Farzad Rahnema; Committee Co-Chair: Eric Elder; Committee Member: C.-K. Chris Wang; Committee Member: Rebecca Howell; Committee Member: Sang Cho. Part of the SMARTech Electronic Thesis and Dissertation Collection.

Radiation therapy has become a very import method for treating cancer patients. Thus, it is extremely important to accurately determine the location of energy deposition during these treatments, maximizing dose to the tumor region and minimizing it to healthy tissue. A Coarse-Mesh Transport Method (COMET) has been developed at the Georgia Institute of Technology in the Computational Reactor and Medical Physics Group for use very successfully with neutron transport to analyze whole-core criticality. COMET works by decomposing a large, heterogeneous system into a set of smaller fixed source problems. For each unique local problem that exists, a solution is obtained that we call a response function. These response functions are pre-computed and stored in a library for future use. The overall solution to the global problem can then be found by a linear superposition of these local problems. This method has now been extended to the transport of photons and electrons for use in medical physics problems to determine energy deposition from radiation therapy treatments. The main goal of this work was to develop benchmarks for testing in order to evaluate the COMET code to determine its strengths and weaknesses for these medical physics applications. For response function calculations, Legendre polynomial expansions are necessary for space, angle, polar angle, and azimuthal angle. An initial sensitivity study was done to determine the best orders for future testing. After the expansion orders were found, three simple benchmarks were tested: a water phantom, a simplified lung phantom, and a non-clinical slab phantom. Three more clinically relevant problems were developed from patient CT scans. Different coarse-mesh sizes and incident energies were tested. The COMET solutions for each case were compared to a reference solution obtained by pure Monte Carlo results from EGSnrc. In most cases, the COMET solutions produced reasonably good agreement with the COMET solutions. It was found that better results were obtained for lower energy incident photon beams as well as for larger mesh sizes. Recommendations were made for future development of COMET and the numerical benchmarks.

A Variational Transport Theory Method for Two-Dimensional Reactor Core Calculations Scott W. Mosher 110 Pages Directed by Dr. Farzad Rahnema It seems very likely that the next generation of reactor analysis methods will be based largely on neutron transport theory, at both the assembly and core levels. Signifi-cant progress has been made in recent years toward the goal of developing a transport method that is applicable to large, heterogeneous coarse-meshes. Unfortunately, the ma-jor obstacle hindering a more widespread application of transport theory to large-scale calculations is still the computational cost. In this dissertation, a variational heterogeneous coarse-mesh transport method has been extended from one to two-dimensional Cartesian geometry in a practical fashion. A generalization of the angular flux expansion within a coarse-mesh was developed. This allows a far more efficient class of response functions (or basis functions) to be employed within the framework of the original variational principle. New finite element equations were derived that can be used to compute the expansion coefficients for an individual coarse-mesh given the incident fluxes on the boundary. In addition, the non-variational method previously used to converge the expansion coefficients was developed in a new and more thorough manner by considering the implications of the fission source treat-ment imposed by the response expansion. The new coarse-mesh method was implemented for both one and two-dimensional (2-D) problems in the finite-difference, multigroup, discrete ordinates approximation. An efficient set of response functions was generated using orthogonal boundary conditions constructed from the discrete Legendre polynomials. Several one and two-dimensional heterogeneous light water reactor benchmark problems were studied. Relatively low-order response expansions were used to generate highly accurate results using both the variational and non-variational methods. The expansion order was found to have a far more significant impact on the accuracy of the results than the type of method. The varia-tional techniques provide better accuracy, but at substantially higher computational costs. The non-variational method is extremely robust and was shown to achieve accurate re-sults in the 2-D problems, as long as the expansion order was not very low.