divergence free flux is modified. alphaEqn is modified