@李东岳 里面的（u）和（y）我想选用特征速度和特征长度，就是物体的运动速度和物体的特征长度，然后du/dy想在流场中提取

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

流体的动能通过涡耗散，而空间网格最小尺寸决定了你数值计算能捕捉到的最小涡的尺寸，因此小于你网格能够捕捉的尺寸的涡，其含有的动能无法通过数值模拟耗散掉。为了解决这个问题，可以理解为增加了一个和“空间尺度”有关的人工耗散，进而对小尺度的涡进行耗散，因此叫做大涡模拟。而刻画这个小涡空间尺度的标准和你的第一层网格有关，因为其决定了小涡的最大尺寸(也就是在空间上被过滤的小涡的最大尺寸)。

来个人回复一下嘛

比如
当用浮点数乘除得到一个整数的时候，很可能是 0.99999999999 1.0000000001。
这个完全无法预测。写函数解决？可是又不知道它本该是整数还是很接近这个整数。

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

@liujuncfd 这恐怕不行吧，比冲是推力除以单位时间消耗燃料的质量（或重量）。

@李东岳

李老师您好。

目前出口是采用如下方程
3.png

这里Uc尝试过：

使用上下层速度的平均值1.5。 在X方向上采集最靠近边界的三个点的速度值进行平均，作为Uc。

如果采用2的话，Uc是一个随着Y坐标以及时间变换的函数Uc(y,t)

Phi是压力泊松方程构造的中间量，这个查了好多文章都说的很模糊，我也就直接用了零梯度。

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

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

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

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

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

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

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

可以设计一个量纲的类，里面包含一个整形数组来记录量纲的分量，如[L, M, T]还要定义量纲的运算。然后设计一个标量的类，里面除了浮点数据之外，再加一个量纲，定义标量的动作时也要注意其中的量纲计算。同理还可以定义向量和张量。比如距离的量纲是[1, 0, 0]，速度的量纲是[1, 0, -1]。

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

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

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

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

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

实际工程上说的压力是什么压，压力表测到的是什么压，还有就是数值模拟计算得到的压力结果我们应该关注什么压呢？或者说从动压静压总压分布图上我们能看出什么东西。谢谢

这个放在他们搞涡那面，估计能讨论一上午 :duang:

我个人也不太熟悉。在我的文章里面，也不太关心涡的情况。

请问燃料燃烧时这种‘理论计算产生的平均总温’应该怎么计算啊？翻了基本燃烧学的书都没这方面的定义。温升燃烧效率.png

@jinlinna 感谢！

@程迪 666啊老铁，听君一席话，胜读百篇论文:146: :146: :146: :146: :146:

@一朵苍凉 我也不太懂，可以看看控制领域的相频特性，也许有帮助:zoule: