The Turbulence Viscosity in the SPS method
Dear All,
Hi.
I have a question about the "prr" parameter in the JSPHCpu.cpp code.
As we know in SPS method the turbulence viscosity is calculated with use of the "rate of the strain tensor (2SS)" function which is named prr in codes and should be equal to 2*(0.5(dw/dx+du/dz))^2 in 2D condition. [S=0.5(dw/dx+du/dz)]
I could not find a relation between the equation of the prr in the code and its equations in the references. Could anyone please give me an explanation about the calculation of the prr in the DualSPHysics code with a reference (In a 3D condition obviously)?
Regrads,
Mojtaba
csph=velocity gradient matrix
const float pow1=csph.xx*csph.xx + csph.yy*csph.yy + csph.zz*csph.zz;
const float prr=pow1+pow1 + csph.xy*csph.xy + csph.xz*csph.xz + csph.yz*csph.yz;
Hi.
I have a question about the "prr" parameter in the JSPHCpu.cpp code.
As we know in SPS method the turbulence viscosity is calculated with use of the "rate of the strain tensor (2SS)" function which is named prr in codes and should be equal to 2*(0.5(dw/dx+du/dz))^2 in 2D condition. [S=0.5(dw/dx+du/dz)]
I could not find a relation between the equation of the prr in the code and its equations in the references. Could anyone please give me an explanation about the calculation of the prr in the DualSPHysics code with a reference (In a 3D condition obviously)?
Regrads,
Mojtaba
csph=velocity gradient matrix
const float pow1=csph.xx*csph.xx + csph.yy*csph.yy + csph.zz*csph.zz;
const float prr=pow1+pow1 + csph.xy*csph.xy + csph.xz*csph.xz + csph.yz*csph.yz;
Comments
We have checked this stuff several times, but we will do it again to be sure.
Regards
Thanks a lot.
I'm looking forward for your reply.
Regards.
Mojtaba
Thanks a lot.
I'm looking forward for your reply.
Regards.
Mojtaba