@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; #}; }

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

@李东岳 可以安装了

@李东岳 李老师，这是我的邮箱 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）要么是有问题

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

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

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

主要是因为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的思路也可以，这个其实也就是个简单验证，也不需要多优雅～

@是嘿嘿牛啊 haha 理解理解

你在机械硬盘跑，数据少的话，不影响跑的速度