|
楼主 |
发表于 2005-5-19 17:28:22
|
显示全部楼层
[求助]请问在Phoenics中可以计算浓度的问题吗?
[这个贴子最后由wuweiyang在 2005/05/19 05:33pm 第 1 次编辑]
谢谢大家对这个贴子的关心和支持,加入变量的问题已经解决了,但由于所计算问题中在部分区域中气体的流动速度很小,几近于零,所以就要考虑到扩散的影响。而扩散系数恰恰为张量,且与速度相关。方程如下:
∂C/∂t + Div•(Cv) − Div•(D•div(C)) − C’W = 0
其中:C为混合气体中组分的浓度,C’W为源项,这是我要计算的方程,
而且D的计算公式如下:D = Dt + Dm (一般来说Dm可忽略不计),其中Dt(i, j) = αl v δ(i, j) + (αl − αt)v(i)v(j)/v
αl和αt表示一些常数,v(i),v(j), v为速度分量和速度,δ(i, j)就是i = j时为1,不相等时为0的那个。
看过原来POLIS里面的例子,里面是不是有个是专门模拟电子原件散热的,有个特殊的物体可以设置横向和纵向不同的散热系数。就想这也可以直接设置吧,但是至今还没找到。
也想过用PLANT来设置,可能是我自己设置的问题吧,将方程中的扩散项先用TERM命令去掉,然后将其离散后扩散项作为源项输入,但是可能自己没有编过程,边界上总是处理不好,还有就是一计算就振荡,就再也没计算。
看了PHOENICS中对浓度的计算公式好像是这样的(∑对j求和):
ρ*∂C/∂t +ρ*∑(∂(C*Uj)/∂xj) =ρ*∑(∂(νL/PrL(C) * ∂C/∂xj)/ ∂xj)
与上面的公式相比,其中扩散中少两项,而且PrL(C)只能通过设置层流普朗特数的PRNDTL来设置,我是没想出来怎样去设置这个值。
下面的图片为上面两个公式的图片,
|
|