CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   boundaryFoam, axisymmetry and turbulence models (http://www.cfd-online.com/Forums/openfoam-solving/88279-boundaryfoam-axisymmetry-turbulence-models.html)

thomas_toulorge May 12, 2011 13:05

boundaryFoam, axisymmetry and turbulence models
 
1 Attachment(s)
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:

Code:

/*---------------------------------------------------------------------------*\
| =========                |                                                |
| \\      /  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.


All times are GMT -4. The time now is 14:52.