Cebeci Smith-Eddy viscosity
Hello! I am trying to implement the Cebeci-Smith model the model calculates the eddy viscosity(turbulent kinematic viscosity) by finding the first y_{m} in which (v_{t})inner=(v_{t})outer and assigning the (v{t})inner in y<ym and outer otherwise. I don't understand how can i assign a value in y, because in my program in matlab, i have to assign a value of v_{t} for each x, not y as the model says! I tried to do this by calculating both inner and outer for all the points in my mesh and then finding the first ym in which they are equal and assigning at first (v_{t})inner and then for y>ym, (v_{t})outer, i post my code:
for j=2:J+1 for i = 2: K+1 if emt_inner(i,j) - emt_outer(i,j) <0 && flag emt(i,j) = emt_inner(i,j); else flag=0; emt(i,j) = emt_outer(i,j); end if emt(i,j) <0 emt(i,j)=0; end end end I turn zero emt where has zero values.I am not sure if this is correct. Please if somebody has any thoughts or has implemented this model i would be pleased if he helps me because i have been stuck for a long time! |
Nobody can help?
|
All times are GMT -4. The time now is 12:14. |