U \frac{∂c}{∂x} + V \frac{∂c}{∂y} − ∇ · (k∇c) + ac = f  (1)^{k = [ \begin{matrix}k11 & k12 \\ k12 & k22 \end{matrix} ]  (2)}^{\frac{∂u}{∂t} + u\frac{∂u}{∂x} + v \frac{∂u}{∂y} = −\frac{∂p}{∂x} + \frac{1}{Re}∇^{2}u  (3)}^{\frac{∂v}{∂t} + u\frac{∂v}{∂x} + v \frac{∂v}{∂y} = −\frac{∂p}{∂y} + \frac{1}{Re}∇^{2}v  (4)}^{∇ · u⃗ = 0  (5)}

1