Function weber_number 

pub fn weber_number(
    density: f64,
    velocity: f64,
    length: f64,
    surface_tension: f64,
) -> f64

Weber number: We = ρv²L / σ (inertia vs surface tension)