Stokes sayısı

Stokes sayısının farklı etkilerini gösteren bir illüstrasyon. Turuncu ve yeşil yörüngeler, sırasıyla küçük ve büyük Stokes sayıları için verilmiştir. Turuncu eğri, birden küçük Stokes sayısına sahip ve akış çizgilerini (mavi) takip eden bir parçacığın yörüngesini, yeşil eğri ise birden büyük Stokes sayısına sahip ve akış çizgilerini takip etmeyen bir parçacığın yörüngesini göstermektedir. Bu parçacık, sarı renkle gösterilen noktada engellerden birine (kahverengi daireler) çarpmaktadır.

Stokes sayısı (Stk), George Gabriel Stokes'un adını taşıyan ve parçacıkların bir akışkan akışı içerisinde süspansiyonda gösterdiği davranışı karakterize eden bir boyutsuz sayıdır. Stokes sayısı, bir parçacığın (veya damlanın) karakteristik zamanı ile akışın veya bir engelin karakteristik zamanı arasındaki oran olarak şu şekilde tanımlanır:

$\mathrm{S}\mathrm{t}\mathrm{k}=\frac{t_0{u}_0}{l_0}$

Burada, $t_0$ parçacığın relaksasyon süresi (İng. relaxation time; sürüklenme nedeniyle parçacık hızının üstel olarak azaldığı zaman sabiti), $u_0$ engelden uzakta akışın akışkan sürati ve $l_0$ engelin karakteristik boyutu (genellikle çapı) veya akıştaki karakteristik bir uzunluk ölçeğidir (örneğin sınır tabakası kalınlığı). Düşük Stokes sayısına sahip bir parçacık akış çizgilerini takip eder (mükemmel adveksiyon), oysa büyük Stokes sayısına sahip bir parçacık ataletinin etkisi altındadır ve başlangıç yörüngesini izlemeye devam eder.

Stokes akışı durumunda, yani parçacık (veya damlacık) Reynolds sayısı yaklaşık olarak birden küçük olduğunda, parçacığın sürükleme katsayısı Reynolds sayısının kendisiyle ters orantılıdır. Bu durumda, parçacığın karakteristik zamanı şu şekilde ifade edilebilir:

$t_0=\frac{{\rho }_p{d}_p^2}{18{\mu }_g}$

Burada, ${\rho }_p$ parçacığın yoğunluğu, $d_p$ parçacık çapı ve ${\mu }_g$ akışkanın dinamik viskozitesidir.

Deneysel akışkan dinamiğinde, Stokes sayısı, parçacık görüntü velocimetrisi (PIV) deneylerinde, çok küçük parçacıkların türbülanslı akışlara katılması ve akışkan hareketinin hızını ve yönünü (aynı zamanda akışkanın hız alanı olarak da bilinir) belirlemek amacıyla optik olarak gözlemlenmesi durumunda bir akış izleyici doğruluğu ölçütüdür. Kabul edilebilir izleme doğruluğu için, parçacık yanıt süresi, akışın en küçük zaman ölçeğinden daha hızlı olmalıdır. Küçük Stokes sayıları daha iyi izleme doğruluğunu gösterir; $\mathrm{S}\mathrm{t}\mathrm{k}\gg 1$ durumunda, parçacıklar özellikle akışın ani bir şekilde yavaşladığı yerlerde akıştan ayrılacaktır. $\mathrm{S}\mathrm{t}\mathrm{k}\ll 1$ durumunda, parçacıklar akış çizgilerini yakından takip eder. Eğer $\mathrm{S}\mathrm{t}\mathrm{k}<0.1$ ise, izleme doğruluğu hataları yüzde 1'in altında kalır.

Parçacık görüntü velocimetrisinde (PIV) relaksasyon süresi ve izleme hatası

PIV için izleme doğruluğu açısından iki farklı parçacık boyutunun karşılaştırılması. Propilen glikolün simüle edilmiş parçacıkları (mavi noktalar) durgunluk noktası akış alanında (gri akış çizgileri) taşınmaktadır. 1 mm parçacıkların durgunluk plakasına çarptığına, 0.1 mm parçacıkların ise akış çizgilerini takip ettiğine dikkat edin.

Stokes sayısı, daha önce tartışıldığı gibi, PIV veri setlerinin kalitesini değerlendirmek için bir ölçüt sunar. Ancak, her uygulamada karakteristik bir hız veya uzunluk ölçeğinin tanımı açık olmayabilir. Bu yüzden, izleme gecikmesinin nasıl ortaya çıktığını daha iyi anlamak için, Stokes rejimindeki bir parçacığın diferansiyel denklemlerinin tanımlanması gereklidir. Belirli bir hızla $v_p(t)$ hareket eden bir parçacık, akışkanın değişken hız alanıyla karşılaşacaktır. Parçacığın Lagrangian referans çerçevesinde akışkanın hızının $v_f(t)$ olduğunu varsayalım. Bu hızlar arasındaki fark, parçacığın yolunu düzeltmek için gerekli sürükleme kuvvetini oluşturur:

$\mathrm{\Delta }v(t)=v_f(t)-v_p(t)$

