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

SOS: the nonpositive k and epsion ,how to get the correct boundary settings

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 6, 2012, 00:36
Default SOS: the nonpositive k and epsion ,how to get the correct boundary settings
  #1
New Member
 
zdong
Join Date: Aug 2011
Posts: 15
Rep Power: 4
abbott.hn is on a distinguished road
hello everybody,when i try to simulation the pantograph wind tunnel experiment,i got an divergent k and epsilon using simpleFoam or pimpleFoam with k-e turbulent modle.

i have looked though the forum a long time for a solution to this kind of divergent probelm, after changing some related settings follow the tips i got,
like checking the model length dimension ,setting the proper magnitude of k and e whenever in the internal field or boundary,setting the relax factor or tolerence
for u and p,changing the numerical scheme for k and e , even modifying the solver for k and e,

however ,all of this didn't change anything

the keypoint should be in the boundary settings,because when i turned off the turbulent modeling, the procedure can run smoothly.

those in the attachment is my cas directory , the grid is not in the constant/polyMesh but the boundary dictionary for the file size consideration.



anyone could help me out?

best regards for you,thanks .

abbott
Attached Files
File Type: gz pantagraphTunnelEx.tar.gz (3.8 KB, 4 views)
abbott.hn is offline   Reply With Quote

Old   March 6, 2012, 02:00
Default
  #2
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 782
Rep Power: 11
Bernhard is on a distinguished road
I cannot run your case without the mesh, so it is difficult to spot the error. Can you give the log file? Does it diverge in the first time-step, or only later on? Is your mesh ok (check with checkMesh)?
Bernhard is offline   Reply With Quote

Old   March 7, 2012, 21:53
Default
  #3
New Member
 
zdong
Join Date: Aug 2011
Posts: 15
Rep Power: 4
abbott.hn is on a distinguished road
Quote:
Originally Posted by Bernhard View Post
I cannot run your case without the mesh, so it is difficult to spot the error. Can you give the log file? Does it diverge in the first time-step, or only later on? Is your mesh ok (check with checkMesh)?
hello Bernhard thanks for your reply ,
the log file is in the attachments ,maybe you can find some hints from that .appreciate for your help
Attached Files
File Type: txt checkMeshLog.txt (4.3 KB, 5 views)
File Type: txt PimpleFoamlog.txt (7.1 KB, 2 views)
File Type: txt simpleFoamlog.txt (8.7 KB, 4 views)
abbott.hn is offline   Reply With Quote

Old   March 10, 2012, 09:11
Default
  #4
New Member
 
zdong
Join Date: Aug 2011
Posts: 15
Rep Power: 4
abbott.hn is on a distinguished road
come on ....
i need you help .
Even when i chang to the rhoSimplecFoam solver,the problem arise the same,and this time ,the "floating points excepting " error came out at the first iterate when solving the pressure eqution,while the simpleFoam or pimpleFoam prompt the bounding turbulent energy and e at the first iterate ,and then the "floating points excepting " error after several iterate.

i will to be crazy ,foamers ,please help me out, thanks so much
abbott.hn is offline   Reply With Quote

Old   March 12, 2012, 03:28
Default
  #5
Senior Member
 
romant's Avatar
 
Roman Thiele
Join Date: Aug 2009
Location: Stockholm, Sweden
Posts: 342
Rep Power: 9
romant is on a distinguished road
Quote:
Originally Posted by abbott.hn View Post
come on ....
i need you help .
Even when i chang to the rhoSimplecFoam solver,the problem arise the same,and this time ,the "floating points excepting " error came out at the first iterate when solving the pressure eqution,while the simpleFoam or pimpleFoam prompt the bounding turbulent energy and e at the first iterate ,and then the "floating points excepting " error after several iterate.

i will to be crazy ,foamers ,please help me out, thanks so much

Hej,

first of all calm down, you cannot expect people to answer within minutes. this is a forum of volunteers not your microsoft helpdesk hotline.

Second, please check your epsilon and k values, inlet and initial. they seem very high. I don't know what you are simulating, but your values are very large. check the tutorial for values that are more appropriate. k = 0.01 and epsilon = 0.5 and then start varying those values towards what ever you expect to have for your k and epsilon values.

try to run the simulation with "standard" values first, as for example the values that are used in the tutorials
__________________
~roman
romant is offline   Reply With Quote

Old   March 12, 2012, 23:49
Default thank for you advice
  #6
New Member
 
zdong
Join Date: Aug 2011
Posts: 15
Rep Power: 4
abbott.hn is on a distinguished road
Quote:
Originally Posted by romant View Post
Hej,

first of all calm down, you cannot expect people to answer within minutes. this is a forum of volunteers not your microsoft helpdesk hotline.

Second, please check your epsilon and k values, inlet and initial. they seem very high. I don't know what you are simulating, but your values are very large. check the tutorial for values that are more appropriate. k = 0.01 and epsilon = 0.5 and then start varying those values towards what ever you expect to have for your k and epsilon values.

try to run the simulation with "standard" values first, as for example the values that are used in the tutorials
hi
thanks for your advice ,i will follow your suggest.
i'a really a lit impatient, however, i will calm myself down and study the tutorials mor carefully.thanks .
abbott.hn 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



All times are GMT -4. The time now is 14:59.