Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA truly Newtonian softening length for disc simulations
437
0
0.0
(
0
)
Authors
Jean-Marc Hure
Ask ChatGPT about the research
No Arabic abstract
The softened point mass model is commonly used in simulations of gaseous discs including self-gravity while the value of associated length lambda remains, to some degree, controversial. This ``parameter is however fully constrained when, in a discretized disc, all fluid cells are demanded to obey Newtons law. We examine the topology of solutions in this context, focusing on cylindrical cells more or less vertically elongated. We find that not only the nominal length depends critically on the cells shape (curvature, radial extension, height), but it is either a real or an imaginary number. Setting lambda as a fraction of the local disc thickness -- as usually done -- is indeed not the optimal choice. We then propose a novel prescription valid irrespective of the disc properties and grid spacings. The benefit, which amounts to 2-3 more digits typically, is illustrated in a few concrete cases. A detailed mathematical analysis is in progress.
rate research
Read More
Gravitational softening length is one of the key parameters to properly set up a cosmological $N$-body simulation. In this paper, we perform a large suit of high-resolution $N$-body simulations to revise the optimal softening scheme proposed by Power et al. (P03). Our finding is that P03 optimal scheme works well but is over conservative. Using smaller softening lengths than that of P03 can achieve higher spatial resolution and numerically convergent results on both circular velocity and density profiles. However using an over small softening length overpredicts matter density at the inner most region of dark matter haloes. We empirically explore a better optimal softening scheme based on P03 form and find that a small modification works well. This work will be useful for setting up cosmological simulations.
This document describes the general process of setting up, running, and analysing disc galaxy simulations using the freely available program Phantom of RAMSES (PoR). This implements Milgromian Dynamics (MOND) with a patch to the RAMSES grid-based $N$-body and hydrodynamical code that uses adaptive mesh refinement. We discuss the procedure of setting up isolated and interacting disc galaxy initial conditions for PoR, running the simulations, and analysing the results. This manual also concisely documents all previously developed MOND simulation codes and the results obtained with them.
We have modeled hydrodynamical tori in the Klu{z}niak-Lee pseudo-Newtonian potential. The tori in equilibrium were perturbed with uniform sub-sonic velocity fields: vertical, radial and diagonal respectively, and allowed to evolve in time. We identify the eigenmodes corresponding to those of slender tori. The results of our simulations are relevant to the investigation of high-frequency quasi-periodic oscillations observed in stellar-mass black hole binaries.
I describe attempts to identify and understand the most isolated galaxies starting from my 1983 Leiden PhD thesis, continuing through a string of graduate theses on various aspects of this topic, and concluding with an up-to-date account of the difficulty to find really isolated objects. The implication of some of the findings revealed on the way and presented here is that the nearby Universe may contain many small dark-matter haloes, and that some such haloes may possibly be accreting intergalactic gas to form dwarf galaxies.
Linear perturbation is used to investigate the effect of gravitational softening on the retrieved two-armed spiral eigenmodes of razor-thin stellar discs. We explore four softening kernels with different degrees of gravity bias, and with/without compact support (compact in the sense that they yield exactly Newtonian forces outside the softening kernel). These kernels are applied to two disc galaxy models with well-known unsoftened unstable modes. We illustrate quantitatively the importance of a vanishing linear gravity bias to yield accurate frequency estimates of the unstable modes. As such, Plummer softening, while very popular amongst simulators, performs poorly in our tests. The best results, with excellent agreement between the softened and unsoftened mode properties, are obtained with softening kernels that have a reduced gravity bias, obtained by compensating for the sub-Newtonian forces at small interparticle distances with slightly super-Newtonian forces at radii near the softening length. We present examples of such kernels that, moreover, are analytically simple and computationally cheap. Finally, these results light the way to the construction of softening methods with even smaller gravity bias, although at the price of increasingly complex kernels.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
