讨论 Fluent计算效率优于OpenFOAM

[Foamer24]

@wwzhao 谢谢，库郎数中的速度指的是特征速度，就用了远场速度，不过用您说的最大速度来预估更合理。还看到东岳老师分享的链接，介绍三维非结构网格的库郎数计算，这个公式中的最小体积可以用画网格软件得到，但不知道通量该如何取值。

[wwzhao]

@Foamer24

一方面，非结构网格的库朗数是通过网格体积和通量（与速度有关）计算的，而一般我们计算库朗数的公式是最简单的一维情况，这两种情况本身就有区别，参考这个链接。

另一方面，有些网格的当地局部速度比特征速度大得多，所以这些网格的库朗数比用远场速度估算得到的库朗数要大。

[李东岳]

@Foamer24 在 讨论 Fluent计算效率优于OpenFOAM 中说：

Mean and max Courant Numbers = 0.001871585 0.5321452

你努力的方向不应该去找库郎数，你的算例只是发散了，库郎数只是表现，不是根本。建议找找发散的原因去去根

[cccrrryyy]

@Foamer24 感谢分享！我没用过fluxScheme，所以这里什么比较好不太清楚。时间步我一般刚开始都是1e-06，之后根据最大CFL数再去调。我个人的感觉是如果需要调到5e-07以下我就会比较注意了，这种量级感觉算的相当慢了。我的网格有内流和外流两部分，内流的话边界层在0.1mm上下，和你差不多。这个就用雷诺数估算了一下。外流（射流）的话就没care边界层了。

我觉得时间步的绝对大小参考意义比较小，毕竟也看你的硬件好坏。可能单个时间步需要的时间（我的大概在一分多钟）和总时间步数（我的大概50，000左右）比较有比较价值。另外，我的ddtScheme一般用backward，因为是二阶的。对流项我也一般不用迎风格式，迎风格式对LES好像不太好。线性求解器里面我的压力方程一直是用PCG在解，没用GAMG。这个点其实我很关心，就是PCG和GAMG到底孰优孰劣，希望有经验的来探讨下！

你这个同样的网格用Fluent只需要1e-06这种量级的时间步吗？这个感觉不科学啊，同样的网格怎么会差这么多呢。

[Foamer24]

@cccrrryyy fluent的网格边界层更密，为1微米，满足y+小于1。

[cccrrryyy]

@Foamer24 意思是Fluent的网格更密然后还能用更大的时间步？

[Foamer24]

@cccrrryyy 是的

[李东岳]

你需要深挖一下

给你举个例子，下面是运算的一个log文件，时间步长非常小，2e-5，算2天之后，结果完全可以是正确的，但我觉得有问题

在自适应调节时间步长的情况下，时间步长会自动跳转满足CFL标准。其实对于我这个，如果在代码里加上几行，输出Ur的话，会发现下图中全场计算域中，就这么一个地方Ur的值特别大，导致局部库郎数过大。

我建议你看看你的算例是不是也是这样，我怀疑跟我差不多

[Foamer24]

@东岳 在 讨论 Fluent计算效率优于OpenFOAM 中说：

Ur

您好！之前测试的时候，壁面附近网格越密，网格高度变低、长宽比也变大，要求的时间步长deltaT就越小。不知道您说的Ur是什么定义？不知道在壁面附近有没有特殊处理，增加它的时间步长。

[李东岳]

@Foamer24 我那个只是举个粒子，你可以看看你的U，不一定看Ur

[nicklj]

@东岳 在 讨论 Fluent计算效率优于OpenFOAM 中说：

你需要深挖一下

给你举个例子，下面是运算的一个log文件，时间步长非常小，2e-5，算2天之后，结果完全可以是正确的，但我觉得有问题

在自适应调节时间步长的情况下，时间步长会自动跳转满足CFL标准。其实对于我这个，如果在代码里加上几行，输出Ur的话，会发现下图中全场计算域中，就这么一个地方Ur的值特别大，导致局部库郎数过大。

我建议你看看你的算例是不是也是这样，我怀疑跟我差不多

想请教一下遇到这种情况应该怎么办呢？特别算vof的时候，界面某些点出经常出现很大的速度，导致时间步长小的无法接受。