前面2e-06,4e-06的数据都能reconstruct,到6e-06时就发生了这个错误，显示找不到/home/mywsl/OpenFOAM/mywsl-6/run/2D_detonation4/processor0/6e-06/polyMesh/pointProcAddressing，但我看了，pointProcAddressing文件存在，有大佬知道这个问题吗

@zzkluck 感谢回复，已经用上windows的paraview了

学到了学到了

@李东岳 多谢李老师鼓励，还得感谢李老师的讨论答疑和提供的交流平台:ok: 。

@李东岳 李老师，我在尝试使用openfoam10当中的multiphaseEuler求解器来模拟不同含气率入口的气液两相离心泵，我想请问一下，这个求解器对网格的要求是不是比较高？我目前使用的是非结构化网格，之前计算过单相流的瞬态和稳态都没有问题。但是在目前计算气液两相时，时间步长会不断减小，从最开始的e-6减小到负十几次，应该是发散了吧？我是不是要使用结构化网格呢？麻烦老师给点意见，谢谢！

如题，数据导入进来是Point Data，尝试用Rotational Extrusiton把Point Data旋转成柱体，但是转出来的是一堆线，渲染时选Surface的结果类似这样：
fig_rotational.png
Screenshot_20241121_234954.png
Representation选Volume会报

ERROR: In vtkOpenGLProjectedTetrahedraMapper.cxx, line 389 vtkOpenGLProjectedTetrahedraMapper (0x55abb08f3810): Encountered non-tetrahedra cell!

同时，也无法画出Contour，无法Extract Surface，但是可以通过Clip筛数据画点云。
如果用Delaunay 2D，可以渲染出2D云图，但是2D Contour也不行。Delaunay 2D -> Rotational Extrusion旋转出来只有旋转起始面和结束面。Delaunay 3D渲染出的结果是空的。
尝试了vtkVolumeOfRevolutionFilter，结果是一样的。所以在此求教各位大佬，有没有办法把2D轴对称数据渲染成3D效果，画出Contour？

也可能需要模版参数特化的某个函数，指定的模版参数有问题，希望有用哈，我遇到的问题是，在自定义基类中把：
BasicTurbulenceModel::correctNut();
换成了TurbulenceModel::correctNut();
解决方法是传入
BasicTurbulenceModel和TurbulenceModel两个模版参数，
或者BasicTurbulenceModel::correctNut();只在子类中出现

@小刘lyw 哈哈哈我主要还是对颗粒数有疑惑，实在不明白应该具体算多少，有好几种公式，可以探讨一下，

遥遥领先

@李东岳 谢谢东子老师

@zzkluck debug大师

@李东岳 模型缺陷吗.....那就有点麻烦了。我不单纯写文章，是老师的项目。我再考虑考虑

@zzkluck :146: :146: :146:

@李东岳 老师我OF-4.1版本算例constant文件夹下得有turbulenceProperties字典，里面写的是laminar，不涉及湍流。好像OF-8以后都不需要turbulenceProperties字典了。

按照文档里的说法，有时换一下meshShrinker也许会有效果，感觉文档里那个例子和你的也有点像。

addLayersControls { ... meshShrinker displacementMotionSolver; solver displacementLaplacian; displacementLaplacianCoeffs { diffusivity quadratic inverseDistance 1(wall); } ... }

相应的还需要调整下fvSolution和fvScheme，具体文档里都有。

话说新版文档虽然变好看了，但图居然是糊的，可能还需要参考下旧版文档。

另外这里也请教下各位老师，按我的理解snappyHexMesh做完snap之后，在layer这步应该是从几何表面出发，向外挤出来一块空间把边界层塞进去。怎么在楼主的这个例子中反而是向内把几何都挤变形了

@Mina-Lee 在 precice tourorials中的受力问题 中说：

因为你编译的openfoam-adapter里，已经对流固耦合交界面的force变量引用了，OF在运行完成后，adapter会引用force变量，precice运行完毕后，会覆盖这个force变量。
precice就是给openfoam施加了一个狄利克雷位移边界条件-理解了这句话就全明白了。

其他的能理解，但是这个“precice运行完毕后，会覆盖这个force变量”指的是
在controldict中定义了“forces”那么openfoam输出的forec是重新计算的，与流固耦合过程中交换的不是一个值的意思吗

调试的时候用ascii，跑算例用binary。如果网格几十万以上，binary会小非常多

你好

PIMPLE: iteration 40 **time step continuity errors : sum local = 20238.39271, global = -1326.825454, cumulative = -1286.292262** **time step continuity errors : sum local = 1965.022984, global = -493.7266286, cumulative = -1780.018891** PIMPLE: iteration 22 GAMG: Solving for p, Initial residual = 0.1839153372, Final residual = 1.265686901e-06, No Iterations 20 GAMG: Solving for p, Initial residual = 0.1398796637, Final residual = 1.397575026e-06, No Iterations 17 GAMG: Solving for p, Initial residual = 0.1463017811, Final residual = 1.437462198e-06, No Iterations 16 GAMG: Solving for p, Initial residual = 0.121320625, Final residual = 6.97291136e-08, No Iterations 20 PIMPLE: iteration 23 .GAMG: Solving for p, Initial residual = 0.1600126604, Final residual = 1.256444652e-06, No Iterations 23 GAMG: Solving for p, Initial residual = 0.09022545056, Final residual = 7.317342264e-07, No Iterations 20 GAMG: Solving for p, Initial residual = 0.08084307981, Final residual = 7.297894683e-07, No Iterations 18 GAMG: Solving for p, Initial residual = 0.06296048891, Final residual = 7.500379861e-08, No Iterations 22

你的问题在于：
1、瞬态项上，你的瞬态项很难收敛，残差爆掉了，这个似乎与你的瞬态项离散方式有关，可以试一下二阶向后差分。

ddtSchemes { default CrankNicolson 0.9; }

2、你的泊松方程收敛性很差，外循环很难收敛，残差很高，你定义的外循环收敛条件根本达不到，泊松方程在非正交修正后，迭代次数反而变多，证明你的网格质量可能较差。
建议：1.更换瞬态项离散方式。2.增加非正交修正次数，提高你的网格质量。
注意：外循环定义50次合理，一般难收敛问题，最多10次就可以收敛了。流固耦合交界面前几个时间步不收敛很正常，但是你这个是流体不收敛，这就不对了，建议你先单独跑一下流体看看，你这个流体设置不太好。
祝好~