Physical Model:
v = 4.0*Q/(pow(d, 2.0)*π)
λ = pow(1.0/(-2.0*log10(2.51/(Re*sqrt(λ)) + (1/3.71)*(k/d))), 2.0)
hr = (λ*L*pow(v, 2.0))/(2.0*d*G)