Інтеграли Гауса

Ми можемо проінтегрувати кривину Ґауса по всьому об'єму $n$ -вимірного многовида:

(1) I_{m} = \int K^{[m]} \sqrt{g} d u^{1} d u^{2} \dots d u^{n}

Назвемо одержану величину інтегралом Ґауса.

Варіація інтеграла Ґауса

Нехай ми маємо замкнутий обмежений многовид, подібний до сфери чи тору. Варіація інтеграла Ґауса дається такою формулою (доведення буде нижче):

(2) δ I_{m} = \frac{1}{2} \int G_{i j}^{[m]} δ g^{i j} \sqrt{g} d u^{1} d u^{2} \dots d u^{n}

де тензор $G_{i j}^{[m]}$ у випадку гіперповерхні визначається наступною формулою (піднімемо один з індексів, щоб формула краще виглядала):

(3) G_{j}^{[m] i} = \frac{(- 1)^{m + 1}}{m!} g_{j p_{1} p_{2} \dots p_{m}}^{i s_{1} s_{2} \dots s_{m}} b_{s_{1}}^{p_{1}} b_{s_{2}}^{p_{2}} \dots b_{s_{m}}^{p_{m}}

В цій формулі використовується тензор метричної матрьошки $g_{j_{1} \dots j_{m}}^{i_{1} \dots i_{m}}$ (дивіться статтю Одиничний антисиметричний тензор), і тензор повної кривини гіперповерхні $b_{i j}$ .

Очевидно, формула (3) задає цілу серію тензорів другого рангу. Назвемо ці тензори тензорами Ейнштейна $m$ -го степеня, оскільки, як легко можна обчислити, тензор другого степеня збігається з класичним тензором Ейнштейна:

(4) G_{j}^{[2] i} = - \frac{1}{2} g_{j r s}^{i p q} b_{p}^{r} b_{q}^{s} = - \frac{1}{4} g_{j r s}^{i p q} R_{p q}^{r s} = R_{j}^{i} - \frac{1}{2} R δ_{j}^{i}

Замітимо, що при $m = n$ тензор метричної матрьошки в формулі (3) дорівнюватиме нулю, оскільки в кожній з двох антисиметричних груп індексів індекси будуть повторюватися. Тому тензор Ейнштейна буде нульовим, а отже і варіація інтеграла Ґауса (формула 2). Це означає, що в цьому випадку інтеграл Ґауса є топологічним інваріантом.

Обчислення варіації інтеграла Ґауса для випадку гіперповерхні

В наступних інтегралах цього пункту, для спрощення запису, ми не будемо писати добутків диференціалів координат $d u^{1} d u^{2} \dots d u^{n}$ . Маємо:

(5) δ I_{m} = δ (\int K^{[m]} \sqrt{g}) = \int (δ K^{[m]}) \sqrt{g} + \int K^{[m]} δ \sqrt{g}

Займемся першим доданком формули (5). Кривина Ґауса обчислюється за формулою:

(6) K^{[m]} = \frac{(- 1)^{m}}{m!} g_{p_{1} p_{2} \dots p_{m}}^{s_{1} s_{2} \dots s_{m}} b_{s_{1}}^{p_{1}} \dots b_{s_{m}}^{p_{m}}

Компоненти тензора метричної матрьошки дорівнюють або нулю, або одиниці, або мінус одиниці, а отже цей тензор можна винести за знак варіації як константу. Також, враховуючи симетричність тензора метричної матрьошки щодо перестановки «вертикальних пар» індексів ${s_{k} p_{k}}$ , одержуємо таку формулу для варіації кривини Ґауса:

(7) δ K^{[m]} = m \cdot \frac{(- 1)^{m}}{m!} g_{j p_{2} \dots p_{m}}^{i s_{2} \dots s_{m}} (δ b_{i}^{j}) b_{s_{2}}^{p_{2}} \dots b_{s_{m}}^{p_{m}} = G_{j}^{[m - 1] i} δ b_{i}^{j}

