ODE Kaku dan Stabilitas

Persamaan diferensial kaku mengandung proses yang berkembang pada skala waktu yang sangat terpisah, sehingga metode eksplisit dipaksa untuk mengambil langkah-langkah yang sangat kecil secara tidak praktis untuk stabilitas; solusi efisiennya memerlukan metode implisit dengan sifat stabilitas yang kuat.

Definition

Persamaan diferensial disebut kaku ketika ia mengakui komponen solusi yang meluruh pada skala waktu yang sangat berbeda, sehingga stabilitas numerik daripada akurasi mendikte ukuran langkah; teori stabilitas menganalisis metode mana yang dapat mengambil langkah besar tanpa pertumbuhan kesalahan.

Scope

Topik ini mencakup fenomena dan definisi informal kekakuan, persamaan uji linier dan wilayah stabilitas absolut, konsep stabilitas-A, stabilitas-A(alpha), dan stabilitas-L, mengapa metode eksplisit gagal pada masalah kaku, dan metode implisit — Runge-Kutta implisit dan rumus diferensiasi mundur — yang menyelesaikannya.

Core questions

• Apa yang membuat suatu masalah kaku, dan mengapa hal itu mengalahkan metode eksplisit?
• Bagaimana wilayah stabilitas absolut didefinisikan melalui persamaan uji linier?
• Apa yang dibutuhkan oleh stabilitas-A dan stabilitas-L, dan mengapa keduanya penting untuk masalah kaku?
• Metode mana yang memberikan stabilitas yang dibutuhkan untuk sistem kaku dan diferensial-aljabar?

Key theories

Stabilitas absolut dan persamaan uji

Menerapkan suatu metode pada persamaan uji linier skalar menghasilkan faktor amplifikasi; himpunan produk ukuran langkah-kali-nilai eigen di mana faktor ini memiliki magnitudo paling banyak satu adalah wilayah stabilitas absolut metode tersebut, yang harus mengandung nilai eigen kaku masalah untuk memungkinkan langkah besar.

Stabilitas-A dan Stabilitas-L

Suatu metode bersifat stabil-A jika wilayah stabilitasnya mencakup seluruh bidang kiri, sehingga stabil untuk semua mode peluruhan terlepas dari ukuran langkah, dan stabil-L jika ia juga meredam mode yang sangat kaku sepenuhnya; sifat-sifat ini memilih metode implisit yang cocok untuk masalah kaku.

Mechanisms

Pada masalah kaku, mode peluruhan tercepat memiliki nilai eigen negatif yang besar; wilayah stabilitas terbatas metode eksplisit memaksa ukuran langkah untuk menyelesaikan mode tersebut bahkan jauh setelah mode tersebut secara fisik telah mati, membuat perhitungan menjadi sangat lambat. Metode implisit seperti metode Euler mundur, skema Runge-Kutta implisit, dan rumus diferensiasi mundur memiliki wilayah stabilitas yang mencakup bidang kiri (atau sebagian besar darinya), sehingga mereka tetap stabil pada langkah besar dan memungkinkan ukuran langkah dipilih hanya berdasarkan akurasi. Setiap langkah kemudian memerlukan penyelesaian sistem aljabar (umumnya nonlinier), biasanya dengan iterasi Newton menggunakan Jacobian.

Clinical relevance

Kekakuan tersebar luas dalam jaringan reaksi kimia, pembakaran, sirkuit listrik, sistem kontrol, dan diskretisasi metode garis dari persamaan diferensial parsial parabolik; mengenali kekakuan dan memilih pemecah implisit yang stabil secara tepat sangat penting untuk mendapatkan hasil dalam waktu yang layak, dan sebagian besar perangkat lunak ODE produksi mencakup deteksi kekakuan otomatis dan pengalihan.

History

Gagasan kekakuan diidentifikasi oleh Curtiss dan Hirschfelder pada tahun 1952, dan teori stabilitas pendukung — stabilitas-A dan batasan orde — dikembangkan oleh Dahlquist; kode rumus diferensiasi mundur Gear dan kemudian metode Runge-Kutta implisit orde tinggi membentuk perangkat praktis untuk masalah kaku dan diferensial-aljabar.

Key figures

• Germund Dahlquist
• C. William Gear
• Ernst Hairer
• Gerhard Wanner

Seminal works

• hairer1996
• iserles2008

Frequently asked questions

Apa sebenarnya yang membuat ODE menjadi kaku?

Kekakuan muncul ketika sistem memiliki komponen yang meluruh jauh lebih cepat daripada evolusi solusi yang diminati. Tidak ada definisi tunggal yang tajam, tetapi ciri praktisnya adalah metode eksplisit dipaksa untuk menggunakan langkah-langkah yang sangat kecil untuk stabilitas bahkan ketika akurasi akan memungkinkan langkah-langkah besar.

Mengapa masalah kaku memerlukan metode implisit?

Metode implisit dapat memiliki wilayah stabilitas yang mencakup seluruh bidang kiri (stabilitas-A), sehingga mereka tetap stabil pada ukuran langkah besar untuk mode yang meluruh dengan cepat. Metode eksplisit memiliki wilayah stabilitas terbatas, yang memaksa langkah-langkah kecil dan membuatnya tidak praktis untuk masalah kaku.