interFoam短周期波浪衰减过于严重
-
我目前正在做波浪模拟相关的模拟，在模拟中我发现OpenFOAM v2106在使用interFoam模拟波浪时，短周期波浪的波高衰减非常严重，在离造波边界约5倍波长的位置波高衰减了至少一半，有人遇到过这种情况吗？
模拟的波浪为StokesI，也就是线性波，但是会有变形；水深在0.68 m和0.8 m之间，波高为0.08 m到0.16 m；严重衰减的波浪周期在1.0 s到1.5 s之间。湍流模型为k-Omega SST（包括可变密度的版本），关闭动量预测，时间格式为Crank-Nicolson （value = 0.9）。造波方法为速度边界造波，对应U的边界为waveVelocity，alpha.water的边界为waveAlpha，在另一侧边界设置消波。时间步长控制在1/1000周期一下，小范围波动不影响结果；x方向网格控制在1/100波长以下，y方向静水位附近网格控制在1/16波高以下。
现在的问题是，造波边界造出来的波形基本上对的上，但是短周期的波衰减过快。当波周期增大时，这种衰减就会相对减少，在2 s左右就和实验对的差不多了。我知道OpenFOAM不太能处理波陡0.05以上的波浪，但通过降低波高减小波陡后并没有很显著的改善。我也看到过文献说OpenFOAM自带的湍流模型没有考虑交界面附近的密度变化导致湍流被高估，但OpenFOAM v2106已经植入了文献作者提供的变密度湍流模型，也没有明显改善。同时，加密网格、降低时间步长、更换离散格式和增加nCorrectors的数量都不适用。
有没有大神遇到过这种情况？如果有，是怎么解决的？
-
@tidedrinker
纯探讨，如不对，请原谅：
是不是破碎了 ？
有文献说，那些理论临界值，自然界有时都很难达到
。。。
模拟上就更难了 -
@tidedrinker
不是经验，但真遇到过：- 破碎比物理实验和理论临界来的”早“
- 一旦破碎，波高会降很多
-
https://github.com/BrechtDevolder/buoyancyModifiedTurbulenceModels
我之前造波也遇到这个问题，层流造波没有衰减，湍流衰减很大，用这个解决了 -
-
@tidedrinker
我没验，你验证下，看到底用哪个理论。。。 深h对T影响有限，我更倾向深水用stokes1还有，请教下：
你熟悉不规则波浪反射的估算吗？ -
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，但依旧迭代不出来，请问您有什么解决办法吗
