Smagorinsky coefficient average on planes parallel to walls
 

Hi,all:
In channel flow LES, the most common way to use dynamic LES is to average Smagorinsky coefficient along the plane parallel to the wall. As OF's structured meshes are not described as i,j,k , it's not so convenient to do this kind of average. In postChannel utility, a typical method to do the average is to set a front face and walk along the normal direction to do the average of each mesh region, so the easiest way to do is like this:
1/ output the cD and cI as OF standard file
2/ call postChannel
3/ set the cD as a function of y, and match the averaged values to all the meshes

My question is： how can I define cD as a function as y and match cD to the mesh at specific y location?
(ps: this method is not efficient at all, I'll appreciate if there is more easier way)
Xianbei

Hi,all:
I think there is no need to write cD to a file, just do the collapse and then match cD at each y location. I'm wondering how can I implement this kind of match? Anyone can give me some advices?
ps:using 'if' is not a efficient way, say: if(abs(y-y_(i))<VSMALL) cD=cD(i), I think this will take a lot of time
Best regards
Xianbei