@hongjiewang 我直接好家伙，这都两年了居然还回复我。 由于种种原因我现在已经不怎么玩OF了，只能给你一些比较通用的建议。首先看这一堆问号大概是segment fault错误，那么我首先猜是产生了异常值导致求解器爆炸,log函数最容易出现的问题就是传入0值和负值，建议你首先检查数值问题。如果不是数值问题，那就单步调试之类的慢慢debug看看怎么回事，你的两个变量是否正确初始化了等等

自己编译求解器很难调试，因为求解的方程都是耦合在一起的，一般发散会互相影响，不好判断问题在哪里

麻烦帮忙改正一下。想要规定thetaBoundary为一个常数值，然后对theta（相分数）进行对比，总是实现不了。谢谢。 const Field<scalar>& theta(); const Field<scalar>& thetaBoundary(); thetaBoundary = {0.5}; if (theta < thetaBoundary)

@zhaodl 这个目前我尚不清楚

算出来的等值面结果如图

Hi，东岳大佬，请教下openfoam里可以实现分区网格计算，比如两个分区，并设置interface boundary condition来coupling这两个区域吗？如果可以，麻烦给点思路，谢谢。 我暂时想到的只有fluid-solid interaction，即FSI toolbox in Foam-extend，改写solid domain为fluid-fluid interaction，来实现，不过有些繁琐了。 BR， Shiyu

Hi, Number44 我遇到和您类似的问题，不过我的case简单不需要并行，但对时间cost要求很严格，我发现我改动的FSI很花费时间，可能与你说的插值函数有关，不知你有改善效率问题吗？ 其次，想问一下，除了foam-extend里面有FSI tool，基金会的openfoam有类似的吗？比较偏向于使用基金会版本，谢谢

@Henry-Zhou 您好，我最近也想模拟固液两相沉积的，这方面有啥进展或者思路吗？

请问如图左边为进口右边为出口，群体平衡模型能否适用于这种结构的模型?

@evensun 能在高校待时间长的，必然都有外快。各有各的路子，有些大家能看见，有些大家看不见。否则待个5 6年必然就走了。

大四学生表示人是在元旦连夜跑的（哈尔滨疫情影响）

@hongjiewang 能帮到你我也很开心，继续加油~

@benqing 我用的2.4.0 设置应该是没有问题的，按手册设置的。好的，我再看看，谢谢

@李东岳 是的，但是我没有在ICEM import Model里面找到相应的*.ansa格式呀，应该用ansa导出什么格式呢？

混个c币

@李东岳 谢谢东岳老师，我还是在win下安装paraview吧~

最近在备课，准备了很多湍流的东西，$\epsilon=\nu \overline{\dfrac{\partial u_{i}^{\prime}}{\partial x_{k}} \dfrac{\partial u_{i}^{\prime}}{\partial x_{k}}}$ 这个定义不是真正的epsilon的定义，真正的定义应该是$\epsilon=2\nu \overline{S_{ij}'S_{ij}'}=2\nu \overline{S': S'}=2\nu |S'|^2$。所以在写代码的时候，大体就是： 通过icoFoam获得U， 通过postProcess获得UMean 二者相减是U' 通过U'获得$S'$ $\epsilon=2\nu |S'|^2$ 对$\epsilon$做平均 大体上，需要做两次平均过程。第一次平均获得uMean，第二次获得$\epsilon$定义里面的平均。这样可以获得解析的$\epsilon$。同时还要加上$\epsilon_{sgs}$。 可以参考Wilcox那本湍流书的109页。 我算了算$Re_\tau=395$的算例，左边是湍流动能耗散率，右边是速度。看起来是符合物理的。就是没有数据可以对。

@lys 用的已有模型，专门针对非牛顿流体的曳力模型还没用过，希望能进一步交流啊

thinkpad t480s i7 8550u 16G内存 平时写论文、做ppt用，也可以开虚拟机改求解器算算小规模问题，真要数值计算还得远程服务器

有用过Foam-extend 4.0 里的interIbFoam的吗？这个求解器是不是不能加湍流模型啊？一加湍流模型边界条件就报错

谢谢李老师，确实是alphatWallFunction出了问题。重新修改了这个函数后，就可以了。再次感谢。

