Формула Остроградского-Гаусса 🌊

Представь себе резервуар с водой, внутри которого работают несколько труб: одни закачивают воду внутрь, другие откачивают её наружу. Снаружи резервуара ты видишь только его поверхность — стенки, через которые вода либо втекает, либо вытекает. Возникает естественный вопрос: можно ли, зная только то, что происходит на самой поверхности резервуара (сколько воды через неё проходит), сказать что-то содержательное о том, что творится внутри — сколько там суммарно работают источники и стоки? И наоборот: зная устройство всех труб внутри, можно ли предсказать суммарный поток через стенки, даже не имея понятия о форме самого резервуара?

Именно на этот вопрос отвечает формула Остроградского-Гаусса — одна из самых мощных теорем математического анализа, связывающая поверхность тела с его объёмом. Она утверждает, что суммарный поток векторного поля через замкнутую поверхность в точности равен тройному интегралу от особой характеристики этого поля — дивергенции — по объёму, который эта поверхность ограничивает. Грубо говоря: «сколько вытекает наружу» всегда равно «сколько производится внутри». Ты уже встречал похожую идею на плоскости — формула Грина связывала циркуляцию поля по замкнутому контуру с двойным интегралом от ротора внутри области. Формула Остроградского-Гаусса — это прямое трёхмерное обобщение той же самой идеи, только вместо контура и площади здесь работают поверхность и объём.

Зачем вообще нужна эта теорема, если её можно принять просто как красивое математическое совпадение? Дело в том, что она даёт совершенно практический вычислительный трюк: если тебе нужно посчитать поток поля через какую-то жутко изогнутую замкнутую поверхность (представь границу турбины или сложную деталь двигателя), прямое вычисление поверхностного интеграла может оказаться кошмаром — придётся аккуратно параметризовать каждый кусок этой поверхности. А тройной интеграл по объёму внутри неё часто вычисляется значительно проще, особенно если дивергенция поля — простая функция координат. И работает эта замена в обе стороны: иногда, наоборот, выгоднее заменить трудный тройной интеграл на поток через удобную поверхность.

А теперь — обещанная связь с искусственным интеллектом, ради которой эта тема особенно важна тем, кто движется в сторону машинного обучения и анализа данных. В современных генеративных моделях (так называемых непрерывных нормализующих потоках, лежащих в основе некоторых диффузионных моделей и нейросетевых дифференциальных уравнений) плотность вероятности данных «течёт» вдоль векторного поля скорости точно так же, как течёт жидкость в резервуаре. И скорость изменения логарифма плотности вероятности вдоль такого потока задаётся именно дивергенцией этого векторного поля — прямое, буквальное применение идеи, которую ты сегодня разберёшь. Более того, сама логика «суммарный поток через границу равен суммарной характеристике внутри» — это глубокий родственник принципа нормировки вероятности: полная вероятность, «втекающая» и «вытекающая» из области значений случайной величины, обязана балансироваться, и уравнения, управляющие этим балансом (уравнение неразрывности, уравнение Фоккера-Планка), — прямые потомки формулы Остроградского-Гаусса.

История: откуда это взялось?

Как это часто бывает с фундаментальными результатами анализа, формула Остроградского-Гаусса не была открыта одним человеком в один момент — её элементы независимо всплывали у нескольких математиков на протяжении полувека. Первым к похожему результату, хотя и в частном случае, подошёл Жозеф Луи Лагранж ещё в 1762 году, изучая задачи гидродинамики. Но именно Карл Фридрих Гаусс, разрабатывая математический аппарат электростатики в 1813 году, впервые применил теорему в узнаваемом виде — для вычисления потока электрического поля через поверхности, окружающие заряды, что впоследствии легло в основу закона Гаусса.

Строгую и полную формулировку теоремы для произвольного векторного поля и произвольной замкнутой поверхности дал русский математик Михаил Васильевич Остроградский в 1826 году. Молодой Остроградский в это время учился в Париже у величайших математиков эпохи — Коши, Лапласа, Фурье и Пуассона, — и представил доказательство общего случая теоремы Парижской академии наук. Полностью его работа была опубликована лишь в 1831 году, но приоритет строгого общего доказательства именно за ним — не случайно во многих европейских и особенно русскоязычных учебниках теорема носит двойное имя. Забавно, что примерно в то же самое время, в 1828 году, английский самоучка Джордж Грин — тот самый, чьё имя носит формула из прошлого урока — независимо вывел двумерный частный случай той же самой идеи, даже не подозревая о работах континентальных коллег.

Судьба этой теоремы оказалась удивительно масштабной: она легла в основу не только классической электродинамики (уравнения Максвелла в интегральной форме — это буквально формула Остроградского-Гаусса, применённая к электрическому и магнитному полям), но и гидродинамики, теории теплопроводности, общей теории относительности и, спустя почти два столетия, современных генеративных моделей машинного обучения. Немногие результаты математики XIX века могут похвастаться настолько широким и до сих пор растущим кругом применений.

Точная формулировка: поток через поверхность и источники внутри

Интуиция: воздушный шар с вентиляторами

Давай разберёмся с самой сутью теоремы на интуитивном примере. Представь прозрачный воздушный шар, внутри которого установлено несколько маленьких вентиляторов — одни нагнетают воздух наружу (это источники), другие, наоборот, всасывают воздух внутрь (это стоки). Если сложить суммарную «производительность» всех вентиляторов внутри шара — с плюсом для тех, что выдувают воздух, и с минусом для тех, что всасывают, — результат должен в точности совпасть с тем, сколько воздуха реально проходит через оболочку шара наружу в единицу времени. Это интуитивно очевидно: воздух не может «взяться из ниоткуда» или «бесследно исчезнуть» внутри шара — весь избыток, произведённый вентиляторами, обязан вытечь через оболочку, а весь недостаток, созданный всасывающими вентиляторами, обязан прийти извне через ту же оболочку.

Теперь переведи эту картинку на язык математики. Роль «вентиляторов» в каждой точке пространства играет локальная характеристика поля, которая называется дивергенцией — она показывает, насколько сильно поле «расходится» из этой точки (положительная дивергенция — точка похожа на источник, отрицательная — на сток). А роль «суммарного потока через оболочку» играет поверхностный интеграл от поля по замкнутой поверхности. Формула Остроградского-Гаусса — это математически строгая версия утверждения «сумма всех источников внутри равна потоку через оболочку».

Теорема

