اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدParticle splitting in smoothed particle hydrodynamics based on Voronoi diagram
109
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
We present a novel method for particle splitting in smoothed particle hydrodynamics simulations. Our method utilizes the Voronoi diagram for a given particle set to determine the position of fine daughter particles. We perform several test simulations to compare our method with a conventional splitting method in which the daughter particles are placed isotropically over the local smoothing length. We show that, with our method, the density deviation after splitting is reduced by a factor of about two compared with the conventional method. Splitting would smooth out the anisotropic density structure if the daughters are distributed isotropically, but our scheme allows the daughter particles to trace the original density distribution with length scales of the mean separation of their parent. We apply the particle splitting to simulations of the primordial gas cloud collapse. The thermal evolution is accurately followed to the hydrogen number density of 10^12 /cc. With the effective mass resolution of ~10^-4 Msun after the multi-step particle splitting, the protostellar disk structure is well resolved. We conclude that the method offers an efficient way to simulate the evolution of an interstellar gas and the formation of stars.
قيم البحث
اقرأ أيضاً
In this paper, we present a new formulation of smoothed particle hydrodynamics (SPH), which, unlike the standard SPH (SSPH), is well-behaved at the contact discontinuity. The SSPH scheme cannot handle discontinuities in density (e.g. the contact disc
ontinuity and the free surface), because it requires that the density of fluid is positive and continuous everywhere. Thus there is inconsistency in the formulation of the SSPH scheme at discontinuities of the fluid density. To solve this problem, we introduce a new quantity associated with particles and density of that quantity. This density evolves through the usual continuity equation with an additional artificial diffusion term, in order to guarantee the continuity of density. We use this density or pseudo density, instead of the mass density, to formulate our SPH scheme. We call our new method as SPH with smoothed pseudo-density (SPSPH). We show that our new scheme is physically consistent and can handle discontinuities quite well.
There has been interest in recent years to assess the ability of astrophysical hydrodynamics codes to correctly model the Kelvin-Helmholtz instability. Smoothed particle hydrodynamics (SPH), in particular, has received significant attention, though t
here has yet to be a clear demonstration that SPH yields converged solutions that are in agreement with other methods. We have performed SPH simulations of the Kelvin-Helmholtz instability using the test problem put forward by Lecoanet et al (2016). We demonstrate that the SPH solutions converge to the reference solution in both the linear and non-linear regimes. Quantitative convergence in the strongly non-linear regime is achieved by using a physical Navier-Stokes viscosity and thermal conductivity. We conclude that standard SPH with an artificial viscosity can correctly capture the Kelvin-Helmholtz instability.
We perform simulations of the Kelvin-Helmholtz instability using smoothed particle hydrodynamics (SPH). The instability is studied both in the linear and strongly non-linear regimes. The smooth, well-posed initial conditions of Lecoanet et al. (2016)
are used, along with an explicit Navier-Stokes viscosity and thermal conductivity to enforce the evolution in the non-linear regime. We demonstrate convergence to the reference solution using SPH. The evolution of the vortex structures and the degree of mixing, as measured by a passive scalar `colour field, match the reference solution. Tests with an initial density contrast produce the correct qualitative behaviour. The L2 error of the SPH calculations decreases as the resolution is increased. The primary source of error is numerical dissipation arising from artificial viscosity, and tests with reduced artificial viscosity have reduced L2 error. A high-order smoothing kernel is needed in order to resolve the initial velocity amplitude of the seeded mode and inhibit excitation of spurious modes. We find that standard SPH with an artificial viscosity has no difficulty in correctly modelling the Kelvin-Helmholtz instability and yields convergent solutions.
We present a thorough numerical study on the MRI using the smoothed particle magnetohydrodynamics method (SPMHD) with the geometric density average force expression (GDSPH). We perform shearing box simulations with different initial setups and a wide
range of resolution and dissipation parameters. We show, for the first time, that MRI with sustained turbulence can be simulated successfully with SPH, with results consistent with prior work with grid-based codes. In particular, for the stratified boxes, our simulations reproduce the characteristic butterfly diagram of the MRI dynamo with saturated turbulence for at least 100 orbits. On the contrary, traditional SPH simulations suffer from runaway growth and develop unphysically large azimuthal fields, similar to the results from a recent study with mesh-less methods. We investigated the dependency of MRI turbulence on the numerical Prandtl number in SPH, focusing on the unstratified, zero net-flux case. We found that turbulence can only be sustained with a Prandtl number larger than $sim$2.5, similar to the critical values of physical Prandtl number found in grid-code simulations. However, unlike grid-based codes, the numerical Prandtl number in SPH increases with resolution, and for a fixed Prandtl number, the resulting magnetic energy and stresses are independent of resolution. Mean-field analyses were performed on all simulations, and the resulting transport coefficients indicate no $alpha$-effect in the unstratified cases, but an active $alphaOmega$ dynamo and a diamagnetic pumping effect in the stratified medium, which are generally in agreement with previous studies. There is no clear indication of a shear-current dynamo in our simulation, which is likely to be responsible for a weaker mean-field growth in the tall, unstratified, zero net-flux simulation.
We present Phantom, a fast, parallel, modular and low-memory smoothed particle hydrodynamics and magnetohydrodynamics code developed over the last decade for astrophysical applications in three dimensions. The code has been developed with a focus on
stellar, galactic, planetary and high energy astrophysics and has already been used widely for studies of accretion discs and turbulence, from the birth of planets to how black holes accrete. Here we describe and test the core algorithms as well as modules for magnetohydrodynamics, self-gravity, sink particles, H_2 chemistry, dust-gas mixtures, physical viscosity, external forces including numerous galactic potentials as well as implementations of Lense-Thirring precession, Poynting-Robertson drag and stochastic turbulent driving. Phantom is hereby made publicly available.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
