[bojiezhang]

hello everyone:
when I run the case groovyWaveTank (http://idisk.mac.com/egpaterson-Publ...ovyWaveTank.gz), it says error
" --> FOAM FATAL IO ERROR:
keyword PIMPLE is undefined in dictionary "/home/host/OpenFOAM/OpenFOAM-2.0.0/swak4Foam/Examples/groovyBC/groovyWaveTank/system/fvSolution".
I have downloaded swak4Foam and compiled it well and I do not know how to solve it, please help me, thank you in advance!
bojiezhang

[Phicau]

Hi

just change word "SIMPLE" to "PIMPLE" in your fvSolution dictionary.

[bojiezhang]

Dear Phicau,
Thank you for your reply! There is no SIMPLE in the fvSolution dictionary.The whole content is:
PISO
{
momentumPredictor no;
nCorrectors 3;
nNonOrthogonalCorrectors 0;
nGammaCorr 1;
nGammaSubCycles 3;
cGamma 1;
pdRefCell 0;
pdRefValue 0;
pRefProbe
{
fields (pd);
probeLocations ((0.51 0.51 0.51));
};
}

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

But when I change PISO to PIMPLE , the error like this:

PIMPLE: Operating solver in PISO mode

Reading field p_rgh



--> FOAM FATAL IO ERROR:
cannot find file

file: /home/host/OpenFOAM/OpenFOAM-2.0.0/swak4Foam/Examples/groovyBC/groovyWaveTank/0/p_rgh at line 0.

From function regIOobject::readStream()
in file db/regIOobject/regIOobjectRead.C at line 73.

FOAM exiting

can you help me?
bojiezhang

[Phicau]

Sorry, I meant PIMPLE instead of PISO...

[bojiezhang]

Originally Posted by Phicau

Sorry, I meant PIMPLE instead of PISO...

But when I change PISO to PIMPLE , the error like this:

PIMPLE: Operating solver in PISO mode

Reading field p_rgh



--> FOAM FATAL IO ERROR:
cannot find file

file: /home/host/OpenFOAM/OpenFOAM-2.0.0/swak4Foam/Examples/groovyBC/groovyWaveTank/0/p_rgh at line 0.

From function regIOobject::readStream()
in file db/regIOobject/regIOobjectRead.C at line 73.

FOAM exiting

[Phicau]

Do you have p_rgh file (which provides the BCs for pressure) in the 0 directory?

[bojiezhang]

Originally Posted by Phicau

Do you have p_rgh file (which provides the BCs for pressure) in the 0 directory?

no, just have the file "pd", and the fvsolution also use the name pd as below:
pd GAMG
{
tolerance 1e-8;
relTol 0.05;

smoother GaussSeidel;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;

cacheAgglomeration false;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
};

pdFinal PCG
{
preconditioner
{
type GAMG;
tolerance 1e-8;
relTol 0;

nVcycles 2;

smoother GaussSeidel;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;

cacheAgglomeration false;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
};

tolerance 1e-8;
relTol 0;
maxIter 20;
};

[Phicau]

Then I guess you will have to rewrite certain things of your case in order for it to work, since it seems it is from a very old OF version.

Take a look at damBreak tutorial case and modify consistently the files.

[bojiezhang]

Originally Posted by Phicau

Do you have p_rgh file (which provides the BCs for pressure) in the 0 directory?

I use the swak4Foam in the openfoam 2.0.0, not in the openfoam-extended, does it matter?

[bojiezhang]

But the file version which written is 2.0
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object pd;
}

[gschaider]

Originally Posted by bojiezhang

I use the swak4Foam in the openfoam 2.0.0, not in the openfoam-extended, does it matter?

The problem you're experiencing have got to do with the face that the waveTank-case was set up for OF and 1.5 and interFoam went through some changes since then. The best thing would be to compare the dambreak-case for 1.5 and 2.0 and apply the same changes you see there to the waveTank

[bojiezhang]

Thank you very much, I know it now! But I have another question, I defined the alpha1 inlet like:
inlet
{
type groovyBC;
valueExpression "(pos().y<=A*cos(-w*time())) ? 1 : 0";
variables "l=1.56;A=0.01;g=vector(0,-9.81,0);k=2*pi/l;w=sqrt(k*mag(g));";
timelines ();
}

