icoFoam的一些细节问题

dyj19901127

@李东岳按照《流动与传热数值计算—若干问题的研究和探讨》第2.8节和Discussion on momentum interpolation method for collocated grids of incompressible flow. Numerical Heat Transfer: Part B: Fundamentals，choi格式与亚松弛因子无关，仍然与时间步长相关。

李东岳

@dyj19901127

参考choi的论文，其引用Majumdar的文章表明"the resulting converged solution is relaxation factor dependent". 具体原理，还需研究。

我刚看了Yu bo的Discussion on momentum interpolation method for collocated grids of incompressible flow. Numerical Heat Transfer: Part B: Fundamentals，他也是说相关的。他们引用的Majumdar和Miller的文章。确认下？

Anyway，多研究研究，看看choi原理是怎么实现的。Choi那个论文很短。Yo bo师从陶院士？不一般啊。

队长别开枪

@李东岳
我觉得关键点在于进行当前时间步Rhie-chow动量插值时上一时间步界面速度$\mathbf{U}_{f}^{n}$的计算，如果没有ddtcorr()，$\mathbf{U}_{f}^{n}$只是简单通过相邻单元基于动量方程系数矩阵主对角元素加权平均得到，如果加了修正项，那么$\mathbf{U}_{f}^{n}$使用的就是上一时间步里基于动量插值得到的结果，所以在phiHbyA的计算里需要修正界面的体积通量，下面的代码是EulerDdtScheme.C527~545行（OpenFOAM 2.4.0）中的一部分

fluxFieldType phiCorr
(
      phi.oldTime() - (mesh().Sf() & fvc::interpolate(U.oldTime()))
);

- (mesh().Sf() & fvc::interpolate(U.oldTime()))用于消去计算HbyA时加上的体中心处的旧值。修正的原因是即便体中心处的速度场满足连续方程，基于动量方程系数矩阵主对角元素加权平均往面中心插值，时间步长越小，越倾向于简单算术平均，再加上非结构不规则网格等因素，很难再保证插值后的面速度场仍然保持0散度条件。
个人浅见，欢迎批评指正。

yhdthu

@李东岳前辈好，我有个问题，这个UEqn.H()分为两个部分，一个是AnUn，另一个是源项S，在同一个时间步内，随着修正次数变化的只是AnUn，但是S是不变的对不？

李东岳

S不一定，如果S和求解的变量有关，那肯定变啊。

yhdthu

@李东岳前辈，关于动量方程中的扩散项，由于其处理方式为半隐半显，即：

tmp<fvVectorMatrix> laminar::divDevRhoReff(volVectorField& U) const 
{
  return
        (
          - fvm::laplacian(muEff(), U)
          - fvc::div(muEff()*dev2(T(fvc::grad(U))))
        );
}

我目前认为这个显式部分是进入了.H()算子中，我想问的就是这部分是不是在压力修正时（即同一个时间步内）的计算不变的？

李东岳

我目前认为这个显式部分是进入了.H()算子中

这个目前我不确定，你可以看一看。:expressionless:

我是河滩

@东岳对流项是不是对UEqn.A()没有贡献？因为单元面上通量phi求和为0，满足连续性方程。

我是河滩

@CFD中文网压力全为零梯度边界，压力泊松方程离散后的系数矩阵奇异，adjustPhi修改了矩阵，让矩阵不再奇异？

李东岳

@我是河滩在 icoFoam的一些细节问题中说：

对流项是不是对UEqn.A()没有贡献？

取决于你的离散格式，大部分格式有贡献

我是河滩

@东岳是的，还和网格有关，均匀网格和非均匀网格产生的结果也不同。