Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userConditions for justifying single-fluid approximation for charged and neutral dust fluids and a smoothed particle magnetohydrodynamics method for dust-gas mixture
72
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We describe a numerical scheme for magnetohydrodynamics simulations of dust-gas mixture by extending smoothed particle magnetohydrodynamics. We employ the single-species particle approach to describe dust-gas mixture with several modifications from the previous studies. We assume that the charged and neutral dusts can be treated as single-fluid and the electro-magnetic force acts on the gas and that on the charged dust is negligible. The validity of these assumption in the context of protostar formation is not obvious and is extensively evaluated. By investigating the electromagnetic force and electric current with terminal velocity approximation, it is found that as the dust size increases, the contribution of dust to them becomes smaller and negligible. We conclude that our assumptions of the electro-magnetic force on the dusts is negligible are valid for the dust size with a d & 10{mu}m. On the other hand, they do not produce the numerical artifact for the dust a d . 10{mu}m in envelope and disk where the perfect coupling between gas and dusts realizes. However, we also found that our assumptions may break down in outflow (or under environment with very strong magnetic field and low density) for the dust a d . 10{mu}m. We conclude that our assumptions are valid in almost all cases where macroscopic dust dynamics is important in the context of protostar formation. We conduct numerical tests of dusty wave, dusty magnetohydrodynamics shock, and gravitational collapse of magnetized cloud core with our simulation code. The results show that our numerical scheme well reproduces the dust dynamics in the magnetized medium.
rate research
Read More
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.
We describe a simple method for simulating the dynamics of small grains in a dusty gas, relevant to micron-sized grains in the interstellar medium and grains of centimetre size and smaller in protoplanetary discs. The method involves solving one extra diffusion equation for the dust fraction in addition to the usual equations of hydrodynamics. This diffusion approximation for dust is valid when the dust stopping time is smaller than the computational timestep. We present a numerical implementation using Smoothed Particle Hydrodynamics (SPH) that is conservative, accurate and fast. It does not require any implicit timestepping and can be straightforwardly ported into existing 3D codes.
We present a fix to the overdamping problem found by Laibe & Price (2012) when simulating strongly coupled dust-gas mixtures using two different sets of particles using smoothed particle hydrodynamics. Our solution is to compute the drag at the barycentre between gas and dust particle pairs when computing the drag force by reconstructing the velocity field, similar to the procedure in Godunov-type solvers. This fixes the overdamping problem at negligible computational cost, but with additional memory required to store velocity derivatives. We employ slope limiters to avoid spurious oscillations at shocks, finding the van Leer Monotonized Central limiter most effective.
The development of smooth particle magnetohydrodynamic (SPMHD) has significantly improved the simulation of complex astrophysical processes. However, the preservation the solenoidality of the magnetic field is still a severe problem for the MHD. A formulation of the induction equation with a vector potential would solve the problem. Unfortunately all previous attempts suffered from instabilities. In the present work, we evolve the vector potential in the Coulomb gauge and smooth the derived magnetic field for usage in the momentum equation. With this implementation we could reproduce classical test cases in a stable way. A simple test case demonstrates the possible failure of widely used direct integration of the magnetic field, even with the usage of a divergence cleaning method.
Supersonic turbulence is believed to be at the heart of star formation. We have performed smoothed particle magnetohydrodynamics (SPMHD) simulations of the small-scale dynamo amplification of magnetic fields in supersonic turbulence. The calculations use isothermal gas driven at rms velocity of Mach 10 so that conditions are representative of star-forming molecular clouds in the Milky Way. The growth of magnetic energy is followed for 10 orders in magnitude until it reaches saturation, a few percent of the kinetic energy. The results of our dynamo calculations are compared with results from grid-based methods, finding excellent agreement on their statistics and their qualitative behaviour. The simulations utilise the latest algorithmic developments we have developed, in particular, a new divergence cleaning approach to maintain the solenoidal constraint on the magnetic field and a method to reduce the numerical dissipation of the magnetic shock capturing scheme. We demonstrate that our divergence cleaning method may be used to achieve $ abla cdot {bf B}=0$ to machine precision, albeit at significant computational expense.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
