题目就是想要问的，有没有大佬也在做这，求交流啊 。
我C++可以说基础等于无，囫囵吞枣看了一本陈志泊的《面向对象的程序设计语言》，明白基本原理，高级的技术都还不会呢。所以做这个代码就等于抄写。要是有也在做这个的，且恰巧没有人交流的，交流交流哈。

@李东岳 delta函数

@李东岳 谢谢李老师，计算之前controlDict已经事先挂载过了

libs ("test1_libIncompressibleTurbulenceModels.so");

@eason 如果是multi-body 问题，然后其中一个body需要IBM，就直接用Gmsh做其他body 的conforming mesh就可以。只有一个body且需要IBM，就可以做一个enclosed的正方体或者圆柱体，画网格就行。

@xpqiu 老师您好，我有一个关于监测多孔介质区域体积力的问题想请教下您，我在波浪场中前后设置两块多孔介质区域porosity1和porosity2，然后用function函数分别输出porosity1和porosity2的受到的体积力，但是两个输出的力是一模一样的，经过计算我发现输出的力是porosity1和porosity2的合力，也就是说openfoam没有办法分别输出porosity1和porosity2各自的受力？谢谢老师！

感谢分享。

在没有独显的服务器上，windows的paraview也是黑屏。

@李东岳 谢谢李老师的回复，这个问题我最后发现是一个很简单的错误导致的。因为decompsePar没有把边界场进行拆分:136: ，导致各个处理器文件夹下的场定义缺失所以发生了错误

@Mania 怀疑fft变换时参数选错了，OF不大可能错，因为你说之前的都没问题。

@Chenwei 这个确实是个问题。主要就是这个边界条件的组合会影响计算的结果。你有尝试过全部totolPressure么。或者进口速度给一个稍微小的值，如1e-3，然后流量吧这个减掉？或者干脆不剪掉，差得多么。

原始算例应该考虑的是稳定性的问题。

请教各位老师，Wilson面的概念和物理意义以及表达形式，搜索到下面的论文可是下载不下来，谢谢:xinxin:
Bohn D , Kerpicci H , Ren J , et al. Homogeneous and heterogeneous nucleation and wet-steam Flow in a nozzle and a LP-turbine cascade. 2002.

@李东岳 我这试了下system下创建fieldAverage文件并复制上述代码，删掉functions和fieldAverage1可以计算完成后处理，controDict不需要挂载，完整的文件内容如下

/*--------------------------------*- C++ -*----------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ / A nd | Version: 9 \\/ M anipulation | ------------------------------------------------------------------------------- Description Calculates and writes the time averages of given list of fields. \*---------------------------------------------------------------------------*/ type fieldAverage; libs ( "libfieldFunctionObjects.so" ); writeControl writeTime; timeStart 0; timeEnd 100; resetOnOutput false; fields ( U { mean on; prime2Mean on; base time; } );

:135: :135:

简单记录一下，两个问题
1.ocean-wave3d这个安装包没有下载好
2.fe40需要gcc5版本，虚拟机里没有统一成一个版本。

气泡群加入到槽道中，有周期性边界，让它们一直输运，我想要考察壁面上的脉动压力。

只有特别微观的槽道，10几个气泡的能算。这个我之前算过，算半个月也没算出来个统计平均值。后来就不整了。目前也没遇到过什么好办法。VOF算这种千万网格打底。除非换欧拉方法，特别快。

看看其他大佬有啥意见

Tryggvason那面算了非常多的这种东西，我估计你肯定知道。

@bestucan
先谢谢老师的回复
老师说的对。。。。

谢谢老师

应该是我的边界条件设置的不合理，我修改了物性和边界条件以后，目前是能继续算下去的。
贴一下我参考的网址：
https://www.researchgate.net/post/Can-anyone-see-this-problem-in-openfoamcan-anyone-solve-it/1

好的好的

我发现这个欧拉场到拉格朗日场，拉格朗日场到欧拉场的映射，有各种不同的方法..

如题，请教各位老师同学，我在LES模拟后，对数据进行后处理时，发现一个问题：定压比热的场值经过平均后被缩小了，但是其他场的平均值没有这个问题，请问各位这是怎么回事呢？平均后的峰值(8000)明显低于真实值(120000)，如图所示：
比热.jpg
比热云图.png

