CFD中文网

我正在模拟三维下的波浪-海流-结构物相互作用，主要用olaFlow，偶尔也用基金会版本的interFoam。然后发现olaFlow的造波边界也好，OF6自带的造波边界也罢，在多面体网格中表现的都不太行，olaFlow的边界尤甚。而且我在网上搜索到的算例，要么是结构网格，要么是笛卡尔网格（切体网格），没怎么看到多面体或者四面体网格。

所以想问一下，是不是OF框架下的边界造波，对多面体网格的支持都不是太好？或者说，相对于横平竖直的结构网格和迪卡尔网格，在多面体网格上应用造波边界需要什么特殊的设置？

问题在于，点在网格单元中的位置不确定，考虑相邻网格间可能有较大速度差异，如果点不在单元正中则单元对应的速度不能代表某个具体点的速度。

这个不太好处理，有限体积法里面REA的average就是对每个网格的速度做了平均出一个体心速度。如果任意空间位置都要做插值，这是高精度有限体积法。

@tidedrinker :146: :146: :146: :146:

Hi https://cfd-china.com/user/tidedrinker

我在使用of2306也出现类似的情况，我的波频是1.3 而且不同网格密度的相位差和衰减率有明显差异

但是我在使用最密网格（单位波高内20个网格）时，总是卡在67 迭代不出来 最终计算发散

Courant Number mean: 1.18882076445e-06 max: 7.50119780436 Interface Courant Number mean: 3.00160678714e-08 max: 3.95428576669 deltaT = 6.18380398924e-152 Time = 67.132 PIMPLE: iteration 1 Updating StokesI wave model for patch inlet smoothSolver: Solving for alpha.water, Initial residual = 5.38442986655e-09, Final residual = 5.38442986655e-09, No Iterations 0 Phase-1 volume fraction = 0.666618301469 Min(alpha.water) = -2.29205431887e-10 Max(alpha.water) = 1.00000001308 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666618301469 Min(alpha.water) = -2.29205431887e-10 Max(alpha.water) = 1.00000001308 Updating shallowWaterAbsorption wave model for patch outlet DICPBiCGStab: Solving for p_rgh, Initial residual = 0.844109801632, Final residual = 0.00558442869736, No Iterations 20 DICPBiCGStab: Solving for p_rgh, Initial residual = 0.742328538632, Final residual = 0.00717435702295, No Iterations 26 time step continuity errors : sum local = 1.13821305653e-09, global = -5.40663033957e-10, cumulative = 2.0897431155e-05 DICPBiCGStab: Solving for p_rgh, Initial residual = 0.823953322827, Final residual = 0.00800680485496, No Iterations 25 DICPBiCGStab: Solving for p_rgh, Initial residual = 0.867377012998, Final residual = 9.38111948679e-10, No Iterations 389 time step continuity errors : sum local = 9.20755023008e-16, global = 7.16963646939e-16, cumulative = 2.08974311557e-05 smoothSolver: Solving for omega, Initial residual = 7.31631987399e-08, Final residual = 4.57762864183e-11, No Iterations 1 smoothSolver: Solving for k, Initial residual = 4.39226774226e-07, Final residual = 2.03583449134e-10, No Iterations 1 ExecutionTime = 29318.96 s ClockTime = 29376 s Courant Number mean: 1.18884127695e-06 max: 7.50169136387 Interface Courant Number mean: 3.0019195664e-08 max: 3.95459850554 deltaT = 8.24321301597e-153 Time = 67.132 PIMPLE: iteration 1 Updating StokesI wave model for patch inlet smoothSolver: Solving for alpha.water, Initial residual = 5.38393271075e-09, Final residual = 5.38393271075e-09, No Iterations 0 Phase-1 volume fraction = 0.666618301469 Min(alpha.water) = -2.29205431887e-10 Max(alpha.water) = 1.00000001308 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666618301469 Min(alpha.water) = -2.29205431887e-10 Max(alpha.water) = 1.00000001308 Updating shallowWaterAbsorption wave model for patch outlet DICPBiCGStab: Solving for p_rgh, Initial residual = 0.844175076514, Final residual = 0.00539645720402, No Iterations 20 DICPBiCGStab: Solving for p_rgh, Initial residual = 0.742024848316, Final residual = 0.00733560856723, No Iterations 27 time step continuity errors : sum local = 1.16277362766e-09, global = -5.78716682763e-10, cumulative = 2.08968524391e-05 DICPBiCGStab: Solving for p_rgh, Initial residual = 0.824114713535, Final residual = 0.00803020015335, No Iterations 25 DICPBiCGStab: Solving for p_rgh, Initial residual = 0.867366339743, Final residual = 9.11722968482e-10, No Iterations 373 time step continuity errors : sum local = 8.95106772158e-16, global = -5.97636368108e-16, cumulative = 2.08968524385e-05 smoothSolver: Solving for omega, Initial residual = 7.30573956466e-08, Final residual = 4.54407875238e-11, No Iterations 1 smoothSolver: Solving for k, Initial residual = 4.38607733913e-07, Final residual = 2.03320386843e-10, No Iterations 1 ExecutionTime = 29358.38 s ClockTime = 29416 s Courant Number mean: 1.18886725019e-06 max: 7.50212848379 Interface Courant Number mean: 3.00241095746e-08 max: 3.95487860459 deltaT = 1.09878323649e-153 Time = 67.132 PIMPLE: iteration 1 Updating StokesI wave model for patch inlet smoothSolver: Solving for alpha.water, Initial residual = 5.38400281089e-09, Final residual = 5.38400281089e-09, No Iterations 0 Phase-1 volume fraction = 0.666618301469 Min(alpha.water) = -2.29205431887e-10 Max(alpha.water) = 1.00000001308 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666618301469 Min(alpha.water) = -2.29205431887e-10 Max(alpha.water) = 1.00000001308 Updating shallowWaterAbsorption wave model for patch outlet -------------------------------------------------------------------------- Primary job terminated normally, but 1 process returned a non-zero exit code. Per user-direction, the job has been aborted. -------------------------------------------------------------------------- -------------------------------------------------------------------------- mpirun detected that one or more processes exited with non-zero status, thus causing the job to be terminated. The first process to do so was: Process name: [[8254,1],7] Exit code: 145 --------------------------------------------------------------------------

