SPHinXsys
alpha version

SPHinXsys (pronunciation: s'finksis) is an acronym from Smoothed Particle Hydrodynamics for industrial compleX systems. It provides C++ APIs for physical accurate simulation and aims to model coupled industrial dynamic systems including fluid, solid, multibody dynamics and beyond with SPH (smoothed particle hydrodynamics), a meshless computational method using particle discretization.
Fluid dynamics, solid dynamics, fluidstructure interactions (FSI), and their coupling to multibody dynamics (with SIMBody library https://simtk.org)
SPH is a fully Lagrangian particle method, in which the continuum media is discretized into Lagrangian particles and the mechanics is approximated as the interaction between them with the help of a kernel, usually a Gaussianlike function. SPH is a mesh free method, which does not require a mesh to define the neighboring configuration of particles, but construct of update it according to the distance between particles. A remarkable feature of this method is that its computational algorithm involves a large number of common abstractions which link to many physical systems inherently. Due to such unique feature, SPH have been used here for unified modeling of both fluid and solid mechanics.
The SPH algorithms are based on the published work of the authors. The algorithms for the discretization of the fluid dynamics equations are based on a weakly compressible fluid formulation, which is suitable for the problems with incompressible flows, and compressible flows with low Mach number (less than 0.3). The solid dynamics equations are discretized by a total Lagrangian formulation, which is suitable to study the problems involving linear and nonlinear elastic materials. The FSI coupling algorithm is implemented in a kinematicforce fashion, in which the solid structure surface describes the phaseinterface and, at the same time, experiences the surface forces imposed by the fluid pressure and friction.
2D models can be built using basic shapes (polygon and circle) and full version of binary operations. 3D models can be generated by simple shapes (brick and sphere), imported from external STL files and processed by applying simple binary operations, e.g. add and substract.
Newtonian fluids with isothermal linear equation of state. Nonnewtonian fluids with OldroydB model. Linear elastic solid, nonlinear elastic solid with NeoHookian model and anisotropic muscle model.
Uniform resolution is used within each fluid or solid bodies. However, it is allowed to use different resolutions for different bodies. For example, one is able to using higher resolution for a solid body which is interacting with a fluid body with lower resolution.
Intel Threading Building Blocks (TBB) is used for the multicore parallelism.
Xiangyu Hu, Luhui Han, Chi Zhang, Shuoguo Zhang, Massoud Rezavand, Yongchuan Yu
Xiangyu Hu (xiang), Department of Mechanical Engineering, Technical University of Munich yu.h u@tum .de
German Research Foundation (Deutsche Forschungsgemeinschaft) DFG HU1527/61, HU1527/101 and HU1527/121.
Luhui Han and Xiangyu Hu, "SPH modeling of fluidstructure interaction", Journal of Hydrodynamics, 2018: 30(1):6269.
Chi Zhang and Massoud Rezavand and Xiangyu Hu, "Dualcriteria time stepping for weakly compressible smoothed particle hydrodynamics", arXiv:1905.12302