\(\frac{\partial}{\partial t}\left(\theta C_j\right)-\nabla\bullet\left(\theta D\nabla C_j\right)+\nabla\bullet\left(\mathbf{q}C_j\right)=R_{j}\)