HPSC Lab at JuliaCon 2024

From July 9th through 13th, JuliaCon 2024 took place at the PSV soccer stadium in Eindhoven, Netherlands. It presented a great opportunity to meet a lot of fellow Julia users and developers and enabled many interesting face-to-face discussions with core Julia developers. It was also nice to meet a number of members of the Trixi Framework development team, such as Hendrik Ranocha (University of Mainz), Erik Faulhaber and Benedict Geihe (both University of Cologne), and Daniel Doehring (RWTH Aachen University).

Niklas (HLRS Stuttgart/HPSC Lab) and Michael (HPSC Lab) were also present and gave talks on particle-based multiphysics simulations with TrixiParticles.jl and on secure numerical computations with fully homomorphic encryption. Overall, it was a great event and we are looking forward to next year!

 

 
The Trixi Framework team at JuliaCon 2024. From left to right: Michael Schlottke-Lakemper, Benedict Geihe, Niklas Neher, Hendrik Ranocha, Daniel Doehring © Daniel Doehring

Together with Erik Faulhaber, Sven Berger, Christian Weißenfels und Gregor Gassner, we have submitted our paper "Robust and efficient pre-processing techniques for particle-based methods including dynamic boundary generation".

 

arXiv:2506.21206 reproduce me!

 

 

Abstract

Obtaining high-quality particle distributions for stable and accurate particle-based simulations poses significant challenges, especially for complex geometries. We introduce a preprocessing technique for 2D and 3D geometries, optimized for smoothed particle hydrodynamics (SPH) and other particle-based methods. Our pipeline begins with the generation of a resolution-adaptive point cloud near the geometry's surface employing a face-based neighborhood search. This point cloud forms the basis for a signed distance field, enabling efficient, localized computations near surface regions. To create an initial particle configuration, we apply a hierarchical winding number method for fast and accurate inside-outside segmentation. Particle positions are then relaxed using an SPH-inspired scheme, which also serves to pack boundary particles. This ensures full kernel support and promotes isotropic distributions while preserving the geometry interface. By leveraging the meshless nature of particle-based methods, our approach does not require connectivity information and is thus straightforward to integrate into existing particle-based frameworks. It is robust to imperfect input geometries and memory-efficient without compromising performance. Moreover, our experiments demonstrate that with increasingly higher resolution, the resulting particle distribution converges to the exact geometry.

Search