@浪迹天大 感谢大佬回答，我是这样设置的，但是不行，具体原因我还在找。

请问一下，有没有基于压力的chemistryModel，现有的openfoam中只有两种化学反应模型，即standardChemistryModel和TDACChemistryModel，我想植入一个基于压力的化学反应模型，有没有提供参考的现有库

@Voynich 你试试，期待你的计算结果

@李东岳 东岳老师，我将代码改了一下，出现了如下错误，您是否知道什么原因呢？1618388889(1).jpg 1618389113(1).jpg

请问各位老师，各位前辈，我想在chtMultiRegionFoam求解器添加热解模型，用来观察大颗粒的热解行为，希望各位不吝赐教。谢谢。

@Samuel-Tu 您好，将您的语句添加在我的case中，其中只修改了patch的名字，结果出现下面的错误，请问应该怎么修改。

#0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigSegv::sigHandler(int) at ??:? #2 ? in /lib/x86_64-linux-gnu/libpthread.so.0 #3 Foam::fvPatch::Cf() const at ??:? #4 Foam::writeCentersFunctionObject::write() at system/controlDict.functions.writeCenters:74 #5 Foam::functionObjectList::execute() at ??:? #6 Foam::Time::run() const at ??:? #7 ? in ~/OpenFOAM/dyfluid-v2006/platforms/linux64GccDPInt32Opt/bin/hybridPorousInterFoam #8 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6 #9 ? in ~/OpenFOAM/dyfluid-v2006/platforms/linux64GccDPInt32Opt/bin/hybridPorousInterFoam Segmentation fault (core dumped)```

gradU输出格式是c34a2688-3643-403a-9b10-1c6d0d510b66-image.png

@veen 感谢您的解答！reconstructParMesh时每个时间步都显示No mesh。
用 touch case.foam && paraview case.foam或者paraFoam -builtin直接处理的时候，在中依然没有物理场信息，不知道我哪里操作有问题了，我用tecplot试试

@wwzhao 好的我明白了，您说的外循环指的是pimple算法的循环，内循环指的是线性方程组的迭代求解，pimple循环中每次都要求解一次压力方程，一共求解2-4次左右，是吗？谢谢老师！

@xpqiu :
谢谢老师您的精彩点评，师之教诲，学生记下来
谢谢老师
@veen : 也谢谢veen

@aiweimo

UEqn.relax() ，即使松弛因子是1，其实也会对 UEqn 的 系数矩阵进行一些调整，以使它更满足对角占优。fvSolution 里面如果没有定义松弛因子，则默认值为1。

@aiweimo

k-w方程是低雷诺数方程，本身不需要壁面函数

严格来说，应该是 $k-\omega$ 方程可以一直算到 $y^+<1$，但是，$k-\omega$ 方程也可以结合壁函数来算$y^+>30$ 的网格。相比之下，标准$k-\varepsilon$模型的$\varepsilon$ 方程则只在log-law region之外才严格有效，不可直接用来算很小 $y^+$ 的网格。
OF 里面的那些湍流壁面边界条件，有一点容易引起误解的是他所有的名字都带 wallFunction，而大部分人看到 wall function，想到的都是标准壁函数，也就是 log-law region的 $U^+=\frac{1}{\kappa}\ln(Ey^+)$ 以及层流底层的 $U^+=y^+$ 。 实际情况下原比这个复杂，因为很难保证整个region 都满足 $y^+ > 30$ 或者 $y^+<1$，所以，对于实际问题，湍流模型的壁面处理需要考虑如何处理$1<y^+<30$ 的区域。文献中有很多相关的研究，openfoam 的 omegaWallFunction 其实就是一种，里面构造了一个blending function 把适合 $y^+ > 30$和 $y^+<1$的 omega 公式blend一下，来让omega 在任意 $y^+$ 下都能算出来一个相对合理的值。
所以，OF 里面那些 wallFunction 的边界，可以理解为定义多种不同的 wallTreatment 方式，具体怎么用，需要结合着实际情况来。2 楼提到的就是一种常见的用法。文献中壁面处理的方式太多太多了。。。

至于 $\omega_{wall}$ 指定为零梯度或者固定值，这种做法不太合理。因为 $\omega$ 在壁面附近需要满足某种渐近特性（asymptotic），零梯度或者固定值都不满足这个条件，虽然可以可以算出来结果，但是不太符合物理。

我最近看了下，

替代文字

上面这个方程应该是稳态算法的方程。里面MRF.DDt(U)用于处理$\Omega\times\bfU_I$。

96bfbd01-3d8c-469a-a303-a79d75b4d80c-image.png

但是上面这个方程不是OpenFOAM植入的方程。在瞬态情况下应该是下面这个：
$$
\frac{\p \bfU_I}{\p t}+\nabla\cdot(\bfU_I\bfU_R)+\Omega\times\bfU_I=...
$$
下面的瞬态代码分别对应：

tmp<fvVectorMatrix> tUEqn ( fvm::ddt(U) //第1项 + fvm::div(phi, U) //第2项 + MRF.DDt(U) //第3项 + turbulence->divDevReff(U) == fvOptions(U) );

:chitang:

@Samuel-Tu 在 local index to global index 中说：

但是写入solver的时候调用不了这个函数。说没有定义此函数，如何调用这个函数呢。

大佬您好，请问您解决这个问题了吗？我也想问一下怎么调用这个函数，感恩！

大佬您好，为啥 那请问局部索引应该怎么显示呢？谢谢！

多谢！已更新

@李东岳 感谢李老师的回复，谢谢

@李东岳 好的，谢谢东岳老师！

@李东岳 岳 东岳老师好，上面公式S其实表达的意思不一样的，我这里用错了，现在存在的困难是怎么求解这个压力方程值p'，感觉好难写出来只关于p'的方程，再进行求解，希望能得到东岳老师和各位前辈的指点，谢谢！
111.png

你这个最好用wedge网格来计算就不会出现这个问题了