pimpleFoam with komegaSST results in error



  • 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.



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



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


Log in to reply
 


CFD中文网 | 东岳流体学术 | 东岳流体商业 | 吉ICP备20003622号-1