
Bu çalışma, menteşeyle birbirine bağlanan iki saldan oluşan bir dalga enerjisi dönüştürücüsünün odaklanmış aşırı dalgalar karşısındaki hareketlerini ve yapısal yüklerini, deneysel ölçümlerle doğrulanan üç boyutlu bir hesaplamalı akışkanlar dinamiği modeliyle incelemiştir. Simcenter STAR-CCM+ kullanılarak dalga genliği, tepe frekansı ve iki sal arasındaki mesafe değiştirilmiş; düşey hareket, yunuslama, menteşe kuvvetleri ve bağlama halatı kuvvetleri karşılaştırılmıştır. En kritik dinamik durum, 0,32 metre odaklanmış dalga genliği ve 0,7 Hz tepe frekansında ortaya çıkmış; her iki salın yunuslama açısı yaklaşık 15 dereceye ulaşırken menteşe kuvveti de incelenen koşullar içindeki en yüksek seviyesine çıkmıştır. Bununla birlikte çalışma hakem değerlendirmesinden geçmemiş bir preprinttir; sonuçlar model ölçekli sayısal benzetimlere dayanmakta, gerçek deniz sahasında tam ölçekli yapısal dayanım veya enerji üretimi kanıtı sunmamaktadır.
Sal aralığı, hareket kabiliyeti ile yapısal dayanabilirlik arasında belirgin bir denge oluşturmaktadır. Sal aralığının 0,2 metreye düşürülmesi iki gövdenin daha eşzamanlı hareket etmesini, daha küçük yunuslama tepkileri ve daha düşük bağlantı yükleri oluşmasını sağlamıştır. Ancak bu durum, enerji dönüşümü için gerekli göreli açısal hareketi de sınırlandırmıştır. En büyük yunuslama hareketi 0,4 metre aralıkta görülürken, en yüksek menteşe ve bağlama sistemi yükleri 0,3 metre aralıkta ortaya çıkmıştır. Aralığın 0,6 metreye yükseltilmesi hidrodinamik etkileşimi zayıflatmış, ancak menteşe üzerindeki statik eğilme momentini artırmıştır.
Araştırmanın temel mühendislik mesajı, dalga enerjisi dönüştürücüsünün yalnızca en yüksek enerji yakalama potansiyeline göre tasarlanamayacağıdır. Menteşe bağlantısı, ön ve arka bağlama sistemleri ile sal aralığı, aşırı deniz durumlarında ayrı ayrı değerlendirilmelidir. Özellikle ön ve arka salın bağlama kuvvetlerinin farklı dalga frekanslarında tepe yapması, iki taraf için tek ve ortak bir tasarım yükünün yeterli olmayabileceğini göstermektedir.
Araştırmanın temel sorusu nedir?
Araştırma, menteşeli çift sal tipi bir dalga enerjisi dönüştürücüsünün ani ve yüksek enerjili dalgalar karşısında nasıl hareket ettiğini ve kritik yüklerin hangi koşullarda oluştuğunu belirlemeye çalışmaktadır. İncelenen temel değişkenler odaklanmış dalganın genliği, dalga spektrumunun tepe frekansı ve ön sal ile arka sal arasındaki mesafedir.
Sal tipi dalga enerjisi dönüştürücülerinde elektrik üretme potansiyeli, büyük ölçüde birbirine bağlı gövdelerin göreli açısal hareketine dayanır. Dalgalar iki salı farklı zamanlarda ve farklı yönlerde döndürdüğünde menteşe veya güç alma sistemi üzerinden mekanik enerji elde edilebilir. Ancak göreli hareketin büyümesi, aynı zamanda menteşede daha yüksek kesme kuvvetleri ve eğilme momentleri oluşmasına neden olabilir. Bu nedenle enerji dönüşümünü artıran bir tasarım, aşırı dalgalardaki yapısal güvenlik bakımından en uygun tasarım olmayabilir.
Çalışma neden önemlidir?
Dalga enerjisi dönüştürücüleri normal çalışma koşullarında enerji üretmekle birlikte, hizmet ömürleri boyunca ani ve olağan dışı dalga olaylarına da dayanmak zorundadır. Araştırmacılar, sal tipi dönüştürücüler üzerine yapılan çalışmaların önemli bir bölümünün enerji yakalama verimine odaklandığını; menteşe bağlantılarında ve bağlama noktalarında oluşan aşırı yüklerin ise yeterince incelenmediğini belirtmektedir.
Çalışmanın hedeflediği literatür boşluğu, iki sal arasındaki mesafenin enerji dönüşümü için yararlı göreli hareket ile yapısal yükleri azaltmak için gerekli hareket sınırlaması arasındaki etkisinin yeterince anlaşılmamış olmasıdır. Araştırma bu boşluğu, doğrulanmış bir viskoz akış modeli üzerinden sistematik dalga ve geometri karşılaştırmaları yaparak ele almaktadır.
Odaklanmış aşırı dalga ne anlama gelmektedir?
Odaklanmış dalga, farklı frekanslara ve başlangıç fazlarına sahip çok sayıda dalga bileşeninin belirli bir konumda ve belirli bir zamanda aynı anda tepe oluşturacak biçimde üst üste bindirilmesidir. Bu yöntem, uzun ve rastlantısal bir düzensiz dalga kaydında aşırı bir dalganın tesadüfen oluşmasını beklemek yerine, yüksek enerjili olayın kontrollü biçimde oluşturulmasını sağlar.
Çalışmada kullanılan yaklaşım NewWave temelli odaklanmış dalga düşüncesine dayanmaktadır. Dalga bileşenlerinin fazları, dalga tepesinin menteşe bölgesine yakın önceden tanımlanmış odak noktasında toplanacağı biçimde ayarlanmıştır.
Sayısal model nasıl kurulmuştur?
Sayısal benzetimler Simcenter STAR-CCM+ yazılımında gerçekleştirilmiştir. Ön sal, arka sal ve aralarındaki bağlantı plakası elastik şekil değiştirmeleri ihmal edilen rijit cisimler olarak modellenmiştir. Hesaplama bakımından sistem, dönel mafsallarla birbirine bağlanan üç rijit cisimden oluşmaktadır.
Durgun su seviyesi üzerinde küresel bir koordinat sistemi tanımlanmış; ön ve arka salın ağırlık merkezlerine ayrıca yerel koordinat sistemleri yerleştirilmiştir. Bu koordinat sistemleri her salın hareketini, iki sal arasındaki göreli hareketi, menteşe kuvvetlerini ve bağlama kuvvetlerini ayrı ayrı izlemek için kullanılmıştır.
Akışın temel korunum denklemleri
Modelde türbülanslı su-hava akışı Reynolds ortalamalı Navier-Stokes denklemleriyle çözülmüştür. Sıkıştırılamaz akış için kütlenin korunumu şu biçimde verilmiştir:
\[ \frac{\partial u_i}{\partial x_i}=0 \]
Burada ui, Kartezyen koordinat sistemindeki hız bileşenlerini metre/saniye cinsinden; xi ise metre cinsinden uzaysal koordinatları göstermektedir. Denklem, bir akışkan hücresine giren ve hücreden çıkan hacimsel akışın dengeli olduğunu ifade eder.
Momentum korunumu çalışmada şu genel biçimde tanımlanmıştır:
\[ \frac{\partial \rho u_i}{\partial t}+\frac{\partial \rho u_j u_i}{\partial x_j}-\frac{\partial}{\partial x_j}\left(\mu_{\mathrm{eff}}\frac{\partial u_i}{\partial x_j}\right)=-\frac{\partial p^{*}}{\partial x_i}-g_i x_i\frac{\partial \rho}{\partial x_i}+\sigma\kappa\frac{\partial\alpha}{\partial x_i} \]
- ρ, akışkan yoğunluğunu kilogram/metreküp cinsinden gösterir.
- t, saniye cinsinden zamandır.
- μeff, moleküler ve türbülanslı momentum taşınımını birlikte temsil eden etkin dinamik viskozitedir.
- p*, hidrostatik basınç bileşeni çıkarılmış değiştirilmiş dinamik basınçtır.
- gi, yerçekimi ivmesinin ilgili yön bileşenidir.
- σ, yüzey gerilim katsayısını; κ, su-hava arayüzünün eğriliğini gösterir.
- α, bir hesaplama hücresinin ne kadarının suyla dolu olduğunu belirleyen hacim oranıdır.
Denklem; zamana bağlı ivmeyi, akışla taşınan momentumu, viskoz ve türbülanslı gerilmeleri, basınç kuvvetlerini, yerçekimi-kaldırma etkisini ve su-hava arayüzündeki yüzey gerilimini aynı çözüm içinde birleştirmektedir.
Su-hava yüzeyi nasıl izlenmiştir?
Serbest yüzey, Akışkan Hacmi yöntemiyle (Volume of Fluid, VOF) izlenmiştir. Her hücreye bir α su hacim oranı atanmıştır. α = 0 hücrenin yalnızca hava, α = 1 yalnızca su içerdiğini; 0 ile 1 arasındaki değerler ise su-hava geçiş yüzeyini göstermektedir.
\[ \frac{\partial\alpha}{\partial t}+\frac{\partial(u_i\alpha)}{\partial x_i}+\frac{\partial\left(u_{c,i}\alpha(1-\alpha)\right)}{\partial x_i}=0 \]
Son terim, su-hava arayüzünün sayısal ağ üzerinde gereğinden fazla yayılmasını önleyen yapay sıkıştırma terimidir. Çalışmada sıkıştırma katsayısı ca = 1 olarak kullanılmıştır.
Karışım yoğunluğu ve etkin viskozite, hücredeki su oranına göre hesaplanmıştır:
\[ \rho=\alpha\rho_{\mathrm{water}}+(1-\alpha)\rho_{\mathrm{air}} \]
\[ \mu_{\mathrm{eff}}=\alpha\mu_{\mathrm{water}}+(1-\alpha)\mu_{\mathrm{air}}+\rho\nu_t \]
Burada νt, türbülanslı girdap viskozitesidir. α değeri suya yaklaştıkça karışımın özellikleri suyun; havaya yaklaştıkça havanın fiziksel özelliklerine yaklaşmaktadır.
Serbest yüzey yüksekliği, dikey doğrultudaki su hacim oranlarının toplanmasıyla elde edilmiştir:
\[ \zeta_{\mathrm{water\ level}}=\sum_{i=0}^{n-1}\alpha_i(z_{i+1}-z_i) \]
zi ve zi+1 hücre yüzlerinin metre cinsinden dikey koordinatlarıdır. Bu işlem, her dikey hücre parçasındaki su miktarını toplayarak anlık su yüzeyi yüksekliğini belirlemektedir.
Türbülans nasıl modellenmiştir?
Türbülans kapanışı için Menter’in SST k-ω modeli kullanılmıştır. Model, duvara yakın bölgelerde k-ω yaklaşımından, serbest kayma akışlarında ise k-ε yaklaşımından yararlanır. Araştırmacılar bu modeli hareketli yüzer cisimlerin çevresindeki akış ayrılması, vorteks saçılması ve viskoz etkileri kabul edilebilir hesaplama maliyetiyle temsil edebilmek amacıyla seçmiştir.
Reynolds gerilmesi şu şekilde ifade edilmiştir:
\[ \tau^{R}_{ij}=\frac{2}{3}k\delta_{ij}-\nu_t\left(\frac{\partial u_i}{\partial x_j}+\frac{\partial u_j}{\partial x_i}\right) \]
Burada k, birim kütle başına türbülans kinetik enerjisini; δij, Kronecker delta fonksiyonunu; νt ise türbülanslı girdap viskozitesini göstermektedir.
Çalışmada k ve ω için taşıma denklemleri de verilmiştir. Ancak PDF’deki Denklem 7’de üretim terimi Pk ile dağılım terimi β*ωk arasında bulunması gereken cebirsel işlem işareti görünür değildir. Bu nedenle söz konusu denklem bu sürümde kesin işareti tahmin edilerek yeniden yazılmamıştır. Denklem grubunun genel işlevi; türbülans üretimini, taşınımını, difüzyonunu ve dağılımını hesaplamaktır.
Çok cisimli hareket denklemi
Üç rijit cisimden oluşan sistemin hareketi Newton-Euler denklemlerinin genelleştirilmiş matris biçimiyle ifade edilmiştir:
\[ \mathbf{M}\ddot{\mathbf{c}}=\mathbf{F}_{ef}+\mathbf{F}_{cf} \]
M, üç cismin kütle ve atalet özelliklerini içeren blok köşegen atalet matrisidir. c, sistemin konum ve dönme koordinatlarını; Fef, dalga, yerçekimi, kaldırma ve bağlama gibi dış kuvvetleri; Fcf ise mafsalların kinematik kısıtlarından kaynaklanan bağlantı kuvvetlerini göstermektedir.
\[ \mathbf{M}=\begin{bmatrix}\mathbf{M}_1&0&0\\0&\mathbf{M}_2&0\\0&0&\mathbf{M}_3\end{bmatrix} \]
Bu yapı, ön salın, arka salın ve bağlantı plakasının atalet özelliklerini ayrı bloklarda tutarken mafsal kuvvetleri aracılığıyla hareketlerin birbirini etkilemesini sağlamaktadır.
Odaklanmış dalga denklemleri
Üç boyutlu tek yönlü dalga alanı, farklı frekanslardaki bileşenlerin üst üste bindirilmesiyle oluşturulmuştur:
\[ \eta(x,y,t)=\sum_{i=1}^{N_f}a_i\cos\left(k_i x\cos\theta+k_i y\sin\theta-\omega_i t+\varphi_i\right) \]
- η, metre cinsinden su yüzeyi yüksekliğidir.
- Nf, dalga bileşeni sayısıdır.
- ai, her bileşenin metre cinsinden genliğidir.
- ki, radyan/metre cinsinden dalga sayısıdır.
- θ, dalganın yayılma yönüdür.
- ωi, açısal frekanstır.
- φi, başlangıç fazıdır.
Her bileşenin belirlenen odak konumunda ve odak zamanında aynı anda tepeye ulaşması için şu faz koşulu uygulanmıştır:
\[ \cos\left(k_i x_j\cos\theta+k_i y_j\sin\theta-\omega_i t_j+\varphi_i\right)=1 \]
Toplam odaklanmış dalga genliği bileşen genliklerinin toplamıdır:
\[ A_f=\sum_{i=1}^{N_f}a_i \]
Her bileşenin genliği normalize edilmiş spektral enerji yoğunluğuna göre dağıtılmıştır:
\[ a_i=A_f S_a(f_i)\,df \]
\[ S_a(f_i)=\frac{S(f_i)}{\sum_{n=1}^{N_f}S(f_n)\,df} \]
Formül tutarlılığı notu: PDF’de açısal frekans ilişkisi ωi = 2π/fi biçiminde yazılmıştır. Bu ifade açısal frekansın birimleri bakımından tutarlı değildir. Açısal frekansın standart ilişkisi, çalışmanın mantığını açıklamak için aşağıdaki biçimde düşünülür:
\[ \omega_i=2\pi f_i \]
Bu düzeltme PDF’deki denklem sessizce değiştirilerek yapılmamış; kaynak sürümündeki olası dizgi veya ifade hatası açıkça belirtilmiştir.
Fiziksel deney ve sayısal tank düzeni
Sayısal modelin temel aldığı fiziksel deneyler Zhejiang Ocean University çekme tankında gerçekleştirilmiştir. Tank 130 metre uzunluğunda, 6 metre genişliğinde ve en fazla 4 metre ayarlanabilir su derinliğine sahiptir. Hidrolik dalga üreteci 0,1-2,0 Hz frekans aralığında dalga oluşturabilmektedir.
Sal hareketleri, gövdelerin üzerine yerleştirilen kızılötesi yansıtıcı işaretleri takip eden optik hareket yakalama kameralarıyla ölçülmüştür. Dalga yüksekliği ölçer 50 Hz, hareket yakalama sistemi 120 Hz ve gerilim sensörleri 100 Hz örnekleme frekansıyla çalışmıştır. Dalga yüksekliği kararlı hâle geldikten sonraki 10-15 periyot analiz için seçilmiştir.
Sayısal dalga tankının toplam uzunluğu 28 metre, su derinliği 3 metre ve genişliği 14 metredir. Tankın giriş ve çıkış uçlarında yedişer metre uzunluğunda dalga zorlama bölgeleri, ortada ise 14 metrelik çalışma bölgesi bulunmaktadır. Giriş bölgesi dalgayı üretmek, çıkış bölgesi ise yansıyan dalgaları sönümlemek için kullanılmıştır.
| Deneysel model parametresi | Değer | Açıklama |
|---|---|---|
| Tek salın karakteristik uzunluğu | 1,2 m | Dalga doğrultusundaki gövde uzunluğu |
| Sal genişliği | 0,6 m | Yatay enine boyut |
| Sal yüksekliği | 0,3 m | Gövdenin toplam düşey boyutu |
| Başlangıç sal aralığı | 0,3 m | Temel doğrulama ve H-serisi koşullarındaki mesafe |
| Ağırlık merkezi koordinatları | X = 0, Y = 0, Z = 0,1 m | Yerel koordinat sistemine göre |
| Tek salın kütlesi | 72,4 kg | İç balastla ayarlanan model kütlesi |
Giriş ve çıkış sınırları dalga üretimi ve emilimini destekleyen hız girişi koşullarıyla tanımlanmıştır. Üst yüzey basınç çıkışı, taban duvar, yan yüzeyler simetri düzlemi ve sal yüzeyleri kaymama koşullu duvar olarak modellenmiştir.
Ağ yapısı ve zaman adımı sonuçları etkiliyor mu?
Hareket eden salların çevresinde örtüşen ağ (overset mesh) yöntemi kullanılmıştır. Böylece sal çevresindeki ağ, arka plan akış ağından bağımsız hareket edebilmiş; gövde hareketi sırasında hücrelerin aşırı bozulması önlenmiştir. Su yüzeyi ve sal çevresi yerel olarak inceltilmiştir.
| Ağ çözünürlüğü | Toplam hücre sayısı | Serbest yüzey inceltmesi | Sal çevresi inceltmesi |
|---|---|---|---|
| Kaba | 975.465 | Δz = Ai/2, Δx = Ai/4, Δy = Ai | Δx = Δy = Δz = Ai |
| Orta | 2.682.123 | Δz = Ai/4, Δx = Ai/8, Δy = Ai/2 | Δx = Δy = Δz = Ai/4 |
| İnce | 13.883.546 | Δz = Ai/8, Δx = Ai/16, Δy = Ai/4 | Δx = Δy = Δz = Ai/8 |
Ağ bağımsızlığı testi 0,32 metre hedef genliği ve 0,4 Hz tepe frekansı kullanılarak yapılmıştır. Kaba ağ dalga tepesini düşük tahmin etmiş ve daha büyük faz hatası üretmiştir. Orta ve ince ağ sonuçları birbirine çok yakın bulunmuş; orta ağdaki tepe genliği hatası hedef değerin yüzde 1,5’i içinde kalmıştır. Hesaplama maliyeti nedeniyle sonraki benzetimlerde 2.682.123 hücreli orta ağ seçilmiştir.
Zaman adımı için Tp/100, Tp/1000 ve Tp/10000 seçenekleri karşılaştırılmıştır. Orta zaman adımı olan Δt = Tp/1000, ince çözümle yakın sonuç verdiği için kullanılmıştır.
Model deneysel verilerle nasıl doğrulanmıştır?
Doğrulamanın ilk aşamasında dalga yüksekliği 0,10 metre ve periyodu 1,67 saniye olan düzenli bir dalga kullanılmıştır. Sayısal modelin tahmin ettiği iki sal arasındaki göreli yunuslama hareketi, fiziksel deneyle hem genlik hem de faz bakımından yakın uyum göstermiştir.
İkinci aşamada odaklanmış dalga üretimi, Birleşik Krallık’taki University of Plymouth COAST Laboratory deneyleriyle karşılaştırılmıştır. Bu deney havuzu 35 metre uzunluğunda ve 15,5 metre genişliğindedir. Dalga üretici tarafında su derinliği 4 metre, çalışma bölgesinde 3 metredir ve sistemde geri beslemeli 24 dalga üreteci bulunmaktadır.
CCP-WSI kör karşılaştırmasının 3BT3 Test 4 koşulu kullanılmıştır. Odaklanmış dalga, 0,101563 Hz ile 2 Hz arasında değişen 244 frekans bileşeninin doğrusal üst üste bindirilmesiyle oluşturulmuştur. Kuramsal odak noktası dalga üretim sınırından 14,8 metre uzaktadır. On üç dalga ölçerden sekizinin kayıtları doğrulamada kullanılmıştır. Sayısal model, odak bölgesindeki dalga tepesi, çukur seviyesi ve faz gelişimini deneylerle yakın biçimde yeniden üretmiştir.
Bununla birlikte aşırı dalga koşullarındaki menteşe ve bağlama kuvvetlerinin tamamı, aynı çift sal sistemi üzerinde doğrudan fiziksel yük ölçümleriyle doğrulanmamıştır. Deneysel doğrulama, göreli hareket ve odaklanmış dalga üretiminin güvenilirliğini desteklemekte; bütün yapısal yük tahminlerini bağımsız olarak doğrulamamaktadır.
Dalga genliği arttığında ne olmaktadır?
H1-H4 koşullarında tepe frekansı 0,4 Hz’de sabit tutulmuş, odaklanmış dalga genliği 0,20 metreden 0,32 metreye yükseltilmiştir. Genlik arttıkça hem ön hem de arka salın düşey hareketi ve yunuslaması büyümüştür. Yunuslama tepe değerleri dalga genliğiyle yaklaşık doğrusal artarken, menteşe ve bağlama kuvvetleri doğrusal olmayan biçimde yükselmiştir.
En yüksek genlikli H4 koşulunda, dalga grubu geçtikten sonra arka salın yunuslama hareketinde giderek küçülen salınımlar görülmüştür. Yazarlar bu davranışı keskin gövde kenarlarında akış ayrılması, vorteks saçılması, doğrusal olmayan viskoz sönüm ve akışkan hafızası etkileriyle ilişkilendirmektedir. Bu mekanizma sayısal akış alanına dayalı bir yorumdur; çalışmada vorteks enerjisi veya sönüm bileşenleri ayrı deneylerle ölçülmemiştir.
Tepe frekansı neden 0,7 Hz’de kritik hâle gelmiştir?
H4-H7 koşullarında dalga genliği 0,32 metrede sabit tutulmuş, tepe frekansı 0,4 Hz’den 0,8 Hz’ye yükseltilmiştir. En büyük hareket H6 koşulunda, yani 0,7 Hz’de oluşmuştur. Ön ve arka salın en yüksek yunuslama tepkileri yaklaşık 15 dereceye ulaşmıştır.
Yazarlar 0,7 Hz koşulunu, sal-menteşe-güç alma sistemi-bağlama bütününün küresel doğal frekansına yakın bir uyarım olarak yorumlamaktadır. Bu durumda gövdelerin, menteşenin ve bağlama sisteminin hareketlerinin faz bakımından birbirini güçlendirdiği ve çok cisimli bağlı rezonans oluşturduğu ileri sürülmektedir. Frekans 0,8 Hz’ye yükseldiğinde uyarım sistemin rezonans bölgesinden uzaklaşmış ve yunuslama tepkisi yeniden azalmıştır.
Çalışmada ayrı bir özdeğer veya modal analiz sunulmadığından, 0,7 Hz’in kesin doğal frekans olduğu doğrudan kanıtlanmış değildir. Bulgular, incelenen dalga koşulları içinde 0,7 Hz civarında güçlü bir dinamik büyütme oluştuğunu göstermektedir.
Tepe frekansı yükseldikçe odaklanmış dalga enerjisinin belirlenen konuma ulaşma zamanı da gecikmiştir. Araştırmacılar bunu derin su dalga dispersiyonu ve yüksek frekans ağırlıklı dalga paketlerinin daha düşük grup hızlarıyla ilişkilendirmektedir.
Menteşe ve bağlama kuvvetleri hangi koşullarda en yüksek olmuştur?
Kuvvetler yoğunluk, yerçekimi ivmesi ve karakteristik sal uzunluğu kullanılarak boyutsuzlaştırılmıştır. Boyutsuz menteşe kuvveti dalga genliği ve frekansla doğrusal olmayan biçimde değişmiştir.
| Koşul | Odaklanmış dalga genliği Af | Tepe frekansı fp | Boyutsuz menteşe kuvveti | Boyutsuz ön bağlama kuvveti | Boyutsuz arka bağlama kuvveti |
|---|---|---|---|---|---|
| H1 | 0,20 m | 0,4 Hz | 0,00339 | 0,00128 | 0,00119 |
| H2 | 0,25 m | 0,4 Hz | 0,00432 | 0,00132 | 0,00116 |
| H3 | 0,30 m | 0,4 Hz | 0,00522 | 0,00136 | 0,00110 |
| H4 | 0,32 m | 0,4 Hz | 0,00576 | 0,00156 | 0,00107 |
| H5 | 0,32 m | 0,6 Hz | 0,00704 | 0,00227 | 0,000927 |
| H6 | 0,32 m | 0,7 Hz | 0,00834 | 0,00480 | 0,00102 |
| H7 | 0,32 m | 0,8 Hz | 0,00789 | 0,00414 | 0,000960 |
H4 koşulundaki menteşe kuvveti H1’in yaklaşık 1,70 katı, H6 koşulundaki kuvvet ise H1’in yaklaşık 2,46 katıdır. Çalışma bunu, yüksek genliklerde anlık ıslak yüzeyin ve basınç dağılımının değişmesiyle büyüyen doğrusal olmayan Froude-Krylov kuvvetlerine bağlamaktadır.
En yüksek menteşe kuvveti ve en yüksek ön bağlama kuvveti H6 koşulunda oluşmuştur. Arka bağlama kuvvetinin en yüksek değeri ise H1 koşulundadır. Bu fark, dalga frekansının yük aktarım yolunu değiştirdiğini ve ön ile arka bağlama sisteminin ayrı tasarım yükleriyle değerlendirilmesi gerektiğini göstermektedir.
Metin, arka bağlama kuvvetinin H1’deki tepesini olası bir “bağlı rezonans” ile ilişkilendirirken çalışmanın ana rezonans koşulunu H6 olarak tanımlamaktadır. Ek bir modal veya yük aktarım analizi verilmediğinden, H1’deki arka bağlama tepesinin ayrıntılı fiziksel mekanizması tam olarak açıklığa kavuşmamıştır.
Düşey kuvvet neden baskındır?
Şekil 15’te düşey Z yönündeki boyutsuz dalga kuvveti, X ve Y yönündeki bileşenlerden daha büyüktür. Z yönündeki değer H1’de 0,03099’dan H6’da 0,04781’e yükselmiştir. Araştırmacılar bu baskınlığı, dalga yüzeyinin yükselmesiyle oluşan anlık kaldırma kuvveti artışına bağlamaktadır.
Yatay kuvvetler daha çok dalga kırınımı ve viskoz sürükleme etkileriyle ilişkilendirilmiştir. Bununla birlikte metin Y yönündeki kuvveti “sıfıra yakın” olarak tanımlarken Şekil 15(e)’de verilen boyutsuz Y kuvvetleri yaklaşık 0,00394-0,01007 aralığındadır ve X yönündeki kuvvetlerle aynı büyüklük mertebesindedir. Bu nedenle “sıfıra yakın” ifadesi grafik değerleriyle tam olarak uyumlu görünmemektedir.
Sal aralığı hareketi ve yük aktarımını nasıl değiştirmiştir?
En kritik dalga koşulu olarak belirlenen 0,32 metre genlik ve 0,7 Hz tepe frekansında, sal aralığı 0,2 metreden 0,6 metreye kadar değiştirilmiştir.
| Koşul | Dalga genliği | Tepe frekansı | Sal aralığı | Temel gözlem |
|---|---|---|---|---|
| T1 | 0,32 m | 0,7 Hz | 0,2 m | En eşzamanlı hareket, en küçük yunuslama ve daha düşük yükler |
| T2 | 0,32 m | 0,7 Hz | 0,3 m | En yüksek menteşe, ön bağlama ve arka bağlama kuvvetleri |
| T3 | 0,32 m | 0,7 Hz | 0,4 m | En büyük yunuslama ve göreli hareket |
| T4 | 0,32 m | 0,7 Hz | 0,5 m | Hidrodinamik etkileşimin azalmaya başladığı ara koşul |
| T5 | 0,32 m | 0,7 Hz | 0,6 m | Daha zayıf hareket bağlantısı ve daha yüksek statik eğilme talebi |
0,2 metrelik küçük aralıkta iki sal arasındaki akış etkileşimi “hidrodinamik gölgeleme” oluşturmuş ve gövdeler neredeyse tek bir rijit cisim gibi eşzamanlı hareket etmiştir. Bu durum yunuslama genliğini ve bağlantı yüklerini azaltmış, ancak enerji dönüşümünde kullanılabilecek göreli açısal hareketi sınırlandırmıştır.
0,4 metrelik aralıkta iki sal daha bağımsız hareket etmiş, aralarında belirgin bir faz farkı oluşmuş ve en büyük yunuslama tepkisi görülmüştür. Araştırmacılar bunu daha yüksek enerji yakalama potansiyeliyle ilişkilendirmektedir. Ancak çalışma gerçek PTO gücünü, üretilen elektriği veya yakalama genişliği oranını hesaplamamıştır. Dolayısıyla 0,4 metrenin enerji üretimi açısından kesin olarak en verimli aralık olduğu değil, göreli hareket bakımından en elverişli koşulu oluşturduğu söylenebilir.
0,3 metrelik aralık, en yüksek menteşe ve bağlama kuvvetlerini üretmiştir. İki salın ivmeleri ve hareket fazları arasındaki fark, menteşeye aktarılan dinamik eğilme momenti ile kesme kuvvetini artırmıştır. Bu nedenle T2 koşulu, incelenen geometri içinde menteşe dayanımı için kritik tasarım durumu olarak değerlendirilmiştir.
0,6 metrelik aralıkta gövdeler arasındaki hidrodinamik bağlantı zayıflamış ve dinamik hareket 0,4 metre koşuluna göre azalmıştır. Buna karşılık daha uzun bağlantı mesafesi menteşe üzerindeki statik eğilme momentini artırmıştır. Bu sonuç, gövdeleri birbirinden uzaklaştırmanın dinamik etkileşimi azaltmasına rağmen yapısal kol uzunluğunu ve statik yük talebini büyütebileceğini göstermektedir.
Şekiller ve grafikler ne göstermektedir?
H1-H4 akış alanı görselleri, dalga genliği arttıkça sal yüzeyleri çevresindeki basınç alanının ve gövde eğimlerinin büyüdüğünü göstermektedir. Yunuslama zaman geçmişlerinde ön ve arka salın tepe ve çukur hareketleri benzer zamanlarda oluşmakta, ancak genlik ve salınımın sönme davranışı iki sal arasında farklılaşmaktadır.
H5-H7 görsellerinde 0,7 Hz koşulunda sal gövdelerinin en belirgin karşıt yönlü eğimleri oluşmuştur. 0,8 Hz’de dalga alanı daha kısa periyotlu olmasına rağmen gövdelerin yunuslama genliği azalmıştır. Bu grafikler, sistem tepkisinin dalga frekansıyla sürekli artmadığını ve belirli bir frekans bölgesinde tepe yaptığını göstermektedir.
T1-T5 görselleri, küçük aralıkta sal hareketlerinin birbirine daha yakın olduğunu; orta aralıklarda faz farkının büyüdüğünü; geniş aralıkta ise iki gövdenin hidrodinamik bağının zayıfladığını göstermektedir. Şekil 18’de menteşe ve her iki bağlama kuvvetinin 0,3 metre aralıkta tepe yaptığı, daha küçük veya daha büyük aralıklarda azaldığı görülmektedir.
Grafiklerde deneysel tekrarlar, standart sapmalar, güven aralıkları veya istatistiksel anlamlılık işaretleri verilmemiştir. Sonuçlar deterministik CFD koşulları arasındaki sayısal karşılaştırmalardır; biyolojik veya gözlemsel araştırmalardaki p değeri türünden bir anlamlılık analizi uygulanmamıştır.
Çalışmanın güçlü yönleri nelerdir?
- Model hem çift salın göreli hareketi hem de bağımsız odaklanmış dalga deneyleriyle iki aşamada doğrulanmıştır.
- Ağ ve zaman adımı bağımsızlığı ayrı ayrı test edilmiş ve seçilen çözünürlüğün gerekçesi sayısal olarak verilmiştir.
- Su-hava arayüzü, viskoz etkiler, türbülans ve hareketli gövdeler üç boyutlu CFD modeli içinde birlikte ele alınmıştır.
- Dalga genliği, frekans ve sal aralığı sistematik koşul serileriyle karşılaştırılmıştır.
- Menteşe kuvveti ile ön ve arka bağlama kuvvetleri ayrı değerlendirilmiştir.
- Enerji dönüşümü için gerekli göreli hareket ile aşırı yüklerden kaçınma arasındaki mühendislik ödünleşimi açık biçimde ortaya konmuştur.
Çalışmanın sınırlılıkları nelerdir?
- Çalışma hakem değerlendirmesinden geçmemiş bir preprinttir.
- Sonuçlar model ölçeğindeki sayısal benzetimlere dayanmaktadır; tam ölçekli deniz prototipi doğrulaması yoktur.
- Sal ve menteşe elemanları rijit kabul edilmiştir; elastik şekil değiştirme, yerel gerilme yığılması ve malzeme hasarı modellenmemiştir.
- Menteşe kuvveti hesaplanmış olsa da ayrıntılı gerilme, yorulma ömrü, çatlak ilerlemesi veya nihai kırılma analizi yapılmamıştır.
- Bağlama sisteminin uzun dönem yorulması, gevşek halat darbeleri ve deniz tabanı bağlantıları incelenmemiştir.
- Aşırı deniz durumu kısa süreli odaklanmış dalgalarla temsil edilmiştir; uzun süreli düzensiz fırtına kayıtları değerlendirilmemiştir.
- Yalnızca belirli genlik, frekans ve aralık değerleri incelenmiştir; elde edilen kritik değerlerin bütün sal geometrilerine genellenmesi mümkün değildir.
- Enerji üretimi, PTO gücü ve gerçek elektrik çıktısı doğrudan hesaplanmamıştır.
- 0,7 Hz rezonans yorumu ayrı bir modal veya özfrekans analiziyle doğrulanmamıştır.
- PDF’de açısal frekans formülü ve bazı kuvvet yorumlarında metin-grafik tutarsızlıkları bulunmaktadır.
- Veriler açık bir arşivde sunulmamış, talep üzerine paylaşılacağı belirtilmiştir.
Çalışma neyi desteklemektedir?
- İncelenen çift sal modeli için 0,7 Hz tepe frekansı güçlü hareket ve menteşe yükü büyütmesi oluşturmuştur.
- Dalga genliğinin artması hareket ve yapısal yükleri büyütmüştür.
- Sal aralığı hem göreli hareketi hem de yük aktarım yolunu değiştirmiştir.
- 0,3 metre aralık en yüksek menteşe ve bağlama yüklerini üretmiştir.
- 0,4 metre aralık en büyük yunuslama hareketini üretmiştir.
- 0,2 metre aralık daha eşzamanlı ve daha düşük yüklü bir hareket sağlamıştır.
- Ön ve arka bağlama sistemlerinin kritik frekansları birbirinden farklı olabilmektedir.
Çalışma neyi kanıtlamamaktadır?
- 0,4 metre aralığın gerçek deniz şartlarında en yüksek elektrik üretimini sağlayacağını kanıtlamamaktadır.
- 0,2 metre aralığın bütün sal tipi dönüştürücüler için en güvenli tasarım olduğunu göstermemektedir.
- Modelin tam ölçekli cihazlarda aynı kuvvet ve hareket değerlerini vereceğini kanıtlamamaktadır.
- Menteşenin belirli bir malzeme veya kesitle yorulmadan hizmet vereceğini göstermemektedir.
- Sistem için uzun dönem fırtına dayanımı veya çalışma ömrü tahmini sunmamaktadır.
- 0,7 Hz’in bütün benzer cihazların doğal frekansı olduğunu göstermemektedir.
- Endüstriyel veya ticari ölçekte başarılı enerji üretimi kanıtı sunmamaktadır.
Geçmiş, bugün ve gelecek açısından önemi
Geçmişte sal tipi dalga enerjisi çalışmalarının önemli bir bölümü normal dalga koşullarındaki hareket ve enerji yakalama performansına yoğunlaşmıştır. Bu araştırma, aynı göreli hareketin aşırı dalgalarda menteşe ve bağlama sistemi için kritik bir yük kaynağına dönüşebileceğini öne çıkarmaktadır.
Bugünkü mühendislik uygulamaları açısından çalışma, cihaz geometrisinin yalnızca enerji üretim hedefiyle değil, hayatta kalabilirlik ve yapısal yük yollarıyla birlikte değerlendirilmesi gerektiğini göstermektedir. Özellikle ön ve arka bağlama hatlarının ayrı yük zarflarıyla tasarlanması ve 0,7 Hz civarındaki dinamik büyütmenin menteşe dayanım hesaplarına dâhil edilmesi önerilmektedir.
Gelecekte aynı yaklaşımın tam ölçekli prototip deneyleri, uzun süreli düzensiz fırtına dalgaları, elastik gövde ve menteşe modelleri, yorulma analizi, ayrıntılı PTO dinamiği ve gerçek enerji üretim ölçümleriyle genişletilmesi gerekmektedir. Böyle bir doğrulama, enerji yakalama potansiyeli ile yapısal güvenlik arasındaki optimum tasarım noktasının daha güvenilir biçimde belirlenmesini sağlayabilir.
Çalışmanın Yöntemi ve Bulguları
Teknik yöntem özeti
| Yöntem bileşeni | Çalışmada kullanılan yaklaşım |
|---|---|
| Çalışma türü | Deneysel verilerle doğrulanan üç boyutlu CFD ve çok cisimli dinamik benzetim |
| Yazılım | Simcenter STAR-CCM+ |
| Akış çözücüsü | Sonlu hacimler yöntemi, zamana bağlı örtük çözüm |
| Akış denklemleri | Reynolds ortalamalı Navier-Stokes denklemleri |
| Türbülans modeli | Menter SST k-ω |
| Serbest yüzey | VOF su-hava hacim oranı yöntemi |
| Hareketli gövde ağı | Overset mesh |
| Cisim modeli | Ön sal, arka sal ve bağlantı plakasından oluşan üç rijit cisim |
| Bağlantı | Kinematik kısıtlı dönel mafsallar |
| Sayısal tank | 28 m uzunluk, 14 m genişlik, 3 m su derinliği |
| Dalga bölgeleri | Girişte ve çıkışta 7 m dalga zorlama/sönümleme bölgesi |
| Seçilen ağ | 2.682.123 hücreli orta çözünürlük |
| Zaman adımı | Δt = Tp/1000 |
| Veri işleme | Dalga yüksekliği, yunuslama zaman geçmişi ve boyutsuzlaştırılmış kuvvet tepeleri |
| İstatistiksel test | Uygulanmamıştır; deterministik koşullar arası CFD karşılaştırması yapılmıştır |
Doğrulama özeti
| Doğrulama aşaması | Koşul | Değerlendirilen ölçüt | Sonuç |
|---|---|---|---|
| Çift sal hareketi | H = 0,10 m, T = 1,67 s düzenli dalga | İki sal arasındaki göreli yunuslama | Deney ve benzetim genlik ve faz bakımından yakın uyum göstermiştir |
| Odaklanmış dalga üretimi | 244 bileşen, 0,101563-2 Hz, odak noktası 14,8 m | Dalga tepesi, çukur seviyesi ve faz gelişimi | Sekiz dalga ölçer konumunda deneyle yakın uyum bildirilmiştir |
| Ağ bağımsızlığı | Af = 0,32 m, fp = 0,4 Hz | Odaklanmış dalga yüksekliği | Orta ağ hatası yüzde 1,5 içinde kalmıştır |
| Zaman adımı bağımsızlığı | Tp/100, Tp/1000 ve Tp/10000 | Dalga yüksekliği ve fazı | Tp/1000 yeterli kabul edilmiştir |
Ana koşullar ve bulgular
| Değişken | İncelenen aralık | Kritik sonuç | Yorum sınırı |
|---|---|---|---|
| Dalga genliği | 0,20-0,32 m; fp = 0,4 Hz | Hareketler ve yükler büyümüş; kuvvet artışı doğrusal olmayan özellik göstermiştir | Yalnızca dört genlik değeri incelenmiştir |
| Tepe frekansı | 0,4-0,8 Hz; Af = 0,32 m | 0,7 Hz’de yaklaşık 15 derece yunuslama ve en yüksek menteşe kuvveti oluşmuştur | Rezonans yorumu bağımsız modal analizle doğrulanmamıştır |
| Sal aralığı | 0,2-0,6 m; Af = 0,32 m, fp = 0,7 Hz | 0,3 m en yüksek bağlantı yükünü, 0,4 m en büyük yunuslamayı üretmiştir | Enerji üretim gücü doğrudan hesaplanmamıştır |
| Ön bağlama sistemi | H1-H7 | En yüksek yük H6 koşulunda oluşmuştur | Boyutsuz tepe kuvveti raporlanmıştır |
| Arka bağlama sistemi | H1-H7 | En yüksek yük H1 koşulunda oluşmuştur | H1 yük tepesinin mekanizması tam açıklanmamıştır |
| Düşey dalga kuvveti | H1-H7 | Z yönü en büyük kuvvet bileşenidir; en yüksek değer H6’dadır | Metindeki Y yönü yorumu grafikle tam uyumlu değildir |
Dayanabilirlik odaklı mühendislik çıkarımı
Çalışmanın sonuçları tek bir “en uygun” sal aralığı tanımlamamaktadır. 0,2 metre aralık yük azaltma bakımından avantajlı, ancak göreli hareket bakımından sınırlayıcıdır. 0,4 metre aralık göreli hareketi büyütmekte, fakat enerji üretimi doğrudan ölçülmediğinden yalnızca enerji yakalama potansiyeli açısından elverişli kabul edilmektedir. 0,3 metre aralık ise menteşe ve bağlama sistemleri açısından kritik yük durumunu oluşturmaktadır. Tasarımın, normal çalışma enerji performansı ile aşırı dalga dayanımı arasında çok amaçlı bir optimizasyon gerektirdiği anlaşılmaktadır.
Kaynak ve Yöntem Notu
Çalışmanın tam özgün adı: Hydrodynamic response and survivability of a hinged dual-raft wave energy converter under focused extreme waves
Yazarlar, PDF’deki sıralarıyla: Peng Xu; Cheng Lu; Yuan Zhang; Changqing Jiang; Xuanyu Chen; Xiwu Gong.
Eş birinci yazar: PDF’de eş katkı veya eş birinci yazarlık bilgisi belirtilmemiştir.
Sorumlu yazar: Yuan Zhang, School of Naval Architecture and Maritime, Zhejiang Ocean University. PDF’de sorumlu yazar e-posta adresi zhangyuan@zjou.edu.cn olarak verilmiştir.
Kurumsal bağlantılar:
- School of Naval Architecture and Maritime, Zhejiang Ocean University, Zhoushan 316022, China.
- Guangxi Key Laboratory of Ocean Engineering Equipment and Technology, Beibu Gulf University, Qinzhou 535011, China.
- Institute for Sustainable and Autonomous Maritime Systems, University of Duisburg-Essen, Duisburg 47057, Germany.
DOI: PDF’nin ön kaynak kimliği sayfasında DOI basılı değildir. Resmî SSRN kaydında DOI 10.2139/ssrn.6947528 olarak yer almaktadır.
Dergi: Hakemli bir dergi sürümü doğrulanmamıştır.
Yayın platformu: SSRN.
Özgün yayınevi: Çalışma bir hakemli dergi makalesi olarak yayımlanmadığından geleneksel anlamda özgün bir dergi yayınevi bulunmamaktadır. İncelenen sürüm SSRN üzerinde sunulan preprinttir.
Yayın/kayıt yılı: 2026. PDF içinde kesin yükleme veya yayımlanma günü belirtilmemiştir.
Hakemlik durumu: Bu çalışma bir preprinttir ve hakem değerlendirmesinden geçmemiştir. PDF’nin her sayfasında “This preprint research paper has not been peer reviewed” uyarısı bulunmaktadır.
Kaynak türü: Deneysel doğrulama içeren üç boyutlu CFD tabanlı sayısal araştırma makalesi ve preprint.
Resmî bağlantı:SSRN kayıt sayfası, özet numarası 6947528.
Finansman: Çalışma, National Natural Science Foundation of China tarafından 52301342, 52271293 ve 52401345 numaralı hibelerle; Natural Science Foundation of Guangxi Zhuang Autonomous Region tarafından 2026GXNSFAA00641300 numaralı hibeyle desteklenmiştir.
Çıkar çatışması: Yazarlar, çalışmayı etkileyebilecek bilinen bir mali çıkar veya kişisel ilişki bulunmadığını beyan etmiştir.
Veri erişimi: Verilerin talep üzerine paylaşılacağı belirtilmiştir.
İçerik hazırlama yöntemi: Bu Türkçe açıklama yüklenen PDF’nin metni, denklemleri, tabloları, hareket grafikleri, kuvvet grafikleri, deney düzeneği görselleri ve sayısal akış görselleri birlikte incelenerek hazırlanmıştır. PDF dışından bilimsel bulgu eklenmemiştir. Dış doğrulama yalnızca DOI ve resmî SSRN kaynak kimliği için kullanılmıştır.
Temel sınırlılıklar: Çalışma hakemlikten geçmemiştir; model ölçeğindedir; gövdeler rijit kabul edilmiştir; tam ölçekli saha doğrulaması, yapısal gerilme ve yorulma analizi, uzun süreli düzensiz fırtına koşulları ve doğrudan PTO enerji üretimi değerlendirmesi bulunmamaktadır. 0,7 Hz’deki bağlı rezonans yorumu ayrı bir modal analizle doğrulanmamıştır. PDF’de açısal frekans denklemi ile Y yönü kuvvet yorumunda açıkça belirtilmesi gereken tutarsızlıklar bulunmaktadır.
Bilimsel uyarı: Bu çalışma hakem değerlendirmesinden geçmemiş bir preprinttir. Sonuçlar sayısal model ve model ölçekli deneysel doğrulama kapsamındadır; gerçek deniz sahasında tam ölçekli dayanım, enerji üretimi veya uzun dönem işletme başarısı kanıtı olarak yorumlanmamalıdır.

Bir yorum bırakın
E-posta adresiniz yayınlanmayacaktır. Gerekli alanlar * ile işaretlenmiştir