甘比斯:一种新的GPU加速N体代码用于双星系统
Maximilian Zimmermann1, Elke Pilat-Lohinger1
1Department of Astrophysics, University of Vienna, Türkenschanzstraße 17, Vienna, 1180 Austria.
概括
我们开发了GANBISS,这是一个GPU加速的N体集成器,用于模拟二进制恒星系统中的行星小盘. 这个代码在NVIDIA GPU上提供高达100倍的加速度,用于复杂的天体物理模拟.
科学领域:
- 计算天体物理学 计算天体物理学
- 行星科学 行星科学
- 双星系统是双星系统.
背景情况:
- 在双星系统中模拟行星小盘的动态进化是计算密集的.
- 现有的方法可能无法有效处理大量的磁盘物体或复杂的引力相互作用.
研究的目的:
- 为了介绍GANBISS,一种新的GPU加速N体集成器,旨在模拟二进制恒星系统中的行星小盘.
- 评估这种新的集成方法的性能和保护性质.
主要方法:
- 使用CUDA C开发了GANBISS,实现了布利尔施-斯托尔集成方法.
- 利用NVIDIA GPU (计算能力≥3.5) 进行加速计算.
- 将GPU性能与传统的CPU计算进行比较.
主要成果:
- 甘比斯可以模拟数千个磁盘物体或多达5000万个不相互作用的无质体的系统.
- 证明了非简单集成器的能量和角动量保存特征.
- 与CPU计算相比,实现了高达100倍的GPU加速度,取决于对象数量.
结论:
- 甘比斯为双星系统动态中的N体模拟提供了显著的性能提升.
- 该代码适用于研究行星盘的演变和其他相关的天体物理问题.
- 突出了GPU加速在解决天体物理学中大规模N体问题的潜力.
相关概念视频
Reduced Mass Coordinates: Isolated Two-body Problem
1.4K
In classical mechanics, the two-body problem is one of the fundamental problems describing the motion of two interacting bodies under gravity or any other central force. When considering the motion of two bodies, one of the most important concepts is the reduced mass coordinates, a quantity that allows the two-body problem to be solved like a single-body problem. In these circumstances, it is assumed that a single body with reduced mass revolves around another body fixed in a position with an...
1.4K
Maxwell-Boltzmann Distribution: Problem Solving
1.6K
Individual molecules in a gas move in random directions, but a gas containing numerous molecules has a predictable distribution of molecular speeds, which is known as the Maxwell-Boltzmann distribution, f(v).
This distribution function f(v) is defined by saying that the expected number N (v1,v2) of particles with speeds between v1 and v2 is given by
This distribution function f(v) is defined by saying that the expected number N (v1,v2) of particles with speeds between v1 and v2 is given by
1.6K
Gravitation Between Spherically Symmetric Masses
952
The gravitational potential energy between two spherically symmetric bodies can be calculated from the masses and the distance between the bodies, assuming that the center of mass is concentrated at the respective centers of the bodies.
952
Second Order systems II
138
In an underdamped second-order system, where the damping ratio ζ is between 0 and 1, a unit-step input results in a transfer function that, when transformed using the inverse Laplace method, reveals the output response. The output exhibits a damped sinusoidal oscillation, and the difference between the input and output is termed the error signal. This error signal also demonstrates damped oscillatory behavior. Eventually, as the system reaches a steady state, the error diminishes to zero.
138
Gravity between Spherical Bodies
8.6K
Newton's law of gravitation describes the gravitational force between any two point masses. However, for extended spherical objects like the Earth, the Moon, and other planets, the law holds with an assumption that masses of spherical objects are concentrated at their respective centers.
This assumption can be proved easily by showing that the expression for gravitational potential energy between a hollow sphere of mass (M) and a point mass (m) is the same as it would be for a pair of extended...
This assumption can be proved easily by showing that the expression for gravitational potential energy between a hollow sphere of mass (M) and a point mass (m) is the same as it would be for a pair of extended...
8.6K
Equilibrium Conditions for a Particle
1.2K
When an object is in equilibrium, it is either at rest or moving with a constant velocity. There are two types of equilibrium: static and dynamic. Static equilibrium occurs when an object is at rest, while dynamic equilibrium occurs when an object is moving with a constant velocity. In both cases, there must be a balance of forces acting on the object.
To understand the concept of equilibrium, let us first consider the forces acting on an object. When different forces act on an object, they can...
To understand the concept of equilibrium, let us first consider the forces acting on an object. When different forces act on an object, they can...
1.2K


