اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدA principal possibility for computer investigation of evolution of dynamical systems with independent on time accuracy
34
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
Extensive N-body simulations are among the key means for the study of numerous astrophysical and cosmological phenomena, so various schemes are developed for possibly higher accuracy computations. We demonstrate the principal possibility for revealing the evolution of a perturbed Hamiltonian system with an accuracy independent on time. The method is based on the Laplace transform and the derivation and analytical solution of an evolution equation in the phase space for the resolvent and using computer algebra.
قيم البحث
اقرأ أيضاً
The recently developed method (Paper 1) enabling one to investigate the evolution of dynamical systems with an accuracy not dependent on time is developed further. The classes of dynamical systems which can be studied by that method are much extended
, now including systems that are; (1) non-Hamiltonian, conservative; (2) Hamiltonian with time-dependent perturbation; (3) non-conservative (with dissipation). These systems cover various types of N-body gravitating systems of astrophysical and cosmological interest, such as the orbital evolution of planets, minor planets, artificial satellites due to tidal, non-tidal perturbations and thermal thrust, evolving close binary stellar systems, and the dynamics of accretion disks.
We present a new parallel code for computing the dynamical evolution of collisional N-body systems with up to N~10^7 particles. Our code is based on the the Henon Monte Carlo method for solving the Fokker-Planck equation, and makes assumptions of sph
erical symmetry and dynamical equilibrium. The principal algorithmic developments involve optimizing data structures, and the introduction of a parallel random number generation scheme, as well as a parallel sorting algorithm, required to find nearest neighbors for interactions and to compute the gravitational potential. The new algorithms we introduce along with our choice of decomposition scheme minimize communication costs and ensure optimal distribution of data and workload among the processing units. The implementation uses the Message Passing Interface (MPI) library for communication, which makes it portable to many different supercomputing architectures. We validate the code by calculating the evolution of clusters with initial Plummer distribution functions up to core collapse with the number of stars, N, spanning three orders of magnitude, from 10^5 to 10^7. We find that our results are in good agreement with self-similar core-collapse solutions, and the core collapse times generally agree with expectations from the literature. Also, we observe good total energy conservation, within less than 0.04% throughout all simulations. We analyze the performance of the code, and demonstrate near-linear scaling of the runtime with the number of processors up to 64 processors for N=10^5, 128 for N=10^6 and 256 for N=10^7. The runtime reaches a saturation with the addition of more processors beyond these limits which is a characteristic of the parallel sorting algorithm. The resulting maximum speedups we achieve are approximately 60x, 100x, and 220x, respectively.
We study the dynamical evolution of the young star cluster Arches and its dependence on the assumed initial stellar mass function (IMF). We perform many direct $N$-body simulations with various initial conditions and two different choices of IMFs. On
e is a standard Kroupa IMF without any mass segregation. The other is a radially dependent IMF, as presently observed in the Arches. We find that it is unlikely for the Arches to have attained the observed degree of mass segregation at its current age starting from a standard non-segregated Kroupa IMF. We also study the possibility of a collisional runaway developing in the first $sim 2-3 rm{Myr}$ of dynamical evolution. We find that the evolution of this cluster is dramatically different depending on the choice of IMF: if a primordially mass segregated IMF is chosen, a collisional runaway should always occur between $2-3 rm{Myr}$ for a broad range of initial concentrations. In contrast, for a standard Kroupa IMF no collisional runaway is predicted. We argue that if Arches was created with a mass segregated IMF similar to what is observed today then at the current cluster age a very unusual, high-mass star should be created. However, whether a collisional runaway leads to the formation of an intermediate-mass black hole (IMBH) depends strongly on the mass loss rate via winds from massive stars. Growth of stellar mass through collisions can be quenched by strong wind mass loss. In that case, the inter-cluster as well as intra-cluster medium are expected to have a significant Helium enrichment which may be observed via Helium recombination lines. The excess amount of gas lost in winds may also be observed via X-ray observations as diffused X-ray sources.
Evolution of a cluster of primordial black holes in the two-body relaxation approximation based on the Fokker-Planck equation is discussed. In our calculation, we consider the self-gravitating cluster with a wide range of black holes masses from $10^
{-4} M_{odot}$ up to $100 M_{odot}$ and the total mass $10^5 M_{odot}$. Moreover, we included a massive black hole in the cluster center which determines the evolution rate of the density profile in its vicinity.
In the last decades we witnessed an increase in studies of open clusters of the Galaxy, especially because of the good determination for a wide range of values of parameters such as age, distance, reddening, and proper motion. The reliable determinat
ion of the parameters strongly depends on the photometry available and especially on the U filter, which is used to obtain the color excess E(B-V) through the color-color diagram (U-B) by (B-V) by fitting a zero age main-sequence. Owing to the difficulty of performing photometry in the U band, many authors have tried to obtain E(B-V) without the filter. But because of the near linearity of the color-color diagrams that use the other bands, combined with the fact that most fitting procedures are highly subjective (many done by eye) the reliability of those results has always been questioned. Our group has recently developed, a tool that performs isochrone fitting in open-cluster photometric data with a global optimization algorithm, which removes the need to visually perform the fits and thus removes most of the related subjectivity. Here we apply our method to a set of synthetic clusters and two observed open clusters (Trumpler 1 and Melotte 105) using only photometry for the BVRI bands. Our results show that, considering the cluster structural variance caused only by photometric and Poisson sampling errors, our method is able to recover the synthetic cluster parameters with errors of less than 10% for a wide range of ages, distances, and reddening, which clearly demonstrates its potential. The results obtained for Trumpler 1 and Melotte 105 also agree well with previous literature values.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
