Üniform Derinlikteki Kanal Akımları ve En Uygun Kesit

6.6. Üniform Derinlikteki Kanal Akımları

Açık kanalların çoğunda, özellikle tarımsal sulama, drenaj hatları, taşkın kontrol sistemleri ve şehirsel su iletimi amacıyla inşa edilen altyapı yapılarında, sıvıların sabit ve üniform bir derinlikte akması genel bir tasarım prensibi olarak tercih edilir. Akış sırasında kanaldaki suyun hızı doğrudan serbest yüzeyde değil, rüzgar sürtünmesi ve yüzey gerilimi gibi etkilerden ötürü serbest yüzeyin hemen altında maksimum değere ulaşmaktadır. Buna karşılık, viskozite ve sınır tabakası etkileri nedeniyle kanal duvarlarında ve tabanında, yani suyun katı cidar ile temas ettiği her noktada akış hızı sıfır (kaymama koşulu) olmaktadır.

Açık kanal akımlarında, mühendislik yaklaşımlarını basitleştirmek adına hız profilini ve cidar (duvar) kayma gerilmesini sabit kabul ederiz. Bu kabuller ışığında, akışkanın debisi (verdi), kesit alanı ve taban eğimi arasındaki matematiksel, fiziksel ve hidrolik ilişkiler ayrıntılı olarak incelenebilir. Akımın üniform olduğu; diğer bir ifadeyle, zamana ve akış yönündeki (x ekseni) mesafeye göre derinliğin, kesit alanının ve ortalama akış hızının değişmediği varsayımı, kanal içerisindeki kuvvet dengelerini kurmamızı sağlar ve hidrodinamik denklemlerin çözümünü mümkün kılar.

6.6.1. Chezy ve Manning Eşitlikleri

Açık kanallarda akışkanın debisinin (verdisinin) doğru bir şekilde hesaplanması hidroloji ve hidrolik mühendisliğinin temelini oluşturur. Bu konu çok eski yıllara dayanmakta olup, tarihsel süreç içerisinde bu konuda pek çok ampirik (deneysel) ve yarı-ampirik formül geliştirilmiştir. Kontrol hacmi yaklaşımı kullanılarak, akış yönü olan x ekseni doğrultusunda Newton'un ikinci hareket kanunundan türetilen momentum eşitliği uygulanarak akış modellenir.

Düzenli (daimi) ve üniform bir akımda, ardışık (1) ve (2) numaralı kesitlerde momentum değişimi sıfırdır, zira akışkanın ivmelenmesi söz konusu değildir ve ortalama kesitsel hızlar birbirine mutlak surette eşittir ($V_1 = V_2$). Bu ideal durumda kontrol hacmine etki eden tüm dış kuvvetlerin toplamı sıfır olmalıdır ($\sum F_x = 0$). Momentum denkleminden yola çıkarak şu temel ifade yazılabilir:

$$ F_1 - F_2 - \tau_w \cdot P \cdot L + W \cdot \sin\theta = 0 $$

Bu formülasyonda yer alan parametrelerin fiziksel karşılıkları şu şekildedir: * $F_1, F_2$: Kontrol hacminin başlangıç ve bitiş uçlarındaki hidrostatik basınç kuvvetleridir. Akım üniform olduğu ve derinlik sabit kaldığı için ($y_1 = y_2$), karşılıklı kesitlerdeki basınç dağılımları aynıdır ve $F_1 = F_2$ olur, böylece birbirlerini yok ederler. * $\tau_w \cdot P \cdot L$: Akışkan tabakaları ile kanalın katı cidar yüzeyi arasındaki etkileşimden (sürtünmeden) doğan toplam direnç veya kayma kuvvetidir. Burada $P$ akışkanın temas ettiği ıslak çevreyi, $L$ kanalın söz konusu kesiminin uzunluğunu ve $\tau_w$ ortalama duvar kayma gerilmesini ifade eder. * $W \cdot \sin\theta$: Kontrol hacmindeki su kütlesinin ağırlığının, kanal eğimi yönündeki (hareket ettirici) bileşenidir. Suyun toplam ağırlığı $W = \rho \cdot g \cdot A \cdot L$ olarak hesaplanır; burada $\rho$ yoğunluğu, $g$ yerçekimi ivmesini ve $A$ kesit alanını temsil eder.

Kuvvet dengesinde hidrostatik basınç kuvvetleri elendiğinde, sistemi hareket ettiren ağırlık bileşeninin tümüyle cidar sürtünmesine harcandığı görülür ve denklem şu şekilde yeniden düzenlenir:

$$ \tau_w = \frac{W \cdot \sin\theta}{P \cdot L} $$

