openfoam中冲蚀公式（erosion）的程序编写

[wzdrfy]

finnie模型，在erosion程序中是
// Finnie’s model
const scalar coeff = p.nParticle()*p.mass()*sqr(magU)/(p_psi_K_);
if (tan(alpha) < K_/6)
{
Q += coeff(sin(2alpha) - 6/K_sqr(sin(alpha)));
}
else
{
Q += coeff(K_*sqr(cos(alpha))/6);
}
但如果是HUSER模型

该如何编写呢？有没有大佬解答一下