but it does not work, and says that
"--> FOAM Warning :
From function groovyBCFvPatchField<Type>::groovyBCFvPatchField(c onst fvPatch& p,const DimensionedField<Type, volMesh>& iF,const dictionary& dict)
in file groovyBCFvPatchField.C at line 119
No value defined for alpha1 on inlet therefore using 40{0}
"

No value is defined in the inlet boundary and the value is set to be zero at last, I do not know how to modify it, can you help me? Thank you in advance!

[gschaider]

Originally Posted by bojiezhang

Thank you very much, I know it now! But I have another question, I defined the alpha1 inlet like:
inlet
{
type groovyBC;
valueExpression "(pos().y<=A*cos(-w*time())) ? 1 : 0";
variables "l=1.56;A=0.01;g=vector(0,-9.81,0);k=2*pi/l;w=sqrt(k*mag(g));";
timelines ();
}

but it does not work, and says that
"--> FOAM Warning :
From function groovyBCFvPatchField<Type>::groovyBCFvPatchField(c onst fvPatch& p,const DimensionedField<Type, volMesh>& iF,const dictionary& dict)
in file groovyBCFvPatchField.C at line 119
No value defined for alpha1 on inlet therefore using 40{0}
"

No value is defined in the inlet boundary and the value is set to be zero at last, I do not know how to modify it, can you help me? Thank you in advance!

Just add a

value uniform 0;

line or so to the BC. This warning is only of consequence if another groovyBC uses alpha1 before that BC is evaluated for the first time. Very unlikely that this will happens, that's why it is only a warning

Bernhard

[bojiezhang]

Originally Posted by gschaider

Just add a

value uniform 0;

line or so to the BC. This warning is only of consequence if another groovyBC uses alpha1 before that BC is evaluated for the first time. Very unlikely that this will happens, that's why it is only a warning

Bernhard

Yeah, it does work, thank you very much! But if I add "value uniform 0" to the inlet boundary, when the wave begins to transmit, it is not like wave and the phenomenon is strange. If I add "value uniform 1" to it, the phenomenon in the boundary is not true! Do you have any advice? Thank you!

[gschaider]

Originally Posted by bojiezhang

Yeah, it does work, thank you very much! But if I add "value uniform 0" to the inlet boundary, when the wave begins to transmit, it is not like wave and the phenomenon is strange. If I add "value uniform 1" to it, the phenomenon in the boundary is not true! Do you have any advice? Thank you!

Not familiar enough with interFoam to say when which BC is evaluated. Maybe setting the value with funkySetBoundaryCondition to be consistent with the evaluated value at t=0 may help

[bojiezhang]

Originally Posted by gschaider

Not familiar enough with interFoam to say when which BC is evaluated. Maybe setting the value with funkySetBoundaryCondition to be consistent with the evaluated value at t=0 may help

OK! Thank you very much! I will try it later!

[idefix]

I have the same problem now. How did you solve yours?

[idefix]

Here´s the error massage:

--> FOAM FATAL IO ERROR:
cannot find file

file: /OpenFOAM/test/0.001/p_rgh at line 0.

From function regIOobject::readStream()
in file db/regIOobject/regIOobjectRead.C at line 73.

FOAM exiting

[JBrake]

I am running a model, similar to the propeller tutorial, where I have a cyclic AMI between a cylinder and the rest of the domain. I get the following error when running:

--> FOAM FATAL ERROR:
cannot find file "/home/jbrakefi/OpenFOAM/olaFlow/My_olaFlow/dynamic_OWSC/0.01/p_rgh"

From function virtual Foam::autoPtr<Foam::ISstream> Foam::fileOperations::uncollatedFileOperation::rea dStream(Foam::regIOobject&, const Foam::fileName&, const Foam::word&, bool) const
in file global/fileOperations/uncollatedFileOperation/uncollatedFileOperation.C at line 505.

FOAM exiting

I have tried searching for answers on this forum and elsewhere with no luck so far. I don't know where to begin to troubleshoot this one.