[kilroy]

I have found OpenFOAM model files of a wigley hull example with an experimental validation report. The model was programmed to be running in version 1.7.x. To be able to make it run in the 2.2.0 version, I had to make the following changes.


In the "fvSchemes" file, "divSchemes" section:

Old code:
divSchemes
{
div(rho*phi,U) Gauss linearUpwind cellLimited Gauss linear 1.0;
div(phi,alpha) Gauss vanLeer;
div(phirb,alpha) Gauss interfaceCompression;
div(phi,k) Gauss upwind;
div(phi,omega) Gauss upwind;
div(phi,R) Gauss upwind;
div(R) Gauss linear;
div(phi,nuTilda) Gauss upwind;
div((nuEff*dev(grad(U).T()))) Gauss linear;
}

New code:
divSchemes
{
div(rho*phi,U) Gauss linearUpwindV grad(U);
div(phi,alpha) Gauss vanLeer;
div(phirb,alpha) Gauss interfaceCompression;
div(phi,k) Gauss upwind;
div(phi,omega) Gauss upwind;
div(phi,R) Gauss upwind;
div(R) Gauss linear;
div(phi,nuTilda) Gauss upwind;
div((muEff*dev(T(grad(U))))) Gauss linear;
div((nuEff*dev(T(grad(U))))) Gauss linear;
}


In the "fvSolution" file, the following section is added:

"(U|k|omega)Final"
{
$U;
tolerance 1e-06;
relTol 0;
}

In the "nut" file, "boundary" section:

Old code
hull
{
type nutWallFunction;
value uniform 0;
}

New code:
hull
{
type nutkWallFunction;
value uniform 0;
}


Now the file is running perfectly on version 2.2.0. But the drag forces I get from it does not match with the values presented in the report. The values I read is 20% over.

Do you think the changes I made can effect the the drag output of the model?

Thanks,

[gaduo]

Hi, kilroy, I hope to get an OpenFOAM model files of a wigley hull example. Do you know where to get it? Thanks