Теорема (Остроградского-Гаусса). Пусть $V$ — ограниченная замкнутая область в пространстве $\mathbb{R}^3$ с кусочно-гладкой границей $\partial V$, а векторное поле $\vec F(x,y,z)=\bigl(P(x,y,z),\,Q(x,y,z),\,R(x,y,z)\bigr)$ непрерывно дифференцируемо (то есть класса $C^1$) во всех точках замкнутой области $V$, включая границу. Тогда суммарный поток поля через границу $\partial V$, взятую с внешней нормалью $\vec n$, равен тройному интегралу от дивергенции поля по объёму $V$:

$$\iint_{\partial V} \vec F\cdot\vec n\,dS=\iiint_V \operatorname{div}\vec F\,dV$$

где дивергенция поля в декартовых координатах вычисляется как

$$\operatorname{div}\vec F=\frac{\partial P}{\partial x}+\frac{\partial Q}{\partial y}+\frac{\partial R}{\partial z}$$

Условия применимости важно держать в голове: поверхность $\partial V$ обязана быть замкнутой (без «дыр» и краёв) и кусочно-гладкой (то есть состоящей из конечного числа гладких кусков, как грани куба), а компоненты поля — непрерывно дифференцируемыми во всей замкнутой области, а не только внутри неё. Нормаль $\vec n$ всегда берётся внешней — направленной наружу из области $V$.

Физическая интерпретация обеих частей формулы стоит того, чтобы проговорить её отдельно. Левая часть, $\iint_{\partial V}\vec F\cdot\vec n\,dS$, — это суммарный поток поля через поверхность: если $\vec F$ — поле скоростей жидкости, это буквально объём жидкости, вытекающий через поверхность за единицу времени. Правая часть, $\iiint_V\operatorname{div}\vec F\,dV$, — это суммарная «мощность источников» внутри объёма: интеграл от локальной плотности источников (дивергенции) по всему объёму. Формула утверждает: эти две величины совпадают всегда, для любого поля класса $C^1$ и любой достаточно приличной замкнутой поверхности.

Разбор примеров

Пример 1 (лёгкий). Найти поток поля $\vec F=(x,y,z)$ через полную поверхность единичного куба $[0,1]^3$, вычислив обе части формулы независимо.

Начнём с правой части — тройного интеграла от дивергенции. Дивергенция поля: $\operatorname{div}\vec F=\dfrac{\partial x}{\partial x}+\dfrac{\partial y}{\partial y}+\dfrac{\partial z}{\partial z}=1+1+1=3$. Поскольку дивергенция постоянна, тройной интеграл равен просто дивергенции, умноженной на объём куба:

Теперь проверим левую часть напрямую, перебрав все шесть граней куба. На грани $x=1$ внешняя нормаль $\vec n=(1,0,0)$, и $\vec F\cdot\vec n=x=1$ (площадь грани равна $1$), вклад — $1$. Аналогично грани $y=1$ и $z=1$ дают по $1$ (там $\vec F\cdot\vec n=y=1$ и $\vec F\cdot\vec n=z=1$ соответственно). На противоположных гранях $x=0$, $y=0$, $z=0$ поле обнуляется в соответствующей компоненте (например, на грани $x=0$ имеем $\vec F=(0,y,z)$, а нормаль $\vec n=(-1,0,0)$, значит $\vec F\cdot\vec n=0$), так что эти три грани дают нулевой вклад.

Ответ: суммарный поток равен $1+1+1+0+0+0=3$, что в точности совпадает с тройным интегралом — теорема подтверждена.

Пример 2 (средний). Вычислить поток поля $\vec F=(x^3,y^3,z^3)$ через сферу радиуса $R$ с центром в начале координат.

Считать этот поток напрямую как поверхностный интеграл по сфере было бы утомительно, но формула Остроградского-Гаусса позволяет свести всё к тройному интегралу. Дивергенция поля: $\operatorname{div}\vec F=3x^2+3y^2+3z^2=3(x^2+y^2+z^2)=3r^2$, где $r$ — расстояние от начала координат (в сферических координатах). Переходя к сферическим координатам ($dV=r^2\sin\varphi\,dr\,d\varphi\,d\theta$), получаем:

Каждый из трёх множителей вычисляется отдельно: $\int_0^{2\pi}d\theta=2\pi$, $\int_0^\pi\sin\varphi\,d\varphi=2$, $\int_0^R r^4\,dr=\dfrac{R^5}5$. Перемножаем всё вместе:

Ответ: поток равен $\dfrac{12\pi R^5}5$. Обрати внимание, насколько это было проще, чем пытаться напрямую параметризовать сферу и интегрировать по ней $x^4+y^4+z^4$ — именно в такой замене и состоит главная практическая ценность теоремы.

Пример 3 (сложный, машинное обучение). В непрерывных нормализующих потоках точки данных перемещаются вдоль траекторий, задаваемых векторным полем скорости $\vec f(x,y,z)$, а плотность вероятности вдоль этого движения меняется по фундаментальному соотношению $\dfrac{d}{dt}\log p(x(t))=-\operatorname{div}\vec f\bigl(x(t)\bigr)$. Пусть поле скорости — простое «растяжение» $\vec f=(x,y,z)$, а начальная плотность вероятности равномерна и нормирована на единичном шаре радиуса $1$ (то есть $\iiint_{\text{шар}}\rho\,dV=1$). Найти суммарный поток «вероятностной массы» через границу шара в начальный момент времени и связать его с изменением логарифма плотности.

Дивергенция поля скорости: $\operatorname{div}\vec f=\dfrac{\partial x}{\partial x}+\dfrac{\partial y}{\partial y}+\dfrac{\partial z}{\partial z}=3$ — постоянна в любой точке пространства. По формуле Остроградского-Гаусса, поток «вероятностной массы» $\rho\vec f$ через границу шара равен тройному интегралу от $\operatorname{div}(\rho\vec f)=\rho\operatorname{div}\vec f=3\rho$ (поскольку $\operatorname{div}\vec f$ — константа, а $\rho$ можно вынести правило произведения дивергенции):

Ответ: суммарный отток вероятностной массы через границу шара равен $3$, что в точности совпадает со значением дивергенции поля скорости. Это не совпадение: для постоянной дивергенции скорость изменения логарифма плотности $\dfrac{d}{dt}\log p=-\operatorname{div}\vec f=-3$ в каждой точке, а теорема Остроградского-Гаусса — это ровно тот математический механизм, который переводит эту локальную, точечную характеристику (дивергенцию) в глобальный, интегральный факт (суммарный поток через границу конечной области). Именно поэтому дивергенция векторного поля скорости — центральная величина при обучении непрерывного нормализующего потока: она напрямую говорит, как быстро «сжимается» или «расширяется» облако вероятности вдоль траектории.

Почему это важно

Формула Остроградского-Гаусса — это, по сути, закон сохранения в самой общей математической форме: то, что происходит на границе области, полностью определяется тем, что происходит внутри неё, и наоборот. В физике эта идея лежит в основе всех уравнений сохранения — массы, заряда, энергии. В машинном обучении та же самая логика управляет тем, как изменяется плотность вероятности вдоль потоков в генеративных моделях: дивергенция векторного поля скорости — это точечная, локальная версия того самого «суммарного потока через границу», о котором говорит теорема. Без ясного понимания этой связи разобраться в математике современных генеративных моделей практически невозможно.

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

Интуиция: кубики, которые обмениваются потоками

А теперь давай разберёмся, почему теорема вообще верна, а не просто «работает на конкретных примерах». Представь, что ты разбил весь объём $V$ на множество очень маленьких кубиков. Для каждого отдельного кубика поток через его шесть граней приближённо равен произведению дивергенции поля в этом кубике на его объём — это, по сути, локальное определение дивергенции (насколько сильно поле «растекается» из бесконечно малого объёма).

Теперь сложи вклады всех кубиков. Ключевое наблюдение: у любых двух соседних кубиков есть общая грань, и поток через эту общую грань считается дважды — один раз как выходящий поток для одного кубика (с нормалью, направленной в одну сторону), и второй раз как входящий поток для соседнего кубика (с нормалью, направленной ровно в противоположную сторону). Поскольку поле непрерывно, эти два вклада почти в точности равны по модулю и противоположны по знаку — они взаимно сокращаются. Единственные вклады, которые не с кем сокращать, — это грани кубиков, лежащие на самой внешней границе всей области $V$. Значит, сумма всех локальных «источников» (дивергенций, умноженных на объёмы маленьких кубиков) в пределе равна ровно суммарному потоку через внешнюю границу — а это и есть формула Остроградского-Гаусса.

Теорема

Схема строгого доказательства (для простой области). Пусть область $V$ можно описать как «$z$-простую»: для каждой точки $(x,y)$ проекции области $V$ на плоскость $xy$ значение $z$ пробегает отрезок от $z=z_1(x,y)$ (нижняя граница) до $z=z_2(x,y)$ (верхняя граница). Тогда для компоненты поля $R(x,y,z)$ (третьей координаты $\vec F=(P,Q,R)$) справедливо:

$$\iiint_V\frac{\partial R}{\partial z}\,dV=\iint_D\left[\int_{z_1(x,y)}^{z_2(x,y)}\frac{\partial R}{\partial z}\,dz\right]dA=\iint_D\Bigl[R\bigl(x,y,z_2(x,y)\bigr)-R\bigl(x,y,z_1(x,y)\bigr)\Bigr]dA$$

где внутренний интеграл по $z$ вычислен по формуле Ньютона-Лейбница, а $D$ — проекция области на плоскость $xy$. Правая часть в точности равна потоку компоненты $(0,0,R)$ через верхнюю и нижнюю части границы $\partial V$ (с учётом знаков внешних нормалей). Аналогичные выкладки для компонент $P$ и $Q$ (при условии, что область также «$x$-простая» и «$y$-простая») дают в сумме полную формулу Остроградского-Гаусса. Для произвольной области доказательство завершается разбиением её на конечное число простых кусков, у которых внутренние границы (общие для соседних кусков) взаимно сокращаются — в точности так, как описано в интуиции выше.

Обрати внимание на структуру этого доказательства: оно не какая-то магия, а прямое следствие формулы Ньютона-Лейбница, применённой «послойно» — сначала интегрируем по одной переменной (получая разность значений на верхней и нижней границе), а затем по оставшимся двум переменным (что даёт поверхностный интеграл). Это тот же самый приём, который лежал в основе доказательства формулы Грина на плоскости, только на одно измерение выше.

Разбор примеров

Пример 1 (лёгкий). Проверить теорему прямым вычислением обеих частей для прямоугольного параллелепипеда $[0,1]\times[0,2]\times[0,3]$ и поля $\vec F=(x,0,0)$.

Дивергенция: $\operatorname{div}\vec F=\dfrac{\partial x}{\partial x}=1$. Тройной интеграл равен объёму параллелепипеда: $\iiint_V 1\,dV=1\cdot2\cdot3=6$.

Теперь поверхностный интеграл. Поле имеет только $x$-компоненту, значит на четырёх гранях, параллельных оси $x$ (то есть $y=0$, $y=2$, $z=0$, $z=3$), нормаль перпендикулярна полю, и вклад равен нулю. Остаются грани $x=0$ и $x=1$. На грани $x=0$: $\vec F=(0,0,0)$, вклад — $0$. На грани $x=1$: нормаль $\vec n=(1,0,0)$, $\vec F\cdot\vec n=x=1$, площадь грани $=2\cdot3=6$, вклад $=1\cdot6=6$.

Ответ: поверхностный интеграл равен $0+6=6$, что точно совпадает с тройным интегралом.

Пример 2 (средний). Вычислить поток поля $\vec F=(0,0,z^2)$ через полную поверхность цилиндра $x^2+y^2\le4$, $0\le z\le3$ двумя способами.

Через тройной интеграл. Дивергенция: $\operatorname{div}\vec F=\dfrac{\partial z^2}{\partial z}=2z$. В цилиндрических координатах:

Через поверхностный интеграл. Боковая поверхность цилиндра имеет нормаль, лежащую в горизонтальной плоскости (её $z$-компонента равна нулю), а поле направлено строго вдоль оси $z$, так что $\vec F\cdot\vec n=0$ на всей боковой поверхности. На нижнем основании $z=0$: поле $\vec F=(0,0,0)$, вклад — $0$. На верхнем основании $z=3$: нормаль $\vec n=(0,0,1)$, $\vec F\cdot\vec n=z^2=9$, площадь основания $=4\pi$, вклад $=9\cdot4\pi=36\pi$.

Ответ: оба способа дают $36\pi\approx113{,}097$ — теорема подтверждена, причём поверхностный интеграл оказался даже проще: не нужно было ничего вычислять на боковой поверхности вовсе.

Пример 3 (сложный, машинное обучение). Проиллюстрируем идею доказательства на двух соседних кубиках — ровно та картина, которая лежит в основе метода накопления локальных эффектов в непрерывных моделях. Рассмотрим объединение кубов $[0,1]\times[0,1]\times[0,1]$ и $[1,2]\times[0,1]\times[0,1]$ (вместе образующих параллелепипед $[0,2]\times[0,1]\times[0,1]$) с полем $\vec F=(x^2,0,0)$. Покажем явно, как поток через внутреннюю общую грань $x=1$ сокращается, оставляя только вклад внешней границы.

Дивергенция: $\operatorname{div}\vec F=2x$. Тройной интеграл по левому кубу: $\int_0^1 2x\,dx=\bigl[x^2\bigr]_0^1=1$. Тройной интеграл по правому кубу: $\int_1^2 2x\,dx=\bigl[x^2\bigr]_1^2=4-1=3$. Сумма по обоим кускам: $1+3=4$ — это должно совпасть с тройным интегралом по всему объединённому параллелепипеду: $\int_0^2 2x\,dx=\bigl[x^2\bigr]_0^2=4$ — совпадает.

Теперь посмотрим на потоки через грани. Для левого куба грань $x=1$ — внешняя (правая грань), с нормалью $\vec n=(1,0,0)$: $\vec F\cdot\vec n=x^2=1$, вклад $=1\cdot1=1$ (площадь грани $=1$). Для правого куба та же самая плоскость $x=1$ — тоже граница, но теперь это его левая грань, с внешней нормалью $\vec n=(-1,0,0)$: $\vec F\cdot\vec n=-x^2=-1$, вклад $=-1$.

Ответ: вклады внутренней грани $x=1$ для левого и правого кубов равны $+1$ и $-1$ — они в точности взаимно сокращаются, как и предсказывает идея доказательства. Оставшиеся вклады внешней границы: грань $x=0$ даёт $0$ (там $\vec F=(0,0,0)$), грань $x=2$ даёт $\vec F\cdot\vec n=x^2=4$, вклад $=4$; остальные четыре грани каждого куба дают $0$ (поле перпендикулярно их нормалям). Суммарный поток через внешнюю границу: $0+4=4$ — совпадает с тройным интегралом по всему объёму. Именно этот механизм — взаимное сокращение вкладов на внутренних гранях и накопление только внешних — лежит в основе того, почему сумма локальных изменений вдоль непрерывной траектории (будь то дивергенция векторного поля в непрерывном нормализующем потоке или ряд последовательных обновлений весов при обучении) сводится к чистому суммарному эффекту на границе всего процесса, а не к запутанной сумме промежуточных вкладов.

Почему это важно

Понимание того, откуда берётся теорема, а не просто умение подставлять числа в готовую формулу, — это именно то, что отличает уверенное владение математическим аппаратом от механического заучивания. Идея взаимного сокращения внутренних потоков — не изолированный трюк для доказательства одной конкретной теоремы: это универсальный принцип телескопирования, который встречается всюду, где сумма локальных вкладов сводится к разности на границах, — от формулы Ньютона-Лейбница до вывода уравнений сохранения в физике и вероятностных моделях.

Вычисление объёма тела через формулу как частный случай

Интуиция: специально подобранное поле превращает поток в объём

Давай разберёмся с ещё одним неожиданным применением теоремы. Что, если специально подобрать векторное поле так, чтобы его дивергенция была равна ровно единице во всех точках пространства? Тогда правая часть формулы Остроградского-Гаусса, тройной интеграл от единичной дивергенции по объёму $V$, автоматически даёт ровно объём этой области: $\iiint_V 1\,dV=V$. А значит, левая часть — поток такого специального поля через границу — тоже численно равна объёму тела!

Это удивительный, на первый взгляд контринтуитивный факт: объём трёхмерного тела можно вычислить, вообще не интегрируя по самому объёму, а лишь пробежавшись по его поверхности. Практическая польза колоссальна: если граница тела задаётся набором явных, легко параметризуемых поверхностей (плоскостей, кусков сфер, цилиндров), а сам объём тела описывается сложной системой неравенств, вычисление потока через границу часто оказывается значительно проще прямого тройного интегрирования.

Теорема

Следствие (вычисление объёма через поток). Если взять поле $\vec F=\left(\dfrac{x}3,\dfrac{y}3,\dfrac{z}3\right)$, его дивергенция равна $\operatorname{div}\vec F=\dfrac13+\dfrac13+\dfrac13=1$ в каждой точке пространства. Тогда объём ограниченной области $V$ равен потоку этого поля через её границу:

$$V=\iiint_V dV=\iint_{\partial V}\left(\frac{x}3,\frac{y}3,\frac{z}3\right)\cdot\vec n\,dS=\frac13\iint_{\partial V}\bigl(x\,dy\,dz+y\,dz\,dx+z\,dx\,dy\bigr)$$

Существуют и другие подходящие поля с единичной дивергенцией (например, $\vec F=(x,0,0)$ или $\vec F=(0,0,z)$), каждое из которых даёт свою, иногда более простую формулу для конкретной геометрии тела — выбор поля обычно определяется тем, какая из граней или частей границы устроена проще для интегрирования.

Разбор примеров

Пример 1 (лёгкий). Найти объём шара радиуса $R=3$, используя формулу потока поля $\vec F=\left(\dfrac{x}3,\dfrac{y}3,\dfrac{z}3\right)$ через сферу.

На сфере радиуса $R$ внешняя единичная нормаль равна $\vec n=\left(\dfrac{x}R,\dfrac{y}R,\dfrac{z}R\right)$. Скалярное произведение:

Это постоянная величина на всей сфере (поскольку на сфере радиуса $R$ выполняется $x^2+y^2+z^2=R^2$). Поток равен этому постоянному значению, умноженному на площадь сферы $4\pi R^2$:

При $R=3$: $V=\dfrac{4\pi\cdot27}3=36\pi$.

Ответ: $V=36\pi\approx113{,}097$ — та самая классическая формула объёма шара, выведенная без единого тройного интеграла.

Пример 2 (средний). Найти объём цилиндра радиуса $R=2$ и высоты $h=5$ через ту же формулу потока, разбивая границу на боковую поверхность и два основания.

На боковой поверхности $x=R\cos\theta$, $y=R\sin\theta$, $z=z$ внешняя нормаль $\vec n=(\cos\theta,\sin\theta,0)$. Тогда:

Площадь боковой поверхности $=2\pi Rh$, вклад $=\dfrac{R}3\cdot2\pi Rh=\dfrac{2\pi R^2h}3$. При $R=2$, $h=5$: $\dfrac{2\pi\cdot4\cdot5}3=\dfrac{40\pi}3$.

На верхнем основании $z=h$, нормаль $(0,0,1)$: $\vec F\cdot\vec n=\dfrac{z}3=\dfrac{h}3$, площадь $=\pi R^2$, вклад $=\dfrac{h}3\pi R^2=\dfrac{\pi R^2h}3$. При тех же числах: $\dfrac{\pi\cdot4\cdot5}3=\dfrac{20\pi}3$. На нижнем основании $z=0$: вклад равен $0$.

