@东岳 谢谢

好像发现原因了，Newtonian没有从singlePhaseTransportMode继承，因此如果在singlePhaseTransportMode中定义虚函数是不行的

@Samuel-Tu 原来这样，简单的重构方法好像确实有一个代价就是网格要比较细，不然不准确。

加油加油

是的，都是绝对坐标系的

最近在家呆的闲，又没有实验室的资源，，想看看磁流体的流固耦合的代码和书最好中文的，，，求推荐啊～

滑移壁面没有剪切力

@bubbleEL 您好，请问您德高斯分布函数已经实现了吗？最近也在实现这种办法，不知道能否指点一下？另外一直在思考超出网格之后发散的数学原理是什么？不知道您这边是不是已经弄清楚了？

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

@程迪 压力都是零梯度边界条件，压力泊松方程的系数矩阵是奇异的，adjustphi让矩阵不奇异了吗？奇异矩阵怎么解的？

可以用吗？如何解决，求分享

@东岳 明白了！谢谢老师

@Exthan 不知道是边界条件设错了还是应该更换求解器

不是的，后面我用komegaSST算了一下这个模型，吧压力求解器改成了PCG之后就收敛了，而且结果比较满意，希望对遇到同样问题的朋友也有参考价值

大家好，我最近在模拟波浪。在现在interfoam底下的case里，对于波的设置条件里，有一项ramptime，我不清楚它的意思，请大家帮忙解答。谢谢🙏

入口边界条件？我有一个问题，最近看到一些模拟射流的把入口延长了一段，加了管内流动，用来避免喷管直接出口处的较大的梯度，东岳老师怎么看@东岳 ，不知道您的入口速度怎么给的

@东岳 我是按这个公式计算的，密度为1，那么u_tau=sqrt(tau_w)，但是我发现，论文里的摩擦速度刚好是我用OF算出来的1.225倍。论文是用fluent算的，fluent里面的空气密度是1.225。我现在想不通 为什么会差这1.225倍。:jingya:

@东岳 感谢🙏🏼 试了一下还是不行，有新的其他函数的错误了😹

回复: 请教个关于编译求解器时报错的问题solidificationFoam.zip :high:

@马乔 在 OF文件读写 中说：

#include "fvCFD.H" // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // class proxyWriter : public regIOobject { private: scalar value_; proxyWriter(const proxyWriter&) = delete; void operator=(const proxyWriter&) = delete; public: proxyWriter(const IOobject& io, const scalar& value) : regIOobject(io), value_(value) { } virtual bool writeData(Ostream& os) const { os<<value_; return os.good(); } proxyWriter& operator++() { ++value_; return *this; } }; int main(int argc, char *argv[]) { #include "setRootCaseLists.H" #include "createTime.H" // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // proxyWriter writer ( IOobject ( "value", runTime.timeName(), runTime, IOobject::NO_READ, IOobject::AUTO_WRITE ), 0. ); Info<< "\nStarting time loop\n" << endl; while (runTime.loop()) { Info<< "Time = " << runTime.timeName() << nl << endl; ++writer; runTime.write(); runTime.printExecutionTime(Info); } Info<< "End\n" << endl; return 0; } :chigua:

@马乔 @Samuel-Tu 请问这个怎么用？这是重新编译一个lib然后link到controlDict吗？那怎么知道是在这一步的什么时候输出？是计算前还是计算后？