pub fn sph_kernel_viscosity_lap(r_len: f32, h: f32) -> f32
Viscosity kernel laplacian: W_lap = (15/(2πh^3))*(-(r^3/2h^3) + r^2/h^2 + h/(2r) - 1)