Positive and negative - Measured in $Coulombs~(C)$
Electron charge:$~e=-1.6021765\times {{10}^{-19}}~C$
Universe is neutral - Default state of a material body is neutral
Point charge (infinitesimally small volume)
Line distribution – linear density ${{\rho }_{\ell }}$ $C/m$
Surface distribution – surface density ${{\rho }_{s}}$ $C/m^{2} $
Volume distribution – volume density ${{\rho }_\text{v}}$ $C/m^{3}$
Flux equals charge producing the Flux $(Coulombs)$
Emerges out of Positive Charge (source)
Lands on Negative Charge (sink)
${{\psi }_{Closed~surface}}=Total~{{Q}_{enclosed~by~the~surface}}$
$D={{\left\{ d\psi /dS \right\}}_{max}}=d\psi /\left( d{{S}_{normal}} \right)~$$C/m^{2}$
Vector direction along flux lines
$d{{\psi }_{dS}}=\vec{D}\cdot \overrightarrow{dS}~~~and~~~{{\psi }_{s}}=\iint_{s}{d\psi }=\iint_{s}{\vec{D}\cdot \overrightarrow{dS}}$
$\left\{ {{\psi }_{closed~surface~S}} \right\}=\unicode{x222F}_s~
\vec{D}\cdot \overrightarrow{dS} =\iiint_{enclosed~volume}{{{\rho }_\text{v}}d\text{v}}$
$\vec{\nabla }\cdot \vec{D}={{\rho }_\text {v}}\left( position \right)$
$\vec{\nabla }\cdot \vec{D}=\left[ \frac{\partial }{\partial x}{{D}_{x}}+\frac{\partial }{\partial y}{{D}_{y}}+\frac{\partial }{\partial z}{{D}_{z}} \right]$
$\vec{\nabla }\cdot \vec{D}=\left[ \frac{1}{\rho }\frac{\partial }{\partial \rho }\left( \rho {{D}_{\rho }} \right)+\frac{1}{\rho }\frac{\partial }{\partial \varphi }{{D}_{\varphi }}+\frac{\partial }{\partial z}{{D}_{z}} \right]$
$\vec{\nabla }\cdot \vec{D}=\left[ \frac{1}{{{r}^{2}}}\frac{\partial }{\partial r}\left( {{r}^{2}}{{D}_{r}} \right)+\frac{1}{r~sin\theta }\frac{\partial }{\partial \theta }\left( sin\theta ~{{D}_{\theta }} \right)+\frac{1}{r~sin\theta }\frac{\partial }{\partial \varphi }{{D}_{\varphi }} \right]$
$\unicode{x222F}_s
~\vec{D}\cdot \overrightarrow{dS} =\iiint_{volume~enclosed~by~S}{\vec{\nabla }\cdot \vec{D}~d \text {v}}$
$\unicode{x222F}_s
~\vec{A}\cdot \overrightarrow{dS} =\iiint_{volume~enclosed~by~S}{\vec{\nabla }\cdot \vec{A}~d \text {v}}$
Total enclosed in a volume, given flux density distribution
${{Q}_{vol}}=\unicode{x222F}_{S~enclosing~vol}
\vec{D}\cdot \overrightarrow{dS} $