@李东岳 only 外力，就是被这个误导了才没转过来弯:136:

@vivian 我也遇到了一样的问题，请问您的问题解决了吗，可以指点一下吗？

大家好，
还有如下问题
/lib64/libstdc++.so.6: version `GLIBCXX_3.4.29' not found
怎么解决？

谢谢

@李东岳 谢谢李老师答复。我现在就在用-oversubscribe方法跑，但是我不知道这个方法具体是怎么调用的核数，比如是不是我的两个cpu都调用了，还是说只调用了一个cpu里的38个线程（一个cpu总共48线程）

三年了。继续膜拜blockMesh大佬

:haqi:

@田畔的风 在 有人在mac上装of么？ 中说：

@一颗鸭蛋 默认情况下，你编译求解器生成的可执行文件会保存到$FOAM_APPBIN，比如我的默认在<上级文件夹>/OpenFOAM-v2206/platforms/darwinARM64ClangDPInt32Opt/bin，但是这个App的加载磁盘是只读的，所以会报错。

解决方案有两个：

修改在Make/files中你的求解器保存路径 EXE = $(FOAM_APPBIN)/xxxFoam -> EXE = <新路径>/xxxFoam

然后在环境变量中为$FOAM_APPBIN附加这个路径，以便让OpenFOAM在执行时能搜索到这个文件

export $FOAM_APPBIN=<新路径>:$FOAM_APPBIN 创建一个可读写的分区，直接编译安装。参考 https://github.com/BrushXue/OpenFOAM-AppleM1。

您的意思是把openfoam安装到创建的可读写分区吗？我用的也是M1的Mac，现在想在Mac上编译自己修改的求解器，大佬能不能给点详细的指导？谢谢🙏

@李东岳 或许是我这边fluent有某个设置错了，杨老师论文中顶面用的slip结果不会出现k衰减。我再试了fluent中顶面设置no-slip，结果是比较符合预期，接近顶面壁面平均速度减小，湍动能增大
Ux剖面.png TKE剖面.png

@GenhouSUN 两台电脑，一个是9.5.0，一个是9.4.0。而且我别的版本安装不上，只在OpenFOAM-v1912版本上用，waves2Foam是2137版本

@Vortex 大佬太厉害了。用语言描述就是，地表速度低，科氏力弱，压力梯度导致风向偏移；高度上升后，气流速度增加，科氏力变强，将偏移拉了回来，并最终与压力梯度产生的压力平衡。

@李东岳 好的，明白了，谢谢李老师

@allanZhong 2年过去了。。。是需要乘以体积。例如我刚写了个科氏力源项：

nergySource { type coded; selectionMode all; field U; codeInclude #{ #}; codeAddSup #{ Pout<< "**codeAddSup**" << endl; const volVectorField& U = eqn.psi(); vectorField& USource = eqn.source(); vector omega(0, 0, 5e-5); const DimensionedField<scalar, volMesh>& V = mesh().V(); USource += (2.0*omega)^U*V; #}; codeAddRhoSup #{ Pout<< "**codeAddRhoSup**" << endl; #}; codeAddAlphaRhoSup #{ Pout<< "**codeAddAlphaRhoSup**" << endl; #}; }

这个方法很好,不过请注意,这个链接中的SurfaceVectors 应该是作为Seed Source 而maskPoints作为Input才对,作者笔误了.亲测可用.原文链接如下 https://zhuanlan.zhihu.com/p/592092340
还有一个带有箭头的流线图链接如下,不过我并没有试过哈,也贴出来.
https://blog.csdn.net/weixin_44571596/article/details/122187429

@李东岳 好的，我先尝试着修改下再去跑跑看

@李东岳 可以安装了

@李东岳 李老师，这是我的邮箱 coolclw67@163.com。之前没研究RANS，但挺感兴趣的，最近我也研究下，谢谢李老师

@李东岳 谢谢东岳老师，看了一下完美满足需求:146:

本人用sprayFaom进行一些模拟，请问如何获得液滴的蒸发速率？
质量方程：1676896511168.jpg
代码：

solve ( fvm::ddt(rho) + fvc::div(phi) - parcels.Srho(rho) );

论文中说Smass是蒸发速率，那我需要把代码中的parcels.Srho(rho)给输出出来，如何实现呢？

前辈，你好！请问sprayFoam中的液滴蒸发率怎么获得？

@wpx http://dyfluid.com/buoyantSimpleFoam.html

@jinlinna 我明天看看

@coolhhh 在 计算y+，Δx+，Δz+的疑问 中说：

但最后还是回归到网格除了法向外，另外两个方向怎么定义x和z方向

是的。一般sci里面怎么弄？我个人怀疑sci就是按照流向方向来了。

@学流体的小明 在 求助：openfoam三维算例中如何沿着某一个方向做空间上的平均 中说：

我算槽道流，是有一个专门的postChannel后处理程序的，可能你得自己编程做了。
可以参考一下这个后处理程序。
https://develop.openfoam.com/Development/openfoam/-/tree/OpenFOAM-v2012/applications/utilities/postProcessing/miscellaneous/postChannel

感谢，我发现了一个叫fluidfoam的后处理python工具，已经植入到openfoam中，需要py3.8以上版本，of的版本估计也得是2106以上，具体可以参照他的官方说明，该工具功能极其强大且便携

@coolhhh 这里的$S_k$指的是k输运方程的源项，应该不是0文件夹中的s

@Leeseungnan 在 openfoam里面的煤燃烧intrinsicRate model方程植入单位不统一 中说：

我发现两个变量的单位不一致，但是在方程植入的时候，程序里面却将两项做的差值处理，这种操作可以吗？

理论上是不可以的。但是还见过有类似的操作。我也很懵。举个例子，下面这个方程，ln()里面不能带单位，但是里面确是有单位的。类似下面这个公式在很多sci里面都有，我只是举例。搞得我很懵。

捕获.PNG

类似你那种两个变量相减，a-b这种，就更需要单位一致了。这个问题我认为你可以问一下openfoam.com那面：1）要么是他们单位写错了，2）要么是有问题

@卡洛 前辈方便的话能加下您的微信吗？想向您具体请教下，谢谢。我的微信是13652081916

@tidedrinker
我没验，你验证下，看到底用哪个理论。。。 深h对T影响有限，我更倾向深水用stokes1

还有，请教下：
你熟悉不规则波浪反射的估算吗？

@杨英狄 在 可凝气体的冷凝数值模拟 中说：

高温水蒸气在液态水中的冷凝及流动

不太好理解。水蒸气在水里冷凝然后流动？

后来这个解决了吗 :-)

主要是因为OpenFOAM那面多相流的工作量比较大，改进也比较多，新老版本都存在，所以就同时有了这么多东西。

twoPhaseEulerFoam: 最老版本的双流体求解器，稳健性一般 multiphaseEulerFoam: 最老版本的多流体求解器，稳健性不好 reactingTwoPhaseEulerFoam: 末期版本的双流体求解器，稳健性比较好 reactingMultiPhaseEulerFoam: 末期版本的多流体求解器，稳健性比较好 twoLiquidMixingFoam 不属于上面几类

想用IATE和类方法计算不含相变和传热的问题

用reacting类求解器

偶尔误入这个帖子，还挺有意思。

@李东岳 十分感谢李老师，后面课题组有培训经费一定要来上您的编程课:huahua:

你好。想问一下在openfoam中离散过程中对于面插值的中心差分格式的源文件代码openfoam/src/finiteVolume/interpolation/surfaceInterpolation/surfaceInterpolationScheme/surfaceInterpolationScheme.C
其中对于点乘面差值格式的代码形式主要由Foam::surfaceInterpolationScheme<Type>::dotInterpolate 函数决定的，我现在的主要想法是想通过修正插值格式在代码sfi[fi] = Sfi[fi] & (lambda[fi]*(vfi[P[fi]] - vfi[N[fi]]) + vfi[N[fi]] ) 后面加体心梯度项点乘偏移量向量，如图所示：
b2b2b5cc-a725-48e9-ad35-ba8bc36d7208-1675996388057.png

我该如何修改代码，或者说大致的思路，请各位指点，拜托了。

@李东岳 是的李老师，这个功能挺有意思的。我前不久写过Pope和Celik Index的类似后处理代码，看到之后新功能就分享给大家了。当然我实现的比较简单，适用性没有那么好，但是足够用了。您说的写到controlDict的思路也可以，这个其实也就是个简单验证，也不需要多优雅～