Function critical_velocity 

pub fn critical_velocity(
    t_e_kev: f64,
    a_i: f64,
    z_i: f64,
    z_eff: f64,
) -> FusionResult<f64>

Critical velocity where electron drag equals ion drag [m/s].