CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > FLUENT

Epsilon residual in FLUENT

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 13, 2018, 13:44
Default Epsilon residual in FLUENT
  #1
New Member
 
Esmaeil Jabbari
Join Date: Dec 2014
Posts: 24
Rep Power: 11
esi1520 is on a distinguished road
Hi Dear Friends,

I have a problem for convergence of turbulence in fluent. In my model mesh quality and aspect ration are 0.27 and 16.1 respectively and use k-epsilon model for turbulence. Epsilon residual is not stable and after some iteration and adding different model like energy, radiation and etc for simulation, epsilon residual rise and i can not control it. Also i use first-order scheme for discretization of trubulaence model. Can anyone help me to set up fluent solver like Multigrid type, under relaxation factor and other way for solution epsilon equation and good simulation.

The residuals graph attached.

Thanks,,,
Attached Images
File Type: jpg Untitled.jpg (130.8 KB, 101 views)
esi1520 is offline   Reply With Quote

Old   May 13, 2018, 14:39
Default
  #2
Senior Member
 
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,667
Rep Power: 65
LuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura about
50 iterations is too few to make any judgement. Does it diverge? If not, keep cranking. When you change models, expect residuals to spike. Don't touch multigrid parameters. If it diverges you can try playing with urf's but you have to know why it diverged. If it doesn't diverge, leave it alone.
LuckyTran is offline   Reply With Quote

Old   December 9, 2019, 03:51
Default
  #3
Member
 
samm
Join Date: Feb 2019
Location: South Korea
Posts: 39
Rep Power: 7
Jegan is on a distinguished road
Hello,

I'm also facing the same situation. Epsilon oscillates too much even though i went through large number of iterations.
Any suggestions to get smooth residuals.

Thanks in Advance

Samm.
Attached Images
File Type: jpg epsiol_oscillates.jpg (92.3 KB, 61 views)
Jegan is offline   Reply With Quote

Old   December 9, 2019, 10:47
Default
  #4
Senior Member
 
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,667
Rep Power: 65
LuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura about
In your case, all your residuals are oscillating. It's not just the epsilon equation. Your solution can be improved everywhere.


There is a warning that the pdf enthalpy is too low, lower than your table, and this means your solution is being clipped. It's only 2 cells so maybe it's not important. Check if your solution makes sense (that you have a flame where you expect it to be). Probably, the solution is completely wrong. If it looks okay... then it's just a matter of better meshing.
LuckyTran is offline   Reply With Quote

Old   December 9, 2019, 10:59
Unhappy
  #5
Member
 
samm
Join Date: Feb 2019
Location: South Korea
Posts: 39
Rep Power: 7
Jegan is on a distinguished road
Thanks for your reply

Yes,You are correct .Still I couldn't get the flame that i expected.
My problem is 2d axisymmetric pressure based k-epsilon model.
I'm using Non-premixed chemical equilibrium model with 12 species

Methane Velocity inlet 4.62 m/s with pressure 4.113
Oxygen Velocity inlet 59.89 m/s with pressure 4.89
temperature at both inlet 290 K.

Hybrid Initialization

My case seems simple but i don't know where i'm making mistake
please help with your valuable suggestions.
Jegan is offline   Reply With Quote

Old   December 9, 2019, 13:03
Default
  #6
Senior Member
 
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,667
Rep Power: 65
LuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura about
If you don't get the flame you expect, then your solution is just wrong. It's meaningless to look at or compare residuals. Your residuals look bounded asymptotic, so it's more-or-less reached iterative convergence. You have to think about why the flame is not where it is. What's wrong with the flame? Do you even have a flame or did it blow out? etc. It could be initialization and/or wrong BC's and/or settings. I see you chose hybrid initialization. Well, a lot can go wrong.
LuckyTran is offline   Reply With Quote

Old   December 10, 2019, 00:04
Unhappy
  #7
Member
 
samm
Join Date: Feb 2019
Location: South Korea
Posts: 39
Rep Power: 7
Jegan is on a distinguished road
Thank you sir

