|
[Sponsors] |
One probable cause of high bounding epsilon & k values |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
April 9, 2015, 05:02 |
One probable cause of high bounding epsilon & k values
|
#1 |
Member
|
Hello friends,
I would like to share some experience regarding the problems which may or may not occur when running OpenFOAM multiphase solvers on sHM meshes with activated turbulence models. This isn't necessarily limited to multiphase solvers but can also occur when using single-phase solvers like simpleFoam etc. Here are some boundary conditions to reproduce the problem
As far as the question about numerical schemes goes I can say: it doesn't matter! I tested every scheme combination there is:
I tested the same case on both meshes as seen below. The result was that the solver always crashed with the left mesh with "bounding" warnings at each time step and the case with the right mesh continued to run, with only one "bounding" warning which was at the beginning of the simulation. As far as I am concerned, the only solution to the described problem is to use one cell size for the complete domain. I am open to other tips and I hope that this post helps someone. Regards, Sebastian |
|
April 9, 2015, 08:23 |
|
#2 |
Member
Join Date: Sep 2014
Location: Germany
Posts: 88
Rep Power: 11 |
Hi Sebastian,
have you checked y+ in the first near wall cell layer? The standard wall functions don't like a y+< 30 to 50. I observed very strange behaviour in the near wall region in cases where the near wall resolution was too high. Do you use wall functions at all? Just a few thoughts about your problem, maybe you have already considered these points. Last edited by TobM; April 9, 2015 at 08:24. Reason: typo |
|
April 9, 2015, 08:37 |
|
#3 | |
Member
|
Quote:
It's acceptable for the air to use the Low-Re wall treatment, since the kinematic viscosity is high. When considering the fluids (liquid metals & slag in my case) the kinematic viscosity drops dramatically and therefore the y+ is rather high. Since I don't assume the gas velocities to be in the range of the liquids I used wall functions (High-Re wall treatment). I guess there is no way to stay out of a specific range in my case. |
||
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Calculation of k and epsilon freezes | Nigirim | OpenFOAM Running, Solving & CFD | 1 | November 14, 2012 08:52 |
epsilon and K blowing up. | sivakumar | OpenFOAM Running, Solving & CFD | 1 | October 25, 2012 05:50 |
Values of epsilon | Gearb0x | OpenFOAM Running, Solving & CFD | 2 | May 20, 2010 13:36 |
Bounding epsilon and K with rasInterFoam | openfoam_user | OpenFOAM Running, Solving & CFD | 0 | October 23, 2008 09:48 |
Multicomponent fluid | Andrea | CFX | 2 | October 11, 2004 06:12 |