Apakah prinsip maksimum / minimum dari persamaan panas dipertahankan oleh diskritisasi Crank-Nicolson?

10

Saya menggunakan skema beda hingga Crank-Nicolson untuk menyelesaikan persamaan panas 1D. Saya bertanya-tanya apakah prinsip maksimum / minimum dari persamaan panas (yaitu bahwa maksimum / minimum terjadi pada kondisi awal atau pada batas-batas) juga berlaku untuk solusi diskritisasi.

Ini mungkin tersirat oleh fakta bahwa Crank-Nicolson adalah skema yang stabil dan konvergen. Tetapi tampaknya Anda mungkin dapat membuktikan ini secara langsung melalui argumen aljabar linier menggunakan matriks yang dibuat dari stensil Crank-Nicolson.

Saya menghargai petunjuk apa pun untuk literatur tentang ini. Terima kasih.

foobarbaz
sumber
Hai foobarbaz, dan selamat datang di scicomp! Saya berasumsi Anda bahwa masalah yang Anda selesaikan tidak memiliki istilah sumber, benar?
Paul

Jawaban:

8

Prinsip maksimum untuk Crank-Nicolson akan berlaku jika untuk timestepkdan spasi gridh. Secara umum, kita dapat mempertimbangkanskemaθdari bentuk un+1=un+μ

μkh21
khθ manaAadalah matriks Laplacian standar dan0θ1. Jikaμ(1-2θ)1
kamun+1=kamun+μ2((1-θ)SEBUAHkamun+θSEBUAHkamun+1)
SEBUAH0θ1 , maka skemanya stabil. (Ini dapat dengan mudah ditunjukkan dengan teknik Fourier.) Namun, kriteria kuat bahwaμ(1-θ)1μ(1-2θ)12 diperlukan agar prinsip maksimum berlaku secara umum.μ(1-θ)12

Untuk bukti, lihat Solusi Numerik dari Persamaan Diferensial Parsial oleh KW Morton . Secara khusus, lihat Bagian 2.10 dan 2.11 dan Teorema 2.2.


Ada juga cara yang baik untuk melihat bahwa prinsip maksimum tidak akan berlaku secara umum untuk Crank-Nicolson tanpa kendala pada .μ

Pertimbangkan persamaan panas pada dengan diskritisasi yang mengandung 3 titik, termasuk batas. Biarkan u k i menunjukkan diskretisasi di timestep k dan titik grid i . Asumsikan batas Dirichlet, sehingga u k 0 = u k 2 = 0 untuk semua k . Kemudian Crank-Nicolson berkurang menjadi ( 1 - μ[0,1]kamusayakksayakamu0k=kamu2k=0k yang selanjutnya dapat direduksi menjadi u n + 1 1 =(1-μ

(1-μ2(-2))kamu1n+1=(1+μ2(-2))kamu1n,
kamu1n+1=(1-μ1+μ)kamu1n.

Jika kita mempertimbangkan kondisi awal , maka kita memiliki u n 1 = ( 1 - μkamu10=1

kamu1n=(1-μ1+μ)n,
kamu1n1kamu1n<0nμ1μ1μ

Menanggapi permintaan foobarbaz, saya telah menambahkan sketsa buktinya.

(1+2θμ)kamujn+1=θμ(kamuj-1n+1+kamuj+1n+1)+(1-θ)μ(kamuj-1n+kamuj+1n)+[1-2(1-θ)μ]kamujn

μ(1-θ)12

kamujn+1kamuj-1n+1kamuj+1n+1kamuj-1nkamuj+1nkamujnkamujn+1kamujn+1

(1+2θμ)kamujn+1>θμ(kamuj-1n+1+kamuj+1n+1)+(1-θ)μ(kamuj-1n+kamuj+1n)+[1-2(1-θ)μ]kamujn=(1+2θμ)kamujn+1

kamujn+1kamu

Ben
sumber
Terima kasih! Apakah Anda mengetahui referensi lain selain Morton? Saya tidak dapat mengakses Bagian atau Teorema itu di pratinjau buku Google. Saya ingin mengerti buktinya.
foobarbaz
@foobarbaz Saya tidak punya referensi lain yang berguna, tetapi saya menambahkan garis besar buktinya. Beri tahu saya jika saya bisa membuatnya lebih jelas.
Ben
0

Stabilitas berarti bahwa gangguan tetap dibatasi waktu. Itu tidak berarti bahwa prinsip maksimum terpenuhi pada tingkat diskrit, itu masalah yang berbeda. Memuaskan prinsip maksimum diskrit sudah cukup tetapi tidak perlu untuk stabilitas.

chris
sumber