基于SPH的风沙运动的数值模拟1-力学学报.doc

基于SPH的风沙运动的数值模拟1) 逯博(((( 买买提明·艾尼* 金阿芳* 徐玉娟* *(新疆大学机械工程学院,乌鲁木齐 830047) 摘要：运用光滑粒子流体动力学（SPH-Smoothed Particles Hydrodynamics）方法对沙粒和气流的相互耦合运动特性进行了分析，研究提出了风沙流的SPH数值方法并进行了数值模拟。首先提出了风沙流的SPH建模方法和基本理论，建立了风沙流动的SPH数值模拟平台。其次通过建立风沙流的SPH模型并施加边界条件，对自然风作用下沙粒的运动情况进行了数值模拟，详细分析了沙粒运动轨迹及特性，最后通过与相关研究成果对比分析，验证了完善后的SPH方法有效性。通过考虑气流场的可变性，在风沙流SPH计算模型中引入了加载（起风）和卸载（停风）方式，观察并对比分析了沙粒的运动轨迹和特性。为进一步研究风沙流的实时动态非线性行为提供了SPH理论基础和数值分析方法。 关键字：SPH，两相流，风沙运动 12-312 引言 沙漠化及其带来的风沙环境问题是人类面临的重大自然灾害之一。那么，找到合理的防止沙漠化的技术就要从研究风沙流动的规律着手。随着计算力学和数值模拟技术的发展，学者们建立了各种各样的数学模型来研究风沙的气固两相流动问题[1~]。本文旨在通过数值模拟的方法对风沙流动的机理和特征进行深入研究。运用一种新的无网格数值模拟方法—SPH（光滑粒子流体动力学）来对风沙流现象进行建模及模拟。气固两相被离散为拉格朗日粒子，流体和固体力学的控制方程则转化为相应的粒子间作用力计算开始后，粒子在这些力的作用下运动演化由粒子携带的信息，通过插值核函数求得流场和固体的运动参数。 1 风沙流动的基本方程 风沙流动属于典型的两相流问题，要考虑风场和沙粒场（气相和固相）两个不同物质构成的场的交互作用。它是一种复杂的、非线性、自组织的二相流。沙粒在空中的移动因其大小不同，受力不同等因素而产生不同的运动轨迹。沙粒在风中的运动形式大致分为三种：跃移、悬移和蠕移。而其中以跃移运动最为重要，跃移颗粒占全部输沙量的75%以上。所以颗粒的跃移运动机制一直是风沙两相流研究的重要内容。 1.1 气相控制方程 本文中气体的运动规律采用经典的黏性不可压缩N-S方程[]： 气相连续方程为： 气相动量方程为： 上式中，为压力；是气相所受的体力；是气体的动力粘。 1.2 沙相控制方程 建立和求解运动轨迹方程的关键确定沙粒在气流中运动的受力情况。研究发现，一个沙粒在气流场中所受到的力主要有[]：气流曳力，重力，气流升力，Magnus力，电力，Saffman力和Basset力。为了简化计算，做如下假设：(1) 忽略气流垂向运动和紊动作用；(2) 沙粒形状为球型体；(3) 暂不考虑电力对风沙流场的影响。根据前面的假设和跃移沙粒的受力分析按照运动学理论，沙粒在气流中运动的轨迹方程可以表示为： () 2 风沙流的SPH描述和离散 流体动力学问题的求解一般都是通过偏微分方程组，但很难求得解析解。因此就要找到一种数值方法来对问题域进行离散化并获取目标点上的变量函数及其导数的近似值，将偏微分方程组转化成一系列只与时间相关的常微分方程。SPH方法就是随之诞生的其中一种数值方法，它的基本思想是：（1）将所要求解的问题域分解成一群任意分布的粒子，这些粒子之间不需要用网格进行连接；（2）用积分表示函数的核近似法对任意函数进行逐步积分；（3）在计算的每一个时间步里，利用支持域里的相邻粒子对应的值叠加求和来取代场函数及其导数的积分表示形式2.1 SPH的基本方程 在SPH方法中，函数的积分表示式为[]： (4) 这里，是包含的积分体积，是狄拉克函数，性质如下：。如用光滑函数来取代狄拉克函数，则函数的积分表示式可写为： () 式中，是光滑函数，是定义光滑函数影响范围的光滑半径。对式()可通过相邻粒子的值来叠加求和进行近似，这样可以得到离散点处函数的粒子近似式。用粒子的体积来取代前面积分式中邻近粒子处的无穷小体元，那么粒子的质量为。最终在粒子处的函数的粒子近似式可表示成： () 式中的是支持域内粒子的总数。同理可得，粒子处的函数空间导数的粒子近似式可写为： (7) 其中，是粒子间距离，。 2.2风沙两相流的SPH耦合问题 本文基于SPH方法给出气固双向耦合作用的相流算法。该算法的基本思想是将固体相离散为与液体相一样的SPH粒子，其物理信息通过模拟气体“颗粒”与固体“颗粒”间的相互耦合等作用，来描述两相流动中的微观特性。 气相控制方程基本理论可以得到SPH中的气流场动量方程： (8) 同理，固体SPH颗粒的运动要满足力学中的基本关系