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

Issue with SimpleFoam with Sharp corners

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 3, 2008, 17:48
Default I have been trying to run a no
  #1
Member
 
kshitij neroorkar
Join Date: Mar 2009
Location: Michigan, USA
Posts: 32
Rep Power: 8
kdneroorkar is on a distinguished road
I have been trying to run a nozzle test case using simpleFoam with FixedValue pressure at inlet and outlet boundaries and zeroGradient velocity at both boundaries. I am getting a strange acceleration zone at the exit of the nozzle due to a very high pressure cell.

I know that a better boundary condition is an outlet pressure and an inlet velocity condition but I donot have this information for this case. Could someone please tell me if I am applying some inconsistent boundary condition or does simpleFoam have a problem with sharp corners.
I am attaching my case folder below for anyne who wishes to run it

kdneroorkar is offline   Reply With Quote

Old   September 4, 2008, 09:37
Default sorry here is the image http:/
  #2
Member
 
kshitij neroorkar
Join Date: Mar 2009
Location: Michigan, USA
Posts: 32
Rep Power: 8
kdneroorkar is on a distinguished road
sorry here is the image
the test case is too big to post on the board so if someone wants to see the case I can email it to them.
kdneroorkar is offline   Reply With Quote

Old   September 10, 2008, 15:24
Default Kshitij, Interesting proble
  #3
Member
 
David P. Schmidt
Join Date: Mar 2009
Posts: 70
Rep Power: 8
schmidt_d is on a distinguished road
Kshitij,

Interesting problem. It might help to list your email address:

Kshitij Neroorkar <kneroork@engin.umass.edu>

-David
schmidt_d is offline   Reply With Quote

Old   September 10, 2008, 16:56
Default Kshitij Neroorkar: Is this a l
  #4
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 208
Rep Power: 9
louisgag is on a distinguished road
Send a message via ICQ to louisgag
Kshitij Neroorkar: Is this a laminar simulation? Otherwise, what turbulence model are you using?


regards
louisgag is offline   Reply With Quote

Old   September 10, 2008, 17:06
Default it is laminar
  #5
Member
 
kshitij neroorkar
Join Date: Mar 2009
Location: Michigan, USA
Posts: 32
Rep Power: 8
kdneroorkar is on a distinguished road
it is laminar
kdneroorkar is offline   Reply With Quote

Old   September 12, 2008, 19:01
Default maybe if you try the same prob
  #6
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 208
Rep Power: 9
louisgag is on a distinguished road
Send a message via ICQ to louisgag
maybe if you try the same problem with a coarser or finer mesh your results will improve
louisgag is offline   Reply With Quote

Old   September 14, 2008, 13:54
Default I would first run checkMesh an
  #7
Senior Member
 
Srinath Madhavan (a.k.a pUl|)
Join Date: Mar 2009
Location: Edmonton, AB, Canada
Posts: 698
Rep Power: 12
msrinath80 is on a distinguished road
I would first run checkMesh and ensure that your mesh passes all checks.
msrinath80 is offline   Reply With Quote

Old   September 15, 2008, 10:05
Default i ran checkMesh and the mesh d
  #8
Member
 
kshitij neroorkar
Join Date: Mar 2009
Location: Michigan, USA
Posts: 32
Rep Power: 8
kdneroorkar is on a distinguished road
i ran checkMesh and the mesh did pass all tests.
Also, the mesh is pretty fine with 140,000 cells
kdneroorkar is offline   Reply With Quote

Old   September 15, 2008, 13:56
Default Is this some kind of axisymmet
  #9
Senior Member
 
Srinath Madhavan (a.k.a pUl|)
Join Date: Mar 2009
Location: Edmonton, AB, Canada
Posts: 698
Rep Power: 12
msrinath80 is on a distinguished road
Is this some kind of axisymmetric case? Can you go to www.rapidshare.com and upload your file there and paste the download link (that it gives you) here. The case might help give a better idea on what could be wrong.
msrinath80 is offline   Reply With Quote

Old   September 15, 2008, 14:09
Default here is the link http://rap
  #10
Member
 
kshitij neroorkar
Join Date: Mar 2009
Location: Michigan, USA
Posts: 32
Rep Power: 8
kdneroorkar is on a distinguished road
here is the link

http://rapidshare.com/files/14553772...le.tar.gz.html
kdneroorkar is offline   Reply With Quote

Old   September 15, 2008, 16:44
Default I just checked your mesh and c
  #11
