WikiDer > Метод гомотопического анализа

Homotopy analysis method
Два пунктирных пути, показанные выше, гомотопны относительно их конечных точек. Анимация представляет одну возможную гомотопию.

В метод гомотопического анализа (ВЕТЧИНА) представляет собой полуаналитический метод решения нелинейный обычный/частичный дифференциальные уравнения. В методе гомотопического анализа используется концепция гомотопия из топология для генерации решения сходящегося ряда для нелинейных систем. Это возможно благодаря использованию гомотопическогоСерия Маклорена иметь дело с нелинейностями в системе.

HAM был впервые разработан в 1992 г. Ляо Шицзюнь из Шанхайский университет Цзяотун в его докторской диссертации[1] и далее модифицированный[2] в 1997 г.[рекламный язык] ввести ненулевой вспомогательный параметр, называемый параметр контроля сходимости, c0, построить гомотопию на дифференциальной системе общего вида.[3] Параметр управления сходимостью - это нефизическая переменная, которая обеспечивает простой способ проверки и обеспечения сходимости ряда решений. Способность HAM естественным образом демонстрировать сходимость решения ряда необычна для аналитических и полуаналитических подходов к нелинейным уравнениям в частных производных.

Характеристики

HAM отличается от других аналитические методы в четырех важных аспектах. Во-первых, это серии метод расширения, не зависящий напрямую от малых или больших физических параметров. Таким образом, он применим не только для слабонелинейных, но и для сильно нелинейных задач, выходя за рамки некоторых присущих стандартным ограничениям. методы возмущения. Во-вторых, HAM - это единый метод для Ляпунов метод искусственного малого параметра, метод дельта-разложения, метод Метод разложения Адомиана,[4] и гомотопический метод возмущения.[5][6] Большая общность метода часто допускает сильную сходимость решения в более крупных пространственных областях и областях параметров. В-третьих, HAM дает отличную гибкость в выражении решения и способах его получения в явном виде. Это дает большую свободу выбора базисные функции искомого решения и соответствующего вспомогательного линейный оператор гомотопии. Наконец, в отличие от других методов аналитической аппроксимации, HAM предоставляет простой способ гарантировать конвергенция серии решений.

Метод гомотопического анализа также может сочетаться с другими методами, используемыми в нелинейных дифференциальных уравнениях, таких как спектральные методы[7] и Аппроксимации Паде. Кроме того, его можно комбинировать с вычислительными методами, такими как метод граничных элементов разрешить линейный метод для решения нелинейных систем. В отличие от численной техники продолжение гомотопии, метод гомотопического анализа представляет собой метод аналитической аппроксимации, в отличие от метода дискретных вычислений. Кроме того, HAM использует параметр гомотопии только на теоретическом уровне, чтобы продемонстрировать, что нелинейная система может быть разбита на бесконечный набор линейных систем, которые решаются аналитически, в то время как методы продолжения требуют решения дискретной линейной системы при изменении параметра гомотопии. решить нелинейную систему.

Приложения

За последние двадцать лет HAM применялся для решения растущего числа нелинейных обычный/уравнения в частных производных в области науки, финансов и техники.[8][9] Например, множественные установившиеся резонансные волны на большой и конечной глубине воды[10] были найдены с волновой резонанс критерий произвольного количества путешествий гравитационные волны; это согласуется с критерием Филлипса для четырех волн малой амплитуды. Далее, единая волновая модель, применяемая с HAM,[11] допускает не только традиционные гладкие прогрессивные периодические / уединенные волны, но также прогрессивные уединенные волны с острым гребнем на конечной глубине воды. Эта модель показывает, что пиковые уединенные волны являются согласованными решениями наряду с известными гладкими. Кроме того, HAM был применен ко многим другим нелинейным задачам, таким как нелинейные теплопередача,[12] то предельный цикл нелинейных динамических систем,[13] Американец пут опцион,[14] точный Уравнение Навье – Стокса,[15] цена опциона согласно стохастическая волатильность,[16] то электрогидродинамический потоки[17] то Уравнение Пуассона – Больцмана. для полупроводниковых приборов,[18] и другие.

Краткое математическое описание

Изотопия чашки кофе в пончик (тор).

Рассмотрим общее нелинейное дифференциальное уравнение

,

куда - нелинейный оператор. Позволять обозначим вспомогательный линейный оператор, ты0(Икс) первоначальное предположение о ты(Икс), и c0 константа (называемая параметром контроля сходимости) соответственно. Использование параметра внедрения q ∈ [0,1] из теории гомотопий можно построить семейство уравнений,

называется уравнением деформации нулевого порядка, решение которого непрерывно изменяется относительно параметра вложения q ∈ [0,1]. Это линейное уравнение

