为什么二维方腔流用pisFoam算出的流线与下方没有涡

[郑澄远]

该怎么解决呢

[dzw05]

你只是在下方没画出来流线，应该是有涡的。

[李东岳]

有可能，你加密试试？不行就降低粘度试试？

[郑澄远]

@东岳 把网格加密后80*80算到time=0.05时就显示

#3  Foam::GaussSeidelSmoother::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::GaussSeidelSmoother::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  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#8  Foam::fvMatrix<double>::solve() at ??:?
#9  Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<double> > const&) at ??:?
#10  Foam::RASModels::kEpsilon<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::correct() at ??:?
#11  ? at ??:?
#12  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13  ? at ??:?

浮点数例外 (核心已转储)

[李东岳]

降低你的时间步长试试

[肖艺]

@东岳
二维顶盖驱动方腔流无滑移时有涡，有滑移时不出现涡，这种情况可能是什么原因造成的？程序是fortran语言自编的