CFD中文网

@fan820 使用DEFINE_DPM_INJECTION_INIT宏可以实现，把TP_VEL(tp)[i]定义成你当前时间步改变的气体流速就行

@cccrrryyy
返回到前面下，Hybrid RANS/LES方法现在还是挺流行的，据我了解，上文提的DES,VLES,以及流体机械文献中常用的PANS都属于混合模型。类似“LES”，应该是指从现象上看获得的结果跟LES类似，但是求解的方程跟LES不同，比如VLES模型（特指南航韩省思教授文章里发展的），就是基于RANS模型框架发展的，当网格较密的时候，便能达到LES的计算精度，但是方程中并没有涉及到过滤速度。

对于URANS，从现象上的比较，对于某些本身具有较强瞬态特性的流动现象，URANS可能会优异些。对于稳态计算（指最终结果与初始值无关），RANS跟URANS应该相差无几，对于非稳态（最终结果与初值有关），那URANS应该由于RANS的。

额外跑个题，对于URANS，多了个瞬态项，那每个时间步计算的时均速度（时均动量）都是从第一个时间步开始累加后平均后获得么？那初值的影响应该有可能很大有可能就没有。对于RANS，初值对结果几乎没啥影响，本身可能就跟一般牛顿法迭代就方程一样，最终获得近似值。

@李东岳 是的，李老师，我在测试不可压缩流中也同样遇到这种情况，但不可压缩流中热力学压力基本不变，只是动态压力变化，密度即使从理想气体方程中获得，感觉整个流场中的压力场（近似热力压力场）跟密度场被强制为常数，可压缩流场中，密度由连续性方程计算，我试过对入口边界的密度边界条件使用NSCBC，但是奇点只是从角点转移到附近的壁面处了，很奇怪，我最近也是一直在搜索相关文献或者论坛，看是否有讨论相关问题的。

@Vortex学习 我没用过Fortran程序计算Rortex，不过我的方法可以借鉴下，就是使用Tecplot打开Fluent结果，并导出基于节点的data文件，里面只包含xyz坐标以及UVW速度，然后我是在matlab中进行处理获得的R值以及矢量值，在替换data文件里的数值，最后再使用Tecplot打开data文件。这样可以处理较大网格量的数据，不过这个网格量得取决于Tecplot软件，我之前有个2000万的网格，Tecplot无法打开，我就很无奈了。

请教大家一个问题，使用可压缩求解器（密度变量是根据连续性方程计算）对泊肃叶管道流动进行模拟时，在入口的角点处（即速度入口边界与壁面边界的交界处）密度值以及压力值提高的较快，并随着时间推荐计算会崩溃，这种情况该如何解决啊？此外，我也是刚接触这部分内容，对于这类情况是否得使用一些特殊的边界条件（比如NSCBC）呢？

@qhwj1990 这几个数值就是单位的差别把，kj=1000j, kgmol=1000mol

发生相变的话，主要的量就是潜热，跟物质本身的标准状态焓无关，但是跟相变的焓差有关，也就是潜热了，所以液相为0，气相设置的大一点，但是焓差也不是代表潜热值，也就是后面红线的那段话了，其实就是单位的问题，这个焓值除摩尔质量，不就是潜热的单位么，以上就是我的理解了。

要不就在ICEM每个面都重新命名下，在重算吧