颠覆“孝顺就是顺从”，建立需求边界模型，颠覆道德绑架，输出既尽孝又不内耗的相处方案。

[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​ρj​Aj​​W(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}}mi​r¨i​j∈Ni​∑​Fijcontact​Fibody​相互作用非局部平滑作用通过核函数W WW在光滑长度h hh内连续加权平均局部接触作用仅当粒子几何重叠/接触时产生离散力需显式接触检测守恒性质天然满足动量/能量守恒对称核函数设计需显式处理碰撞中的能量耗散恢复系数e ee控制典型应用自由表面流、溃坝、大变形固体、天体物理颗粒流、粉体混合、岩石破碎、料仓卸料关键区别SPH粒子 ≠ 物理颗粒1个SPH粒子代表连续介质的一小块如1 mm³水可穿透/重叠DEM粒子 物理颗粒1个DEM粒子对应1个真实颗粒如1颗沙粒不可穿透需接触力阻止

方法论对比详表特性SPHDEM框架类型拉格朗日无网格方法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^

O(N

理论最坏实际O ( N ) O(N)O(N)通过空间分割优化并行扩展性优秀纯拉格朗日GPU友好良好接触检测需通信但商业软件如EDEM已支持大规模并行开源实现DualSPHysics, PySPHLIGGGHTS, YADE

典型应用场景对比✅SPH 优势场景自由表面流动溃坝、波浪冲击、液滴飞溅无需界面捕捉算法大变形固体金属成型、弹道冲击避免网格畸变多相流气-液界面通过密度/状态方程区分相天体物理星系形成、超新星爆发GADGET系列✅DEM 优势场景颗粒系统沙土流动、药片混合、谷物流动破碎/断裂岩石爆破、脆性材料失效通过bond breaking建模摩擦主导过程料仓堵塞、颗粒筛分工业设备搅拌器、输送带、粉碎机EDEM广泛用于工程设计⚠️易混淆场景问题正确方法常见误区模拟沙子流动DEM沙粒为离散实体误用SPH将沙视为连续体丢失颗粒尺度效应模拟水冲击沙堆SPH-DEM耦合水用SPH沙用DEM单独用SPH无法捕捉颗粒级渗流或单独用DEM无法模拟流体模拟泥石流SPH-DEM耦合泥浆用SPH石块用DEM单独用SPH丢失大石块惯性效应

SPH-DEM 耦合优势互补当问题同时涉及连续流体和离散颗粒时需耦合两种方法耦合机制Drag ForcePressure GradientMomentum ExchangeVolume DisplacementDEM_Phase接触检测法向/切向力滚动摩擦SPH_Phase核函数近似动量方程状态方程SPH Fluid ParticlesDEM Solid Particles典型耦合策略耦合类型描述应用案例单向耦合流体→颗粒忽略颗粒对流体反作用低浓度悬浮颗粒 1 % 1\%1%体积分数双向耦合流体↔颗粒动量交换泥石流、高浓度浆料输送解析尺度颗粒尺寸 ≫ 流体网格/粒子颗粒周围流场显式解析计算昂贵非解析尺度颗粒尺寸 ≈ 流体粒子通过阻力模型如Gidaspow近似流体作用力代表性研究泥石流模拟SPH模拟泥浆基质DEM模拟大石块成功复现抛射、爬高现象 [[2]]抛石加固DEM模拟落石运动SPH模拟水动力研究管道防护 [[5]]磨料水射流SPH模拟水射流DEM模拟磨料颗粒预测材料侵蚀 [[4]]气泡-颗粒相互作用SPH模拟气泡上升DEM模拟颗粒沉降研究三相流 [[7]]

选择指南何时用SPH何时用DEM是否否是低浓度高浓度是否是否问题类型涉及流体/大变形固体?涉及离散颗粒?考虑FEM/MPM纯SPH颗粒浓度?SPH 单向耦合SPH-DEM双向耦合问题类型仅涉及颗粒/刚体?颗粒间接触主导?考虑刚体动力学纯DEM考虑CFD-DEM流体用网格法决策矩阵问题特征推荐方法理由自由表面水流SPH天然处理界面无需VOF/Level Set沙土流动干燥DEM颗粒摩擦/碰撞主导连续假设失效湿沙流动SPH-DEM耦合水连续 沙离散需同时建模金属3D打印熔池SPH高温流体大变形网格法易畸变药片混合机DEM颗粒形状/摩擦关键工业标准方法水下爆炸冲击结构SPH-DEM水SPH 破碎结构DEM bonded particles[[10]]

前沿挑战与发展趋势挑战SPHDEM耦合方法多尺度问题微观湍流需LES/DES颗粒团聚需多尺度接触模型流体-颗粒尺度分离如CFD-DEM vs SPH-DEM计算效率邻域搜索瓶颈KD-Tree/Hash Grid优化接触检测瓶颈GPU加速接触算法耦合接口通信开销需异步更新策略精度提升δ-SPH抑制振荡、Riemann-SPH提升界面非球形颗粒Superquadric, Polyhedra全解析耦合Resolved couplingvs 非解析Unresolved新兴方向机器学习增强diffSPH参数优化数字孪生实时DEM仿真多物理场耦合SPH-DEM-FEM热-力-流耦合[[2]]实践建议纯流体问题→ 优先选 SPH自由表面或传统CFD内部流纯颗粒问题→ 优先选 DEMLIGGGHTS/YADE开源EDEM商业流-固耦合→ 根据颗粒浓度选择低浓度CFD-DEMOpenFOAM LIGGGHTS via CFDEMcoupling高浓度/自由表面SPH-DEMDualSPHysics内置DEM模块参考资源SPH-DEM耦合综述Trujillo-Vela et al. (

,Computers and Geotechnics[[35]]开源耦合框架DualSPHysics内置DEM模块CFDEMcouplingOpenFOAM LIGGGHTS验证案例库SPHERIC Benchmark含SPH-DEM耦合测试

www.crm.88888-www.crm.88888最新版N.27.82.12-2285安卓网应用