CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   Phoenics (https://www.cfd-online.com/Forums/phoenics/)
-   -   EP Anomalies and Discretisation Schemes (https://www.cfd-online.com/Forums/phoenics/51874-ep-anomalies-discretisation-schemes.html)

Kevin Chow March 22, 2005 09:15

EP Anomalies and Discretisation Schemes
 
Hi chaps,

I'm simulating large (VERY LARGE) tank flows, which consists of a high velocity jet (~40-50m/s) blowing into a tank of around 25,000m^3 volume. To sort out turbulence, the bog-standard K-Epsilon model is used.

Sometimes after solution, the calculation of EP shows an extremely small minimum or extremely large maximum somewhere in the grid; for some reason, a single cell in the flowfield has a completely out-of-bounds EP value, yet this cell doesn't seem to affect anything. Other times, the length scale (and thus ENUT) seems to be very strangely computed, leading to large fields of large length scale and ENUT. This behaviour seems to be more predominant using 2nd order discretisation schemes, but using the hybrid/UDS still results in some strange single-cell EP results (unsurprisingly, EP doesn't converge well at all, and I get the feeling that if this behaviour isn't halted after more than the first 5 iterations, it gets increasingly worse with no recourse).

Can anyone shed any light or possible ways around this?

PattiMichelle April 3, 2005 15:24

Re: EP Anomalies and Discretisation Schemes
 
Do you have a pressure relief in the domain? What's the domain fluid setting?

Kevin Chow April 7, 2005 06:36

Re: EP Anomalies and Discretisation Schemes
 
The domain has a single fixed inlet and 3 other fixed pressure outlets; fluid is air (although I've done a bit of multi-phase work involving slightly heavier gases).

I'm thinking that setting a higher initial condition might stop the formation of extremely low EP (and KE, of order 1e-5 to -10) but even if one out of bounds cell does appear, it seems to propagate, even using a scheme like Van Leer.

Patti Sheaffer April 7, 2005 16:01

Re: EP Anomalies and Discretisation Schemes
 
I've had problems with higher-order discreetization scheems in PHOENICS. It's amazing how unstable things can become! More relaxation and gridsize (bigger) seem to help often.

Kevin Chow April 7, 2005 18:05

Re: EP Anomalies and Discretisation Schemes
 
I was afraid you would mention those two things. I can under-relax the flow (0.01 on EP and KE) and usually it would work okay, but as soon as I go above some threshold in the relaxation it'll start to go funny, and I don't know why. It's actually quite irritating, since I'm running a 540-timestep simulation that will take around 22-26 days, so I don't really want to increase the mesh density or do anything that will take even more time.

Can anyone shed any more light on this phenomenon?

Patti Sheaffer April 7, 2005 18:31

Re: EP Anomalies and Discretisation Schemes
 
My current understanding is that the solution matrix becomes unstable and I have seen that happen in other types of simulations - the existence of a relaxation threshold. I believe it's sometimes very difficult to get around this problem. I believe PHOENICS should have some built-in code to pre-detect when relaxation settings will cause numerical instabilities.

Kevin Chow April 8, 2005 04:13

Re: EP Anomalies and Discretisation Schemes
 
This is just a thought, but do you think that you can limit the "damage" by setting much tighter VARMAX and VARMIN? I know it won't reduce the propagation of the anomalies (well, it might, I don't know) but it should reduce the damage caused.

nasreddine April 8, 2005 11:03

Re: EP Anomalies and Discretisation Schemes
 
discretisation k-epsilon

nasreddine May 10, 2005 08:19

k-epsilon model in turbulent flow
 
pouvez vous s'il vous plait m'envoyer 3 à 4 thèses et quelques articles qui ont pour sujet les écoulements turbulents leur modélisation en modèle k-epsilon ; discrétisation spatiale simple et simpler ; méthode des projections ; discrétisation des équations de dissipation et celle de l'énergie cinétique turbulente. voici mes coordonnées: tèl: 0021698374840 adresse: 2 rue de tozeur cité elward oued ellil 2021 tunis tunisie merci pour votre collaboration

timur May 10, 2005 08:46

Re: k-epsilon model in turbulent flow
 
Ты бы еще на руÑÑком напиÑал :)

Ali Safavi November 7, 2005 04:41

Re: k-epsilon model in turbulent flow
 
Hi I am MS student of Chemical engineering in " IRAN UNIVERSITY OF SCIENCE AND TECHNOLOGY" I need to know mor about "k-epsilon model in turbulent flow".Please help me .

Sincerely yours Ali Safavi MS Student Chem. Eng. Dept. Iran University of Science and Technology, Tehran, Iran

GAURANG September 24, 2006 05:16

Re: k-epsilon model in turbulent flow
 
details

wiercipieta September 18, 2007 17:15

Re: k-epsilon model in turbulent flow
 
I would like to know how about so-called shear stresses. How do they look like in that model.

tushar December 19, 2007 05:35

Re: k-epsilon model in turbulent flow *NM*
 


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