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++

@ansel 不需要让VMware找到wsl啊，可以把wsl理解成一个没有GUI的虚拟机，它可以用物理机的全部CPU资源(虽然读写性能会弱一点)，但是VMware是一个有GUI的虚拟机，它只能用你允许范围内的CPU/内存资源. 这两者是相互独立且互不影响的....

目前想要用多孔介质代替小孔计算浮动owc的重叠网格模型

@Aaron 有回流应该是正常现象，跟毛细数和粘度比有关吧，参考这篇Advances in Water Resources 170 (2022) 104324

@RQLI_CFD 在 请教：第三代涡识别方法Liutex的代码在linux中运行的问题 中说：

@lizhisongsjtu 感觉这个liutex没有那么好啊，我计算的是槽道流，liutex识别的结果很混乱（liutex的等值面），感觉还没Q准则显示的涡结构清晰，但是看刘老师文章里的效果挺好的。。。

您好，我也遇到这个问题了，Q准则识别的槽道流涡结构比Liutex好太多了，Liutex识别出来的涡都是断断续续的，怎么调阈值都不行。请问您解决了吗？

这个进口水面的高度可以自行调节，请参考无痛笔记，variableHeightFlowRate边界条件的讨论

不用啊，按照大佬2楼说的那么写就能用