[noodleluvmay]

Hello,

I am very new to this forum and CFD simulation on OpenFoam in general, this question might be slightly dumb but hopefully someone could shine some light.

OpenFoam Version: 18.10 installed from cfdsupport.

I am trying the OpenFoam motorBike tutorial, have finished meshing using snappyHexMesh, trying to solve first using potentialFoam then simpleFoam as I was told that this practice would increase stability in solving the case. However, while I was executing the potentialFoam solver through command:

$ potentialFoam.exe -noFunctionObjects -writep

the potential flow solved fine, then right after the line
"Calculating approximate pressure field", OpenFoam outputs error then exits:

-----

--> FOAM FATAL IO ERROR:
keyword div(div(phi,U)) is undefined in dictionary "D:/Program_Files/OpenFOAM/18.10/liuni-dev/run/tutorials/motorBike/system/fvSchemes.divSchemes"

file: D:/Program_Files/OpenFOAM/18.10/liuni-dev/run/tutorials/motorBike/system/fvSchemes.divSchemes from line 30 to line 35.

From function const Foam::entry& Foam::dictionary::lookupEntry(const Foam::word&, bool, bool) const
in file db/dictionary/dictionary.C at line 566.

FOAM exiting

//Screenshot:



//////////////////////////////////////////////////////////////

Checking in the fvSchemes dictionary, I've not changed anything from the default tutorial file, div(phi,U) seems to be defined okay?

FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
default steadyState;
}

gradSchemes
{
default Gauss linear;
grad(U) cellLimited Gauss linear 1;
}

divSchemes
{
default none;
div(phi,U) bounded Gauss linearUpwindV grad(U);
div(phi,k) bounded Gauss upwind;
div(phi,omega) bounded Gauss upwind;
div((nuEff*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
default Gauss linear corrected;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}

wallDist
{
method meshWave;
}


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

//Screenshot:


////////////////////////////////////////////////////

I'm not sure if I have described the question detailed enough, please let me know if any other information needed in order to pinpoint the problem!

Thanks in advance!

[clapointe]

Well the error appears to be telling you that a line with "div(div(phi,U))" in your div schemes does not exist (2 divs). So I would try adding one (P.s. -- the tutorials for potentialFoam include a line like "div(div(phi,U)) Gauss linear" in fvSchemes).

Caelan

[noodleluvmay]

Hi Caelan,

Thank you very much for the reply. For some reason my new version of OpenFOAM did not include the div(div(phi, U)) in the tutorial. I have added a line of

div(div(phi,U)) Gauss linear;

It works fine now.

Thanks!

Nian