CFD与DEM耦合的FLUENT实现.docVIP

CFD与DEM耦合的FLUENT实现

2007中国ANSYS-CFD年会论文集

CFD与DEM耦合的FLUENT实现

罗志国～邹宗树

东北大学材料与冶金学院～辽宁沈阳110004

摘要:本文利用商业软件FLUENT为平台，将DEM程序利用UDF编译进FLUENT，从而实现了DEM与CFD的耦合。本工作充分发挥了FLUENT的现有功能，同时通过二次开发对其功能进行了拓展。最后利用其对水平管的气力输送过程进行了模拟。

关键词:离散单元方法计算流体力学颗粒-流体流动

引言

单个颗粒的信息(速度、空间位置、颗粒之间的作用力)对深入了解输送过程是十分必要。目前最先进的技术PEPT(PositronEmissionPositionTracking)也只能获得单个颗粒的运动轨迹。利用连续介质力学方法和离散介质单元方法的耦合方法，即CFD-DEM(computationalfluiddynamicsanddiscreteelementmethod)可获得每个颗粒在任意时刻的所有信息，并可进行定量分析。[1-5]CFD-DEM是近几年发展起来的一种用以模拟气固两相流动的计算技术。该方法把气相当作连续介质处理，通过求解连续介质控制方程获得气相运动信息，即所说的计算流体力学方法;而固相则当作离散介质处理，通过求解牛顿第二定律控制方程获得单个颗粒的运动轨迹，即离散单元方法。两相之间的相互作用通过牛顿第三定律进行耦合，也就是两相之间的相互作用力应该大小相等，方向相反。

前人的CFD-DEM程序基本上是自编程序。本文利用商业软件FLUENT为平台，将DEM[6]程序利用UDF编译进FLUENT。本工作充分发挥了FLUENT的现有功能，同时通过二次开发对其功能进行了拓展。最后利用其对水平管的气力输送过程进行了模拟。

1数学控制模型

1.1颗粒相-离散单元法

颗粒一般有两种运动形式，平移和旋转，由牛顿第二运动定律确定。在运动过程中，颗粒之间、颗粒与壁面之间法发生碰撞以及受周围流体的曳力作用，动量和能量均发生改变。任意时刻t在两相流体系中的颗粒平移及旋转运动公式如下。

根据牛顿第二定律，颗粒运动方程可表示为:

平动方程:

kidvi，，m,f，mg，f，f,ipf,iic,ijd,ijdt,1j

(1)

转动方程:

kid,iI，，TT,，,,,icijrijdt,1j

(2)

其中m，I，k，v，,分别代表颗粒i的质量、转动惯量、接触颗粒数目、以及该颗粒iiiii

的平动和转动速度。f，mg，f，f分别代表流体对颗粒的曳力、自身重力、颗粒之间pf,iic,ijd,ij

的弹性作用力、阻尼力，T，T分别为切向力产生的力矩和转动摩擦力矩。c,ijr,ij

罗志国，1971年出生，男，满族，辽宁省沈阳市人，讲师，博士，研究方向为颗粒体系的模拟与仿真，冶金过程的模拟与仿真，铁水预处理等。

135

2007中国ANSYS-CFD年会论文集

1.2气相-连续介质力学方法

流场由基于流体网格的连续性方程和纳维-斯托克斯方程来确定。质量方程:

,,，，，,,,u,0t,动量方程:(3)

kc

fpfi,,,,，，,ufi,1，，p，，，,,,,uu,,,,，,,,,，,,gfftΔV,

(4)

式中P代表压力，τ，ε，?V指粘性应力张量、空隙率及流体网格体积，具体各项力的计算公式列下表。

表,各项力的计算公式

力及力矩符号公式

34**2f正向弹性力,ERn,cn,ijn3

1**2阻尼力fdn,ij,,，，8mER,vnijnn,ij

3，，2,,,,,,,min,f,,t,ijt,ij,maxscn,ij,,,1,1,δf切向弹性力ct,iijt,ij,,,,,δt,ij,maxt,ij,,,,，，

12,,,,1,t,ijt,ij,max,,f阻尼力dt,ij,,mF,6vtijscn,ijt,ij,,,t,ij,max,,

21,,，，0.5cρπRu,vu,vεf曳力f,id,ifiiiiii0

R，(F，F)T力矩i,jict,ijdt,ij

fMg重力g,ii

Em，m111R*2,,iji其中:,,,,,E,,，m,,,,,n,ijt,ij,maxsn2*22(1,)2(1,,),RRRRiji

,，，,，，,v,v,v，ω，R,ω，Rv,v,n,nv,v，n，nijjijjiin,ijijt,ijij

22，，1.5,logRe2ρRεu,v，，,,fiiii4.810p,i,,,,,,Re,,3.7,0.65exp,,c,0.63，p,i0,di0.5,,μ2f,,Re,pi，