运用离散相模型（DPM）模拟蒸发器中气液两相流

chenhuan1115

我现在做的是关于一个蒸发器内部流场和温度场的模拟，模拟制冷剂R22的气液两相流动，经过计算制冷剂液体的体积分数占2.5%，采用离散相模型来模拟制冷剂的两相流动，其中还伴有制冷剂液滴的蒸发，气态的制冷剂设为连续相，液态的设置成离散相，蒸发温度-10度，壁面温度设置成定壁温0度，边界调节是进出口escape  壁面为trap，调试过先计算连续相在加入离散相，但是在连续相收敛之后加入离散相残差曲线震荡的同样很大，振荡的周期是连续相和离散相耦合计算的步数相关的，在计算过程中还出现很多低温区域，会出现温度限制的提示：temperature limited to 1.000000e+000 in 202 cells on zone 2 in domain 1，调试了很长时间没有得到收敛的结果，所以来这里请教下版主，麻烦帮我看看，先谢谢啦！再问下版主在初始化的时候有没有必要点击initialize—reset DPM sources啊？

[ 本帖最后由 chenhuan1115 于 2011-4-12 09:16 编辑 ]

chenhuan1115

其中mumber of continuous phase iterations per DPM iteration调试过的值有5 ，10 ，20 ，50，震荡都比较大，得不到收敛的结果

chenhuan1115

选择的粒径分布式是rosin-rammler，最小直径1e-5，最大是0.0001，平均直径是5e-5，spread parameter是3.5 ，mumber of Diameters是10

chenhuan1115

随机跟踪模型

chenhuan1115

离散相物性设置

chenhuan1115

饱和蒸汽压是线性变化的

chenhuan1115

二阶迎风差分

chenhuan1115

离散相模拟的高手 版主xrs333帮我看看啊，谢谢啦

xrs333

物性的问题：
密度用常数不妥。
沸点比蒸发温度还低（当然，后者并非真的是物性参数）。
饱和蒸汽压很重要，你用分段线性，分段要足够多，数据一定要符合实际，最好用多项式，且要覆盖流场的温度范围。

chenhuan1115

原帖由 xrs333 于 2011-4-12 11:06 发表

                               
登录/注册后可看大图

物性的问题：
密度用常数不妥。
沸点比蒸发温度还低（当然，后者并非真的是物性参数）。
饱和蒸汽压很重要，你用分段线性，分段要足够多，数据一定要符合实际，最好用多项式，且要覆盖流场的温度范围。

谢谢版主xrs333，我先按版主说的调试调试
1.我发现离散相液滴的密度只有常数这一个选项，不知是哪里设置错啦
2.我设置的沸点是制冷剂r22的标准沸点，我的理解是沸点就是饱和蒸汽压下的温度，是不是沸点也不应设置成常数呢而是根据压力变化的呢？
3.饱和蒸汽压我分成10段，温度范围为-50度到40度，进口温度为-10度，定壁温为0度，这个范围应该可以覆盖的我现在用多项式拟合调试下
4.问下版主每次初始化完毕后用不用rest DPM sources？
再次谢谢版主，调试完毕后有问题再向你请教！

xrs333

1 我看错了，以为是气体密度。
2 那么蒸发温度要比沸点低，它是一个计算参数，高于该温度才计蒸发。沸点也应是压力的函数。
3 温度范围不但要覆盖实际温度范围，还要考虑迭代过程中温度波动。所以你给的温度范围肯定不够。
4 rest DPM sources是指把颗粒传递给连续相的动量、质量和能量全清零。初次初始化rest与否都一样。但是在计算中如果停下来修改了什么设置，然后你希望分散相计算重新来过，那么就rest一下。如果是为了获得稳定解，只要不发散，则是否rest也无所谓。

chenhuan1115

现在自己的疑惑明朗了很多，已经在计算啦，准备算下看看结果，谢谢你的指点啊

chenhuan1115

版主，经过这几天的调试，我调大啦饱和蒸汽压的范围，温度从163.15K到353.15k，设置了20个点，线性变化的，沸点也是线性变化的，蒸发温度仍然为263.15K，把连续相的速度进口改成了质量进口，这样可以输入表压354800pa，但是结果还是不太理想，不知道哪里还没调好，请版主指教下

xrs333

把蒸发温度降低。

mifengdatou

马克一下，正在学习dpm