Version working for Poiseuille flow. Dimensionless beta modified and recovering max velocity = 1.5.