Stokes sürükleme kuvveti şu şekilde tanımlanır:

$F_D=3\pi \mu d_p\mathrm{\Delta }v$

Parçacığın kütlesi ise şu şekilde hesaplanır:

$m_p={\rho }_p\frac{4}{3}\pi (\frac{d_p}{2}{)}^3={\rho }_p\frac{\pi d_p^3}{6}$

Bu bağlamda, Newton'un ikinci yasasına göre parçacığın ivmesi şu şekilde bulunabilir:

$\frac{d{v}_p(t)}{dt}=\frac{F_D}{m_p}=\frac{18\mu }{{d_p}^2{\rho }_p}\mathrm{\Delta }v(t)$

Relaksasyon süresi olan $t_0=\frac{{\rho }_p{d}_p^2}{18{\mu }_g}$ yerine konularak, ifade şu hale getirilir:

$\frac{d{v}_p(t)}{dt}=\frac{1}{t_0}\mathrm{\Delta }v(t)$

Yukarıdaki birinci derece diferansiyel denklem, Laplace dönüşümü yöntemi ile şu şekilde çözülebilir:

$t_{0}s{v}_p(s)=v_f-v_p(s)$

$\frac{v_p(s)}{v_f(s)}=\frac{1}{t_{0}s+1}$

Bu çözüm, frekans alanında, karakteristik zamanı $t_0$ olan birinci dereceden bir sistemi tanımlar. Bu nedenle, −3 dB kazanç (kesim) frekansı şu şekilde olacaktır:

$f_{-3\text{ dB}}=\frac{1}{2\pi t_0}$

Kesim frekansı ve parçacık transfer fonksiyonu, yan panelde çizilen grafikte, dengesiz akış uygulamalarında PIV hatasının ve bu hatanın türbülans spektral nicelikleri ve kinetik enerji üzerindeki etkisinin değerlendirilmesini sağlar.

Havada farklı parçacık çapları için propilen glikol parçacığının Bode grafiği.

Şok dalgası üzerinden geçen parçacıklar

Önceki bölümde tartışılan parçacık izleme yanılma hatası, frekans alanında belirgindir, ancak parçacık hareketinin akış alanı ölçümleri yapmak amacıyla izlendiği durumlarda (örneğin, parçacık görüntü velocimetrisinde) bunu değerlendirmek zor olabilir. Yukarıda belirtilen diferansiyel denklemin basit fakat anlamlı bir çözümü, zorlama fonksiyonu $v_f(t)=V_u-\mathrm{\Delta }VH(t)$ bir Heaviside basamak fonksiyonu olduğunda mümkündür; bu, bir şok dalgası üzerinden geçen parçacıkları temsil eder. Bu durumda, $V_u$ şokun yukarısındaki akış hızıdır; $\mathrm{\Delta }V$ ise şok boyunca meydana gelen hız düşüşüdür.

Bir parçacığın basamak yanıtı basit bir üstel fonksiyondur:

$v_p(t)=(V_u-\mathrm{\Delta }V)+\mathrm{\Delta }V{e}^{-t/t_0}$

Hızı zamanın bir fonksiyonu olarak, parçacık hız dağılımını mesafenin bir fonksiyonu olarak dönüştürmek amacıyla, $x$ yönünde 1 boyutlu bir hız sıçraması olduğunu varsayalım. Şok dalgasının $x=0$ konumunda olduğunu varsayalım ve ardından önceki denklemi entegre edelim:

$x_{\text{particle}}=\underset{0}{\overset{\mathrm{\Delta }t}{\int }}{v}_p(t)dt=\underset{0}{\overset{\mathrm{\Delta }t}{\int }}(V_u-\mathrm{\Delta }V)dt+\underset{0}{\overset{\mathrm{\Delta }t}{\int }}\mathrm{\Delta }V{e}^{-t/t_0}dt$

$x_{\text{particle}}=\mathrm{\Delta }t(V_u-\mathrm{\Delta }V)+\mathrm{\Delta }t\mathrm{\Delta }V(1-e^{-\mathrm{\Delta }t/t_0})$

$\mathrm{\Delta }t=3t_0$ (hızın \% 95 değişim süresi) relaksasyon süresini dikkate alarak, şu sonucu elde ederiz:

$x_{\text{particle},95\mathrm{\%}}=3t_0(V_u-\mathrm{\Delta }V)+3t_0\mathrm{\Delta }V(1-e^{-3})$

$x_{\text{particle},95\mathrm{\%}}=3t_0(V_u-0.05\mathrm{\Delta }V)$

Bu durum, parçacık hızının, şok dalgasından itibaren $x_{\text{particle},95\mathrm{\%}}$ mesafede aşağı akış hızının $\mathrm{\%}5$'ine yerleşeceği anlamına gelir. Pratikte bu, bir PIV sistemi için şok dalgasının yaklaşık olarak bu $x_{\text{particle},95\mathrm{\%}}$ mesafesi kadar bulanık görüneceği anlamına gelir.

