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

simpleFoam pipe flow case doesn't converge

Register Blogs Community New Posts Updated Threads Search

Like Tree1Likes
  • 1 Post By alexeym

 
 
LinkBack Thread Tools Search this Thread Display Modes
Prev Previous Post   Next Post Next
Old   February 7, 2014, 23:31
Default simpleFoam pipe flow case doesn't converge
  #1
Senior Member
 
Join Date: Nov 2010
Location: USA
Posts: 1,232
Rep Power: 24
me3840 is on a distinguished road
Hello all, I have a simple case of pipe flow to get me introduced to openFOAM. I've only used commercial codes in my life, so this is an interesting experience!

Alright, so I have a pipe r=0.1m and L=0.5m. I built the grid in Salome, it's ~150k cells, tet+prisms. An image of it is attached. It isn't the best, but there aren't any major issues and checkMesh agrees (I think).
http://www.dropbox.com/s/5ex5w0g9saqjqn3/checkMesh.log
http://www.dropbox.com/s/fe1ujn42eathiqm/mesh.png


Here are the residuals from the current solution I have computed. Not really much going on besides tke bottoming out:
http://www.dropbox.com/s/3o10ds4gbo6bvlf/residuals.png

I'm hoping to get just a 5m/s flow of air through this guy. My current fvSchemes and fvSolution are attached. I also attached my U and p files from /0, since I have some concerns with them.
http://www.dropbox.com/s/rxvmt6vvy0317al/p
http://www.dropbox.com/s/fg7kpwfzkv73fq5/U
https://www.dropbox.com/s/fot4xa5rwu8fha7/fvSolution
https://www.dropbox.com/s/wnv8nyuqnjtv2my/fvSchemes

Okay, here are my questions:
1. The syntax for fvSchemes was a little confusing to me, but I'm pretty sure I have selected all 2nd order schemes. Some of the syntax for the upwind ones were giving me errors. Could someone, for example, give me an example for the syntax of div(phi,U) for a 2nd order bounded upwind scheme?
2. In commercial codes I'm pretty used to starting off with 2nd order immediately. Is that too harsh a condition for openFOAM which lacks all the stability algorithms and accelerators and whatnot commercial vendors use? Even for flow in a pipe?
3. Can someone tell me what the point is behind setting a patch type in /constant/polyMesh/boundary when I'm just going to define the BC in the /0 files?
4. I'm not really convinced I'm applying the BCs I want to, I'm hoping I haven't over-specified the system. Let's start with U:
4A. I went with fixedValue everywhere. For commercial codes U is not normally specified explicitly on the outlet boundary, so I'm questioning whether or not I should set it to be the same as my inlet. I'm hoping that setting a fixed (0 0 0) at the wall boundary should be sufficient for a no-slip condition.
4B. For p, I'm pretty sure I set the BC for the inlet as zeroGradient, which is what I'd like for a constant-velocity inlet (I think..). Then I specified a constant-pressure outlet condition at the gauge value. Do these (and the conditions in U) make the problem well-posed?
5. Does anything else strike someone as strange that might be making this case not converge?

Thanks for any and all comments, I'm looking forward to doing more OpenFOAM, it's quite fun!

Last edited by me3840; February 8, 2014 at 11:03.
me3840 is offline   Reply With Quote

 


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
solving a case on flow through duct hariya03 OpenFOAM Running, Solving & CFD 0 August 10, 2013 16:30
Cannot validate simple pipe flow using simpleFoam inf.vish OpenFOAM Running, Solving & CFD 4 August 6, 2013 07:32
Pipe flow with pressure-inlet lummz FLUENT 3 October 13, 2012 13:29
About Turbulence Intensity (Pipe flow assimilated) gRomK13 Main CFD Forum 1 July 10, 2009 03:11
Terrible Mistake In Fluid Dynamics History Abhi Main CFD Forum 12 July 8, 2002 09:11


All times are GMT -4. The time now is 15:18.