标记网格法的基本原理.docxVIP

标记网格法的基本原理

））2．3．2（（2告＋告＋宅＝－岳＋古（扫＋告）＋如gx皂＋笠＋宁＝考咕（＋启导

）

）

2．3．

2（

（2

式(2,I)、(2.2)、｛2.3)中u,u代表流体流动盯在水平(xi?]入竖直方向上y向）

的速度分歪，p为无呈纲化后的压力值，即P=P/(p 沪），其中P为压力值，单位为Pa.p为流体密度，单位为K护n3.Re为滋体雷诺指数，且Re=p--lJ·L/rJ=ll?l./v: Fr为弗劳德数，且Fr=UJ.{L-了，其中g为重力加速度．u、L分别为特征参考速度及特征参考长戌：9x、gy分别为x1年y向上的顶力加速度分昆．

将上述两个动虽方程(2,2)、(2,3)分别对．x,y求偏导．再将求导后的方程相加，

即可得压力泊松方程，如式(2.4)所示

vip=-

沪a2 a,uu 护斗 aJJ 1 沪D 护D

- 2笃 一歹 飞 七 如 ＋动 ＝Sp

(2.4)

式(2.4)中速度散昼D为

D=·竺＋竺

扣 ay·

(2.5)

根据连续性方程

,(2

可知，速度散盘D应等于O，即（2..4)式的源项Sp中不应

该再含有与D有关的项，但是研究人员Harlow和Welch在1965年发现，式(2.4}中与D有关的项是非常或要的，因为若给出的初始条件不准确，或压力泊松方程迭代解不完菩，在求解压力泊松方程时将导致迭代误差的逐渐积累，段终导致数值解的不棕定，但是可以通过限制时间步长，、空间步长等的大小来使数值解在迭代过程中趋千收敛，并大大减少迭代误差的积系，以使每一个时间步下获得稳定的数值解．

2.1.2 交错网格系统

由于采用有限差分法来求捎N-S方程组及连续性方程，因此首先要对待求解的区域进行网格划分，邓么这里将采用怎样的一种网格系统来标记各待求的未知物理：：：呢？接下来将详细介绍．

I.交错网格系统的结构

＾里 吵 只

,R

、

j+1

...

P心.,

P

＂二认二l! 卢

、~ J

心，，..

如口A 尸d 尸切

三产｀

己

犁

,flf

。心,,

Pv P心

llI

J一l.

心比

p克，

心沁

，一1

氏1:

tP f

t

P f 凸

小1 PaUI

| i I 计i

恤贮 水平方向 (X方向）

图2.2标记网格法中采用的网格结构

对控制方程组采用有限差分法求销，其网格系统结构如图,22所示，其中压力设宜在单个矩形网格单元的中心，速度设没在相应网格单元边线的中点上，接下来将详细讨论这种设置方法的好处．

标记网格法的基本特征是在离散化控制方程组时采用固定、纵横交铅的刚性网格系统，其是不随流体一起运动的矩形网格系统凹，这种交错的网格系轨的茎本思想为，假设在处理：：：维不可压缩粘性流体问题时，流体流动的方向为水平方向（设为X方向），垂直于流动方向的方向为Y方向，则可以对X.?Y横栽面的矩形区域进行网格划分，这样在水平、竖丑方向上均均匀分布者多个扣同的小矩形网格单元，然后在计江流休的速度、压力分鱼时可以取不同的网格节点，即速度分复u,v和压力p分别设窟在面一个网格单元的竖豆边、水平边中心及网格单元中心节点上（箭头方向水平向右、竖直向上的节点分别用来计算速度分量u,v)，如图2,2所示，这样做的目的是为了克服在丰交错的网格系统中偶节点和奇节点上压力值收敛千不同值的困难．

交错网格系统的优缺点

值彴注送的是，交错网格系统有一些其他方法（比如VOF法）无法比拟的优点比如：（I) 计饵的网格中心节点压力值是花定的，不会 出现不梯定的 悖况 (2)不需要丽先确定在刚性壁面上的压力条件：（3)计笲物理往u、V, p时，各自在网格单元的不同节点上进行求觯，每一个压力值不会同时收敛干两

个值(4)利用标记物粒子来虚拟地追踪流体的流动行为，其只随流体一起运动，并不参与纳维并i托克斯方程的求解过程，段后根据每一个粒子所在位堂的流体流动速度来计算其下一时刻的位登;(S)通过连接不同时刻所有标记物粒

子的位贡，可以动态观察流体流动时自由表面的位翌和形貌的变化；(6)由于该法使用了标记物粒子，通过求得每一个标记物粒子在不同时刻不同位置的溢度和应力（即标记物粒子的温度史和应力史，标记物粒子代表的就是流体本身），再结合聚合物结晶动力学方程，便可以求得最终注塑制品中不同位悝的结晶度分布．

当然，