[gerritgroot]

Hi,

I'm completely new to OpenFoam and I got stuck with fvSchemes.

I'm trying to run a rotational periodic case with SRFSimpleFoam. Because of bad convergence, I'd like to start with first order schemes
However, I'm having problems with the scheme specs of "div((nuEff*dev(T(grad(Urel)))))" which only seems to work with Gauss linear;

Furthermore I have problems running the laplacian schemes.

I specified the fvSchemes file as:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.0.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default         steadyState;
}

gradSchemes
{
    default         none;
    grad(p)         Gauss linear;
    grad(Urel)      Gauss linear;
    grad(omega)     Gauss linear;
    grad(k)         Gauss linear;
}

divSchemes
{
    default                         none;
    div(phi,Urel)                   bounded Gauss upwind; // bounded Gauss linear;
    div(phi,k)                      bounded Gauss upwind; // bounded Gauss linear;
    div(phi,epsilon)                bounded Gauss upwind; // bounded Gauss linear;
    div(phi,omega)                  bounded Gauss upwind; // bounded Gauss linear;
    div(phi,R)                      bounded Gauss upwind; // bounded Gauss linear;
    div(R)                          bounded Gauss upwind; // bounded Gauss linear;
    div(phi,nuTilda)                bounded Gauss upwind; // bounded Gauss linear;
    div((nuEff*dev(T(grad(Urel))))) Gauss linear; // bounded Gauss upwind; did not work
}

laplacianSchemes
{
    default                        none;
    laplacian(nuEff,Urel)          Gauss linear corrected; //Gauss upwind uncorrected; did not work
    laplacian((1|A(Urel)),p)       Gauss linear corrected; //Gauss upwind uncorrected; did not work
    laplacian(DkEff,k)             Gauss linear corrected; //Gauss upwind uncorrected; did not work
    laplacian(DepsilonEff,epsilon) Gauss linear corrected; //Gauss upwind uncorrected; did not work
    laplacian(DomegaEff,omega)     Gauss linear corrected; //Gauss upwind uncorrected; did not work
    laplacian(DREff,R)             Gauss linear corrected; //Gauss upwind uncorrected; did not work
    laplacian(DnuTildaEff,nuTilda) Gauss linear corrected; //Gauss upwind uncorrected; did not work
}

interpolationSchemes
{
    default           linear;
    interpolate(Urel) linear;
}

snGradSchemes
{
    default         corrected;
}

fluxRequired
{
    default         no;
    p               ;
}


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

div((nuEff*dev(T(grad(Urel))))) Gauss linear;

Works fine, but

div((nuEff*dev(T(grad(Urel))))) bounded Gauss upwind;

Gives the error:

[2]
[2]
[2] --> FOAM FATAL IO ERROR:
[2] unknown div scheme bounded

Valid div schemes are :
1(Gauss)
[2]
[2] file: IOstream.divSchemes.div((nuEff*dev(T(grad(Urel))))) at line 0.
[2]
[2]     From function divScheme<Type>::New(const fvMesh&, Istream&)
[2]     in file /opt/OpenFOAM/OpenFOAM-2.3.x/src/finiteVolume/lnInclude/divScheme.C at line 82.
[2]
FOAM parallel run exiting
[2]

While: div((nuEff*dev(T(grad(Urel))))) Gauss upwind;

Gives the error:

[2]
[2]
[2] --> FOAM FATAL IO ERROR:
[2] attempt to read beyond EOF
[2]
[2] file: IOstream.divSchemes.div((nuEff*dev(T(grad(Urel))))) at line 0.
[2]
[2]     From function ITstream::read(token&)
[2]     in file db/IOstreams/Tstreams/ITstream.C at line 83.
[2]
FOAM parallel run exiting
[2]

Furthermore, while trying to run the first laplacian(nuEff,Urel) with Gauss upwind uncorrected;

Gives the error:

[2]
[2]
[2] --> FOAM FATAL ERROR:
[2]
    request for surfaceScalarField uncorrected from objectRegistry region0 failed
    available objects of type surfaceScalarField are
1(phi)
[2]
[2]     From function objectRegistry::lookupObject<Type>(const word&) const
[2]     in file /opt/OpenFOAM/OpenFOAM-2.3.x/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line
[0]
198

.
[0]


[2]
[0]

FOAM parallel run aborting
--> FOAM FATAL ERROR:
[2]


[0]
    request for surfaceScalarField uncorrected from objectRegistry region0 failed
    available objects of type surfaceScalarField are
1(phi)
[0]
[0]     From function objectRegistry::lookupObject<Type>(const word&) const
[0]     in file /opt/OpenFOAM/OpenFOAM-2.3.x/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 198.
[0]
FOAM parallel run aborting
[0]

Anyone any idea? How do I set everything to first order to avoid bad convergence in the first iterations?
Thanks in advance!