CFD中文网

@李东岳 这个案例我用的商软，发现空隙率0.01和0.5差别不大，影响最大的是阻力系数，我试了这两个空隙率，水都能很顺畅的流出来

想开个贴讨论下多相流界面追踪类方法目前的前沿水平，例如能算到的最高参数，各类方法都有什么限制等等。目前多相流算法里，计算界面追踪的多相动力学，比如droplet collision, droplet splashing, droplet impact on solid wall，这几个典型案例，目前有用VOF，level-set，phase field，LBM, SPH, 欧拉拉格朗日耦合其它模型的。这几个案例在不同Re和We下会产生许多种演化形态，实验都做出来了，但感觉计算还没跟上。

调研过不下百篇的文献，基本上关键词相关的都翻了一遍，在文献中看到droplet collision大概最好的也就是Re=6500,We=1500，密度比DR=1000左右的，大致是这个级别，像droplet impact有篇08年的用TransAT算到了RE=11000，而最近16年左右的用levelset这些也算到Re=8000，We=1500左右。大多数都是采用了一些改进的数值技巧来稳定模拟，而且网格也都比较密，计算资源消耗也都比较大，估计都得上超算，当然和实验也都对的上，基本都是到了prompt splattering这样的形态。

想讨论下大家都做过或者见过的这类模拟最高参数是多少? 或者往上算各类方法的限制在哪？ 因为多相流在算到高RE和We的时候会数值发散，稳定性应该是主要限制，还有例如VOF和level set会在薄液层处由于数值格式产生人工破碎的非物理情况，需要单独使用人工识别薄液层再稳定的方式，对不规则碰撞并不实用。

解决办法一：关闭3d加速，但会卡成ppt，其它办法暂未发现

用于求解不可压的pseudo-compressible method, 被认为有更好的收敛特性，见A. J. Chorin, JCP, 1967, pp. 12-26

@李东岳 LaTeX党遇到这个真的心累。只好把它留到最后一口气吃掉。

@东岳 我最近用reactingMultiphaseEulerFoam做算例，也是设置的时候某一相为0，湍流模型一直出现divide by zero，找了半天没发现错误，更改之后搞定了，谢谢东岳老师

最后自己修改了mapFields把映射内部场的代码删掉了

@肖艺 暂时还没研究这个问题，不过想了一下，网格运动的控制方程只能是一个，考虑平面4个点，运动应该是参考距离/角度进行同幅度运动？好像是这样

@东岳 谢谢教授

好像是编译的方法才能正常使用

您好，询问一下在ALM的使用中你有没有考虑分布因子的设置，可以交流一下吗，感谢

哦，明白了。

新人发帖，还不懂如何编辑帖子，请各位老师见谅。
捕获.JPG
该公式是安德森的空气动力学基础一书中对等熵条件下准一维流动推导出的速度面积关系式。根据该关系式，当管道流通面积A不变时，速度u是不变的，则根据质量守恒关系可以知道密度和压力都是不变的。但是我在观察真实条件下的抽气管道流动中，在管道流通面积不变的情况下，上游管道静压大于下游管道静压（下游安装有抽气机组）。请问这是否说明等熵条件下准一维流动速度面积关系式不适用于真实条件下的均匀管道流动？
流体力学小白，还希望有老师能指点一下。

$\frac{\mathrm{d} A}{A}=\left ( M^{2}-1 \right )\frac{\mathrm{d} u}{u}$

@TINGHAIK 谢谢您，问题已解决

2021-09-17 15-46-22屏幕截图.png

@李东岳 在 OpenFOAM里面的dnsFoam为什么这么菜？ 中说：

@ZifeiYin

Hi,

请参考下述内容，把第一行第二句话自行删去：

鉴于CFD圈子的小众性，“CFD中国”执行最严厉的友善管理制度：所有用户请以专业的态度对待所有问题。不得在所有版块发表轻蔑性言论；
例如，下面的语句被认为是具有轻蔑性的：“楼主回去看看书吧”，“这种问题你也问？”，“不会用百度？”
“CFD中国”建议用户对每个问题严谨友善处理。例如上述内容可以这样表述：“建议楼主参考《数值传热学》第48页，有对xxx问题的详细表述”，“个人认为这个问题是一个非常普适性的问题，在大部分的教材中都有涉及，楼主可以参考xxx”，“这个问题和CFD的关系不大，建议楼主去百度。google去搜索相关内容”

哈哈哈哈，东岳老师干得漂亮

@东岳 李老师，我没推出来

标准的限制方法是保证插值后网格顶点（vertice）的值不会出现极大值或极小值（比周围网格的最大值大或比周围网格的最小值小）因为梯度插值最容易出现越界值的就是在顶点处，从而用一个标量来scale梯度gradient；另外一种方法是保证网格的面心midponts of cell faces处不出现极大值或极小值即可，就认为满足了插值后不存在越界的问题。我想的是关于第二种方法，面心处插值后不出现最大值或最小值不能代表顶点处不出现最值，对吗

@李东岳 在 请问UprimeMean是如何计算的 中说：

所以你要计算得是$\overline{{\bfU '}^2}=\overline{\bfU^2}-{\overline{\bfU}}^2$？$\bfU '$是脉动速度，$\bfU$是速度，$\overline{\bfU}$是平均速度

老师，请问openFoam中可压缩求解器（rhoreactingFoam）使用Les模型，使用uprimemean计算的是favre平均的值吗还是雷诺平均的值，如果想要分别计算favre平均和雷诺平均可以做到吗

@肖艺 STAR CCM+
例子： https://www.cfd-schuck.de/en/cfd/anwendungsbeispiele/simulation-of-gear.html