ترغب بنشر مسار تعليمي؟ اضغط هنا

Numerical aspects of Giant Impact simulations

80   0   0.0 ( 0 )
 نشر من قبل Christian Reinhardt
 تاريخ النشر 2017
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

In this paper we present solutions to three short comings of Smoothed Particles Hydrodynamics (SPH) encountered in previous work when applying it to Giant Impacts. First we introduce a novel method to obtain accurate SPH representations of a planets equilibrium initial conditions based on equal area tessellations of the sphere. This allows one to imprint an arbitrary density and internal energy profile with very low noise which substantially reduces computation because these models require no relaxation prior to use. As a consequence one can significantly increase the resolution and more flexibly change the initial bodies to explore larger parts of the impact parameter space in simulations. The second issue addressed is the proper treatment of the matter/vacuum boundary at a planets surface with a modified SPH density estimator that properly calculates the density stabilizing the models and avoiding an artificially low density atmosphere prior to impact. Further we present a novel SPH scheme that simultaneously conserves both energy and entropy for an arbitrary equation of state. This prevents loss of entropy during the simulation and further assures that the material does not evolve into unphysical states. Application of these modifications to impact simulations for different resolutions up to $6.4 cdot 10^6$ particles show a general agreement with prior result. However, we observe resolution dependent differences in the evolution and composition of post collision ejecta. This strongly suggests that the use of more sophisticated equations of state also demands a large number of particles in such simulations.



قيم البحث

اقرأ أيضاً

173 - A. Gabbana , D. Simeoni , S. Succi 2019
We present an analytical derivation of the transport coefficients of a relativistic gas in (2+1) dimensions for both Chapman-Enskog (CE) asymptotics and Grads expansion methods. Moreover, we develop a systematic calibration method, connecting the rel axation time of relativistic kinetic theory to the transport parameters of the associated dissipative hydrodynamic equations. Comparison between the analytical results and numerical simulations, shows that the CE method correctly captures dissipative effects, while Grads method does not. The resulting calibration procedure based on the CE method opens the way to the quantitative kinetic description of dissipative relativistic fluid dynamics under fairly general conditions, namely flows with strongly non-linearities, in non-ideal geometries, across both ultra-relativistic and near-non-relativistic regimes.
The mass balance of mountain glaciers is of interest for several applications (local hydrology or climate projections), and turbulent fluxes can be an important contributor to glacier surface mass balance during strong melting events. The underlying complex terrain leads to spatial heterogeneity and non-stationarity of turbulent fluxes. Due to the contribution of thermally-induced flows and gravity waves, exchange mechanisms are fully three-dimensional, instead of only vertical. Additionally, glaciers have their own distinct microclimate, governed by a down-glacier katabatic wind, which protects the glacier ice and interacts with the surrounding flows on multiple scales. In this study, we perform large-eddy simulations with the WRF model with dx=48 m to gain insight on the boundary-layer processes over an Alpine valley glacier, the Hintereisferner (HEF). We choose two case studies from a measurement campaign (August 2018) with different synoptic wind directions (South-West and North-West). Model evaluation with an array of eddy-covariance stations on the glacier tongue and surroundings reveals that WRF is able to simulate the general glacier boundary-layer structure. Under southwesterly airflow, the down-glacier wind is supported by the South-Western synoptic wind direction, a stable boundary layer is present over the ice surface, and local processes govern the turbulence kinetic energy production. Under northwesterly airflow, a cross-glacier valley flow and a breaking gravity wave lead strong turbulent mixing and to the subsequent erosion of the glacier boundary layer. Stationarity analyses of the sensible heat flux suggest non-stationary behaviour for both case study days, while non-stationarity is highest on the NW day during the gravity-wave event. These results suggest that the synoptic wind direction has, in addition to upstream topography and the atmospheric stability, a strong impact on whether a local glacier boundary layer can form or not, influencing whether a glacier is able to maintain its own microclimate.
The rheology of a three-dimensional granular jet during an impact is investigated numerically. The cone-like scattering pattern and the sheet-like pattern observed in an experiment [X. Cheng, et al. Phys. Rev. Lett. 99, 188001 (2007)] can be reproduc ed through our calculation. We discuss the constitutive equation for granular jet impact in terms of our simulation. From the analysis of an effective friction constant, which is the ratio between the shear stress and the pressure the assumption of the zero yield stress would be natural in our setup and the shear visocity is not small in contrast to the suggestion by the experiment.
The purpose of this article is numerical verification of the thory of weak turbulence. We performed numerical simulation of an ensemble of nonlinearly interacting free gravity waves (swell) by two different methods: solution of primordial dynamical e quations describing potential flow of the ideal fluid with a free surface and, solution of the kinetic Hasselmann equation, describing the wave ensemble in the framework of the theory of weak turbulence. Comparison of the results demonstrates pretty good applicability of the weak turbulent approach.
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 pointe d 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.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا