@zwenzhi 在 snappyHexMesh网格质量控制 中说： 多面体网格和切割体网格在同样网格数量上，精度方面哪个更好呀？谢谢～～ 这个问题不好说

@吴优 谢谢答复，最近也研究了一下，答主回答的方法是可行的，可是现在好像是没办法做到紧耦合，如果采用UDF的写结构求解程序的话，因为我们没有宏命令控制求解器在一个时间步内暂停，然后进行流体与结构求解的多次迭代求解，收敛后再向下一时间步推进。

完了，竟然现在不能翻墙了，你直接下载paraview-5.6.0的source包，拷贝到TP-7文件夹下试试？

@东岳 我没有找到李老师说的“tutorial: decayIsoTurb”。不过在OF里找到了这个前处理程序 正在看这篇文章，还看到了一篇中文的“一种改进的均匀各向同性湍流初始化方法” 我看的比较慢，有些基本概念还不太清楚。

@南风 他这个意思是群泡传质速率，不是单泡的。另外这个值与你其他参数取值也有关系

怎么样，这个后来算的

@D-Benjamin 具体一共都有哪些函数我也不是很清楚，具体函数和怎么用通常都是看看openfoam自带的边界或者求解器或者其它的应用有类似你想获得功能，就去翻源代码。gsum应该也可以对边界上的网格体积或者面进行求和，应该就能实现你说的功能。

@Peng 大神能否分享一下fluent 的udf实现方法，openfoam里我在低版本里可以用了(做了一些修改)，非常感谢

edit -> find data 如果数据只在网格或者点，可以先用filter point data to cell 或者cell data to point筛一下。

感谢分享

Really? 那得谢谢王博士！

CFDer的旗帜 4600！

插值只对patch进行，和是否并行无关。并行对local patch进行插值。

@cfd_fly 您好，请问文档还有吗？刚才下载看分享过期了

@TINGHAIK 完全可以，没有问题的。

@ibelief @ibelief 在 fluent dpm模型如何统计壁面有多少质量的粒子被吸附。 中说： @xulei-zju Report > Discrete phase > SAMPLE-> append file 会生成一个dpm 的文件 可以用excel txt 等工具打开，统计

@Izumi 这个没添加，我回去试试

@东岳 结果.tar 这个是发散的相分数，速度，压力图

@J 把你最近学习的成果总结一下，发上来嘛

为什么这个要交opeartor splitting呢？没理解。

@东岳 谢谢东岳大神

@东岳 偶然翻看到之前的讨论，结合东岳老师提到的声速有关的内容和最近的一些新体会，贴上来和大家讨论。 对于可压缩流动的定义好像存在不同的理解，我现在倾向于认为正确的是，是否严格意义的可压缩要看密度是否随压力变化。如果密度是随压力和温度、组分等变化，那么是完全的可压缩。如果密度不随压力变化，而只随温度、组分等变化，也是可压缩（因为密度会变化），但不是完全的可压缩。后者在很多传热、燃烧应用中是广泛存在的（密度受压力影响很小，密度差主要由温度和组分变化导致），在处理的时候大多采用所谓的低马赫数近似（low Mach number approximation）。其本质是假设密度不随压力变化，不一定是字面意思上的马赫数很低才可以使用，更多的应该是指代一套数值方法（这点不是很确定 ）。FOAM里面有一种状态方程是incompressiblePerfectGas，和普通的perfectGas区别在于，前者的压力用一个参考压力（常数），后者用的则是计算得到的压力（非常数）。 感觉FOAM里面之所以有基于psi的热物性和基于rho的热物性就是为了能够处理不同类型的可压缩流动。传热求解器（solver/heatTransfer）比如buoyantPimpleFoam采用了基于rho的热物性，而燃烧求解器（solver/combustion）比如reactingFoam，fireFoam采用了基于psi的热物性。燃烧求解器也有用基于rho的热物性的，目前还不清楚背后的原因。可压缩求解器（solver/compressible）中，像rhoPimpleFoam本来是采用基于psi的热物性，但高版本开始也把基于rho的热物性揉进去了，很好奇如果用rhoPimpleFoam采用同一套状态方程去求解同一个问题，但一个用psiThermo一个用rhoThermo，结果是否会一样？ 啰嗦了两大段，还是回到之前的问题。很明显reactingFoam和rhoReactingFoam两个求解器调用的压力方程形式不同是因为密度更新的方式不同，而密度更新方式的不同很可能来源于是否采用低马赫数假设这一套计算方法。这一点现在还不明确，所以发上来希望有大佬能给讲讲 另外有个问题，如果采用低马赫数近似，密度随压力不变，那是不是就相当于认为声速无穷大？这样假设会导致什么呢？

个人理解，如果不改transport property的话，引用的nu应该要加上phase1的关键字，比如可能要输出的应该是nu.air。关键是看nu在你用的solver里注册成了什么名字，这个名字可能和输入的名字完全不一样（注：是注册的名字不是变量名）

