跟进一下，这是不同边界条件下的压力云图。可以直观的看到入口处压力场出现震荡。
c908d44b-b2db-44fb-91fb-276fddda4220-image.png

我也尝试从数值角度说一下我的看法。计算求解的是NS方程，对于入口处速度和压力均为定值的初始条件，可能会导致结果不满足NS方程。造成所谓的过度约束。而入口出口分别指定速度和压力为固定值，则会避免过度约束的现象发生。

欢迎各位大佬批评指正！

chemkin原来使用的输入文件中的传输数据chemkinToFOAM识别不了，需要将组分的传递数据用openFOAM的格式自己写一个transportProperties

@YuPeng 好的，感谢！

@Joann 在 什么原因使得“on signal 9 (Killed)” 中说：

大家好，计算后发现求解器的log在算第一个时间的时候就出现了这样的情况：

Primary job terminated normally, but 1 process returned a non-zero exit code. Per user-direction, the job has been aborted. -------------------------------------------------------------------------- [joann:06935] Read -1, expected 5328, errno = 3 [joann:06936] Read -1, expected 5400, errno = 3 [joann:06926] Read -1, expected 5400, errno = 3 [joann:06931] Read -1, expected 10368, errno = 3 [joann:06932] Read -1, expected 5616, errno = 3 [joann:06933] Read -1, expected 6984, errno = 3 -------------------------------------------------------------------------- mpirun noticed that process rank 0 with PID 0 on node joann exited on signal 9 (Killed).

对于这个问题有点迷茫，特来请教各位一下如何解决这个问题？
你好，请问您解决了吗？是什么原因呢

@Hope 18230146705@163.com，谢谢师兄，祝师兄科研顺利！

如果是偏应用的东西，一般来说也不需要人指导，自己跑几个算例就行，虽然学习周期要更长一些（相对商软）

之前我讲课的时候，OpenFOAM应用还是重头戏。http://dyfluid.com/class.html 这个下面有课程演变记录。另外自己每年CFD功力都在增长，7 8年前的自己跟现在完全不能比。2025年，OpenFOAM资料天花乱坠。基本所有的流体课题组都在做CFD。OpenFOAM应用也顺势的要更加简单（因为资料多了，玩的人也多了，交流也多了）。 OpenFOAM应用也不会作为我以后讲课的重头戏（今年只放了2个下午）。

因此类似咨询OpenFOAM应用的问题，我也没有详细回答。不过总体来说，现在做OpenFOAM应用的门槛应该比20年前要小很多。学习周期也要比之前缩短。但还是要比fluent长，毕竟一个点点点，一个要自己改文本文件。

另外关于OpenFOAM应用的问题，我本想直接放弃不讲，直接上算法编程。但是90%的人真是0基础。不讲直接要废。所以今年就放了2个下午。也只有2个下午，如果会还可以不听。

现在做算法多了，已经对OpenFOAM应用没有什么大的兴趣，之前写的洋洋洒洒的这种文章，现在已经写不出来了 OpenFOAM教程 初级入门建议

李老师，感谢您的修改，对比后我发现我对fvsolutions的设置还不太明确。用您改后的案例，openfoamv9跑出来了，但用openfoam-v2206还是报原来的错误。目前，我准备尝试用clion调试找一下壁面函数中wallDist的问题

@Joe 我也是小白还在折腾这个三相的求解器😂，OF似乎没有现成的三相空化求解器，所以想着只能在两相模型的基础上推广到三相然后自己改一个求解器出来了。

这个案例中，看stl文件，模型并没有闭合，为啥可以mode inside。还有是不是进行STL文件进行网格划分时需要一个整体的STL文件?但是看这个案例中并没有给出。

直接mapFields就可以吧？我用过，是可以把自适应加密之后的网格，映射到没有加密的规整网格上的。
如果是说机理上面可不可以，那其实是损失了加密之后那些网格的信息。你用另外一个计算模型也是要自适应网格吗？那可以直接把自适应加密之后的网格也拷贝过去。

@李东岳 在 interfoam计算实际河道 中说：

我也不知道了。类似这种我都是让甲方提供 ..

谢谢东哥，我慢慢找吧，实在弄不出来可能就换东西做了

1000w规模

想debug看看，就是跑不动啊..

最近在搞反应ODE看了一下代码，做个记录。

在chemistryModel.C中有定义：

forAll(RR_, fieldI) { RR_.set ( fieldI, new DimensionedField<scalar, volMesh> ( IOobject ( "RR." + Y_[fieldI].name(), mesh.time().timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE ), mesh, dimensionedScalar("zero", dimMass/dimVolume/dimTime, 0.0) ) ); }

RR并不是体场。所以在求解器层面，可以定义一个volScalarField，然后把RR的值复制给这个volScalarField然后进行输出，并在paraview查看。

也可以吧上面的NO_WRITE改成AUTO_WRITE，就是不能再paraview查看。

@Amadeus 谢谢您 我这边换了个喷油的方式 可以算了 等手头事情做完了我试试您说的方法！

@HITSC30 感谢回答，两篇文章写得很详细！

33095a1e-6724-4493-b8a0-376e8d4e19a6-图片.png
请问各位老师，这是我做的颗粒对壁面的磨损，但是从paraview中看磨损在碰到壁面时为贯穿式，其中壁面通过SW建模，fluent meshing画网格。有没有解决方法，十分感谢。下面是liggghts代码：
221c09f4-6c30-4abb-9817-98cea3a76378-image.png
e4ce9c6d-8813-406b-a6ba-36e600ff8f26-image.png

好的，感谢李博

你的这个算例是并行计算吗? 如果是并行，那么建议把网格划分方式先改成

numberOfSubdomains 6; method simple; simpleCoeffs { n (3 2 1); delta 0.001; }

这种比较少的核心数跑一次.

如果串行也会出现这种情况就确实很麻烦....

66fe34e5-2b9e-4952-ab9b-9c798f0e4be5-1740230382354.png

#include "kinematicMomentumTransportModels.H" #include "transportModel.H" #include "addToRunTimeSelectionTable.H" #include "RASModel.H" // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // makeBaseMomentumTransportModel ( geometricOneField, geometricOneField, incompressibleMomentumTransportModel, IncompressibleMomentumTransportModel, transportModel ); #define makeRASModel(Type) makeTemplatedTurbulenceModel (transportModelIncompressibleTurbulenceModel, RAS, Type) #include "myRAS" makeRASModel(myRAS); 以上是models.C 单独编译myRAS是没有报错的，编入的时候出错了

@李东岳 是的李老师。此外我还发现一个问题，我通过paraview中，切片+fliters+Alphabetical+plot on short lines的方式得到的结构物表面的压力曲线，除以（1/2U^2）之后得到的Cpn曲线，发现与其他论文中发布的结果之间幅值仍然差一个密度的结果，但是在postprocessing中读取的升力结果是匹配的，差异很小。想咨询是什么原因导致的，是读取表面压力的时候读取错了吗