Function drag_force 

pub fn drag_force(
    fluid_density: f32,
    velocity: f32,
    drag_coeff: f32,
    area: f32,
) -> f32

Drag force: F_d = 0.5 * rho * v^2 * Cd * A