Source author record

R. Valdarnini

R. Valdarnini appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph astro-ph.IM physics.comp-ph physics.flu-dyn astro-ph.CO

Catalog footprint

What is connected

5works

5topics

3close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2016arXiv

Improved Performances in Subsonic Flows of an SPH Scheme with Gradients Estimated using an Integral Approach

In this paper we present results from a series of hydrodynamical tests aimed at validating the performance of a smoothed particle hydrodynamics (SPH) formulation in which gradients are derived from an integral approach. We specifically investigate the code behavior with subsonic flows, where it is well known that zeroth-order inconsistencies present in standard SPH make it particularly problematic to correctly model the fluid dynamics. In particular we consider the Gresho-Chan vortex problem, the growth of Kelvin-Helmholtz instabilities, the statistics of driven subsonic turbulence and the cold Keplerian disc problem. We compare simulation results for the different tests with those obtained, for the same initial conditions, using standard SPH. We also compare the results with the corresponding ones obtained previously with other numerical methods, such as codes based on a moving-mesh scheme or Godunov-type Lagrangian meshless methods. We quantify code performances by introducing error norms and spectral properties of the particle distribution, in a way similar to what was done in other works. We find that the new SPH formulation exhibits strongly reduced gradient errors and outperforms standard SPH in all of the tests considered. In fact, in terms of accuracy we find good agreement between the simulation results of the new scheme and those produced using other recently proposed numerical schemes. These findings suggest that the proposed method can be successfully applied for many astrophysical problems in which the presence of subsonic flows previously limited the use of SPH, with the new scheme now being competitive in these regimes with other numerical methods.

preprint2012arXiv

Hydrodynamic capabilities of an SPH code incorporating an artificial conductivity term with a gravity-based signal velocity

This paper investigates the hydrodynamic performances of an SPH code incorporating an artificial heat conductivity term in which the adopted signal velocity is applicable when gravity is present. In accordance with previous findings it is shown that the performances of SPH to describe the development of Kelvin-Helmholtz instabilities depend strongly on the consistency of the initial condition set-up and on the leading error in the momentum equation due to incomplete kernel sampling. An error and stability analysis shows that the quartic B-spline kernel (M_5) possesses very good stability properties and we propose its use with a large neighbor number, between ~50 (2D) to ~ 100 (3D), to improve convergence in simulation results without being affected by the so-called clumping instability. SPH simulations of the blob test show that in the regime of strong supersonic flows an appropriate limiting condition, which depends on the Prandtl number, must be imposed on the artificial conductivity SPH coefficients in order to avoid an unphysical amount of heat diffusion. Results from hydrodynamic simulations that include self-gravity show profiles of hydrodynamic variables that are in much better agreement with those produced using mesh-based codes. In particular, the final levels of core entropies in cosmological simulations of galaxy clusters are consistent with those found using AMR codes. Finally, results of the Rayleigh-Taylor instability test demonstrate that in the regime of very subsonic flows the code has still several difficulties in the treatment of hydrodynamic instabilities. These problems being intrinsically due to the way in which in standard SPH gradients are calculated and not to the implementation of the artificial conductivity term.

preprint2010arXiv

The impact of numerical viscosity in SPH simulations of galaxy clusters

A SPH code employing a time-dependent artificial viscosity scheme is used to construct a large set of N-body/SPH cluster simulations for studying the impact of artificial viscosity on the thermodynamics of the ICM and its velocity field statistical properties. Spectral properties of the gas velocity field are investigated by measuring for the simulated clusters the velocity power spectrum E(k). The longitudinal component E_c(k) exhibits over a limited range a Kolgomorov-like scaling k^{-5/3}, whilst the solenoidal power spectrum component E_s(k) is strongly influenced by numerical resolution effects. The dependence of the spectra E(k) on dissipative effects is found to be significant at length scales 100-300Kpc, with viscous damping of the velocities being less pronounced in those runs with the lowest artificial viscosity. The turbulent energy density radial profile E_{turb}(r) is strongly affected by the numerical viscosity scheme adopted in the simulations, with the turbulent-to-total energy density ratios being higher in the runs with the lowest artificial viscosity settings and lying in the range between a few percent and ~10%. These values are in accord with the corresponding ratios extracted from previous cluster simulations realized using mesh-based codes. At large cluster radii, the mass correction terms to the hydrostatic equilibrium equation are little affected by the numerical viscosity of the simulations, showing that the X-ray mass bias is already estimated well in standard SPH simulations. Finally, simulations in which the gas can cool radiatively are characterized by the presence in the cluster inner regions of high levels of turbulence, generated by the interaction of the compact cool gas core with the ambient medium.

preprint1995arXiv

Three-point correlation function of galaxy clusters in cosmological models: a strong dependence on triangle shape

In this paper, we use large $\pppm$ N-body simulations to study the three-point correlation function $ζ$ of clusters in two theoretical models. The first model (LCDM) is a low-density flat model of $Ω_0=0.3$, $Λ_0=0.7$ and $h=0.75$, and the second model (PIM) is an Einstein-de-Sitter model with its linear power spectrum obtained from observations. We found that the scaled function $Q(r,u,v)$, which is defined as the ratio of $ζ(r, ru, ru+rv)$ to the hierarchical sum $ξ(r)ξ(ru)+ ξ(ru) ξ(ru+rv) +ξ(ru+rv)ξ(r)$ (where $ξ$ is the two-point correlation function of clusters), depends weakly on $r$ and $u$, but very strongly on $v$. $Q(r,u,v)$ is about 0.2 at $v=0.1$ and 1.8 at $v=0.9$. A model of $Q(r,u,v)=Θ10^{1.3v^2}$ can fit the data of $ζ$ very nicely with $Θ\approx 0.14$. This model is found to be universal for the LCDM clusters and for the PIM clusters. Furthermore, $Q(r,u,v)$ is found to be insensitive to the cluster richness. We have compared our N-body results with simple analytical theories of cluster formation, like the peak theories or the local maxima theories. We found that these theories do not provide an adequate description for the three-point function of clusters. We have also examined the observational data of $ζ$ presently available, and have not found any contradiction between the observations and our model predictions. The $v$-dependence of $q$ in a projected catalogue of clusters is shown to be much weaker than the $v$-dependence of $Q$ found in the three-dimensional case. It would be important to search for the $v$-dependence of $Q$ in redshift samples of rich clusters.

preprint1994arXiv

The 1-Point Cluster Distribution Function and its Moments

We derive the 1-point probability density function of the smoothed 3-D Abell-ACO cluster density field and we compare it with that of artificial cluster samples, generated as high peaks of a Gaussian field in such a way that they reproduce the low-order (2- and 3-point) correlation functions and the observed cluster selection functions. We find that both real and simulated {\em pdf}'s are well approximated by a log-normal distribution even when the Gaussian smoothing radius is as large as 40 $h^{-1}$ Mpc. Furthermore the low-order moments of the {\em pdf} are found to obey a relation $γ= S_{3} σ^4$, with $γ$ being the skewness and $S_{3}\approx 1.8$. These results are consistent with clusters being high-peaks of an underlying initial Gaussian density field. A by-product of our analysis is that when we rescale the {\em pdf} cluster moments to those of the QDOT-IRAS galaxies, using linear biasing with $b_{cI}\sim 4.5$ and for the common smoothing radius of 20 $h^{-1}$ Mpc, we find them to be significantly smaller than those directly estimated from the QDOT data by Saunders et al. (1991).

Institution

Affiliation not imported yet

This author record came from a source that does not expose affiliation metadata. Once the author claims the profile or we enrich the record from another provider, this section will link to the concrete institution.

Topic footprint