@东岳 多谢李老师了，是时间步长的原因，我开始按库伦数取得 听了您的建议往小了试，没有出现发散

感谢回复 李老师 但应该不是时间步长的问题。无论时间步长取得多短，总会在第五步出现发散

Time = 6e-05

Courant Number mean: 1.96719e+25 max: 9.01107e+27
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
#4  Foam::symGaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
#5  Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6  ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/icoFoam"
#7  ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/icoFoam"
#8  ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/icoFoam"
#9  ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/icoFoam"
#10  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11  ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/icoFoam"

使用icoFoam求解的时候在第五步出现发散，用paraview查看时发现第五步出现涡流。速度，粘度均很小，基本可以排除是湍流。请问这种现象是求解器使用不当还是其他原因，如果是的话，请问出现涡流的情况该使用什么求解器。