CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   p relaxationFactor in twoPhaseEulerFoam (http://www.cfd-online.com/Forums/openfoam/91927-p-relaxationfactor-twophaseeulerfoam.html)

hkhosravi August 26, 2011 11:42

p relaxationFactor in twoPhaseEulerFoam
 
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>::prevIter() const
in file /home/hamed/OpenFOAM/OpenFOAM-2.0.0/src/OpenFOAM/lnInclude/GeometricField.C at line 844.

FOAM aborting

#0 Foam::error::printStack(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??? :confused:

wyldckat August 26, 2011 12:51

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 August 26, 2011 14:16

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 August 26, 2011 14:39

Quote:

Originally Posted by hkhosravi (Post 321822)
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 ;)

Quote:

Originally Posted by hkhosravi (Post 321822)
i`m sure the problem relate to "p" relaxation factor, because when i delete it, everything is correct !!

:eek: 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 August 26, 2011 15:53

:eek: 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 August 26, 2011 16:11

After googling it a bit... :eek: 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 August 27, 2011 02:18

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

thanks Bruno

Regards


All times are GMT -4. The time now is 09:02.