CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

p relaxationFactor in twoPhaseEulerFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   August 26, 2011, 11:42
Default p relaxationFactor in twoPhaseEulerFoam
  #1
New Member
 
H Kh
Join Date: Dec 2010
Location: Tehran , Iran
Posts: 16
Rep Power: 6
hkhosravi is on a distinguished road
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???
hkhosravi is offline   Reply With Quote

Old   August 26, 2011, 12:51
Default
  #2
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,488
Blog Entries: 34
Rep Power: 86
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
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
wyldckat is offline   Reply With Quote

Old   August 26, 2011, 14:16
Default
  #3
New Member
 
H Kh
Join Date: Dec 2010
Location: Tehran , Iran
Posts: 16
Rep Power: 6
hkhosravi is on a distinguished road
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 !!
hkhosravi is offline   Reply With Quote

Old   August 26, 2011, 14:39
Default
  #4
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,488
Blog Entries: 34
Rep Power: 86
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Quote:
Originally Posted by hkhosravi View Post
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 View Post
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
wyldckat is offline   Reply With Quote

Old   August 26, 2011, 15:53
Default
  #5
New Member
 
H Kh
Join Date: Dec 2010
Location: Tehran , Iran
Posts: 16
Rep Power: 6
hkhosravi is on a distinguished road
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.
hkhosravi is offline   Reply With Quote

Old   August 26, 2011, 16:11
Default
  #6
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,488
Blog Entries: 34
Rep Power: 86
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
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
wyldckat is offline   Reply With Quote

Old   August 27, 2011, 02:18
Default
  #7
New Member
 
H Kh
Join Date: Dec 2010
Location: Tehran , Iran
Posts: 16
Rep Power: 6
hkhosravi is on a distinguished road
yes, it`s a bug and solved in OF 2.0.1.

thanks Bruno

Regards
hkhosravi is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Something wrong in UEqns.H within twoPhaseEulerFoam cheng1988sjtu OpenFOAM 2 June 24, 2011 10:48
twoPhaseEulerFoam freemankofi OpenFOAM 0 May 23, 2011 16:24
stratified horizontal two phase flow usinfg twoPhaseEulerFoam karthik1414 OpenFOAM 0 April 12, 2011 09:57
problems in Two Phase flow using twoPhaseEulerFoam with OpenFoam 1.6 raagh77 OpenFOAM Running, Solving & CFD 0 March 6, 2010 06:11
TwoPhaseEulerFoam bed tutorial case stable in 1.5, crashes in 1.6 hemph OpenFOAM 3 December 5, 2009 05:19


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