Теорема Гаусса - Остроградского

Цели урока

• Понять физический смысл дивергенции как плотности источников векторного поля
• Формулировать и применять теорему Гаусса-Остроградского для вычисления потоков
• Знать закон Гаусса в электродинамике и его связь с теоремой дивергенции
• Видеть теорему как частный случай обобщённой теоремы Стокса в размерности 3

Предварительные знания

• Теорема Грина
• Поверхностные интегралы
• Дивергенция и ротор

• Теорема Грина

Как вычислить полный выходящий поток жидкости через поверхность ракурса, не интегрируя по каждому квадратному сантиметру поверхности?

• NVIDIA PhysX в 4000 игровых движков: теорема дивергенции конвертирует дорогие поверхностные интегралы в объёмные, ускоряя физику жидкостей в 10-100x
• Закон Гаусса в электродинамике: поток E через любую поверхность = Q_enclosed / epsilon_0 - прямое следствие div E = rho/epsilon_0
• OpenFOAM и CFD: метод конечных объёмов для расчёта аэродинамики Airbus и NASA строится на дискретной версии теоремы Гаусса
• Астрофизика: масса внутри орбиты определяет гравитационную силу независимо от распределения - следствие теоремы Гаусса для гравитации

Тройное независимое открытие

Михаил Остроградский (1801-1862) доказал теорему в 1826 году, представил Парижской академии наук в 1828-м. Карл Фридрих Гаусс опубликовал версию в 1813 году (для частного случая). Джордж Грин доказал 2D-случай в 1828-м. Три человека, три разные постановки, одна идея. Классический пример математической конвергенции: фундаментальные истины переоткрываются независимо. Сегодня теорему называют по-разному в разных странах: теорема Гаусса (в физике), теорема Остроградского (в России), теорема дивергенции (в математике). Под всеми именами - один результат.

Теорема Гаусса - Остроградского

1831 год. Остроградский доказывает теорему в Петербурге. Гаусс публикует то же самое в Германии. Грин доказал двумерный случай ещё в 1828-м. Три человека, три страны, три года разницы - одна теорема. NVIDIA PhysX использует её в 4000 игровых движков для симуляции жидкости: вместо интегрирования по всей поверхности - дивергенция по объёму, экономия 10-100x.

Для поля E в электростатике: div E = rho/epsilon_0 (уравнение Максвелла). Теорема Гаусса немедленно даёт: поток E через любую поверхность = полный заряд внутри / epsilon_0. Именно так работает закон Гаусса в физике.

Физический смысл дивергенции и уравнение неразрывности

Уравнение неразрывности в гидродинамике: d(rho)/dt + div(rho*v) = 0. Для несжимаемой жидкости (rho = const): div v = 0. Это значит: поток через любую замкнутую поверхность внутри потока равен нулю - нет ни источников, ни стоков. Авиационные CFD-симуляторы Boeing и Airbus решают это уравнение на сетках из 10-100 миллионов ячеек.

Теорема Гаусса требует, чтобы F имело непрерывные частные производные в замкнутом объёме V. Если поле имеет особенность внутри (как гравитация точечной массы в r=0) - нужно исключить особую точку малым шаром и взять предел.

Теорема дивергенции в вычислительной физике

FVM (Finite Volume Method) - метод конечных объёмов - строится целиком на теореме Гаусса. Расчётная область разбивается на маленькие ячейки. В каждой ячейке интеграл дивергенции заменяется суммой потоков через грани. OpenFOAM - открытый CFD-солвер - решает уравнения Навье-Стокса на сетках из миллионов ячеек. Airbus использует его для расчёта аэродинамики. NASA - для симуляции сопел ракетных двигателей.

Идея FVM: вместо дифференциального уравнения div F = q решаем его интегральный аналог. Для каждой ячейки: сумма потоков через грани = источники внутри. Это именно теорема Гаусса в дискретной форме.

Связь с теоремой Стокса и обобщение

Теоремы Грина, Стокса и Гаусса - три лика одной теоремы. Это станет ясно с дифференциальными формами: int_M d(omega) = int_{partial M} omega. При dim M = 1: формула Ньютона-Лейбница. При dim M = 2: теорема Грина. При dim M = 2 в R^3: теорема Стокса. При dim M = 3: теорема Гаусса.

Четыре уравнения Максвелла в компактной записи через дифференциальные формы: dF = 0 и d(*F) = J. Здесь F - 2-форма электромагнитного поля на 4-мерном пространстве-времени. Теоремы Гаусса и Стокса спрятаны внутри оператора d.

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

Теорема Гаусса пронизывает всю математическую физику - от уравнений Максвелла до уравнений Навье-Стокса.

• Электродинамика — Связанная тема
• Гидродинамика — Связанная тема
• Гравитация — Связанная тема
• Численные методы — Связанная тема

Итоги

• Теорема Гаусса: oiint_{dV} F dS = iiint_V div F dV - поток через поверхность = интеграл дивергенции по объёму
• Дивергенция в точке = предел потока через малую сферу / объём сферы - координатно-независимое определение
• div B = 0 для магнитного поля: нет монополей, поток через любую замкнутую поверхность = 0
• Тождество div(rot F) = 0 следует из d^2 = 0 - фундаментальное алгебраическое свойство
• FVM: дискретная теорема Гаусса на сетке - основа CFD-симуляций в авиации и энергетике

Вопросы для размышления

• Почему теорема Гаусса позволяет вычислить гравитацию от сферически-симметричного тела как от точечной массы - что здесь физически происходит?
• Как уравнение неразрывности div(rho*v) = 0 для несжимаемой жидкости связано с теоремой Гаусса - откуда берётся это уравнение?
• В чём принципиальное сходство между теоремой Грина (2D) и теоремой Гаусса (3D) - есть ли единый принцип, из которого они оба следуют?

Связанные уроки

• calc-25-green-theorem — Теорема Грина - 2D-предшественник
• calc-27-diff-forms — Формы - современный язык теоремы дивергенции
• calc-28-manifold-int — Интегрирование на многообразиях обобщает этот результат
• la-04-matrix-ops