“…We first deal with the weak solutions, which are known (see [1], [10]) to exist globally in time for all $a\in L_{\sigma}^{2}$ , satisfying the identity : $\langle u(t), \varphi\rangle=\langle e^{-tA}a, \varphi\rangle+\int_{0}^{t}\langle u\otimes u, \nabla e^{-(t-s)A}\varphi\rangle ds$ given in (i) is constructed via approximate solutions $\{u_{N}\}$ as given in [1], [7], [13], which satisfy $\lim_{Narrow\infty}\int_{0}^{\infty}||u_{N}(t)-u(t)||_{2}^{2}dt=0$ .…”