N-body simulation for self-gravitating collisional systems with a new SIMD instruction set extension to the x86 architecture, Advanced Vector eXtensions

We present a high-performance N-body code for self-gravitating collisional systems accelerated with the aid of a new SIMD instruction set extension of the x86 architecture: Advanced Vector eXtensions (AVX), an enhanced version of the Streaming SIMD Extensions (SSE). With one processor core of Intel Core i7-2600 processor (8 MB cache and 3.40 GHz) based on Sandy Bridge micro-architecture, we implemented a fourth-order Hermite scheme with individual timestep scheme (Makino and Aarseth, 1992), and achieved the performance of ∼20 giga floating point number operations per second (GFLOPS) for double-precision accuracy, which is two times and five times higher than that of the previously developed code implemented with the SSE instructions (Nitadori et al., 2006b), and that of a code implemented without any explicit use of SIMD instructions with the same processor core, respectively. We have parallelized the code by using so-called NINJA scheme (Nitadori et al., 2006a), and achieved ∼90 GFLOPS for a system containing more than N = 8192 particles with 8 MPI processes on four cores. We expect to achieve about 10 tera FLOPS (TFLOPS) for a self-gravitating collisional system with N ∼ 10⁵ on massively parallel systems with at most 800 cores with Sandy Bridge micro-architecture. This performance will be comparable to that of Graphic Processing Unit (GPU) cluster systems, such as the one with about 200 Tesla C1070 GPUs (Spurzem et al., 2010). This paper offers an alternative to collisional N-body simulations with GRAPEs and GPUs.     

Effect of a uniform magnetic field on nonlinear magnetocenvection in a rotating fluid spherical shell

Abstract

Dynamic interaction between magnetic field and fluid motion is studied through a numerical experiment of nonlinear three-dimensional magnetoconvection in a rapidly rotating spherical fluid shell to which a uniform magnetic field parallel to its spin axis is applied. The fluid shell is heated by internal heat sources to maintain thermal convection. The mean value of the magnetic Reynolds number in the fluid shell is 22.4 and 10 pairs of axially aligned vortex rolls are stably developed. We found that confinement of magnetic flux into anti-cyclonic vortex rolls was crucial on an abrupt change of the mode of magnetoconvection which occurred at Δ = 1 ～ 2, where A is the Elsasser number. After the mode change, the fluid shell can store a large amount of magnetic flux in itself by changing its convection style, and the magnetostrophic balance among the Coriolis, Lorentz and pressure forces is established. Furthermore, the toroidal/poloidal ratio of the induced magnetic energy becomes less than unity, and the magnetized anti-cyclones are enlarged due to the effect of the magnetic force. Using these key ideas, we investigated the causes of the mode change of magnetoconvection. Considering relatively large magnetic Reynolds number and a rapid rotation rate of this model, we believe that these basic ideas used to interpret the present numerical experiment can be applied to the dynamics in the Earth's and other planetary cores.     

Unexpected formation modes of the first hard binary in core collapse

The conventional wisdom for the formation of the first hard binary in core collapse is that three-body interactions of single stars form many soft binaries, most of which are quickly destroyed, but eventually one of them survives. We report on direct N-body simulations to test these ideas, for the first time. We find that the assumptions are incorrect in the majority of the cases: (1) quite a few three-body interactions produce a hard binary from scratch; (2) in many cases there are more than three bodies directly and simultaneously involved in the production of the first binary. The main reason for the discrepancies is that the core of a star cluster, at the first deep collapse, contains typically only five or so stars. Therefore, the homogeneous background assumption, which still would be reasonable for, say, 25 stars, utterly breaks down. There have been some speculations in this direction, but we demonstrate this result here explicitly, for the first time.     

Linear Magnetoconvection in Rotating Fluid Spheres Permeated by a Uniform Axial Magnetic Field

A linear analysis of thermally driven magnetoconvection is carried out with emphasis on its application to convection in the Earth's core. We consider a rotating and self-gravitating fluid sphere (or spherical shell) permeated by a uniform magnetic field parallel to the spin axis. In rapidly rotating cases, we find that five different convective modes appear as the uniform field is increased; namely, geostrophic, polar convective, magneto-geostrophic, fast magnetostrophic and slow magnetostrophic modes. The polar convective (P) and magneto-geostrophic (E) modes seem to be of geophysical interest. The P mode is characterized by such an axisymmetric meridional circulation that the fluid penetrates the equatorial plane, suggesting that generation of quadrapole from dipole fields could be explained by a linear process. The E mode is characterized by a few axially aligned columnar rolls which are almost two-dimensional due to a modified Proudman-Taylor theorem.