Знайдемо тепер варіацію тензора повної кривини гіперповерхні.

(8) δ b_{i}^{j} = δ (g^{j s} b_{s i}) = b_{s i} δ g^{j s} + g^{j s} δ b_{s i}

Підставляючи (8) в (7), знаходимо:

(9) δ K^{[m]} = R_{i j}^{[m]} δ g^{i j} + G^{[m - 1] i j} δ b_{i j}

де введено позначення тензора Річчі $m$ -го степеня:

(10) R_{j}^{[m] i} = \frac{(- 1)^{m}}{(m - 1)!} g_{j p_{2} \dots p_{m}}^{s_{1} s_{2} \dots s_{m}} b_{s_{1}}^{i} b_{s_{2}}^{p_{2}} \dots b_{s_{m}}^{p_{m}}

(аналогічно до тензора Ейнштейна, тензор Річчі другого степеня збігається з класичним тензором Річчі)

Далі, тензор повної кривини гіперповерхні дорівнює скалярному добутку вектора нормалі $𝐧$ до гіперповерхні і другої похідної радіус-вектора $𝐫_{i j}$ :

(11) b_{i j} = (𝐧 \cdot 𝐫_{i j})

Знаходимо варіацію:

(12) δ b_{i j} = (δ 𝐧 \cdot 𝐫_{i j}) + (𝐧 \cdot δ 𝐫_{i j})

Варіація одиничного вектора нормалі ортогональна до самого вектора нормалі:

(13) (δ 𝐧 \cdot 𝐧) = 0

а тому лежить в дотичній до многовида гіперплощині, і її можна розкласти за базисними векторами гіперповерхні:

(14) δ 𝐧 = α^{i} 𝐫_{i}

де $α^{i}$ — деякі, поки що невідомі коефіцієнти розкладу. Оскільки вектор нормалі завжди ортогональний до базисних векторів:

(15) (𝐧 \cdot 𝐫_{j}) = 0

то маємо також:

(16) δ (𝐧 \cdot 𝐫_{j}) = (δ 𝐧 \cdot 𝐫_{j}) + (𝐧 \cdot δ 𝐫_{j}) = α^{i} g_{i j} + (𝐧 \cdot δ 𝐫_{j}) = 0

Із (14) ш (16) маємо:

(17) δ 𝐧 = - g^{i j} (𝐧 \cdot δ 𝐫_{j}) 𝐫_{i}

Підставимо тепер (17) в (12). Знаходимо:

(18) δ b_{i j} = - g^{s p} (𝐧 \cdot δ 𝐫_{s}) (𝐫_{p} \cdot 𝐫_{i j}) + (𝐧 \cdot δ 𝐫_{i j}) = - Γ_{i j}^{s} (𝐧 \cdot δ 𝐫_{s}) + (𝐧 \cdot δ 𝐫_{i j})

Тепер ми готові обчислювати перший інтеграл формули (5). Підставляючи вирази (9) і (18), знаходимо:

(19) \int (δ K^{[m]}) \sqrt{g} = \int R_{i j}^{[m]} δ g^{i j} \sqrt{g} + \int G^{[m - 1] i j} ((𝐧 \cdot δ 𝐫_{i j}) - Γ_{i j}^{s} (𝐧 \cdot δ 𝐫_{s})) \sqrt{g}

Другий інтеграл в правій частині формули (19) розбивається на два, із них перший інтеграл ми можемо проінтегрувати частинами:

(20) \int (\sqrt{g} G^{[m - 1] i j} 𝐧 \cdot \frac{\partial}{\partial u^{j}} δ 𝐫_{i}) = \int \frac{\partial}{\partial u^{j}} (\sqrt{g} G^{[m - 1] i j} (𝐧 \cdot δ 𝐫_{i})) d u^{1} d u^{2} \dots d u^{n} - \int (\frac{\partial}{\partial u^{j}} (\sqrt{g} G^{[m - 1] i j} 𝐧) \cdot δ 𝐫_{i})

