CFD中文网

@Yung 这个 CAD 可以导出的格式有点少呀，0b3894d9-c7b1-4e5a-b38a-a1907754d22a-image.png ，输出-其他格式里面我这个版本只有这几个，你可以试一下 *.stl,*.igs不过我都没用过，不确定到底可不可以。

@jeffx 我也是报这个错误

利用fluent和rocky双向耦合时，fluent报错提示微信图片_20240703191405.png 微信图片_20240703191509.jpg

好的师兄我再去试试，之前用的是simple算法，然后今天上午试了一下simplec残差曲线又没出现振荡了，我把残差监视调一调再试试。

这些辐射系数都是和 wavenumber 相关, 在RT E中需要这些带量纲参数.
Thanks

@allanzhong 看这个意思，应该和你说的差不多，理论上肯定为，但实际上从理论转换为数值算法的时候总会有些偏差，算法在具体实现的时候又会有一些偏差。如果是这样，这些偏差我认为是可以接受的。

@两月三年 是的，这样显示的时候是表示选上了。注意选的位置就好。

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

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

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

还没有专门看过这方面的研究

@hxy970122 扫了一眼Fluent的组分输运方程，对于各个组分就没有重力的影响，所以如果组分分离就肯定是因为扩散率不同。但是在binary mixture里面，似乎两个组分的扩散系数是相同的？

重力的影响是在质量方程和动量方程里面，所以就是所有组分受到的重力作用是平均重力……

2e37529c-113a-4d19-a63f-e23bf85c3593-image.png

@evensun 求解器只负责求解数学式（控制方程），为什么要求解这个控制方程而不是别的控制方程是用户的选择，或者说取决于用户需要的是什么。
如果你需要百分百精确的描述一个非定常的物理问题，那你用稳态求解器原则上来讲是不合理的，但这个不合理不是说解方程这一步不对，而是你没有求解与你的需求相匹配的方程。如果你只需要部分描述这个物理问题，比如你关注时间平均量，那这时候用RANS去解这个问题，这一套下来是合理的。因为你需要的只是时均量，只是这个物理问题的一部分。
从这个角度去思考，你用稳态求解器计算了一个非定常的过程，求解器没有发散，甚至可能还体现出了一些有规律的震荡。这个现象本身其实是有意义的，但你在分析它的时候要非常小心。比如你认为“算得的结果呈现出了涡的发展变化过程”，这个表述是值得商榷的。

再啰嗦一点，关于RANS我一直觉得有个很有意思的角度。我们都知道雷诺数是决定流动的重要因素，我们可以认为控制方程也有一个“数值雷诺数”，它和普通的雷诺数定义类似，只是粘度项应该替换为分子粘度+数值粘度。在RANS中，数值粘度一般高出分子粘度2个量级甚至3个量级，所以尽管物理问题可能雷诺数很大，但数值雷诺数可能只在10的3次方或者4次方量级。10的3次方已经和我们通常认为的临界雷诺数接近了，导致RANS结果看起来是稳定的，就好像层流一样。
说白了，RANS模型对湍流的耗散非常强，把很多非定常的结构都给耗散掉了，留下了定常的部分。但这个耗散能力是有限度的，如果你求解的这个物理问题经过RANS模型后它的数值雷诺数正好处于临界雷诺数附近，那可能就是会有一些小的震荡。
这个湍流模型是自己开发的吗？你可以通过尝试不同的湍流粘度来验证一下，比如量级直接调大10倍，你说的这些震荡还存在吗？

@东岳 “......在计算域内部添加一个吹气孔进行吹气”
现在是要实现类似的吹气过程，将网格作为温度、质量（化学组分）源项，且网格赋予固定的速度以及压力，求解结果中反查这个网格的速度值和设定值是有偏差的（怀疑和组分浓度扩散相关，扩散系数存在导致网格有一定的附加“扩散速度”），结果的压力值误差更大（大于100%，猜测openfoam中应该用特殊的方式去限定这个压力约束）。

......BTW真正关心的问题来了，在fireFoam application源程序中没有找到质量运移的扩散系数D，也没找到施密特数Sc的设定位置（施密特数（Schmidt number, Sc）是一个无量纲的标量，定义为动黏滞系数和扩散系数的比值，用来描述同时有动量扩散及质量扩散的流体。施密特数的命名是为了纪念德国工程师 Ernst Heinrich Wilhelm Schmidt (1892-1975)。）

请问下会了吗

@李东岳 我发了帖子之后才发现你之前发过通知了，不是催更，抱歉抱歉~:chitang:

感谢感谢！

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

对于不可压缩，大涡模型采用的亚格子模型采用的雷诺时均的形式，因此可以用雷诺时均计算大涡湍动能和亚格子湍动能，两者的和就是总湍动能。
$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}$
请问各位大佬，我的理解对么

@小火人 为什么不直接用绿色版paraview呢

ICEM用block画的多个区域，怎么各自设置成单独的part?

来个人回复一下嘛