CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > Siemens > STAR-CCM+

LES simulation - Not able to bring down M value of every cell below 0.1

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

LinkBack Thread Tools Search this Thread Display Modes
Old   February 5, 2018, 14:08
Default LES simulation - Not able to bring down M value of every cell below 0.1
New Member
Join Date: Feb 2018
Posts: 1
Rep Power: 0
shettysds is on a distinguished road
"We tried doing a LES simulation on Tundish (3.18 m*0.8 m*0.78 m). We are refining the cell based on ratio of sgs_tke and total_tke (M = (k_sgs)/(k_sgs + k_res)).

Converting the RANS simulation to a LES simulation we used a field function to refine the mesh at important areas "(${M} > 0.1) ? max(${cellSize}/40, 0.000002) :0" to find the cells with high values of M and splitting them by 40 times i.e. making the cells finer, which eventually resulted in billion cells. Even though cells were as small as 1.6e-5m, there were always some of these cells having M value closer to 1 or 1. when we back traced, we found out that these cells have resolved TKE(k_res) value =0. which basically means M = 1 = (k_sgs)/(k_sgs + k_res) . So no matter how much finer we make the mesh, these cells will always have M value of 1.

Our question would be; Is it fine to have some cells with value of M=1? if yes, then what percentage of these cells are ideal?"

i am getting k_res = 0 in these cells. thus making TKE = K_sgs. which would basically mean that these cells are badly resolved to have M =1. if we had some values for K_res instead of zero, we would have understood it is less than K_sgs so it has to be made finer. with k_res coming as zero at these cells, we are not able to understand why we are not getting any values for K_res at these cells. and if we can ignore these cells.
shettysds is offline   Reply With Quote

Old   February 7, 2018, 19:39
Senior Member
Join Date: Nov 2010
Location: USA
Posts: 1,232
Rep Power: 24
me3840 is on a distinguished road
According to Pope it's acceptable to have the total resolved energy be 80% of the total energy, or 0.2 by your metric. Whether or not you want to follow that is your choice.

However your case is suspicious. Are you using the hybrid upwind discretization? I'm not sure how it's implemented to compute the SGS KE, but the code may lock certain poor cells using the upwind scheme with the SGS, perhaps. Are those areas poor quality?
me3840 is offline   Reply With Quote


les model

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
[Other] refineWallLayer Error Yuby OpenFOAM Meshing & Mesh Conversion 2 November 11, 2021 11:04
Neighboring cells in tetrahedral mesh vishwesh OpenFOAM Programming & Development 9 November 10, 2017 07:06
How to use "translation" in solidBodyMotionFunction in OpenFOAM rupesh_w OpenFOAM Running, Solving & CFD 5 August 16, 2016 04:27
Fluent RPvar list for molar-weight pilou FLUENT 6 February 21, 2011 03:08
boundary condition scalarField in parallel nlc OpenFOAM 4 November 5, 2010 19:26

All times are GMT -4. The time now is 06:19.