@yuan您好！我觉得是坐标值而非绝对值。

@李东岳 在 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中读取的升力结果是匹配的，差异很小。想咨询是什么原因导致的，是读取表面压力的时候读取错了吗

@Amadeus 在 preCICE软件在CentOS7超算集群上无root权限从源码编译安装记录贴 中说：

但是如果碰到了如下和_MT相关的错误

此处错误的原因在于spooles库编译时没有开启多线程(Multi-Thread)，如果你的calculix-preCICE不需要启用多线程，那么就无需往下看了.

spooles多线程编译

@Amadeus 在 preCICE软件在CentOS7超算集群上无root权限从源码编译安装记录贴 中说：

spooles可以按照指南编译

按照preCICE官方给出的指南是默认没有打开spooles库的多线程编译选项的.
因此在make lib之后，需要进入 spooles库所在路径/spooles.2.2/MT/src 中，然后执行

make makeLib

注意大小写敏感. 然后就会把多线程支持库合并编译进之前生成的spooles.a文件.

最后回到calculix-adapter文件夹，make clean清理之前生成的文件，把makefile文件第55行的CFLAGS改回来即可

CFLAGS = -Wall -O3 -fopenmp $(INCLUDES) -DARCH="Linux" -DSPOOLES -DARPACK -DMATRIXSTORAGE -DUSE_MT

最后 make -j 重新编译

@ShaneHEEE 设置的是从x方向流入 我再检查一下

压力速度耦合那个我当时咋处理的有点忘了 :135: ，可能是松弛因子，要不就是换成SIMPLE或者SIMPLEC

不过fixedMean看起来更加符合物理，我在讲课的时候都让他们改fixedMean

http://dyfluid.com/openfoam.html#id7 看看这个对不对

TRPjet是啥，估计应该是fluid:592
solid2
solid1 之一

@李东岳 有的 我这边一个算例大概要4-6h，网格不算太多

兄弟 研究的咋样啦

解决了，在网格合并命令mergeMeshes之前执行setWaveParameters和setWaveFields命令就行了

可以考虑试试Nektar++