phi k_a*(1-epsilon)*c Local reaction rate alpha (eta/(Da*L^2))*q*(1-epsilon)/(q+epsilon) Drag-force coefficient