Source author record

G. Lanzafame

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

6works
4topics
2close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

6 published item(s)

preprint2012arXiv

Implicit integrations for SPH in Semi-Lagrangian approach: application to the accretion disc modelling in a microquasar

Current explicit integration techniques in fluid dynamics are deeply limited by the Courant-Friedrichs-Lewy condition of the time step progression, based on the adopted spatial resolution coupled with the maximum value between the kinetic velocity or the signal transmission speed in the computational domain. Eulerian implicit integration techniques, even though more time consuming, can allow to perform stable computational fluid dynamics paying the price of a relatively larger inaccuracy in the calculations, without suffering such a strict temporal limitation. In this paper, we present a simple and effective scheme to perform Free Lagrangian Smooth Particle Hydrodynamics (SPH) implicit integrations in Semi-Lagrangian approach without any Jacobian matrix inversion operations for viscous Navier-Stokes flows. Applications to SPH accretion disc simulation around a massive black hole (MBH) in a binary stellar system are shown, together with the comparison to the same results obtained according to the traditional explicit integration techniques. Some 1D and 2D critical tests are also discussed to check the validity of the technique.

preprint2011arXiv

An approach for solving the boundary free edge difficulties in SPH modelling: application to a viscous accretion disc in close binaries

In this work, we propose a SPH interpolating Kernel reformulation suitable also to treat free edge boundaries in the computational domain. Application to both inviscid and viscous stationary low compressibility accretion disc models in Close Binaries (CB) are shown. The investigation carried out in this paper is a consequence of the fact that a low compressibility modelling is crucial to check numerical reliability. Results show that physical viscosity supports a well-bound accretion disc formation, despite the low gas compressibility, when a Gaussian-derived Kernel (from the Error Function) is assumed, in extended particle range - whose Half Width at Half Maximum (HWHM) is fixed to a constant $h$ value - without any spatial restrictions on its radial interaction (hereinafter GASPHER). At the same time, GASPHER ensures adequate particle interpolations at the boundary free edges. Both SPH and adaptive SPH (hereinafter ASPH) methods lack accuracy if there are not constraints on the boundary conditions, in particular at the edge of the particle envelope: Free Edge (FE) conditions. In SPH, an inefficient particle interpolation involves a few neighbour particles; instead, in the second case, non-physical effects involve both the boundary layer particles themselves and the radial transport. Either in a regime where FE conditions involve the computational domain, or in a viscous fluid dynamics, or both, a GASPHER scheme can be rightly adopted in such troublesome physical regimes. Despite the applied low compressibiity condition, viscous GASPHER model shows clear spiral pattern profiles demonstrating the better quality of results compared to SPH viscous ones. Moreover a successful comparison of results concerning GASPHER 1D inviscid shock tube with analytical solution is also reported.

preprint2011arXiv

An approach to the Riemann problem in the light of a reformulation of the state equation for SPH inviscid ideal flows: a highlight on spiral hydrodynamics in accretion discs

In physically inviscid fluid dynamics, "shock capturing" methods adopt either an artificial viscosity contribution or an appropriate Riemann solver algorithm. These techniques are necessary to solve the strictly hyperbolic Euler equations if flow discontinuities (the Riemann problem) are to be solved. A necessary dissipation is normally used in such cases. An explicit artificial viscosity contribution is normally adopted to smooth out spurious heating and to treat transport phenomena. Such a treatment of inviscid flows is also widely adopted in the Smooth Particle Hydrodynamics (SPH) finite volume free Lagrangian scheme. In other cases, the intrinsic dissipation of Godunov-type methods is implicitly useful. Instead "shock tracking" methods normally use the Rankine-Hugoniot jump conditions to solve such problems. A simple, effective solution of the Riemann problem in inviscid ideal gases is here proposed, based on an empirical reformulation of the equation of state (EoS) in the Euler equations in fluid dynamics, whose limit for a motionless gas coincides with the classical EoS of ideal gases. The application of such an effective solution to the Riemann problem excludes any dependence, in the transport phenomena, on particle smoothing resolution length $h$ in non viscous SPH flows. Results on 1D shock tube tests, as well as examples of application for 2D turbulence and 2D shear flows are here shown. As an astrophysical application, a much better identification of spiral structures in accretion discs in a close binary (CB), as a result of this reformulation is also shown here.

