该项目实现了基于平滑粒子流体动力学(SPH)方法的经典二维溃坝数值模拟。项目通过在MATLAB环境中构建完整的拉格朗日粒子体系,能够精确捕捉流体在突发性溃堤后的非线性大变形过程以及自由表面的演化细节。其核心功能涵盖了流体粒子与边界粒子的初始化配置、基于三次样条核函数的空间场变量估算、以及针对纳维-斯托克斯方程的无网格化数值求解。
实现方法上,系统采用显式时间积分方案更新粒子的位置与速度。通过求解连续性方程更新各粒子的密度,并利用状态方程(Equation of State)计算压力,进而获得压力梯度力。为