[shock77]

Hi,


I am trying to use the kOmegaSSTSAS model with rhoCentralFoam. Unfortunaletely I have some troubles starting a simulation with it. kOmegaSST works fine.


When I define my turbulence model like:

simulationType  RAS;

RAS
{
    RASModel        kOmegaSSTSAS;

    turbulence      on;

    printCoeffs     on;


kOmegaSSTSASCoeffs
    {
        // Default SST coefficients
        alphaK1     0.85;
        alphaK2     1.0;
        alphaOmega1 0.5;
        alphaOmega2 0.856;
        beta1       0.075;
        beta2       0.0828;
        betaStar    0.09;
        gamma1      5/9;
        gamma2      0.44;
        a1          0.31;
        b1          1.0;
        c1          10.0;
        F3          no;

        // Default SAS coefficients
        Cs          0.11;
        kappa       0.41;
        zeta2       3.51;
        sigmaPhi    2.0/3.0;
        C           2;


        delta cubeRootVol;

        cubeRootVolCoeffs
        {
         deltaCoeff      1;
    }
    }
}

I get the error Message:

--> FOAM FATAL IO ERROR: 
keyword delta is undefined in dictionary "/home/case_sas/constant/turbulenceProperties.RAS.kOmegaSSTCoeffs"

Which is weird, because thats the way its definied in the sourcecode? So I switch kOmegaSSTSASCoeffs to kOmegaSSTCoeffs and the simulation starts, but crashes after the first iteration:

Mean and max Courant Numbers = 0.01507015757285037026 0.1734610019636567413
deltaT = 1.172882946281961179e-08
Time = 1.172882946281961179e-08

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 0.9983413068152621195, Final residual = 4.523050796451504211e-09, No Iterations 1000
smoothSolver:  Solving for Uy, Initial residual = 0.999903823391497748, Final residual = 2.612725370460363529e-10, No Iterations 4
smoothSolver:  Solving for Uz, Initial residual = 0.9982737144596987022, Final residual = 4.656112886042951805e-09, No Iterations 1000
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for e, Initial residual = 0.1588687188325532407, Final residual = 2.735319357247494125e-09, No Iterations 1000
#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::sqrt(Foam::Field<double>&, Foam::UList<double> const&) at ??:?
#4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
#5  Foam::RASModels::kOmegaSSTSAS<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > >::Qsas(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const at ??:?
#6  Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > > >, Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > >::correct() at ??:?
#7  ? at ??:?
#8  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9  ? at ??:?
Floating point exception (core dumped)

Its my first time running a hybrid model. I hope someone more expericenced can help.

[shock77]

This is how it is described in the .H file of the turbulence model.

    References:   34     \verbatim
   35         Egorov, Y., & Menter F.R. (2008).
   36         Development and Application of SST-SAS Model in the DESIDER Project.
   37         Advances in Hybrid RANS-LES Modelling,
   38         Notes on Num. Fluid Mech. And Multidisciplinary Design,
   39         Volume 97, 261-270.
   40     \endverbatim
   41 
   42     The model coefficients are
   43     \verbatim
   44         kOmegaSSTSASCoeffs
   45         {
   46             // Default SST coefficients
   47             alphaK1     0.85;
   48             alphaK2     1.0;
   49             alphaOmega1 0.5;
   50             alphaOmega2 0.856;
   51             beta1       0.075;
   52             beta2       0.0828;
   53             betaStar    0.09;
   54             gamma1      5/9;
   55             gamma2      0.44;
   56             a1          0.31;
   57             b1          1.0;
   58             c1          10.0;
   59             F3          no;
   60 
   61             // Default SAS coefficients
   62             Cs          0.11;
   63             kappa       0.41;
   64             zeta2       3.51;
   65             sigmaPhi    2.0/3.0;
   66             C           2;
   67 
   68             // Delta must be specified for SAS e.g.
   69             delta cubeRootVol;
   70 
   71             cubeRootVolCoeffs
   72             {}
   73         }
   74     \endverbatim
   75 
   76 SourceFiles
   77     kOmegaSSTSAS.C

But kOmegaSSTSASCoeffs is not accepted.

[Brodenson]

Hello!, I had a similar problem but trying to run with simpleFoam. I get a
[1] --> FOAM FATAL IO ERROR:
[1] keyword delta is undefined in dictionary "IOstream.RAS"