My boundary conditions seems well as it was taken from the experimental work.so maybe no problem with that. Most probably initialization issue, if i use standard initialization,will it be good or something else i have to change.

Currently in this case residuals oscillating too much and flame seems completely not OK..no flame in the chamber,near inlet it was looking less temperature complete domain is same max temp say(3050 k). this may imbalance stuck in my domain or i don't know how to sort it out.

help me with some suggestions
Jegan is offline   Reply With Quote

Old   December 10, 2019, 10:36
Default
  #8
Senior Member
 
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,667
Rep Power: 65
LuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura about
Hybrid initialization is nice and automatic and in most cases but stupid.


Try either:
1) Standard initialization. You have to do some thinking to figure out the initial conditions.

2) Patch the progress variable in the combustion chamber. This forces there to be a flame and hopefully it anchors itself in the right place. You can also patch the temperature where it doesn't make sense. This is more-or-less the same as standard initialization, but you get to keep your presumably nice flow field that you've already achieved over 100k iterations.
LuckyTran is offline   Reply With Quote

Old   December 10, 2019, 23:34
Default
  #9
Member
 
samm
Join Date: Feb 2019
Location: South Korea
Posts: 39
Rep Power: 7
Jegan is on a distinguished road
Thanks again

I started the case again with standard initialization (atmospheric condition(pressure and temp) and other values zero). But i couldn't get the desired flame.

I attached residuals and temperature contour.

how can i tune this case in right way.Help me
Attached Images
File Type: jpg residualserrro.jpg (106.1 KB, 34 views)
File Type: jpg temperature_error.jpg (106.5 KB, 27 views)
Jegan is offline   Reply With Quote

Old   December 12, 2019, 09:01
Smile
  #10
New Member
 
SANKET BHATT
Join Date: Nov 2019
Posts: 10
Rep Power: 6
sanket2309 is on a distinguished road
Hi Jegan,

Use FMG initialization after standard initialization. Hopefully it will work. Give initial value of epsilon wisely in standard initialization. For achieve that value you have to do run model for 3-4 times. But, you will get experience from that.

Thanks.

Regards,
Sanket
sanket2309 is offline   Reply With Quote

Old   December 12, 2019, 09:22
Unhappy
  #11
Member
 
samm
Join Date: Feb 2019
Location: South Korea
Posts: 39
Rep Power: 7
Jegan is on a distinguished road
Thank you for your kind response

I think FMG initialization is unavailable for Non-premixed combustion
In standard initialization i simply gave atmospheric condition (pressure and temperature) i have no idea about other values so i put zero for all...now the case achieved the same flame as i attached in my previous reply.
Completely helpless situation now ..don't know what to alter to get on right way.
Give your valuable suggestions

Thanks
Jegan is offline   Reply With Quote

Old   December 25, 2019, 02:05
Default
  #12
Member
 
samm
Join Date: Feb 2019
Location: South Korea
Posts: 39
Rep Power: 7
Jegan is on a distinguished road
Any help please..why the heavy imbalance occurs. Unable to achieve the desired flame. constant maximum temperature after some iterations how to sort out this issue. where and how to resolve this problem
Jegan is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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
chtMultiRegionSimpleFoam: maximum number of iterations excedeed. Nkl OpenFOAM Running, Solving & CFD 19 October 10, 2019 02:42
Suppress twoPhaseEulerFoam energy AlmostSurelyRob OpenFOAM Running, Solving & CFD 33 September 25, 2018 17:45
chtMultiRegionSimpleFoam turbulent case Aditya Patil OpenFOAM Running, Solving & CFD 6 April 24, 2017 22:13
Simulation seems to converge but crashes suddenly xxxx OpenFOAM 16 September 12, 2014 08:07
pisoFoam with k-epsilon turb blows up - Some questions Heroic OpenFOAM Running, Solving & CFD 26 December 17, 2012 03:34


All times are GMT -4. The time now is 22:26.