pub fn peclet_number(velocity: f64, length: f64, diffusivity: f64) -> f64
Peclet number: Pe = vL/α (advection vs diffusion)