Navier-Stokes denklemleri

Navier-Stokes denklemleri, ismini Claude-Louis Navier ve George Gabriel Stokes'tan almış olan, sıvılar ve gazlar gibi akışkanların hareketini tanımlamaya yarayan bir dizi denklemden oluşmaktadır.Bu denklemler; akışkan içerisindeki birim kütleye etki eden momentum (ivmelenme) değişimlerinin, basınç değişimleri ve sürtünme kayıplarına neden olan viskoz kuvvetlerin (sürtünmeye benzer) toplamına eşit olduğunun doğruluğunu ortaya koymaktadır. Bu viskoz kuvvetler moleküller arası etkileşimlerden meydana gelmekte ve akışkanın akmaya ne kadar dirençli (viskoz) olduğunu göstermektedir. Böylece, Navier-Stokes denklemlerinin, verilen akışkanın herhangi bir bölgesindeki kuvvetler dengesinin dinamik ifadesi olduğu söylenebilir.Bu denklemler en kullanışlı denklemlerin başında gelmektedirler. Çünkü, gerek akademik gerekse ekonomik birçok fenomenin fiziğini açıklamaktadır. Hava akımları ve okyanus akıntılarının, boru içindeki su akışının, galaksideki yıldız hareketlerinin, kanat etrafındaki hava akımlarının modellenmesinde ve hesaplarında sıkça kullanılırlar.Navier-Stokes denklemlerinin detayına girmeden önce, akışkanlar hakkında bazı kabuller yapılması gereklidir. Öncelikle akışkanın sürekli olduğu kabul edilir. Yani akışkanın tamamının aynı özellikte olduğu içinde farklı biçimler (formlar) bulunmadığı kabul edilir. Bir başka gerekli kabul de konu ile ilgili tüm alanların basınç, hız, yoğunluk, sıcaklık vs., diferansiyel olduğudur. (faz değişimleri olmadan)Bunun için, bazı hallerde kontrol hacmi adı verilen, rastgele seçilmiş sonlu bir hacim belirlemek gereklidir, bu hacim üzerinde bu prensipler kolayca uygulanabilir. Bu sonlu hacim  ile gösterilir ve yüzeyi sınırlandırılır . Kontrol hacmi, sabit kalabilir veya akışkan ile hareket edebilir. Temel kabuller bunlardır, bununla beraber, farklı uygulamalarda özel kabuller de yapılabilir.Hareket eden akışkanın özelliklerinin değişiminin ölçülebilmesi için iki yol vardır. Örneğin dünya atmosferindeki rüzgar hızının değişimleri ele alınacak olursa; bu değişiklikler bir meteoroloji istasyonu ölçüm cihazı (anemometre) veya bir hava balonu yolu ile ölçülebilir. Şüphesiz, ilk durumdaki anemometre boşlukta sabit bir nokta boyunca geçiş yapan tüm hareketli parçacıkların hızını ölçerken, ikinci durumda bahsedilen aygıt akışkan ile beraber hareket ederken hızdaki değişimi ölçer.Aynı durumda, yoğunluk, sıcaklık vb. değişimler de ölçümü etkileyecektir. Bu nedenle, bu iki hal için bir ayrım yapılmalıdır. Bir alanın boşluktaki sabit bir pozisyona göre türevi uzaysal (spatial) veya Euleryen türev (Eulerian derivative) olarak adlandırılır. Hareketli bir parçacığın izlenmesi türevi gerçek (substantive), Lagrangyan (Lagrangian) veya maddi (material) türev olarak adlandırılır.Burada  akışkanın hızıdır. Denklemin sağ tarafındaki ilk terim alışılmış Euleryen türevi (sabit bir referans üzerindeki türev) iken, ikinci terim akışkan hareketi ile oluşan değişiklikleri ifade eder. Bu etki adveksiyon olarak adlandırılır.En genel biçimde, bir korunum kanunu şunu ifade eder, bir kontrol hacmi üzerinde tanımlanmış hacim özelliği (bulk property) değişiminin oranı  hacim sınırları boyunca hareket eden akışkanın dışarı taşıdığı kayıp ve artı kontrol hacminin iç tarafındaki kazançlar ve kayıplara eşit kabul edilir. Bu, aşağıdaki integral denklemi ile ifade edilir.Eğer kontrol hacmi boşluk içinde sabitlenmiş ise bu integral denkleminden aşağıdaki şekilde bir ifade yazılabilir.Ayrıca, kontrol hacminin içinde, bu son denklemde elde edilmiş olan sağ taraftaki ilk terimin ifade edilmesi için diverjans teoremi kullanılmıştır. Böylece:Yukarıdaki ifade boşlukta sabit kalan bir kontrol hacminde  için geçerlidir. Çünkü  zaman içinde sabittir, değişmez. Bu sayede "" ve "" ifadeleri birbirinin yerine yazılabilir. Böylece ifade tüm alanlar için geçerli olur, ve integral çıkartılabilir.Momentumun korunumu, yoğunluk yerine momentumun vektör bileşenleri ve akışkan üzerine etkiyen kuvvetler ile, süreklilik denklemine benzer bir yaklaşım yapılarak ifade edilir. Süreklilik denkleminde  yerine belirli bir yönde birim hacim başına net momentum yazılır, , burada  hızın  bileşenidir. (hız x, y veya z yönleri boyunca olmak üzere), akışkan üzerine etkiyen kuvvetin  bileşenidir (her birim hacim başına gerçek kuvvet). Genel kuvvetler yerçekimi ve basınç gradyenlerini kapsar. Bu şu şekilde de ifade edilebilir:tensörü, akışkan parçacığı üzerine uygulanmış yüzey kuvvetleri olarak tanımlanır (gerilme tensörü). Akışkan girdap gibi bağımsız bir eğme bükme hareketi yapmadıkça,  simetrik bir tensördür. Genel olarak, biçim:Burada  normal gerilmeler,  teğetsel gerilmeler (kesme gerilmeleri) ve p gerilme tensörünün izotropik parçası ile birleştirilmiş statik basınçtır.matris izi (İng. trace) akışkanın dengede olup, olmadığı mutlaka tanımlanması (hacim vizkozitesi (bulk viscosity) olmadıkça) ile daima -3p'dir.Bu denklemler hala tamamlanmamıştır. Tamamlamak için, 'nin şekli üzerinde bir varsayım yapılmalıdır, şöyle ki, gerilme tensörü için aşağıda gösterildiği gibi bir süreklilik kanununa ihtiyaç vardır.Akış, sürekli ve diferansiyel kabul edilmiş ve korunum kanunları çerçevesinde kısmi diferansiyel denklemler ile ifade edilmiştir. Akışın sıkıştırılamaz (sabit yoğunluk) olduğu durumda, değişkenler, basınç ve hız bileşenleri için çözülmüştür. Bu değişkenler, Navier-Stokes denklemlerinin üç bileşeni, kütlenin korunumu (süreklilik denklemi) ilave edilerek, kapalı bir sistem için kısmi diferansiyel denklemler ile, sınır şartlarına uygun olarak çözülebilir. Sıkıştırılamaz akış durumunda, yoğunluk sistem için diğer bir bilinmeyen haline gelir, sistem için bir durum denklemi ilavesi ile saptanır. Durum denkleminde genelde akışkanın sıcaklığı işin içine girer, o yüzden denklem enerjinin korunumu için de mutlaka çözülmelidir. Bu denklemler non-lineer'dir (yani lineer değildir) ve kapalı formdaki analitik çözümleri sadece çok basit sınır şartları için bilinir.Denklemler, akım ve girdap fonksiyonu ikinci değişkenleri için Wilkinson denklemlerine dönüştürülebilirdir. Çözüm akışkan özelliklerine (viskozite, özgül ısı ve ısıl iletkenlik gibi) ve çalışma alanındaki sınır şartlarına bağlıdır.Denklem akışkanlarla ilgili problemlerin çözümü için, genel bazı durumlar için sadeleştirilip, genelleştirilerek kullanılabilir.Buradan denklemi türetebilmek için, öncelikle denge hali ifade edilir, pij=-pδij. Newtonyen bir akışkan için, bu denge değerinden gerilim tensörünün sapması, hızın gradyeni içinde lineerdir. Galile sabiti (Galilean covariance) nedeni ile açık şekilde hız üzerinde bağımlı değildir. Diğer bir ifade ile, pij+pδij,  de lineerdir. Akışkanların dönme sabiti belirlenir (sıvı kristal (liquid crystal) olmayanlar). pij+pδij izli ve izsiz simetrik tensörlerine ayrılır. Benzer olarak  izli, izsiz simetrik ve antisimetrik tensorlere ayrılır. Antisimetrik parça sıfıra gider, izli parça ve izsiz simetrik parçaya uygun iki katsayı vardır.  nin izsiz simetrik parçası,  dir, burada d uzaysal ölçü sayısıdır ve izli parça  dır. Bu nedenle, en genel lineer dönme sabiti şu şekilde verilir;μ ve μB bazı katsayılardır. μ kesme vizkozitesi (shear viscosity) ve μB hacim vizkozitesi (bulk viscosity) olarak adlandırılır. Bu ampirik (deneysel) bir incelemedir, hacim vizkozitesi çoğu akışkan için ihmal edilebilirdir, bu nedenle çoğu zaman ihmal edilir.Denklem içinde −2/3 ile çarpım görünmesi bununla açıklanır. Bu çarpım, 1 veya 2 uzaysal boyut içinde değiştirilebilir.yüksek süpersonik ve hipersonik uçuşlar gibi sıra dışı örnekler hariç, çoğunlukla ihmal edilebilirdir.Altı bilinmeyen (u, v, w, T, e and ) ve altı denklemden oluşan yukarıdaki gibi bir çözüm sistemi elde edilmiş olur.Bunlar, akış başlamadan önce bir miktar kesme dayanım kabiliyetleri olan akışkanlardır. Örnek olarak, diş macunu verilebilir.(Eğer  fakat akışkan sıkıştırılabilir ise sonuçta Euler denklemleri olarak bilinen denklemler elde edilir; burada, önemli olan sıkıştırılabilir akış ve akış içindeki şok dalgalarıdır.)Şunu ifade etmek gerekir ki, Navier-Stokes denklemleri akışkan akışını sadece yaklaşık olarak tanımlayabilir ve çok küçük ölçeklerde veya sıra dışı şartlarda, gerçek akışkanlar diğer maddeleri ve molekülleri içeren karışımlardır, Navier-Stokes denklemleri ile homojen ve sürekli akışlar modellenmiş ve bunun üzerinden sonuçlar elde edilmiştir. Bununla beraber Navier-Stokes denklemleri pratikteki problemlerin çözümü için, geniş bir aralıkta faydalı olur.