CFD中文网

试了一下关闭超线程和关闭E核（小核）的性能差别。关闭超线程对8核并行速度无影响，还是157s左右；超线程和小核都关闭的话，8核并行速度会变慢，计算时间200s以上了。

想问一下李老师这个计划还有吗！！很想要这本湍流模型:mianmo:

@bingningmeng45 你好，请问您做圆柱绕流划分网格的话，要满足什么样的条件才是DNS啊:143:

@linhan-ge 在 多相流VOF 中说：

VOF-DEM Study of Solid Distribution Characteristics in Slurry Taylor Flow-Based Multiphase Microreactors
强啊前辈，我这进展还缓慢着呢
刚看了您“VOF-DEM Study of Solid Distribution Characteristics in Slurry Taylor Flow-Based Multiphase Microreactors”的模型部分，请问您是在interFoam上加了Laplacian filter就可以抑制假速度、算出液膜了吗？因为我现在就做气液流动，速度也跟您差不多（1cm/s~10cm/s），所以想请教一下这个方案可行性。我直接用interFoam算的液膜厚度非常厚，不知道问题在哪里。
另外我看到您是设置接触角90°，那最开始设置气相区域的时候需要留出一个液膜厚度吗？还是直接把整个管径宽度都设置成气体都可以？我比较担心如果把整段管都设置成气体会让气体“粘”在壁面上。
谢谢！

@ZhangMengYi 没有，好像只见过Deen那面搞过，不知道他们能不能给你

@李东岳 老师，DPMFoam的解析目前您是否已经整理完成了呢，在哪里可以学习到呢

@coolhhh 感谢！

@队长别开枪 好的 谢谢前辈！

我是南京航空航天大学航空学院2022级硕士研究生吴祥清，跟随张老师链接文本从事CFD研究，本人的主要研究方向为高速出入水，初步学习商软入门CFD现在准备编程写算法。查找学习资料了解到CFD中文网以及老师，所以想申请2022CFD青年成长支持计划，我目前正在使用LS-DYNA软件进行垂直破冰仿真，开始着手发相关论文。在本科期间接触流体力学便想深入了解这个专业，也申请了这个方向的研究生，如果论文写作顺利，便会在研究生二年级申请硕博连读，为国家军事建设贡献自己的一份力量。目前在上课之余也开始着手C++入门。但是CFD学习之路道阻且长，希望可以得到老师资助来更好的学习CFD。

@turbulen 没有见过你说的AR线性滤波， RFG 算是一种谐波合成法吧， RFG DSRFG 以及G.T. Bitsuamlak "A divergence-free inflow turbulence generator using spectral representation method for large-eddy simulation of ABL flows"这个方法，都应该算是一类的。做LES还有用Synthetic Eddy Method 和Digital Filter Method的，这两种方法在openfoam.com中有，如果你在用openfoam.org,可以去找一下这个网址：https://github.com/NHERI-SimCenter/TurbulenceInflowTool/tree/master/openFOAM_code/OpenFOAM6/turbulentInflow

@Foamer24 +1

@李东岳 在 如果用DPMFoam求解稀相流会怎么样？误差大么 中说：

@Alvin 不知道你的问题解决没有。我在我自己写的求解器中，通过

将p方程换为p_rgh方程 将g去掉

解决了。因为轻颗粒的浮力在方程\eqref{gra}中包含了，重颗粒也如此。就像 @dzw05 说的，连续相可以认为是单向不可压缩流，类似simpleFoam。 @Alvin 你试试将DPMFoam中的phig去掉，运算看看？

这个1,2两项是二选一即可吗，phig指的是不是phicForces里面的rAUcf*（g&mesh.Sf()）

@cav 你好，请问关于气泡群直接数值模拟的气泡合并问题有进展了么？想了解一下~

@红豆沙 fluent meshing画的多面体of是支持的，多面体和六面体核心结合的新型网格形式of不支持，不知道新版的有没有更新

@random_ran 您好，想咨询一下，这么多的网格您是如何进行后处理的呢？paraview 会崩溃的吧？

刚接触sprayfoam，请教三个个问题：

初始场指定区域怎么布置液滴？setFields？ 液滴的喷射位置和喷射速度怎么设置为跟随inlet边界？ 我模拟的是预气化的液滴，比如C7H16既有液滴也有气体，我在算例里面该怎么区分他们？

谢谢前辈，希望解答或者提供一些解决问题的方向:140:

这个是我用openfoam10跑的。这个数据还是不太好。我看了一下原始文献以及论坛里面的一些文章。这个跟湍流模型很敏感。我尝试了一下调节进口的湍流强度，也很敏感。这个我后续不会再调了。我觉得应该留给专门做这方面的人去调试一下

impingingJet10.tar.xz

@falltime nice！一针见血。厉害厉害！:146: :146:

@李东岳 走一步看一步吧..先把前辈们的东西学会:xiezuoye:

@Ryo 前辈好！ 请问Gmesh能够对非结构网格划分block吗?

@李东岳 跟湍流项有关的参数，在发散的时候，在进口是正常的。这个case我也特意关注yplus+， 优化到<10左右。

