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

@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，网格不算太多

兄弟 研究的咋样啦