A divergência de alguns campos podem ser adaptadas utilizando o delta de dirac, de forma que o teorema da divergência não é violado. Assista a aula sobre função delta de dirac a alguns videos atrás

Inclusive, por causa disso, nos seus resultados, o fluxo total envolvendo uma carga pontual é a carga/(epsilon 0). Já segundo meu material didático, é igual à carga total...

A densidade de fluxo elétrico é uma quantidade escalar, definida como o módulo do campo elétrico gerado por uma carga elétrica livre, e corresponde, na visão de Faraday, ao número de linhas de campo que atravessam uma área determinada. Dito de outra forma, é a quantidade de fluxo que passa por uma área definida que seja perpendicular à direção do campo. Não há, portanto, incompatibilidade entre esta formulação e a que é apresentada nesta área. Nossa formulação é mais sofisticada matematicamente e define o fluxo como a integral de superfície do produto escalar entre o vetor campo elétrico e o vetor área diferencial.

Isto não tem nada a ver com a noção de densidade de fluxo, e sim com as unidades utilizadas. No sistema SI aparece o fator epsilon 0, no CGS não aparece.