CFD中文网

fluent初始化显示mass diffusivity：invalid diffusivity please check material。用的组分输运涡耗散模型

@东岳 光看流场的矢量图可能确实看着还凑合对称，但是假如我要在出口得到RTD曲线（停留时间分布曲线），就会得到对称出口明显不同的RTD曲线？

或许CATIA的多截面曲面可以实现这个的建模。

描述热扩散的无量纲物理量热扩散因子究竟是什么？有无针对不同情形下的具体公式？
1.png 2.png 3.png

后来这个怎么样了

感谢分享！！！:huahua:

@李东岳 确实是我驱动段给的压力和温度不够，我错把激波后的压力和温度当成驱动段的压力和温度了，现在重新设置驱动段的压力和温度后，确实能够达到高马赫数了

这个其实是个好话题，很多大佬在sci里面也明确说了，LES之类的模拟，如果对于这种类似的圆管，施加cyclic边界条件，会对湍流产生抑制。但是我也没见过大佬们进行对比。

另一方面，感觉边界条件，直觉感觉，单单的一个patch引起的差异，不至于那么大。LES模型参数的影响更大。如果谁有类似的文章，可以发来看看。

你可以读取openFOAM处理好的网格文件啊，它的网格文件还是很有规律的，至少我现在已经实现读取和处理openFOAM的网格数据了，自己在编写求解器。

请问各位大佬，对于可压缩流大涡模拟的总湍动能如何定义呢？

对于不可压缩，大涡模型采用的亚格子模型采用的雷诺时均的形式，因此可以用雷诺时均计算大涡湍动能和亚格子湍动能，两者的和就是总湍动能。
$K=0.5*(\overline{uu}-\bar{u}\bar{u}+\overline{vv}-\bar{v}\bar{v}+\overline{ww}-\bar{w}\bar{w})+\bar{k}_{sgs}$ 但是对于可压缩流的大涡模拟，亚格子模型采用的为法夫尔时均形式，此时如何定义总湍动能呢？我的理解是
$\widetilde{K}=0.5*(\overline{\rho uu}-\bar{\rho}\widetilde{u}\widetilde{u}+\overline{\rho vv}-\bar{\rho}\widetilde{v}\widetilde{v}+\overline{\rho ww}-\bar{\rho}\widetilde{w}\widetilde{w})+\widetilde{k}_{sgs}$
请问各位大佬，我的理解对么

想尝试与一篇文献中实验结果作对比，但根据文中所提供的条件算出的入口速度很小（0.18m/s），以此作为入口的边界条件；同时入口边界的alpha.water也固定为1，计算一段时间后如图所示，发现入口处几乎没有水相的存在了。这是为什么？

ae0a7d49-651d-4636-a27b-c8a0ff4b9334-image.png

入口处的边界条件如下：

U { type fixedValue; value uniform (0.18 0 0); } p_rgh { type zeroGradient; } alpha.water { type fixedValue; value uniform 1; } alpha.air { type zeroGradient; }

也许再过一个世纪，湍流就不需要上帝来解答了，哈哈

@李东岳 李老师有没有兴趣看下风工程这块的有个商业软件WT，里面是速度-压力耦合求解器，真的非常快非常快。我看了下，软件底层是从PHOENICS来的，然后搞了一套叫MIGAL的算法。

现在CFD用的SIMPLE这种segregated solver确实太慢了，城市这种尺度下的流动本身大概率又比较简单，不可压，也没有化学反应之类的，可以暂时不考虑温度？可能会有颗粒流的问题，但是不是也可以先只考虑流动对颗粒的影响。这种情况下感觉流场本身的求解其实就看求解器的效率了。

我搜了下，OpenFOAM这边其实已经有一些coupled solver 的工作，不知道成不成熟。

JCP
让大家膜拜一下。:sunglasses:

@先先贝 我的方法都在之前写的帖子里面分享出来了。算几个时间步很快的。而且串行计算三个时间步就差不多了，你设置的自适应加密等级也基本到位了。这个时候再运行一次setFields，给加密后的较密的网格赋值，再去分块。
分块的话也得想想加密的部分会不会在你流场里面移动，你的不同核之间的负载要基本平衡，可能还需要你手动去划分区域。

@小八 ICER，文章已经找不到了。化工方向的:+1:

@xpqiu

十分感谢老师的指点！

我检查了所有异常算例，发现问题的确都与0文件中的初始设置有关。

在这个算例里，命令subsetMesh -overwrite frameCellSet -patch frame会将frameCellSet网格区域与frame链接，但是由于没有显式划分frame，因此，subsetMesh会默认生成一个empty类型的边界。

六楼的分析存在大问题，没有参考的必要了。

再次感谢两位老师的指点:146: ~

@东岳 谢谢老师，我已经找到网格生成规律了，从左至右，由下至上地规律依次生词。满足误之前需求工作，reNumberMesh有需要我再去琢磨吧，再谢！

很遗憾，木有啊:sad:

这个是不是用到了基于物理量+动网格技术。。。