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

High viscosity ratio in fire simulation

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 13, 2020, 07:35
Default High viscosity ratio in fire simulation
  #1
New Member
 
Daniel Herranz
Join Date: Feb 2019
Posts: 5
Rep Power: 7
Danhehe is on a distinguished road
Hello, I am doing right now a transient simulation of an underground fire. The fire starts inside a train and the species and heat has to move arround the whole station.


I need to simulate like 1200(s) and my mesh has about 400k elements, which are polyhedra with maximum skewness 0.54 and minimum orthogonal quality 0.2.


I have done an steady case without the fire, just with natural flow inside the station to start the transient case from this solution.


For my transient simulation, I have set up a CFL=10 and an initial time step of 5e-4 s. This time step grows continualy by a factor of 2 every time step. It reaches something like 0.15 (meaning a CFL=7) and then, suddenly, at second 5 or 6 of simulation, viscosity ratio exceeds 10^5 and time step starts to fall.


My boundary conditions are ok as I checked them with my partners and the values are right.



Which could be my problem? Mesh? Schemes?...

Thanks
Danhehe is offline   Reply With Quote

Old   March 13, 2020, 07:49
Default Time-step
  #2
Senior Member
 
vinerm's Avatar
 
Vinerm
Join Date: Jun 2009
Location: Nederland
Posts: 2,946
Blog Entries: 1
Rep Power: 36
vinerm will become famous soon enough
Fire spread rate is usually quite high. Hence, even with a very coarse mesh, 0.15 s is quite large value for time-step. You have to use a smaller time-step. You can certainly increase the time-step by a factor of 2 but doing that every time-step is not a good idea. Use an increase factor of 1.1 or maximum 1.2 if you want it to increase with each time-step. A reduction factor of 0.9 would do good. Fix the maximum value such that the CFL does not go beyond 10.
__________________
Regards,
Vinerm

PM to be used if and only if you do not want something to be shared publicly. PM is considered to be of the least priority.
vinerm is offline   Reply With Quote

Old   March 13, 2020, 08:14
Default
  #3
New Member
 
Daniel Herranz
Join Date: Feb 2019
Posts: 5
Rep Power: 7
Danhehe is on a distinguished road
Quote:
Originally Posted by vinerm View Post
Fire spread rate is usually quite high. Hence, even with a very coarse mesh, 0.15 s is quite large value for time-step. You have to use a smaller time-step. You can certainly increase the time-step by a factor of 2 but doing that every time-step is not a good idea. Use an increase factor of 1.1 or maximum 1.2 if you want it to increase with each time-step. A reduction factor of 0.9 would do good. Fix the maximum value such that the CFL does not go beyond 10.

So you are basically saying that my maximum time steps has to be in the order of 10^-2? Or are you saying that a growth rate of 2 for time step is too high?


I do not think the growth rate is the problem,as it goes to 0.15 and stay there for more than 6 or 7steps, I think it is more related to 0.15 being too high, which is a real problem for me.


Thanks for your answer
Danhehe is offline   Reply With Quote

Old   March 13, 2020, 08:22
Default Growth Rate and Final Time-Step
  #4
Senior Member
 
vinerm's Avatar
 
Vinerm
Join Date: Jun 2009
Location: Nederland
Posts: 2,946
Blog Entries: 1
Rep Power: 36
vinerm will become famous soon enough
If the objective is not to observe the initial spreading rate but the spreading extent over 1200 s, then the growth rate is not that important. However, if that is of concern, then a growth rate of 2 is rather high. Final time-step is certainly high and that's what causes the trouble with numerics.
__________________
Regards,
Vinerm

PM to be used if and only if you do not want something to be shared publicly. PM is considered to be of the least priority.
vinerm 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
Turbulence Viscosity Ratio issue in EXtended domain with Pressure Far Field BCs Muneeb FLUENT 0 December 6, 2018 15:48
Erroneous eddy viscosity ratio for pipe flow preis OpenFOAM Running, Solving & CFD 1 May 11, 2018 19:58
3D CFD simulation of propeller blade - Turbulent viscosity limited to viscosity ratio Quentin_C FLUENT 1 April 4, 2017 05:32
High curvature 90 degree bend using SST, eddy viscosity drop?? tinab0binuh Main CFD Forum 0 August 22, 2014 17:26
turbulent viscosity ratio on a large duct newbie FLUENT 4 March 7, 2008 16:05


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