\vec{D}\cdot grad\varphi=div(\varphi\vec{D})-\varphi div\vec{D}=div(\varphi\vec{D})-\varphi\rho