《Computers & Geosciences》：Enhanced rapid simulation of debris flow using game engine-based SPH-DEM methodology
编辑推荐：
理解泥石流的形成与演化对于有效灾害防治至关重要。虚拟现实提供了直观的可视化方式展示其启动与发展过程，然而现有虚拟仿真研究通常采用SPH–SPH模型来同时表示固相和液相。但研究人员已知SPH不适用于固体建模，因为其核函数缺乏断裂的物理基础，难以捕捉结构破坏和颗粒
理解泥石流的形成与演化对于有效灾害防治至关重要。虚拟现实提供了直观的可视化方式展示其启动与发展过程，然而现有虚拟仿真研究通常采用SPH–SPH模型来同时表示固相和液相。但研究人员已知SPH不适用于固体建模，因为其核函数缺乏断裂的物理基础，难以捕捉结构破坏和颗粒尺度破碎。使用SPH描述两相还需要复杂的多相扩展、额外的本构定律以及专门的接触-摩擦处理，这显著增加了建模难度。同时，实时三维可视化解决方案仍然有限。为应对这些挑战，本研究提出了一种紧密集成物理建模与可视化的虚拟仿真框架。在力学模拟（可视化之前的物理计算）方面，采用SPH–DEM耦合模型模拟泥石流运动。为克服邻域搜索的性能瓶颈，设计了一种混合哈希网格–八叉树加速结构，在密度计算中效率提升38.7%，在压力和粘度模块中效率提升34.4%，尤其在密集粒子区域。为缓解SPH边界密度低估问题，引入了一种基于实时密度和速度梯度的自适应补偿策略。采用经典两相溃坝实验验证所提模型，结果显示前缘位移平均误差为8.552%。敏感性分析表明，在关键参数扰动下，前缘传播保持稳定并表现出物理上合理的响应特性。验证结果表明，实现的SPH–DEM求解器可在游戏引擎环境中稳定运行，并显示出初步的数值可行性。未来工作将纳入更复杂的实验场景和实际案例，以进一步评估模拟精度并增强计算稳定性。在可视化方面，该框架集成了屏幕空间流体渲染（SSF），并通过双边滤波增强表面质量。实验在配备NVIDIA RTX 3070 Ti GPU的平台上进行。结果表明，在多达20,000个粒子的情况下，系统保持实时性能，同时在物理精度和视觉逼真度之间实现了合理平衡。总体而言，本研究为在游戏引擎环境中集成泥石流模拟与可视化的可行性提供了初步证据。
**论文解读文章**
**研究背景与问题**
泥石流是一种突发性高、破坏性极强的地质灾害，对山区基础设施和居民安全构成严重威胁。深入理解其形成机制与演化过程，是揭示成灾机理、提高预测精度及制定有效防控策略的基础。传统研究主要依赖水槽实验和野外观测，但这些方法受环境条件、高昂成本和安全风险限制，难以全面捕捉泥石流的动态特征。为克服这些局限，研究者转向数值模拟方法探索泥石流动力学过程，并在历史灾害重建和减灾评估中取得显著成果。当前，多种流固耦合数值模型被提出，按离散化方式分为欧拉-欧拉、欧拉-拉格朗日和拉格朗日-拉格朗日方法。其中，拉格朗日-拉格朗日模型（如SPH–FEM、SPH–SPH和SPH–DEM）在处理大变形、自由表面和复杂流-粒相互作用方面更具灵活性。然而，数值模拟固有地缺乏直观的三维表示、交互性有限且对实时外部数据响应弱。虚拟仿真技术整合计算机图形学和交互技术，可在三维环境中再现流动动力学、碰撞及相关过程，并实现实时交互。现有虚拟仿真研究多采用SPH–SPH模型描述固液相，但SPH不适合固体建模，因其核函数缺乏断裂物理基础，难以捕捉结构破坏和颗粒尺度破碎；同时两相描述需要复杂多相扩展和特殊接触摩擦处理，增加建模难度。此外，实时三维可视化方案仍然有限。
**研究目标与意义**
为应对上述挑战，本研究提出一种紧密集成物理建模与可视化的虚拟仿真框架。物理模拟采用SPH–DEM耦合模型更准确捕获泥石流中固液相互作用；可视化引入屏幕空间流体渲染（SSF）技术。该研究初步证明了在游戏引擎环境中集成泥石流模拟与可视化的可行性，为虚拟仿真在灾害防治中的应用提供了技术参考。论文发表在《Computers》。
**关键技术与方法**
1. **混合哈希网格–八叉树加速结构**：针对邻域搜索性能瓶颈，设计混合哈希网格–八叉树结构，在密度计算中效率提升38.7%，在压力和粘度模块提升34.4%，尤其适用于密集粒子区域。
2. **自适应边界密度补偿策略**：基于实时密度和速度梯度，对边界粒子进行识别和密度补偿，缓解SPH边界密度低估问题，提升系统稳定性和边界精度。
3. **屏幕空间流体渲染（SSF）与双边滤波**：利用SSF实现流体相实时渲染，并对深度图应用双边滤波进行边缘保持平滑，增强复杂边界特征的视觉表现。
4. **游戏引擎集成**：基于Unreal Engine 5.3.2平台，利用Niagara粒子系统实现SPH–DEM求解器，结合静态网格渲染模块处理固相。
**研究结果**
**2.1 SPH-DEM动态模型**
建立了SPH与DEM的双向力学耦合机制：SPH累积流体对固体的总力，DEM固体粒子插值到邻域流体粒子的SPH方程中计算反作用力；同时处理两相间的接触力和摩擦力，以及固体间的碰撞力。
**2.2 SPH-DEM模型的实现**
在Unreal Engine的Niagara粒子系统中实现SPH–DEM模型，利用其数据驱动、基于节点的流程实现每个粒子属性的精细控制和更新规则定制。
**2.3 泥石流可视化实现**
可视化包括四个核心模块：固体渲染、深度纹理计算、法线纹理计算和材质生成。固体模块通过实例化预设岩石模型进行网格渲染；流体模块使用SSF技术，结合双边滤波优化深度图，最终融合两相渲染结果。
**2.4 加速结构性能评估**
混合哈希网格–八叉树结构在密集粒子区域显著提升效率：密度计算加速38.7%，压力和粘度模块加速34.4%。
**2.5 密度补偿验证**
自适应密度补偿策略有效缓解边界密度低估，通过溃坝实验验证前缘位移平均误差为8.552%，敏感性分析表明关键参数扰动下前缘传播稳定且物理响应合理。
**2.6 整体系统性能**
在NVIDIA RTX 3070 Ti GPU平台上，20,000个粒子时系统保持实时性能，物理精度与视觉逼真度取得合理平衡。
**讨论与结论**
**讨论部分**：研究人员指出，Liu（2022）在SPH-DEM框架中引入结构动力学模块并开发SPH-DEM-FEM三相耦合模型，可详细模拟泥石流与防护结构相互作用，为失效机制和设计优化提供重要见解。相比之下，本研究尚未纳入真实的泥石流-结构相互作用，限制了模型在物理精度和系统完整性方面的表现。未来需扩展耦合模型以支持结构交互，并进一步评估更复杂场景下的模拟精度和计算稳定性。
**研究结论**（翻译自Conclusion部分，原文未单独列出，根据讨论和摘要提炼）：本研究提出了一种基于游戏引擎的SPH-DEM泥石流模拟与可视化集成框架，通过混合邻域搜索加速结构和自适应边界密度补偿策略，实现了高效稳定的物理模拟；结合屏幕空间流体渲染技术，在实时性能下提供了合理的视觉表现。实验验证了模型在经典两相溃坝中的前缘位移误差为8.552%，且对参数扰动表现出稳定性。该工作为在游戏引擎环境中进行泥石流虚拟仿真提供了初步可行性证据，未来将纳入更多复杂实验场景和实际案例以进一步验证精度和稳定性。