Ответ: суммарный поток $=\dfrac{40\pi}3+\dfrac{20\pi}3+0=\dfrac{60\pi}3=20\pi$, что точно совпадает с прямой формулой объёма цилиндра $V=\pi R^2h=\pi\cdot4\cdot5=20\pi$.

Пример 3 (сложный, машинное обучение). Стандартный трёхмерный симплекс $\{w_1\ge0,\ w_2\ge0,\ w_3\ge0,\ w_1+w_2+w_3\le a\}$ (при $a=2$) — геометрический образ допустимой области параметров, суммарно не превосходящих некоторого предела (например, область до нормировки в распределении, где $w_i$ — ненормированные веса категорий дискретного распределения). Найти его объём через формулу потока.

Граница симплекса состоит из четырёх плоских граней: трёх координатных ($w_1=0$, $w_2=0$, $w_3=0$) и одной наклонной ($w_1+w_2+w_3=a$). На грани $w_1=0$ внешняя нормаль $\vec n=(-1,0,0)$ (область лежит в полупространстве $w_1\ge0$), и первая компонента поля там равна $\dfrac{w_1}3=0$, так что $\vec F\cdot\vec n=0$. По той же логике грани $w_2=0$ и $w_3=0$ тоже дают нулевой вклад.

На наклонной грани $w_1+w_2+w_3=a$ внешняя единичная нормаль равна $\vec n=\left(\dfrac1{\sqrt3},\dfrac1{\sqrt3},\dfrac1{\sqrt3}\right)$. Тогда:

— постоянная величина на всей этой грани. Сама грань представляет собой равносторонний треугольник со стороной $a\sqrt2$ (расстояние между вершинами $(a,0,0)$ и $(0,a,0)$), а его площадь равна $\dfrac{\sqrt3}4\bigl(a\sqrt2\bigr)^2=\dfrac{\sqrt3}4\cdot2a^2=\dfrac{a^2\sqrt3}2$. Вклад этой грани:

При $a=2$: $V=\dfrac86=\dfrac43$.

Ответ: $V=\dfrac43\approx1{,}333$, что совпадает с известной формулой объёма координатного тетраэдра $V=\dfrac{a^3}6$. Содержательный вывод для машинного обучения: объём допустимой области параметров (в нашем примере — область, где сумма ненормированных весов не превышает заданного предела) можно вычислить, интегрируя не по самому объёму, а по его границе, — а это может быть значительно удобнее, если границы области заданы простыми линейными или иными явными ограничениями, что типично для областей допустимых значений в задачах с ограничениями.

Почему это важно

Возможность посчитать объём через поток — это не изолированный математический курьёз, а рабочий инструмент, который регулярно используется там, где прямое интегрирование по объёму затруднено, а граница тела устроена значительно проще. Более того, сам факт, что объём — это частный случай потока некоторого специального поля, лишний раз подчёркивает, насколько универсальна формула Остроградского-Гаусса: она не просто связывает поток и дивергенцию, а даёт единый язык, на котором можно выражать самые разные геометрические и физические величины.

Упрощение потоков через сложные поверхности

Интуиция: заткни дырявую пещеру простой заглушкой

Представь, что ты — спелеолог, изучающий сложную систему пещер, и хочешь измерить, сколько воды суммарно протекает через все входы и выходы этой запутанной пещерной системы. Вместо того чтобы аккуратно измерять поток через каждый извилистый проход по отдельности, ты можешь просто «заткнуть» всю систему одной простой воображаемой заглушкой (скажем, плоским диском) — и если внутри пещеры нет источников или стоков воды (вода просто протекает насквозь, не появляясь и не исчезая), суммарный поток через заглушку будет в точности равен суммарному потоку через настоящий сложный вход.

Формально эта идея опирается на теорему Остроградского-Гаусса напрямую: если дивергенция поля равна нулю всюду внутри некоторого объёма, ограниченного двумя разными поверхностями с общим краем, то потоки через обе эти поверхности (при согласованной ориентации) обязаны совпадать — ведь их разность равна тройному интегралу от нулевой дивергенции, то есть нулю.

Теорема

Следствие (замена поверхности). Пусть поверхности $S_1$ и $S_2$ имеют общий край (границу) и вместе с этим краем ограничивают объём $V$, в котором $\operatorname{div}\vec F=0$ всюду. Тогда при согласованной ориентации (обе нормали направлены «в одну и ту же сторону» относительно объёма $V$) потоки через $S_1$ и $S_2$ равны:

$$\iint_{S_1}\vec F\cdot\vec n\,dS=\iint_{S_2}\vec F\cdot\vec n\,dS$$

В частности, если поле имеет особую точку (например, точечный заряд), но $\operatorname{div}\vec F=0$ везде, кроме этой единственной точки, то поток через любую замкнутую поверхность, охватывающую особую точку, один и тот же — независимо от формы и размера поверхности, лишь бы она охватывала ровно ту же особенность.

Разбор примеров

Пример 1 (лёгкий). Вычислить поток постоянного поля $\vec F=(0,0,1)$ через верхнюю полусферу радиуса $R=3$ (часть сферы $x^2+y^2+z^2=R^2$ при $z\ge0$), заменив её на плоский диск с той же границей.

Дивергенция постоянного поля равна нулю всюду ($\operatorname{div}\vec F=0$), значит поток через полусферу (с нормалью, направленной вверх/наружу) равен потоку через диск $x^2+y^2\le R^2$ в плоскости $z=0$ с той же ориентацией (нормаль вверх, $\vec n=(0,0,1)$), поскольку вместе они образуют замкнутую поверхность (полусфера сверху и диск снизу), ограничивающую полушарие без источников поля внутри.

На диске: $\vec F\cdot\vec n=1$ повсюду, площадь диска $=\pi R^2=\pi\cdot9=9\pi$.

Ответ: поток через диск, а значит и через полусферу, равен $9\pi\approx28{,}274$ — при этом прямое вычисление по полусфере потребовало бы явной параметризации через сферические углы, а замена на плоский диск свела задачу к элементарному умножению.

Пример 2 (средний, закон Гаусса). Точечный заряд величины $q$ находится в начале координат, создавая поле $\vec F=q\dfrac{\vec r}{|\vec r|^3}$ (здесь $\vec r=(x,y,z)$). Вычислить поток этого поля через сферу радиуса $R$ и объяснить, почему тот же ответ получится для любой другой замкнутой поверхности, охватывающей заряд.

На сфере радиуса $R$ поле имеет модуль $|\vec F|=q/R^2$ и направлено строго вдоль внешней нормали (радиально), значит $\vec F\cdot\vec n=q/R^2$ — постоянная величина на всей сфере. Площадь сферы $=4\pi R^2$, поток:

