@李东岳 李老师，这是我的邮箱 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 理解理解

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

@李东岳 有两个问题，第一，计算稳态流场的情况下buoyantBoussinesqSimpleFoam可以采用这种方式吗？第二，这种方法是不是没计算离散相方程。

用多重网格算法求解线性方程组Ax=b，算法求解过程中按理只需要给定系数粗细网格的系数矩阵就可以了，但在看算法理论时，都涉及到粗细网格的位置关系的判定来确定插值矩阵、限定矩阵。这个就需要在算法上导入网格具体的物理位置了吧。这样好像不对（按理是只要给定系数矩阵A，b,就能求解Ax=b）对于这点一直想不明白:mihu:

@hy1112006 是的李老师，在cfd-online上找到试了一个方法，虽然有系统警告，但暂且管用：
cafe1fe5-612f-45bc-adac-914ea277d477-image.png

警告如下，目前不影响运行：
7997c7d3-2fcc-47c1-a44b-c92d5a4c78f1-image.png