JOSEPHINE: A parallel SPH code for free-surface flows

Cherfils, Jean-Marc; Pinon, Grégory; Rivoalen, Elie

doi:10.1016/j.cpc.2012.02.007

Cited by 48 publications

(17 citation statements)

References 26 publications

Supporting

Mentioning

Contrasting

Unclassified

Order By: Relevance

“…The most straightforward method of spatial decomposition is a partitioning of a computational domain into equal rectangular subdomains [6,7]. In the context of spatial nonuniformity of particle concentration, the partitioning can be performed recursively, using the orthogonal recursive bisection (ORB) method.…”

Section: Introductionmentioning

confidence: 99%

Parallel SPH modeling using dynamic domain decomposition and load balancing displacement of Voronoi subdomains

Егорова

Dyachkov

Паршиков

et al. 2019

Computer Physics Communications

View full text Add to dashboard Cite

A highly adaptive load balancing algorithm for parallel simulations using particle methods, such as molecular dynamics and smoothed particle hydrodynamics (SPH), is developed. Our algorithm is based on the dynamic spatial decomposition of simulated material samples between Voronoi subdomains, where each subdomain with all its particles is handled by a single computational process which is typically run on a single CPU core of a multiprocessor computing cluster.The algorithm displaces the positions of neighbor Voronoi subdomains in accordance with the local load imbalance between the corresponding processes. It results in particle transfers from heavy-loaded processes to less-loaded ones. Iteration of the algorithm puts into alignment the processor loads. Convergence to a well-balanced decomposition from imbalanced one is improved by the usage of multi-body terms in the balancing displacements.The high adaptability of the balancing algorithm to simulation conditions is illustrated by SPH modeling of the dynamic behavior of materials under extreme conditions, which are characterized by large pressure and velocity gradients, as a result of which the spatial distribution of particles varies greatly in time. The higher parallel efficiency of our algorithm in such conditions is demonstrated by comparison with the corresponding static decomposition of the computational domain. Our algorithm shows almost perfect strong scalability in tests using from tens to several thousand processes.

show abstract

Section: Introductionmentioning

confidence: 99%

Parallel SPH modeling using dynamic domain decomposition and load balancing displacement of Voronoi subdomains

Егорова

Dyachkov

Паршиков

et al. 2019

Computer Physics Communications

View full text Add to dashboard Cite

show abstract

“…As for the treatment of polygon walls, mirror particle is a choice. Mirror particle can obtain smooth pressure distribution in two dimensions.…”

Section: Introductionmentioning

confidence: 99%

Improvement of pressure distribution to arbitrary geometry with boundary condition represented by polygons in particle method

Zhang

Koshizuka

Murotani

et al. 2017

Numerical Meth Engineering

View full text Add to dashboard Cite

The boundary condition represented by polygons in the moving particle semi-implicit method can accurately represent geometries and treat complex geometry with high efficiency. However, inaccurate wall contribution to the Poisson's equation leads to drastic numerical oscillation. To address this issue, in this research, we analyzed the problems of the Poisson's equation used in the boundary condition represented by polygons. The new Poisson's equation is proposed based on the improved source term (Tanaka and Masunaga, Trans Jpn Soc Comput Eng Sci, 2008). The asymmetric gradient model (Khayyer and Gotoh, Coastal Engineering Journal, 2008) is also adopted to further suppress the numerical oscillation of fluid particles. The proposed method can dramatically improve the pressure distribution to arbitrary geometry in three dimensions and keep the efficiency. Four examples including the hydrostatic simulation, dam break simulation, and two complex geometries are verified to show the general applicability of the proposed method. Figure 17. Snapshots of the pressure distribution in a water tank with a wedge at t = 0.5, 1.0, 1.5, 2.0, 2.5, 3.0 s. Figure 20. Snapshots of the complex geometry at t = 0.05, 0.09, 0.14, 0.16, and 0.18 s.

show abstract

“…C'est sur ce point que nous focalisons nos travaux de modélisation. L'outil numérique JOSEPHINE (CHERFILS et al, 2012) que nous utilisons est un code de type SPH. C'est une méthode lagrangienne, sans maillage, bien adaptée au suivi des mouvements d'interface.…”

Section: Introductionunclassified

Modélisation par la méthode SPH de l'impact de la houle sur une paroi verticale

Lu¹,

Cherfils²,

Pinon³

et al. 2014

XIIIèmes JNGCGC, Dunkerque

View full text Add to dashboard Cite

Résumé :Un outil numérique, de type SPH (Smoothed Particle Hydrodynamics), a été développé dans le but de simuler l'interaction hydrodynamique entre la houle et des ouvrages portuaires ou côtiers, dotés de systèmes de récupération d'énergie marine. Nous nous intéressons plus particulièrement à la prédiction des pics de pression survenant lors de l'impact de vagues sur un ouvrage, afin de calculer les efforts auxquels la structure sera soumise. L'impact d'un soliton sur un mur vertical, adossé à une plage inclinée, a donc été simulé dans un canal à houle numérique, montrant ainsi la capacité du code à traiter ce type de problème. En faisant varier les paramètres d'entrée du batteur, nous avons retrouvé les différentes formes d'impacts mis en évidence par les expériences : aerated impact, flip through, air pocket et slosh impact. Afin d'étudier le lien entre la forme de l'interface et la dynamique de l'impact, un cas académique de "jet triangulaire" a été étudié de manière similaire à l'étude de . Plusieurs variantes de la méthode SPH (Rusanov-SPH, Delta-SPH, …) ont été testées en termes de stabilité et précision. Cela a conduit à une nette amélioration du calcul de la pression sur la paroi et de la trajectoire du jet tout en préservant le caractère lagrangien de la méthode. D'autres variantes de la méthode SPH, basées sur les schémas type ALE (Arbitraire Lagrange Euler) et les solveurs de Riemann permettent d'obtenir de bons résultats (GUILCHER et al., 2010, RAFIEE et al., 2013, moyennant des temps de calculs supérieurs. Les résultats obtenus ici sont comparés avec ceux d'une autre méthode numérique, de type VOF, ainsi qu'avec un modèle analytique.

show abstract

JOSEPHINE: A parallel SPH code for free-surface flows

Cited by 48 publications

References 26 publications

Parallel SPH modeling using dynamic domain decomposition and load balancing displacement of Voronoi subdomains

Parallel SPH modeling using dynamic domain decomposition and load balancing displacement of Voronoi subdomains

Improvement of pressure distribution to arbitrary geometry with boundary condition represented by polygons in particle method

Modélisation par la méthode SPH de l'impact de la houle sur une paroi verticale

Contact Info

Product

Resources

About