Skema beda hingga untuk "persamaan gelombang", metode karakteristik

10

Pertimbangkan masalah berikut mana istilah pemaksaan dapat bergantung pada u , v (lihat Edit 1 di bawah untuk formulasi), dan W dan turunan pertamanya. Ini adalah persamaan gelombang 1 + 1 dimensi. Kami memiliki data awal yang ditentukan pada { u + v = 0 } .

Wuv=F
u,vW{u+v=0}

Saya tertarik pada solusi di dalam domain ketergantungan dari suatu interval dan sedang mempertimbangkan skema beda hingga berikut.

{u+v=0,u[uM,uM]}
  • Tujuannya adalah untuk berkembang oleh W u ( u , v + 1 ) - W u ( u , v ) = F ( u , v ) dan juga W v ( u + 1 , v ) - W v ( u , v ) = F ( u , v ) . Skema ini terintegrasi dalam arti bahwa WWuWu(u,v+1)Wu(u,v)=F(u,v)Wv(u+1,v)Wv(u,v)=F(u,v)
    W(kamu,v)+Wkamu(kamu,v)+Wv(kamu+1,v)=W(kamu+1,v+1)=W(kamu,v)+Wv(kamu,v)+Wkamu(kamu,v+1)
    jadi saya dapat secara konsisten menghitung dari data awal dengan mengintegrasikan ke atas; maka saya hanya benar-benar perlu melihat persamaan evolusi untuk W v dan W u .WWvWkamu
  • Untuk data awal, kita memerlukan kondisi kompatibilitas . Yang menunjukkan bahwa saya dapat menghitung data awal dengan menggunakan maju (dalam u ) beda hingga dari W pada waktu awal dengan nilai-nilai yang diberikan W tWkamu(kamu,v)-Wv(kamu+1,v-1)=W(kamu+1,v-1)-W(kamu,v)kamuWWtpada titik setengah bilangan bulat .(kamu+0,5,v-0,5)

Pertanyaan :

  1. Apakah ini skema yang terkenal? Secara khusus, di mana saya dapat menemukan analisis skema ini?
  2. Adakah hal yang jelas harus saya perhatikan?

Latar Belakang : Berpura-pura tidak tahu apa-apa (yang mungkin benar, karena saya seorang ahli matematika murni yang mencoba mempelajari sedikit mesin perhitungan).


Edit 1 : Hanya untuk memperjelas (untuk mengatasi beberapa komentar): persamaan di t koordinat akan W t t - W x x = F dan u dan v batal koordinat yang diberikan oleh (hingga beberapa faktor renormalising dari 2) u = t + x dan v = t - x . Jadi data awal di { u + v = 0 } sebenarnya di { t = 0 } .x t

Wtt-Wxx=F
kamuvkamu=t+xv=t-x{kamu+v=0}{t=0}

Jadi alih-alih sebuah mesh disesuaikan dengan saya menganggap mesh disesuaikan dengan ( u , v ) yang ¨ diproteksi 45 derajat ¨. Dibandingkan dengan ( t , x ) di mana t , x mengambil nilai integer, orang dapat menganggap u , v mesh memiliki titik tambahan di mana keduanya (tetapi tidak hanya satu) t dan x mengambil nilai setengah bilangan bulat.(t,x)(kamu,v)(t,x)t,xkamu,vtx

Willie Wong
sumber
Saya agak bingung dengan subskrip Anda, tetapi bagi saya ini sepertinya semacam rumusan domain waktu hingga-perbedaan . . . mungkin dengan formulasi mesh terhuyung-huyung (setengah indeks?).
meawoppl
1
@meawoppl: Dia hanya memanggil variabelnya bukan x , t seperti yang biasa dilakukan. (Dalam formulasi u , v yang biasa , mereka juga dirotasi oleh 45 dalam bidang ruang-waktu melawan x , t , tapi itu masalah yang terpisah.)kamu,vx,tkamu,v45x,t
Wolfgang Bangerth
Saya telah mengedit untuk mengklarifikasi (penjelasan Wolfgang Bangerth adalah apa yang ada dalam pikiran saya).
Willie Wong

Jawaban:

6

Pasti ada literatur tentang skema seperti ini. Dua kata kunci tersebut

  • Metode karakteristik yang dimodifikasi
  • Skema semi-Lagrangian

