圆柱绕流 高雷诺数10e5
-
圆柱绕流在低雷诺数（600左右）的时候使用icoFoam能正常计算（nu=0.08），把nu改为10e-5后，加入了湍流模型，用pisoFoam 求解，但在计算一开始就会直接爆掉，最大库朗数1600多，怎么改都算不了，不知道为啥，请大家指点一下？是网格的原因吗还是我的设置问题？
圆柱直径4m，入流10m/s，网格量大概190000。用的k-epsilon，时间步长原本用的是0.005，后面改成0.0005了也不行，应该是设置问题，但不知道是哪。
贴一下信息
logCourant Number mean: 1.64478e+27 max: 6.7649463e+34 Time = 0.006 PIMPLE: iteration 1 smoothSolver: Solving for Ux, Initial residual = 0.99999999, Final residual = 6.1529517e-10, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 0.99731182, Final residual = 1.7252366e-11, No Iterations 1 [2] #0 Foam::error::printStack(Foam::Ostream&) at ??:? [2] #1 Foam::sigFpe::sigHandler(int) at ??:? [2] #2 ? in /lib/x86_64-linux-gnu/libpthread.so.0 [2] #3 Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:? [2] #4 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:? [2] #5 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? [2] #6 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:? [2] #7 Foam::fvMatrix<double>::solveSegregatedOrCoupled(Foam::dictionary const&) at ??:? [2] #8 Foam::fvMesh::solve(Foam::fvMatrix<double>&, Foam::dictionary const&) const at ??:? [2] #9 ? in ~/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/pimpleFoam [2] #10 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6 [2] #11 ? in ~/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/pimpleFoam [hipeson-X12DPi-N-T-6:358676] *** Process received signal *** [hipeson-X12DPi-N-T-6:358676] Signal: Floating point exception (8) [hipeson-X12DPi-N-T-6:358676] Signal code: (-6) [hipeson-X12DPi-N-T-6:358676] Failing at address: 0x3e800057914 [hipeson-X12DPi-N-T-6:358676] [ 0] /lib/x86_64-linux-gnu/libpthread.so.0(+0x14420)[0x7f9c482c6420] [hipeson-X12DPi-N-T-6:358676] [ 1] /lib/x86_64-linux-gnu/libpthread.so.0(raise+0xcb)[0x7f9c482c62ab] [hipeson-X12DPi-N-T-6:358676] [ 2] /lib/x86_64-linux-gnu/libpthread.so.0(+0x14420)[0x7f9c482c6420] [hipeson-X12DPi-N-T-6:358676] [ 3] /home/hipeson/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5scaleERNS_5FieldIdEES3_RKNS_9lduMatrixERKNS_10FieldFieldIS1_dEERKNS_8UPtrListIKNS_17lduInterfaceFieldEEERKS2_h+0x242)[0x7f9c48dc4792] [hipeson-X12DPi-N-T-6:358676] [ 4] /home/hipeson/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver6VcycleERKNS_7PtrListINS_9lduMatrix8smootherEEERNS_5FieldIdEERKS8_S9_S9_S9_S9_S9_RNS1_IS8_EESD_h+0x8ea)[0x7f9c48dc634a] [hipeson-X12DPi-N-T-6:358676] [ 5] /home/hipeson/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5solveERNS_5FieldIdEERKS2_h+0x668)[0x7f9c48dc8658] [hipeson-X12DPi-N-T-6:358676] [ 6] /home/hipeson/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x625)[0x7f9c4c7200d5] [hipeson-X12DPi-N-T-6:358676] [ 7] /home/hipeson/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE24solveSegregatedOrCoupledERKNS_10dictionaryE+0x498)[0x7f9c4bea8728] [hipeson-X12DPi-N-T-6:358676] [ 8] /home/hipeson/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZNK4Foam6fvMesh5solveERNS_8fvMatrixIdEERKNS_10dictionaryE+0x28)[0x7f9c4be412e8] [hipeson-X12DPi-N-T-6:358676] [ 9] pimpleFoam(+0x2fc8c)[0x55fd9630cc8c] [hipeson-X12DPi-N-T-6:358676] [10] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf3)[0x7f9c480e4083] [hipeson-X12DPi-N-T-6:358676] [11] pimpleFoam(+0x32aae)[0x55fd9630faae] [hipeson-X12DPi-N-T-6:358676] *** End of error message *** -------------------------------------------------------------------------- Primary job terminated normally, but 1 process returned a non-zero exit code. Per user-direction, the job has been aborted. -------------------------------------------------------------------------- -------------------------------------------------------------------------- mpirun noticed that process rank 2 with PID 0 on node hipeson-X12DPi-N-T-6 exited on signal 8 (Floating point exception). --------------------------------------------------------------------------U
dimensions [0 1 -1 0 0 0 0]; internalField uniform (9 0 0); boundaryField { INLET { type fixedValue; value uniform (9 0 0); } OUTLET { type zeroGradient; } WALL { type slip; } CYLINDER { type fixedValue; value uniform (0 0 0); } frontAndBackPlanes { type empty; } }P
dimensions [0 2 -2 0 0 0 0]; internalField uniform 0; boundaryField { INLET { type zeroGradient; } OUTLET { type fixedValue; value uniform 0; } WALL { type zeroGradient; } CYLINDER { type zeroGradient; } frontAndBackPlanes { type empty; } }k
dimensions [0 2 -2 0 0 0 0]; internalField uniform 0.375; boundaryField { INLET { type fixedValue; value uniform 0.375; } OUTLET { type zeroGradient; } WALL { type kqRWallFunction; value uniform 0.375; } CYLINDER { type kqRWallFunction; value uniform 0.375; } frontAndBackPlanes { type empty; } }nut
dimensions [0 2 -1 0 0 0 0]; internalField uniform 0; boundaryField { INLET { type calculated; value uniform 0; } OUTLET { type calculated; value uniform 0; } WALL { type nutkWallFunction; value uniform 0; } CYLINDER { type nutkWallFunction; value uniform 0; } frontAndBackPlanes { type empty; } }epsilon
dimensions [0 2 -3 0 0 0 0]; internalField uniform 0.07; boundaryField { INLET { type fixedValue; value uniform 0.07; } OUTLET { type zeroGradient; } WALL { type epsilonWallFunction; value uniform 0.07; } CYLINDER { type epsilonWallFunction; value uniform 0.07; } frontAndBackPlanes { type empty; } }感谢大家帮助！
-
谢谢，我试了一下，还是不行，我网格用ICEM画的，网格质量很高，就是不知道哪出的问题。还请大家指教一下
fvsolutionsolvers { p { solver PCG; preconditioner DIC; tolerance 1e-06; relTol 0.05; } pFinal { $p; relTol 0; } "(U|k|epsilon)" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-05; relTol 0.1; } "(U|k|epsilon)Final" { $U; relTol 0; } } PIMPLE { nNonOrthogonalCorrectors 0; nCorrectors 2; }fvscheme
ddtSchemes { default backWard; } gradSchemes { default Gauss linear; } divSchemes { default none; div(phi,U) Gauss linearUpwind grad(U); turbulence Gauss upwind; div(phi,k) $turbulence; div(phi,epsilon) $turbulence; div(phi,nuTilda) $turbulence; div(phi,R) $turbulence; div(R) Gauss linear; div((nuEff*dev2(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear corrected; } interpolationSchemes { default linear; } snGradSchemes { default corrected; } -
-
你这个算例有几处不太合理的地方：
- 估计你的网格是边界层首层厚度很小的那种吧，那么你用 standard k epsilon 这个湍流模型是不合适的，建议使用 kOmegaSST。
- 抛开 k-epsilon 模型的合理性不谈，你这个 k epsilon nut 的 Inlet 边界条件设置也是不合理的。首先，你的入口速度是 9m/s，k 的 Inlet 值为 0.375，意味着入口湍流度为约 5.5%，这个值是偏大的。而 epsilon 的如何值为 0.07，意味着入口的 nut 值为 0.09*0.375^2/0.07=0.1808，超过流体粘度的 10000倍，这个是严重偏大的值。
-
@李东岳 李老师，因为不知道怎么设初始值合适，所以我当时是参照foam里面其他案例设置的nut，用的piso，这应该属于URANS吧，我看论坛以前讨论过RANS是用稳态模拟，URANS是瞬态模拟。
nut：
dimensions [0 2 -1 0 0 0 0]; internalField uniform 0; boundaryField { INLET { type calculated; value uniform 0; } OUTLET { type calculated; value uniform 0; } WALL { type nutkWallFunction; value uniform 0; } CYLINDER { type nutkWallFunction; value uniform 0; } frontAndBackPlanes { type empty; } }
