@沐风 您好! 我最近刚开始学习OF。关于您的 文献1的validation，我也想试一下，请问在什么标准求解器上修改比较好？想要应用boussinesq 假设，但看buoyantBoussinesqPimpleFoam 求解器描述为湍流的。但该算例的自然对流Ra=6*10^5,即为层流。希望您指教一下，非常感谢！！

感觉是不是fluent有什么黑科技，在fluent里能算，在of里发散了。。

补充一下~算例用的求解器是自己在原求解器基础上修改的，这个求解器已经可以确定没有语句或者逻辑上的问题了。之前用这个求解器已经正确跑过算例了，结果也是和物理相符合的。只是前两天给虚拟机扩容，操作失败，然后所有的算例文件还没来得及保存。这两天新设的这个算例总是发散，我觉得是边界条件的问题，但是又不太懂，所以还是得麻烦各位老师~

大家好，最近我在用双流体模型做一个沉积案例，其中一相由于重力会沉积到壁面，并且附着在上边，想请教一下各位大神，如何设置壁面类型或者某一相的参数，才能出现附着的效果。谢谢！！！

@李东岳 谢谢李老师，我试过mesh(),还是过不了。这相当于一种以cell为计算单位的简化。在一些粒子数量级较大的情况下，可以应用。

@veen 好的。等你的教程，嘿嘿

@xpqiu 谢谢邱老师的热心指点！我明白了！再次感谢:xinxin:

//- Interactions between real and referred (off processor) particles void realReferredInteraction();

请问，在pairCollision中提到的referred particles 具体指什么？谢谢。

@jyang
在decomposeParDict文件最后，加入：
constraints
{
constraint1
{
type preservePatches;
patches (region1_to_region2 region2_to_region1);
}
}

@bestucan 和论文原作者确认了，说我这样是等价的，调整了一下其他地方，目前算出来对上了。。

@xpqiu 对，nutW2是分配了空间了的，大小就是cellID的大小。
造成这个的原因我怀疑是不是因为的我边界有转折导致的。。我是提取的roof这个边界上对应的第一层网格的值的。
121.png
因为当时我做计算的时候除以了nutW1之后计算报错了，我才发现nutW1[0]=0，不能被除。然后才改成使用了nutW2的

OpenFOAM有个快捷命令

src

进入源代码所在文件夹，然后

grep -R 'ParcelType::trackingData' ./

就把所有文件抓住来了，如果找定义，只看.H文件

preCICE很amazing，用好了应该比自己改求解器简单高效 https://www.bilibili.com/video/BV1iT4y1P7e7/
（现在发帖竟然需要间隔300s:jingya: ）

@西湖冷月 借帖求教大佬，狭窄裂缝中的湍流用哪种方法模拟好呢？

:xinxin2:

OpenFOAM v2006中，adiabaticPerfectFluid状态方程如下：

$\rho (p)=\rho_{ref}(\frac{p+B}{p_{ref}+B})^{\frac{1}{n_{T}}}$

根据熵的计算公式：

$dS=c_{p}\frac{dT}{T}-(\frac{\partial v}{\partial T})_{p}dp$

求解器中将这一公式分为两部分进行积分计算，其中第一部分$c_{p}\frac{dT}{T}$ 的积分计算在thermophysicalModels/specie/thermo/hConst/hConstThermoI.H中实现，即$c_{p}*log(\frac{T}{T_{std}})$
但第二部分的积分，由于公式是绝热状态的，没有温度作为参数，想请教一下官方求解器中为什么会进行以下处理：

通过源代码公式的求导反推，官方求解器在计算第二部分积分时，在状态方程的分母上直接添加了温度T，如下式：
CodeCogsEqn (1).gif
公式源代码详见adiabaticPerfectFluidI.H文件的S函数：

template<class Specie> inline Foam::scalar Foam::adiabaticperfectFluid<specie>::S ( scalar p, scalar T ) const { scalar n=1-1.0/gamma_; return -pow(p0_+B_,1.0/gamma_)*(pow((p+B_),n)-pow((Pstd+B_),n)) /(rho0_*T*n); }

除上述绝热状态方程外，线性变化的状态方程linear也采用了同样的处理方式，想知道为什么可以这样处理？
另外一个问题是，如果采用上述处理之后，对于气体的$c_{p}-c_{v}$，根据公式:
CodeCogsEqn (2).gif 该如何处理？

那创新性不是拉满了，厉害了！:146:

做一个点线图，横轴是孔隙率，纵轴是最大液相分数，看看到哪里突变了。二分法找这个突变点，应该很快。

出图的时候用同一标尺（最大值最小值相等），这样不同图之间也有可比性。看看那个刚好导致突变点的孔隙率的case，液相分数是怎么消失的，是跑出边界了还是在流场内被计算过程吃了

不知道你对QBMM这面的研究进展怎样了:135:

@Samuel-Tu :
谢兄之玉言，弟牢记

然兄亦暗示：
有的参数需调整，以保证物理意义一致。

另外，如兄及各位大学士不弃，望再赐一教：
OF能自动算雷诺数吗？

感激不尽