我最近看了下，

替代文字

上面这个方程应该是稳态算法的方程。里面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网格来计算就不会出现这个问题了

@王慧博 在 矩阵求解器 中说：

LUsolve

您好。请问LUsolve是OF自带的函数吗？ 还有一个问题，为啥LUsolve可以直接调用而不是通过对象访问成员函数的方式调用呢？难道LUsolve是一个全局函数吗？谢谢大神！！！

主要输OF中的残差是线性方程组求解的残差，不是时间步迭代的残差，您说对吗？李老师。另外，是不是只有稳态问题我们才会关注线性方程组求解的残差残差，瞬态问题的residual 有意义吗？

已解决。。。先给出ADD的表达式，然后对其进行限制。修改代码为：
屏幕截图 2021-04-07 112950.png

@李东岳 WALE不需要vanDriest吧？

@huntersong 哇，好详细，我去试试， 十分感谢！

@wangfei9088 哇，太谢谢了。我也考虑过是单位的问题，但还未在代码中找到数据验证猜想，非常感谢大佬！

在dargmodel.C这个衍生类中没有mesh这个成员对象，不知道该如何引入

@xpqiu 感谢您的答疑
原来如此，是两个box相当于框定了一个水域区范围间的意思。
zoneID场的值设置明白了，还想问一下zoneID这个场是不是因为重叠网格所以要设置的，具体作用是什么呢？

我大部分情况下模拟结果都是这种
上面呢种是因为激波前移叠加了？
不太懂啊
Snipaste_2021-04-02_21-04-11.png

@bestucan 好的，谢谢老师的耐心解答！

@sr0922msc 我没有尝试过运行4.0之前的版本，即使是4.0版本，我在运行tutorials案例的时候，也遇到了一些问题，有些可以运行，有些不可以

因为这不是C:mihu:
这是C++。
你在这里修改，但是其他过程在其他地方修改了。而且这个赋值八成也是不成功的。给类成员赋值一般都由成员函数完成。

改起来还是挺麻烦的。:chitang: