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 } .
Saya tertarik pada solusi di dalam domain ketergantungan dari suatu interval dan sedang mempertimbangkan skema beda hingga berikut.
- 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 Wjadi 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 .
- 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 tpada titik setengah bilangan bulat .
Pertanyaan :
- Apakah ini skema yang terkenal? Secara khusus, di mana saya dapat menemukan analisis skema ini?
- 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 } .
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.
sumber
Jawaban:
Pasti ada literatur tentang skema seperti ini. Dua kata kunci tersebut
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
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= 0 F F= 0 F F 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).
sumber
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.
sumber