以计算流体动力学模型（CFD）模拟的戈壁地表风沙两相流运动特征

尘暴过程中沙砾质戈壁地表的风沙两相流运动，一直是风沙物理学关注的问题。但由于技术限制，野外观测及风洞实验难以捕捉大规模颗粒运动。为此，本文通过构建流体-颗粒碰撞互馈模型进行三维场域数值模拟，分析了不同砾石盖度下尘暴来流中的风沙两相流运动规律。结果表明：基于三维场流体-颗粒碰撞互馈模型修正后的求解结果能够较为准确地模拟戈壁地表风沙两相流运动过程。(1)不同来流尘暴过程，风相速度均会随砾石盖度的增大而减小，并在砾石盖度大于40%的戈壁近地表形成拜格诺焦点。相较于风相，沙相降速不明显。此外，距砾石床面4 cm以上，风沙两相速度基本保持一致。(2)当尘暴速度大于10 m·s^-1时，在距地1.5～2.1 m高度形成风沙两相流速度交点，交点上下表现为“风快沙慢”及“风慢沙快”现象。(3)尘暴过程中近地表沙尘浓度沿均值线呈“脉动”规律，表现为地表粉尘沉降与再释放循环过程，同时砾石盖度大小会直接影响沙尘沉积分布，当盖度大于40%时能够有效抑制沙尘释放。