Некоторые векторные поля являются суммой безвихревого и соленоидального векторного поля.
В физике и математике теорема о разложении Гельмгольца или основная теорема векторного исчисления [1] [2] [3] [4] [5] [6] [7] утверждает, что некоторые дифференцируемые векторные поля могут быть разложены в сумму безвихревого ( безроторного ) векторного поля и соленоидального ( бездивергентного ) векторного поля. В физике часто обсуждается только разложение достаточно гладких , быстро затухающих векторных полей в трех измерениях. Она названа в честь Германа фон Гельмгольца .
Определение
Для векторного поля , определенного на домене , разложение Гельмгольца представляет собой пару векторных полей и таких, что:
Здесь — скалярный потенциал , — его градиент , а — дивергенция векторного поля . Безвихревое векторное поле называется градиентным полем и называется соленоидальным полем или полем вращения . Это разложение не существует для всех векторных полей и не является единственным . [8]
История
Разложение Гельмгольца в трех измерениях было впервые описано в 1849 году [9] Джорджем Габриэлем Стоксом для теории дифракции . Герман фон Гельмгольц опубликовал свою статью о некоторых основных гидродинамических уравнениях в 1858 году, [10] [11], которая была частью его исследования теорем Гельмгольца, описывающих движение жидкости вблизи вихревых линий. [11] Их вывод требовал, чтобы векторные поля достаточно быстро затухали на бесконечности. Позже это условие удалось ослабить, и разложение Гельмгольца можно было распространить на более высокие измерения. [8] [12] [13] Для римановых многообразий было выведено разложение Гельмгольца-Ходжа с использованием дифференциальной геометрии и тензорного исчисления . [8] [11] [14] [15]
Разложение стало важным инструментом для решения многих задач теоретической физики [11] [14] , но также нашло применение в анимации , компьютерном зрении , а также робототехнике [15] .
Трехмерное пространство
Многие учебники по физике ограничивают разложение Гельмгольца трехмерным пространством и ограничивают его применение векторными полями, которые достаточно быстро затухают на бесконечности, или функцией выпуклости , которая определена на ограниченной области . Затем можно определить векторный потенциал , так что поле вращения задается как , используя ротор векторного поля. [16]
Пусть будет векторным полем на ограниченной области , которая дважды непрерывно дифференцируема внутри , и пусть будет поверхностью, которая охватывает область . Тогда может быть разложена на компоненту без завихрений и компоненту без дивергенции следующим образом: [17]
где
и является оператором набла относительно , а не .
Если и, следовательно, неограниченно и исчезает быстрее, чем при , то имеем [18]
Это справедливо, в частности, если дважды непрерывно дифференцируемо по и имеет ограниченный носитель.
Вывод
ДоказательствоПредположим, что у нас есть векторная функция , ротор которой нам известен , и дивергенция, в области и полях на границе. Записывая функцию с использованием дельта-функции в виде,
где — векторный оператор Лапласа , имеем
где мы использовали векторное тождество Лапласа:
дифференцирование/интегрирование по и в последней строке линейность аргументов функции:
Затем, используя векторные тождества
мы получаем
Благодаря теореме о дивергенции уравнение можно переписать как
с нормалью к внешней поверхности .
Определение
мы наконец получаем
Пространство решений
Если — разложение Гельмгольца , то — другое разложение тогда и только тогда, когда
- и
- где
- — гармоническое скалярное поле ,
- является векторным полем, которое удовлетворяет
- скалярное поле.
Доказательство: Положим и . Согласно определению разложения Гельмгольца, условие эквивалентно
- .
Расхождение каждого члена этого уравнения дает , следовательно, является гармоническим.
Наоборот, при любой гармонической функции , является соленоидальной, поскольку
Таким образом, согласно предыдущему разделу, существует векторное поле такое, что .
Если — другое такое векторное поле, то
выполняется , следовательно,
для некоторого скалярного поля .
Поля с заданной дивергенцией и ротором
Термин «теорема Гельмгольца» может также относиться к следующему. Пусть C — соленоидальное векторное поле , а d — скалярное поле на R 3 , которые достаточно гладкие и которые исчезают быстрее, чем 1/ r 2 на бесконечности. Тогда существует векторное поле F такое, что
если дополнительно векторное поле F исчезает при r → ∞ , то F является уникальным. [18]
Другими словами, векторное поле может быть построено как с заданной дивергенцией, так и с заданным ротором, и если оно также исчезает на бесконечности, оно однозначно определяется своей дивергенцией и ротором. Эта теорема имеет большое значение в электростатике , поскольку уравнения Максвелла для электрических и магнитных полей в статическом случае имеют именно такой тип. [18] Доказательство проводится с помощью конструкции, обобщающей приведенную выше: положим
где представляет собой ньютоновский потенциальный оператор. (При действии на векторное поле, такое как ∇ × F , он определяется как действующий на каждый компонент.)
Слабая формулировка
Разложение Гельмгольца можно обобщить, уменьшив предположения о регулярности (необходимость существования сильных производных). Предположим, что Ω — ограниченная, односвязная, липшицева область . Каждое квадратично-интегрируемое векторное поле u ∈ ( L 2 (Ω)) 3 имеет ортогональное разложение: [19] [20] [21]
где φ находится в пространстве Соболева H 1 (Ω) квадратично интегрируемых функций на Ω , частные производные которых, определенные в смысле распределения , квадратично интегрируемы, и A ∈ H (rot, Ω) — пространство Соболева векторных полей, состоящее из квадратично интегрируемых векторных полей с квадратично интегрируемым ротором.
Для несколько более гладкого векторного поля u ∈ H (rot, Ω) справедливо аналогичное разложение:
где φ ∈ H 1 (Ω), v ∈ ( H 1 (Ω)) d .
Вывод из преобразования Фурье
Обратите внимание, что в теореме, сформулированной здесь, мы наложили условие, что если не определено на ограниченной области, то будет затухать быстрее, чем . Таким образом, преобразование Фурье от , обозначенное как , гарантированно существует. Мы применяем соглашение
Преобразование Фурье скалярного поля является скалярным полем, а преобразование Фурье векторного поля является векторным полем той же размерности.
Теперь рассмотрим следующие скалярные и векторные поля:
Следовательно
Продольные и поперечные поля
Терминология, часто используемая в физике, относится к компоненту без завихрений векторного поля как к продольной компоненте , а к компоненте без расхождения как к поперечной компоненте . [22] Эта терминология исходит из следующей конструкции: вычислить трехмерное преобразование Фурье векторного поля . Затем разложить это поле в каждой точке k на два компонента, один из которых указывает продольно, т. е. параллельно k , другой из которых указывает в поперечном направлении, т. е. перпендикулярно k . До сих пор мы имели
Теперь применим обратное преобразование Фурье к каждому из этих компонентов. Используя свойства преобразований Фурье, выводим:
Так как и ,
мы можем получить
так что это действительно разложение Гельмгольца. [23]
Обобщение на более высокие измерения
Матричный подход
Обобщение на размерности невозможно выполнить с помощью векторного потенциала, поскольку оператор вращения и векторное произведение определены (как векторы) только в трех измерениях.
Пусть — векторное поле в ограниченной области , которое убывает быстрее, чем при и .
Скалярный потенциал определяется аналогично трехмерному случаю как:
где ядро интегрирования снова является фундаментальным решением уравнения Лапласа , но в d-мерном пространстве:
с объемом d-мерных единичных шаров и гамма -функцией .
Для , просто равно , что дает тот же префактор, что и выше. Вращательный потенциал является антисимметричной матрицей с элементами:
Над диагональю находятся элементы, которые снова появляются в зеркальном отображении на диагонали, но с отрицательным знаком. В трехмерном случае элементы матрицы просто соответствуют компонентам векторного потенциала . Однако такой матричный потенциал может быть записан как вектор только в трехмерном случае, поскольку действителен только для .
Как и в трехмерном случае, градиентное поле определяется как
Поле вращения, с другой стороны, определяется в общем случае как расходимость строк матрицы:
В трехмерном пространстве это эквивалентно вращению векторного потенциала. [8] [24]
Тензорный подход
В -мерном векторном пространстве с , можно заменить соответствующей функцией Грина для Лапласа , определяемой как
, где для индекса используется соглашение Эйнштейна о суммировании . Например, в 2D.
Следуя тем же шагам, что и выше, мы можем записать,
где есть дельта Кронекера (и снова используется соглашение о суммировании). Вместо определения векторного Лапласа, использованного выше, мы теперь используем тождество для символа Леви-Чивиты ,
которое справедливо в размерностях, где есть -компонентный мультииндекс . Это дает
Поэтому мы можем записать,
где
Обратите внимание, что векторный потенциал заменяется ранговым тензором в размерностях.
Поскольку является функцией только , можно заменить , что дает Интегрирование по частям затем можно использовать для определения
того, где находится граница . Эти выражения аналогичны приведенным выше для трехмерного пространства.
Для дальнейшего обобщения на многообразия см. обсуждение разложения Ходжа ниже.
Дифференциальные формы
Разложение Ходжа тесно связано с разложением Гельмгольца, [25] обобщающим векторные поля на R 3 на дифференциальные формы на римановом многообразии M . Большинство формулировок разложения Ходжа требуют, чтобы M было компактным . [26] Поскольку это не так для R 3 , теорема о разложении Ходжа не является строго обобщением теоремы Гельмгольца. Однако ограничение компактности в обычной формулировке разложения Ходжа можно заменить подходящими предположениями о распаде на бесконечности для задействованных дифференциальных форм, что дает надлежащее обобщение теоремы Гельмгольца.
Расширения полей, не затухающих на бесконечности
Большинство учебников имеют дело только с векторными полями, затухающими быстрее, чем с на бесконечности. [16] [13] [27] Однако Отто Блюменталь в 1905 году показал, что адаптированное ядро интегрирования может быть использовано для интегрирования полей, затухающих быстрее, чем с , что существенно менее строго. Чтобы добиться этого, ядро в интегралах свертки должно быть заменено на . [28]
С еще более сложными ядрами интегрирования решения могут быть найдены даже для расходящихся функций, которые не должны расти быстрее, чем полиномиальные. [12] [13] [24] [29]
Для всех аналитических векторных полей, которые не должны стремиться к нулю даже на бесконечности, методы, основанные на частичном интегрировании и формуле Коши для повторного интегрирования [30], могут быть использованы для вычисления замкнутых решений вращательного и скалярного потенциалов, как в случае многомерных полиномиальных , синусоидальных , косинусных и экспоненциальных функций . [8]
Уникальность решения
В общем случае разложение Гельмгольца не определено однозначно. Гармоническая функция — это функция, которая удовлетворяет . Добавляя к скалярному потенциалу , можно получить другое разложение Гельмгольца:
Для векторных полей , затухающих на бесконечности, вполне правдоподобно, что скалярные и вращательные потенциалы также затухают на бесконечности. Поскольку является единственной гармонической функцией с этим свойством, что следует из теоремы Лиувилля , это гарантирует уникальность градиентных и вращательных полей. [31]
Эта уникальность не распространяется на потенциалы: в трехмерном случае скалярный и векторный потенциалы совместно имеют четыре компонента, тогда как векторное поле имеет только три. Векторные поля инвариантны к калибровочным преобразованиям, а выбор соответствующих потенциалов, известный как фиксация калибровки, является предметом калибровочной теории . Важными примерами из физики являются калибровочное условие Лоренца и калибровка Кулона . Альтернативой является использование полоидально-тороидального разложения .
Приложения
Электродинамика
Теорема Гельмгольца представляет особый интерес в электродинамике , поскольку ее можно использовать для записи уравнений Максвелла в потенциальном изображении и более простого их решения. Разложение Гельмгольца можно использовать для доказательства того, что при заданной плотности электрического тока и плотности заряда можно определить электрическое поле и плотность магнитного потока . Они являются уникальными, если плотности исчезают на бесконечности и предполагается то же самое для потенциалов. [16]
Динамика жидкости
В динамике жидкости проекция Гельмгольца играет важную роль, особенно для теории разрешимости уравнений Навье-Стокса . Если проекция Гельмгольца применяется к линеаризованным несжимаемым уравнениям Навье-Стокса, получается уравнение Стокса . Оно зависит только от скорости частиц в потоке, но больше не от статического давления, что позволяет свести уравнение к одному неизвестному. Однако оба уравнения, Стокса и линеаризованные уравнения, эквивалентны. Оператор называется оператором Стокса . [32]
Теория динамических систем
В теории динамических систем разложение Гельмгольца может быть использовано для определения «квазипотенциалов», а также для вычисления функций Ляпунова в некоторых случаях. [33] [34] [35]
Для некоторых динамических систем, таких как система Лоренца ( Эдвард Н. Лоренц , 1963 [36] ), упрощенная модель атмосферной конвекции , можно получить замкнутое выражение разложения Гельмгольца:
Разложение Гельмгольца со скалярным потенциалом задается как:
Квадратичный скалярный потенциал обеспечивает движение в направлении начала координат, что отвечает за устойчивую точку фиксации для некоторого диапазона параметров. Для других параметров поле вращения обеспечивает создание странного аттрактора , заставляя модель проявлять эффект бабочки . [8] [37]
Медицинская визуализация
В магнитно-резонансной эластографии , варианте МРТ, где механические волны используются для исследования вязкоупругости органов, разложение Гельмгольца иногда используется для разделения измеренных полей смещения на его компонент сдвига (без расходимости) и компонент сжатия (без завитков). [38] Таким образом, комплексный модуль сдвига может быть рассчитан без учета вкладов волн сжатия.
Компьютерная анимация и робототехника
Разложение Гельмгольца также используется в области компьютерной инженерии. Это включает в себя робототехнику, реконструкцию изображений, а также компьютерную анимацию, где разложение используется для реалистичной визуализации жидкостей или векторных полей. [15] [39]
Смотрите также
Примечания
- ↑ Дэниел Александр Мюррей : Элементарный курс интегрального исчисления . American Book Company, 1898. С. 8.
- ^ Дж. У. Гиббс , Эдвин Бидвелл Уилсон : Векторный анализ . 1901, стр. 237, ссылка из интернет-архива .
- ↑ Оливер Хевисайд : Электромагнитная теория . Том 1, издательская компания «Электрик», 1893.
- ^ Уэсли Стокер Баркер Вулхаус : Элементы дифференциального исчисления . Уил, 1854.
- ^ Уильям Вулси Джонсон : Элементарный трактат по интегральному исчислению: основанный на методе скоростей или флюксий . John Wiley & Sons, 1881.
См. также: Метод флюксий . - ↑ Джеймс Берни Шоу: Вектор исчисления: с приложениями к физике . Д. Ван Ностранд, 1922, стр. 205.
См. также: Теорема Грина . - ↑ Джозеф Эдвардс: Трактат об интегральном исчислении . Том 2. Chelsea Publishing Company, 1922.
- ^ abcdef Эрхард Глётцль, Оливер Рихтерс: Разложение Гельмгольца и потенциальные функции для n-мерных аналитических векторных полей . В: Журнал математического анализа и приложений 525(2), 127138, 2023, doi :10.1016/j.jmaa.2023.127138, arXiv :2102.09556v3. Рабочий лист Mathematica в doi :10.5281/zenodo.7512798.
- ↑ Джордж Габриэль Стокс : О динамической теории дифракции . В: Труды Кембриджского философского общества 9, 1849, стр. 1–62. doi :10.1017/cbo9780511702259.015, см. стр. 9–10.
- ^ Герман фон Гельмгольц : Über Integrale der Hydrodynamischen Gleichungen, welche den Wirbelbewegungen entsprechen . В: Journal für die reine und angewandte Mathematik 55, 1858, стр. 25–55, номер документа : 10.1515/crll.1858.55.25 (sub.uni-goettingen.de, digizeitschriften.de). На странице 38 компоненты скорости жидкости ( u , v , w ) выражены через градиент скалярного потенциала P и ротор векторного потенциала ( L , M , N ).
- ^ abcd Альп Кустепели: О теореме Гельмгольца и ее обобщении для многослойных сред . В: Электромагнетизм 36.3, 2016, стр. 135–148, doi :10.1080/02726343.2016.1149755.
- ^ ab Ton Tran-Cong: О теореме разложения Гельмгольца и уравнении Пуассона с бесконечной областью определения . В: Quarterly of Applied Mathematics 51.1, 1993, стр. 23–35, JSTOR 43637902.
- ^ abc D. Petrascheck, R. Folk: Теорема разложения Гельмгольца и расширение Блюменталя с помощью регуляризации . В: Физика конденсированных сред 20(1), 13002, 2017, doi :10.5488/CMP.20.13002.
- ^ ab Вольфганг Шпрессиг: О разложениях Гельмгольца и их обобщениях – Обзор . В: Математические методы в прикладных науках 33.4, 2009, стр. 374–383, doi :10.1002/mma.1212.
- ^ abc Харш Бхатия, Грегори Норгард, Валерио Паскуччи, Пир-Тимо Бремер: Разложение Гельмгольца-Ходжа – Обзор . В: Труды IEEE по визуализации и компьютерной графике 19.8, 2013, стр. 1386–1404, doi :10.1109/tvcg.2012.316.
- ^ abc Дитмар Петрашек: Повторное рассмотрение разложения Гельмгольца . В: European Journal of Physics 37.1, 2015, Статья 015201, doi : 10.1088/0143-0807/37/1/015201.
- ^ "Теорема Гельмгольца" (PDF) . Университет Вермонта. Архивировано из оригинала (PDF) 2012-08-13 . Получено 2011-03-11 .
- ^ abc Дэвид Дж. Гриффитс : Введение в электродинамику . Prentice-Hall, 1999, стр. 556.
- ^ Шериф Амруш, Кристин Бернарди , Моник Дож , Виветт Жиро : Векторные потенциалы в трехмерных негладких областях . В: Математические методы в прикладных науках 21(9), 1998, стр. 823–864, doi :10.1002/(sici)1099-1476(199806)21:9<823::aid-mma976>3.0.co;2-b, Bibcode :/abstract 1998MMAS...21..823A .
- ^ Р. Дотрей и Ж.-Л. Лионс. Спектральная теория и ее приложения, том 3 Математического анализа и численных методов для науки и техники. Springer-Verlag, 1990.
- ^ V. Girault , PA Raviart: Методы конечных элементов для уравнений Навье–Стокса: теория и алгоритмы. Springer Series in Computational Mathematics. Springer-Verlag, 1986.
- ^ AM Stewart: Продольные и поперечные компоненты векторного поля . В: Sri Lankan Journal of Physics 12, стр. 33–42, 2011, doi :10.4038/sljp.v12i0.3504 arXiv :0801.0335
- ^ Роберт Литтлджон: Гамильтониан классического электромагнитного поля. Онлайн-лекции, berkeley.edu.
- ^ ab Эрхард Глётцль, Оливер Рихтерс: Разложение Гельмгольца и вращательные потенциалы в n-мерных декартовых координатах . 2020, arXiv :2012.13157.
- ^ Фрэнк У. Уорнер: Теорема Ходжа . В: Основания дифференцируемых многообразий и групп Ли . (= Graduate Texts in Mathematics 94). Springer, Нью-Йорк, 1983, doi :10.1007/978-1-4757-1799-0_6.
- ^ Кантарелла, Джейсон; ДеТурк, Деннис; Глюк, Герман (2002). «Векторное исчисление и топология доменов в 3-пространстве». The American Mathematical Monthly . 109 (5): 409–442. doi :10.2307/2695643. JSTOR 2695643.
- ^ Р. Дуглас Грегори: Теорема Гельмгольца, когда область бесконечна и когда поле имеет особые точки . В: The Quarterly Journal of Mechanics and Applied Mathematics 49.3, 1996, стр. 439–450, doi :10.1093/qjmam/49.3.439.
- ^ Отто Блюменталь : Über die Zerlegung unendlicher Vektorfelder . В: Mathematische Annalen 61.2, 1905, стр. 235–250, doi : 10.1007/BF01457564.
- ^ Мортон Э. Гуртин: О теореме Гельмгольца и полноте функций напряжения Папковича-Нейбера для бесконечных областей . В: Архив рациональной механики и анализа 9.1, 1962, стр. 225–233, doi :10.1007/BF00253346.
- ^ Коши, Огюстен-Луи (1823). «Тренте-Чинквием Лесон». Резюме уроков, проводимых в Королевской политехнической школе по исчислению бесконечно малых величин (на французском языке). Париж: Королевская империя. стр. 133–140.
- ^ Шелдон Акслер, Пол Бурдон, Уэйд Рэми: Ограниченные гармонические функции . В: Теория гармонических функций (= Graduate Texts in Mathematics 137). Springer, Нью-Йорк, 1992, стр. 31–44, doi :10.1007/0-387-21527-1_2.
- ^ Александр Дж. Хорин, Джерролд Э. Марсден: Математическое введение в механику жидкости (= Тексты по прикладной математике 4). Springer US, Нью-Йорк 1990, doi :10.1007/978-1-4684-0364-0.
- ^ Томохару Суда: Построение функций Ляпунова с использованием разложения Гельмгольца–Ходжа . В: Дискретные и непрерывные динамические системы – A 39.5, 2019, стр. 2437–2454, doi :10.3934/dcds.2019103.
- ^ Томохару Суда: Применение разложения Гельмгольца–Ходжа к изучению некоторых векторных полей . В: Journal of Physics A: Mathematical and Theoretical 53.37, 2020, стр. 375703. doi :10.1088/1751-8121/aba657.
- ^ Джозеф Сюй Чжоу, MDS Алию, Эрик Аурелл, Суй Хуан: Квазипотенциальный ландшафт в сложных мультистабильных системах . В: Журнал интерфейса Королевского общества 9.77, 2012, стр. 3539–3553, doi :10.1098/rsif.2012.0434.
- ^ Эдвард Н. Лоренц : Детерминированный непериодический поток . В: Журнал атмосферных наук 20.2, 1963, стр. 130–141, doi :10.1175/1520-0469(1963)020<0130:DNF>2.0.CO;2.
- ^ Хайнц-Отто Пейтген, Хартмут Юргенс, Дитмар Саупе: Странные аттракторы: Локус хаоса . В: Хаос и фракталы . Спрингер, Нью-Йорк, стр. 655–768. дои : 10.1007/978-1-4757-4740-9_13.
- ^ Армандо Мандука: МР-эластография: принципы, руководства и терминология . В: Магнитный резонанс в медицине , 2021, doi :10.1002/mrm.28627 PMID 33296103.
- ^ Херш Бхатия, Валерио Паскуччи, Пир-Тимо Бремер: Естественное разложение Гельмгольца-Ходжа для анализа потоков с открытой границей . В: Труды IEEE по визуализации и компьютерной графике 20.11, ноябрь 2014 г., стр. 1566–1578, ноябрь 2014 г., doi : 10.1109/TVCG.2014.2312012.
Ссылки
- Джордж Б. Арфкен и Ханс Дж. Вебер, Математические методы для физиков , 4-е издание, Academic Press: Сан-Диего (1995) стр. 92–93
- Джордж Б. Арфкен и Ханс Дж. Вебер, Математические методы для физиков – Международное издание , 6-е издание, Academic Press: Сан-Диего (2005) стр. 95–101
- Резерфорд Арис , Векторы, тензоры и основные уравнения механики жидкости , Prentice-Hall (1962), OCLC 299650765, стр. 70–72