Turbulence viscosity SPS

edited February 26 in Old versions

Dear all,


I am trying to modify the DSPH code implementing a Dynamic Smagorinsky viscosity model.

I see in the literature that the kernel function serves as a grid filter. At the standard Smagorinsky model

in the DSPH code I see:

 const float pow1=gradvel.xx*gradvel.xx + gradvel.yy*gradvel.yy + gradvel.zz*gradvel.zz;

const float prr=pow1+pow1 + gradvel.xy*gradvel.xy + gradvel.xz*gradvel.xz + gradvel.yz*gradvel.yz;

Where is the influence of the kernel function in this? When I search definition for the gradvel.* I do not find the influence of the kernel.

And I need this, to take an influence of 2h for the test filter.

Any help would be appreciated. 

Sign In or Register to comment.