附件：发散的各参数截图
https://drive.google.com/file/d/1U7dGKk-jv6OYMdOf8Ukiy1Gak648iJEU/view?usp=sharing

我用子系统上的OpenFOAM，总感觉并行计算没用，和我不开并行计算时间差不多。。。而且还decomposePar不了，好像是安装的时候少了啥东西。。。

@知行 5.p63、p96~p107 每页的转置符号$\rightarrow$ $\bf T$

@夜阑烟寒 从昨天开始domecross用不了了，你那是不是也是？我用了得一年半了吧，第一次开始这么长时间用不了

@Zhy2022 :high:

还有一个myFoam，github可以搜到，里面有个lusgsFoam求解器

@coolhhh 看起来这个应该不是数值上的问题，而是从物理到数值模型的改进的问题。这个问题如果能处理掉，我觉得可以发一篇sci，名字就叫“An improved Yang turblence model for large scale atmosphere simulation”

:chitang:

@Calf-Z-DNS 你好，看到你提的这个问题，请问你解决了吗？
我也想计算整个流动区域近壁面好几层网格的yPlus，而不是仅仅贴壁面的一层网格的yPlus，暂时我看到的wallShearStress和yPlus的计算都是贴壁面一层的值；
或者像东岳老师的示例代码那样直接给定了不变的uTau，而其实2维流动中，沿着流线方向uTau是变化的，所以我觉得比较复杂的是，怎么用某一壁面网格上的wallShearStress（或uTau）确定 对应法线方向上的多层网格yPlus，这个对应关系不知道如何寻找。
希望我表达清楚了，非常感谢！

@李东岳 好的李老师

595269271@qq.comqq.com

@东岳 https://www.linkedin.com/jobs/view/postdoc-in-cfd-of-turbulent-environmental-flows-at-university-of-houston-1628249421
我离青千条件远着呢，先在这边混着，等我老婆博士毕业了再说:chigua:

@Caijinjin 如果可以提供下述信息：CPU型号、节点数、操作系统、价格这些参数（如果不方便公开，可以发我邮箱），我可以远程给你向日葵看一下。另外，你们为什么换IB交换机？我看你们去年6月份买的集群，这半年一直没用么

@尚善若水 边界条件里面有inject vorton的时间我记得，应该没有很久吧。我觉得对计算速度应该没啥影响。

@zhe 直接设置水体区域，没有另外处理

@zhe Pout 是processer输出的 mpirundebug 这样的话 就并行debug了 你试试~就知道
还是info最直接

@李东岳 李老师，早上好。我翻了翻我以前做的文章，以弦长为Re=6000的矩形case，和 @gengchen 发的图有点像，网格尺寸是0.0125弦长，用的是笛卡尔网格，yPlus可以小于10。 如果这位同学的网格量估计在5~10亿以下，我应该是可以提供这个算力的。

gengchen

@星星星星晴 明白了.我这个方法确实直接套用不上.并且碰撞模型从同一个cell扩展到周围cells之间的碰撞,还会有其他的问题.我再想想,谢谢了.以后常联络!

我测试了一下，单核计算，对于这个150万网格的算例，每个网格的操作要占用30ms来计算，150万网格就是45000秒。

因为我对Liutex这个计算过程不熟悉，不清楚这个是否能做简化。

@星星星星晴 不好意思，一直在工作和学习两边忙，没有怎么来论坛。都还好，谢谢。至于openfoam的粒子追踪模型，都是不考虑粒子体积的，只是一个点。但在欧拉下计算NS方程，是考虑了粒子占网格的比值（这也不是考虑了体积）这只是考虑了由于粒子存在在网格，占据了原本流体的一部分体积，对计算流体造成了影响。总的来说，还是算的流，粒子只是一个源项。

@李东岳 多谢李老师

调换完全没问题。水湍流、空气层流也没问题。

@散漫守望2016 多谢您的指导，用自己的foam算过软件自带的例子，结果和标准的例子结果比较吻合，应该不是程序的问题。根据您的分析，我觉得很可能还是我模型结构的问题。十分感谢。

@aaaachenjian 模拟是在微重力环境下，上下壁面有个温差，在自由面上加个热毛细力，已经放弃用VOF方法了，像这种变形较小的直接挪网格，这是算法天生问题不太好解决

@李东岳 李老师，对于周期边界下的压力场，piso循环中有pEqn.setReference(pRefCell,pRefValue);用来确定压力场。那在计算周期边界下的被动标量时，应该用什么来对温度场进行修正呢？
:135:

@星星星星晴 在 LES与RANS的结果不相同 中说：

大佬求下面链接的文档啊。。
http://www.dyfluid.com/LES.html

http://www.dyfluid.com/theory.pdf 我最近更新在这里了 你可以看下 51页

我整个进口是开放的，不是东岳大佬说的有壁面的情况。

如果你不是这个边界条件，是开放的，那存在卷吸的情况，会很不一样

@King_RVM 方便问您一些具体的网格导入和修改边界条件的细节吗，因为我现在导入网格改边界条件总是失败