|
[Sponsors] |
January 11, 2007, 17:18 |
Hi all,
I'm having fits wit
|
#1 |
Member
Michael Rangitsch
Join Date: Mar 2009
Location: Midland, Michigan, USA
Posts: 31
Rep Power: 17 |
Hi all,
I'm having fits with a non-reacting case mixing fluids with different temperatures. I'm continually getting a temperature out of range error in the janafThermo routines. The domain is a quarter of a pipe, with zero gradient boundaries on the walls, and quarter symmetry used (2 symmetry planes). Looking at the data, the boundary value at the input has the correct value (300K). The temperatures at the wall of the inlet plane and the symmetry edges of the inlet plane go down as the solution progresses, eventually giving me the low temperature error. I've tried restricting the courant number to 0.1 (and lower), but it just delays the problem. I'm using OpenFoam 1.3 on a linux box... Any hints would be appreciated... Mike |
|
January 11, 2007, 19:02 |
Sounds like a boundry conditio
|
#2 |
Senior Member
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,907
Rep Power: 33 |
Sounds like a boundry condition problem to me... Let's investigate:
First, check the mesh: is the non-orthogonality bad. Running checkMesh should tell you all you need to know. Second, check if you are satisfying continuity. Have a look at local and global continuity error: they should be of the same order of magnitude as the pressure solver convergence tolerance. If this is too big, try tightening the pressure tolerance. If all is well, let's try first order discretisation: switch the convection discretisation to upwind and diffusion to Gauss linear limited 0.5. If the error does not go away, I'm pretty sure it's your boundary conditions. Keep us posted, Hrv
__________________
Hrvoje Jasak Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk |
|
January 12, 2007, 13:01 |
Hi Hrvoje,
The grid checked
|
#3 |
Member
Michael Rangitsch
Join Date: Mar 2009
Location: Midland, Michigan, USA
Posts: 31
Rep Power: 17 |
Hi Hrvoje,
The grid checked out just fine. OK's on all parameters (it's a mixed grid, but the part with problems was nice, uniform hex's, very orthogonal). I tightened down on the pressure convergence tolerance (10e-12 instead of 10e-9), as you suggested. It seems to be working. The walls and and symmetry boundaries at the inlet are picking up the correct (inlet) temperature, and there are no excursions below the inlet value anywhere in the domain. Thanks for the help! Mike |
|
January 12, 2007, 13:11 |
Hi Mike,
The tolerance soun
|
#4 |
Senior Member
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,907
Rep Power: 33 |
Hi Mike,
The tolerance sounds awfully tight (might be necessary though). Is this by any chance a very large hex+tet grid with lots of flow in the hex part and nothing much happening in a big tet box? You bay be better off trying the AMG solver on the pressure equation - that has got better smoothing properties and may let you get away without such tight convergence. In any case, it should be faster... Hrv
__________________
Hrvoje Jasak Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk |
|
January 12, 2007, 13:39 |
Hi again Hrvoje,
The domain
|
#5 |
Member
Michael Rangitsch
Join Date: Mar 2009
Location: Midland, Michigan, USA
Posts: 31
Rep Power: 17 |
Hi again Hrvoje,
The domain is two cocentric cylinders connected by a couple holes (the inner cylinder is capped off). The tet grids are near the holes, where the two hex grids meet. Where the flow goes through the hole, the grid is actually triangular wedges. I've imported the grid from Gambit. I'm trying to initialize the flow to simulate an ignition. I'll try switching to the AMG solver for pressure and see what happens. Thanks again for the help... Mike |
|
February 28, 2007, 16:22 |
Hello again all,
I've been
|
#6 |
Member
Michael Rangitsch
Join Date: Mar 2009
Location: Midland, Michigan, USA
Posts: 31
Rep Power: 17 |
Hello again all,
I've been doing other things, but the problem didn't go away, it just took longer to become evident. The temperature at the symmetry planes continually decreases until I get the janaf bound error. The flow field looks correct, but |
|
March 1, 2007, 06:55 |
Hello,
I had a maybe simila
|
#7 |
New Member
Ulrich Brammer
Join Date: Mar 2009
Location: Berlin, Germany
Posts: 1
Rep Power: 0 |
Hello,
I had a maybe similar problem with reactingFoam. I did not define for all reactants the conditions at startTime. I thought this is done with the definition of ydefault. The simulations later showed that the symmetryplane was treated like a wall regarding this one reactant. The problem disappeared with explicitly defining the conditions for this reactant at startTime. Ulrich |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
zero temperatures | Georg | CFX | 1 | December 1, 2008 17:19 |
wall and symmetry boundary conditions | malik | FLUENT | 2 | April 22, 2008 08:48 |
unrealistic temperatures | Sudipto Mukhopadhyay | FLUENT | 3 | April 1, 2007 12:10 |
Nastran Temperatures as Wall boundary condition | Holger | Siemens | 3 | September 10, 2005 12:59 |
symmetry and rotating wall | eric | FLUENT | 1 | October 29, 2003 10:39 |