Какова цель использования интеграции по частям при получении слабой формы для дискретизации FEM?

24

При переходе от сильной формы PDE к форме FEM кажется, что всегда следует делать это, сначала указав вариационную форму. Для этого вы умножаете сильную форму на элемент в некотором (соболевском) пространстве и интегрируете по своему региону. Это я могу принять. Я не понимаю, почему нужно использовать формулу Грина (один или несколько раз).

Я в основном работал с уравнением Пуассона, так что если мы возьмем это (с однородными граничными условиями Дирихле) в качестве примера, то есть

2u=f,uΩu=0,uΩ

тогда утверждается, что правильный способ формирования вариационной формы

Ωfvdx=Ω2uvdx=ΩuvdxΩnuvds=Ωuvdx.

Но что мешает мне использовать выражение в первой строке, разве это не вариационная форма, которую можно использовать для получения формы FEM? Разве не соответствующий билинейной и линейных форм b(u,v)=(2u,v) и ? Проблема здесь в том, что если я использую линейные базисные функции (функции формы), то у меня будут проблемы, потому что моя матрица жесткости будет нулевой (не обратимой)? Но что, если я использую нелинейные функции формы? Я все еще должен использовать формулу Грина? Если мне не нужно: это целесообразно? Если нет, то есть ли у меня вариационная, но не слабая формулировка?l(v)=(f,v)

Теперь, допустим, у меня есть PDE с производными более высокого порядка, означает ли это, что существует много возможных вариационных форм, в зависимости от того, как я использую формулу Грина? И все они приводят к (различным) приближениям FEM?

Кристиан
источник
Связанный: scicomp.stackexchange.com/questions/667/…
Пол

Ответы:

18

Краткий ответ:

Нет, вам не нужно делать интеграцию для определенных FEM. Но в вашем случае вы должны это сделать.


Длинный ответ:

  • Допустим, - решение с конечными элементами. Если в качестве основы вы выберете кусочно-линейный многочлен, то, взяв за него Δ , вы получите распределение 1-го порядка (подумайте о производной по ступенчатой ​​функции Хевисайда), а интегрирование - Δ u hH - 1, умноженное на v, даст только имеет смысл, когда вы воспринимаете это как дуальную пару, а не как L 2 -иннер. Вы также не получите нулевую матрицу, теорема Рисса о представлении говорит, что в φ - Δ u hH 1 0 есть элементUчасΔ-ΔUчасЧАС-1vL2φ-ΔUчасЧАС01может характеризовать пару двойственности с помощью скалярного произведения в : - Δ у ч , v Н - 1 , Н 1 0 = Омф - Δ U чv скалярное произведение в  H 1 . Интегрирование по частям элемент за элементом для u h пролит свет на эту пару двойственности: для T элемент в этой триангуляции Ωu hЧАС1

    -ΔUчас,vЧАС-1,ЧАС01знак равноΩφ-ΔUчасvвнутренний продукт в ЧАС1,
    UчасT это говорит, что - Δ у ч должны включать межэлементный скачок потока в его двойственности паре представления, обратите вниманиеинтегрирование на границе каждого элемента также двойственность пара между H 1 / 2 и Н - 1 / 2 , Даже если вы используете квадратичный базис, который имеет неисчезающий Δ на каждом элементе, вы все равно не сможете записать ( Δ u , v ) как внутренний продукт из-за присутствия этого межэлементного скачка потока.
    ΩUчасvзнак равно-ΣT(TΔUчасv+TUчасNvdS),
    -ΔUчасЧАС1/2ЧАС-1/2Δ(ΔU,v)
  • WК,пWК,пЧАС1ЧАС2ЧАС2ΔUчас

  • {σзнак равно-U,σзнак равное,
    J(v)знак равно| |σ+U| |L2Ω2+| |σ-е| |L2Ω2,
    В соответствии с тем же духом, что и функционал Ритца-Галеркина, формулировка конечного элемента минимизации над функционалом в пространстве конечных элементов не требует интегрирования по частям.
Шухао Цао
источник
17

ЧАС2ЧАС2

Reid.Atcheson
источник
1
ЧАС2L2
1
То, что вы говорите, по сути правильно. Что касается PDE более высокого порядка, вам не обязательно использовать пространства с более высокой регулярностью, поскольку запись смешанной формулировки (см. Ответ Shuhao) может помочь. Вы также можете использовать другие методы, такие как штраф за прыжок, чтобы избежать этой трудности. Для классического ответа FEM, хотя, да, вам нужна более высокая регулярность.
Reid.Atcheson
2
Позвольте мне подчеркнуть важность симметрии. Если дифференциальный оператор самосопряженный, я ожидаю, что в итоге получится симметричная матрица. Без интеграции по частям это не будет иметь место.
Стефано М
1
Вычислительная выгода была моей главной мыслью при добавлении этого, но есть ли сильные теоретические преимущества симметрии (кроме более простых доказательств фактов, которые, вероятно, все еще сохраняются в эллиптическом случае, даже если дискретизация несимметрична)?
Reid.Atcheson
15

Прекрасные ответы уже есть на этой странице, но все еще есть (маленький) недостаток.

ОП спросил:

Теперь, допустим, у меня есть PDE с производными более высокого порядка, означает ли это, что существует много возможных вариационных форм, в зависимости от того, как я использую формулу Грина? И все они приводят к (различным) приближениям FEM?

Интегрирование по частям ( правильно ) важно при наличии граничных условий типа Неймана. На самом деле именно по IBP вы принимаете во внимание Neumann bc в своей вариационной формулировке. Форма Неймана до н.э. зависит от того, как вы интегрируете по частям, ср. этот ответ на интегрирование по частям в линейной упругости. Таким образом, даже для эллиптических уравнений в частных производных второго порядка интегрирование по частям должно выполняться заданным образом, чтобы восстановить вариационную формулировку, действительную для неймановских или смешанных граничных условий. (И это, конечно, независимо от того, что вы дискретизируете с помощью FEM).

В математической физике, где Нейман до н.э. имеет четко определенное значение (тепловой поток, напряжение ...), интеграция по частям важна для обеспечения правильной интерпретации результатов. Даже для однородных условий Дирихле и МКЭ это верно, поскольку, если мы используем метод множителей Лагранжа для наложения БК, множители становятся физическими величинами, такими как концентрированные потоки или силы.

Стефано М
источник