物理模型和数值方法

  • 格式:doc
  • 大小:143.74 KB
  • 文档页数:4

物理模型和数值方法

1.1 基本物理模型

图1.1.1

坐标原点位于槽道中心,X-,Y-,Z-分别是流向,展向和垂向;在X和Y方向上,作用周期性边界条件,考虑壁面无滑移条件。计算域是22hhh,对应的网格数是129129128;用欧拉网格来离散流体区域和时间步长t,同时在X,Y,Z方向上最少的网格点能保证最小网格长度比流体尺度小。

1.2 物理方程

/0iiux

(1)2.21Reiiijijjjiuuuputxxx (2)

iu是无量纲速度矢量;ix是笛卡尔坐标;P是脉动运动压力;ij是无量纲压力梯度;Re是雷诺数。

34PDPpppduCuuuudtd (3)

其中,pu是颗粒的速度,u是流体的速度,pd是颗粒的直径。

0.6872410.15ReReDppC (4)

iidyvdt (5)

,11/NifpDjjfVmf (6)

现在开始考虑颗粒与流场相互耦合时,对N-S方程进行改造以适应变化后的流场;在我们的研究中,我们只是加入了单向的颗粒程序,并未对N-S进行改造,现在简单介绍一下对N-S方程进行改造的方程,方便以后更深的研究。

首先公式(1)是一个大家非常熟悉的流体连续方程;颗粒对流体的影响是在公式(2)中是作为一个额外的动量项出现的,在公式(2)中p是压力,f和分别代表流体密度和运动粘度,if是颗粒对每一个流体单元的作用力,额外的压力梯度使得流量速率保持常数。其中公式(1)和(2)依据一个半隐式分步法,时间集成对扩散项采用的是隐式二阶Crank-Nicolson方法,对对流项和外部力采用一个二阶Adams-Bashforth方法,在空间交错网格上采用中心差分格式;由于大规模的连续性,对于泊松方程中的拟压力项采用傅立叶变换。公式(3)是得到简化的颗粒方程,其中iv是颗粒速度,iu颗粒位置上无扰动的流体速度,pd和p分别代表颗粒的直径和颗粒密度,DC是修正在颗粒惯性影响下的拖曳力,系数公式如方程式(4)所示;方程(5)中iy是获取颗粒距壁面的距离,对其求导求得颗粒速度;if是一个控制体流体对于颗粒总的作用力力,pm是颗粒质量,N是一个控制体积V中的颗粒数目,Djf是方程(3)中的拖曳力。

带入方程(1),(2),(3)中得到粒子所在位置流体的速度,进而可以与自由流体进行比较,但就目前来说,我还不太了解对于此程序的实现,故只是对自由流体程序基础上加上了单行的颗粒程序进行加入颗粒后的流体模拟。在以后工作中对其进行实现。

DC是阻力系数,Rep是颗粒雷诺数,Re/pppduu;/ppfSt,反映的是粒子反应时间,其中22/;/18fpppud,(是流体动力粘度)。

图1.1.2

粒子反应时间对粒子轨迹的影响,小颗粒与流体同步,大的颗粒轨迹则会变化不同于流场;Hussain很好的验证了这一现象,气溶胶颗粒在重力场下,随着流场的分布,同时能很好的看出流场的流动。

根据以上方程式,用合适的插值方法离散,来获得颗粒的速度,以及颗粒所在位置流体的速度,从而与自由流体进行比较,分析颗粒对流体的影响。

1.3 数值方法

首先对解析方法进行一个简单的介绍;

(1)流体求解方法

PS(伪谱法)伪谱方法,是谱方法中常用于求解高度非线性方程的一类方法。它又被称作“离散变量表示法(Discrete Variable Representation method)”。它能很方便的处理正交网格,提高计算精度。

FD有限差分方法(Finite Difference Method):是计算机数值模拟最早采用的方法,至今仍被广泛运用.该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域.有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组.该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法.对于有限差分格式,从格式的精度来划分,有一阶格式、二阶格式和高阶格式.从差分的空间形式来考虑,可分为中心格式和逆风格式.考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等.目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式.差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定.构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法.其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度.通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式.

FV有限体积法(Finite Volume Method):又称为控制体积法.其基本思路是:将计算区域划分为一系列不重复的控制体积,并使每个网格点周围有一个控制体积;将待解的微分方程对每一个控制体积积分,便得出一组离散方程.其中的未知数是网格点上的因变量的数值.为了求出控制体积的积分,必须假定值在网格点之间的变化规律,即假设值的分段的分布的分布剖面.从积分区域的选取方法看来,有限体积法属于加权剩余法中的子区域法;从未知解的近似方法看来,

有限体积法属于采用局部近似的离散方法.简言之,子区域法属于有限体积发的基本方法.

(2)粒子时间集成方法

AB2:二阶Adams-Bashforth插值方法

CN:Crank-Nicolson插值方法

RK:Runge-Kutta插值方法

(3)流体速度插值

L6:六阶Lagrange插值方法

LH4:四阶Lagrange-Hermite多项式

H3:三阶Hermite插值多项式

TL:Tri-linear方法(控制法向网格点)

(4)垂向网格点的控制

HT:离散双曲正切拉伸

Chebyshev:切比雪夫距离

在我们研究中,用到的数值方法分别是:

① 流体求解方法:有限差分法-Crank-Nicolson method(有限差分方法的一种,它在时间方向上是隐世的二阶方法)

② 流体时间集成(非线性+粘性条件):AB2(二阶Adams-Bashforth插值方法+ CN:Crank-Nicolson插值方法)

③ 粒子时间集成:RK4(四阶Runge-Kutta插值方法)

④ 流体速度插值:H3(三阶Hermite插值多项式)

⑤ 网格解析度:128*129*128

⑥ 垂向离散点的控制:切比雪夫多项式

根据这些方法,以程序语言实现,从而对问题进行计算,得出计算结果。