|
[Sponsors] |
February 7, 2014, 23:31 |
simpleFoam pipe flow case doesn't converge
|
#1 |
Senior Member
Join Date: Nov 2010
Location: USA
Posts: 1,232
Rep Power: 24 |
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. |
|
|
|
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 |