Senior Member
 
Srinath Madhavan (a.k.a pUl|)
Join Date: Mar 2009
Location: Edmonton, AB, Canada
Posts: 698
Rep Power: 12
msrinath80 is on a distinguished road
I just checked your mesh and case setup. I could not identify any obvious errors. However, your log snippet worries me a lot.

Listed below are the snippets for time step continuity errors at the 1st, 10th, 20th and 100th, 500th, 2000th and 4000th iteration. Your cumulative error does not look good, neither does the 'sum local'.

time step continuity errors : sum local = 1.02549, global = -0.0188398, cumulative = -0.0188398

time step continuity errors : sum local = 0.0699843, global = -0.000584316, cumulative = 0.0335601

time step continuity errors : sum local = 0.00333286, global = 0.000201054, cumulative = 0.0354029

time step continuity errors : sum local = 0.000304726, global = -3.63236e-06, cumulative = 0.0340224

time step continuity errors : sum local = 0.000303121, global = 3.93966e-06, cumulative = 0.0318931

time step continuity errors : sum local = 0.000308258, global = -1.41759e-05, cumulative = 0.0297984

time step continuity errors : sum local = 0.00030529, global = 1.66529e-05, cumulative = 0.0286185

I recommend going over the first few posts of this thread[1] to clarify what each error term means.

Things you may want to try:

1. Increase the number of non-orthogonal correctors (at least 1 corrector to start with).

2. Like Hrv recommends in that other[1] post, converge the pressure more tightly (at least in the beginning). Set the tolerance in fvSolution to (say) 1e-08 instead of 1e-06.

I am no expert in CFD, so you may want to recheck your Boundary conditions as well.

References:

[1] http://www.cfd-online.com/OpenFOAM_D...es/1/1671.html
msrinath80 is offline   Reply With Quote

Old   September 15, 2008, 17:06
Default Oh, and one more thing. You ca
  #12
Senior Member
 
Srinath Madhavan (a.k.a pUl|)
Join Date: Mar 2009
Location: Edmonton, AB, Canada
Posts: 698
Rep Power: 12
msrinath80 is on a distinguished road
Oh, and one more thing. You can always use icoFoam to run your laminar simpleFoam problem in a transient manner. The transient simulation should converge to steady state if the problem is indeed steady. I recommend maintaining a Max. Courant number of at most 0.75 even if you are not interested in resolving any transient behavior.
msrinath80 is offline   Reply With Quote

Old   September 15, 2008, 17:57
Default Hi Srinath I have tried to ru
  #13
Member
 
kshitij neroorkar
Join Date: Mar 2009
Location: Michigan, USA
Posts: 32
Rep Power: 8
kdneroorkar is on a distinguished road
Hi Srinath
I have tried to run icoFoam using the same boundary conditions and i got a very similar result. In this case the continuity errors look alright. I am just pasting the last few time step continuity errors below

time step continuity errors : sum local = 3.42216e-13, global = -1.07768e-14, cumulative = 8.9745e-10

time step continuity errors : sum local = 4.29076e-12, global = -8.72646e-15, cumulative = 8.97432e-10

time step continuity errors : sum local = 3.58808e-13, global = -8.47658e-15, cumulative = 8.97417e-10

this was using 2 non-orthogonal correctors.

Thanks a lot for your help
Kshitij
kdneroorkar is offline   Reply With Quote

Old   September 15, 2008, 18:12
Default How many time steps did you ru
  #14
Senior Member
 
Srinath Madhavan (a.k.a pUl|)
Join Date: Mar 2009
Location: Edmonton, AB, Canada
Posts: 698
Rep Power: 12
msrinath80 is on a distinguished road
How many time steps did you run and what was your Max. Co? Try plotting the time variation of velocity or pressure at a particular point in the domain to see if the variable has indeed converged. See the oodles/pitzDaily tutorial to see how to add probes to your simulation.
msrinath80 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
Courant no%23 and sharp corner in 3D jam OpenFOAM Running, Solving & CFD 6 November 18, 2008 09:21
Avoiding Skew in Tight Corners Marc FLUENT 4 July 23, 2007 15:21
Flow near sharp corners Harish Main CFD Forum 4 February 21, 2007 22:55
How to treat the 2D corners when using NSCBC? leaf Main CFD Forum 0 May 26, 2006 04:12
why the zones gets split at sharp corners/bends Laxminarayana FLUENT 0 February 16, 2005 03:13


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