[Calmly]

Hello!

Let's say I want to set different schemes to all of the interpolation "values" and not to set it with the "default ..."

How can I find in which "values" I need to define schemes?

[Joshua14]

You should be able to get rid of the default entry and run the case via whatever your solver is (simpleFoam?). The case should crash and the log file will say it is missing scheme blah. Now you know what scheme to add. Do this till the case doesn't crash anymore.

[Calmly]

i am using the rhoCentralFoam solver and when I do as you said I get this:

--> FOAM FATAL IO ERROR:
[2] keyword flux(rhoU) is undefined in dictionary "IOstream.interpolationSchemes"
[2]
[2] file: IOstream.interpolationSchemes from line 0 to line 0.
[2]
[2] From function const Foam::entry& Foam::dictionary::lookupEntry(const Foam::word&, bool, bool) const
[3]
[3]
[3] --> FOAM FATAL IO ERROR:
[3] keyword flux(rhoU) is undefined in dictionary "IOstream.interpolationSchemes"
[3]
[3] file: IOstream.interpolationSchemes from line 0 to line 0.
[3]
[3] From function const Foam::entry& Foam::dictionary::lookupEntry(const Foam::word&, bool, bool) const
[3] in file [0]
[0]
[0] --> FOAM FATAL IO ERROR:
[0] keyword flux(rhoU) is undefined in dictionary "/home/themistoklis/OpenFOAM/themistoklis-4.1/VKI/2D/openField_ground_2D/Turbulent/ras/blast_2D_lU_1500_newMesh/system/fvSchemes.interpolationSchemes"
[0]
[0] file: /home/themistoklis/OpenFOAM/themistoklis-4.1/VKI/2D/openField_ground_2D/Turbulent/ras/blast_2D_lU_1500_newMesh/system/fvSchemes.interpolationSchemes from line 46 to line 50.
[0]
[0] From function const Foam::entry& Foam::dictionary::lookupEntry(const Foam::word&, bool, bool) const
[0] in file db/dictionary/dictionary.C at line 441.[2] in file db/dictionary/dictionary.C at line 441.
[2]
FOAM parallel run exiting
[2]

[0]
FOAM parallel run exiting
[0]
db/dictionary/dictionary.C at line 441.
[3]
FOAM parallel run exiting
[3]

[1]
[1] --> FOAM FATAL IO ERROR:
[1] keyword flux(rhoU) is undefined in dictionary "IOstream.interpolationSchemes"
[1]
[1] file: IOstream.interpolationSchemes from line 0 to line 0.
[1]
[1] From function const Foam::entry& Foam::dictionary::lookupEntry(const Foam::word&, bool, bool) const
[1] in file db/dictionary/dictionary.C at line 441.
[1]
FOAM parallel run exiting
[1]
[themistoklis-444:13538] 3 more processes have sent help message help-mpi-api.txt / mpi-abort

I don't understand how to define the flux(rhoU) in the interpolation schemes...

I have managed to define the schemes I want to some other values like rho, U, T but this...I can't..

[Joshua14]

What does your fvSchemes file look like? I'm guessing you should have something like the following in that file.

gradSchemes
{
.
.
.
}

divSchemes
{
.
.
.
}

interpolationSchemes
{
flux(rhoU) linear;
}

[Calmly]

Well, yeah that worked. Thanks!

But what if I want to use a different scheme like Minmod? It doesn't work.. I get this error if I put anything else except "linear"

[3] --> FOAM FATAL IO ERROR:
[3] attempt to read beyond EOF
[3]
[3] file: IOstream.interpolationSchemes.dotInterpolate(S,tau MC) at line 0.
[3]
[3] From function virtual Foam::Istream& Foam::ITstream::read(Foam::token&)
[3] in file db/IOstreams/Tstreams/ITstream.C at line 82.
[3]
FOAM parallel run exiting
[3]
[1]
[1]
[1] --> FOAM FATAL IO ERROR:
[1] attempt to read beyond EOF
[1]
[1] file: IOstream.interpolationSchemes.dotInterpolate(S,tau MC) at line 0.
[1]
[1] From function virtual Foam::Istream& Foam::ITstream::read(Foam::token&)
[1] in file db/IOstreams/Tstreams/ITstream.C at line 82.
[1]
FOAM parallel run exiting
[1]

[2]
[2] --> FOAM FATAL IO ERROR:
[2] attempt to read beyond EOF
[2]
[2] file: IOstream.interpolationSchemes.dotInterpolate(S,tau MC) at line 0.
[2]
[2] From function virtual Foam::Istream& Foam::ITstream::read(Foam::token&)
[2] in file db/IOstreams/Tstreams/ITstream.C at line 82.
[2]
FOAM parallel run exiting
[2]
[0]
[0]
[0] --> FOAM FATAL IO ERROR:
[0] attempt to read beyond EOF
[0]
[0] file: /home/themistoklis/OpenFOAM/themistoklis-4.1/VKI/2D/openField_ground_2D/Turbulent/ras/blast_2D_lU_baloon_minmod/system/fvSchemes.interpolationSchemes.dotInterpolate(S,ta uMC) at line 51.
[0]
[0] From function virtual Foam::Istream& Foam::ITstream::read(Foam::token&)
[0] in file db/IOstreams/Tstreams/ITstream.C at line 82.
[0]
FOAM parallel run exiting
[0]
[themistoklis-444:24232] 3 more processes have sent help message help-mpi-api.txt / mpi-abort
[themistoklis-444:24232] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages

The above error occurs when I replace the linear with the Minmod scheme.

But, like I told you it works with this:

interpolationSchemes
{
default none;
reconstruct(rho) Minmod;
reconstruct(U) MinmodV;
reconstruct(T) Minmod;
flux(rhoU) linear;
dotInterpolate(S,tauMC) linear;
interpolate(muEff) linear;
interpolate(rho) linear;
}

[Joshua14]

Unfortunately I don't know much about only being able to use the linear scheme. It may be the only scheme available for the flux. Maybe someone else knows more.

[ykanani]

Hi,

When you get "attempt to read beyond EOF" it means that OpenFOAM is trying to read something else beyond the point you ended the line with ";"
So perhaps you are not providing all information needed to use that scheme. It could be a word or number. You can check the header file (.H) corresponding to the scheme that you are trying to use. It will give you useful information regarding the syntax etc. Unfortunately, I have not used minmod myself so I cannot help you with that specific scheme.


Regards,

[Calmly]

Originally Posted by ykanani

Hi,

When you get "attempt to read beyond EOF" it means that OpenFOAM is trying to read something else beyond the point you ended the line with ";"
So perhaps you are not providing all information needed to use that scheme. It could be a word or number. You can check the header file (.H) corresponding to the scheme that you are trying to use. It will give you useful information regarding the syntax etc. Unfortunately, I have not used minmod myself so I cannot help you with that specific scheme.


Regards,

You are right. By adding "neg" or "pos" after the "Minmod" it works but I don't understand what these mean (neg and pos)

[ykanani]

Originally Posted by Calmly

You are right. By adding "neg" or "pos" after the "Minmod" it works but I don't understand what these mean (neg and pos)

It seems they have explained it in the following paper:

https://www.researchgate.net/publica...rk_of_OpenFOAM


Regards,