# boundary-repulsion

**Repository Path**: fortran-sph/boundary-repulsion

## Basic Information

- **Project Name**: boundary-repulsion
- **Description**: 🕵️‍♂️ 对 SPH 的固壁边界排斥力进行研究
- **Primary Language**: FORTRAN
- **License**: MIT
- **Default Branch**: master
- **Homepage**: None
- **GVP Project**: No

## Statistics

- **Stars**: 0
- **Forks**: 0
- **Created**: 2022-07-04
- **Last Updated**: 2022-07-04

## Categories & Tags

**Categories**: Uncategorized

**Tags**: None

## README

# 边界排斥力

- 如果固壁边界存在排斥力，若没有与流体粒子相耦合计算，则初始时刻，流体粒子将受到流体域的压力合力，导致粒子向固壁移动，
再受到极大的排斥力，将有可能向流体域快速发射。

- 如果固壁边界没有排斥力，而是采用固壁虚粒子与流体粒子耦合计算，固壁边界粒子初始化间距若明显小于流体粒子（即更密），
则固壁边界对流体粒子的合力将大于流体域粒子对流体粒子的合力，流体将有可能有向内流动的趋势。

以上情况，都将对粒子初始分布进行不必要的扰动！这是本库的**研究背景**。

```text
# 二维粒子分布，假设 1:5 是固壁虚粒子，6:9 是流体实粒子，进行粒子 6 的压强求和
  5     9
  4
  3  6  8
  2
  1     7
```

```powershell
>> fpm run  # 使用 fpm 运行代码
           1   1.6550342035820654E-003   1.6550342035820654E-003
           2   4.7456099237337653E-003   2.3728049618668826E-003
           3   6.8209261325098006E-003   0.0000000000000000
           4   4.7456099237337653E-003  -2.3728049618668826E-003
           5   1.6550342035820654E-003  -1.6550342035820654E-003
           6   0.0000000000000000        0.0000000000000000
           7  -1.6550342035820654E-003   1.6550342035820654E-003
           8  -6.8209261325098006E-003   0.0000000000000000
           9  -1.6550342035820654E-003  -1.6550342035820654E-003
   9.4912198474675306E-003   0.0000000000000000
```

**研究结果：**

可能的方案是，既设立边界排斥力（但会增加计算量），也设立实虚粒子耦合计算，实虚粒子的粒子间距大致相同，固壁虚粒子可以设置 2 层粒子。