Чи дотримується принцип максимального / мінімального рівняння тепла за допомогою дискретизації Кранка-Ніколсона?


10

Я використовую схему кінцевих різниць Кранка-Ніколсона для вирішення 1D рівняння тепла. Мені цікаво, чи справедливий принцип максимального / мінімального рівняння тепла (тобто, що максимум / мінімум виникає при початковій умові або на межах) також для дискретного рішення.

Це, мабуть, має на увазі той факт, що Кранк-Нікольсон є стійкою та конвергентною схемою. Але здається, що ви могли б довести це безпосередньо через аргумент лінійної алгебри, використовуючи матриці, створені з трафарету Кранка-Ніколсона.

Я вдячний за будь-які вказівки до літератури на це. Дякую.


Привіт foobarbaz, і ласкаво просимо до scicomp! Я вважаю, що проблема, яку ви вирішуєте, не має джерела, правда?
Павло

Відповіді:


8

Максимальний принцип для Кранка-Нікольсона буде дотримуватися, якщо для кроку часуkта відстані сіткиh. Загалом, ми можемо вважатиθ-схему виду un+1=un+μ

μkh21
khθ деA- стандартна матриця Лаплаціа, а0θ1. Якщоμ(1-2θ)1
un+1=un+μ2((1θ)Aun+θAun+1)
A0θ1 , тоді схема стабільна. (Це легко показати методами Фур'є.) Однак, більш сильний критерій, щоμ(1-θ)1μ(12θ)12 потрібно для принципу максимуму в цілому.μ(1θ)12

Для доказів див. Числові рішення часткових диференціальних рівнянь К. В. Мортона . Зокрема, подивіться на розділи 2.10 та 2.11 та теорему 2.2.


Також є хороший спосіб побачити, що принцип Максиму взагалі не буде дотримуватися для Кранка-Ніколсона без обмеження на .μ

[0,1]uikkiu0k=u2k=0k

(1μ2(2))u1n+1=(1+μ2(2))u1n,
u1n+1=(1μ1+μ)u1n.

u10=1

u1n=(1μ1+μ)n,
u1n1u1n<0nμ1μ1μ

У відповідь на прохання foobarbaz я додав ескіз доказу.

(1+2θμ)ujn+1=θμ(uj1n+1+uj+1n+1)+(1θ)μ(uj1n+uj+1n)+[12(1θ)μ]ujn

μ(1θ)12

ujn+1uj1n+1uj+1n+1uj1nuj+1nujnujn+1ujn+1

(1+2θμ)ujn+1>θμ(uj1n+1+uj+1n+1)+(1θ)μ(uj1n+uj+1n)+[12(1θ)μ]ujn=(1+2θμ)ujn+1

ujn+1u


Дякую! Чи трапляється вам знати про іншу довідку, окрім Мортона? Я не можу отримати доступ до цих розділів або теореми в попередньому перегляді книг Google. Я хотів би зрозуміти доказ.
foobarbaz

@foobarbaz У мене немає ще однієї посилання, але я додав контур доказу. Дайте мені знати, чи можу я зробити це зрозумілішим.
Бен

0

Стійкість означає, що збурення залишаються обмеженими у часі. Це не означає, що принцип максимального задовольняється на дискретному рівні, це вже інше питання. Задоволення принципу дискретного максимуму достатньо, але не потрібно для стабільності.

Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.