Перший інтеграл за формулою Остроградського-Ґауса перетворюється в інтеграл по межі області інтегрування:

(21) \oint G^{[m - 1] i j} (𝐧 \cdot δ 𝐫_{i}) d S_{j}

яка у випадку замкнутого многовида вироджується в точку (і відповідно цей інтеграл стає рівним нулю). Але ми можемо також розглядати і інтеграл Ґауса по частині многовиду, в цьому випадку інтеграл (21) треба залишити.

Другий інтеграл в правій частині формули (20) розбивається на два інтеграла:

(22) \int \frac{\partial}{\partial u^{j}} (\sqrt{g} G^{[m - 1] i j}) (𝐧 \cdot δ 𝐫_{i}) + \int \sqrt{g} G^{[m - 1] i j} (𝐧_{j} \cdot δ 𝐫_{i})

Збираючи формули (20-22) в (19), одержуємо:

(23) \int (δ K^{[m]}) \sqrt{g} = \int R_{i j}^{[m]} δ g^{i j} \sqrt{g} + \oint G^{[m - 1] i j} (𝐧 \cdot δ 𝐫_{i}) d S_{j} -

- \int (\frac{1}{\sqrt{g}} \frac{\partial}{\partial u^{j}} (\sqrt{g} G^{[m - 1] i j}) + Γ_{s j}^{i} G^{[m - 1] s j}) (𝐧 \cdot δ 𝐫_{i}) \sqrt{g} - \int \sqrt{g} G^{[m - 1] i j} (𝐧_{j} \cdot δ 𝐫_{i})

Розберемося з третім інтегралом в правій частині формули (23). Вираз у дужках цього інтеграла дорівнює дивергенції тензора Ейнштейна:

(24) \frac{1}{\sqrt{g}} \frac{\partial}{\partial u^{j}} (\sqrt{g} G^{[m - 1] i j}) + Γ_{s j}^{i} G^{[m - 1] s j} = \partial_{j} G^{[m - 1] i j} + Γ_{s j}^{i} G^{[m - 1] s j} + Γ_{s j}^{j} G^{[m - 1] i s} = \nabla_{j} G^{[m - 1] i j}

Але, як легко показати (використовуючи рівняння Петерсона-Кодацці), дивергенція тензора Ейнштейна тотожно дорівнює нулю.

В останньому інтегралі формули (24) похідні вектора нормалі до гіперповерхні дорівнюють:

(25) 𝐧_{j} = - b_{j}^{s} 𝐫_{s}

Отже формула (23) стає простішою:

(26) \int (δ K^{[m]}) \sqrt{g} = \int R_{i j}^{[m]} δ g^{i j} \sqrt{g} + \oint G^{[m - 1] i j} (𝐧 \cdot δ 𝐫_{i}) d S_{j} + \int \sqrt{g} G^{[m - 1] i j} b_{j}^{s} (𝐫_{s} \cdot δ 𝐫_{i})

Часткова згортка тензора Ейнштейна з тензором повної кривини дає тензор Річчі:

(27) G^{[m - 1] i j} b_{j}^{s} = R^{[m] i s}

і внаслідок симетрії тензора Річчі маємо:

(28) R^{[m] i j} (𝐫_{i} \cdot δ 𝐫_{j}) = \frac{1}{2} R^{[m] i j} ((𝐫_{i} \cdot δ 𝐫_{j}) + (δ 𝐫_{i} \cdot 𝐫_{j})) = \frac{1}{2} R^{[m] i j} δ g_{i j} = - \frac{1}{2} R_{i j}^{[m]} δ g^{i j}

Щодо останнього перетворення формули (28), дивіться статтю Прості обчислення диференціальної геометрії. Отже формула (26) ще раз спрощується:

(29) \int (δ K^{[m]}) \sqrt{g} = \frac{1}{2} \int R_{i j}^{[m]} δ g^{i j} \sqrt{g} + \oint G^{[m - 1] i j} (𝐧 \cdot δ 𝐫_{i}) d S_{j}

