pimpleFoam with komegaSST results in error

drink water

The error is as follows. However, if I input command "unset FOAM_SIGFPE" before "pimpleFoam", everything is OK. the OpenFOAM could run as usual. Could anybody help me for this stupid question?

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 1


PIMPLE: No convergence criteria found


PIMPLE: Operating solver in transient mode with 1 outer corrector
PIMPLE: Operating solver in PISO mode


Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kOmegaSST
Selecting patchDistMethod meshWave
bounding k, min: 0 max: 0.00375 average: 0.00375
bounding omega, min: 0 max: 0.559 average: 0.559
RAS
{
    model           kOmegaSST;
    turbulence      on;
    printCoeffs     on;
    alphaK1         0.85;
    alphaK2         1;
    alphaOmega1     0.5;
    alphaOmega2     0.856;
    gamma1          0.555556;
    gamma2          0.44;
    beta1           0.075;
    beta2           0.0828;
    betaStar        0.09;
    a1              0.31;
    b1              1;
    c1              10;
    F3              false;
}

No MRF models present

No finite volume options present
#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::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
#5  Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel<Foam::IncompressibleMomentumTransportModel<Foam::transportModel> > >, Foam::IncompressibleMomentumTransportModel<Foam::transportModel> >::F2() const at ??:?
#6  Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel<Foam::IncompressibleMomentumTransportModel<Foam::transportModel> > >, Foam::IncompressibleMomentumTransportModel<Foam::transportModel> >::F23() const at ??:?
#7  Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel<Foam::IncompressibleMomentumTransportModel<Foam::transportModel> > >, Foam::IncompressibleMomentumTransportModel<Foam::transportModel> >::correctNut() at ??:?
#8  ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
#9  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10  ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
Floating point exception (core dumped)

李东岳

The minimum value of your $k$ and $\omega$ is 0, which looks problematic to me. Double check your boundary conditions. kOmega is not as stable as kEpsilon model, you need to be careful. kOmegaSST model is even worse.

drink water

@李东岳谢谢李老师！我估计是边界条件的事，有些边界条件用的fixedValue，0，但是我改为10e-6后仍旧报错，我再试试。ps：帖子用的ubuntu系统，没中文输入法，无奈只能用英文发帖了

drink water

@李东岳确实如李老师所说，将固定壁面边界条件Omega的值设置为1e-6后，能算了！！！

CFD中文网

pimpleFoam with komegaSST results in error