|
[Sponsors] |
March 13, 2020, 07:35 |
High viscosity ratio in fire simulation
|
#1 |
New Member
Daniel Herranz
Join Date: Feb 2019
Posts: 5
Rep Power: 7 |
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 |
|
March 13, 2020, 07:49 |
Time-step
|
#2 |
Senior Member
|
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. |
|
March 13, 2020, 08:14 |
|
#3 | |
New Member
Daniel Herranz
Join Date: Feb 2019
Posts: 5
Rep Power: 7 |
Quote:
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 |
||
March 13, 2020, 08:22 |
Growth Rate and Final Time-Step
|
#4 |
Senior Member
|
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. |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |