[hkhosravi]

dear foamers

i`m using twoPhaseEulerFoam for my simuation.
relaxationfactor for my simulation is:
relaxationFactors
{
Ua 0.7;
Ub 0.7;
p 0.3;
alpha 0.2;
beta 0.2;
Theta 0.2;
k 0.4;
epsilon 0.4;
}

but an error was occure only for p (pressure relaxation factor) !!
the error is:
Courant Number mean: 0.00035655 max: 0.0004
Max Ur Courant Number = 0.0004
deltaT = 1.1999e-05
Time = 1.1999e-05

PIMPLE: iteration 1
DILUPBiCG: Solving for alpha, Initial residual = 1.08695e-06, Final residual = 1.52185e-22, No Iterations 1
Dispersed phase volume fraction = 0.11 Min(alpha) = 0 Max(alpha) = 0.55
DILUPBiCG: Solving for alpha, Initial residual = 8.6948e-07, Final residual = 1.52191e-22, No Iterations 1
Dispersed phase volume fraction = 0.11 Min(alpha) = 0 Max(alpha) = 0.55
kinTheory: max(Theta) = 1e-05
kinTheory: min(nua) = 2.94999e-08, max(nua) = 2.97152e-06
kinTheory: min(pa) = 0, max(pa) = 1.20197e-10
GAMG: Solving for p, Initial residual = 1, Final residual = 0.0810062, No Iterations 4


--> FOAM FATAL ERROR:
previous iteration field
IOobject: volScalarField p "/home/hamed/OpenFOAM/hamed-2.0.0/mycase/al96-H0=10-03/0"

not stored. Use field.storePrevIter() at start of iteration.

From function GeometricField<Type, PatchField, GeoMesh>:revIter() const
in file /home/hamed/OpenFOAM/OpenFOAM-2.0.0/src/OpenFOAM/lnInclude/GeometricField.C at line 844.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) in "/home/hamed/OpenFOAM/OpenFOAM-2.0.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/hamed/OpenFOAM/OpenFOAM-2.0.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2
in "/home/hamed/OpenFOAM/OpenFOAM-2.0.0/platforms/linuxGccDPOpt/bin/twoPhaseEulerFoam"
#3
in "/home/hamed/OpenFOAM/OpenFOAM-2.0.0/platforms/linuxGccDPOpt/bin/twoPhaseEulerFoam"
#4
in "/home/hamed/OpenFOAM/OpenFOAM-2.0.0/platforms/linuxGccDPOpt/bin/twoPhaseEulerFoam"
#5
in "/home/hamed/OpenFOAM/OpenFOAM-2.0.0/platforms/linuxGccDPOpt/bin/twoPhaseEulerFoam"
#6 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#7
in "/home/hamed/OpenFOAM/OpenFOAM-2.0.0/platforms/linuxGccDPOpt/bin/twoPhaseEulerFoam"
Aborted

any solution or idea???

[wyldckat]

Greetings hkhosravi,

OpenFOAM's take on folder and file names is that the file system reflects the program variables... with the exception of time folders which should be properly formatted numbers.
Therefore "al96-H0=10-03" is a very bad folder name! I'm stuned how OpenFOAM didn't stop right at the beginning telling you that the name "al96-H0=10-03" is invalid...

Try again with another folder name for your case! Something more... simple! You could try with "al96_H0__10_03".

Best regards,
Bruno

[hkhosravi]

Hi Bruno

thanks for quick reply.
I changed the folder name to "al96_03" also "al96", but there are the same error.

i`m sure the problem relate to "p" relaxation factor, because when i delete it, everything is correct !!

[wyldckat]

Originally Posted by hkhosravi

I changed the folder name to "al96_03" also "al96", but there are the same error.

Well, at least that's a relief... it would be very annoying that one's liberty to create crazy folder names (programmatically-wise) would be restricted as well

Originally Posted by hkhosravi

i`m sure the problem relate to "p" relaxation factor, because when i delete it, everything is correct !!

How could I have not spotted this before... Does the file "0/p" exist and have the necessary boundaries and field? That's what the error message is talking about!
In case you can't define it yourself, I think you can use the following instructions to write the p field: http://openfoamwiki.net/index.php/Ti...gisteredObject

[hkhosravi]

How could I have not spotted this before... Does the file "0/p" exist and have the necessary boundaries and field? That's what the error message is talking about![/QUOTE]

The file "0/p" exist and have the correct BC, because I can run the case without using relaxation factor for "p".

also, In the stored time directory, "p" field exist and I can see pressure field in paraview.

[wyldckat]

After googling it a bit... It's a bug! See http://www.openfoam.com/mantisbt/view.php?id=245

It should be fixed in OpenFOAM 2.0.1, so you should upgrade!

Best regards,
Bruno

[hkhosravi]

yes, it`s a bug and solved in OF 2.0.1.

thanks Bruno

Regards