CFD中文网

近壁面网格长宽比本来就很大，如果用RANS计算，有的模型对近壁面y+要求很小，比如SST k-omega要求y+≈1，但是对x+或者z+基本没有明确的限制，为了减少网格数量，一般流向和展向的网格无量纲参数可以在百千量级，这样的话近壁面网格长宽比就会很大，也是可以正常计算的。

Examinination of bubble collisions and coalescence in bubbly flows
D. Broder and M. Sommerfeld
还有这一篇，这好像都是小的那个做的，大的具体尺寸我也不知道，不过什么时候能在做实验不好说，实验仪器快被转移了，大的那个。。

@东岳 多谢提示，我再研究研究。

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

= =fvOptions(U)是用于处理源项的，暂时可以忽略。

SimpleFoam动量方程参考 http://dyfluid.com/icoFoam.html 中的方程1，只不过没有时间项。我建议你看看这个文档，simpleFoam和icoFoam区别非常小。

@FluidGao 嗯嗯，感谢。

在做单气泡尺寸的时候我都是用db = C_PHASE_DIAMETER(cell, thread_g)来获取气泡直径，最近想尝试双气泡直径，但是在参考的曳力计算式中有大气泡分数（fraction of large bubbles）这一项，不知道这项的值应该用什么程序句子来获取呢？

还请大佬能帮忙解惑，万分感谢！

前几天OKS的时候有个同学问欧拉方程是不是速度U和压力p按照声速传播。当时我给的答案是不是。在这里详细说一下这个问题。对于双曲系统
\begin{equation}
q_t+Aq_x=0
\end{equation}
其中$q$是求解的矢量，比如$[p,u]$；$_t$是对变量求时间导数，$_x$是对空间求导数，$A$是一个矩阵，但不是对角阵。双曲系统中的$A$是可以对角化的：
\begin{equation}
A=R \Lambda R^{-1}
\end{equation}
其中$R$是特征向量组建的矩阵，$\Lambda$是特征值组件的对角阵，如果定义$w=R^{-1}q$，最原始的方程就可以变成
\begin{equation}
w_t+\Lambda w_x=0
\end{equation}
这个方程式一系列解耦的双曲方程，可以所有的$w$都按照特征值的方向进行传播。因此对于欧拉方程，p和u并不是按照声速传播，按照声速传播的是$w$，$w$可通过$w=R^{-1}q$求出。

$w$也就是特征变量哦。如果感觉困惑，或许可以看看这个：http://homepages.see.leeds.ac.uk/~amt6xw/Distance Learning/CFD5020COMP/node17.html

adept

@李东岳
谢谢东岳老师！明白是什么意思了！

相变，多相气液固 、传质、传热 、
要是能把天然气水合物这玩意 模拟了，，那也差不多千古流芳了 哈哈哈

一个小白 小白 不能再白的人 在一群大佬环绕中 悄摸摸的说到

@李东岳 感谢东岳大神的翻牌

参考你的帮助，刚刚自我捅咕解决了。。。3q

@jeffx 我也是报这个错误

@ibelief 在 颗粒表面化学反应的文献 中说：

基于简单碰撞理论煤粉燃烧动力学模型的研究

非常感谢

@东岳 谢谢！

各位老师好：
我用Fluent计算流道进出口之间的总压降，发现几个问题如下：

我设置压力入口总压为1MPa，压力出口静压为0，计算结束后发现出口的总压反而比入口的还高，在没有外界做功的情况下我认为这是不合理的，总压平白无故增加了，为什么？ 我详细看了看出口的压力分布，采用Mass-Weighted Average表面积分方法，发现出口处的总压和动压差距很大，按理来说，出口静压设置为0以后，总压和静压应该是相等的，但结果显示并不相等（见图1）。换用Area-Weighted Average方法或Facet Average 方法，问题依然存在，进一步导出出口的ASCII文件，发现出口处各单元的总压和动压也不相等（见图2）。总压≠静压+动压这个问题仅存在于计算结果中的出口处，入口完全没问题。 改变边界条件，换用速度入口和outflow出口的组合时上述问题2得以解决。但只要使用的是压力出口，不管是速度入口压力出口的组合，还是压力入口压力出口的组合，都存在这个问题。

Fluent的具体设置如图3所示，因为模型比较简单，收敛的也很快，计算在87次迭代后收敛，收敛情况及监控物理量如图4所示。
我想知道为什么会出现上述问题，出口的总压为什么会比入口高？出口处计算出来的压力和我预想的不一致的原因是什么？恳请各位老师指教。

d2425a1d-5980-4190-952f-e3b2e246058f-image.png
图 1 采用Mass-Weighted Average表面积分方法时的出口压力结果

9a171ad1-dde7-4ed6-a22e-5f3958cdd438-image.png
图 2 ASCII文件中的出口压力结果

7163498f-eed6-450d-8255-6d377f2c2630-image.png
图3 Fluent的具体设置

a703bcd8-ae07-4b9f-a66f-4ca91b040653-image.png
图4 计算收敛情况及监控物理量的变化情况

还有为什么在相同边界条件下，NIXTURE模型加表面张力和不加表面张力，收敛性差距很大？有人说MIXTURE可以不用开启表面张力正确吗？请问各位大佬！

:wolaile: :wolaile: :wolaile: :wolaile: :wolaile:

:xiabanle: :xiabanle: :xiabanle: :xiabanle:

:tishizi:

各位安好，
想请教一下，对于离心叶轮来说，他不同叶高的截面不是一个圆柱，而是一个曲面。这种情况下，怎么提取不同叶高面。我们知道，在CFD-Post里面有turbo模板，可以生成对应的turboSurface，但是换了后处理软件，就没有办法生成了。
因此，想请教一下，不同叶高曲面表达是什么（图片中蓝色曲面）？应该是在笛卡尔坐标系和柱坐标的转换。查阅了CFD-Post的帮助文档和Fluent的帮助文档里面，都没有这一块内容。所以想发个帖咨询下大家
各位如果有什么线索也可以一起讨论，谢谢！
14da7003-9457-43c8-a206-c5bde261a26f-image.png
8672775e-2671-4021-b7ed-35fe6e11b1fa-image.png

[相关资料]
https://www.youtube.com/watch?v=C6xUFZD6wik - 创建柱坐标截面

但这样你就无法附加除了zeroGradient的边界条件了？