Результат вообще не зависит от $R$ — что не случайность: можно проверить прямым дифференцированием, что $\operatorname{div}\vec F=0$ всюду, кроме самого начала координат (где поле не определено). Значит, для любой другой замкнутой поверхности $S$, охватывающей заряд, объём между $S$ и маленькой сферой вокруг заряда не содержит источников поля, и по следствию о замене поверхности поток через $S$ обязан совпадать с потоком через сферу.

Ответ: поток равен $4\pi q$ независимо от радиуса $R$ и вообще от формы охватывающей поверхности — это в точности закон Гаусса в электростатике.

Пример 3 (сложный, машинное обучение). В некоторых архитектурах нормализующих потоков (например, объём-сохраняющих преобразованиях) требуется, чтобы преобразование не искажало плотность вероятности — формально это условие $\operatorname{div}\vec f=0$ для поля скорости $\vec f$. Рассмотрим поле $\vec f=(y,-x,0)$ (чисто «вихревое» вращение вокруг оси $z$) и покажем, что оно объём-сохраняющее, вычислив его поток через две совершенно разные замкнутые поверхности: сферу радиуса $2$ и куб $[-1,1]^3$.

Проверим дивергенцию: $\operatorname{div}\vec f=\dfrac{\partial y}{\partial x}+\dfrac{\partial(-x)}{\partial y}+\dfrac{\partial0}{\partial z}=0+0+0=0$. По формуле Остроградского-Гаусса поток через сферу радиуса $2$ сразу равен $\iiint 0\,dV=0$.

Проверим это независимо и для куба, вычислив вклад одной грани явно, чтобы убедиться в согласованности. На грани $x=1$: нормаль $\vec n=(1,0,0)$, $\vec f\cdot\vec n=y$. Интегрируем по квадрату $y\in[-1,1]$, $z\in[-1,1]$:

(интеграл нечётной функции $y$ по симметричному отрезку равен нулю). По аналогичной причине (нечётность подынтегрального выражения относительно соответствующей координаты) все остальные пять граней куба тоже дают нулевой вклад.

Ответ: поток через обе поверхности — и сферу, и куб — равен $0$, что подтверждает теорему и одновременно демонстрирует главное практическое следствие для генеративных моделей: если поле скорости потока имеет нулевую дивергенцию (объём-сохраняющее преобразование), суммарный «отток» вероятностной массы через любую воображаемую замкнутую границу в пространстве состояний всегда равен нулю — облако вероятности меняет форму, но не «теряет» и не «приобретает» суммарную массу ни через одну поверхность, какую бы причудливую форму она ни имела.

Почему это важно

Приём замены сложной поверхности на удобную — это, пожалуй, самое практически ценное следствие теоремы Остроградского-Гаусса. Закон Гаусса в электростатике, расчёт объём-сохраняющих преобразований в генеративных моделях, оценка потоков в сложных инженерных конструкциях — во всех этих задачах прямое вычисление поверхностного интеграла по «настоящей» сложной границе было бы практически невозможным, а замена на удобную поверхность (при условии отсутствия источников между ними) превращает задачу в элементарную.

Практика: 30 заданий

Базовые задания (1–10)

Задание 1: Найти дивергенцию поля $\vec F=(x^2,y^2,z^2)$ в точке $(1,2,3)$.

Задание 2: Вычислить поток поля $\vec F=(x^2,y^2,z^2)$ через сферу радиуса $3$ с центром в начале координат.

Задание 3: Для поля $\vec F=(2x,3y,z)$ и куба $[0,2]^3$ вычислить $\operatorname{div}\vec F$ и поток через полную поверхность куба.

Задание 4: Вычислить поток поля $\vec F=(yz,xz,xy)$ через сферу радиуса $1$.

Задание 5: Проверить теорему на параллелепипеде $[0,1]\times[0,2]\times[0,3]$ с полем $\vec F=(x,0,0)$, вычислив поток напрямую.

Задание 6: Найти дивергенцию поля $\vec F=(e^x\sin y,\,e^x\cos y,\,z^2)$ и её значение в точке $(0,0,1)$.

Задание 7: Найти объём шара радиуса $R=2$, используя формулу потока с полем $\vec F=\left(\dfrac{x}3,\dfrac{y}3,\dfrac{z}3\right)$.

Задание 8: Вычислить поток поля $\vec F=(x,y,z)$ через полную поверхность цилиндра радиуса $1$ и высоты $5$.

Задание 9: Вычислить поток поля $\vec F=(x^3,y^3,z^3)$ через сферу радиуса $1$.

Задание 10: Найти объём тетраэдра $w_1\ge0$, $w_2\ge0$, $w_3\ge0$, $w_1+w_2+w_3\le2$ через формулу потока.

Средние задания (11–20)

Задание 11: Вычислить поток поля $\vec F=(x^2,y^2,z^2)$ через полную поверхность единичного куба $[0,1]^3$, используя тройной интеграл.

Задание 12: Вычислить поток поля $\vec F=(0,0,z^2)$ через полную поверхность цилиндра радиуса $2$ и высоты $3$ двумя способами.

Задание 13: Найти объём конуса с вершиной в начале координат, высотой $h=3$ вдоль оси $z$ и радиусом основания $R=2$, вычислив тройной интеграл в цилиндрических координатах.

Задание 14: Точечный заряд $q=5$ помещён внутри произвольной замкнутой поверхности. Найти поток поля $\vec F=q\vec r/|\vec r|^3$ через сферу радиуса $10$ вокруг заряда.

Задание 15: Заряд $q=7$ находится внутри произвольной (не обязательно сферической) замкнутой поверхности $S$. Найти поток поля через $S$.

Задание 16: Вычислить поток поля $\vec F=(0,0,1)$ через верхнюю полусферу радиуса $R=3$, заменив её на диск.

Задание 17: Найти объём симплекса $w_1,w_2,w_3\ge0$, $w_1+w_2+w_3\le3$ через формулу потока.

Задание 18 (машинное обучение): В непрерывном нормализующем потоке поле скорости $\vec f=(x,y,z)$, начальная (ненормированная) плотность равна $1$ на единичном шаре. Найти суммарный поток "вероятностной массы" через границу шара.

Задание 19: Вычислить поток поля $\vec F=(y,-x,0)$ через полную поверхность куба $[-1,1]^3$.

Задание 20: Вычислить поток того же поля $\vec F=(y,-x,0)$ через сферу радиуса $5$.

Продвинутые задания (21–30)

