Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSimulations of planet-disc interactions using Smoothed Particle Hydrodynamics
206
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We have performed Smoothed Particle Hydrodynamics (SPH) simulations to study the time evolution of one and two protoplanets embedded in a protoplanetary accretion disc. We investigate accretion and migration rates of a single protoplanet depending on several parameters of the protoplanetary disc, mainly viscosity and scale height. Additionally, we consider the influence of a second protoplanet in a long time simulation and examine the migration of the two planets in the disc, especially the growth of eccentricity and chaotic behaviour. One aim of this work is to establish the feasibility of SPH for such calculations considering that usually only grid-based methods are adopted. To resolve shocks and to prevent particle penetration, we introduce a new approach for an artificial viscosity, which consists of an additional artificial bulk viscosity term in the SPH-representation of the Navier-Stokes equation. This allows for an accurate treatment of the physical kinematic viscosity to describe the shear, without the use of an artificial shear viscosity.
rate research
Read More
In this study we investigated the capabilities of the mesh-free, Lagrangian particle method (Smoothed Particle Hydrodynamics, SPH) to simulate the detailed hydrodynamic processes generated by both spilling and plunging breaking waves within the surf zone. The weakly-compressible SPH code DualSPHysics was applied to simulate wave breaking over two distinct bathymetric profiles (a plane beach and fringing reef) and compared to experimental flume measurements of waves, flows, and mean water levels. Despite the simulations spanning very different wave breaking conditions (including an extreme case with violently plunging waves on an effectively dry reef slope), the model was able to reproduce a wide range of relevant surf zone hydrodynamic processes using a fixed set of numerical parameters. This included accurate predictions of the nonlinear evolution of wave shapes (e.g., asymmetry and skewness properties), rates of wave dissipation within the surf zone, and wave setup distributions. By using this mesh-free approach, the model was able to resolve the critical crest region within the breaking waves, which provided robust predictions of the wave-induced mass fluxes within the surf zone responsible for the undertow. Within this breaking crest region, the model results capture how the potential energy of the organized wave motion is initially converted to kinetic energy and then dissipated, which reproduces the distribution of wave forces responsible for wave setup generation across the surf zone. Overall, the results reveal how the mesh-free SPH approach can accurately reproduce the detailed wave breaking processes with comparable skill to state-of-the-art mesh-based Computational Fluid Dynamics (CFD) models, and thus can be applied to provide valuable new physical insight into surf zone dynamics.
The gravitational interaction between a protoplanetary disc and planetary sized bodies that form within it leads to the exchange of angular momentum, resulting in migration of the planets and possible gap formation in the disc for more massive planets. In this article, we review the basic theory of disc-planet interactions, and discuss the results of recent numerical simulations of planets embedded in protoplanetary discs. We consider the migration of low mass planets and recent developments in our understanding of so-called type I migration when a fuller treatment of the disc thermodynamics is included. We discuss the runaway migration of intermediate mass planets (so-called type III migration), and the migration of giant planets (type II migration) and the associated gap formation in the disc. The availability of high performance computing facilities has enabled global simulations of magnetised, turbulent discs to be computed, and we discuss recent results for both low and high mass planets embedded in such discs.
We develop a formalism for the identification and accurate estimation of the strength of structure formation shocks during cosmological smoothed particle hydrodynamics simulations. Shocks not only play a decisive role for the thermalization of gas in virialising structures but also for the acceleration of relativistic cosmic rays (CRs) through diffusive shock acceleration. Our formalism is applicable both to ordinary non-relativistic thermal gas, and to plasmas composed of CRs and thermal gas. To this end, we derive an analytical solution to the one-dimensional Riemann shock tube problem for a composite plasma of CRs and thermal gas. We apply our methods to study the properties of structure formation shocks in high-resolution hydrodynamic simulations of the LCDM model. We find that most of the energy is dissipated in weak internal shocks with Mach numbers M~2 which are predominantly central flow shocks or merger shock waves traversing halo centres. Collapsed cosmological structures are surrounded by external shocks with much higher Mach numbers up to M~1000, but they play only a minor role in the energy balance of thermalization. We show that after the epoch of cosmic reionisation the Mach number distribution is significantly modified by an efficient suppression of strong external shock waves due to the associated increase of the sound speed of the diffuse gas. Invoking a model for CR acceleration in shock waves, we find that the average strength of shock waves responsible for CR energy injection is higher than that for shocks that dominate the thermalization of the gas. When combined with radiative dissipation and star formation, our formalism can also be used to study CR injection by supernova shocks, or to construct models for shock-induced star formation in the interstellar medium. (abridged)
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 discontinuity 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.
At present, the giant impact (GI) is the most widely accepted model for the origin of the Moon. Most of the numerical simulations of GI have been carried out with the smoothed particle hydrodynamics (SPH) method. Recently, however, it has been pointed out that standard formulation of SPH (SSPH) has difficulties in the treatment of a contact discontinuity such as a core-mantle boundary and a free surface such as a planetary surface. This difficulty comes from the assumption of differentiability of density in SSPH. We have developed an alternative formulation of SPH, density independent SPH (DISPH), which is based on differentiability of pressure instead of density to solve the problem of a contact discontinuity. In this paper, we report the results of the GI simulations with DISPH and compare them with those obtained with SSPH. We found that the disk properties, such as mass and angular momentum produced by DISPH is different from that of SSPH. In general, the disks formed by DISPH are more compact: while formation of a smaller mass moon for low-oblique impacts is expected with DISPH, inhibition of ejection would promote formation of a larger mass moon for high-oblique impacts. Since only the improvement of core-mantle boundary significantly affects the properties of circumplanetary disks generated by GI and DISPH has not been significantly improved from SSPH for a free surface, we should be very careful when some conclusions are drawn from the numerical simulations for GI. And it is necessary to develop the numerical hydrodynamical scheme for GI that can properly treat the free surface as well as the contact discontinuity.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
