Physics-driven complex relaxation for multi-body systems of SPH method

Chenxi Zhao,Yongchuan Yu,Oskar J. Haidn,Xiangyu Hu

2024-09-26

Abstract:In the smoothed particle dynamics (SPH) method, the characteristics of a target particle are interpolated based on the information from its neighboring particles. Consequently, a uniform initial distribution of particles significantly enhances the accuracy of SPH calculations. This aspect is particularly critical in Eulerian SPH, where particles are stationary throughout the simulation. To address this, we introduce a physics-driven complex relaxation method for multi-body systems. Through a series of two-dimensional and three-dimensional case studies, we demonstrate that this method is capable of achieving a globally uniform particle distribution, especially at the interfaces between contacting bodies, and ensuring improved zero-order consistency. Moreover, the effectiveness and reliability of the complex relaxation method in enhancing the accuracy of physical simulations are further validated.

Computational Engineering, Finance, and Science

What problem does this paper attempt to address?

This paper attempts to solve the problem of how to obtain a uniform and isotropic initial particle distribution when using the Smoothed Particle Hydrodynamics (SPH) method in multi - body systems, especially in multi - body systems in the Eulerian framework. Specifically, the paper focuses on the following points:

1. **Ensuring zero - order consistency**: At the interface of multi - body systems, especially for complex geometries, ensuring zero - order consistency of particle distribution (i.e., the balance of forces between particles and the position no longer changes with time) is crucial. This is very important for improving the accuracy and reliability of physical simulations.

2. **Eliminating gaps between multi - body systems**: When importing different geometries, due to different parsers and discrete origin settings, gaps may occur between multi - body systems. These gaps will affect the accuracy of the simulation and may even lead to calculation crashes, especially in the Eulerian framework.

3. **Generating high - quality particle distributions**: Many existing methods have difficulty ensuring the quality of particle distribution near the boundary and the overall isotropy when dealing with complex geometries. For example, methods using virtual particles or manually arranging boundary particles have limitations.

To solve these problems, the author proposes a physically - driven complex relaxation method. The main features of this method include:

- **Boolean operations to ensure no gaps**: Through Boolean operations based on the level - set field, ensure that there are no gaps between multi - body systems.
- **Simultaneously relaxing solid and fluid particles**: In each time step, perform a physically - driven relaxation process on the internal solid particles and the external fluid particles simultaneously.
- **Improved boundary treatment**: Combine the surface bounding method and the static confinement method to ensure the position and distribution quality of boundary particles.

Through these improvements, this method can achieve better zero - order consistency and uniform particle distribution at the interface of multi - body systems, thereby improving the accuracy and efficiency of SPH simulations in the Eulerian framework.

### Formula Summary

1. **Definition of level - set field**:
\[
\Gamma=\{(x, y)\mid\phi(x, y, t) = 0\}
\]
\[
N=\frac{\nabla\phi}{|\nabla\phi|}
\]

2. **Physically - driven particle relaxation**:
\[
\frac{dv}{dt}=F_p
\]
\[
F_{p,i}=-\frac{2p_0 V_i}{m_i}\sum_j\nabla_i W_{ij}V_j
\]
\[
\Delta t\leq0.25\frac{s}{h|\frac{dv}{dt}|}
\]
\[
r^{n + 1}=r^n+\frac{1}{2}F_p^n(\Delta t)^2
\]

3. **Surface bounding method**:
\[
r_i =
\begin{cases}
r_i-\left(\phi_i+\frac{1}{2}\Delta x\right)N_i & \text{if }\phi_i\geq-\frac{1}{2}\Delta x \\
r_i & \text{otherwise}
\end{cases}
\]

4. **Static confinement method**:
\[
F_{p,i}=-\frac{2p_0 V_i}{m_i}\left(\sum_{j,\phi_j < 0}\nabla W_{ij}V_j+\sum_{k,\phi_k > 0}H(\phi_k,\epsilon)l_f^m\nabla_i W_{ik}\right)
\]
\[
H(\phi,