Laplace Dönüşümleri ve Diferansiyel Denklemlerin Sayısal Çözümleri ile ilgili Genel Bilgiler

Özet

Diferansiyel denklemler, mühendislik ve temel bilimlerde fiziksel olayların matematiksel olarak modellenmesinde temel bir araçtır. Ancak bu denklemlerin doğrudan çözümü, özellikle yüksek mertebeden, başlangıç değerli veya parçalı tanımlanmış fonksiyonlar içeren durumlarda, cebirsel denklemlere kıyasla oldukça zordur. Bu nedenle, diferansiyel denklemleri daha kolay çözülebilir bir biçime dönüştürmenin yollarının aranması doğaldır. Bu amaçla geliştirilen yöntemlerden biri, diferansiyel denklemi uygun bir integral dönüşümü yardımıyla cebirsel bir denkleme dönüştürmektir. Aslında bir diferansiyel denklemi cebirsel bir denkleme dönüştürmenin birden fazla yolu vardır. Bu tür yöntemlerin ortak özelliği, diferansiyel denklemdeki her terimi belirli bir çekirdek (kernel) fonksiyon ile çarparak, denklemin tanım bölgesi boyunca integre etmeleridir. Bu işlem sonucunda, bilinmeyen fonksiyonun türevlerini içermeyen, yalnızca dönüşüm uzayında tanımlı cebirsel bir denklem elde edilir. Bu yaklaşımın en önemli avantajı şudur: Diferansiyel denklemin çözümü, dönüşüm uzayında cebirsel olarak bulunur ve daha sonra ters dönüşüm uygulanarak zaman uzayındaki gerçek çözüm elde edilir. Bu nedenle bu yöntemler genel olarak integral dönüşümleri olarak adlandırılır. Her integral dönüşümü, kendine özgü integral sınırlarına ve çekirdek fonksiyonlara sahiptir ve belirli tip problemlerde etkilidir. Bu makalede ele alınan Laplace dönüşümü, integral dönüşümleri arasında en yaygın ve en güçlü olanlardan biridir. Laplace dönüşümü, özellikle zaman değişkeni için tanımlanmış başlangıç değer problemlerinde ve yarı sonsuz geometrilerde (bağımsız değişkenin sıfırdan sonsuza uzandığı problemler) son derece etkilidir. Bu özellik, Laplace dönüşümünü mühendislik uygulamaları için ideal hale getirir. Laplace dönüşümü, başta sabit katsayılı lineer diferansiyel denklemler ve bu tür denklemlerden oluşan denklem sistemlerinin çözümünde kullanılacaktır. Özellikle, diferansiyel denklemin sağ tarafında homojen olmayan terimler, ani değişimler, sıçrama süreksizlikleri veya parçalı tanımlı fonksiyonlar bulunduğu durumlarda, Laplace dönüşümü klasik çözüm yöntemlerine göre çok daha büyük bir basitlik sağlar. Diferansiyel denklemler, doğa olaylarını, mühendislik sistemlerini ve fiziksel süreçleri matematiksel olarak ifade etmenin en güçlü araçlarından biridir. Bu denklemler; hareket, ısı transferi, akışkanlar mekaniği, elektrik devreleri ve kontrol sistemleri gibi birçok alanda karşımıza çıkar. Şimdiye kadar diferansiyel denklemler için ele alınan yöntemler genellikle analitik çözüm yöntemleri olmuştur ve bu yöntemlerle elde edilen çözümler analitik çözüm olarak adlandırılmıştır. Analitik çözümler, çözüm fonksiyonunun bağımsız değişkene açık veya kapalı biçimde bağlı olduğu durumlardır. Örneğin 𝑦 = 𝑥2 gibi açık biçimde ifade edilebilen çözümler, bağımsız değişkenin herhangi bir değerinde doğrudan hesap yapılmasına olanak tanır. Bu tür çözümler, hem yoruma açık olmamaları hem de kesin sonuçlar vermeleri nedeniyle matematiksel açıdan son derece değerlidir. Bununla birlikte bazı analitik çözümler kapalı biçimde ifade edilir ve bu durumda çözüm, açık bir fonksiyon şeklinde yazılamaz. Böyle çözümlerde, fonksiyonun davranışını incelemek için ek işlemler yapılması gerekir. Ancak analitik çözüm yöntemleri, pratikte karşılaşılan diferansiyel denklemlerin büyük bir kısmı için yetersiz kalmaktadır. Özellikle lineer olmayan, değişken katsayılı ya da karmaşık sınır ve başlangıç koşullarına sahip denklemler çoğunlukla analitik olarak çözülemez. Uygulamada karşılaşılan gerçek sistemlerin büyük bir bölümü bu tür denklemlerle ifade edildiğinden, analitik çözümlerin sadece istisnai durumlarda mümkün olduğu görülmektedir. Bu noktada, diferansiyel denklemlerin çözümü için yaklaşık çözüm yöntemleri devreye girer. Yaklaşık çözümler, çözümün tam ifadesini vermek yerine, belirli noktalar için çözüm değerlerinin hesaplanmasını amaçlar. Bu yöntemlerin temelinde, diferansiyel denklemin sürekli yapısının küçük adımlara bölünmesi ve çözümün adım adım ilerletilmesi fikri yer alır. Sayısal çözüm yöntemleri sayesinde, analitik olarak çözülemeyen birçok problem yeterli doğrulukta çözülebilir hale gelir. Sayısal çözümde elde edilen sonuçlar genellikle sayısal tablolar veya grafikler şeklinde sunulur. Bu yaklaşım, özellikle mühendislik uygulamalarında son derece kullanışlıdır. Çünkü çoğu zaman bir sistemin belirli zamanlardaki veya belirli konumlardaki davranışı, çözümün genel ifadesinden daha büyük önem taşır. Sayısal yöntemler bu ihtiyacı doğrudan karşılar. Bu bölümde ele alınan sayısal çözüm yöntemlerinin diferansiyel denklemlerle yakından ilişkili olması nedeniyle, konuya sayısal integral alma işlemi ile başlanmıştır. Bir diferansiyel denklemin çözümü, özünde bu denklemin integrasyonuna eşdeğer olduğundan, integralin sayısal olarak nasıl hesaplandığının anlaşılması büyük önem taşır. Bu nedenle dikdörtgen şerit yöntemi, trapez kuralı ve Simpson kuralı gibi temel sayısal integrasyon yöntemleri tanıtılarak, diferansiyel denklemlerin sayısal çözümüne bir zemin hazırlanmıştır. Sayısal çözüm yöntemleri arasında en basit olanı Euler yöntemidir. Bu yöntem, her ne kadar yüksek doğruluk sağlamasa da, sayısal çözüm mantığını kavramak açısından temel bir öneme sahiptir. Euler yöntemini takiben, daha yüksek doğruluk sağlayan düzeltilmiş Euler yöntemi ele alınmış ve ardından mühendislik uygulamalarında yaygın olarak kullanılan Runge–Kutta yöntemleri tanıtılmıştır. Özellikle dördüncü mertebeden Runge–Kutta yöntemi, doğruluk ve hesaplama maliyeti arasındaki dengesi nedeniyle öne çıkmaktadır. [1-53]