[thomas_toulorge]

Hello,

I am trying to obtain a fully developed velocity profile in a duct at low Reynolds number (Re ~ 6000) by using an axisymmetric geometry with boundaryFoam. I would like to use the k-omega SST turbulence model, but I get a strange error that seems to occur at initialization of the turbulence model object:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.7.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 1.7.x-3776603e4c6c
Exec   : /opt/openfoam171/applications/bin/linux64GccDPOpt/boundaryFoam
Date   : May 12 2011
Time   : 18:34:16
Host   : PMA-08-003
PID    : 10174
Case   : /home/u0056407/DG/Tests/SideCavity/CFD/Duct/Duct_Test
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

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

Create mesh for time = 0

Reading field U

Creating face flux

Selecting incompressible transport model Newtonian
Selecting RAS turbulence model kOmegaSST


--> FOAM FATAL ERROR: 
Maximum number of iterations reached. Increase maxIter.
    maxIter:20
    nChangedCells:0
    nChangedFaces:4


    From function FaceCellWave<Type>::FaceCellWave(const polyMesh&, const labelList&, const List<Type>, UList<Type>&, UList<Type>&, const label maxIter)
    in file /home/opencfd/OpenFOAM/OpenFOAM-1.7.x/src/OpenFOAM/lnInclude/FaceCellWave.C at line 862.

FOAM exiting

I also get this error for other low-Re turbulence models (LamBremhorstKE, LienCubicKELowRe, LienLeschzinerLowRe), while the case runs when "RASModel" is set to "laminar". However, boundaryFoam works with the komegaSST model for a channel geometry as in the tutorial.

Am I doing something wrong? Is it even possible at all to use boundaryFoam in an axisymmetric setting? The files are attached.

Thanks in advance!

Thomas.