ubuntu cent都用。不支持也没关系。我这也没见得需要他们支持什么..

@李东岳 谢谢李老师的回答。

说明前期没收敛啊，网格无关性是要做的

这个文章没有讨论【一维压力波求解】

@bestucan 谢谢大佬

不错不错，文章挺狠，大佬的文章就是到位

@bestucan 万分感谢！感谢赵老师，感谢您！:xinxin: 已经向您邮箱发了邮件zhou_hli@163.com

@chszkc 谢谢大佬！！学到啦！！！

曲线救国：建一个矢量场函数grad($$Velocity[0])，再去用这玩意儿的第2个分量..

@李东岳 里面的（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模型参数的影响更大。如果谁有类似的文章，可以发来看看。