[jaason]

Hi FOAMers,
I want to implement a 2nd order upwind scheme. However i get the error message below whenever I try to use a scheme like QUICK or linearUpwind for my InterpolationSchemes, gradSchemes or my LaplacianSchemes. I can only implement a second order upwind scheme in the divergence section. My fvSchemes file is pasted below

attempt to read beyond EOF

file: /home/owner/OpenFOAM/owner-2.2.0/run/Thesis/sonicFoam/Sept12/system/fvSchemes.interpolationSchemes.interpolate(U_0) at line 59.

    From function ITstream::read(token&)
    in file db/IOstreams/Tstreams/ITstream.C at line 83.

FOAM exiting

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         none;
    grad(U)        Gauss upwind phi;
    grad(p)        Gauss upwind phi;
    grad(e)        Gauss upwind phi;
}

divSchemes
{
    default         none;
    div(phi,U)      Gauss upwind phi;
    div(phid,p)     Gauss upwind phi;
    div(phi,e)      Gauss upwind phi;
    div(phi,K)      Gauss upwind phi;
    div(phiv,p)     Gauss upwind phi;
    div((muEff*dev2(T(grad(U))))) Gauss linear 1;
}

laplacianSchemes
{
    default                 none;
    laplacian(muEff,U)        Gauss linear corrected;
    laplacian(alphaEff,e)    Gauss linear corrected;
    laplacian(Dp,p)        Gauss linear corrected;
}

interpolationSchemes
{
    default                         none;
    interpolate(rho)                 upwind phi;
    interpolate(rho_0)                 upwind phi;
    interpolate((((1|A(U))*rho_0)*U_0))        upwind phi;
    interpolate(HbyA)                upwind phi;
    interpolate(((1|A(U))*rho_0))            upwind phi;
    interpolate(thermo:psi)            upwind phi;
    interpolate(U_0)                 QUICK;
}

snGradSchemes
{
    default         corrected;
}

fluxRequired
{
    default         no;
    p               ;
}

what I want to do is substitute linearUpwind for all upwind terms in the code as to make it second order and therefore have to deal with less diffusion type inaccuracies. any suggestions??
Am I going about implementing a second order upwind scheme incorrectly?

I am using sonicFoam as a solver but I had the same problems with rhoCentralFoam

any help would be greatly appreciated.

[Yahoo]

Hey,
Were you able to solve this error. I am getting the same one?

Originally Posted by jaason

Hi FOAMers,
I want to implement a 2nd order upwind scheme. However i get the error message below whenever I try to use a scheme like QUICK or linearUpwind for my InterpolationSchemes, gradSchemes or my LaplacianSchemes. I can only implement a second order upwind scheme in the divergence section. My fvSchemes file is pasted below

attempt to read beyond EOF

file: /home/owner/OpenFOAM/owner-2.2.0/run/Thesis/sonicFoam/Sept12/system/fvSchemes.interpolationSchemes.interpolate(U_0) at line 59.

    From function ITstream::read(token&)
    in file db/IOstreams/Tstreams/ITstream.C at line 83.

FOAM exiting

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         none;
    grad(U)        Gauss upwind phi;
    grad(p)        Gauss upwind phi;
    grad(e)        Gauss upwind phi;
}

divSchemes
{
    default         none;
    div(phi,U)      Gauss upwind phi;
    div(phid,p)     Gauss upwind phi;
    div(phi,e)      Gauss upwind phi;
    div(phi,K)      Gauss upwind phi;
    div(phiv,p)     Gauss upwind phi;
    div((muEff*dev2(T(grad(U))))) Gauss linear 1;
}

laplacianSchemes
{
    default                 none;
    laplacian(muEff,U)        Gauss linear corrected;
    laplacian(alphaEff,e)    Gauss linear corrected;
    laplacian(Dp,p)        Gauss linear corrected;
}

interpolationSchemes
{
    default                         none;
    interpolate(rho)                 upwind phi;
    interpolate(rho_0)                 upwind phi;
    interpolate((((1|A(U))*rho_0)*U_0))        upwind phi;
    interpolate(HbyA)                upwind phi;
    interpolate(((1|A(U))*rho_0))            upwind phi;
    interpolate(thermo:psi)            upwind phi;
    interpolate(U_0)                 QUICK;
}

snGradSchemes
{
    default         corrected;
}

fluxRequired
{
    default         no;
    p               ;
}

what I want to do is substitute linearUpwind for all upwind terms in the code as to make it second order and therefore have to deal with less diffusion type inaccuracies. any suggestions??
Am I going about implementing a second order upwind scheme incorrectly?

I am using sonicFoam as a solver but I had the same problems with rhoCentralFoam

any help would be greatly appreciated.

[Bernhard]

Without looking I would say you have to supply "QUICK phi" instead of just "QUICK", because you need to tell the solver with respect to which variable it has to take the upwind values. This is a feature. The error message shown it the one you typically get when you have too little arguments in a line in a sub-dictionary.

[ron_OFuser]

Hello everyone,

I am a newbie..
What is the difference between the following categories.

(1)

divSchemes
{
default none;
div(phi,U) Gauss upwind phi;
div(phid,p) Gauss upwind phi;
div(phi,e) Gauss upwind phi;
div(phi,K) Gauss upwind phi;
div(phiv,p) Gauss upwind phi;
div((muEff*dev2(T(grad(U))))) Gauss linear 1;
}

(2)

interpolationSchemes
{
default none;
interpolate(rho) upwind phi;
interpolate(rho_0) upwind phi;
interpolate((((1|A(U))*rho_0)*U_0)) upwind phi;
interpolate(HbyA) upwind phi;
interpolate(((1|A(U))*rho_0)) upwind phi;
interpolate(thermosi) upwind phi;
interpolate(U_0) QUICK phi;
}

Any help will be greatly appreciated.

Thank you

[sharonyue]

Originally Posted by ron_OFuser

Hello everyone,

I am a newbie..
What is the difference between the following categories.

(1)

divSchemes
{
default none;
div(phi,U) Gauss upwind phi;
div(phid,p) Gauss upwind phi;
div(phi,e) Gauss upwind phi;
div(phi,K) Gauss upwind phi;
div(phiv,p) Gauss upwind phi;
div((muEff*dev2(T(grad(U))))) Gauss linear 1;
}

(2)

interpolationSchemes
{
default none;
interpolate(rho) upwind phi;
interpolate(rho_0) upwind phi;
interpolate((((1|A(U))*rho_0)*U_0)) upwind phi;
interpolate(HbyA) upwind phi;
interpolate(((1|A(U))*rho_0)) upwind phi;
interpolate(thermosi) upwind phi;
interpolate(U_0) QUICK phi;
}

Any help will be greatly appreciated.

Thank you

Hello,

the entry in interpolationSchemes is for some terms which need to get values from cell centres to faces. It includes all the terms in it.
In divSchemes, it only include the schemes for convections.

eg,

interpolate(HbyA) upwind phi;

This will call upwind scheme when u wanna get HbyA's face value from the centers. U can only define it here instead of defining it in divSchemes, cause HbyA is not a convection term

Best,