SPH 与 DEM 方法系统对比
@[toc0]
SPH 与 DEM 方法系统对比：异同与耦合应用
一、核心本质差异
| 维度 | SPH (Smoothed Particle Hydrodynamics) | DEM (Discrete Element Method) |
|---|---|---|
| 物理本质 | 连续介质近似：粒子代表连续流体/固体的"质量点"，通过核函数重构场变量 | 离散实体模拟：粒子代表真实物理颗粒（沙粒、药丸、矿石等），具有明确几何边界 |
| 数学基础 | 基于积分-微分近似： ⟨ A ( r i ) ⟩ = ∑ j m j A j ρ j W ( r i − r j , h ) \displaystyle \langle A(\mathbf{r}_i) \rangle = \sum_j m_j \frac{A_j}{\rho_j} W(\mathbf{r}_i - \mathbf{r}_j, h)⟨A(ri)⟩=j∑mjρjAjW(ri−rj,h) | 基于牛顿第二定律 + 接触力学： m i r ¨ i = ∑ j ∈ N i F i j contact + F i body \displaystyle m_i \ddot{\mathbf{r}}_i = \sum_{j \in \mathcal{N}_i} \mathbf{F}_{ij}^{\text{contact}} + \mathbf{F}_i^{\text{body}}mir¨i=j∈Ni∑Fijcontact+Fibody |
| 相互作用 | 非局部平滑作用：通过核函数W WW在光滑长度h hh内连续加权平均 | 局部接触作用：仅当粒子几何重叠/接触时产生离散力（需显式接触检测） |
| 守恒性质 | 天然满足动量/能量守恒（对称核函数设计） | 需显式处理碰撞中的能量耗散（恢复系数e ee控制） |
| 典型应用 | 自由表面流、溃坝、大变形固体、天体物理 | 颗粒流、粉体混合、岩石破碎、料仓卸料 |
💡关键区别：
- SPH粒子 ≠ 物理颗粒：1个SPH粒子代表连续介质的一小块（如1 mm³水），可穿透/重叠
- DEM粒子 = 物理颗粒：1个DEM粒子对应1个真实颗粒（如1颗沙粒），不可穿透（需接触力阻止）
二、方法论对比详表
| 特性 | SPH | DEM |
|---|---|---|
| 框架类型 | 拉格朗日无网格方法（Meshless） | 离散粒子动力学（Discrete） |
| 控制方程 | Navier-Stokes / 弹性动力学方程的粒子离散化 | 牛顿运动方程 + 接触力模型（Hertz-Mindlin, Linear Spring-Dashpot） |
| 时间积分 | 通常显式（Leapfrog, Verlet） | 通常显式（Velocity-Verlet），受接触刚度限制（Δ t ∝ m / k n \Delta t \propto \sqrt{m/k_n}Δt∝m/kn） |
| 边界处理 | 虚粒子法、边界积分法、镜像粒子 | 刚性/柔性壁面建模（同颗粒接触逻辑） |
| 计算复杂度 | O ( N log N ) O(N \log N)O(NlogN)（需邻域搜索） | O ( N 2 ) O(N^2)O(N2)理论最坏，实际O ( N ) O(N)O(N)（通过空间分割优化） |
| 并行扩展性 | 优秀（纯拉格朗日，GPU友好） | 良好（接触检测需通信，但商业软件如EDEM已支持大规模并行） |
| 开源实现 | DualSPHysics, PySPH | LIGGGHTS, YADE |
三、典型应用场景对比
✅SPH 优势场景
- 自由表面流动：溃坝、波浪冲击、液滴飞溅（无需界面捕捉算法）
- 大变形固体：金属成型、弹道冲击（避免网格畸变）
- 多相流：气-液界面（通过密度/状态方程区分相）
- 天体物理：星系形成、超新星爆发（GADGET系列）
✅DEM 优势场景
- 颗粒系统：沙土流动、药片混合、谷物流动
- 破碎/断裂：岩石爆破、脆性材料失效（通过bond breaking建模）
- 摩擦主导过程：料仓堵塞、颗粒筛分
- 工业设备：搅拌器、输送带、粉碎机（EDEM广泛用于工程设计）
⚠️易混淆场景
| 问题 | 正确方法 | 常见误区 |
|---|---|---|
| 模拟"沙子流动" | DEM（沙粒为离散实体） | 误用SPH（将沙视为连续体，丢失颗粒尺度效应） |
| 模拟"水冲击沙堆" | SPH-DEM耦合（水用SPH，沙用DEM） | 单独用SPH（无法捕捉颗粒级渗流）或单独用DEM（无法模拟流体） |
| 模拟"泥石流" | SPH-DEM耦合（泥浆用SPH，石块用DEM） | 单独用SPH（丢失大石块惯性效应） |
四、SPH-DEM 耦合：优势互补
当问题同时涉及连续流体和离散颗粒时，需耦合两种方法：
耦合机制
典型耦合策略
| 耦合类型 | 描述 | 应用案例 |
|---|---|---|
| 单向耦合 | 流体→颗粒（忽略颗粒对流体反作用） | 低浓度悬浮颗粒（< 1 % <1\%<1%体积分数） |
| 双向耦合 | 流体↔颗粒动量交换 | 泥石流、高浓度浆料输送 |
| 解析尺度 | 颗粒尺寸 ≫ 流体网格/粒子 | 颗粒周围流场显式解析（计算昂贵） |
| 非解析尺度 | 颗粒尺寸 ≈ 流体粒子 | 通过阻力模型（如Gidaspow）近似流体作用力 |
代表性研究
- 泥石流模拟：SPH模拟泥浆基质，DEM模拟大石块，成功复现抛射、爬高现象 [[2]]
- 抛石加固：DEM模拟落石运动，SPH模拟水动力，研究管道防护 [[5]]
- 磨料水射流：SPH模拟水射流，DEM模拟磨料颗粒，预测材料侵蚀 [[4]]
- 气泡-颗粒相互作用：SPH模拟气泡上升，DEM模拟颗粒沉降，研究三相流 [[7]]
五、选择指南：何时用SPH？何时用DEM？
决策矩阵
| 问题特征 | 推荐方法 | 理由 |
|---|---|---|
| 自由表面水流 | SPH | 天然处理界面，无需VOF/Level Set |
| 沙土流动（干燥） | DEM | 颗粒摩擦/碰撞主导，连续假设失效 |
| 湿沙流动 | SPH-DEM耦合 | 水（连续）+ 沙（离散）需同时建模 |
| 金属3D打印熔池 | SPH | 高温流体+大变形，网格法易畸变 |
| 药片混合机 | DEM | 颗粒形状/摩擦关键，工业标准方法 |
| 水下爆炸冲击结构 | SPH-DEM | 水（SPH）+ 破碎结构（DEM bonded particles）[[10]] |
六、前沿挑战与发展趋势
| 挑战 | SPH | DEM | 耦合方法 |
|---|---|---|---|
| 多尺度问题 | 微观湍流需LES/DES | 颗粒团聚需多尺度接触模型 | 流体-颗粒尺度分离（如CFD-DEM vs SPH-DEM） |
| 计算效率 | 邻域搜索瓶颈（KD-Tree/Hash Grid优化） | 接触检测瓶颈（GPU加速接触算法） | 耦合接口通信开销（需异步更新策略） |
| 精度提升 | δ-SPH抑制振荡、Riemann-SPH提升界面 | 非球形颗粒（Superquadric, Polyhedra） | 全解析耦合（Resolved coupling）vs 非解析（Unresolved） |
| 新兴方向 | 机器学习增强（diffSPH参数优化） | 数字孪生（实时DEM仿真） | 多物理场耦合（SPH-DEM-FEM热-力-流耦合）[[2]] |
🔍实践建议：
- 纯流体问题→ 优先选 SPH（自由表面）或传统CFD（内部流）
- 纯颗粒问题→ 优先选 DEM（LIGGGHTS/YADE开源，EDEM商业）
- 流-固耦合→ 根据颗粒浓度选择：
- 低浓度：CFD-DEM（OpenFOAM + LIGGGHTS via CFDEMcoupling）
- 高浓度/自由表面：SPH-DEM（DualSPHysics内置DEM模块）
参考资源
- SPH-DEM耦合综述：Trujillo-Vela et al. (2020),Computers and Geotechnics[[35]]
- 开源耦合框架：
- DualSPHysics（内置DEM模块）
- CFDEMcoupling（OpenFOAM + LIGGGHTS）
- 验证案例库：SPHERIC Benchmark（含SPH-DEM耦合测试）