Örneğin, durgunluk sıcaklığı 298 K ve Mach sayısı $M=2$ olan normal bir şok dalgasını düşünelim. $d_p=1\mu \text{m}$ çapındaki bir propilen glikol parçacığı, akışı $x_{\text{particle},95\mathrm{\%}}=5\text{ mm}$ kadar bulanıklaştıracaktır; oysa $d_p=10\mu \text{m}$ çapındaki bir parçacık, akışı $x_{\text{particle},95\mathrm{\%}}=500\text{ mm}$ kadar bulanıklaştıracaktır (bu, çoğu durumda kabul edilemez PIV sonuçları verecektir).

Bir şok dalgası, akışın ani yavaşlamasının en kötü senaryosu olmasına rağmen, PIV'deki parçacık izleme hatasının etkisini ve bu hatanın hız alanlarının $x_{\text{particle},95\mathrm{\%}}$ mertebesindeki uzunluk ölçeklerinde bulanıklaşmasına neden olduğunu gösterir.

Non-Stokesian sürükleme rejimi

Önceki analiz, ultra-Stokesian rejiminde doğru olmayacaktır, yani parçacık Reynolds sayısı birden çok daha büyük olduğunda geçerli değildir. Mach sayısının birlikten çok daha küçük olduğunu varsayarak, Israel ve Rosner tarafından genelleştirilmiş bir Stokes sayısı formu gösterilmiştir.

$$\text{Stk}\text{e}=\text{Stk}\frac{24}{\text{Re}o}\underset{0}{\overset{{\text{Re}}_o}{\int }}\frac{d{\text{Re}}^{\prime }}{C_D({\text{Re}}^{\prime }){\text{Re}}^{\prime }}$$

Burada ${\text{Re}}_o$, "parçacık serbest-akış Reynolds sayısı" olarak tanımlanır,

$${\text{Re}}_o=\frac{{\rho }_g|𝐮|d_p}{{\mu }_g}$$

Ek olarak tanımlanan fonksiyon $\psi ({\text{Re}}_o)$ ise; non-Stokesian sürükleme düzeltme faktörünü ifade eder,

$$\text{Stk}e=\text{Stk}\cdot \psi (\text{Re}o)$$

Bu fonksiyon şu şekilde tanımlanır,

$\psi$, küresel bir parçacık için non-Stokesian sürükleme düzeltme faktörünü tanımlar

$$\psi (\text{Re}o)=\frac{24}{\text{Re}o}\underset{0}{\overset{{\text{Re}}_o}{\int }}\frac{d{\text{Re}}^{\prime }}{C_D({\text{Re}}^{\prime }){\text{Re}}^{\prime }}$$

Sınırlandırıcı parçacık serbest-akış Reynolds sayıları dikkate alındığında, $\text{Re}o\to 0$ olduğunda $C_D(\text{Re}o)\to 24/{\text{Re}}_o$ ve dolayısıyla $\psi \to 1$ olur. Bu nedenle, beklenildiği gibi Stokesian sürükleme rejiminde düzeltme faktörü birdir. Wessel & Righi $\psi$'yi Schiller & Naumann'ın bir küre üzerindeki sürükleme için ampirik korelasyonundan $C_D(\text{Re})$ kullanarak değerlendirmiştir.

$$\psi (\text{Re}o)=\frac{3(\sqrt{c}\text{Re}{o}^{1/3}-\arctan (\sqrt{c}\text{Re}{o}^{1/3}))}{c^{3/2}\text{Re}o}$$

Burada sabit $c=0.158$ olarak alınmıştır. Geleneksel Stokes sayısı, büyük parçacık serbest-akış Reynolds sayıları için sürükleme kuvvetini önemli ölçüde küçümseyecektir. Bu da parçacıkların akış yönünden sapma eğilimini fazla tahmin edecektir. Bu durum, sonraki hesaplamalarda veya deneysel karşılaştırmalarda hatalara yol açacaktır.

Anizokinetik parçacık örneklemesinin uygulanması

Örneğin, Belyaev ve Levin tarafından tanımlanan hizalanmış, ince duvarlı dairesel bir nozul aracılığıyla parçacıkların seçici olarak yakalanması şu şekilde ifade edilmektedir:

$c/c_0=1+(u_0/u-1)(1-\frac{1}{1+\mathrm{S}\mathrm{t}\mathrm{k}(2+0.617u/u_0)})$

burada $c$ parçacık konsantrasyonu, $u$ hız ve alt simge 0, nozulun oldukça yukarısındaki koşulları belirtir. Karakteristik mesafe nozulun çapıdır. Bu durumda Stokes sayısı şu şekilde hesaplanır,

$\mathrm{S}\mathrm{t}\mathrm{k}=\frac{u_0{V}_s}{dg}$

burada $V_s$ parçacığın çökelme hızı, $d$ örnekleme tüpünün iç çapı ve $g$ yerçekimi ivmesidir.

Ayrıca bakınız

• Stokes yasası – Reynolds sayısı birden küçük olan parçacıklar üzerindeki akışkanlarda sürükleme kuvvetini ifade eder.

Diğer okumalar

• Şablon:Kitap kaynağı
• Şablon:Kitap kaynağı
• Şablon:Dergi kaynağı
• Şablon:Dergi kaynağı