[ray.hill4@ntlworld.com]

This a bug that I have researched and has been reqularly reported over the years. I have openFoam 4.1 on Xubuntu 17.04. After about 12 attempts ate apply various older recommendations I have made no progress, but think it may be due to a missing function in a complied library somewhere. I believe most of the software to be working OK since I successfully completed the airfoil2D example earlier today. I only have about 10 hours experience of using this.

Thanks in advance for advice that may follow.

I have tried numerous variations of the div((nrEff.. statement) my latest fvSchemes file follows:

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

ddtSchemes
{
default steadyState;
}

gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
grad(U) Gauss linear;
}

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


laplacianSchemes
{
default Gauss linear orthogonal ;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
p;
}
wallDist
{
method meshWave;
}

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

[vatavuk]

Hi Ray,

Apparently your file was prepared for OpenFoam 2.1.0 and you are trying to run it on OF 4.1.
I suggest that you find a tutorial prepared for OF 4.1, that is similar to your case and use the fvSchemes file from the tutorial in your case.
Best regards,
Paulo

[ray.hill4@ntlworld.com]

Yes I would look there if I could find one but none seem to exist! I did originally start from the tutorial stuff bundled with the 4.1 build. The official documentation doesn't seem to have much to offer on the subject, as yet. The fact there is some backward history makes me think , as you say documentation (some faults were fixed this way) or a repetitive omission from a rolling forward set of original build statements.

Update: have run foamDictionary system/fvSolution and foamDictionary system/fvSchemes and cut and pasted the output of each into the fvShemes and fvSolution files in case\system it now runs a bit further but now fails looking for rho, I will submit this as a separate issue. I found that it was important that any variables in divSchemes sections needed to be declared in the solvers section of fvSolution.

[vatavuk]

Hi Ray,

Which solver are you using?
Try this: Delete the line with the keyword and run the solver. Do you receive an error message complaining the keyword is not present?

Best regards,
Paulo

[ray.hill4@ntlworld.com]

Tried your suggestion knocked out all rho references and tried with both simpleFoam and icoFoam solvers. Glad no one is paying me to do this, I would be fired by now! I was actually trying to get a rocket shape (from an .STL) and then try treating it as airfoil because I was curious how much (if any) aerodynamic lift was generated by the rocket at flight angles >10deg (initially using the incompressible velocity range). At least the nuEff issue got fixed (maybe that helps others), but I have spent far too much time on this so have to put it aside and move onto more productive things. Thanks again for your input.