Русская Википедия:Тензор напряжений Максвелла
Тензор напряжений Максвелла (назван в честь Джеймса Клерка Максвелла) представляет собой симметричный тензор второго порядка, используемый в классическом электромагнетизме для представления взаимодействия между электромагнитными силами и механическим импульсом. В простых случаях, таких как точечный заряд, свободно движущийся в однородном магнитном поле, легко рассчитать силы, действующие на заряд, согласно силе Лоренца. В более сложных случаях такая обычная процедура может стать непрактично сложной с уравнениями, охватывающими несколько строк. Поэтому удобно собрать многие из этих членов в тензоре напряжений Максвелла и использовать тензорную арифметику, чтобы найти ответ на поставленную задачу.
В релятивистской формулировке электромагнетизма тензор Максвелла появляется как часть электромагнитного тензора энергии-импульса, который является электромагнитной составляющей полного тензора энергии-импульса. Последний описывает плотность и поток энергии и импульса в пространстве-времени.
Обоснование
Ниже показано, что электромагнитная сила записывается параметрами E и B. Используя векторное исчисление и уравнения Максвелла, ищется симметрия в выражениях, содержащих E и B, а введение тензора напряжений Максвелла упрощает результат.
| Имя | Дифференциальная форма |
|---|---|
| Закон Гаусса (в вакууме) | <math>\nabla \cdot \mathbf{E} = \frac {\rho} {\varepsilon_0}</math> |
| Закон Гаусса для магнетизма | <math>\nabla \cdot \mathbf{B} = 0</math> |
| Уравнение Максвелла – Фарадея (закон индукции Фарадея) |
<math>\nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}} {\partial t}</math> |
| Круговой закон Ампера (в вакууме) (с поправкой Максвелла) |
<math>\nabla \times \mathbf{B} = \mu_0 \mathbf{J} + \mu_0 \varepsilon_0 \frac{\partial \mathbf{E}} {\partial t}\ </math> |
Шаблон:Ordered listВ приведённом выше соотношении сохранения импульса, <math>\nabla \cdot \boldsymbol{\sigma}</math> является плотностью потока импульса и играет роль, аналогичную <math>\mathbf{S}</math> в теореме Пойнтинга.
Приведённый выше вывод предполагает полное знание параметров ρ и J (как свободных, так и ограниченных зарядов и токов). В случае нелинейных материалов (таких как магнитное железо с BH-кривой (кривой плотности магнитного потока)) необходимо использовать нелинейный тензор напряжений Максвелла. [1]
Уравнение
В физике тензор напряжений Максвелла является тензором напряжений электромагнитного поля. Как указано выше в единицах СИ, это определяется как:
- <math>\sigma_{ij} =
\varepsilon_0 E_i E_j + \frac{1}{\mu_0}B_i B_j - \frac{1}{2}\left(\varepsilon_0 E^2 + \frac{1}{\mu_0}B^2\right)\delta_{ij},
</math>
где ε0 — электрическая постоянная, μ0 — магнитная постоянная, E — электрическое поле, B — магнитное поле, а δij — дельта Кронекера . В гауссовых единицах СГС это определяется как:
- <math>\sigma_{ij} = \frac{1}{4\pi}\left(E_i E_j + H_i H_j - \frac{1}{2}\left(E^2 + H^2\right)\delta_{ij}\right),</math>
где H — намагничивающее поле.
Альтернативный способ выражения этого тензора:
- <math>
\overset{\leftrightarrow}{\boldsymbol{\sigma}} = \frac{1}{4\pi} \left[ \mathbf{E} \otimes \mathbf{E} + \mathbf{H} \otimes \mathbf{H} - \frac{E^2 + H^2}{2}\mathbb{I} \right], </math>
где ⊗ — диадическое произведение, а последний тензор <math>\mathbb{I}</math> — единичная диада:
- <math>\mathbb{I} \equiv \begin{pmatrix}
1 & 0 & 0 \\
0 & 1 & 0 \\
0 & 0 & 1
\end{pmatrix} =
\left(\mathbf{\hat x} \otimes \mathbf{\hat x} + \mathbf{\hat y} \otimes \mathbf{\hat y} + \mathbf{\hat z} \otimes \mathbf{\hat z}\right)
</math>
Элемент ij тензора напряжений Максвелла имеет единицы импульса на единицу площади в единицу времени и даёт поток импульса, параллельный i-й оси, пересекающий поверхность, перпендикулярную j-й оси (в отрицательном направлении) в единицу времени.
Эти единицы также можно рассматривать как единицы силы на единицу площади (отрицательное давление), а элемент ij тензора также можно интерпретировать как силу, параллельную оси i, действующую на поверхность, перпендикулярную оси j, на единицу. площади. Действительно, диагональные элементы задают натяжение (напряжение, вытягивание), действующее на дифференциальный элемент площади по нормали к соответствующей оси. В отличие от сил, вызванных давлением идеального газа, элемент площади в электромагнитном поле также испытывает силу, направленную не по нормали к элементу. Этот сдвиг задается недиагональными элементами тензора напряжений.
Только магнетизм
Если поле является только магнитным (что в значительной степени верно, например, для двигателей), некоторые члены выпадают, и уравнение в единицах СИ принимает вид:
- <math>\sigma_{ij} = \frac{1}{\mu_0} B_i B_j - \frac{1}{2\mu_0} B^2 \delta_{ij} \,.</math>
Для цилиндрических объектов, таких как ротор двигателя, это выражение упрощается до:
- <math>\sigma_{rt} = \frac{1}{\mu_0} B_r B_t - \frac{1}{2\mu_0} B^2 \delta_{rt} \,.</math>
где r — сдвиг в радиальном (наружу от цилиндра) направлении, t — сдвиг в тангенциальном (вокруг цилиндра) направлении. Это тангенциальная сила, которая вращает двигатель. Br — плотность потока в радиальном направлении, а Bt — плотность потока в тангенциальном направлении.
В электростатике
В электростатике эффекты магнетизма отсутствуют. В этом случае магнитное поле исчезает, <math>\mathbf{B} = \mathbf{0}</math>, и мы получаем тензор электростатических напряжений Максвелла . Он дается в виде компонентов
- <math>
\sigma_{ij} = \varepsilon_0 E_i E_j - \frac{1}{2}\varepsilon_0 E^2\delta_{ij} </math>
и в символической форме
- <math>
\boldsymbol{\sigma} = \varepsilon_0\mathbf{E} \otimes \mathbf{E} - \frac{1}{2}\varepsilon_0(\mathbf{E} \cdot \mathbf{E})\mathbf{I}, </math>
где <math>\mathbf{I}</math> является подходящим тождественным тензором (обычно <math>3\times3</math>).
Собственное значение
Собственные значения тензора напряжений Максвелла определяются выражением:
- <math>\{\lambda\} = \left\{-\left(\frac{\varepsilon_0}{2}E^2 + \frac{1}{2\mu_0}B^2\right),~ \pm\sqrt{\left(\frac{\varepsilon_0}{2}E^2 - \frac{1}{2\mu_0}B^2\right)^2 + \frac{\varepsilon_0}{\mu_0}\left(\boldsymbol{E} \cdot \boldsymbol{B}\right)^2}\right\}.</math>
Эти собственные значения получаются итеративным применением леммы об определителе матрицы в сочетании с формулой Шермана-Моррисона .
Отмечая, что матрица характеристического уравнения, <math>\overleftrightarrow{\boldsymbol{\sigma}} - \lambda\mathbf{\mathbb{I}}</math> может записываться как
- <math>\overleftrightarrow{\boldsymbol{\sigma}} - \lambda\mathbf{\mathbb{I}} = -\left(\lambda + V\right)\mathbf{\mathbb{I}} + \varepsilon_0\mathbf{E}\mathbf{E}^\textsf{T} + \frac{1}{\mu_0}\mathbf{B}\mathbf{B}^\textsf{T},</math>
где
- <math>V = \frac{1}{2}\left(\varepsilon_0 E^2 + \frac{1}{\mu_0}B^2\right),</math>
мы устанавливаем
- <math>\mathbf{U} = -\left(\lambda + V\right)\mathbf{\mathbb{I}} + \varepsilon_0\mathbf{E}\mathbf{E}^\textsf{T}.</math>
Применяя лемму об определителе матрицы один раз, мы получаем:
- <math>\det{\left(\overleftrightarrow{\boldsymbol{\sigma}} - \lambda\mathbf{\mathbb{I}}\right)} = \left(1 + \frac{1}{\mu_0}\mathbf{B}^\textsf{T}\mathbf{U}^{-1}\mathbf{B}\right)\det{\left(\mathbf{U}\right)}.</math>
Применение его снова даёт,
- <math>\det{\left(\overleftrightarrow{\boldsymbol{\sigma}} - \lambda\mathbf{\mathbb{I}}\right)} = \left(1 + \frac{1}{\mu_0}\mathbf{B}^\textsf{T}\mathbf{U}^{-1}\mathbf{B}\right) \left(1 - \frac{\varepsilon_0\mathbf{E}^\textsf{T}\mathbf{E}}{\lambda + V}\right) \left(-\lambda - V\right)^3.</math>
Из последнего множимого из правой части выражения сразу видно, что <math>\lambda = -V</math> является одним из собственных значений.
Чтобы найти обратную <math>\mathbf{U}</math>, воспользуемся формулой Шермана-Моррисона:
- <math>\mathbf{U}^{-1} = -\left(\lambda + V\right)^{-1} - \frac{\varepsilon_0 \mathbf{E}\mathbf{E}^\textsf{T}}{\left(\lambda + V\right)^2 - \left(\lambda + V\right) \varepsilon_0\mathbf{E}^\textsf{T}\mathbf{E}}.</math>
Факторизовав <math>\left(-\lambda - V \right)</math> членом определителя, нам осталось найти нули рациональной функции:
- <math>\left(-\left(\lambda + V\right) - \frac{\varepsilon_0\left(\mathbf{E} \cdot \mathbf{B}\right)^2}{\mu_0\left(-\left(\lambda + V\right) + \varepsilon_0\mathbf{E}^\textsf{T}\mathbf{E}\right)}\right) \left(-\left(\lambda + V\right) + \varepsilon_0\mathbf{E}^\textsf{T} \mathbf{E}\right).</math>
Таким образом, как только мы решим
- <math>-\left(\lambda + V\right) \left(-\left(\lambda + V\right) + \varepsilon_0 E^2\right) - \frac{\varepsilon_0}{\mu_0}\left(\mathbf{E} \cdot \mathbf{B}\right)^2 = 0,</math>
мы получаем два других собственных значения.
Смотрите также
- Исчисление Риччи
- Плотность энергии электрических и магнитных полей
- Вектор пойнтинга
- Электромагнитный тензор энергии-импульса
- Магнитное давление
- Сила магнитного натяжения
Ссылки
- Дэвид Дж. Гриффитс, «Введение в электродинамику», стр. 351–352, Benjamin Cummings Inc., 2008 г.
- Джон Дэвид Джексон, «Классическая электродинамика, 3-е изд.», John Wiley & Sons, Inc., 1999.
- Ричард Беккер, «Электромагнитные поля и взаимодействия», Dover Publications Inc., 1964.
