拟涡位移模型及其气粒两相流动数值模拟.pdfVIP

中国工程热物理学会 编号：986047 多相流学术会议 拟涡位移模型及气粒两相流动数值模拟} 葛满初 宋伟东 中科院工程热物理研究所 [摘要]本文首先用离散涡方法对绕圆柱和半无限长平板的流动作了数值模拟。在讨论粒 子和涡团耦合运动的基础上首次提出了粒子在流场中运动的拟涡位移模型，并对几个典型的 气粒两相流动作了数值模拟。文章最后对流场中粒子运动的瞬时状态作了统计处理，得到了 关于粒子在流场中分布的概率密度函数和粒子脉动速度的分布函数并用拟涡位移模型对壁面 附近粒子的运动一沉积和反弹作了数值模拟。计算结果是令人满意的。 关键词： 计算流体力学 气粒两相流动 拟涡位移 一、前言 含有大量固体颗粒或液滴的气体流动称为气粒流动，或多相流动。气粒两相流动的复杂 性给研究这一现象带来了很大的困难，要想给出一个能够精确描述这一过程的数学模型是相 当不容易的。在长期的研究实践中，人们发展了各种简化的数学模型，并在应用中取得了较 好的效果。目前，多相流系统的研究方法主要有三种：宏观连续介质理论的欧拉方法，微观 的运动理论和粒子运动’的拉格朗日方法，经过众多学者的不懈努力，气粒多相流研究领域已 经建立起一套较为完整而相当有效的理论和方法，许多模型也在实际工程应用发挥了相当大 的作用。但是对于这一复杂现象内部机理认识的不足仍然制约着气粒多相流研究的发展。本 文在研究中引入了离散涡方法并首次提出了模拟粒子在流场中运动的拟涡位移模型。对气粒 两相流动在涡团的尺度上进行了直接模拟，并试图对这一复杂现象内部机理作初步的探讨。 二、基本方程、模型和求解方法 离散涡方法 方程中的非线性项用一组尺度很小而又具有有限涡核的涡元(”涡团”)间的无粘性相互诱导作 用来模拟，(2)N-S方程中的粘性项可以通过在涡元移动中附加一个高斯分布的随机位移来模 拟，(3)固体壁面上无滑移边界条件可以通过在壁面上生成新涡来实现。涡方法对二维不可压 缩，非定常，高Reynolds数流动的直接数值模拟有不少优越性：例如，数值耗散小，具有自适 应性，无需区分层流或湍流流态，且适合于模拟绕复杂几何外形的分离流动． 二维黏性流动的涡量动力学方程有如下形式 旦竺+矿．V∞：oV2国 a (1) 涡方法在求解方程时采用了算子分裂的办法．方程被分成对流方程 丝+矿．V∞：0 a (2) ·国家基金和国家攀登B项目资助课题 皿一73 和扩散方程 丝：uv2∞ 曰． (3) 并顺序求解．在求解方程(2)中。连续的涡量场离散被为一系列涡团，即 甜(z，f)=∑rJy。[z—z，O)】 “ (4) 这里N是涡团的个数，y，是一个表示一个位于Z，处的涡团内涡量分布的函数它满足归一条 弘i)d孑=1 (5) 如果假定涡团内不同的涡量分布仅仅依赖于一个参数仃，，那么y一可有如下形式： y，(z—z，)=(1／盯?)／(1z—z，I／or。) ，m 其中，称为形函数或分布函数，它对于所有涡团都是一样的。盯·称为截断参数，它反映了 涡团的尺寸大小，y·称为光滑函数． 根据Biot—Savart定律，由～个给定涡量场可以求得速度场 v(z，t)=JK(x-Z。)∞(z’，t)dV。