Function bernoulli_pressure_kernel 

pub fn bernoulli_pressure_kernel<R>(
    p1: &Pressure<R>,
    v1: &Speed<R>,
    h1: &Length<R>,
    v2: &Speed<R>,
    h2: &Length<R>,
    density: &Density<R>,
) -> Result<Pressure<R>, PhysicsError>
where
    R: RealField + FromPrimitive,

Calculates pressure $P_2$ using Bernoulli’s principle.

$P_1 + \frac{1}{2}\rho v_1^2 + \rho g h_1 = P_2 + \frac{1}{2}\rho v_2^2 + \rho g h_2$

Solves for $P_2$: $P_2 = P_1 + \frac{1}{2}\rho(v_1^2 - v_2^2) + \rho g(h_1 - h_2)$