preprint2010arXiv

An approach to the Riemann problem for SPH inviscid ideal flows

In the physically non viscous fluid dynamics, "shock capturing" methods adopt either an artificial viscosity contribution or an appropriate Riemann solver algorithm. These techniques are necessary to solve the strictly hyperbolic Euler equations if flow discontinuities (the Riemann problem) must be solved. A necessary dissipation is normally used in such cases. An explicit artificial viscosity contribution is normally adopted to smooth out spurious heating and to treat transport phenomena. Such a treatment of inviscid flows is also widely adopted in the Smooth Particle Hydrodynamics (SPH) finite volume free Lagrangian scheme. In other cases, the intrinsic dissipation of Godunov - type methods is implicitly useful. Instead "shock tracking" methods normally use the Rankine - Hugoniot jump conditions to solve such problem. A simple, effective solution of the Riemann problem in inviscid ideal gases is here proposed, based on an empirical reformulation of the equation of state (EoS) in the Euler equations in fluid dynamics, whose limit for a motionless gas coincides with the classical EoS of ideal gases. The application of such effective solution of the Riemann problem excludes any dependence, in the transport phenomena, on particle resolution length $h$ in non viscous SPH flows. Results on 1D shock tube tests are here shown.

preprint2010arXiv

An approach to the Riemann problem for SPH inviscid ideal flows: consequences for the state equation

In the non viscous fluid dynamics, Smooth Particle Hydrodynamics (SPH), as a free Lagrangian "shock capturing" method adopts either an artificial viscosity contribution or an appropriate Riemann solver technique. An explicit or an implicit dissipation, introduced in such techniques, is necessary to solve the Euler equations to solve flow discontinuities (the Riemann problem). Dissipation is useful to smooth out spurious heating and to treat transport phenomena. A simple, effective solution of the Riemann problem is here proposed, based on an empirical reformulation of the equation of state (EoS) in the Euler equations in fluid dynamics, whose limit for a motionless gas coincides with the classic EoS of ideal gases. Results on 1D shock tube tests are here shown, as well as a 3D transport application on accretion discs in close binaries (CBs).

preprint2010arXiv

Improvements in the X-ray luminosity function and constraints on the Cosmological parameters from X-ray luminous clusters

We show how to improve constraints on Ω_m, σ_8, and the dark-energy equation-of-state parameter, w, obtained by Mantz et al. (2008) from measurements of the X-ray luminosity function of galaxy clusters, namely MACS, the local BCS and the REFLEX galaxy cluster samples with luminosities L> 3 \times 10^{44} erg/s in the 0.1--2.4 keV band. To this aim, we use Tinker et al. (2008) mass function instead of Jenkins et al. (2001) and the M-L relationship obtained from Del Popolo (2002) and Del Popolo et al. (2005). Using the same methods and priors of Mantz et al. (2008), we find, for a Λ$CDM universe, Ω_m=0.28^{+0.05}_{-0.04} and σ_8=0.78^{+0.04}_{-0.05}$ while the result of Mantz et al. (2008) gives less tight constraints $Ω_m=0.28^{+0.11}_{-0.07}$ and σ_8=0.78^{+0.11}_{-0.13}. In the case of a wCDM model, we find Ω_m=0.27^{+0.07}_{-0.06}, $σ_8=0.81^{+0.05}_{-0.06}$ and $w=-1.3^{+0.3}_{-0.4}$, while in Mantz et al. (2008) they are again less tight Ω_m=0.24^{+0.15}_{-0.07}, σ_8=0.85^{+0.13}_{-0.20} and w=-1.4^{+0.4}_{-0.7}. Combining the XLF analysis with the f_{gas}+CMB+SNIa data set results in the constraint Ω_m=0.269 \pm 0.012, σ_8=0.81 \pm 0.021 and w=-1.02 \pm 0.04, to be compared with Mantz et al. (2008), Ω_m=0.269 \pm 0.016, σ_8=0.82 \pm 0.03 and w=-1.02 \pm 0.06. The tightness of the last constraints obtained by Mantz et al. (2008), are fundamentally due to the tightness of the $f_{gas}$+CMB+SNIa constraints and not to their XLF analysis. Our findings, consistent with w=-1, lend additional support to the cosmological-constant model.