В оригинале у Зорича формула выглядит иначе:

У Зорича (Формула (6') рядом с (12)) определение дивергенции выглядит аналогично:

и там также применимы ваши утверждения
1) Явно указано, что к нулю стремится не просто величина объёма

а диаметр этого объёма

2) Точка, в которой берётся градиент (здесь обозначенная

) не входит как аргумент в подынтегральную функцию, а входит как параметр в выбираемый объём

таким образом, что этот объём содержит эту точку.
Однако используя предложенный вами метод (см цитату из топика
topic77894.html ниже) для

где

,
который я и скопировал для градиента, получаем верный результат для дивергенции:

.
Смотрите
topic77894.html Цитирую вас:
"Когда

направлен только по радиусу, то он имеет вид

где

- произвольная функция.
Теперь, зная, что

можно взять тонкий сферический слой

и проинтегрировать по нему дивергенцию:

Теперь, поскольку всё сферически-симметрично, разделив на площадь сферы, можно найти и саму

".
У вас

(мое

) стремиться к нулю, а не диаметр этого объема (

или

).
Почему для дивергенции ваш подход дает правильное выражение, а мое повторение вашего подхода для градиента дает ошибку, хотя определения градиента и дивергенции у Зорича аналогичны.
Может неверен ваш метод для дивергенции.
Тогда как его следует изменить, чтобы получать более корректно дивергенцию из теоремы Гаусса для

, а за одно и для градиента ?