@东岳 非惯性坐标系啊，体积力与每个网格的具体位置有关，这样就不用动网格了

并不是这么绝对。就像是密度基求解器和压力基求解器一样。

理论上，rhoCentralFoam是一个密度基求解器，原本用来求解欧拉方程，方程具有双曲型，不同的特征值。问题在于求解亚音速的时候，方程特征会变化，波动速度传输趋于无穷大。rhoCentralFoam调用的格式还取决于波动速读会导致时间布长太小。可以看看中心迎风格式的原始资料，讨论都是基于双曲型欧拉方程的（具备不同的特征值，非不可压缩） https://www.jianguoyun.com/p/De-T2z4Q9s3ZBhik7bkB

不过，我并没有测试过使用密度基中心迎风格式求解亚音速问题，从理论上来看，这种收敛性很差。不过，我觉得你们用一个一维案例测试下算法，比如试试用rhoCentralFoam求解1D喷管，逐渐降低速度看看收敛性，再应用你们的实际问题更好，欢迎反馈测试结果！This is scientific research :ok2:

@东岳 :chigua:

:wocao:

Screenshot_20190429_094156.jpg 这两行后面缺斜杠

@浪迹天大

你这个问题很有意思，从你的推导分析来看，SuSp 确实是削弱了对角占优。

但其中有一处错误：V30 应该为 V-30。不过这改变的是源项，而与对角元素无关。我们一起分析下是怎么回事。

如果将方程按照这种形式写：$\frac{\partial k}{\partial t}=-(-A)-(-B)$

那么代码可以写成：== -fvm::SuSp(-A/k, k) - fvm::SuSp(-B/k, k)。

实际计算时，假设 A=-30，B=50。
那么 $-A/k>0$，$-B/k<0$，fvm前面的负号与 == 抵消，代码实际执行的是

fvm.diag() += V*30 fvm.source() -= V*(-50)*k

这又的确增加了对角占优，显然 fvm::SuSp 不能随心所欲地用！

-fvm::SuSp(-A/k, k) 和 fvm::SuSp(A/k, k) 的结果不一样，而前一种写法才是我们想要的结果。

我翻了 OpenFOAM 的代码后发现，凡是在 == 右边使用的 Sp 和 SuSp 前面必有负号，凡是在 == 左边使用的 Sp 前面必为正号 (代码中基本没有出现 fvm::Su)。

看来 SuSp 是个深坑，用起来要非常小心才行。

提供点具体的错误信息老铁 :laile:

多谢老师

哈哈哈，没事，谢谢东岳老师！

谢谢教授:xinxin:

@东岳 是的，我要用到坐标三个方向的具体数值，标量

@史浩 我不是大佬。最近在做一个VAR的项目，考虑温度场-电磁场-流场-组织耦合计算，需要自编程序，有程序验证的问题。对OpenFoam，我没做过这方面的案例，但我想考察2个方面：一是看一下OpenFoam与自编程序的结果对比，二是看一下OpenFoam与自编程序的计算速度对比。
ESR和VAR比较类似，有相通的地方，你如果做了OpenFoam算例的话可以交流一下。

@史浩 xiexie

在这里讨论buoyantPimpleFoam，rhoPimpleFoam和buoyantPimpleFoam大同小异

首先，buoyantPimpleFoam解析中的方程(3)中的$A_{\mathrm{P}}$包含密度$\rho$；

1/Ap也是除掉密度$\rho$的

不可压缩算法和可压缩算法的压力方程不一样，可压缩算法的压力方程需要乘以$\rho$，不可压缩不需要，buoyantPimpleFoam解析中的方程(15)都需要乘以$\rho$，你看那个$\frac{\rho^{n}}{A}$，所以，

你会发现surfaceScalarField rhorAUf("rhorAUf", fvc::interpolate(rho*rAU))中出现了rho*rAU

@桎梏 我没有报错:quwan:

@Dingcy 谢谢

@夜阑烟寒 您好，我是把这个曲面展开成这个样子的，之前的在tecplot里面变换了坐标，还是曲面，想变成这种平面
QQ截图20190423005626.png

@Dingcy 哈哈我是在用低雷诺数湍流模型算超临界CO2流动传热的例子，物性变化很剧烈，不好收敛。

@东岳 就是一个二维方腔的自然对流，但是不太清楚这个温度边界条件该如何设定

入门不久，尝试着答一下。
我认为openFoam中的边界条件大致分两类：

强制性边界条件，直接规定数值，如fixedValue，fixedMean，这类边界条件的信息会一直存在于求解所得的结果文件，算很多步以后查看文件就会发现边界条件并没有变化； 参考性边界条件，这类边界条件大多是第二类边界条件，如出口定义inletOutlet， 对于流出的气体而言，该边界其实是zeroGradient，这时定义的value只是初始迭代的时候给的一个初值，该边界的数值由流场内部的点插值得到