Задание 21: Заряд $q=2$ смещён в точку $(0,0,0)$ и находится внутри куба $[-2,2]^3$. Найти поток поля $\vec F=q\vec r/|\vec r|^3$ через поверхность куба.

Задание 22: Найти объём эллипсоида $\dfrac{x^2}{1^2}+\dfrac{y^2}{2^2}+\dfrac{z^2}{3^2}\le1$ через обобщение формулы потока (с заменой $x=au$, $y=bv$, $z=cw$, сводящей эллипсоид к единичному шару с якобианом $abc$).

Задание 23 (машинное обучение): Поле скорости объём-сохраняющего потока $\vec f=(\sin y+z,\ \cos x-z,\ x-y)$. Проверить, что $\operatorname{div}\vec f=0$, и найти поток через сферу радиуса $2$.

Задание 24: Вычислить поток поля $\vec F=(x^3,y^3,z^3)$ через сферу радиуса $2$.

Задание 25: Найти объём шарового слоя между сферами радиуса $1$ и радиуса $2$ (с общим центром) через формулу потока поля $\vec F=\left(\dfrac{x}3,\dfrac{y}3,\dfrac{z}3\right)$.

Задание 26: Найти дивергенцию поля $\vec F=(x^2y,\,-xy^2,\,z)$ и вычислить поток через полную поверхность цилиндра радиуса $1$ и высоты $2$.

Задание 27 (машинное обучение): Равномерная плотность $\rho=1$ задана в шаре радиуса $R=1$, поле скорости $\vec v=\dfrac{(x,y,z)}\tau$ с характерным временем $\tau=2$. Найти поток "вероятностной массы" $\rho\vec v$ через границу шара и сравнить с тройным интегралом от $\operatorname{div}(\rho\vec v)$.

Задание 28: Найти объём конуса с радиусом основания $R=3$ и высотой $h=4$ через формулу тройного интеграла в цилиндрических координатах.

Задание 29: Заряд $q=3$ смещён в точку $(1,1,1)$ внутри куба $[-3,3]^3$. Найти поток поля $\vec F=q(\vec r-\vec r_0)/|\vec r-\vec r_0|^3$ через поверхность куба, где $\vec r_0=(1,1,1)$.

Задание 30 (машинное обучение): В непрерывном нормализующем потоке поле скорости $\vec f=(x,y,z)$. Найти скорость изменения логарифма плотности вероятности $\dfrac{d}{dt}\log p$ в произвольной точке и суммарный поток $\vec f$ через границу единичного шара.

Частые ошибки

Давай пройдёмся по типичным промахам, которые встречаются у тех, кто впервые применяет формулу Остроградского-Гаусса.

Первая и самая частая ошибка — перепутать внешнюю и внутреннюю нормаль. Теорема справедлива именно для нормали, направленной наружу из объёма $V$; если по невнимательности взять нормаль внутрь, весь поток поменяет знак на противоположный, а результат окажется неверным без каких-либо признаков ошибки в самих вычислениях.

Вторая ошибка — применить формулу к незамкнутой поверхности. Формула Остроградского-Гаусса работает только тогда, когда поверхность $\partial V$ реально ограничивает некоторый объём целиком, без «дыр» и открытых краёв — для одиночного куска поверхности (например, одной полусферы без диска-основания) нужно либо мысленно замкнуть её дополнительной поверхностью, либо считать поверхностный интеграл напрямую, без применения теоремы.

Третья ошибка — забыть проверить, что поле непрерывно дифференцируемо во всей замкнутой области, включая границу. Если поле имеет особую точку внутри объёма (как электростатическое поле точечного заряда в начале координат), теорему в её базовой форме напрямую применять нельзя — нужно либо исключить малую окрестность особой точки из области интегрирования, либо использовать следствие о равенстве потоков через поверхности с общим краем.

Четвёртая ошибка — при вычислении дивергенции перепутать, какая частная производная берётся по какой переменной: $\operatorname{div}(P,Q,R)=\partial P/\partial x+\partial Q/\partial y+\partial R/\partial z$, а не, скажем, $\partial P/\partial y+\partial Q/\partial x+\partial R/\partial z$ (это была бы уже совсем другая величина, не имеющая отношения к теореме).

Пятая ошибка — при использовании формулы вычисления объёма через поток забыть, что поле должно иметь именно единичную дивергенцию, а не произвольную. Если взять, скажем, поле $\vec F=(x,y,z)$ (с дивергенцией $3$, а не $1$), поток через границу даст не объём тела, а его утроенное значение — легко забыть про этот масштабирующий множитель.

Шестая ошибка, особенно частая в задачах на замену сложной поверхности удобной: применить приём замены, не проверив, что дивергенция поля действительно равна нулю во всём объёме между двумя поверхностями. Если внутри этого объёма есть хотя бы точка с ненулевой дивергенцией, потоки через две поверхности, вообще говоря, различаются — и разница как раз равна тройному интегралу от дивергенции по объёму между ними.

Главное запомнить

Прежде чем двигаться дальше, зафиксируем ключевые идеи этого урока в сжатом виде.

• Формула Остроградского-Гаусса связывает поток векторного поля через замкнутую поверхность $\partial V$ с тройным интегралом от дивергенции поля по объёму $V$: $\iint_{\partial V}\vec F\cdot\vec n\,dS=\iiint_V\operatorname{div}\vec F\,dV$.

• Дивергенция поля $\vec F=(P,Q,R)$ вычисляется как $\operatorname{div}\vec F=\partial P/\partial x+\partial Q/\partial y+\partial R/\partial z$ и физически означает локальную «плотность источников» поля в данной точке.

• Условия применимости: область $V$ ограничена, граница $\partial V$ кусочно-гладкая и замкнутая, поле непрерывно дифференцируемо во всей замкнутой области; нормаль всегда берётся внешней.

• Идея доказательства строится на разбиении объёма на малые кубики: потоки через внутренние общие грани соседних кубиков взаимно сокращаются, остаётся только вклад внешней границы всего объёма.

• Формула Остроградского-Гаусса — это прямое трёхмерное обобщение формулы Грина: там циркуляция по контуру связывалась с ротором внутри плоской области, здесь поток через поверхность связывается с дивергенцией внутри объёма.

• Специально подобранное поле с единичной дивергенцией (например, $\vec F=(x/3,y/3,z/3)$) превращает поток через границу тела в его объём — неожиданный, но очень практичный частный случай формулы.

• Если дивергенция поля равна нулю в объёме между двумя поверхностями с общим краем, потоки через эти поверхности совпадают — это позволяет заменять сложные для интегрирования поверхности удобными (закон Гаусса в электростатике — прямое следствие этого факта).

