Katı Adi Diferansiyel Denklemler ve Kararlılık
Katı diferansiyel denklemler, geniş ölçüde ayrılmış zaman ölçeklerinde gelişen süreçleri içermektedir; bu nedenle açık yöntemler, kararlılık için pratik olmayan derecede küçük adımlar atmak zorunda kalmaktadır. Bunların etkin çözümü, güçlü kararlılık özelliklerine sahip örtük yöntemler gerektirmektedir.
Tanım
Bir diferansiyel denklem, çok farklı zaman ölçeklerinde bozunan çözüm bileşenlerini kabul ettiğinde katı olarak adlandırılmaktadır; böylece adım boyutunu doğruluktan ziyade sayısal kararlılık belirlemektedir. Kararlılık teorisi, hangi yöntemlerin hata büyümesi olmadan büyük adımlar atabileceğini analiz etmektedir.
Kapsam
Bu konu, katılık olgusunu ve gayri resmi tanımını, doğrusal test denklemini ve mutlak kararlılık bölgesini, A-kararlılığı, A(alfa)-kararlılığı ve L-kararlılığı kavramlarını, açık yöntemlerin katı problemlerde neden başarısız olduğunu ve bunları çözen örtük yöntemleri — örtük Runge-Kutta ve geri fark formüllerini — ele almaktadır.
Temel sorular
- Bir problemi katı yapan nedir ve neden açık yöntemleri başarısız kılar?
- Mutlak kararlılık bölgesi, doğrusal test denklemi aracılığıyla nasıl tanımlanır?
- A-kararlılığı ve L-kararlılığı ne gerektirir ve katı problemler için neden önemlidir?
- Katı ve diferansiyel-cebirsel sistemler için gerekli kararlılığı hangi yöntemler sağlar?
Temel kuramlar
- Mutlak kararlılık ve test denklemi
- Bir yöntemin skaler doğrusal test denklemine uygulanması bir amplifikasyon faktörü üretir; bu faktörün büyüklüğünün en fazla bir olduğu adım boyutu-özdeğer çarpımları kümesi, yöntemin mutlak kararlılık bölgesidir ve bu bölge, büyük adımlara izin vermek için problemin katı özdeğerlerini içermelidir.
- A-kararlılığı ve L-kararlılığı
- Bir yöntem, kararlılık bölgesi tüm sol yarı düzlemi içeriyorsa A-kararlıdır, böylece adım boyutundan bağımsız olarak tüm bozunan modlar için kararlıdır; ve ek olarak çok katı modları tamamen sönümlüyorsa L-kararlıdır. Bu özellikler, katı problemlere uygun örtük yöntemleri belirlemektedir.
Mekanizmalar
Katı bir problemde, en hızlı bozunan mod, büyük bir negatif özdeğere sahiptir. Açık bir yöntemin sınırlı kararlılık bölgesi, fiziksel olarak sönümlendikten çok sonra bile o modu çözmek için adım boyutunu zorlamakta, bu da hesaplamayı umutsuzca yavaşlatmaktadır. Geriye doğru Euler yöntemi, örtük Runge-Kutta şemaları ve geri fark formülleri gibi örtük yöntemler, sol yarı düzlemi (veya çoğunu) kapsayan kararlılık bölgelerine sahiptir; bu nedenle büyük adımlarda kararlı kalmakta ve adım boyutunun yalnızca doğruluk tarafından seçilmesine izin vermektedir. Her adım daha sonra (genellikle doğrusal olmayan) bir cebirsel sistemin çözümünü gerektirmekte olup, bu genellikle Jacobi matrisi kullanılarak bir Newton iterasyonu ile yapılmaktadır.
Klinik önem
Katılık, kimyasal reaksiyon ağları, yanma, elektrik devreleri, kontrol sistemleri ve parabolik kısmi diferansiyel denklemlerin çizgi yöntemi ayrıklaştırmalarında yaygındır. Katılığı tanımak ve uygun şekilde kararlı bir örtük çözücü seçmek, sonuçları makul bir sürede elde etmek için esastır ve çoğu üretim ODE yazılımı, otomatik katılık tespiti ve geçişi içermektedir.
Tarihçe
Katılık kavramı, 1952 yılında Curtiss ve Hirschfelder tarafından tanımlanmıştır ve destekleyici kararlılık teorisi — A-kararlılığı ve mertebe engelleri — Dahlquist tarafından geliştirilmiştir. Gear'ın geri fark formülü kodları ve daha sonraki yüksek mertebeli örtük Runge-Kutta yöntemleri, katı ve diferansiyel-cebirsel problemler için pratik araç setini oluşturmuştur.
Öne çıkan isimler
- Germund Dahlquist
- C. William Gear
- Ernst Hairer
- Gerhard Wanner
İlgili konular
Temel eserler
- hairer1996
- iserles2008
Sıkça sorulan sorular
- Bir adi diferansiyel denklemi tam olarak ne katı yapar?
- Katılık, sistemin, ilgilenilen çözümün evrimleşmesinden çok daha hızlı bozunan bileşenlere sahip olması durumunda ortaya çıkar. Tek bir kesin tanımı yoktur, ancak pratik göstergesi, doğruluk büyük adımlara izin verse bile açık yöntemlerin kararlılık için çok küçük adımlar kullanmak zorunda kalmasıdır.
- Katı problemler neden örtük yöntemler gerektirir?
- Örtük yöntemler, tüm sol yarı düzlemi (A-kararlılığı) kapsayan kararlılık bölgelerine sahip olabilir, böylece hızlı bozunan modlar için büyük adım boyutlarında kararlı kalırlar. Açık yöntemler, küçük adımları zorlayan ve onları katı problemler için pratik olmayan sınırlı kararlılık bölgelerine sahiptir.