@李东岳 谢谢李老师，我试过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能自动算雷诺数吗？

感激不尽

@bestucan 对 用tecplot成功了 tecplot里 reverse distance的插值方法可以把无序的数据插值到有序的网格上 不过有几百组数据 还得弄个宏什么的

这个用interFoam+动网格不就可以么

@浪迹天大 非常感谢！今天仔细看了下，终于搞定了。

对于重叠网格运动网格方面做了新的尝试，在算例tut/multiphase/interFoam/laminar/sloshingTank3D6Dof中拷贝gen6Dof文件夹进入自己的background模型文件夹，/constant中的dynamicMeshDict也要拷贝到自己的/constant中，如果只想重叠网格中Zone1（模型）网格运动 需要把dynamicMeshDict改成如下，当然之前的topoSetDict也要增加movingZone的部分

dynamicFvMesh dynamicOversetFvMesh; solver multiSolidBodyMotionSolver; movingZone { solidBodyMotionFunction tabulated6DoFMotion; CofG (0 0 0）； timeDataFileName "<constant>/6Dof.dat"; }

下面说下gen6Dof文件问题：

先进入/Make中

gen6Dof.C EXE = ./gen6Dof

回到上阶文件夹，更改gen6Dof.C文件 ，更改所需运动，6个自由度随便改改，然后终端输入wmake编译，完成之后终端输入 ./gen6Dof,全部没有问题之后把生成的6Dof.dat拷贝到constant文件夹中，然后回到上阶文件夹运行 movingDynamicMesh（别忘了改controDict文件的时常/6Dof时常,然后就可以观察 paraFoam了。
大致思路就这样，我也试了的确可行。
如果有相似重叠网格中运动网格类问题可以参考下