在粗网格和中等网格却计算很稳定。根据https://cfd-china.com/topic/5462/mpirun-detected-that-one-or-more-processes-exited-with-non-zero-status，尝试修改了fvSolution，但依旧迭代不出来，请问您有什么解决办法吗

收藏一下，感觉有用

我最近出现了这个情况。还有人有这个情况么，很奇怪

@tidedrinker 前辈你好，我在使用olaFlow做一个二维的狭缝共振，但是把两个浮箱加上之后造出的波浪不向前传递是什么原因呀，详细的我放在了这个问题里面问题地址

@tidedrinker 在 functions里的forces在二相流里能否输出基于p_rgh的压力和转矩？ 中说：

@李东岳 这个forces是写在controlDict/functions里的，在求解过程中每个（或每若干个）时间步输出一次。这个压力和转矩是随时间变化的，不是根据结果文件用postProcessor -func那一套程序得到的。所以修改结果文件的文件名应该没什么用处？

请问一下forces在计算完成后，如果想要得到其他面上的forces，有什么办法后处理提取吗？看见您说不是根据结果文件用postProcessor -func那一套程序得到的,请教一下，这有什么比较好的处理办法吗？

@tidedrinker 谢谢大佬！

@bestucan 我决定直接放图上去了，不用GCI了，主要是因为GCI到底应该限制在多少以下没有定论，同时GCI也需要求解结果满足单调性，偏偏我这个满足不了……

确实有很大改善，可能上方空气部分需要设置一个比较大的高度

@bestucan 头文件里说的是如果非保守通量（ non-conservative fluxes，应该是翻译成非保守通量吧？）对算法的预测部分有影响的话是需要correctPhi的，但是什么情况下会影响，我完全没有概念。至于后半说多相流的部分，4.1多相流里不开correctPhi的还是有几个的。所以我完全没提取出什么有效信息。

@tidedrinker 感谢分享！:xiexie:

@Hungryandfool 您说的我也想过，但我主要是应付不可预测的爆发性计算需求，超算中心这种流程繁琐，而且灵活性不足。

@悬铃神木 谢谢。

我在用interFoam时，笛卡尔网格可以得到合理的计算结果。但换成非结构网格（prism之类）后，计算虽没有发散，但界面有问题，请问有什么需要改进的地方吗？@东岳 @tidedrinker

关于OpenFOAM中在采用K-Omega SST模型时，使用/不使用壁面函数的Y+及边界设置问题，可以参考这个链接链接文本，解释的很详细。

@tidedrinker 0_1541396115574_8a12b935-4737-407d-8a75-79e3bcf4fe98-image.png

@李东岳 我最后还是把这个打开了。我模拟的是潮流能水轮机，开之前功率系数不过30%出头，开了之后到了45%左右（TSR=6），功率系数-叶尖速比曲线基本上与实验对上了。而且开了之后对时间步长就不敏感了，时间步从0.0001s到0.002s，CoN从40多一直到700多，功率系数变化在5%以内。

而且，OpenFOAM/tutorials/multiphase/interPhaseChangeDyMFoam/propeller这个算例也是开启了动量预测，我考虑是不是固定转速的旋转机械都需要开这个选项。

@李东岳
工程上有意义的是达到指定误差的计算代价，一般理论上认为误差可以从某种定义的残差推算出的不确定来衡量，时间以taubench用25000网格迭代10步作为1个时间单位来衡量。

这样的假设下做横向比较是可能的。但是，具体情况嘛，呵呵。湍流模型的不确定度就说不清楚。