@东岳 多谢东岳老师指点 我再改一改边界条件算一下，搞定了再发帖，用这套软件的道友太难找了

会不会是安装破解出问题了，或者电脑显卡不太够？或者改一下背景颜色为纯白试试？

@东岳 好的，谢谢东岳老师的回答，另外想问一下老师，DPM中颗粒碰撞力的PairMOdel有没有什么参考资料？？谢谢

感谢分享！alias确实更方便！

@CiCi 最近怎么样

@xinxinzhou1995 可以交流一下吗，关于动态接触角模型

残差图

这个叫函数的重载，在面向对象的编程语言中是个普遍应用的特性。

@队长别开枪 大佬，教程有了吗

是欧拉-欧拉，确实如果引入边界条件的话不知道压力场怎么处理。谢谢分析，我再想想。加多孔介质主要就是物理上好解释，我现在觉得可以在wall的区域再附加一个指向wall外面的源项（根据具体问题可以指定方向），类似重力，以此抵消重力以及其它数值问题造成的穿透问题，相当于在这个源项下，wall区域的固体会被很缓慢地排出wall。我估计一开始如果就加上的话，可能就可以阻止固相进入wall，而且有多孔介质巨大的阻力系数在，排出的固相应该速度非常小，对wall外面的计算应该影响有限。

@东岳 那看来目前这个版本不支持ccmToFoam这个命令，@veen 的方法仍然可行，谢谢大家，祝大家节日快乐。

@mr-zhang 同问，楼主问题怎么处理的呢？分享一下呗~

@wwzhao 嗯嗯我已经找啊

@dyb 你可以参考自带的案例——airfoil2D，里面有升阻力系数监测的方法。

这个是不是用到了基于物理量+动网格技术。。。

@chengan-wang 没做过分段式的，单个翼型使用blockMesh生成结构化网格是可以的，网上有可以找到python或matlab用来生成blockMeshDict字典文件的脚本，你可以试试多段翼型的网格设计。

@chengan-wang 或许你可以搞一个局部细加密。

谢谢！！已好

这个估计很难了。你这个涉及流固耦合的问题。CFD这边做不鸟。但是你可以这么做，流场获得你需要的压强时程曲线，然后你在固体场中把压强加上去。这样是可以的。

波浪得求解器很多得。一般三个，waves2Foam，olaFlow, ihFoam。我现在可能开始用olaFlow

@东岳 感謝东岳老師的回覆.首先說個題外話,最近在英國這邊一個公衆號上看到您發的一片勸解博後選擇的文章.因爲是轉發的,起初不知道是您發的,後來越看越像您的口吻,最後在下面看到"轉自李东岳"的時候,心裏很是激動.謝謝东岳老師對後輩們的勸解! 我又做了一個用MPPICInterFoam在沒有粒子的情況下,繞柱(光滑壁面)的模擬.這次用了LES湍流模型.效果很好.並且並行運算時,LES比K-omega SST運行更快.大概是因爲這種情況更符合實際情況吧.跟實驗數據比較,甚至比之前interFoam做的還要好.圖中我做了MPPICInterFoam和interFoam的對比.上面的是MPPICInterFoam的.在柱後水流形成的湍流旋渦非常明顯.在不同高度的對比,效果也很好.至少說明,在水動力這塊兒,這個求解器配合LES,可以得到不錯的效果.當然,邊界條件也很重要. 有問題的大家一起討論啊.

@Samuel-Tu 非常感谢您的回复。讲的非常明白，我之前没有想到这些。其实，我想实现的就是一个热源点慢慢加热最终提高出口温度，或许我应该试试改成热源强度再试一试看。非常感谢您。

感谢分享！！！收录在此 https://www.cfd-china.com/topic/3499/openfoam小代码

@xpqiu 感谢，我去查一下!不过目前我用break也能顺利终止程序，并行和串行都可以

没有实际做过，但是这个链接提供的方法看起来很靠谱。https://www.physixfan.com/如何计算湍流能谱/

@东岳 非常感谢东岳老师！按照您的指导，我设置了新的边界条件immersedboundaryEnthalpyFvPatchField，其他部分的设置目前都正常，但是想要使用浸没温度边界上的函数来更新焓值时会报错，我的程序部分如下所示： const label patchi = patch().index(); wordList type = thermo.T().boundaryField().types(); if(isA<immersedBoundaryFvPatchScalarField>(type[patchi])) { immersedBoundaryFvPatchScalarField& TIb = thermo.T().boundaryField()[patchi]; scalarField & Tref = TIb.refValue(); } 之后程序就会报错， error: invalid initialization of reference of type ‘Foam::immersedBoundaryFvPatchScalarField& {aka Foam::immersedBoundaryFvPatchField<double>&}’ from expression of type ‘const Foam::fvPatchField<double>’ thermo.T().boundaryField()[patchi]; 应该是无法将基类转换为派生类，想和老师请教一下遇到这个问题要怎么解决呢，麻烦老师了！！！