اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدAssessing the Hierarchical Hamiltonian Splitting Integrator for Collisionless N-body Simulations
59
0
0.0
(
0
)
نشر من قبل
Gabriela Aguilar-Arg\\\"uello Miss
تاريخ النشر
2020
مجال البحث
فيزياء
والبحث باللغة
English
تأليف
G. Aguilar-Arguello
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
The N-body problem has become one of the hottest topics in the fields of computational dynamics and cosmology. The large dynamical range in some astrophysical problems led to the use of adaptive time steps to integrate particle trajectories, however, the search of optimal strategies is still challenging. We quantify the performance of the hierarchical time step integrator Hamiltonian Splitting (HamSp) for collisionless multistep simulations. We compare with the constant step Leap-Frog (LeapF) integrator and the adaptive one (AKDK). Additionally, we explore the impact of different time step assigning functions. There is a computational overhead in HamSp however there are two interesting advantages: choosing a convenient time-step function may compensate and even turn around the efficiency compared with AKDK. We test both reversibility and time symmetry. The symmetrized nature of the HamSp integration is able to provide time-reversible integration for medium time scales and overall deliver better energy conservation for long integration times, and the linear and angular momentum are preserved at machine precision. We address the impact of using different integrators in astrophysical systems. We found that in most situations both AKDK and HamSp are able to correctly simulate the problems. We conclude that HamSp is an attractive and competitive alternative to AKDK, with, in some cases, faster and with better energy and momentum conservation. The use of recently discussed Bridge splitting techniques with HamSp may allow to reach considerably high efficiency.
قيم البحث
اقرأ أيضاً
We present a general scheme for constructing Monte Carlo realizations of equilibrium, collisionless galaxy models with known distribution function (DF) f_0. Our method uses importance sampling to find the sampling DF f_s that minimizes the mean-squar
e formal errors in a given set of projections of the DF f_0. The result is a multi-mass N-body realization of the galaxy model in which ``interesting regions of phase-space are densely populated by lots of low-mass particles, increasing the effective N there, and less interesting regions by fewer, higher-mass particles. As a simple application, we consider the case of minimizing the shot noise in estimates of the acceleration field for an N-body model of a spherical Hernquist model. Models constructed using our scheme easily yield a factor ~100 reduction in the variance in the central acceleration field when compared to a traditional equal-mass model with the same number of particles. When evolving both models with a real N-body code, the diffusion coefficients in our model are reduced by a similar factor. Therefore, for certain types of problems, our scheme is a practical method for reducing the two-body relaxation effects, thereby bringing the N-body simulations closer to the collisionless ideal.
We present sixth- and eighth-order Hermite integrators for astrophysical $N$-body simulations, which use the derivatives of accelerations up to second order ({it snap}) and third order ({it crackle}). These schemes do not require previous values for
the corrector, and require only one previous value to construct the predictor. Thus, they are fairly easy to implemente. The additional cost of the calculation of the higher order derivatives is not very high. Even for the eighth-order scheme, the number of floating-point operations for force calculation is only about two times larger than that for traditional fourth-order Hermite scheme. The sixth order scheme is better than the traditional fourth order scheme for most cases. When the required accuracy is very high, the eighth-order one is the best. These high-order schemes have several practical advantages. For example, they allow a larger number of particles to be integrated in parallel than the fourth-order scheme does, resulting in higher execution efficiency in both general-purpose parallel computers and GRAPE systems.
Many problems encountered in plasma physics require a description by kinetic equations, which are posed in an up to six-dimensional phase space. A direct discretization of this phase space, often called the Eulerian approach, has many advantages but
is extremely expensive from a computational point of view. In the present paper we propose a dynamical low-rank approximation to the Vlasov--Poisson equation, with time integration by a particular splitting method. This approximation is derived by constraining the dynamics to a manifold of low-rank functions via a tangent space projection and by splitting this projection into the subprojections from which it is built. This reduces a time step for the six- (or four-) dimensional Vlasov--Poisson equation to solving two systems of three- (or two-) dimensional advection equations over the time step, once in the position variables and once in the velocity variables, where the size of each system of advection equations is equal to the chosen rank. By a hierarchical dynamical low-rank approximation, a time step for the Vlasov--Poisson equation can be further reduced to a set of six (or four) systems of one-dimensional advection equations, where the size of each system of advection equations is still equal to the rank. The resulting systems of advection equations can then be solved by standard techniques such as semi-Lagrangian or spectral methods. Numerical simulations in two and four dimensions for linear Landau damping, for a two-stream instability and for a plasma echo problem highlight the favorable behavior of this numerical method and show that the proposed algorithm is able to drastically reduce the required computational effort.
We describe the numerical code N-MODY, a parallel particle-mesh code for collisionless N-body simulations in modified Newtonian dynamics (MOND). N-MODY is based on a numerical potential solver in spherical coordinates that solves the non-linear MOND
field equation, and is ideally suited to simulate isolated stellar systems. N-MODY can be used also to compute the MOND potential of arbitrary static density distributions. A few applications of N-MODY indicate that some astrophysically relevant dynamical processes are profoundly different in MOND and in Newtonian gravity with dark matter.
Initial conditions for (Newtonian) cosmological N-body simulations are usually set by re-scaling the present-day power spectrum obtained from linear (relativistic) Boltzmann codes to the desired initial redshift of the simulation. This back-scaling m
ethod can account for the effect of inhomogeneous residual thermal radiation at early times, which is absent in the Newtonian simulations. We analyse this procedure from a fully relativistic perspective, employing the recently-proposed Newtonian motion gauge framework. We find that N-body simulations for LambdaCDM cosmology starting from back-scaled initial conditions can be self-consistently embedded in a relativistic space-time with first-order metric potentials calculated using a linear Boltzmann code. This space-time coincides with a simple N-body gauge for z<50 for all observable modes. Care must be taken, however, when simulating non-standard cosmologies. As an example, we analyse the back-scaling method in a cosmology with decaying dark matter, and show that metric perturbations become large at early times in the back-scaling approach, indicating a breakdown of the perturbative description. We suggest a suitable forwards approach for such cases.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