Pratik mühendislik uygulamalarında kanalın taban eğimi genelde çok küçük (örneğin %1'in altında) olduğundan, trigonometrik bir yaklaşım ile $\sin\theta \approx \tan\theta = S_0$ (taban eğimi) alınabilir. İlaveten, hidrolik yarıçap $R = A/P$ (kesit alanının ıslak çevreye oranı) ve akışkanın özgül ağırlığı $\gamma = \rho \cdot g$ tanımları denkleme yerleştirildiğinde, kayma gerilmesi ile kanal geometrisi ve eğimi arasında temel bir ilişki kurulur:

$$ \tau_w = \gamma \cdot R \cdot S_0 $$

Açık kanallardaki akışlar doğası gereği genellikle türbülanslıdır. Yüksek Reynolds ($Re$) sayılarında duvar kayma gerilmesinin, ortalama akış hızının karesini içeren kinetik enerji terimiyle orantılı olduğu deneysel olarak kanıtlanmıştır. Yani, $\tau_w = k \cdot \rho \cdot V^2 / 2$ ilişkisi mevcuttur ($k$ burada sürtünme faktörüdür). Bu iki eşitlik birbirine eşitlendiğinde Chezy denkleminin çekirdeği elde edilir:

$$ \frac{k \cdot \rho \cdot V^2}{2} = \rho \cdot g \cdot R \cdot S_0 $$ $$ V = \left(\frac{2g}{k}\right)^{1/2} \cdot \sqrt{R \cdot S_0} $$

Bu genel hız denkleminde yer alan kareköklü sabit ifade, yani $(2g/k)^{1/2}$ terimi, ünlü Fransız mühendis Antoine de Chézy'ye ithafen Chezy katsayısı ($C$) olarak tanımlanmıştır. Böylece hidrolik biliminin en temel eşitliklerinden biri olan Chezy Eşitliği ortaya çıkar:

$$ V = C \cdot \sqrt{R \cdot S_0} $$

Chezy formülü, yapısının basitliği sayesinde oldukça pratik ve kullanışlı bir araçtır. Ne var ki, en büyük handikapı (sakıncası) $C$ katsayısının sabit bir sayı değil, boyutlu (metrik sistemde $m^{1/2}/s$) ve kanal pürüzlülüğü ile hidrolik yarıçapa bağımlı değişken bir parametre olmasıdır. Tarihsel gelişim içinde $C$ katsayısını daha güvenilir bir şekilde belirlemek maksadıyla Bazin, Kutter ve Ganguillet gibi pek çok bilim insanı uzun yıllar süren denemeler yaparak çeşitli katsayı tabloları, Bazin katsayısı ($Bz$) ve Kutter katsayısı ($m$) gibi yeni değişkenler geliştirmişlerdir.

Ancak hidrolik biliminde açık kanallar için bugün tüm dünyada standart haline gelmiş ve en yaygın olarak kullanılan formül, İrlandalı mühendis Robert Manning (1889) tarafından önerilen Manning Eşitliği'dir. Manning, deneysel veri setlerini incelediğinde, $C$ Chezy katsayısının, yüzey pürüzlülük katsayısının ($n$) tersi ile hidrolik yarıçapın $1/6$ kuvvetinin çarpımına eşit bir trend izlediğini ampirik olarak kanıtlamıştır:

$$ C = \frac{1}{n} R^{1/6} $$

Bu keşfi Chezy denkleminde yerine koyarsak, açık kanallarda en çok başvurulan ve endüstri standardı kabul edilen ortalama hız denklemi olan Manning Hız Formülü'nü elde ederiz:

$$ V = \frac{1}{n} R^{2/3} S_0^{1/2} $$

Manning hız ifadesinden yola çıkarak hacimsel akış debisi (verdi) formülüne ($Q = V \cdot A$) geçildiğinde ise akım şiddetini veren nihai denklem ortaya çıkar:

$$ Q = \frac{1}{n} A \cdot R^{2/3} \cdot S_0^{1/2} $$

Bu formüllerdeki birim ve anlamlar şu şekildedir: * $V$: Akışkanın ortalama kesit hızı (m/s) * $n$: Manning pürüzlülük katsayısı (s/m$^{1/3}$). Pürüzlülük derecesini yansıtır. * $R$: Hidrolik yarıçap (m) * $S_0$: Kanal taban eğimi (boyutsuz veya m/m) * $Q$: Debisi ya da verdi (m³/s) * $A$: Kanal kesit alanı (m²)

6.6.2. En Uygun Kesitin Saptanması (Hidrolik Açıdan En Ekonomik Kesit)

Chezy ve özellikle Manning eşitlikleri detaylı olarak incelendiğinde; eğer inşa edilecek bir kanalın taban eğimi ($S_0$) ile yüzey malzemesinin pürüzlülük katsayısı ($n$) sabit tutulursa, hızın ve dolayısıyla kanaldan iletilecek debinin ($Q$) doğrudan hidrolik yarıçapa ($R$) bağlı olduğu açıkça görülür. Belirli, sabit bir kazı alanı ($A$) kısıtlaması altında debinin maksimum yapılabilmesi için, ıslak çevrenin ($P$) mutlak surette minimum değerine indirilmesi gerekir. Zira çevre azaldıkça sürtünme yüzeyi azalır ve akışkan daha az dirençle karşılaşır.

Sıvı iletim kapasitesini maksimuma çıkartan (aynı pürüzlülük ve eğimdeki alternatiflerine göre en az enerji kaybı yaratan) kesite hidrolik literatüründe "en uygun kesit" veya "en ekonomik kesit" adı verilmektedir. Minimum ıslak çevre, kazı miktarını ve beton (veya diğer malzeme) astarlama maliyetini minimize edeceği için, hidrolik optimizasyon aynı zamanda devasa bir ekonomik optimizasyon demektir.

Matematiksel ve geometrik olarak bilindiği üzere, düzlemde belirli bir alan için etrafını saran çevresi en küçük olan şekil kusursuz bir dairedir. Bu nedenle açık kanallar açısından da teorik olarak en yüksek debiyi iletebilecek en ideal yapı, su ile tam sınır hizasında olan yarım daire şeklindeki kesitlerdir. Fakat gerek hafriyat makinelerinin (ekskavatörlerin) geometrisi, gerekse kalıp döküm işçiliğindeki zorluklar nedeniyle uygulamada çoğunlukla inşaası ve bakımı kolay olan dikdörtgen, yamuk (trapez) ve dairesel boru kesitleri tercih edilmektedir.

Dikdörtgen Kanallarda En Uygun Kesit

Kanal taban genişliği $b$, akışkan (su) derinliği $y$ olan standart bir dikdörtgen kesitli kanalda hidrolik parametreler oldukça basittir. Kesit alanı $A = b \cdot y$ ve akışkanın temas ettiği ıslak çevre $P = b + 2y$ olarak ifade edilir. Hidrolik yarıçap $R$ ise alanın çevreye bölünmesiyle $R = \frac{b \cdot y}{b + 2y}$ olur.

Amacımız alanı ($A$) sabit kabul ederek çevreyi ($P$) minimize etmektir. Bunun için alan denkleminden $b$ çekilerek çevre denklemine yerleştirilir: $b = A/y \implies P = A/y + 2y$. Çevre fonksiyonunun su derinliğine ($y$) göre türevi alınıp minimum noktasını bulmak için sıfıra eşitlendiğinde:

$$ \frac{dP}{dy} = -\frac{A}{y^2} + 2 = 0 \implies A = 2y^2 $$

Bulunan optimizasyon şartı $A = b \cdot y$ temel bağıntısında yerine konulursa; $b \cdot y = 2y^2 \implies b = 2y$ sonucu şüpheye yer bırakmaksızın elde edilir.

Sonuç olarak; dikdörtgen kanallarda en uygun kesit formu, kanal taban genişliğinin, içinden akan su derinliğinin tam iki katı ($b = 2y$) olarak tasarlandığı durumdur. Bu kritik koşulda kanalın hidrolik yarıçapı $R = y/2$ olmaktadır. Alan sabitken maksimum debi bu geometrik kural gözetilerek elde edilir.

Yamuk (Trapez) Kanallarda En Uygun Kesit

Mühendislikte toprak veya stabilize kanallarda, yan duvarların dik yapılması cidarın (şev) zamanla çökmesine neden olur. Göçmeleri engellemek için eğimli kenarlar uygulanır ve ortaya yamuk (trapez) kesit çıkar. Yanal kenarların eğimi (şev açısı), yatayla yapılan açının tanjantı bağlamında $1/m = \tan\theta$ olarak gösterilir. Simetrik bir yamuk kesitte temel hidrolik formüller şunlardır: * Kesit alanı: $A = y(b + m \cdot y)$ * Islak çevre: $P = b + 2y\sqrt{1+m^2}$

Yine belirli bir $A$ alanı için cidar sürtünmesini karakterize eden ıslak çevreyi ($P$) minimum yapan koşulu bulmak için ileri türev işlemleri uygulandığında, hedeflenen sıvı derinliği ($y$) şu şekilde elde edilir:

$$ y = \left( \frac{A}{2\sqrt{1+m^2} - m} \right)^{1/2} $$

Tasarım mühendisi hem su derinliğini hem de şev eğimini (açıyı) seçme özgürlüğüne sahipse, hidrolik olarak kusursuz olan mükemmel yamuk kesitin bulunması için $m$ bağımsız değişkenine göre ikinci bir optimizasyon yapılır. Bu matematiksel analiz sonucunda tabiatta bal peteklerinde de görülen altıgenin bir yarısı olan şekle, yani en iyi trapez açısının $\theta = 60^\circ$ olduğu değere ulaşılır. Bu durumda $m$ değeri $m = \cot(60^\circ) = 1/\sqrt{3}$ olmaktadır.

Şev açısının 60 derece alındığı ideal ve en uygun yamuk kesitin nihai optimizasyon özellikleri şu bağıntılara dönüşür: $$ A = y^2 \sqrt{3} $$ $$ b = \frac{2y}{\sqrt{3}} $$ $$ P = 2y\sqrt{3} $$ ve her ideal açık kanal geometrisinde olduğu gibi hidrolik yarıçap yine $R = y/2$ bulunur.

Dairesel Kesitli Kanallarda Optimum Akış Koşulları

Şehir altyapılarında yaygın olarak kullanılan atıksu (kanalizasyon) kanalları, tam kapalı (dairesel) boru profilleridir. Ancak bu borular basınçlı değil, cazibe (yerçekimi) etkili serbest yüzeyli akım prensibiyle çalışırlar. Yaygın bir yanılgı, bu kanallarda debi ve hızın maksimum olduğu derinlik seviyesinin, borunun tamamen dolu olduğu (%100 doluluk) duruma denk geldiğinin düşünülmesidir.

Gerçekte ise boru tepe noktasına doğru doldukça, kesit alanı artmasına rağmen ıslak çevre (sürtünme alanı) oran olarak çok daha şiddetli, keskin bir şekilde büyür. Tavan sürtünmesi su akışını ciddi şekilde baskılar. Çözümlemeler ve grafiksel analizler dairesel bir kanalda, tamamen dolu kesitten kaynaklanan aşırı sürtünmenin olmadığı ve maksimum debinin elde edildiği eşik noktasının, boru çapının (%100'ü) değil, su seviyesinin tam olarak $y = 0.938D$ (çapın %93.8'i) veya merkeze göre $303^\circ$ açı taradığı durumda gerçekleştiğini matematiksel olarak kanıtlamıştır.

6.6.3. Chezy ve Manning Formüllerinin Mühendislik Uygulamalarındaki Pratik Önemi

Mühendislik projelerinde, sulama şebekelerinden, devasa hidroelektrik santrallerinin deşarj tünellerine kadar çok geniş bir spektrumda bu eşitlikler aktif olarak kullanılmaktadır. Hidrolik tasarım safhasında doğru yüzey pürüzlülük katsayısı ($n$) seçimi hayati önem arz eder. Yanlış, iyimser bir $n$ katsayısı seçimi (örneğin pürüzlü ve otlu bir yüzeyin düzgün perdahlı beton gibi kabul edilmesi), sahadaki uygulamanın kağıt üzerindeki verileri tutmamasına, suyun çok daha düşük bir hızda akmasına, kanal kapasitesinin yetersiz kalmasına ve nihayetinde felaket niteliğinde çevresel taşkınlara sebep olabilir.

Manning pürüzlülük katsayısı $n$; çok temiz ve pürüzsüz cam/plastik yüzeylerde $0.010$ mertebelerindeyken; pürüzlü moloz taş veya kaba beton duvarlarda $0.035 - 0.040$, çok otlanmış bakımsız doğal dere yataklarında ve sel yayılım düzlemlerinde ise bitki örtüsünün direncinden ötürü $0.150$ gibi devasa değerlere ulaşabilir. Bu, su akış hızının sadece pürüzlülük sebebiyle teorik maksimuma kıyasla 10 kattan fazla azaldığı anlamına gelmektedir.

Tasarımlarda sadece pürüzlülük değil kanal optimizasyonu da kritiktir. Su doğada kendi yolunu açarken hidro-morfolojik olarak yatağını yarım daire veya hidrolik açıdan optimum yamuk profiline ulaştırma eğilimindedir. İnşaat mühendislerinin ve tasarımcıların yapay betonarme kanallarda optimum dikdörtgen profili ($b=2y$) veya optimum yamuk profili (şev açısı $60^\circ$) şiddetle tercih etmeleri yalnızca kağıt üzerindeki bir matematiksel formül başarısı değil, aynı zamanda devasa altyapı yatırımlarında hafriyat, beton kalıp işçiliği ve kaplama malzeme sarfiyatlarını milyonlarca lira ölçeğinde minimize eden muazzam bir mühendislik ekonomisi harikasıdır.