підставимо одержаний вираз в (5), також обчисливши варіацію кореня з визначника:

(30) δ \sqrt{g} = - \frac{1}{2} \sqrt{g} g_{i j} δ g^{i j}

В результаті одержуємо такий вираз для варіації інтеграла Ґауса:

(31) δ I_{m} = \frac{1}{2} \int (R_{i j}^{[m]} - K^{[m]} g_{i j}) \sqrt{g} δ g^{i j} + \oint G^{[m - 1] i j} (𝐧 \cdot δ 𝐫_{i}) d S_{j}

Для одержання формули (2) лишається тільки скористатися основним зв'язком між тензорами Ейнштейна та Річчі:

(32) G_{i j}^{[m]} = R_{i j}^{[m]} - K^{[m]} g_{i j}

Обчислення варіації інтеграла Ґауса для випадку парного степеня

Якщо степінь кривини Ґауса парний, то як кривина Ґауса так і тензор Ейнштейна виражаються через тензор Рімана, і ми можемо робити всі обчислення, користуючись тільки внутрішньою геометрією.

Як і раніше, варіація інтеграла Ґауса розкладається на два інтеграла:

(33) δ I_{2 m} = \int (δ K^{[2 m]}) \sqrt{g} d u^{1} \dots d u^{n} - \frac{1}{2} \int K^{[2 m]} g_{i j} δ g^{i j} \sqrt{g} d u^{1} \dots d u^{n} = J_{1} + J_{2}

Обчислимо спочатку перший інтеграл, що стоїть в правій частині рівності (33). Кривина Ґауса обчислюється за формулою:

(34) K^{[2 m]} = \frac{1}{2^{m} (2 m)!} g_{k_{1} l_{1} \dots k_{m} l_{m}}^{i_{1} j_{1} \dots i_{m} j_{m}} R_{i_{1} j_{1}}^{k_{1} l_{1}} \dots R_{i_{m} j_{m}}^{k_{m} l_{m}}

Оскільки коефіцієнти тензора метричної матрьошки є постійними числами, і «вертикальні пари» індексів тензора метричної матрьошки можна переставляти, то варіація кривини Ґауса (34) дорівнює:

(35) δ K^{[2 m]} = \frac{\partial K^{[2 m]}}{\partial R_{i j}^{k l}} δ R_{i j}^{k l} = c_{k l}^{i j} δ R_{i j}^{k l}

де тензор $c_{k l}^{i j}$ знаходиться за такою формулою (кількість множників на одиницю менша ніж у формулі (34)):

(36) c_{k l}^{i j} = \frac{m}{2^{m} (2 m)!} g_{k l k_{2} l_{2} \dots k_{m} l_{m}}^{i j i_{2} j_{2} \dots i_{m} j_{m}} R_{i_{2} j_{2}}^{k_{2} l_{2}} \dots R_{i_{m} j_{m}}^{k_{m} l_{m}}

Нам треба обчислити інтеграл:

(37) J_{1} = \int (δ K^{[2 m]}) d τ = \int c_{k l}^{i j} (δ R_{i j}^{k l}) d τ

Підінтегральний вираз можна розбити на два доданки таким чином:

(38) c_{k l}^{i j} δ R_{i j}^{k l} = c_{k l}^{i j} δ (g^{l s} R_{s i j}^{k}) = c_{k l}^{i j} g^{l s} δ R_{s i j}^{k} + c_{k l}^{i j} R_{s i j}^{k} δ g^{l s} =

= c_{s}^{i j k} δ R_{i j k}^{s} + c_{k_{1} i}^{i_{1} j_{1}} R_{j i_{1} j_{1}}^{k_{1}} δ g^{i j}

Останній доданок можна виразити через тензор Річчі $2 m$ -того степеня:

