DLR-Logo -> http://www.dlr.de
DLR Portal Home | Imprint | Privacy Policy | Contact | Deutsch
Fontsize: [-] Text [+]

cuSwift - A suite of numerical integration methods for modelling planetary systems implemented in C/CUDA

Hellmich, S. and Mottola, S. and Hahn, G. and Kührt, E. and Hlawitschka, M. (2014) cuSwift - A suite of numerical integration methods for modelling planetary systems implemented in C/CUDA. Asteroids, Comets, Meteors, June 30th – July 4th 2014, Helsinki.

Full text not available from this repository.


Simulations of dynamical processes in planetary systems represent an important tool for studying their orbital evolution. Using modern numerical integration methods, it is possible to model systems containing many thousands of objects over time scales of several hundred million years. However, in general supercomputers are needed to get reasonable simulation results in acceptable execution times. To exploit the ever growing computation power of Graphics Processing Units (GPUs) in modern desktop computers we implemented cuSwift, a library of numerical integration methods for studying long- term dynamical processes in planetary systems. cuSwift can be seen as a re-implementation of the famous SWIFT integrator package written by Hal Levison and Martin Duncan. cuSwift is written in C/CUDA and contains different integration methods for various purposes. So far, we have implemented three algorithms: a 15th order Radau integrator, the Wisdom-Holman Mapping (WHM) integrator and the Regularized Mixed Variable Symplectic (RMVS) Method. These algorithms treat only the planets as mutually gravitationally interacting bodies whereas asteroids and comets (or other minor bodies of interest) are treated as massless test particles which are gravitationally influenced by the massive bodies but do not affect each other or the massive bodies. The main focus of this work is on the symplectic methods (WHM and RMVS) which use a larger time step and thus are capable of integrating many particles over a large time span. As an additional feature, we implemented the non-gravitational Yarkovsky effect as described by M. Brož. With cuSwift we show that the use of modern GPUs makes it possible to speed up these methods by more than one order of magnitude compared to the single-core CPU implementation, thereby enabling modest workstation computers to perform long-term dynamical simulations. We use these methods to study the influence of the Yarkovsky effect on resonant asteroids. We present first results and compare them with integrations done with the original algorithms implemented in SWIFT in order to assess the numerical precision of cuSwift and to demonstrate the speedup we achieved using the GPU.

Item URL in elib:https://elib.dlr.de/90955/
Document Type:Conference or Workshop Item (Speech)
Title:cuSwift - A suite of numerical integration methods for modelling planetary systems implemented in C/CUDA
AuthorsInstitution or Email of AuthorsAuthor's ORCID iDORCID Put Code
Hlawitschka, M.Department of computer science, Leipzig University, Augustusplatz 10, 04109 LeipzigUNSPECIFIEDUNSPECIFIED
Date:July 2014
Refereed publication:No
Open Access:No
Gold Open Access:No
In ISI Web of Science:No
Keywords:numerical integration methods, Solar system dynamics, GPGPU, Jupiter Trojans
Event Title:Asteroids, Comets, Meteors
Event Location:Helsinki
Event Type:international Conference
Event Dates:June 30th – July 4th 2014
Organizer:University of Helsinki
HGF - Research field:Aeronautics, Space and Transport
HGF - Program:Space
HGF - Program Themes:Space Exploration
DLR - Research area:Raumfahrt
DLR - Program:R EW - Space Exploration
DLR - Research theme (Project):R - Exploration of the Solar System
Location: Berlin-Adlershof
Institutes and Institutions:Institute of Planetary Research
Institute of Planetary Research > Asteroids and Comets
Deposited By: Hellmich, Stephan
Deposited On:19 Nov 2014 07:29
Last Modified:19 Nov 2014 07:29

Repository Staff Only: item control page

Help & Contact
electronic library is running on EPrints 3.3.12
Website and database design: Copyright © German Aerospace Center (DLR). All rights reserved.