CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

First order in fvSchemes does not seem to work

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 30, 2015, 20:06
Question First order in fvSchemes does not seem to work
  #1
Member
 
Anonymous
Join Date: Jul 2009
Posts: 44
Rep Power: 15
gerritgroot is on a distinguished road
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:
Code:
/*--------------------------------*- 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:
Code:
[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:
Code:
[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:
Code:
[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!
gerritgroot is offline   Reply With Quote

Reply

Tags
fvschemes first order

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
simpleFoam: switch so second order fvSchemes okroud OpenFOAM Running, Solving & CFD 11 September 11, 2013 12:52
First order simulation better than second order Torque_Converter CFX 7 January 8, 2013 04:07
How to compute UDS fourth order gradient Emma66 FLUENT 1 April 27, 2010 05:44
OpenFOAM - 2nd order of discretization? makaveli_lcf OpenFOAM Running, Solving & CFD 0 November 12, 2009 07:15
Unstable flow simpleFoam 2nd order Valle OpenFOAM 0 August 26, 2009 08:12


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