BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Rubén M. Cabezón

Rubén M. Cabezón contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.IMDistributed, Parallel, and Cluster Computingastro-ph.HEastro-ph.SRData Structures and Algorithms

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

6 published item(s)

preprint2022arXiv

Axisymmetric magneto-hydrodynamics with SPH

Many interesting terrestrial and astrophysical scenarios involving magnetic fields can be approached in axial geometry. Even though the Lagrangian smoothed particle hydrodynamics (SPH) technique has been successfully extended to handle magneto-hydrodynamic (MHD) problems, a well-verified, axisymmetric MHD scheme based on the SPH technique does not exist. In this work, we propose and check a new axisymmetric MHD hydrodynamic code that can be applied to astrophysical and engineering problems which display an adequate geometry. We show that a hydrodynamic code built on these axisymmetric premises is able to produce similar results to standard 3D-SPHMHD codes but with much lesser computational effort.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Conservative, density-based smoothed particle hydrodynamics with improved partition of the unity and better estimation of gradients

The correct evaluation of gradients is at the cornerstone of the smoothed particle hydrodynamics (SPH) technique. Using an integral approach to estimate gradients has proven to enhance accuracy substantially. Such approach retains the Lagrangian structure of SPH equations and is fully conservative. But, in practice, it is difficult to make the Lagrangian formulation totally consistent to an exact partition of the unity. In this paper we study, among other things, the connection between the choice of the volume elements (VEs), which enters in the SPH summations, and the accuracy in the gradient estimation within the integral approach scheme (ISPH). A new variant of VEs are proposed which improve the partition of the unity and is fully compatible with the Lagrangian formulation of SPH, including the grad-h corrections. Using analytic considerations, simple static toy models in 1D, and a few full 3D test cases, we show that any improvement in the partition of the unity also leads to a better calculation of gradients when the integral approach is used jointly. Additionally, we propose an easy-to-implement modification of the ISPH scheme, which makes it more flexible and better suited to handle sharp density contrasts. The ISPH code built with the proposed scheme has been validated with a good number of standard tests, some of them involving contact discontinuities. The performance of the code was excellent in all of them, showing that an improvement in the partition of the unity is not detrimental of the good conservation of energy, momentum, and entropy typical of Lagrangian schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Finding Neighbors in a Forest: A b-tree for Smoothed Particle Hydrodynamics Simulations

Finding the exact close neighbors of each fluid element in mesh-free computational hydrodynamical methods, such as the Smoothed Particle Hydrodynamics (SPH), often becomes a main bottleneck for scaling their performance beyond a few million fluid elements per computing node. Tree structures are particularly suitable for SPH simulation codes, which rely on finding the exact close neighbors of each fluid element (or SPH particle). In this work we present a novel tree structure, named \textit{$b$-tree}, which features an adaptive branching factor to reduce the depth of the neighbor search. Depending on the particle spatial distribution, finding neighbors using \tree has an asymptotic best case complexity of $O(n)$, as opposed to $O(n \log n)$ for other classical tree structures such as octrees and quadtrees. We also present the proposed tree structure as well as the algorithms to build it and to find the exact close neighbors of all particles. We assess the scalability of the proposed tree-based algorithms through an extensive set of performance experiments in a shared-memory system. Results show that b-tree is up to $12\times$ faster for building the tree and up to $1.6\times$ faster for finding the exact neighbors of all particles when compared to its octree form. Moreover, we apply b-tree to a SPH code and show its usefulness over the existing octree implementation, where b-tree is up to $5\times$ faster for finding the exact close neighbors compared to the legacy code.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Self-gravitating barotropic equilibrium configurations of rotating bodies with SPH

We present a novel relaxation method to build three-dimensional rotating structures of barotropic bodies using the SPH technique. The method is able to relax gaseous structures in rigid as well as differential rotation. The relaxation procedure strongly relies on the excellent conservation of angular momentum that characterizes the SPH technique. The method has been successfully applied to a variety of zero-temperature white dwarfs and polytropic self-gravitating structures. Our SPH results have been validated by comparing the main features (energies, central densities and the polar to equatorial radius ratio) to those obtained with independent, albeit grid-based methods, as for example, the self-consistent field method, showing that both methods agree within few percents.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Detection of Silent Data Corruptions in Smoothed Particle Hydrodynamics Simulations

Silent data corruptions (SDCs) hinder the correctness of long-running scientific applications on large scale computing systems. Selective particle replication (SPR) is proposed herein as the first particle-based replication method for detecting SDCs in Smoothed particle hydrodynamics (SPH) simulations. SPH is a mesh-free Lagrangian method commonly used to perform hydrodynamical simulations in astrophysics and computational fluid dynamics. SPH performs interpolation of physical properties over neighboring discretization points (called SPH particles) that dynamically adapt their distribution to the mass density field of the fluid. When a fault (e.g., a bit-flip) strikes the computation or the data associated with a particle, the resulting error is silently propagated to all nearest neighbors through such interpolation steps. SPR replicates the computation and data of a few carefully selected SPH particles. SDCs are detected when the data of a particle differs, due to corruption, from its replicated counterpart. SPR is able to detect many DRAM SDCs as they propagate by ensuring that all particles have at least one neighbor that is replicated. The detection capabilities of SPR were assessed through a set of error-injection and detection experiments and the overhead of SPR was evaluated via a set of strong-scaling experiments conducted on an HPC system. The results show that SPR achieves detection rates of 91-99.9%, no false-positives, at an overhead of 1-10%.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Equalizing resolution in smoothed-particle hydrodynamics calculations using self-adaptive sinc kernels

The smoothed-particle hydrodynamics (SPH) technique is a numerical method for solving gas-dynamical problems. It has been applied to simulate the evolution of a wide variety of astrophysical systems. The method has a second-order accuracy, with a resolution that is usually much higher in the compressed regions than in the diluted zones of the fluid. In this work, we propose and check a scheme to balance and equalize the resolution of SPH between high- and low-density regions. This method relies on the versatility of a family of interpolators called Sinc kernels, which allows increasing the interpolation quality by varying only a single parameter (the exponent of the Sinc function). The scheme is checked and validated through a number of numerical tests, from standard one-dimensional Riemann problems in shock tubes, to multidimensional simulations of explosions, hydrodynamic instabilities and the collapse of a sun-like polytrope. The analysis of the hydrodynamical simulations suggests that the scheme devised to equalizing accuracy improves the treatment of the post-shock regions and, in general, of the rarefacted zones of fluids while causing no harm to the growth of hydrodynamic instabilities. The method is robust and easy to implement with a low computational overload. It conserves mass, energy, and momentum and reduces to the standard SPH scheme in regions of the fluid that have smooth density gradients.

0 citations0 reviewsNo code linkedOpen access