【求助】模拟文氏管汽蚀的一些问题

enigmafwj

恳请各位前辈指教~！先谢谢了！

基本情况是这样的，我想计算文氏管汽蚀引起的压力脉动：
solver选择了基于压力的分离求解器；
几何构型为二维轴对称，采用结构网格，近壁面加密；
采用混合模型；
假设为均匀多相流；
综合考虑几何构型和流场压力梯度变化，采用realizable k-e和非平衡壁面函数；
流体为液氧，物性值由《低温工程》得到；
主相为液氧，次相为气氧，采用FLUENT自带的cavitation模型，fg=1e-9；
操作压力设为0；
采用速度入口\压力出口，进口湍流参数的计算以喉部直径为特征尺寸，假设入口全为液氧，且出口若有回流，回流也全为液氧，出入口温度在80K左右，入口压力大概有100个大气压，实验中出口无压力测点，流体出文氏管后流经一个阀和集合器后，压力大概为80个大气压，但由于电脑硬件原因，只模拟到文氏管，此处出口压力即设为大概80个大气压（这里不知道是不是造成倒流的原因，因为实际中文氏管出口压力并非80个大气压）；
压力速度耦合采用couple方案，松弛因子调整为：pressure=0.2，density=0.3，其它全为0.1；
差分格式：pressure=PRESTO!，vapor=QUICK，其它采用二阶迎风；
初始化：compute from inlet，Gau-pressure设为实验测得的入口静压减去0.5*density*v^2的值，v由入口质量流率得到。

先用常物性算稳态流场时，湍流粘性比一开始无超限，随着迭代，超限发生在出口边界处，迭代一些步数后超限消失；常物性稳态计算过程中，始终有倒流(reversed flow)在边界处发生，2-5faces；
壁面网格y+用自适应经多次调整为30-300，网格倾斜度在GAMBIT中已校正；
而后各物性参数(包括饱和蒸汽压和表面张力系数)写UDF用温度拟合，并给能量方程添加汽化潜热源项(UDF源程序见晚辈另一个帖子)，转入非稳态计算，时间精度为二阶，时间步长由公式算得。
单机3核运算（不知对UDF的写法有什么特别要求？）

问题是：
1。若把密度也用温度拟合，迭代没几步连续方程就发散了，请问该如何解决？
2。在解释UDF时，对那个Stack Size堆栈数有何要求和限制？
3。物性参数设置时，Material面板中的标准状态焓对计算是否有影响？好像只有焓差才有物理意义吧？
4。若采用FLUENT自带的cavitation模型，是否需要为各相质量方程添加液化\汽化输运源项？若要添加，是在phase Interation里添加还是在Boundary condition里的fluid分别选择各相来添加？
5。一开始的初始化时，点击Init后，是否还需要点击apply？点与不点有什么区别呢？
6。湍流粘性比超限的问题除了在稳态计算时有上述发生过程外，在各物性用温度拟合，并添加了潜热源项，转入瞬态计算后，在汽蚀区尾部、轴线附近有200-300cells超限且无消失趋势，这个超限在刚转入瞬态计算时
   是没有的，是随着迭代逐渐增多，最后稳定在300+cells，此时结果还可信吗？
7。瞬态计算过程中，始终有回流出现，且回流面积有赠多趋势，请问该如何解决？
8。出口质量流率在稳态计算时符合实际情况，为负，且与入口质量流率相差在1-3%；但是如上述转入瞬态计算后，随着迭代，出口质量流率最后变为正数，即出口全为倒流，此时结果不用想也知道不可信了，那么该
   如何解决呢？
9。瞬态计算迭代到一定时间步，即上述倒流、超限比较严重时，FLUENT无法继续迭代，信息显示如下：
   received suspend command for a pmi context that doesn't exist:unmatched id=2
      unable to read the cmd header on the pmi context,Undefined dynamic error code.
           The parallel FLUENT process could not be started.
      job aborted:
      rank:node:exit code[:error message]
      0:***:-1073741819:process 0 exited without calling finalize
      1:***:-1073741819:process 1 exited without calling finalize
      2:***:123
      注：***为计算机名

恳请各位前辈高人为我解惑，感激不尽！