|
[Sponsors] |
April 15, 2013, 08:24 |
problem from k epsilon to k omega
|
#1 |
Member
Nick Gutschow
Join Date: Jan 2013
Posts: 36
Rep Power: 13 |
Foamers,
I have been running a modified convective heat transfer solver (buoyanrBoussinesqPimpleFoam with the addition of CO2 mass fraction) to model a person breathing out in a closed room. Basically the person-like pillar breathes out hot CO2 rich air at a constant rate and his "ears" are constant pressure BCs that allow a mass balance in the room. I started with a k epsilon based approach (using k and epsilon values at the inlet I borrowed from another model without much thought, just to kick the tires a bit, which where 0.015 in both cases) and I was getting good results consistent with previous work we have done on this. Note the two k epsilon attached photos. "Profile" shows the person profile and you can see the breathing out effect and "backview" is a view from the back so you can see the "ears" doing their mass balance thing. This was all good, but over time the k epsilon scheme showed instability at the "ears". My k and epsilon values would sometimes suddenly exploded to huge values, crashing the model. So, I switch to k omega to try and fix this, and so far I have had no instability issues. BUT, the behavior of my plume is much different. Here I actually calculated a k and omega based on the breath out velocity of 0.06 m/s, turbulent intensity of 10%, and the length scale of the "mouth" which is 0.04 m (so I have k 5.55E-5 and omega 0.186). You can see the results of this in my k omega pictures attached below. Note how much more diffuse the plume is. Also, note how my back view shows a disturbing lack of symmetry on the effects of my ears. So now my questions: 1) It seems to me, if all I changed was from k epsilon to k omega, through the selection of better values of k and omega, I should be able to produce the exact same results. Is this not correct? 2) Why would I be getting this non-symmetric behavior if I have a fully symmetric domain and BCs and had fully symmetric result in a k epsilon scheme? 3) If anyone has any tips on how to perhaps just make k epsilon more stable that approach might be better as k omega is not giving me good results. 4) Any other insight into this that anyone has would be most appreciated. -NG |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Bounding epsilon or bounding omega | Stylianos | OpenFOAM | 8 | February 23, 2018 13:41 |
___Minus epsilon problem___ | maysmech | OpenFOAM | 0 | January 3, 2011 03:10 |
Can I solve this problem by Fluent? | Kai_kc | FLUENT | 1 | October 27, 2010 05:29 |
Problem with epsilon | joshK | OpenFOAM | 1 | October 8, 2010 01:55 |
Adiabatic and Rotating wall (Convection problem) | ParodDav | CFX | 5 | April 29, 2007 19:13 |