We provide a novel and efficient algorithm for computing accelerations in the periodic large-N-body problem that is at the same time significantly faster and more accurate than previous methods. Our representation of the periodic acceleration is precisely mathematically equivalent to that determined by Ewald summation and is computed directly as an infinite lattice sum using the Newtonian kernel (|r|-1). Retaining this kernel implies that one can (i) extend the standard open boundary numerical algorithms and (ii) harness the tremendous computational speed possessed by Graphics Processing Units (GPUs) in computing Newtonian kernels straightforwardly to the periodic domain. The precise form of our direct interactions is based upon the adaptive softening length methodology introduced for open boundary conditions by Price and Monaghan. Furthermore, we describe a new Fast Multipole Method (FMM) that represents the multipoles and Taylor series as collections of pseudoparticles. Using these techniques we have computed forces to machine precision throughout the evolution of a 1 billion particle cosmological simulation with a price/performance ratio more than 100 times that of current numerical techniques operating at much lower accuracy.
|Advisor:||Pinto, Philip A.|
|Commitee:||Arnett, W. D., Dave, Romeel A., Eisenstein, Daniel J., Strittmatter, Peter A.|
|School:||The University of Arizona|
|School Location:||United States -- Arizona|
|Source:||DAI-B 70/10, Dissertation Abstracts International|
|Keywords:||Cosmology, Gravitational N-body problem, Periodic accelerations|
Copyright in each Dissertation and Thesis is retained by the author. All Rights Reserved
The supplemental file or files you are about to download were provided to ProQuest by the author as part of a
dissertation or thesis. The supplemental files are provided "AS IS" without warranty. ProQuest is not responsible for the
content, format or impact on the supplemental file(s) on our system. in some cases, the file type may be unknown or
may be a .exe file. We recommend caution as you open such files.
Copyright of the original materials contained in the supplemental file is retained by the author and your access to the
supplemental files is subject to the ProQuest Terms and Conditions of use.
Depending on the size of the file(s) you are downloading, the system may take some time to download them. Please be