[TomWaits]

Hi FOAMers,

I am new to OpenFOAM and I am having an issue with running a case with kOmegaSST on v2006 - I don't have any issues with k-epsilon. I get the following cryptic error:

Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kOmegaSST
Selecting patchDistMethod meshWave
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in /lib64/libpthread.so.0
#3  Foam::patchWave::getValues(Foam::MeshWave<Foam::wallPoint, int> const&) at ??:?
#4  Foam::patchWave::correct() at ??:?
#5  Foam::patchWave::patchWave(Foam::polyMesh const&, Foam::HashSet<int, Foam::Hash<int> > const&, bool) at ??:?
#6  Foam::patchDistMethods::meshWave::correct(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&) at ??:?
#7  Foam::wallDist::wallDist(Foam::fvMesh const&, Foam::HashSet<int, Foam::Hash<int> > const&, Foam::word const&) at ??:?
#8  Foam::wallDist::wallDist(Foam::fvMesh const&, Foam::word const&) at ??:?
#9  Foam::wallDist const& Foam::MeshObject<Foam::fvMesh, Foam::UpdateableMeshObject, Foam::wallDist>::New<>(Foam::fvMesh const&) at ??:?
#10  Foam::kOmegaSSTBase<Foam::eddyViscosity<Foam::RASModel<Foam::IncompressibleTurbulenceModel<Foam::transportModel> > > >::kOmegaSSTBase(Foam::word const&, Foam::geometricOneField const&, Foam::geometricOneField const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel const&, Foam::word const&) at ??:?
#11  Foam::RASModels::kOmegaSST<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::kOmegaSST(Foam::geometricOneField const&, Foam::geometricOneField const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel const&, Foam::word const&, Foam::word const&) at ??:?
#12  Foam::RASModel<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::adddictionaryConstructorToTable<Foam::RASModels::kOmegaSST<Foam::IncompressibleTurbulenceModel<Foam::transportModel> > >::New(Foam::geometricOneField const&, Foam::geometricOneField const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel const&, Foam::word const&) at ??:?
#13  Foam::RASModel<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::New(Foam::geometricOneField const&, Foam::geometricOneField const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel const&, Foam::word const&) at ??:?
#14  Foam::TurbulenceModel<Foam::geometricOneField, Foam::geometricOneField, Foam::incompressibleTurbulenceModel, Foam::transportModel>::adddictionaryConstructorToTable<Foam::RASModel<Foam::IncompressibleTurbulenceModel<Foam::transportModel> > >::NewTurbulenceModel(Foam::geometricOneField const&, Foam::geometricOneField const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel const&, Foam::word const&) at ??:?
#15  ? at ??:?
#16  ? at ??:?
#17  __libc_start_main in /lib64/libc.so.6
#18  ? at ??:?
Floating point exception (core dumped)

Here is my turbulenceProperties file:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2006                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      turbulenceProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

simulationType RAS;

RAS
{
    RASModel            kOmegaSST;

    turbulence          on;

    printCoeffs         on;
}

// ************************************************************************* //

my 0/k:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1812                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

kInlet          0.015; // approx k = 1.5*(I*U)^2

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform $kInlet;

boundaryField
{
    inlet
    {
        type            turbulentIntensityKineticEnergyInlet;
        intensity       0.1;
        value           uniform $kInlet;
    }

    outlet
    {
        type            inletOutlet;
        inletValue      uniform $kInlet;
        value           uniform $kInlet;
    }

    slip1
    {
        type            slip;
    }

    slip2
    {
        type            slip;
    }

    symmetry1
    {
        type            symmetryPlane;
    }

    symmetry2
    {
        type            symmetryPlane;
    }
}

// ************************************************************************* //

my 0/nut

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1812                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -1 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    inlet
    {
        type            calculated;
        value           uniform 0;
    }

    outlet
    {
        type            calculated;
        value           uniform 0;
    }

    slip1
    {
        type            slip;
    }

    slip2
    {
        type            slip;
    }

    symmetry1
    {
        type            symmetryPlane;
    }

    symmetry2
    {
        type            symmetryPlane;
    }

}

// ************************************************************************* //

and my 0/omega:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1812                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      omega;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

omegaInlet      2.23607; // omega = Cu^-0.25*sqrt(k)/L , where Cu \approx 0.09

dimensions      [0 0 -1 0 0 0 0];

internalField   uniform $omegaInlet;

boundaryField
{
    inlet
    {
        type            turbulentMixingLengthFrequencyInlet;
        mixingLength    0.1;
        value           uniform $omegaInlet;
    }

    outlet
    {
        type            inletOutlet;
        inletValue      uniform $omegaInlet;
        value           uniform $omegaInlet;
    }

    slip1
    {
        type            slip;
    }

    slip2
    {
        type            slip;
    }

    symmetry1
    {
        type            symmetryPlane;
    }
    
    symmetry2
    {
        type            symmetryPlane;
    }
}

// ************************************************************************* //

Could anybody please give me some direction as to what the source of the error may be? I would be grateful for any help.

Many thanks in advance,
Tom Waits

[tomf]

Hi Tom,

Not sure if this would be the reason, but it might be that if no value is being set explicitly a value of 0 is used. For omega this can result in a division by 0, which would blow up. As such it is maybe best to specify also a uniform value for the slip and symmetry boundaries:

    slip2
    {
        type            slip;
        value           uniform $omegaInlet;
    }

    symmetry1
    {
        type            symmetryPlane;
        value           uniform $omegaInlet;
    }

and similar for k and the other boundaries. It should not be necessary for nut, but no problems if you specify it there either.

If this does not help it is probably somewhere else, but I have no suggestion where to look.

Cheers,
Tom