“…Each point emission source f j (r, t) can be described through its emission rate Q j (t) and emission site r j , that is, f j (r, t) = Q j (t)δ(r − r j ), where δ(r−r j ) is the Dirac delta centered at r j ∈ D. The domain of function f k (r, t) is restricted to a line Γ k ⊂ D in the case of a line source, and to a two-dimensional set A l ⊂ D in the case of an area source f l (r, t). It is important to note that each linearly distributed source, as well as each source distributed over an area, can be approximated by the sum of point sources [16]. However, the formulation of the control problem does not require such transformation, because such details are part of the numerical scheme used to solve the dispersion model (1)- (9).…”