Jika kita menggunakan Metode Garis untuk diskritisasi (waktu dan ruang terpisah diskritisasi) PDE hiperbolik yang kita peroleh setelah diskritisasi spasial dengan metode numerik favorit kita (fx. Metode Volume Hingga) dalam praktiknya pemecah ODE mana yang kita gunakan untuk diskretisasi temporal (TVD / SSP / dll)?
Beberapa informasi tambahan ditambahkan: Masalah akurasi dapat menjadi masalah untuk masalah yang tidak mulus. Telah diketahui bahwa PDE hiperbolik nonlinear dapat menyebabkan guncangan dalam waktu yang terbatas meskipun solusi awalnya halus di mana akurasi kasus dapat menurun ke urutan pertama untuk metode tingkat tinggi.
Analisis Stabilitas ODE biasanya dilakukan berdasarkan linierisasi untuk mendapatkan sistem ODE semi-diskrit linear dari bentuk q_t = J q (dengan vektor perturbasi qa), di mana nilai eigen J harus ditingkatkan dalam wilayah stabilitas absolut waktu yang dipilih- metode melangkah. Strategi alternatif adalah dengan menggunakan pseudospectra atau mungkin metode energi untuk analisis stabilitas.
Saya memahami bahwa motivasi untuk metode TVD / SSP adalah untuk menghindari osilasi palsu yang disebabkan oleh metode yang melangkah waktu yang dapat mengakibatkan perilaku yang tidak fisik. Pertanyaannya adalah apakah pengalaman menunjukkan jenis metode penginjilan waktu ini lebih unggul dibandingkan, misalnya, kuda kerja klasik sebagai Metode Runge-Kutta eksplisit atau yang lainnya. Jelas, mereka harus memiliki sifat yang lebih baik untuk kelas masalah di mana solusi dapat menunjukkan guncangan. Karena itu orang dapat berargumen bahwa kita hanya boleh menggunakan metode-metode jenis ini untuk integrasi waktu.
sumber
Ya itu penting. Dua hal biasa yang harus diperhatikan:
Ketepatan. Beberapa skema ODE lebih akurat daripada yang lain, urutan lebih tinggi, dan sebagainya. Aturan praktisnya adalah memilih metode dengan urutan akurasi yang mirip dengan diskritisasi spasial Anda.
Stabilitas. Untuk masalah hiperbolik Anda mengharapkan operator memiliki nilai eigen imajiner murni, jadi Anda menginginkan pemecah ODE yang mencakup beberapa bagian dari akses imajiner dalam domain stabilitasnya. Lihat misalnya Lampiran G dalam Fornberg, Panduan Praktis untuk Metode Pseudospectral.
Dengan persamaan hiperbolik, beberapa orang ingin memastikan bahwa solusi mereka selalu positif, sehingga ada berbagai jenis filter dan trik untuk memastikan ini. Tapi saya hampir tidak tahu apa-apa tentang ini.
Saya jauh dari ahli, tetapi saya pikir saya akan mencoba menjawab karena pertanyaannya sudah ada di sini untuk sementara waktu.
sumber