• В непрерывных нормализующих потоках скорость изменения логарифма плотности вероятности вдоль траектории равна минус дивергенции векторного поля скорости — прямое, дифференциальное применение той же идеи, что и в теореме.

• Объём-сохраняющие преобразования в генеративных моделях характеризуются нулевой дивергенцией поля скорости — а значит, суммарный поток вероятностной массы через любую замкнутую границу в пространстве состояний для таких преобразований всегда равен нулю.

• Теорема работает в обе стороны: сложный поверхностный интеграл иногда выгоднее заменить простым тройным интегралом, а иногда, наоборот, сложный тройной интеграл — простым поверхностным.

Связь с другими темами курса

Формула Остроградского-Гаусса опирается на целую цепочку понятий, пройденных ранее, и продолжает главную линию курса — построение мостов между локальными и глобальными характеристиками поля. Само понятие дивергенции ты уже встречал в уроке про теорию поля, где разбирался с дивергенцией и ротором как точечными характеристиками векторного поля; сегодняшняя теорема как раз и показывает, зачем эти точечные характеристики вообще нужны — они интегрируются в глобальные, измеримые на практике величины. Вычисление тройных интегралов, лежащее в основе правой части формулы, напрямую использует технику, отработанную ещё в уроке про тройные интегралы, включая переход к сферическим и цилиндрическим координатам.

Особенно тесная связь — с формулой Грина из предыдущего урока. Формула Остроградского-Гаусса — это её прямое трёхмерное обобщение: там, где формула Грина связывала циркуляцию поля по замкнутому плоскому контуру с двойным интегралом от ротора внутри области, здесь формула Остроградского-Гаусса связывает поток поля через замкнутую поверхность с тройным интегралом от дивергенции внутри объёма. Обе теоремы — частные случаи одной и той же более общей идеи «интеграл по границе равен интегралу характеристики поля по внутренности», и обе, в свою очередь, оказываются частными случаями ещё более общей теоремы Стокса, к которой курс перейдёт в следующем уроке — там мы наконец увидим, как формула Грина, формула Остроградского-Гаусса и классическая теорема Стокса связаны единой формулировкой на языке дифференциальных форм.

Интересные факты

Несколько любопытных деталей, которые стоит знать помимо основной теории.

• Уравнения Максвелла в их интегральной форме — это, по сути, четыре применения формулы Остроградского-Гаусса и теоремы Стокса к электрическому и магнитному полям; закон Гаусса для электрического поля, который мы разобрали в примере про точечный заряд, — лишь первое из этих четырёх уравнений.

• Формула для вычисления объёма тела через поверхностный интеграл — это математическая основа алгоритма, основанного на теореме о дивергенции, который до сих пор используется в компьютерной графике и системах автоматизированного проектирования для быстрого вычисления объёма трёхмерных полигональных моделей: вместо разбиения модели на тетраэдры программа просто суммирует вклады всех треугольных граней поверхности.

• Уравнение неразрывности в гидродинамике и его вероятностный аналог — уравнение Фоккера-Планка, описывающее эволюцию плотности вероятности случайного процесса во времени, — оба выводятся из формулы Остроградского-Гаусса, применённой к бесконечно малому объёму: это буквально дифференциальная (точечная) форма той же самой теоремы, которую мы изучали сегодня в её интегральной форме.

• Термин «дивергенция» (от латинского divergere — «расходиться») был предложен именно в контексте векторного анализа XIX века и с тех пор используется совершенно неизменно во всех областях, где применяется идея «локального источника поля», — от классической электродинамики до современных архитектур нейросетей, вычисляющих дивергенцию векторных полей автоматическим дифференцированием.

Лайфхаки и полезные трюки

Несколько практических приёмов, которые экономят время и снижают риск ошибок при работе с формулой Остроградского-Гаусса.

• Прежде чем вычислять сложный поверхностный интеграл по замкнутой поверхности напрямую, всегда сначала проверь, не проще ли посчитать дивергенцию поля и свести задачу к тройному интегралу — особенно если дивергенция получается заметно проще исходных компонент поля (как в примере с $e^x\sin y$ и $e^x\cos y$, где слагаемые взаимно уничтожаются).

• Если тело имеет высокую симметрию (шар, куб, симметричные относительно начала координат области) и поле содержит нечётные степени координат, сразу проверь, не обнулится ли тройной интеграл по симметрии — это часто избавляет от долгих вычислений в лоб.

• Для вычисления объёма тела со сложной, но кусочно-плоской или кусочно-параметризуемой границей, попробуй применить формулу потока со специально подобранным полем единичной дивергенции — часто это оказывается быстрее прямого тройного интегрирования, особенно если некоторые грани дают нулевой вклад из-за перпендикулярности поля и нормали.

• Если нужно вычислить поток поля через поверхность с особой точкой внутри (например, электростатическое или гравитационное поле точечного источника), сразу вспомни о приёме замены поверхности: замени сложную поверхность на маленькую сферу вокруг особой точки — результат гарантированно совпадёт, если дивергенция вне особой точки равна нулю.

• Для проверки правильности решения всегда, когда это возможно, вычисляй обе части формулы независимо — поверхностный интеграл напрямую и тройной интеграл от дивергенции; совпадение результатов служит надёжной защитой от арифметических ошибок в длинных вычислениях.

• При работе с полями, связанными с генеративными моделями или нормализующими потоками, сразу проверь дивергенцию поля скорости: нулевая дивергенция сигнализирует об объём-сохраняющем преобразовании, а её конкретное значение прямо говорит о скорости локального «сжатия» или «расширения» плотности вероятности вдоль траектории.

Формула Остроградского-Гаусса закрывает важный шаг в понимании того, как локальные, точечные характеристики поля — дивергенция в каждой отдельной точке пространства — складываются в глобальные, измеримые на практике величины: суммарный поток через границу целого объёма. Ты увидел, что эта идея работает не только как абстрактная теорема, но и как рабочий вычислительный инструмент — позволяющий находить объёмы тел через интегрирование по их поверхности, упрощать расчёт потоков через сложные конфигурации заменой на удобные поверхности и выводить закон Гаусса в электростатике почти без усилий. А главное — ты увидел, что та же самая логика лежит в основе того, как современные генеративные модели отслеживают изменение плотности вероятности вдоль непрерывных траекторий: дивергенция векторного поля скорости в этих моделях — не абстрактная математическая деталь, а буквально та же самая величина, которую ты сегодня научился вычислять и интерпретировать. В следующем уроке мы познакомимся с теоремой Стокса — общей формулировкой, которая объединит сегодняшнюю формулу Остроградского-Гаусса и формулу Грина из прошлого урока в единую конструкцию.