Correct answer is (a) Gauss divergence theorem
To explain: The general diffusion term is div(Γ gradΦ). Integrating for the finite volume method, it becomes
∫CV div(Γ gradΦ)dV
Applying the Gauss divergence theorem,
∫A\(\vec{n}.\)(Γ gradΦ)dA
This is the boundary based integration as the boundaries will be areas.