с известным исходным предположением U(Икс; 0) = ты0(Икс) когда q = 0, но эквивалентно исходному нелинейному уравнению , когда q = 1, т.е. U(Икс; 1) = ты(Икс)). Следовательно, как q увеличивается от 0 до 1, решение U(Икс; q) уравнения деформации нулевого порядка изменяется (или деформируется) от выбранного первоначального предположения ты0(Икс) к решению ты(Икс) рассматриваемого нелинейного уравнения.

Расширение U(Икс; q) в сериале Тейлора о q = 0, имеем гомотопический ряд Маклорена

Предполагая, что так называемый параметр контроля сходимости c0 уравнения деформации нулевого порядка правильно выбрано так, чтобы указанный ряд сходился при q = 1, имеем решение гомотопического ряда

Из уравнения деформации нулевого порядка можно напрямую вывести основное уравнение тым(Икс)

называется мthуравнение деформации порядка, где и за k > 1, а правая часть рм зависит только от известных результатов ты0, ты1, ..., тым − 1 и может быть легко получен с помощью программного обеспечения компьютерной алгебры. Таким образом, исходное нелинейное уравнение преобразуется в бесконечное количество линейных, но без предположения о каких-либо малых / больших физических параметрах.

Поскольку HAM основан на гомотопии, у каждого есть большая свобода выбора начального предположения. ты0(Икс) вспомогательный линейный оператор , а параметр контроля сходимости c0 в уравнении деформации нулевого порядка. Таким образом, HAM предоставляет математикам свободу выбора типа уравнения для уравнения деформации высокого порядка и базовых функций его решения. Оптимальное значение параметра контроля сходимости c0 определяется минимумом квадрата остаточной ошибки определяющих уравнений и / или граничных условий после решения общей формы для выбранного начального предположения и линейного оператора. Таким образом, параметр контроля сходимости c0 является простым способом гарантировать сходимость решения гомотопического ряда и отличает HAM от других аналитических методов аппроксимации. В целом метод дает полезное обобщение концепции гомотопии.

Радиолюбитель и компьютерная алгебра

HAM - это метод аналитической аппроксимации, разработанный для компьютерной эры с целью «вычислений с функциями вместо чисел». В сочетании с системой компьютерной алгебры, такой как Mathematica или же Клен, можно получить аналитические аппроксимации сильно нелинейной задачи сколь угодно высокого порядка с помощью HAM всего за несколько секунд. Вдохновленный недавними успешными применениями HAM в различных областях, пакет Mathematica на основе HAM, названный BVPh, стал доступен в Интернете для решения нелинейных краевых задач. [4]. BVPh - это пакет решателя для сильно нелинейных ОДУ с особенностями, множественными решениями и многоточечными граничными условиями в конечном или бесконечном интервале, и включает поддержку некоторых типов нелинейных уравнений в частных производных.[8] Другой код Mathematica, основанный на HAM, APOh, был создан для решения явного аналитического приближения оптимальной границы исполнения американского пут-опциона, который также доступен в Интернете. [5].

Анализ частотной характеристики нелинейных осцилляторов

Недавно сообщалось, что HAM можно использовать для получения аналитических решений нелинейных уравнений частотной характеристики. Такие решения способны улавливать различные нелинейные поведения, такие как тип упрочнения, смягчения или смешанное поведение осциллятора.[19][20] Эти аналитические уравнения также полезны для предсказания хаоса в нелинейных системах.[21]

