interFoam error
 

Hello foamers ;

I'm trying to solve free surface problem using interFoam solver

after setting fields and specify the boundary conditions alpha1 , p , and U

i executed the interFoam solver and got that error ;

Create time

Create mesh for time = 0


Reading g
Reading field p

Reading field alpha1

Reading field U

Reading/calculating face flux field phi

Reading transportProperties

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian


--> FOAM FATAL IO ERROR:
Unable to set reference cell for field p
Please supply either pRefCell or pRefPoint


file: /home/openfoam1/Desktop/damBreakTest/system/fvSolution::PISO from line 55 to line 60.

From function void Foam::setRefCell
(
const volScalarField&,
const dictionary&,
label& scalar&,
bool
)
in file cfdTools/general/findRefCell/findRefCell.C at line 112.

FOAM exiting


thank you very much

From the error u can say that you have mentioned pRefValue and pRefPoint in the PISO loop in the system/fvSolution dictionary. adding both in the PISO loop u can solve the problem.

hello thank you for your reply


i did , i put pRefPoint 0 ; pRefValue 0 in the PISO loop

and got again an error

Create time

Create mesh for time = 0


Reading g
Reading field p

Reading field alpha1

Reading field U

Reading/calculating face flux field phi

Reading transportProperties

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar
time step continuity errors : sum local = 3.00873e-19, global = -1.87935e-19, cumulative = -1.87935e-19


--> FOAM FATAL ERROR:
incompatible dimensions for operation
[pcorr[-1 3 -1 0 0 0 0] ] == [div(phi)[0 0 -1 0 0 0 0] ]

From function checkMethod(const fvMatrix<Type>&, const GeometricField<Type, fvPatchField, volMesh>&)
in file /home/openfoam1/OpenFOAM/OpenFOAM-1.6.x/src/finiteVolume/lnInclude/fvMatrix.C at line 1219.

FOAM aborting

#0 Foam::error::printStack(Foam::Ostream&) in "/home/openfoam1/OpenFOAM/OpenFOAM-1.6.x/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/openfoam1/OpenFOAM/OpenFOAM-1.6.x/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 void Foam::checkMethod<double>(Foam::fvMatrix< double> const&, Foam::DimensionedField<double, Foam::volMesh> const&, char const*) in "/home/openfoam1/OpenFOAM/OpenFOAM-1.6.x/applications/bin/linuxGccDPOpt/interFoam"
#3 Foam::tmp<Foam::fvMatrix<double> > Foam::operator==<double>(Foam::tmp<Foam:: fvMatrix<double> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/home/openfoam1/OpenFOAM/OpenFOAM-1.6.x/applications/bin/linuxGccDPOpt/interFoam"
#4 main in "/home/openfoam1/OpenFOAM/OpenFOAM-1.6.x/applications/bin/linuxGccDPOpt/interFoam"
#5 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#6 _start at /build/buildd/glibc-2.9/csu/../sysdeps/i386/elf/start.S:122
Aborted

sorry, it is pRefCell not pRefPoint. You can check the variable declarations in createFields.C in interFoam source directory.

Regarding pcorr error, have a look at the "continuityErrs.H" file in source directory. What I feel from the error, you might messed up with dimension, try to see one-to-one matching with tutorial case, if you are doing for the first time.

All the best.

Regards
Santosh...

yes i do that for the first time ,, I'm trying to simulate the behavior of a water droplet on a hydrophobic flat plate under wind drag and the plate have to rotate

how can i make the flat plat rotates, is it possible in openfoam ?

thank you

Yeah u can model plate rotation using derived BC, i.e, rotatingWallVelocity.

Santosh.

oh , if it can do ,, does it put into the consideration the centrifugal force (caused by the rotation ) , if it do ,that sound very nice

Hello,

I am kind of new user in openfoam. I want to simulate such problem that you mentioned in forum( simulating a droplet on hydrophobic surface). would you please help to simulate this.or can u send me your setup..I would be so grateful if you help me through this .

thanks

Mehran

Regarding the error:

The solver needs to know what and where the reference pressure is. Add the following to your fvsolution file:

PISO
{
...
pRefPoint (-0.081 -0.0257 8.01);
pRefValue 1e5;
}