CFD中文网

@bestucan 谢谢老师的手把手指导！原来是可以沿着基类往上找的，之前一看到继承了好几个基类就感觉太麻烦找不出来了。

使用grep在src中搜还可以搜到用法，也很不错，但是函数定义换行了就不好搜了。比如这个要是用 grep -rP "::average\n\(" . 啥也找不到。

编译本地doxygen也学到了，openfoam2406版本的方法整理一下

sudo apt-get install doxygen graphviz cd /usr/share/doc/openfoam2406-common/doc/Doxygen/ gzip -dk Doxyfile.gz #解压doxygen配置文件 ./Allwmake firefox html/index.html

本地查看foam下的函数不会卡了，不过网页中没有的函数原型本地也没有

十分感谢您的分享！

@Mina-Lee 在 大步长下变形网格发散 中说：

换句话说，是不是overset边界，根本不影响变形。fixed边界，才会影响变形。

确实是啊！displacementLaplacianFvMotionSolver类中组建cellDisplacement场时，只对pointDisplacement为fixedValue类型的patch进行了处理（改为cellMotion类型，并对face各点求平均），其它类型的没做任何处理；最后cellDisplacement再插值回pointDisplacement时，点位移场只有原本的fixedValue边界没被修改，其它的都被求解结果修改了。

我看代码时没深想到这一层，只在用的时候发现有时远场使用fixedValue和slip效果不同，我在这个贴子中的回复属于是瞎猫碰上死老鼠了，看来远场只要不是fixedValue，给任意都行。

最后我的问题还是用的单一区域网格以及求扩散方程变形成功实现的，为了验证还用StarCCM算了比过，吻合得很好。倒是在StarCCM中试过重叠网格，overset交界面网格难匹配，常常挖洞出错，算出来的结果锯齿也很厉害。至于rbf，有些懒得折腾了，以后遇到非用不可的情况再说。

这句话应该是您的“泣血之言”，指给后来的人看看:146: :146:

@Mina-Lee 在 大步长下变形网格发散 中说：

如果你要用重叠网格，记住，不要用逆距离，OF的这个方法存在BUG！

@李东岳 我用的v2006，Phi是速度势，用-grad(Phi)计算的速度和输出的U有些差异，尤其是在远场边界。

@李东岳 感谢老师回复！动量预测是关掉了的。

打开试了试，步长极小，残差也一直保持在一个值不下降，cumulative则很大，如下：

Courant Number mean: 7.953414508e-05 max: 0.4993036972 deltaT = 4.979675437e-35 Time = 6.00001e-05 PIMPLE: iteration 1 DICPCG: Solving for pcorr, Initial residual = 1, Final residual = 0.04941709914, No Iterations 29 DICPCG: Solving for pcorr, Initial residual = 0.03122316034, Final residual = 0.001556737496, No Iterations 121 DICPCG: Solving for pcorr, Initial residual = 0.006804280077, Final residual = 0.0003314734696, No Iterations 359 DICPCG: Solving for pcorr, Initial residual = 0.002308369425, Final residual = 9.709852492e-11, No Iterations 911 time step continuity errors : sum local = 1.793718362e-20, global = 1.178518414e-23, cumulative = 1.037419916e+19 smoothSolver: Solving for Ux, Initial residual = 3.985498772e-05, Final residual = 2.227155701e-11, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 6.302253167e-05, Final residual = 2.385619185e-11, No Iterations 1 GAMG: Solving for p, Initial residual = 0.004038713667, Final residual = 0.0001722133227, No Iterations 18 GAMG: Solving for p, Initial residual = 0.0005583071565, Final residual = 2.673668049e-05, No Iterations 34 GAMG: Solving for p, Initial residual = 9.485975434e-05, Final residual = 4.53317788e-06, No Iterations 33 GAMG: Solving for p, Initial residual = 2.596708738e-05, Final residual = 1.231466526e-06, No Iterations 38 time step continuity errors : sum local = 7.929103562e-15, global = -8.302368146e-16, cumulative = 1.037419916e+19 GAMG: Solving for p, Initial residual = 0.003717218882, Final residual = 0.000164557501, No Iterations 17 GAMG: Solving for p, Initial residual = 0.0005378364039, Final residual = 2.53035408e-05, No Iterations 34 GAMG: Solving for p, Initial residual = 9.265726415e-05, Final residual = 4.535601234e-06, No Iterations 36 GAMG: Solving for p, Initial residual = 2.556961887e-05, Final residual = 9.918214454e-11, No Iterations 548 time step continuity errors : sum local = 6.482264028e-19, global = 6.603865839e-20, cumulative = 1.037419916e+19 PIMPLE: iteration 2 ... PIMPLE: iteration 50 smoothSolver: Solving for Ux, Initial residual = 1.226945138e-10, Final residual = 6.6691269e-16, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 1.651845516e-10, Final residual = 3.177048458e-16, No Iterations 1 GAMG: Solving for p, Initial residual = 0.000949721577, Final residual = 4.516299707e-05, No Iterations 41 GAMG: Solving for p, Initial residual = 0.0004457930337, Final residual = 2.206483819e-05, No Iterations 23 GAMG: Solving for p, Initial residual = 7.977007558e-05, Final residual = 3.747481321e-06, No Iterations 36 GAMG: Solving for p, Initial residual = 2.418915618e-05, Final residual = 1.162060008e-06, No Iterations 48 time step continuity errors : sum local = 7.482194647e-15, global = -7.221446192e-16, cumulative = 1.037419916e+19 GAMG: Solving for p, Initial residual = 0.003705959955, Final residual = 0.0001646823634, No Iterations 17 GAMG: Solving for p, Initial residual = 0.0005377574684, Final residual = 2.535476181e-05, No Iterations 34 GAMG: Solving for p, Initial residual = 9.241219827e-05, Final residual = 4.525895839e-06, No Iterations 36 GAMG: Solving for p, Initial residual = 2.542969083e-05, Final residual = 9.774400044e-11, No Iterations 548 time step continuity errors : sum local = 6.388632178e-19, global = 6.502735317e-20, cumulative = 1.037419916e+19 PIMPLE: not converged within 50 iterations ExecutionTime = 21923.2 s ClockTime = 21963 s Courant Number mean: 7.964991508e-05 max: 0.4989139005 deltaT = 4.99051583e-35 Time = 6.00001e-05

看到站内其它的帖子，试了试有量纲/无量纲计算（此前圆柱直径0.05，nu=50；改为圆柱直径0.0005，nu=0.5），也无法收敛。

走投无路试了一下StarCCM，多低雷诺数都能算，和OpenFOAM那些能收敛的算例比得还贼好