@队长别开枪
内柱：Cylinder z1(0 0 -20) z2 (0 0 15) 半径 6.25
外环：z1(0 0 -7.25) z2(0 0 5.25) 内半径 7.5 外半径 12.5
外环原先是关于水面对称，下降一米让其上浮
麻烦大佬给看看
据另一位做六自由度双浮箱的同学说，他也遇到了类似的问题 只要overset之间相互接触就会发散，不知问题到底出在哪
H
hachikoi
@hachikoi
帖子
-
overset 小间隙-多浮体算例 计算结果明显有误 -
动网格和overset物体速度快速震荡导致崩溃的问题老哥们，你们是如何解决这个问题
我在使用两套overset做浮体自衰减的时候也遇到类似的问题 -
overset 小间隙-多浮体算例 计算结果明显有误大佬有没有什么建议呀 @队长别开枪
-
SHM尖锐直角边界层添加@xpqiu 感谢大佬回复，问题已解决 请问大佬有遇到过多重overset计算错误的问题吗
-
overset 小间隙-多浮体算例 计算结果明显有误背景网格的插值区，计算区和hole
中心柱overset的插值区，计算区和hole
外环柱overset的插值区，计算区和hole
-
overset 小间隙-多浮体算例 计算结果明显有误Rigid-body motion of the dynamicWECMiddle Centre of rotation: (0 0 -0.999946925817) Orientation: (1 0 0 0 1 0 0 0 1) Linear velocity: (0 0 0.00133640649968) Angular velocity: (0 0 0) inverseDistance : detected 3 mesh regions zone:0 nCells:3302728 voxels:(92 92 92) bb:(-1250.00076811 -1250.00076811 -800.000768115) (-749.999231885 -749.999231885 -499.999231885) zone:1 nCells:1329243 voxels:(92 92 92) bb:(-17.5000567663 -17.5000567663 -25.0007917663) (17.5000557663 17.5000547663 -1.56244523367) zone:2 nCells:1783736 voxels:(92 92 92) bb:(-4.41261149239 5.32858310761 -12.2503419182) (11.2436789924 17.5000299924 10.2504540666) Overset analysis : nCells : 6415707 calculated : 5795510 interpolated : 159835 (from local:606 mixed local/remote:100 remote:159129) hole : 455019 Updating shallowWaterAbsorption wave model for patch inlet Updating shallowWaterAbsorption wave model for patch outlet Updating shallowWaterAbsorption wave model for patch rightSide Updating shallowWaterAbsorption wave model for patch leftSide Execution time for mesh.update() = 41.86 s smoothSolver: Solving for alpha.water, Initial residual = 8.27297434572e-10, Final residual = 8.27297434572e-10, No Iterations 0 Phase-1 volume fraction = 0.666552338974 Min(alpha.water) = -0.000288674615659 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666552338827 Min(alpha.water) = -0.000288619975557 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 0.000659631763981, Final residual = 6.27878556034e-06, No Iterations 129 DILUPBiCGStab: Solving for p_rgh, Initial residual = 0.00981519206159, Final residual = 6.50074384625e-05, No Iterations 95 time step continuity errors : sum local = 5.76258113924e-11, global = 3.76596506058e-11, cumulative = 2.12788820726e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 0.000194046206327, Final residual = 1.08198900443e-06, No Iterations 137 DILUPBiCGStab: Solving for p_rgh, Initial residual = 5.43945734601e-05, Final residual = 8.53757949514e-10, No Iterations 262 time step continuity errors : sum local = 3.60665856703e-15, global = -2.65475508792e-15, cumulative = 2.1278855525e-09 smoothSolver: Solving for omega, Initial residual = 7.30361462036e-07, Final residual = 3.97781225689e-10, No Iterations 1 smoothSolver: Solving for k, Initial residual = 3.49789830215e-05, Final residual = 3.4609685557e-09, No Iterations 2 PIMPLE: iteration 2 smoothSolver: Solving for alpha.water, Initial residual = 8.16330553355e-10, Final residual = 8.16330553355e-10, No Iterations 0 Phase-1 volume fraction = 0.666552338827 Min(alpha.water) = -0.000288619975557 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666552338696 Min(alpha.water) = -0.00028861608681 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 3.64320126049e-05, Final residual = 3.53085923555e-07, No Iterations 139 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1.6562573215e-05, Final residual = 1.63348260399e-07, No Iterations 126 time step continuity errors : sum local = 6.65506742196e-13, global = 6.65113179978e-13, cumulative = 2.12855066568e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 7.05305597495e-06, Final residual = 4.86049614563e-08, No Iterations 143 DILUPBiCGStab: Solving for p_rgh, Initial residual = 3.0062961958e-06, Final residual = 3.82850298196e-10, No Iterations 202 time step continuity errors : sum local = 1.20177987263e-15, global = 2.18943311055e-16, cumulative = 2.12855088463e-09 smoothSolver: Solving for omega, Initial residual = 5.22965286878e-10, Final residual = 5.22965286878e-10, No Iterations 0 smoothSolver: Solving for k, Initial residual = 8.764653366e-09, Final residual = 8.764653366e-09, No Iterations 0 PIMPLE: iteration 3 smoothSolver: Solving for alpha.water, Initial residual = 8.34145135352e-10, Final residual = 8.34145135352e-10, No Iterations 0 Phase-1 volume fraction = 0.666552338696 Min(alpha.water) = -0.00028861608681 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666552338578 Min(alpha.water) = -0.000288612326428 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 2.26546928799e-06, Final residual = 2.26376259356e-08, No Iterations 133 DILUPBiCGStab: Solving for p_rgh, Initial residual = 7.73727028489e-07, Final residual = 7.58193596374e-09, No Iterations 142 time step continuity errors : sum local = 3.32815134752e-14, global = 3.32775001386e-14, cumulative = 2.12858416213e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 3.24877961418e-07, Final residual = 2.65040195196e-09, No Iterations 151 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1.36702609294e-07, Final residual = 9.73556700389e-10, No Iterations 142 time step continuity errors : sum local = 4.50942836308e-15, global = 3.24662375376e-15, cumulative = 2.12858740875e-09 smoothSolver: Solving for omega, Initial residual = 5.24948390428e-10, Final residual = 5.24948390428e-10, No Iterations 0 smoothSolver: Solving for k, Initial residual = 8.92548068881e-09, Final residual = 8.92548068881e-09, No Iterations 0 PIMPLE: iteration 4 smoothSolver: Solving for alpha.water, Initial residual = 8.39608653245e-10, Final residual = 8.39608653245e-10, No Iterations 0 Phase-1 volume fraction = 0.666552338578 Min(alpha.water) = -0.000288612326428 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666552338472 Min(alpha.water) = -0.000288608472331 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 2.67244877187e-07, Final residual = 2.53807729896e-09, No Iterations 156 DILUPBiCGStab: Solving for p_rgh, Initial residual = 8.94717131277e-08, Final residual = 8.39554535244e-10, No Iterations 127 time step continuity errors : sum local = 3.85428841404e-15, global = 2.11240953212e-15, cumulative = 2.12858952116e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1.8066547318e-08, Final residual = 6.21581241679e-10, No Iterations 95 DILUPBiCGStab: Solving for p_rgh, Initial residual = 6.6868460246e-09, Final residual = 9.95083309543e-10, No Iterations 14 time step continuity errors : sum local = 3.31155029848e-15, global = 2.63481371504e-16, cumulative = 2.12858978464e-09 smoothSolver: Solving for omega, Initial residual = 5.25158396946e-10, Final residual = 5.25158396946e-10, No Iterations 0 smoothSolver: Solving for k, Initial residual = 8.94269523662e-09, Final residual = 8.94269523662e-09, No Iterations 0 PIMPLE: iteration 5 smoothSolver: Solving for alpha.water, Initial residual = 8.49256382234e-10, Final residual = 8.49256382234e-10, No Iterations 0 Phase-1 volume fraction = 0.666552338472 Min(alpha.water) = -0.000288608472331 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666552338375 Min(alpha.water) = -0.000288604542788 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1.63410214157e-07, Final residual = 1.58076469785e-09, No Iterations 115 DILUPBiCGStab: Solving for p_rgh, Initial residual = 9.91057467543e-09, Final residual = 9.88772998935e-10, No Iterations 27 time step continuity errors : sum local = 4.32969109332e-15, global = 3.26463782778e-15, cumulative = 2.12859304928e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 6.98750307525e-09, Final residual = 9.88518391293e-10, No Iterations 14 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1.44234940454e-09, Final residual = 9.8025747263e-10, No Iterations 4 time step continuity errors : sum local = 3.90688904253e-15, global = 3.55096694681e-15, cumulative = 2.12859660025e-09 smoothSolver: Solving for omega, Initial residual = 5.25212396357e-10, Final residual = 5.25212396357e-10, No Iterations 0 smoothSolver: Solving for k, Initial residual = 8.94951618755e-09, Final residual = 8.94951618755e-09, No Iterations 0 ExecutionTime = 2032.89 s ClockTime = 2056 s Courant Number mean: 6.42988902653e-07 max: 0.749137867548 Interface Courant Number mean: 7.49720100137e-09 max: 0.0185446989986 deltaT = 0.000122067132065 Time = 0.0125595 PIMPLE: iteration 1 forces forces: rho: rho Not including porosity effects forces forces: rho: rho Not including porosity effects forces forces: rho: rho Not including porosity effects Rigid-body motion of the dynamicWECMiddle Centre of rotation: (0 0 -0.999947421749) Orientation: (1 0 0 0 1 0 0 0 1) Linear velocity: (0 0 -0.00946196445189) Angular velocity: (0 0 0) inverseDistance : detected 3 mesh regions zone:0 nCells:3302728 voxels:(92 92 92) bb:(-1250.00076811 -1250.00076811 -800.000768115) (-749.999231885 -749.999231885 -499.999231885) zone:1 nCells:1329243 voxels:(92 92 92) bb:(-17.5000567663 -17.5000567663 -25.0007917663) (17.5000557663 17.5000547663 -1.56244523367) zone:2 nCells:1783736 voxels:(92 92 92) bb:(-4.41261149239 5.32858310761 -12.2503424141) (11.2436789924 17.5000299924 10.2504535706) Overset analysis : nCells : 6415707 calculated : 5795510 interpolated : 159835 (from local:606 mixed local/remote:100 remote:159129) hole : 455019 Updating shallowWaterAbsorption wave model for patch inlet Updating shallowWaterAbsorption wave model for patch outlet Updating shallowWaterAbsorption wave model for patch rightSide Updating shallowWaterAbsorption wave model for patch leftSide Execution time for mesh.update() = 43.03 s smoothSolver: Solving for alpha.water, Initial residual = 5.62693014026e-10, Final residual = 5.62693014026e-10, No Iterations 0 Phase-1 volume fraction = 0.66655233832 Min(alpha.water) = -0.000288604542788 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666552338241 Min(alpha.water) = -0.000288600000297 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 0.00610895566778, Final residual = 1.76062946381e-05, No Iterations 187 DILUPBiCGStab: Solving for p_rgh, Initial residual = 0.000464626734755, Final residual = 4.31868001907e-06, No Iterations 135 time step continuity errors : sum local = 8.24477020959e-11, global = 6.57112660894e-11, cumulative = 2.19430786634e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 0.000453800900833, Final residual = 3.94151614699e-06, No Iterations 147 DILUPBiCGStab: Solving for p_rgh, Initial residual = 0.000147155197186, Final residual = 5.92856342839e-10, No Iterations 256 time step continuity errors : sum local = 7.18291390647e-15, global = -4.28128690922e-15, cumulative = 2.19430358505e-09 smoothSolver: Solving for omega, Initial residual = 8.2278501676e-07, Final residual = 3.51637527097e-10, No Iterations 1 smoothSolver: Solving for k, Initial residual = 3.84521991436e-05, Final residual = 9.69284044515e-09, No Iterations 1 PIMPLE: iteration 2 smoothSolver: Solving for alpha.water, Initial residual = 5.18632871574e-10, Final residual = 5.18632871574e-10, No Iterations 0 Phase-1 volume fraction = 0.666552338241 Min(alpha.water) = -0.000288600000297 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666552338169 Min(alpha.water) = -0.000288587367539 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 8.26591595287e-05, Final residual = 7.90344473439e-07, No Iterations 112 DILUPBiCGStab: Solving for p_rgh, Initial residual = 3.51930809664e-05, Final residual = 3.22380409956e-07, No Iterations 136 time step continuity errors : sum local = 3.9496037615e-12, global = 3.83624615083e-12, cumulative = 2.1981398312e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1.71819636959e-05, Final residual = 1.70602992744e-07, No Iterations 169 DILUPBiCGStab: Solving for p_rgh, Initial residual = 7.29985597309e-06, Final residual = 3.6900051247e-10, No Iterations 226 time step continuity errors : sum local = 4.54190924639e-15, global = -2.53540221376e-15, cumulative = 2.1981372958e-09 smoothSolver: Solving for omega, Initial residual = 5.7544716098e-10, Final residual = 5.7544716098e-10, No Iterations 0 smoothSolver: Solving for k, Initial residual = 2.66232952193e-08, Final residual = 3.81501056948e-09, No Iterations 1 PIMPLE: iteration 3 smoothSolver: Solving for alpha.water, Initial residual = 5.38627943033e-10, Final residual = 5.38627943033e-10, No Iterations 0 Phase-1 volume fraction = 0.666552338169 Min(alpha.water) = -0.000288587367539 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666552338103 Min(alpha.water) = -0.000288576915369 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 3.90201537023e-06, Final residual = 3.77572757143e-08, No Iterations 139 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1.67020370565e-06, Final residual = 1.49001158022e-08, No Iterations 133 time step continuity errors : sum local = 1.86485260925e-13, global = 1.62184431696e-13, cumulative = 2.19829948023e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 8.50239811927e-07, Final residual = 3.72560315988e-09, No Iterations 135 DILUPBiCGStab: Solving for p_rgh, Initial residual = 3.74883768517e-07, Final residual = 5.88727468347e-10, No Iterations 143 time step continuity errors : sum local = 6.87566012445e-15, global = 2.61279932851e-15, cumulative = 2.19830209303e-09 smoothSolver: Solving for omega, Initial residual = 5.82349200687e-10, Final residual = 5.82349200687e-10, No Iterations 0 smoothSolver: Solving for k, Initial residual = 5.00557298894e-09, Final residual = 5.00557298894e-09, No Iterations 0 PIMPLE: iteration 4 smoothSolver: Solving for alpha.water, Initial residual = 5.67042728486e-10, Final residual = 5.67042728486e-10, No Iterations 0 Phase-1 volume fraction = 0.666552338103 Min(alpha.water) = -0.000288576915369 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666552338043 Min(alpha.water) = -0.000288566342977 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 2.27362568217e-07, Final residual = 2.24680912268e-09, No Iterations 148 DILUPBiCGStab: Solving for p_rgh, Initial residual = 9.01785649227e-08, Final residual = 9.24088431816e-10, No Iterations 131 time step continuity errors : sum local = 1.10699528045e-14, global = 6.27749363094e-15, cumulative = 2.19830837052e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 4.51230608273e-08, Final residual = 9.6696452405e-10, No Iterations 117 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1.92996272797e-08, Final residual = 7.57883021183e-10, No Iterations 102 time step continuity errors : sum local = 8.10115202691e-15, global = 8.02059882017e-15, cumulative = 2.19831639112e-09 smoothSolver: Solving for omega, Initial residual = 5.82669063188e-10, Final residual = 5.82669063188e-10, No Iterations 0 smoothSolver: Solving for k, Initial residual = 5.04746641437e-09, Final residual = 5.04746641437e-09, No Iterations 0 PIMPLE: iteration 5 smoothSolver: Solving for alpha.water, Initial residual = 5.98420767334e-10, Final residual = 5.98420767334e-10, No Iterations 0 Phase-1 volume fraction = 0.666552338043 Min(alpha.water) = -0.000288566342977 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666552337988 Min(alpha.water) = -0.000288555539222 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 2.69188395856e-08, Final residual = 9.84231165019e-10, No Iterations 51 DILUPBiCGStab: Solving for p_rgh, Initial residual = 5.26446402637e-09, Final residual = 9.89036434404e-10, No Iterations 4 time step continuity errors : sum local = 1.00675233117e-14, global = 4.92300885785e-15, cumulative = 2.19832131413e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1.20396438674e-09, Final residual = 9.80998673529e-10, No Iterations 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 9.89421021316e-10, Final residual = 9.89421021316e-10, No Iterations 0 time step continuity errors : sum local = 1.01490028677e-14, global = 4.93624069931e-15, cumulative = 2.19832625037e-09 smoothSolver: Solving for omega, Initial residual = 5.82678554279e-10, Final residual = 5.82678554279e-10, No Iterations 0 smoothSolver: Solving for k, Initial residual = 5.04987373213e-09, Final residual = 5.04987373213e-09, No Iterations 0 ExecutionTime = 2526.09 s ClockTime = 2556 s From virtual bool Foam::regIOobject::readIfModified() in file db/regIOobject/regIOobjectRead.C at line 271 Re-reading object fvSchemes from file "/mnt/d/Tutorials/UMaine/decay/decayRASkOmegaSST/decayWECMiddleSpar/decayHeave1V0/background/system/fvSchemes" From virtual bool Foam::regIOobject::readIfModified() in file db/regIOobject/regIOobjectRead.C at line 271 Re-reading object fvSolution from file "/mnt/d/Tutorials/UMaine/decay/decayRASkOmegaSST/decayWECMiddleSpar/decayHeave1V0/background/system/fvSolution" Courant Number mean: 4.99784351636e-07 max: 0.569730783454 Interface Courant Number mean: 5.92138108728e-09 max: 0.0209198035939 deltaT = 0.00014349253782 Time = 0.012703 PIMPLE: iteration 1 forces forces: rho: rho Not including porosity effects forces forces: rho: rho Not including porosity effects forces forces: rho: rho Not including porosity effects Rigid-body motion of the dynamicWECMiddle Centre of rotation: (0 0 -0.999948522777) Orientation: (1 0 0 0 1 0 0 0 1) Linear velocity: (0 0 -0.00588418016272) Angular velocity: (0 0 0) inverseDistance : detected 3 mesh regions zone:0 nCells:3302728 voxels:(92 92 92) bb:(-1250.00076811 -1250.00076811 -800.000768115) (-749.999231885 -749.999231885 -499.999231885) zone:1 nCells:1329243 voxels:(92 92 92) bb:(-17.5000567663 -17.5000567663 -25.0007917663) (17.5000557663 17.5000547663 -1.56244523367) zone:2 nCells:1783736 voxels:(92 92 92) bb:(-4.41261149239 5.32858310761 -12.2503435152) (11.2436789924 17.5000299924 10.2504524696) Overset analysis : nCells : 6415707 calculated : 5795510 interpolated : 159835 (from local:606 mixed local/remote:100 remote:159129) hole : 455019 Updating shallowWaterAbsorption wave model for patch inlet Updating shallowWaterAbsorption wave model for patch outlet Updating shallowWaterAbsorption wave model for patch rightSide Updating shallowWaterAbsorption wave model for patch leftSide Execution time for mesh.update() = 41.31 s smoothSolver: Solving for alpha.water, Initial residual = 4.31650457418e-10, Final residual = 4.31650457418e-10, No Iterations 0 Phase-1 volume fraction = 0.666552337956 Min(alpha.water) = -0.000288555539222 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666552337912 Min(alpha.water) = -0.000288542572206 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 0.00236165217602, Final residual = 2.3271010307e-05, No Iterations 107 DILUPBiCGStab: Solving for p_rgh, Initial residual = 0.000412146953105, Final residual = 4.00471828788e-06, No Iterations 119 time step continuity errors : sum local = 5.53488316472e-11, global = -3.70246236732e-11, cumulative = 2.1613016267e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 5.71606087261e-05, Final residual = 5.21097796183e-07, No Iterations 138 DILUPBiCGStab: Solving for p_rgh, Initial residual = 2.35125965116e-05, Final residual = 9.88359006305e-10, No Iterations 215 time step continuity errors : sum local = 1.51332173524e-14, global = 4.52701018796e-15, cumulative = 2.16130615371e-09 smoothSolver: Solving for omega, Initial residual = 9.65248956233e-07, Final residual = 3.75230694262e-10, No Iterations 1 smoothSolver: Solving for k, Initial residual = 4.42745824449e-05, Final residual = 4.43358495994e-09, No Iterations 2 PIMPLE: iteration 2 smoothSolver: Solving for alpha.water, Initial residual = 5.14219964389e-10, Final residual = 5.14219964389e-10, No Iterations 0 Phase-1 volume fraction = 0.666552337912 Min(alpha.water) = -0.000288542572206 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666552337866 Min(alpha.water) = -0.000288456349238 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1.58024318209e-05, Final residual = 1.50927926767e-07, No Iterations 135 DILUPBiCGStab: Solving for p_rgh, Initial residual = 3.73850293211e-06, Final residual = 5.74727092462e-09, No Iterations 165 time step continuity errors : sum local = 6.69721692409e-14, global = -3.03450383446e-14, cumulative = 2.16127580867e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1.49876301523e-06, Final residual = 1.39670720646e-08, No Iterations 123 DILUPBiCGStab: Solving for p_rgh, Initial residual = 6.76715177075e-07, Final residual = 7.98008738186e-10, No Iterations 148 time step continuity errors : sum local = 1.22344536329e-14, global = -5.78068108562e-15, cumulative = 2.16127002799e-09 smoothSolver: Solving for omega, Initial residual = 6.32837304089e-10, Final residual = 6.32837304089e-10, No Iterations 0 smoothSolver: Solving for k, Initial residual = 1.53271768735e-08, Final residual = 2.46952790962e-09, No Iterations 1 PIMPLE: iteration 3 smoothSolver: Solving for alpha.water, Initial residual = 5.76591398138e-10, Final residual = 5.76591398138e-10, No Iterations 0 Phase-1 volume fraction = 0.666552337866 Min(alpha.water) = -0.000288456349238 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666552337822 Min(alpha.water) = -0.000288366845858 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 5.77537562021e-07, Final residual = 5.48225168085e-09, No Iterations 102 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1.56677431916e-07, Final residual = 1.4141887376e-09, No Iterations 140 time step continuity errors : sum local = 2.23087432416e-14, global = -1.32089190452e-15, cumulative = 2.1612687071e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 6.59080063594e-08, Final residual = 9.9295509487e-10, No Iterations 132 DILUPBiCGStab: Solving for p_rgh, Initial residual = 3.25757162506e-08, Final residual = 9.37883485476e-10, No Iterations 103 time step continuity errors : sum local = 1.46786743791e-14, global = -1.4653311533e-14, cumulative = 2.16125405379e-09 smoothSolver: Solving for omega, Initial residual = 6.49111444261e-10, Final residual = 6.49111444261e-10, No Iterations 0 smoothSolver: Solving for k, Initial residual = 4.0118490653e-09, Final residual = 4.0118490653e-09, No Iterations 0 PIMPLE: iteration 4 smoothSolver: Solving for alpha.water, Initial residual = 6.1973337999e-10, Final residual = 6.1973337999e-10, No Iterations 0 Phase-1 volume fraction = 0.666552337822 Min(alpha.water) = -0.000288366845858 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.66655233778 Min(alpha.water) = -0.000288274846298 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 5.74398048665e-08, Final residual = 9.82136138533e-10, No Iterations 120 DILUPBiCGStab: Solving for p_rgh, Initial residual = 7.02405804242e-09, Final residual = 9.98479118902e-10, No Iterations 9 time step continuity errors : sum local = 1.27292028318e-14, global = -7.9193313664e-16, cumulative = 2.16125326185e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 2.30850519151e-09, Final residual = 9.51027282588e-10, No Iterations 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 9.75182912812e-10, Final residual = 9.75182912812e-10, No Iterations 0 time step continuity errors : sum local = 1.26570808677e-14, global = -4.64721120315e-16, cumulative = 2.16125279713e-09 smoothSolver: Solving for omega, Initial residual = 6.49171974338e-10, Final residual = 6.49171974338e-10, No Iterations 0 smoothSolver: Solving for k, Initial residual = 4.0055169094e-09, Final residual = 4.0055169094e-09, No Iterations 0 PIMPLE: iteration 5 smoothSolver: Solving for alpha.water, Initial residual = 6.63999797793e-10, Final residual = 6.63999797793e-10, No Iterations 0 Phase-1 volume fraction = 0.66655233778 Min(alpha.water) = -0.000288274846298 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666552337738 Min(alpha.water) = -0.00028818033445 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1.61804964183e-08, Final residual = 9.64384311072e-10, No Iterations 8 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1.73777456394e-09, Final residual = 9.66855572087e-10, No Iterations 1 time step continuity errors : sum local = 1.11810878112e-14, global = 1.93801378722e-15, cumulative = 2.16125473514e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1.14956370447e-09, Final residual = 9.54095995662e-10, No Iterations 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 9.55825543332e-10, Final residual = 9.55825543332e-10, No Iterations 0 time step continuity errors : sum local = 1.1122693211e-14, global = 1.9981078044e-15, cumulative = 2.16125673325e-09 smoothSolver: Solving for omega, Initial residual = 6.49156613018e-10, Final residual = 6.49156613018e-10, No Iterations 0 smoothSolver: Solving for k, Initial residual = 4.00623653355e-09, Final residual = 4.00623653355e-09, No Iterations 0 ExecutionTime = 2914.54 s ClockTime = 2959 s -
overset 小间隙-多浮体算例 计算结果明显有误计算结果
/*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2306 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : _fbf00d6bf2-20230626 OPENFOAM=2306 version=v2306 Arch : "LSB;label=32;scalar=64" Exec : overInterDyMFoam -parallel Date : Nov 24 2023 Time : 21:42:20 Host : DESKTOP-KVQEVMV PID : 617 I/O : uncollated Case : /mnt/d/Tutorials/UMaine/decay/decayRASkOmegaSST/decayWECMiddleSpar/decayHeave1V0/background nProcs : 8 Hosts : ( (DESKTOP-KVQEVMV 8) ) Pstream initialized with: floatTransfer : false nProcsSimpleSum : 0 nonBlockingExchange: 0 (tuning: 0) exchange algorithm : -1 commsType : nonBlocking polling iterations : 0 trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20) allowSystemOperations : Allowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 Selecting dynamicFvMesh dynamicOversetFvMesh Selecting motion solver: multiSolidBodyMotionSolver Applying motion to entire mesh Selecting solid-body motion function linearMotion Applying solid body motion linearMotion to 1449278 points of cellZone zonePlatform Selecting motion solver: rigidBodyMotion Applying motion to cellSet: cellWECMiddle Selecting rigidBodySolver Newmark PIMPLE: no residual control data found. Calculations will employ 5 corrector loops Reading field p_rgh Reading field U Reading/calculating face flux field phi Creating cellMask field to block out hole cells --> FOAM Warning : From bool Foam::oversetPolyPatch::master() const in file oversetPolyPatch/oversetPolyPatch.C at line 145 The master overset patch is not the first patch. Generally the first patch should be an overset patch to guarantee consistent operation. --> FOAM Warning : From bool Foam::oversetPolyPatch::master() const in file oversetPolyPatch/oversetPolyPatch.C at line 145 The master overset patch is not the first patch. Generally the first patch should be an overset patch to guarantee consistent operation. Creating interpolatedCells field Reading transportProperties Selecting incompressible transport model Newtonian Selecting incompressible transport model Newtonian Selecting turbulence model type RAS Selecting RAS turbulence model kOmegaSST Selecting patchDistMethod meshWave RAS { RASModel kOmegaSST; turbulence on; printCoeffs on; alphaK1 0.85; alphaK2 1; alphaOmega1 0.5; alphaOmega2 0.856; gamma1 0.555555555556; gamma2 0.44; beta1 0.075; beta2 0.0828; betaStar 0.09; a1 0.31; b1 1; c1 10; F3 false; decayControl false; kInf 0; omegaInf 0; } Reading g Reading hRef Calculating field g.h No MRF models present No finite volume options present Reading/calculating face velocity Uf Courant Number mean: 0 max: 0 Starting time loop Courant Number mean: 0 max: 0 Interface Courant Number mean: 0 max: 0 deltaT = 0.012 Time = 0.012 PIMPLE: iteration 1 forces forces: rho: rho Not including porosity effects forces forces: rho: rho Not including porosity effects forces forces: rho: rho Not including porosity effects Rigid-body motion of the dynamicWECMiddle Centre of rotation: (0 0 -0.999950388828) Orientation: (1 0 0 0 1 0 0 0 1) Linear velocity: (0 0 0.00826852870788) Angular velocity: (0 0 0) inverseDistance : detected 3 mesh regions zone:0 nCells:3302728 voxels:(92 92 92) bb:(-1250.00076811 -1250.00076811 -800.000768115) (-749.999231885 -749.999231885 -499.999231885) zone:1 nCells:1329243 voxels:(92 92 92) bb:(-17.5000567663 -17.5000567663 -25.0007917663) (17.5000557663 17.5000547663 -1.56244523367) zone:2 nCells:1783736 voxels:(92 92 92) bb:(-4.41261149239 5.32858310761 -12.2503453812) (11.2436789924 17.5000299924 10.2504506036) Overset analysis : nCells : 6415707 calculated : 5795510 interpolated : 159835 (from local:606 mixed local/remote:100 remote:159129) hole : 455019 Selecting waveModel shallowWaterAbsorption Wave model: patch inlet Type : shallowWaterAbsorption Velocity field name : U Phase fraction field name : alpha.water Transformation from local to global system : (1 -7.77501308186e-32 7.17951403026e-23 7.77501308186e-32 1 0 -0 -0 1) Number of paddles: 1 Reference water depth : 199.999999317 Active absorption: 1 Updating shallowWaterAbsorption wave model for patch inlet Selecting waveModel shallowWaterAbsorption Wave model: patch outlet Type : shallowWaterAbsorption Velocity field name : U Phase fraction field name : alpha.water Transformation from local to global system : (-1 -3.16385813721e-32 7.17951403181e-23 3.16385813721e-32 -1 0 -0 -0 1) Number of paddles: 1 Reference water depth : 199.999999317 Active absorption: 1 Updating shallowWaterAbsorption wave model for patch outlet Selecting waveModel shallowWaterAbsorption Wave model: patch rightSide Type : shallowWaterAbsorption Velocity field name : U Phase fraction field name : alpha.water Transformation from local to global system : (-1.61558713389e-32 1 1.05282388693e-18 -1 -1.61558713389e-32 -0 -0 -0 1) Number of paddles: 1 Reference water depth : 199.999999317 Active absorption: 1 Updating shallowWaterAbsorption wave model for patch rightSide Selecting waveModel shallowWaterAbsorption Wave model: patch leftSide Type : shallowWaterAbsorption Velocity field name : U Phase fraction field name : alpha.water Transformation from local to global system : (2.15411617852e-32 -1 -5.88975123258e-18 1 2.15411617852e-32 0 -0 -0 1) Number of paddles: 1 Reference water depth : 199.999999317 Active absorption: 1 Updating shallowWaterAbsorption wave model for patch leftSide Execution time for mesh.update() = 31.64 s smoothSolver: Solving for alpha.water, Initial residual = 2.61985281096e-07, Final residual = 9.3161532224e-09, No Iterations 23 Phase-1 volume fraction = 0.66655233285 Min(alpha.water) = 0 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.66655233262 Min(alpha.water) = -4.06575814682e-23 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1, Final residual = 0.00934508785126, No Iterations 132 DILUPBiCGStab: Solving for p_rgh, Initial residual = 3.93805455326e-05, Final residual = 9.81032812434e-08, No Iterations 157 time step continuity errors : sum local = 1.31199443232e-09, global = 5.31630376358e-10, cumulative = 5.31630376358e-10 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1.05167563797e-05, Final residual = 1.04029750059e-07, No Iterations 144 DILUPBiCGStab: Solving for p_rgh, Initial residual = 3.3336587163e-06, Final residual = 6.33395425015e-10, No Iterations 233 time step continuity errors : sum local = 6.1751513569e-12, global = 5.28990371798e-12, cumulative = 5.36920280076e-10 smoothSolver: Solving for omega, Initial residual = 0.000102569828088, Final residual = 2.86343248998e-09, No Iterations 3 smoothSolver: Solving for k, Initial residual = 0.999999998453, Final residual = 9.33528038491e-09, No Iterations 66 PIMPLE: iteration 2 smoothSolver: Solving for alpha.water, Initial residual = 1.27402311839e-07, Final residual = 8.81590646774e-09, No Iterations 2 Phase-1 volume fraction = 0.666552379444 Min(alpha.water) = -3.83779065969e-28 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666552381689 Min(alpha.water) = -1.621967531e-06 Max(alpha.water) = 1.00000000013 DILUPBiCGStab: Solving for p_rgh, Initial residual = 9.61393401599e-05, Final residual = 2.95529990406e-07, No Iterations 137 DILUPBiCGStab: Solving for p_rgh, Initial residual = 9.32733975555e-05, Final residual = 8.98829957804e-07, No Iterations 58 time step continuity errors : sum local = 2.97727638562e-09, global = 1.12545136858e-09, cumulative = 1.66237164866e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 8.25486710901e-05, Final residual = 7.9528642603e-07, No Iterations 111 DILUPBiCGStab: Solving for p_rgh, Initial residual = 2.1598896674e-05, Final residual = 2.69740404024e-10, No Iterations 215 time step continuity errors : sum local = 8.87609898546e-13, global = 2.19947839059e-13, cumulative = 1.6625915965e-09 smoothSolver: Solving for omega, Initial residual = 1.05706805697e-06, Final residual = 3.06165398729e-09, No Iterations 2 smoothSolver: Solving for k, Initial residual = 2.17813199337e-05, Final residual = 8.83104952045e-09, No Iterations 11 PIMPLE: iteration 3 smoothSolver: Solving for alpha.water, Initial residual = 1.08796870261e-07, Final residual = 6.88901360938e-09, No Iterations 3 Phase-1 volume fraction = 0.666552350255 Min(alpha.water) = -9.78968376404e-18 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666552348822 Min(alpha.water) = -4.03306452835e-06 Max(alpha.water) = 1.00000000001 DILUPBiCGStab: Solving for p_rgh, Initial residual = 0.000260128355201, Final residual = 2.38866841716e-06, No Iterations 80 DILUPBiCGStab: Solving for p_rgh, Initial residual = 3.97760690631e-05, Final residual = 3.74084676736e-07, No Iterations 98 time step continuity errors : sum local = 6.45059786074e-10, global = 6.27376936037e-10, cumulative = 2.28996853253e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 5.18893127841e-05, Final residual = 4.93880984564e-07, No Iterations 140 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1.00033869064e-05, Final residual = 9.88435829622e-10, No Iterations 215 time step continuity errors : sum local = 1.79260009404e-12, global = 8.66436780833e-14, cumulative = 2.29005517621e-09 smoothSolver: Solving for omega, Initial residual = 2.53541232093e-07, Final residual = 3.5196014854e-09, No Iterations 2 smoothSolver: Solving for k, Initial residual = 2.94650853902e-05, Final residual = 7.53900830036e-09, No Iterations 11 PIMPLE: iteration 4 smoothSolver: Solving for alpha.water, Initial residual = 7.59424341536e-08, Final residual = 8.47687199918e-09, No Iterations 1 Phase-1 volume fraction = 0.666552350369 Min(alpha.water) = -1.7236282171e-09 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666552350436 Min(alpha.water) = -0.000416460503671 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 0.000167246351824, Final residual = 1.40182075997e-06, No Iterations 107 DILUPBiCGStab: Solving for p_rgh, Initial residual = 5.16507599585e-05, Final residual = 4.37468479875e-07, No Iterations 64 time step continuity errors : sum local = 6.00063571489e-10, global = -2.63397570915e-10, cumulative = 2.0266576053e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 8.29697246248e-05, Final residual = 8.11491883613e-07, No Iterations 151 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1.44922341191e-05, Final residual = 6.85329760483e-10, No Iterations 217 time step continuity errors : sum local = 1.22184650642e-12, global = -9.18162980737e-13, cumulative = 2.02573944232e-09 smoothSolver: Solving for omega, Initial residual = 1.52703236419e-07, Final residual = 7.73512960224e-09, No Iterations 1 smoothSolver: Solving for k, Initial residual = 8.02738033721e-06, Final residual = 9.37089259035e-09, No Iterations 23 PIMPLE: iteration 5 smoothSolver: Solving for alpha.water, Initial residual = 5.62034902794e-08, Final residual = 8.1720263608e-09, No Iterations 1 Phase-1 volume fraction = 0.666552346355 Min(alpha.water) = -3.66106998989e-10 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666552344225 Min(alpha.water) = -0.000289714267251 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 0.000179222117898, Final residual = 1.69745284627e-06, No Iterations 123 DILUPBiCGStab: Solving for p_rgh, Initial residual = 2.19874219811e-05, Final residual = 1.89807848936e-07, No Iterations 148 time step continuity errors : sum local = 3.71990181733e-10, global = 2.12074694904e-11, cumulative = 2.04694691181e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 3.41850750136e-05, Final residual = 3.01352648671e-07, No Iterations 126 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1.45426477763e-05, Final residual = 9.45301356671e-10, No Iterations 178 time step continuity errors : sum local = 1.78223171988e-12, global = -2.78459422187e-13, cumulative = 2.04666845238e-09 smoothSolver: Solving for omega, Initial residual = 1.60789297939e-07, Final residual = 8.69888544542e-09, No Iterations 1 smoothSolver: Solving for k, Initial residual = 4.85804248814e-06, Final residual = 9.50862177947e-09, No Iterations 7 ExecutionTime = 557.59 s ClockTime = 560 s Courant Number mean: 1.63749449277e-05 max: 49.9913602954 Interface Courant Number mean: 6.95403940679e-07 max: 3.55036474521 deltaT = 0.000240041477749 Time = 0.01224 PIMPLE: iteration 1 forces forces: rho: rho Not including porosity effects forces forces: rho: rho Not including porosity effects forces forces: rho: rho Not including porosity effects Rigid-body motion of the dynamicWECMiddle Centre of rotation: (0 0 -0.999948351772) Orientation: (1 0 0 0 1 0 0 0 1) Linear velocity: (0 0 0.00870400248066) Angular velocity: (0 0 0) inverseDistance : detected 3 mesh regions zone:0 nCells:3302728 voxels:(92 92 92) bb:(-1250.00076811 -1250.00076811 -800.000768115) (-749.999231885 -749.999231885 -499.999231885) zone:1 nCells:1329243 voxels:(92 92 92) bb:(-17.5000567663 -17.5000567663 -25.0007917663) (17.5000557663 17.5000547663 -1.56244523367) zone:2 nCells:1783736 voxels:(92 92 92) bb:(-4.41261149239 5.32858310761 -12.2503433442) (11.2436789924 17.5000299924 10.2504526406) Overset analysis : nCells : 6415707 calculated : 5795510 interpolated : 159835 (from local:606 mixed local/remote:100 remote:159129) hole : 455019 Updating shallowWaterAbsorption wave model for patch inlet Updating shallowWaterAbsorption wave model for patch outlet Updating shallowWaterAbsorption wave model for patch rightSide Updating shallowWaterAbsorption wave model for patch leftSide Execution time for mesh.update() = 40.12 s smoothSolver: Solving for alpha.water, Initial residual = 8.36900182683e-09, Final residual = 8.36900182683e-09, No Iterations 0 Phase-1 volume fraction = 0.666552343446 Min(alpha.water) = -0.000289714267251 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666552342656 Min(alpha.water) = -0.000289584326027 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 0.0228028713103, Final residual = 0.000158123063473, No Iterations 103 DILUPBiCGStab: Solving for p_rgh, Initial residual = 0.000476023249144, Final residual = 4.19682603135e-06, No Iterations 132 time step continuity errors : sum local = 8.93614961201e-11, global = 7.31703459824e-11, cumulative = 2.11983879837e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 0.000192656249381, Final residual = 1.36166337649e-06, No Iterations 163 DILUPBiCGStab: Solving for p_rgh, Initial residual = 6.70954837961e-05, Final residual = 9.678157846e-10, No Iterations 265 time step continuity errors : sum local = 1.58991434156e-14, global = -7.99361856081e-15, cumulative = 2.11983080475e-09 smoothSolver: Solving for omega, Initial residual = 1.65654710594e-06, Final residual = 9.11101653564e-10, No Iterations 1 smoothSolver: Solving for k, Initial residual = 9.84944919853e-05, Final residual = 7.48999418113e-09, No Iterations 2 PIMPLE: iteration 2 smoothSolver: Solving for alpha.water, Initial residual = 6.06658567684e-09, Final residual = 6.06658567684e-09, No Iterations 0 Phase-1 volume fraction = 0.666552342656 Min(alpha.water) = -0.000289584326027 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666552341998 Min(alpha.water) = -0.000289373569159 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 6.09071653377e-05, Final residual = 5.86498406611e-07, No Iterations 135 DILUPBiCGStab: Solving for p_rgh, Initial residual = 2.67956672479e-05, Final residual = 2.01508297651e-07, No Iterations 149 time step continuity errors : sum local = 3.3067714853e-12, global = 3.30648753556e-12, cumulative = 2.12313729228e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1.18188041175e-05, Final residual = 7.73661385616e-08, No Iterations 156 DILUPBiCGStab: Solving for p_rgh, Initial residual = 5.55315884803e-06, Final residual = 6.90107537102e-10, No Iterations 213 time step continuity errors : sum local = 1.07514813361e-14, global = 1.01727697345e-14, cumulative = 2.12314746505e-09 smoothSolver: Solving for omega, Initial residual = 1.82307346832e-09, Final residual = 1.82307346832e-09, No Iterations 0 smoothSolver: Solving for k, Initial residual = 6.07724465011e-08, Final residual = 3.64698977772e-09, No Iterations 1 PIMPLE: iteration 3 smoothSolver: Solving for alpha.water, Initial residual = 4.95615350895e-09, Final residual = 4.95615350895e-09, No Iterations 0 Phase-1 volume fraction = 0.666552341998 Min(alpha.water) = -0.000289373569159 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666552341433 Min(alpha.water) = -0.00028921797597 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 5.34942515551e-06, Final residual = 2.89216757541e-08, No Iterations 153 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1.71613055767e-06, Final residual = 1.69845429434e-08, No Iterations 139 time step continuity errors : sum local = 2.85144328781e-13, global = 2.82486392054e-13, cumulative = 2.12342995145e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 7.57358817177e-07, Final residual = 4.68743757682e-09, No Iterations 148 DILUPBiCGStab: Solving for p_rgh, Initial residual = 3.30816263038e-07, Final residual = 7.98377737781e-10, No Iterations 155 time step continuity errors : sum local = 1.24636818568e-14, global = 1.09792420598e-14, cumulative = 2.12344093069e-09 smoothSolver: Solving for omega, Initial residual = 1.92226875551e-09, Final residual = 1.92226875551e-09, No Iterations 0 smoothSolver: Solving for k, Initial residual = 1.34277181319e-08, Final residual = 2.5244507753e-09, No Iterations 1 PIMPLE: iteration 4 smoothSolver: Solving for alpha.water, Initial residual = 4.14180619703e-09, Final residual = 4.14180619703e-09, No Iterations 0 Phase-1 volume fraction = 0.666552341433 Min(alpha.water) = -0.00028921797597 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.66655234095 Min(alpha.water) = -0.000289070037242 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 9.69585590378e-07, Final residual = 8.74528894708e-09, No Iterations 140 DILUPBiCGStab: Solving for p_rgh, Initial residual = 2.81105228656e-07, Final residual = 2.72909179203e-09, No Iterations 126 time step continuity errors : sum local = 4.56821656802e-14, global = 1.54082829784e-14, cumulative = 2.12345633897e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 9.25226124652e-08, Final residual = 6.67713245408e-10, No Iterations 124 DILUPBiCGStab: Solving for p_rgh, Initial residual = 4.2978743237e-08, Final residual = 9.74587014884e-10, No Iterations 109 time step continuity errors : sum local = 1.66374816607e-14, global = 1.64971601989e-14, cumulative = 2.12347283613e-09 smoothSolver: Solving for omega, Initial residual = 1.95628819863e-09, Final residual = 1.95628819863e-09, No Iterations 0 smoothSolver: Solving for k, Initial residual = 5.62125672229e-09, Final residual = 5.62125672229e-09, No Iterations 0 PIMPLE: iteration 5 smoothSolver: Solving for alpha.water, Initial residual = 3.60189757848e-09, Final residual = 3.60189757848e-09, No Iterations 0 Phase-1 volume fraction = 0.66655234095 Min(alpha.water) = -0.000289070037242 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666552340533 Min(alpha.water) = -0.000288909525708 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 3.67847330072e-07, Final residual = 1.98395076727e-09, No Iterations 149 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1.06975221441e-07, Final residual = 1.02026137309e-09, No Iterations 136 time step continuity errors : sum local = 1.74929356065e-14, global = 1.46083619252e-14, cumulative = 2.12348744449e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 4.31574420273e-08, Final residual = 9.65472946334e-10, No Iterations 119 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1.49455168275e-08, Final residual = 8.73609993614e-10, No Iterations 102 time step continuity errors : sum local = 1.40323018836e-14, global = 2.66664498002e-15, cumulative = 2.12349011114e-09 smoothSolver: Solving for omega, Initial residual = 1.95616404291e-09, Final residual = 1.95616404291e-09, No Iterations 0 smoothSolver: Solving for k, Initial residual = 5.67015345545e-09, Final residual = 5.67015345545e-09, No Iterations 0 ExecutionTime = 1116.64 s ClockTime = 1124 s Courant Number mean: 1.23753094205e-06 max: 2.67476268045 Interface Courant Number mean: 1.22084714043e-08 max: 0.040419192207 deltaT = 8.97430936596e-05 Time = 0.0123298 PIMPLE: iteration 1 forces forces: rho: rho Not including porosity effects forces forces: rho: rho Not including porosity effects forces forces: rho: rho Not including porosity effects Rigid-body motion of the dynamicWECMiddle Centre of rotation: (0 0 -0.999947523286) Orientation: (1 0 0 0 1 0 0 0 1) Linear velocity: (0 0 0.00975950626086) Angular velocity: (0 0 0) inverseDistance : detected 3 mesh regions zone:0 nCells:3302728 voxels:(92 92 92) bb:(-1250.00076811 -1250.00076811 -800.000768115) (-749.999231885 -749.999231885 -499.999231885) zone:1 nCells:1329243 voxels:(92 92 92) bb:(-17.5000567663 -17.5000567663 -25.0007917663) (17.5000557663 17.5000547663 -1.56244523367) zone:2 nCells:1783736 voxels:(92 92 92) bb:(-4.41261149239 5.32858310761 -12.2503425157) (11.2436789924 17.5000299924 10.2504534691) Overset analysis : nCells : 6415707 calculated : 5795510 interpolated : 159835 (from local:606 mixed local/remote:100 remote:159129) hole : 455019 Updating shallowWaterAbsorption wave model for patch inlet Updating shallowWaterAbsorption wave model for patch outlet Updating shallowWaterAbsorption wave model for patch rightSide Updating shallowWaterAbsorption wave model for patch leftSide Execution time for mesh.update() = 39.64 s smoothSolver: Solving for alpha.water, Initial residual = 1.89506393044e-09, Final residual = 1.89506393044e-09, No Iterations 0 Phase-1 volume fraction = 0.666552340293 Min(alpha.water) = -0.000288909525708 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666552339977 Min(alpha.water) = -0.000288846889594 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 0.00695412581841, Final residual = 6.58530505703e-05, No Iterations 148 DILUPBiCGStab: Solving for p_rgh, Initial residual = 0.000208251042114, Final residual = 2.05743333992e-06, No Iterations 125 time step continuity errors : sum local = 5.895261321e-11, global = -3.50991476391e-11, cumulative = 2.0883909635e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 0.00010380483611, Final residual = 9.76125551947e-07, No Iterations 112 DILUPBiCGStab: Solving for p_rgh, Initial residual = 4.2687177802e-05, Final residual = 9.50115732418e-10, No Iterations 213 time step continuity errors : sum local = 2.05769533909e-14, global = 1.23631828244e-14, cumulative = 2.08840332668e-09 smoothSolver: Solving for omega, Initial residual = 6.1051758878e-07, Final residual = 4.82103625934e-10, No Iterations 1 smoothSolver: Solving for k, Initial residual = 3.01714859593e-05, Final residual = 9.37540222535e-09, No Iterations 1 PIMPLE: iteration 2 smoothSolver: Solving for alpha.water, Initial residual = 1.63714691561e-09, Final residual = 1.63714691561e-09, No Iterations 0 Phase-1 volume fraction = 0.666552339977 Min(alpha.water) = -0.000288846889594 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666552339706 Min(alpha.water) = -0.000288806485034 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 2.23761401106e-05, Final residual = 2.08661810278e-07, No Iterations 134 DILUPBiCGStab: Solving for p_rgh, Initial residual = 9.2229647521e-06, Final residual = 8.84903438452e-08, No Iterations 146 time step continuity errors : sum local = 1.93794820853e-12, global = 1.93793640713e-12, cumulative = 2.09034126309e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 4.31258175956e-06, Final residual = 3.09958221974e-08, No Iterations 129 DILUPBiCGStab: Solving for p_rgh, Initial residual = 2.16173305739e-06, Final residual = 7.61858990072e-10, No Iterations 182 time step continuity errors : sum local = 1.79845063561e-14, global = -7.97949298118e-15, cumulative = 2.0903332836e-09 smoothSolver: Solving for omega, Initial residual = 6.86324608392e-10, Final residual = 6.86324608392e-10, No Iterations 0 smoothSolver: Solving for k, Initial residual = 2.59199319034e-08, Final residual = 3.72155986805e-09, No Iterations 1 PIMPLE: iteration 3 smoothSolver: Solving for alpha.water, Initial residual = 1.464298727e-09, Final residual = 1.464298727e-09, No Iterations 0 Phase-1 volume fraction = 0.666552339706 Min(alpha.water) = -0.000288806485034 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666552339468 Min(alpha.water) = -0.000288763304213 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1.11368632855e-06, Final residual = 9.12572869945e-09, No Iterations 145 DILUPBiCGStab: Solving for p_rgh, Initial residual = 4.51282511339e-07, Final residual = 4.39016946545e-09, No Iterations 148 time step continuity errors : sum local = 9.74366513151e-14, global = -9.73777380055e-14, cumulative = 2.09023590586e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 2.01376950811e-07, Final residual = 9.97729577993e-10, No Iterations 132 DILUPBiCGStab: Solving for p_rgh, Initial residual = 9.03614929149e-08, Final residual = 6.77515000163e-10, No Iterations 146 time step continuity errors : sum local = 1.65113572261e-14, global = -1.60186025025e-14, cumulative = 2.09021988725e-09 smoothSolver: Solving for omega, Initial residual = 6.93072732889e-10, Final residual = 6.93072732889e-10, No Iterations 0 smoothSolver: Solving for k, Initial residual = 4.6403143789e-09, Final residual = 4.6403143789e-09, No Iterations 0 PIMPLE: iteration 4 smoothSolver: Solving for alpha.water, Initial residual = 1.34365518119e-09, Final residual = 1.34365518119e-09, No Iterations 0 Phase-1 volume fraction = 0.666552339468 Min(alpha.water) = -0.000288763304213 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.66655233926 Min(alpha.water) = -0.000288719352487 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 7.66014318067e-08, Final residual = 8.71802548845e-10, No Iterations 148 DILUPBiCGStab: Solving for p_rgh, Initial residual = 2.2919023473e-08, Final residual = 9.70813176141e-10, No Iterations 97 time step continuity errors : sum local = 2.14135270487e-14, global = -1.99862439702e-14, cumulative = 2.09019990101e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1.08856570606e-08, Final residual = 9.76022773547e-10, No Iterations 58 DILUPBiCGStab: Solving for p_rgh, Initial residual = 4.07307604617e-09, Final residual = 9.29257136903e-10, No Iterations 21 time step continuity errors : sum local = 1.67526610845e-14, global = 8.81445888723e-15, cumulative = 2.09020871547e-09 smoothSolver: Solving for omega, Initial residual = 6.93390357342e-10, Final residual = 6.93390357342e-10, No Iterations 0 smoothSolver: Solving for k, Initial residual = 4.67499439824e-09, Final residual = 4.67499439824e-09, No Iterations 0 PIMPLE: iteration 5 smoothSolver: Solving for alpha.water, Initial residual = 1.2520957575e-09, Final residual = 1.2520957575e-09, No Iterations 0 Phase-1 volume fraction = 0.66655233926 Min(alpha.water) = -0.000288719352487 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666552339079 Min(alpha.water) = -0.000288674615659 Max(alpha.water) = 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1.50640815229e-08, Final residual = 9.99499573606e-10, No Iterations 31 DILUPBiCGStab: Solving for p_rgh, Initial residual = 2.01438500186e-09, Final residual = 9.86510742554e-10, No Iterations 10 time step continuity errors : sum local = 1.8597559213e-14, global = 9.98907756168e-15, cumulative = 2.09021870455e-09 DILUPBiCGStab: Solving for p_rgh, Initial residual = 1.37064254853e-09, Final residual = 9.87017098349e-10, No Iterations 1 DILUPBiCGStab: Solving for p_rgh, Initial residual = 9.9024984935e-10, Final residual = 9.9024984935e-10, No Iterations 0 time step continuity errors : sum local = 1.88156236059e-14, global = 9.85210599887e-15, cumulative = 2.09022855665e-09 smoothSolver: Solving for omega, Initial residual = 6.93401046968e-10, Final residual = 6.93401046968e-10, No Iterations 0 smoothSolver: Solving for k, Initial residual = 4.67636350062e-09, Final residual = 4.67636350062e-09, No Iterations 0 ExecutionTime = 1551.51 s ClockTime = 1567 s Courant Number mean: 7.44659378168e-07 max: 0.33617200472 Interface Courant Number mean: 5.92610839418e-09 max: 0.0133063005972 deltaT = 0.000107691712392 Time = 0.0124375 PIMPLE: iteration 1 forces forces: rho: rho Not including porosity effects forces forces: rho: rho Not including porosity effects forces forces: rho: rho Not including porosity effects Rigid-body motion of the dynamicWECMiddle Centre of rotation: (0 0 -0.999946925817) Orientation: (1 0 0 0 1 0 0 0 1) Linear velocity: (0 0 0.00133640649968) Angular velocity: (0 0 0) inverseDistance : detected 3 mesh regions zone:0 nCells:3302728 voxels:(92 92 92) bb:(-1250.00076811 -1250.00076811 -800.000768115) (-749.999231885 -749.999231885 -499.999231885) zone:1 nCells:1329243 voxels:(92 92 92) bb:(-17.5000567663 -17.5000567663 -25.0007917663) (17.5000557663 17.5000547663 -1.56244523367) zone:2 nCells:1783736 voxels:(92 92 92) bb:(-4.41261149239 5.32858310761 -12.2503419182) (11.2436789924 17.5000299924 10.2504540666) Overset analysis : nCells : 6415707 calculated : 5795510 interpolated : 159835 (from local:606 mixed local/remote:100 remote:159129) hole : 455019 -
overset 小间隙-多浮体算例 计算结果明显有误最近在做同轴双浮体装置，目前是内浮体不动（后续会动，所以做成了overset），外浮体初始位移为-1做自衰减
内浮体的overset较大，将外浮体overset包裹，四面消波处理
目前计算结果，外浮体在最初的几秒存在负速度，这是完全不合理的，之前尝试过单独外浮体/单独内浮体/外浮体overset+内浮体背景扣除都是没有问题的。怀疑是overset插值引起的错误结果，请大家帮忙看看是什么问题引起的呢？
fvSchemes：
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2306 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object fvSchemes; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // ddtSchemes { default Euler; } gradSchemes { default Gauss linear; } divSchemes { div(rhoPhi,U) Gauss limitedLinearV 1; div(U) Gauss linear; div(phi,alpha) Gauss vanLeer; div(phirb,alpha) Gauss linear; div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; div(phi,k) Gauss upwind; div(phi,epsilon) Gauss upwind; div(phi,omega) Gauss upwind; } laplacianSchemes { default Gauss linear corrected; } interpolationSchemes { default linear; } snGradSchemes { default corrected; } oversetInterpolation { method inverseDistance; } fluxRequired { default no; p_rgh; pcorr; alpha.water; } wallDist { method meshWave; //method Poisson; } oversetInterpolationSuppressed { grad(p_rgh); surfaceIntegrate(phiHbyA); } /* oversetInterpolationRequired { k; epsilon; omega; nut; alpha.water; } */ // ************************************************************************* //fvSolution:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2306 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object fvSolution; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // solvers { "cellDisplacement.*" { solver PCG; preconditioner DIC; tolerance 1e-06; relTol 0; maxIter 100; } "alpha.water.*" { nAlphaCorr 2; nAlphaSubCycles 1; cAlpha 1; icAlpha 0; MULESCorr yes; nLimiterIter 5; alphaApplyPrevCorr no; solver smoothSolver; smoother symGaussSeidel; tolerance 1e-08; relTol 0; } "pcorr.*" { solver PCG; preconditioner DIC; tolerance 1e-09; relTol 0; } p_rgh { solver PBiCGStab; preconditioner DILU; tolerance 1e-09; relTol 0.01; } p_rghFinal { $p_rgh; relTol 0; } "(U|k|omega|epsilon).*" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-08; relTol 0; } yPsi { solver PBiCGStab; preconditioner DILU; tolerance 1e-6; relTol 0.0; } yPsiFinal { $yPsi; relTol 0; } } PIMPLE { momentumPredictor no; nOuterCorrectors 5; nCorrectors 2; nNonOrthogonalCorrectors 1; ddtCorr yes; correctPhi no; moveMeshOuterCorrectors no; turbOnFinalIterOnly no; oversetAdjustPhi yes; } relaxationFactors { fields { } equations { ".*" 1; } } cache { grad(U); } // ************************************************************************* //dynamicMeshDict
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2306 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object dynamicMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dynamicFvMesh dynamicOversetFvMesh; solvers { dynamicPlatform { motionSolverLibs (libfvMotionSolvers); motionSolver multiSolidBodyMotionSolver; multiSolidBodyMotionSolverCoeffs { zonePlatform { solidBodyMotionFunction linearMotion; linearMotionCoeffs { velocity (0 0 0); } } } } dynamicWECMiddle { motionSolverLibs (librigidBodyMeshMotion); motionSolver rigidBodyMotion; report on; cellSet cellWECMiddle; cOfGdisplacement CofG; bodyIdCofG 1; solver { type Newmark; } accelerationRelaxation 0.8; accelerationDamping 0.9; nIter 3; bodies { dynamicWECMiddle { type rigidBody; parent root; mass 2.0126e+06; centreOfMass (0 0 0); // COG in the local coordinate system, based on transform inertia (1.3312e+08 0 0 1.3312e+08 0 2.1384e+08); // (Ixx Ixy Ixz Iyy Iyz Izz) transform (1 0 0 0 1 0 0 0 1) (0 0 -1); joint { type Pz; } patches ( WECMiddle ); innerDistance 100; outerDistance 101; } } restraints { } /* motionSolverLibs (sixDoFRigidBodyMotion); motionSolver sixDoFRigidBodyMotion; cellSet cellWECMiddle; patches (WECMiddle); innerDistance 100; outerDistance 101; centreOfMass (0 0 -1); cOfGdisplacement CofG; mass 2.0126e+06; momentOfInertia (1.3312e+08 1.3312e+08 2.1384e+08); report on; accelerationRelaxation 0.8; accelerationDamping 0.9; solver { type Newmark; } constraints { zAxis { sixDoFRigidBodyMotionConstraint line; direction (0 0 1); } } */ } } // ************************************************************************* //0:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2306 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volScalarField; object alpha.water; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 0 0 0 0 0 0]; internalField uniform 0; boundaryField { #includeEtc "caseDicts/setConstraintTypes" inlet { type zeroGradient; } outlet { type zeroGradient; } "(rightSide|leftSide)" { type zeroGradient; } lowerSide { type zeroGradient; } atmosphere { type inletOutlet; inletValue $internalField; value $internalField; } oversetWECMiddle { type overset; } WECMiddle { type zeroGradient; } oversetPlatform { type overset; } platform { type zeroGradient; } } // ************************************************************************* // /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2306 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volScalarField; object p_rgh; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 -1 -2 0 0 0 0]; internalField uniform 0; boundaryField { #includeEtc "caseDicts/setConstraintTypes" inlet { type fixedFluxPressure; value $internalField; } outlet { type fixedFluxPressure; value $internalField; } "(rightSide|leftSide)" { type fixedFluxPressure; value $internalField; } lowerSide { type fixedFluxPressure; value $internalField; } atmosphere { type totalPressure; p0 $internalField; } oversetWECMiddle { type overset; } WECMiddle { type fixedFluxPressure; } oversetPlatform { type overset; } platform { type fixedFluxPressure; } } // ************************************************************************* // /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2306 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class pointVectorField; object pointDisplacement; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 0 0 0 0 0]; internalField uniform (0 0 0); boundaryField { #includeEtc "caseDicts/setConstraintTypes" inlet { type fixedValue; value $internalField; } outlet { type fixedValue; value $internalField; } "(rightSide|leftSide)" { type fixedValue; value $internalField; } lowerSide { type fixedValue; value $internalField; } atmosphere { type fixedValue; value $internalField; } oversetWECMiddle { patchType overset; type zeroGradient; } WECMiddle { type calculated; value $internalField; } oversetPlatform { patchType overset; type zeroGradient; } platform { type fixedValue; value $internalField; } } // ************************************************************************* // /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2306 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volVectorField; object U; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { #includeEtc "caseDicts/setConstraintTypes" inlet { type waveVelocity; value $internalField; } outlet { type waveVelocity; value $internalField; } "(rightSide|leftSide)" { type waveVelocity; value $internalField; } lowerSide { type fixedValue; value $internalField; } atmosphere { type pressureInletOutletVelocity; value $internalField; } oversetWECMiddle { type overset; } WECMiddle { type movingWallVelocity; value $internalField; } oversetPlatform { type overset; } platform { type fixedValue; value $internalField; } } // ************************************************************************* //请大家帮忙看看是不是哪里设置有问题
-
openfoam模拟船舶运动参考multiphase\overInterDyMFoam\rigidBodyHull
-
OpenFOAM 浮体铰接可以计算吗可以实现 主要还是动网格方面 可以参考rigidBodyMotion方法
-
ofv2212版本里的rigidbodyhull案例的一些疑问你好 还有一个问题不知道你有没有注意到：rigidBodyHull算例似乎是把螺旋桨overset嵌在船体overset中，这样不会导致背景网格在与船体overset插值时认为螺旋桨内存在水体吗
-
mpirun detected that one or more processes exited with non-zero statushttps://cfd-china.com/user/香柏树
请问您是重新计算 还是只计算最新几步的结果
我删除了最后几个时间步 重新计算 依旧迭代不出来PIMPLE: iteration 1 Updating StokesI wave model for patch inlet smoothSolver: Solving for alpha.water, Initial residual = 2.85603060095e-08, Final residual = 1.44420404532e-10, No Iterations 1 Phase-1 volume fraction = 0.66661826176 Min(alpha.water) = -5.22905791777e-10 Max(alpha.water) = 1.00014678682 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.66661826176 Min(alpha.water) = -5.22905791774e-10 Max(alpha.water) = 1.00014678682 Updating shallowWaterAbsorption wave model for patch outlet GAMGPCG: Solving for p_rgh, Initial residual = 0.833496905521, Final residual = 5.36596999747e-06, No Iterations 9 GAMGPCG: Solving for p_rgh, Initial residual = 0.748404838772, Final residual = 4.322497918e-06, No Iterations 9 time step continuity errors : sum local = 6.40365228609e-12, global = -3.66804749067e-13, cumulative = 2.06361704946e-05 GAMGPCG: Solving for p_rgh, Initial residual = 0.804959354481, Final residual = 4.4360097837e-06, No Iterations 9 GAMGPCG: Solving for p_rgh, Initial residual = 0.834220247626, Final residual = 4.86674559877e-06, No Iterations 9 time step continuity errors : sum local = 2.88518833016e-11, global = -4.33468036047e-13, cumulative = 2.06361700612e-05 smoothSolver: Solving for omega, Initial residual = 8.63899230563e-07, Final residual = 4.28995095404e-10, No Iterations 2 smoothSolver: Solving for k, Initial residual = 4.09430928447e-06, Final residual = 1.30320142023e-09, No Iterations 2 ExecutionTime = 232.9 s ClockTime = 234 s Courant Number mean: 5.31856748554e-06 max: 6.80956597246 Interface Courant Number mean: 4.94232067864e-07 max: 6.80956597246 deltaT = 2.06579061066e-06 Time = 67.1303 PIMPLE: iteration 1 Updating StokesI wave model for patch inlet smoothSolver: Solving for alpha.water, Initial residual = 1.69468041505e-08, Final residual = 5.40084448828e-11, No Iterations 1 Phase-1 volume fraction = 0.666618261841 Min(alpha.water) = -4.88357838579e-10 Max(alpha.water) = 1.00014678669 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666618261841 Min(alpha.water) = -4.88070715976e-10 Max(alpha.water) = 1.00014678669 Updating shallowWaterAbsorption wave model for patch outlet GAMGPCG: Solving for p_rgh, Initial residual = 0.856084241578, Final residual = 3.91327711378e-06, No Iterations 9 GAMGPCG: Solving for p_rgh, Initial residual = 0.729556404756, Final residual = 4.52724813344e-06, No Iterations 9 time step continuity errors : sum local = 4.23132086572e-12, global = -1.39652678751e-13, cumulative = 2.06361699215e-05 GAMGPCG: Solving for p_rgh, Initial residual = 0.824421265103, Final residual = 5.74632165117e-06, No Iterations 9 GAMGPCG: Solving for p_rgh, Initial residual = 0.858181339411, Final residual = 6.32297321249e-06, No Iterations 9 time step continuity errors : sum local = 2.37161385914e-11, global = -7.86308500369e-13, cumulative = 2.06361691352e-05 smoothSolver: Solving for omega, Initial residual = 5.2565171788e-07, Final residual = 3.49278374392e-09, No Iterations 1 smoothSolver: Solving for k, Initial residual = 2.66135012542e-06, Final residual = 7.33859698627e-11, No Iterations 2 ExecutionTime = 257.81 s ClockTime = 259 s Courant Number mean: 3.23080344243e-06 max: 4.63861488263 Interface Courant Number mean: 2.53529659665e-07 max: 4.63861488263 deltaT = 4.45346458291e-07 Time = 67.1303 PIMPLE: iteration 1 Updating StokesI wave model for patch inlet smoothSolver: Solving for alpha.water, Initial residual = 1.28408799391e-08, Final residual = 1.08364680021e-11, No Iterations 1 Phase-1 volume fraction = 0.666618261848 Min(alpha.water) = -4.66529651265e-10 Max(alpha.water) = 1.00014678621 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.666618261848 Min(alpha.water) = -4.66080169839e-10 Max(alpha.water) = 1.00014678621 Updating shallowWaterAbsorption wave model for patch outlet GAMGPCG: Solving for p_rgh, Initial residual = 0.876222984413, Final residual = 4.82480308588e-06, No Iterations 9 GAMGPCG: Solving for p_rgh, Initial residual = 0.76986420467, Final residual = 5.12962944087e-06, No Iterations 9 time step continuity errors : sum local = 4.52648844462e-12, global = -1.38478426638e-13, cumulative = 2.06361689967e-05 GAMGPCG: Solving for p_rgh, Initial residual = 0.845759428026, Final residual = 5.01065970515e-06, No Iterations 9 GAMGPCG: Solving for p_rgh, Initial residual = 0.874705259624, Final residual = 4.63805995333e-06, No Iterations 9 time step continuity errors : sum local = 1.00680418755e-11, global = -1.63755700805e-13, cumulative = 2.0636168833e-05 smoothSolver: Solving for omega, Initial residual = 2.70184710423e-07, Final residual = 2.50452634116e-10, No Iterations 1 smoothSolver: Solving for k, Initial residual = 1.28659289894e-06, Final residual = 1.3294530894e-09, No Iterations 1 ExecutionTime = 282.65 s ClockTime = 284 s -
interFoam短周期波浪衰减过于严重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，但依旧迭代不出来，请问您有什么解决办法吗
-
SHM尖锐直角边界层添加环形浮体的边角周围做了加了处理，之前以为是过渡不均匀造成的，但改善效果不明显
浮体外半径12.5m，内半径7.5m 高度12.5m 向上提升5m做自衰减 -
SHM尖锐直角边界层添加最近在使用OF计算环形浮体的自衰减，采用k-omega SST模型，网格基于SHM离散，但是边角处的边界层始终贴合不上，请大神支支招。
具体问题如下：
-
边界贴合质量：尝试通过提升stl质量（SW和犀牛加密）/增加surfaceFeatureExtractDict-includedAngle/提升snappyHexMeshDict-features和refinementSurfaces level均没能很好地将边角贴合地很好。远看质量还行，但是边角网格质量细看却很差，请问有何提升方法
-
基于k-omega SST计算的yplus在1000左右，首层边界层设置为1.75mm，发现是在内环面拐角处边界层未添加上，尝试好多办法均无法添加，大家有没有添加的办法
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2306 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object snappyHexMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // castellatedMesh true; snap true; addLayers true; geometry { WECMiddle.stl { type triSurfaceMesh; name WECMiddle; regions { OBJECT // Defined in .stl file { name WECMiddle; } } } WECMiddleWaveElevationRefineZone { type searchableBox; min (-17.5 -17.5 -1.25); max (17.5 17.5 11.25); } WECMiddleCornerRefineZone1 { type searchableCone; point1 (0 0 -1.5); radius1 12.75; innerRadius1 7.25; point2 (0 0 -1); radius2 12.75; innerRadius2 7.25; } WECMiddleCornerRefineZone2 { type searchableCone; point1 (0 0 11); radius1 12.75; innerRadius1 7.25; point2 (0 0 11.5); radius2 12.75; innerRadius2 7.25; } }; castellatedMeshControls { maxLocalCells 200000000; maxGlobalCells 300000000; minRefinementCells 20; nCellsBetweenLevels 3; maxLoadUnbalance 0.1; allowFreeStandingZoneFaces true; resolveFeatureAngle 30; features ( { file "WECMiddle.eMesh"; level 5; } ); refinementSurfaces { WECMiddle { level (3 3); // 1.25/2/2^3 0.0781 patchInfo { type wall; } regions { WECMiddle { level (3 3); patchInfo { type wall; } } } } } refinementRegions { WECMiddleWaveElevationRefineZone { mode inside; levels ((2 2)); } WECMiddleCornerRefineZone1 { mode inside; levels ((4 4)); } WECMiddleCornerRefineZone2 { mode inside; levels ((4 4)); } } locationInMesh (0 0 0); // Generate meshes outside the body } snapControls { tolerance 2; implicitFeatureSnap false; explicitFeatureSnap true; multiRegionFeatureSnap true; detectNearSurfacesSnap true; nSmoothPatch 5; nSolveIter 300; nRelaxIter 8; nFeatureSnapIter 10; nSmoothInternal 5; nFaceSplitInterval 5; } addLayersControls { layers { WECMiddle { nSurfaceLayers 15; } } relativeSizes false; //thickness 0.1; firstLayerThickness 0.00175; //finalLayerThickness 0.1; expansionRatio 1.25; featureAngle 360; slipFeatureAngle 0; nGrow 0; nBufferCellsNoExtrude 0; minMedialAxisAngle 90; maxFaceThicknessRatio 0.5; maxThicknessToMedialRatio 0.3; minThickness 1e-06; nLayerIter 60; nRelaxIter 5; nSmoothSurfaceNormals 10; nSmoothNormals 3; nSmoothThickness 10; nRelaxedIter 20; nMedialAxisIter 10; } meshQualityControls { minVol 1e-20; minTetQuality 1e-20; minArea 1e-20; minTwist 0.05; minDeterminant 1e-06; minFaceWeight 0.02; minVolRatio 0.01; minTriangleTwist -1; minFlatness 0.5; maxNonOrtho 70; maxBoundarySkewness 20; maxInternalSkewness 4; maxConcave 80; nSmoothScale 4; errorReduction 0.75; relaxed { minVol 1e-30; minTetQuality 1e-30; minArea 1e-30; minTwist 0.001; minDeterminant 1e-06; minFaceWeight 1e-06; minVolRatio 0.01; minTriangleTwist -1; minFlatness 0.5; maxNonOrtho 75; maxBoundarySkewness 20; maxInternalSkewness 8; maxConcave 80; nSmoothScale 4; errorReduction 0.75; } } mergeTolerance 1e-08; debug 0; // ************************************************************************* //
- 首层边界层目前为1.75mm是调试很多次后得到的最小厚度，但仍然不满足要求，请问大家有没有进一步降低的方法
先谢谢各位了
-
-
snappyHexMesh定义了多域网格（MultiRegion）失败。@jpzhao 在 snappyHexMesh定义了多域网格（MultiRegion）失败。 中说：
@李子橙 所以你的意思是如果是多个域计算传热，就每个域单独画网格，然后再将这几个网格合并?
请问这个问题解决了吗 mergeMeshes的关键问题是子域网格的边界如何设定或者如何取消呢？
-
在mergeMesh后删除边界各位大佬，最近想用OF计算固定的多体装置，由于在大域内生成的结构过于粗糙 需在小域里生成 然后用mergeMesh合并到大域里去，但是小域的边界没办法删除
请问各位大佬，有没有什么设置建议？删除边界或指定特殊的边界？ -
如何在mergeMesh后删除边界各位大佬，最近想用OF计算固定的多体装置，由于在大域内生成的结构过于粗糙 需在小域里生成 然后用mergeMesh合并到大域里去，但是小域的边界没办法删除
请问各位大佬，有没有什么设置建议？删除边界或指定特殊的边界？ -
rigidBodyMotion和sixDoFRigidBodyMotion的区别各位大佬，本人刚入坑不久，想用OF计算浮体运动，请问大家sixDoFRigidBodyMotion和rigidBodyMotion两个求解器的区别究竟是什么 翻找好多资料似乎都没有讲清楚这个事
从功能来看，rigidBodyMotion似乎支持更多约束类型，其他似乎没有区别了\multiphase\interFoam\RAS\floatingObject算例提供了sixDoFRigidBodyMotion和rigidBodyMotion两种方法计算，似乎也没什么差别
那么计算浮体在波浪中的运动，应该选用哪种求解器呢？ -
Wave2FOAM 和 OpenFOAM自带造波各位大佬 最近刚入手OF 发现OF里自带了造波和消波功能，想请教一下OF自带能否消除浮体和造波板之间的反射波？我看到/tutorials/multiphase\interFoam/laminar/waves/mangroveInteraction/constant/waveProperties文件启用了activeAbsorption yes选项，该选项是否可用于消除消除浮体和造波板之间的反射波
另外 注意到Wave2FOAM大概已经一年没有更新了 是不是因为OF的造波已经很强大了
小白发言 祝各位大佬paper全中~~~~