В математике уравнение в частных производных ( УЧП ) — это уравнение, которое вычисляет функцию между различными частными производными функции многих переменных .
Функцию часто считают «неизвестной», которую нужно решить, аналогично тому, как x считается неизвестным числом, которое нужно решить в алгебраическом уравнении, таком как x 2 − 3 x + 2 = 0 . Однако записать явные формулы для решений уравнений в частных производных обычно невозможно. Соответственно, существует огромное количество современных математических и научных исследований методов численной аппроксимации решений некоторых уравнений в частных производных с использованием компьютеров. Уравнения в частных производных занимают также большой сектор чисто математических исследований , в которых обычные вопросы заключаются, вообще говоря, в выявлении общих качественных особенностей решений различных уравнений в частных производных, таких как существование, единственность, регулярность и устойчивость. [1] Среди многих открытых вопросов - существование и гладкость решений уравнений Навье-Стокса , названных одной из проблем Премии тысячелетия в 2000 году.
Уравнения с частными производными повсеместно распространены в математически ориентированных научных областях, таких как физика и техника . Например, они лежат в основе современного научного понимания звука , тепла , диффузии , электростатики , электродинамики , термодинамики , гидродинамики , упругости , общей теории относительности и квантовой механики ( уравнение Шредингера , уравнение Паули и т. д.). Они также возникают из многих чисто математических соображений, таких как дифференциальная геометрия и вариационное исчисление ; среди других заметных применений они являются основным инструментом доказательства гипотезы Пуанкаре на основе геометрической топологии .
Частично из-за такого разнообразия источников существует широкий спектр различных типов уравнений в частных производных, и были разработаны методы работы со многими из возникающих отдельных уравнений. По сути, обычно признается, что не существует «общей теории» уравнений в частных производных, а специальные знания в некоторой степени разделены между несколькими существенно различными подобластями. [2]
Обыкновенные дифференциальные уравнения образуют подкласс уравнений в частных производных, соответствующих функциям одной переменной. Стохастические уравнения в частных производных и нелокальные уравнения по состоянию на 2020 год являются особенно широко изученными расширениями понятия «PDE». Более классические темы, по которым все еще ведется много активных исследований, включают эллиптические и параболические уравнения в частных производных, механику жидкости , уравнения Больцмана и дисперсионные уравнения в частных производных. [3]
Функция u ( x , y , z ) трех переменных является « гармонической » или «решением уравнения Лапласа », если она удовлетворяет условию
Природу этой неудачи можно увидеть более конкретно в случае следующего УЧП: для функции v ( x , y ) двух переменных рассмотрим уравнение
Характер этого выбора варьируется от PDE к PDE. Чтобы понять это для любого данного уравнения, важными организационными принципами обычно являются теоремы существования и единственности . Во многих вводных учебниках роль теорем существования и единственности ОДУ может быть несколько неясной; половина существования обычно не нужна, поскольку можно напрямую проверить любую предлагаемую формулу решения, в то время как половина уникальности часто присутствует только в фоновом режиме, чтобы гарантировать, что предлагаемая формула решения является максимально общей. Напротив, для PDE теоремы существования и единственности часто являются единственным средством, с помощью которого можно ориентироваться в множестве различных решений. По этой причине они также имеют основополагающее значение при проведении чисто численного моделирования, поскольку необходимо понимать, какие данные должен задать пользователь, а какие оставить для расчета компьютеру.
Чтобы обсудить такие теоремы существования и единственности, необходимо точно определить область определения «неизвестной функции». В противном случае, говоря только в терминах типа «функция двух переменных», невозможно осмысленно сформулировать результаты. То есть область определения неизвестной функции следует рассматривать как часть структуры самого УЧП.
Ниже приведены два классических примера таких теорем существования и единственности. Несмотря на то, что два рассматриваемых УЧП очень похожи, существует разительная разница в поведении: для первого УЧП имеется свободное предписание одной функции, тогда как для второго УЧП имеется свободное предписание двух функций.
Возможны и другие явления. Например, следующее УЧП , естественным образом возникающее в области дифференциальной геометрии , иллюстрирует пример, где существует простая и вполне явная формула решения, но со свободным выбором только трёх чисел и даже не одной функции.
В отличие от предыдущих примеров, это УЧП является нелинейным из-за квадратных корней и квадратов. Линейное УЧП — это такое уравнение, что, если оно однородно, сумма любых двух решений также является решением, и любое постоянное кратное любому решению также является решением .
Правильность относится к общему схематическому пакету информации о PDE. Чтобы сказать, что PDE корректен, необходимо иметь:
Из-за необходимости применимости к нескольким различным PDE это несколько расплывчато. Требование «непрерывности», в частности, неоднозначно, поскольку обычно существует множество неэквивалентных средств, с помощью которых его можно строго определить. Однако несколько необычно изучать УЧП без указания того, каким образом оно корректно.
Энергетический метод представляет собой математическую процедуру, которую можно использовать для проверки корректности начально-краевых задач (ИБВП). [4] В следующем примере энергетический метод используется для решения, где и какие граничные условия следует наложить, чтобы полученный IBVP был корректным. Рассмотрим одномерное гиперболическое УЧП, заданное формулой
где – константа, – неизвестная функция с начальным условием . Умножение и интегрирование по области дает
Используя это
Здесь обозначает стандартную норму . Для корректности мы требуем, чтобы энергия решения не возрастала, т. е. была равна , что достигается заданием if и at if . Это соответствует лишь наложению граничных условий на притоке. Корректность допускает рост данных (начальных и граничных), и, таким образом, достаточно показать, что это справедливо, когда все данные установлены на ноль.
Теорема Коши–Ковалевского для задач Коши с начальными значениями по существу утверждает, что если все члены уравнения в частных производных состоят из аналитических функций и выполняется определенное условие трансверсальности (гиперплоскость или, в более общем смысле, гиперповерхность, на которой задаются начальные данные, должна быть нехарактеристичны относительно оператора в частных производных), то на некоторых областях обязательно существуют решения, которые также являются аналитическими функциями. Это фундаментальный результат в изучении аналитических уравнений в частных производных. Удивительно, но теорема не справедлива в случае гладких функций; пример , открытый Гансом Леви в 1957 году, состоит из линейного уравнения в частных производных, коэффициенты которого гладкие (т. е. имеют производные всех порядков), но не аналитические, для которого не существует решения. Таким образом, теорема Коши-Ковалевского обязательно ограничивается областью применения аналитическими функциями.
При написании УЧП принято обозначать частные производные с помощью индексов. Например:
Греческая буква Δ обозначает оператор Лапласа ; если u — функция n переменных, то
УЧП называется линейным, если оно линейно относительно неизвестного и его производных. Например, для функции u от x и y линейное УЧП второго порядка имеет вид
Тремя основными типами нелинейных УЧП являются полулинейные УЧП, квазилинейные УЧП и полностью нелинейные УЧП.
Ближайшими к линейным УЧП являются полулинейные УЧП, в которых только производные высшего порядка выступают в виде линейных членов с коэффициентами, которые являются функциями независимых переменных. Младшие производные и неизвестная функция могут появляться произвольно. Например, общее полулинейное УЧП второго порядка с двумя переменными имеет вид
В квазилинейном УЧП производные высшего порядка также появляются только как линейные члены, но с коэффициентами, возможно, функциями неизвестных и производных низшего порядка:
УЧП без каких-либо свойств линейности называется полностью нелинейным и обладает нелинейностью в одной или нескольких производных высшего порядка. Примером может служить уравнение Монжа–Ампера , возникающее в дифференциальной геометрии . [5]
Эллиптические , параболические и гиперболические уравнения в частных производных второго порядка широко изучаются с начала двадцатого века. Однако существует много других важных типов УЧП, таких как нелинейное уравнение Кортевега – де Фриза третьего порядка . Существуют также гибриды, такие как уравнение Эйлера-Трикоми , которые варьируются от эллиптического до гиперболического для разных областей области. Существуют также важные расширения этих базовых типов для PDE более высокого порядка, но такие знания более специализированы.
Эллиптическая/параболическая/гиперболическая классификация дает представление о соответствующих начальных и граничных условиях , а также о гладкости решений. Полагая u xy = u yx , общее линейное УЧП второго порядка с двумя независимыми переменными имеет вид
Точнее, замена ∂x на X , и аналогично для других переменных (формально это делается преобразованием Фурье ), преобразует УЧП с постоянным коэффициентом в многочлен той же степени , с членами высшей степени (однородный многочлен , здесь квадратичная форма ), являющаяся наиболее значимой для классификации.
Точно так же, как конические сечения и квадратичные формы классифицируются на параболические, гиперболические и эллиптические на основе дискриминанта B 2 − 4 AC , то же самое можно сделать и для УЧП второго порядка в данной точке. Однако дискриминант в УЧП определяется как B 2 − AC , поскольку принято считать, что член xy равен 2 B , а не B ; формально дискриминант (соответствующей квадратичной формы) равен (2 B ) 2 − 4 AC = 4( B 2 − AC ) с коэффициентом 4, опущенным для простоты.
Если имеется n независимых переменных x 1 , x 2 , …, x n , общее линейное уравнение в частных производных второго порядка имеет вид
Классификация зависит от подписи собственных значений матрицы коэффициентов a i , j .
Теория эллиптических, параболических и гиперболических уравнений изучалась на протяжении веков, в основном сосредоточенная вокруг или основанная на стандартных примерах уравнения Лапласа , уравнения теплопроводности и волнового уравнения .
Классификацию уравнений в частных производных можно распространить на системы уравнений первого порядка, где неизвестное u теперь является вектором с m компонентами, а матрицы коэффициентов A ν представляют собой матрицы размером m на m для ν = 1, 2, …, n. . Уравнение в частных производных принимает вид
Геометрическая интерпретация этого условия следующая: если данные для u заданы на поверхности S , то можно определить нормальную производную u на S из дифференциального уравнения. Если данные о S и дифференциальное уравнение определяют нормальную производную от u на S , то S нехарактеристична. Если данные о S и дифференциальное уравнение не определяют нормальную производную от u на S , то поверхность является характеристической , а дифференциальное уравнение ограничивает данные о S : дифференциальное уравнение является внутренним для S.
Линейные УЧП можно свести к системам обыкновенных дифференциальных уравнений с помощью важного метода разделения переменных. Этот метод основан на особенности решений дифференциальных уравнений: если можно найти какое-либо решение, которое решает уравнение и удовлетворяет граничным условиям, то это решение (это также относится и к ОДУ). В качестве анзаца мы предполагаем , что зависимость решения от параметров пространства и времени можно записать как произведение членов, каждый из которых зависит от одного параметра, а затем посмотрим, можно ли это сделать для решения проблемы. [8]
В методе разделения переменных УЧП сводится к УЧП с меньшим количеством переменных, которое представляет собой обыкновенное дифференциальное уравнение, если с одной переменной - их, в свою очередь, легче решить.
Это возможно для простых УЧП, которые называются разделимыми уравнениями в частных производных , а область определения обычно представляет собой прямоугольник (произведение интервалов). Разделимые УЧП соответствуют диагональным матрицам : если рассматривать «значение фиксированного x » как координату, каждую координату можно понимать отдельно.
Это обобщает метод характеристик , а также используется в интегральных преобразованиях .
В особых случаях можно найти характеристические кривые, на которых уравнение сводится к ОДУ — изменение координат в области для выпрямления этих кривых позволяет разделить переменные и называется методом характеристик .
В более общем плане можно найти характерные поверхности. Для решения уравнения в частных производных второго порядка см. метод Шарпита .
Интегральное преобразование может преобразовать УЧП в более простое, в частности, в разделимое УЧП. Это соответствует диагонализации оператора.
Важным примером этого является анализ Фурье , который диагонализует уравнение теплопроводности, используя собственный базис синусоидальных волн.
Если область конечная или периодическая, подходит бесконечная сумма решений, таких как ряд Фурье , но для бесконечных областей обычно требуется интеграл решений, такой как интеграл Фурье . Решение для точечного источника приведенного выше уравнения теплопроводности является примером использования интеграла Фурье.
Часто УЧП можно привести к более простой форме с известным решением путем подходящей замены переменных . Например, уравнение Блэка – Шоулза
Неоднородные уравнения [ необходимы пояснения ] часто можно решить (для УЧП с постоянными коэффициентами всегда решаются), найдя фундаментальное решение (решение для точечного источника), а затем выполнив свертку с граничными условиями, чтобы получить решение.
В обработке сигналов это аналогично пониманию фильтра по его импульсной характеристике .
Принцип суперпозиции применим к любой линейной системе, включая линейные системы УЧП. Распространенной визуализацией этой концепции является взаимодействие двух волн в фазе, которые объединяются, что приводит к большей амплитуде, например sin x + sin x = 2 sin x . Тот же принцип можно наблюдать в PDE, где решения могут быть действительными или сложными и аддитивными. Если u 1 и u 2 являются решениями линейного УЧП в некотором функциональном пространстве R , то u = c 1 u 1 + c 2 u 2 с любыми константами c 1 и c 2 также являются решением этого УЧП в том же функциональном пространстве.
Общеприменимых методов решения нелинейных уравнений в уравнениях не существует. Тем не менее, результаты о существовании и единственности (такие как теорема Коши-Ковалевского ) часто возможны, как и доказательства важных качественных и количественных свойств решений (получение этих результатов является основной частью анализа ). Вычислительное решение нелинейных УЧП, метод разделенных шагов , существует для конкретных уравнений, таких как нелинейное уравнение Шредингера .
Тем не менее, некоторые методы можно использовать для нескольких типов уравнений. H - принцип является наиболее мощным методом решения недоопределенных уравнений. Теория Рикье-Жане является эффективным методом получения информации о многих аналитических переопределенных системах.
Метод характеристик может быть использован в некоторых весьма частных случаях для решения нелинейных уравнений в частных производных. [10]
В некоторых случаях УЧП можно решить с помощью анализа возмущений , при котором решением считается поправка к уравнению с известным решением. Альтернативой являются методы численного анализа , от простых схем конечных разностей до более зрелых многосеточных методов и методов конечных элементов . Многие интересные задачи в науке и технике решаются таким образом с помощью компьютеров , иногда высокопроизводительных суперкомпьютеров .
Начиная с 1870 года работы Софуса Ли поставили теорию дифференциальных уравнений на более удовлетворительную основу. Он показал, что теории интеграции старых математиков могут быть отнесены к общему источнику путем введения того, что сейчас называется группами Ли ; и что обыкновенные дифференциальные уравнения, допускающие одни и те же бесконечно малые преобразования, представляют сравнимые трудности при интегрировании. Он также остановился на теме трансформаций контакта .
Общий подход к решению УЧП использует свойство симметрии дифференциальных уравнений — непрерывные бесконечно малые преобразования решений в решения ( теория Ли ). Непрерывная теория групп , алгебры Ли и дифференциальная геометрия используются для понимания структуры линейных и нелинейных уравнений в частных производных, для создания интегрируемых уравнений, для нахождения их пар Лакса , операторов рекурсии, преобразования Беклунда и, наконец, для поиска точных аналитических решений УЧП.
Методы симметрии получили признание при изучении дифференциальных уравнений, возникающих в математике, физике, технике и многих других дисциплинах.
Метод разложения Адомиана , [11] искусственный метод малого параметра Ляпунова и его метод гомотопического возмущения — все это частные случаи более общего метода гомотопического анализа . [12] Это методы разложения в ряд, и, за исключением метода Ляпунова, они не зависят от малых физических параметров по сравнению с хорошо известной теорией возмущений , что придает этим методам большую гибкость и общность решения.
Тремя наиболее широко используемыми численными методами для решения PDE являются метод конечных элементов (FEM), методы конечных объемов (FVM) и методы конечных разностей (FDM), а также другие методы, называемые бессеточными методами , которые были созданы для решения задач, в которых вышеупомянутые методы ограничены. МКЭ занимает видное место среди этих методов, особенно его исключительно эффективная версия высшего порядка hp-FEM . Другие гибридные версии FEM и бессеточных методов включают обобщенный метод конечных элементов (GFEM), расширенный метод конечных элементов (XFEM), спектральный метод конечных элементов (SFEM), бессеточный метод конечных элементов , разрывный метод конечных элементов Галеркина (DGFEM), элементный метод. свободный метод Галеркина (EFGM), интерполяционный безэлементный метод Галёркина (IEFGM) и др.
Метод конечных элементов (МКЭ) (его практическое применение, часто известное как анализ конечных элементов (FEA)) представляет собой численный метод поиска приближенных решений уравнений в частных производных (PDE), а также интегральных уравнений. [13] [14] Подход к решению основан либо на полном исключении дифференциального уравнения (стационарные задачи), либо на преобразовании УЧП в аппроксимирующую систему обыкновенных дифференциальных уравнений, которые затем численно интегрируются с использованием стандартных методов, таких как метод Эйлера, Рунге-Кутта и др.
Методы конечных разностей — это численные методы аппроксимации решений дифференциальных уравнений с использованием конечно-разностных уравнений для аппроксимации производных.
Подобно методу конечных разностей или методу конечных элементов, значения вычисляются в дискретных местах сетчатой геометрии. «Конечный объем» относится к небольшому объему, окружающему каждую узловую точку сетки. В методе конечного объема поверхностные интегралы в уравнении в частных производных, содержащие член дивергенции, преобразуются в объемные интегралы с использованием теоремы о дивергенции . Эти члены затем оцениваются как потоки на поверхностях каждого конечного объема. Поскольку поток, входящий в данный объем, идентичен потоку, выходящему из соседнего объема, эти методы по своей конструкции сохраняют массу.
Нечеткое дифференциальное уравнение основано на нечеткой логике.
Метод нечеткого дифференциального включения использует нечеткую логику с преобразованием Лапласа для решения УЧП.
Уравнения в частных производных используются во многих областях: Астрономия , Космология , Квантовая механика , Теплопередача , Электромагнетизм , Гидродинамика , Упругость (физика) , Тензор упругости , Тензорный оператор Аналитическая геометрия , Искусственный интеллект , Глубокое обучение , Языковая модель
Некоторые распространенные PDE
Типы граничных условий
Различные темы
{{citation}}
: CS1 maint: multiple names: authors list (link).