各位老师好，关于openFoam模拟颗粒相关的问题： （1）open foam是否可以模拟颗粒-颗粒以及颗粒-壁面的黏结效应，颗粒在壁面的黏结系数是哪个参数？算例文档中只有颗粒与壁面的反弹与损耗系数。 （2）颗粒类型在哪里选择，目前至看到parcelTypeID 1，有其他形状的颗粒类型吗

感谢各位大佬！

@zhe 因为最后一步你重新给p.U()速度了啊。。

@freeazma 如果我想在 fluent 里面同时查看/获取/输出固体和流体温度，能实现吗？

这几天备课OKS，今年大幅度增加湍流的内容。前几天也在朋友圈发了个图。简单说一下。LES处理壁面几种方法： LES-RANS混合方程：用RANS求解壁面区域，LES求解主流区。也就是你需要调用DES、DDES等之类 Wall-modelled LES：对壁面区域进行模化，主流区采用LES。这个进一步还区分不同的模型，比如代数类、PDE类、ODE类等。 http://www.dyfluid.cn/theory.pdf 混合方程我在这里更新了一点内容。 另外没理解你说的壁面雷诺数什么意思。 Ref: Journal of Computational Physics 305 (2016) 589–603

各位大佬好，本人是OF小白，毕设使用OF研究氢氧燃烧现象，求解器采用reactingFoam，在thermophysicalProperties设定中需要加入惰性气体N2，OF提供的reactingFoam燃烧案例均只设定inertSpecie N2即可，但在跑的过程中出现keyword N2 is undefine的报错，试用banaba法得到以下反馈,，似乎是需要找到N2的关键词，有大佬知道应该如何修改嘛``` Attemp to return primitive entry ITstream:/home/user/OenFOAM/run/EX/cinstant/thermophysicalProperties/N2,line 32 ,IO stream: Vertion 2.0 ,format ASCII ,primitiveEntry 'N2' comprises on line 32 the word 'banana' as a sub-dictionary From function virtual const Foam:dictionary &Foam:primitiveEntry : :dic() const in file db/dictionary/primitiveEntry/primitiveEntry.C at line 190. FOAM aborting

http://dyfluid.com/docs/install.html 按照这个安装 不会有任何问题

我用interIbFoam求解有结构物的明渠流问题，在应用了k-epsilon湍流模型后，nut文件的浸没边界条件出现如下报错： --> FOAM FATAL ERROR: evaluate() cannot be called for a genericFvPatchField (actual type immersedBoundaryWallFunction) on patch vegetation of field nut in file "/mnt/g/foam/thomaschi-4.0/run/cylBumpInterIbFoam/0/nut" You are probably trying to solve for a field with a generic boundary condition. From function genericFvPatchField<Type>::evaluate(const Pstream::commsTypes) in file fields/fvPatchFields/basic/generic/genericFvPatchField.C at line 760. FOAM exiting 有人知道是什么原因吗？该怎么解决？

@荷兰后腰阿贾克斯 嗯嗯，谢谢老哥了

Thanks a lot. I wonder if the partial elimination algorithm mentioned in Passalaqcua's paper was implemented in the current version?

Fluent模拟相变问题，固液密度变化或者让液体密度随温度变化，就会出现压力线性修正，然后回流发散。请问这有解决办法吗，

parcels.theta()().write(); 后处理通过theta场和网格体积求parcel的体积。 总质量不就是书橱文件中的mass introduced么

@我是河滩 在 指定cellzone区域添加动量源项的问题 中说： 遍历了整个网格区域进行搜索 前辈您好，您是如何遍历整个网格区域并获取网格坐标的呢，是利用mesh.C().component()语句吗

感谢解读，学习到了@李东岳

Error: received a fatal signal (Segmentation fault). UDF 代码： #include "udf.h" #include <math.h> DEFINE_MASS_TRANSFER(liq_gas_source, cell, thread, from_index, from_species_index, to_index, to_species_index) { real m_lg, K, R; real YI_SAT = 1.5e-6; Thread *gas, *liq; gas = THREAD_SUB_THREAD(thread, from_index); liq = THREAD_SUB_THREAD(thread, to_index); m_lg = 0.0; if((C_YI(cell,liq,1)>YI_SAT)&&(C_VOF(cell,liq)>0.)&&(C_VOF(cell,gas)<=1)) { R=45.833e-6*sqrt(CURRENT_TIME-30.0)+0.833e-6; K= 7.38e-9*(1/(R)+1/sqrt(3.14157.38e-9(CURRENT_TIME-30.0))); m_lg =-1000K(C_YI(cell,liq,1)- YI_SAT); } return (m_lg); }

@veen 从源码编译of是个挺痛苦的事情，建议不把of安装在root，可以考虑试一试easybuild，它是基于module的环境管理，应该比较成熟了。我从0开始，最开始在PC上面编of-3.0，到后来在实验室的机群上home编译of-4.0, foam-extend等等，多版本的管理用module还是方便些。印象中，gcc（如果你用gcc的话）和mpi的版本很重要。

已解决：使用foam-extend4.1版本中的compressibleInterFoam求解器，模型修改为球对称模型，计算正确。 P.S. foam-extend4.1中的compressibleInterFoam没有耦合温度方程，计算过程中没有出现上述问题。 参考文献：Numerical modeling of laser generated cavitation bubbles with the finite volume and volume of fluid method, using OpenFOAM

@李东岳 感谢老师关心！已经在家啦，趁着有空补补基础hhh

@Samuel-Tu 感谢

最近对两个概念非常混淆，即摩擦系数和切线回弹系数，这两个系数是否就是同一个概念？ 比如在以下的颗粒pairCollision参数设置中，法向碰撞的回弹系数根据alpha设定。但是切向的回弹系数，是不是就是这个mu? pairSpringSliderDashpotCoeffs { useEquivalentSize no; alpha 0.01; b 1.5; mu 0.1; cohesionEnergyDensity 0; collisionResolutionSteps 1; };

@波流力 上传了所有案例文件，需要的请自行下载。见桥梁水动力： https://mp.weixin.qq.com/s?__biz=MzU5NTgwODcyMg==&mid=2247483815&idx=1&sn=a0ef629459736bde6b94ebf0d31f68b1&chksm=fe6d191ac91a900cd2a2c2d45a6aafb711da8aab3f3fc7ad7b1d40969df1cc1f0f741080bc10&token=1944678271&lang=zh_CN#rd

压力速度迭代感觉国际上这面玩的少了 就刚出那时候搞出来一些SIMPLEC之类的

回档了 这是什么意思..

@同学博 好的 十分感谢！

受教了，谢谢

@Zhy2022 不好意思，8中也有非2D的定义-SymmTensorI.H

感谢各位老师和同行的鼎力相助，目前入口最高压力能够达到30MPa，和实验匹配效果较好。 最后再一次感谢大家！

各位朋友们好！ 我最近在看lagrangian库中的solidParticle库代码，对于里面的 类模板的使用 不太懂，solidParticleCloud.C部分代码如下： void Foam::solidParticleCloud::move(const dimensionedVector& g) //在solidParticleCloud类中定义一个move函数 { const volScalarField& rho = mesh_.lookupObject<const volScalarField>("rho"); const volVectorField& U = mesh_.lookupObject<const volVectorField>("U"); const volScalarField& nu = mesh_.lookupObject<const volScalarField>("nu"); interpolationCellPoint<scalar> rhoInterp(rho); interpolationCellPoint<vector> UInterp(U); interpolationCellPoint<scalar> nuInterp(nu); solidParticle::trackingData td(*this, rhoInterp, UInterp, nuInterp, g.value()); //在solidParticleCloud类下面定义了一个trackingData类的对象td？？ Cloud<solidParticle>::move(td, mesh_.time().deltaTValue()); } 代码中的这句话： solidParticle::trackingData td(*this, rhoInterp, UInterp, nuInterp, g.value()); 目前自己很困惑，为什么在solidParticleCloud类的函数move中可以 定义一个 solidParticle类中的 trackingData类的 对象td呢？ 自己查找了solidParticleCloud.H文件发现 solidParticleCloud的基类是 Cloud<solidParticle>类，这个Cloud<solidParticle>类 中solidParticle是实参，对Cloud <ParticleType> 类模板进行了实例化. 我怀疑难道是因为Cloud<solidParticle>类中的 实参是solidParticle，所以其子类solidParticleCloud中函数 就可以调用/定义 solidParticle类中的trackingData类？ 请哪位高手能够指点一下呢？感激不尽