Рекомендации

  1. ^ Ляо, С.Дж. (1992), Предлагаемая методика гомотопического анализа для решения нелинейных задач., Кандидатская диссертация, Шанхайский университет Цзяо Тонг
  2. ^ Ляо, С.Дж. (1999), "Явная, полностью аналитическая аппроксимация задач вязкого течения Блазиуса", Международный журнал нелинейной механики, 34 (4): 759–778, Bibcode:1999IJNLM..34..759L, Дои:10.1016 / S0020-7462 (98) 00056-0
  3. ^ Ляо, С.Дж. (2003), За пределами возмущения: введение в метод гомотопического анализа, Бока-Ратон: Chapman & Hall / CRC Press, ISBN 978-1-58488-407-1[1]
  4. ^ Адомян Г. (1994). Решение пограничных задач физики: метод декомпозиции. Kluwer Academic Publishers.
  5. ^ Лян, Сунсинь; Джеффри, Дэвид Дж. (2009), "Сравнение метода гомотопического анализа и метода гомотопических возмущений через уравнение эволюции", Коммуникации в нелинейной науке и численном моделировании, 14 (12): 4057–4064, Bibcode:2009CNSNS..14.4057L, Дои:10.1016 / j.cns.2009.02.016
  6. ^ Sajid, M .; Хаят, Т. (2008), "Сравнение методов HAM и HPM в нелинейных уравнениях теплопроводности и конвекции", Нелинейный анализ: приложения в реальном мире, 9 (5): 2296–2301, Дои:10.1016 / j.nonrwa.2007.08.007
  7. ^ Motsa, S.S .; Сибанда, П .; Awad, F.G .; Шатей, С. (2010), "Новый метод спектрально-гомотопического анализа для задачи МГД Джеффри – Хамеля", Компьютеры и жидкости, 39 (7): 1219–1225, Дои:10.1016 / j.compfluid.2010.03.004
  8. ^ а б Ляо, С.Дж. (2012), Метод гомотопического анализа в нелинейных дифференциальных уравнениях., Берлин и Пекин: Springer & Higher Education Press, ISBN 978-7-04-032298-9 [2]
  9. ^ Ваджравелу, К .; Ван Гордер (2013), Явления нелинейного течения и гомотопический анализ, Берлин и Пекин: Springer & Higher Education Press, ISBN 978-3-642-32102-3 [3]
  10. ^ Xu, D.L .; Lin, Z.L .; Liao, S.J .; Стиассние, М. (2012), «О стационарных полностью резонансных прогрессивных волнах в воде конечной глубины», Журнал гидромеханики, 710: 379–418, Bibcode:2012JFM ... 710..379X, Дои:10.1017 / jfm.2012.370
  11. ^ Ляо, С.Дж. (2013), «Действительно ли существуют пиковые уединенные волны на воде?», Коммуникации в нелинейной науке и численном моделировании, 19 (6): 1792–1821, arXiv:1204.3354, Bibcode:2014CNSNS..19.1792L, Дои:10.1016 / j.cnsns.2013.09.042
  12. ^ Аббасбанди, С. (2006), "Применение метода гомотопического анализа к нелинейным уравнениям, возникающим при теплопередаче", Письма о физике A, 360 (1): 109–113, Bibcode:2006ФЛА..360..109А, Дои:10.1016 / j.physleta.2006.07.065
  13. ^ Chen, Y.M .; Лю, J.K. (2009), "Равномерно допустимое решение предельного цикла уравнения Дуффинга – Ван дер Поля", Сообщения об исследованиях в области механики, 36 (7): 845–850, Дои:10.1016 / j.mechrescom.2009.06.001
  14. ^ Чжу, С.П. (2006), "Точное и явное решение для оценки американских пут-опционов", Количественные финансы, 6 (3): 229–242, Дои:10.1080/14697680600699811
  15. ^ Туркилмазоглу, М. (2009), «Чисто аналитические решения потока сжимаемого пограничного слоя, обусловленного пористым вращающимся диском с теплопередачей», Физика жидкостей, 21 (10): 106104–106104–12, Bibcode:2009ФФЛ ... 21дж6104Т, Дои:10.1063/1.3249752
  16. ^ Пак, Сан Хён; Ким, Чон-Хун (2011), «Метод гомотопического анализа для определения цены опционов при стохастической волатильности», Письма по прикладной математике, 24 (10): 1740–1744, Дои:10.1016 / j.aml.2011.04.034
  17. ^ Мастроберардино, А. (2011), "Применение метода гомотопического анализа к электрогидродинамическому потоку", Commun. Нелинейный. Sci. Нумер. Simulat., 16 (7): 2730–2736, Bibcode:2011CNSNS..16.2730M, Дои:10.1016 / j.cns.2010.10.004
  18. ^ Нассар, Кристофер Дж .; Ревелли, Джозеф Ф .; Боуман, Роберт Дж. (2011), "Применение метода гомотопического анализа к уравнению Пуассона – Больцмана для полупроводниковых устройств", Коммунальное нелинейное научное численное моделирование, 16 (6): 2501–2512, Bibcode:2011CNSNS..16.2501N, Дои:10.1016 / j.cns.2010.09.015
  19. ^ Таджаддодианфар, Фарид (2017). «Нелинейная динамика резонаторов MEMS / NEMS: аналитическое решение методом гомотопического анализа». Микросистемные технологии. 23 (6): 1913–1926. Дои:10.1007 / s00542-016-2947-7.
  20. ^ Таджаддодианфар, Фарид (март 2015 г.). «О динамике бистабильных микро / нанорезонаторов: аналитическое решение и нелинейное поведение». Коммуникации в нелинейной науке и численном моделировании. 20 (3): 1078–1089. Bibcode:2015CNSNS..20.1078T. Дои:10.1016 / j.cnsns.2014.06.048.
  21. ^ Таджаддодианфар, Фарид (январь 2016 г.). «Прогнозирование хаоса в дуговых микро-нанорезонаторах с электростатическим возбуждением: аналитический подход». Коммуникации в нелинейной науке и численном моделировании. 30 (1–3): 182–195. Дои:10.1016 / j.cns.2015.06.013.

внешняя ссылка