Setelah 20 menit googling: beberapa makalah yang mungkin penting adalah http://dx.doi.org/10.1137/0719063 dan http://dx.doi.org/10.1137/0728024 (cari dari sana). Itu mungkin bukan referensi terbaik di luar sana, tetapi itu harus menjadi titik awal untuk membawa Anda ke literatur yang tepat.

Saya menganggap ini sebagai metode garis yang diputar dengan pemisahan dimensional. Agaknya Anda sangat menyadari kesetaraan persamaan dan bentuk biasa dari gelombang persamaan di bawah transformasi u = t + x , v = t - x . Bagi saya, berguna untuk memikirkan skema Anda dalam bentuk persamaan gelombang tradisional ini. Apa yang dilakukan skema adalah mengintegrasikan pertama di sepanjang satu set karakteristik, kemudian di sepanjang yang lain. Integrasi dilakukan dengan menggunakan pemisahan dimensional dan metode Euler

Wtt-Wxx=F
kamu=t+x,    v=t-x.
, keduanya merupakan urutan pertama yang akurat.

Tentu saja, karena Anda mengintegrasikan sepanjang karakteristik, skema Anda akan tepat dalam kasus . Artinya, kesalahan numerik dalam skema Anda hanya akan disebabkan oleh integrasi numerik F (ini mungkin jelas, tetapi mungkin berguna untuk menunjukkan kepada mereka yang terbiasa dengan metode numerik yang lebih tradisional). Selanjutnya, skema Anda stabil tanpa syarat untuk kasus F = 0 . Tidak ada lagi yang bisa dikatakan tentang stabilitas tanpa mengetahui beberapa sifat dari F . Secara umum, skema akan stabil hanya di bawah beberapa pembatasan ukuran langkah hingga (karena metode Euler eksplisit). Jika Jacobian dari FF=0FF=0FF memiliki nilai eigen murni imajiner, skema ini tidak akan stabil.

Pendekatan diskritisasi umum untuk mengurangi PDE ke sistem ODE (seperti dalam metode Anda) dikenal sebagai metode garis. Seperti halnya metode diskritisasi garis apa pun, Anda dapat meningkatkan urutan keakuratan dengan menggunakan pemecah ODE tingkat tinggi dan Anda dapat meningkatkan stabilitas dengan menggunakan pemecah ODE implisit yang sesuai (dengan peningkatan biaya komputasi per langkah yang menyertainya).

David Ketcheson
sumber
"tetapi Google akan membantu Anda lebih banyak" Sebenarnya itu adalah salah satu masalah besar. Saya tidak yakin untuk apa Google (saya menduga literatur numerik dapat menggunakan beberapa istilah berbeda dari literatur murni). Jika Anda dapat menyarankan beberapa kata kunci yang harus saya cari, saya akan berterima kasih. ("Metode garis", misalnya, mengarahkan saya ke banyak informasi [mungkin bahkan sedikit bagi saya untuk dapat menyaring melalui :-)].)
Willie Wong
@ WillieWong - Salah satu referensi untuk persamaan hiperbolik yang biasa kita kutip adalah Metode Volume Hingga LeVeque untuk Masalah Hiperbolik . Saya tidak yakin apakah ini referensi yang tepat untuk Anda mulai, tetapi setidaknya akan memberi Anda pengantar istilah dan teknik di lapangan.
Aron Ahmadia
Oke, saya menambahkan beberapa kata kunci dan referensi. Saya harap mereka membantu.
David Ketcheson
Terima kasih banyak untuk referensi! Itu membuat saya awal yang baik.
Willie Wong
5

Berawal dari tempat David Ketcheson meninggalkan saya dalam jawabannya, sedikit lebih banyak pencarian mengungkapkan beberapa catatan sejarah.

Skema yang saya uraikan di atas dianggap sudah kembali pada tahun 1900 oleh J. Massau, di Mémoire sur l'intégration graphique des équations aux dérivées partielles . Karya ini diterbitkan ulang pada tahun 1952 oleh G. Delporte, Mons.

Analisis modern pertama (walaupun singkat) tentang konvergensi dan semacamnya diberikan oleh Courant, Friedrichs, dan Lewy dalam makalah klasik mereka tahun 1928 di Math. Ann.

Willie Wong
sumber
Wow, saya tidak percaya saya tidak menyadari bahwa ini ada di koran CFL ...
David Ketcheson