(39) c_{k_{1} i}^{i_{1} j_{1}} R_{i_{1} j_{1}}^{k_{1} j} = \frac{1}{2} \cdot \frac{1}{2^{m} (2 m - 1)!} g_{k_{1} i k_{2} l_{2} \dots k_{m} l_{m}}^{i_{1} j_{1} i_{2} j_{2} \dots i_{m} j_{m}} R_{i_{1} j_{1}}^{k_{1} j} R_{i_{2} j_{2}}^{k_{2} l_{2}} \dots R_{i_{m} j_{m}}^{k_{m} l_{m}} = \frac{1}{2} R_{i}^{[2 m] j}

Підставимо формули (37-39) в (33), маємо:

(40) δ I_{2 m} = J_{1} + J_{2} = \int (c_{s}^{i j k} δ R_{i j k}^{s} + \frac{1}{2} R_{i j}^{[2 m]} δ g^{i j} - \frac{1}{2} K^{[2 m]} g_{i j} δ g^{i j}) d τ

= \int c_{s}^{i j k} δ R_{i j k}^{s} d τ + \frac{1}{2} \int G_{i j}^{[2 m]} δ g^{i j} d τ

Ми знову маємо два інтеграла. Покажемо, що перший з них перетворюється в інтеграл по межі області інтегрування інтеграла Ґауса (а тому дорівнює нулю, якщо інтеграл Ґауса береться по всьому об'єму замкненого многовида).

Ці обчислення доволі складні, тому для зручності допоміжні обчислення винесено в окрему статтю. Скористаємося результатом цих обчислень:

(41) = \int_{Ω} c_{s}^{i j k} δ R_{i j k}^{s} d τ = \oint_{S} w_{s}^{i j k} δ Γ_{i j}^{s} n_{k} d S - \int_{Ω} (\nabla_{k} w_{s}^{i j k}) δ Γ_{i j}^{s} d τ

де позначено:

(42) w_{s}^{i j k} = c_{s}^{i k j} - c_{s}^{i j k} = - 2 c_{s}^{i j k}

(останнє перетворення записано з огляду на антисиметрію тензора $c_{s}^{i j k}$ по останній парі індексів)

Можна легко показати (аналогічно тому, як це було зроблено для тензора Ейнштейна), що дивергенція тензора $c_{s}^{i j k}$ , а отже і $w_{s}^{i j k}$ дорівнює нулю.

Дійсно, беручи дивергенцію від формули (36) наприклад за індексом $j$ (тобто $\nabla_{j} c_{k l}^{i j}$ ), одержуємо доданки виду:

(43) \frac{m}{2^{m} (2 m)!} g_{k l k_{2} l_{2} \dots k_{m} l_{m}}^{i j i_{2} j_{2} \dots i_{m} j_{m}} (\nabla_{j} R_{i_{2} j_{2}}^{k_{2} l_{2}}) \dots R_{i_{m} j_{m}}^{k_{m} l_{m}}

Якщо розглянути в формулі (43) три індекси $j, i_{2}, j_{2}$ то можна помітити, що тензор метричної матрьошки не змінюється при їх циклічній перестановці. Оскільки при додаванні за цими індексами зустрічаються усі три циклічні перестановки, то внаслідок диференціальної тотожності Біанкі сума буде дорівнювати нулю.

Остаточно маємо для варіації інтеграла Ґауса по деякій області $Ω$ многовиду:

(44) δ \int_{Ω} K^{[2 m]} d τ = \frac{1}{2} \int_{Ω} G_{i j}^{[2 m]} δ g^{i j} d τ - 2 \oint_{S} c_{s}^{i j k} δ Γ_{i j}^{s} d τ

Посилання

Інтеграли Гауса

Зміст

Варіація інтеграла Ґауса

Обчислення варіації інтеграла Ґауса для випадку гіперповерхні

Обчислення варіації інтеграла Ґауса для випадку парного степеня

Посилання

Навігаційне меню

Інтеграли Гауса

Варіація інтеграла Ґауса

Обчислення варіації інтеграла Ґауса для випадку гіперповерхні

Обчислення варіації інтеграла Ґауса для випадку парного степеня

Посилання

Навігаційне меню

Пошук