Bu çalışma, nükleer reaktörlerin borular, pompalar, bağlantılar, ısı kaynakları ve hidrolik kayıplardan oluşan bütün devrelerini OpenFOAM ortamında bir boyutlu olarak modellemek için yeni bir termal-hidrolik çerçeve geliştirmektedir. Yöntem, foamForNuclear platformunun mevcut tek fazlı kütle, momentum ve enerji denklemlerini değiştirmek yerine; bir boyutlu bileşenleri kütle, momentum ve enerjiyi koruyarak birbirine bağlayan yeni sınır koşulları, çok kollu bağlantı modelleri, basınç sıçraması hesapları, dört bölgeli pompa modeli ve otomatik ağ oluşturma araçları eklemektedir. Analitik doğrulama problemlerinde akış, basınç ve sıcaklık sonuçları %1’in altında bağıl hatalar vermiş; basitleştirilmiş iki çevrimli bir reaktör devresinde TRACE ile yapılan karşılaştırmada kararlı durum çekirdek debileri yaklaşık %3 içinde uyuşmuştur. Bununla birlikte çalışma yalnızca tek fazlı akışları kapsamaktadır ve deneysel bir reaktör sistemiyle değil, analitik çözümler ve başka bir sistem koduyla karşılaştırılmıştır.
Yeni çerçevenin temel amacı, düşük hesaplama maliyetli bir boyutlu sistem modelleriyle yüksek çözünürlüklü CFD veya gözenekli ortam modellerini aynı OpenFOAM altyapısında birleştirebilecek ortak bir temel oluşturmaktır. Geliştirilen branchConnector1D sınır koşulu, bağlı boruların kesitleri farklı olduğunda hızları eşitlemek yerine kütlesel debiyi korumaktadır. pressureCoupling1D ise dinamik basınç, yerel hidrolik kayıp, pompa basıncı ve hidrostatik farkları tek bir basınç sıçraması denkleminde birleştirmektedir. Dört bölgeli homolog pompa modeli, pompa durması sırasında debinin sıfıra yaklaşmasını, tersine dönmesini ve pompanın ataletle yavaşlamasını temsil edebilmektedir.
Türkiye açısından değerlendirme: Çalışma; Türkiye’de nükleer enerji mühendisliği, reaktör güvenlik analizi, araştırma reaktörleri, küçük modüler reaktör çalışmaları, enerji sistemleri eğitimi ve açık kaynaklı mühendislik yazılımı geliştirme faaliyetleri açısından yöntemsel değer taşımaktadır. Yaklaşımın Türkiye’deki araştırma ve tasarım projelerinde kullanılabilmesi için hedef reaktörün gerçek geometri, malzeme, pompa, vana, ısı değiştirici ve kontrol verileriyle yeniden kurulması; deneysel devreler ve uluslararası karşılaştırma problemleriyle doğrulanması; yazılım kalite güvencesi, sürüm kontrolü ve belirsizlik analizlerinin tamamlanması gerekir. Bu araştırmadan geliştirilen çerçevenin lisanslama analizlerine hazır olduğu, iki fazlı kazaları modelleyebildiği, TRACE veya RELAP5’in yerini doğrudan alabileceği ya da gerçek bir nükleer tesisin güvenliğini kanıtladığı sonucu çıkarılamaz.
Neden bir boyutlu reaktör sistem kodlarına ihtiyaç duyulmaktadır?
Nükleer santrallerin termal-hidrolik davranışı; reaktör çekirdeği, sıcak ve soğuk kollar, pompalar, buhar üreteçleri, basınç kapları, bağlantı boruları ve güvenlik sistemlerinin birlikte çalışmasıyla belirlenmektedir. Bütün tesisin yüksek çözünürlüklü üç boyutlu CFD ile modellenmesi, rutin güvenlik analizleri için çok yüksek hesaplama maliyeti oluşturabilir.
Bir boyutlu sistem kodları, boru ve ekipmanları eksen boyunca sıralanmış kontrol hacimleri olarak temsil eder. Böylece basınç, debi, sıcaklık ve enerji dağılımı bütün reaktör devresinde görece düşük hesaplama maliyetiyle izlenebilir. TRACE, RELAP5 ve ATHLET gibi araçlar bu yaklaşımın yerleşik örnekleridir.
Ancak yeni reaktör tasarımlarında doğal dolaşım, büyük havuzlar, termal tabakalaşma ve karmaşık karışım bölgeleri gibi üç boyutlu olaylar da önem kazanmıştır. Bu nedenle tek bir çözünürlük düzeyi yerine, bazı bölgelerin bir boyutlu sistem modeliyle, bazı bölgelerin gözenekli ortam yaklaşımıyla, kritik bölgelerin ise CFD ile çözüldüğü çok doğruluklu yöntemler geliştirilmektedir.
foamForNuclear içindeki eksiklik neydi?
foamForNuclear, OpenFOAM tabanlı çok fizikli ve çok doğruluklu bir nükleer analiz platformudur. Platform, GeN-Foam ile OFFBEAT çözücülerinin geliştirilmesi ve birleştirilmesiyle; akış, ısı transferi, nötronik ve yakıt davranışı gibi farklı fiziklerin ortak bir ortamda çalıştırılmasını amaçlamaktadır.
Çalışmaya göre foamForNuclear, CFD ve gözenekli ortam tabanlı termal-hidrolik bölgeleri aynı simülasyon içinde çalıştırabilmesine rağmen bütün santral devresini kurmak için gereken yerel bir bir boyutlu sistem altyapısına sahip değildi. Pompa, vana, dirsek, T-bağlantı, çap değişimi ve yerel kayıp gibi bileşenler için haricî bir sistem koduna veya OpenModelica bağlantısına ihtiyaç duyuluyordu.
Araştırmacılar bu bağımlılığı kaldırmak için OpenFOAM içinde çalışan yerel bir çerçeve geliştirmiştir. Çalışmanın “OpenFOAM içine doğrudan gömülmüş ilk uygulama” olduğu yönündeki ifade araştırmacıların kendi öncelik iddiasıdır; çalışma içinde kapsamlı ve bağımsız bir yazılım önceliği incelemesi sunulmamıştır.
Geliştirilen çerçeve neyi değiştirmektedir?
Çalışma, foamForNuclear’ın mevcut tek fazlı gözenekli ortam termal-hidrolik denklemlerini yeniden yazmamaktadır. Kütle, momentum ve entalpi korunumunu çözen alan denklemleri aynı kalmaktadır. Yeni çalışma esas olarak şu altyapıyı eklemektedir:
- Bir boyutlu boruları farklı yön ve kesitlerde birbirine bağlayan sınır koşulları,
- Kesit değişimlerinde kütlesel debiyi koruyan hız eşleştirmesi,
- Üç veya daha fazla kolun birleştiği bağlantılar,
- Dinamik basınç, yerel kayıp, pompa ve hidrostatik etkileri içeren basınç sıçraması,
- İleri ve geri akışı kapsayan pompa karakteristikleri,
- Pompa durması ve ataletle yavaşlama modeli,
- Bağlantılarda kütlesel debi ağırlıklı entalpi karışımı,
- Python üzerinden otomatik bir boyutlu boru ağı oluşturma araçları.
OpenFOAM’ın mevcut bir boyutlu yaklaşımı neden yeterli değildi?
OpenFOAM, ince bir üç boyutlu ağın bazı yüzeylerinde empty sınır koşulunu kullanarak iki veya bir boyutlu çözümler oluşturabilmektedir. Ancak empty koşulu, çözümün yalnızca ilgili yüzeye teğet doğrultularda değişmesine izin verdiği için yön değiştiren boru hatlarında sorun oluşturmaktadır. Örneğin 90° dirsek içeren bir devre, tek doğrultulu düz bir ağ gibi temsil edilemez.
Bir çözüm olarak boruların dış yüzeylerinde kaymalı duvar, bağlantılarda ise cyclic veya uyumsuz ağlar için cyclicAMI kullanılabilir. Bu yöntem boruları farklı açılarla bağlayabilse de çalışma üç temel sorun belirlemiştir:
- Kütle korunumu: Kesit alanları farklı iki borunun sınırında aynı hızın uygulanması, aynı kütlesel debinin geçtiğini garanti etmez.
- Bütün alanlara aynı sınır tipinin uygulanması: Basınç, hız ve entalpi bağlı olmalıdır; ancak hidrolik çap gibi yardımcı bir alan sınır boyunca sürekli olmak zorunda değildir.
- Karmaşık ağ oluşturma: Her bağlantı için öteleme vektörü, dönme açısı ve dönme merkezinin elle tanımlanması büyük devrelerde hata riskini artırmaktadır.
branchConnector1D ne yapmaktadır?
Araştırmacılar farklı alan değişkenlerine farklı davranışlar atayabilmek için branchConnector1D adlı bir bağlantı sınır koşulu geliştirmiştir. Basınç, hız, entalpi ve seçilen türbülans değişkenleri gibi fiziksel olarak bağlanması gereken alanlar çevrimsel bağlantıyla eşleştirilirken, hidrolik çap gibi bağlantı boyunca farklı kalabilen alanlarda sıfır gradyan uygulanmaktadır.
Bu ayrım, yalnızca korunması veya sürekli olması gereken değişkenlerin karşı boruyla güçlü biçimde bağlanmasını sağlamaktadır. Böylece bütün alanların zorunlu olarak aynı değere getirildiği yapay bir bağlantı önlenmektedir.
Farklı çaptaki borularda kütle nasıl korunmaktadır?
İki borunun kesitleri, yoğunlukları veya etkin akış alanları farklı olduğunda hızların eşitlenmesi yerine kütlesel debilerin eşitlenmesi gerekir. Çalışmanın kullandığı temel ilişki şöyledir:
\[ \alpha_P \rho_P A_P U_P = \alpha_N \rho_N A_N U_N \]
Burada \(\alpha\) fazın kontrol hacminde kapladığı oranı, \(\rho\) yoğunluğu, \(A\) kesit alanını ve \(U\) hızı göstermektedir. \(P\) ve \(N\), bağlantının iki tarafındaki kontrol hacimlerini ifade etmektedir.
Komşu taraftan gelen hız katkısı kesit, yoğunluk ve faz oranı kullanılarak ölçeklendirilmektedir. Böylece küçük kesitli boruda hız yükselirken büyük kesitli boruda düşebilmekte, fakat bağlantıdan geçen toplam kütlesel debi aynı kalmaktadır.
Çok kollu bağlantılar nasıl modellenmiştir?
Standart çevrimsel OpenFOAM bağlantısı iki yüzey arasındaki ikili ilişki için tasarlanmıştır. T-bağlantıda ise bir ana boru iki kola ayrılabilir veya iki kol tek boruda birleşebilir. Çalışmanın 7. sayfasındaki Şekil 1, ana borunun ucunda iki ayrı bağlantı yüzeyi oluşturularak üst ve alt kollara bağlandığı düzeni göstermektedir.
Üç kollu bağlantıdaki kütle korunumu şu şekilde yazılmıştır:
\[ \alpha_a\rho_aA_aU_a+ \alpha_b\rho_bA_bU_b+ \alpha_c\rho_cA_cU_c=0 \]
Hızlar bağlantıya doğru olduğunda pozitif, bağlantıdan uzaklaştığında negatif kabul edilmektedir. Bir sınır yüzeyinin etkili komşu hızı, yalnızca doğrudan eşleştirildiği yüzeyden değil, bağlantıdaki diğer bütün kollardan gelen kütlesel debilerden de hesaplanmaktadır.
Diğer kollardan gelen katkı açık bir kaynak terimi olarak eklendiği için formülasyon tamamen değil, yarı örtüktür. Dengesiz debi dağılımlarında yakınsamayı iyileştirmek amacıyla bu açık terime isteğe bağlı gevşetme uygulanabilmektedir.
Bir hücrenin ucunda birden fazla bağlantı yüzeyi oluşturulması geometrik hücre hacmini bozabileceğinden, araştırmacılar bağlantıya komşu hücrelerin hacimlerini ve basınç-hız bağlaşımındaki yüzey katkılarını yeniden hesaplamıştır.
Basınç sıçraması hangi etkileri içeriyor?
Bağlantıdan geçen kütlesel debinin korunması, bağlantının iki tarafındaki basınçların mutlaka eşit olması gerektiği anlamına gelmez. Kesit değişimi, dirsek, vana, pompa ve yükseklik farkı basınçta fiziksel bir sıçrama meydana getirebilir.
Bu amaçla pressureCoupling1D adlı bir sınır koşulu geliştirilmiştir. Kaldırma kuvveti içeren OpenFOAM çözücülerinde taşınan değiştirilmiş basınç:
\[ p_{rgh}=p-\rho gz \]
olarak tanımlanmaktadır. İki taraf arasındaki basınç sıçraması çalışmada şu bileşenlerle ifade edilmiştir:
\[ p_{rgh,1}-p_{rgh,2} = \frac{1}{2}\rho\left(U_1^2-U_2^2\right) -\frac{1}{2}K\rho_uU_u^2 +S_p +gz\left(\rho_2-\rho_1\right) \]
Birinci terim dinamik basınç farkını, ikinci terim yerel hidrolik kaybı, \(S_p\) pompanın oluşturduğu basınç artışını ve son terim yoğunluk ile yükseklikten kaynaklanan hidrostatik düzeltmeyi temsil etmektedir.
Yerel kayıp katsayısı \(K\) ve pompa kaynağı zamanla değişen fonksiyonlar olarak verilebildiğinden vana hareketi veya pompa durması gibi geçici olaylar modellenebilmektedir. Akış yön değiştirdiğinde kayıp terimi yeni giriş tarafına göre yeniden yönlendirilmektedir.
Pompa için hangi iki model geliştirilmiştir?
İlk model, pompa yüksekliğini hacimsel debinin polinomu olarak tanımlamaktadır:
\[ h(Q)=\sum_i c_iQ^i \]
Pompa hızı değiştiğinde katsayılar benzerlik yasalarına göre ölçeklendirilmektedir. Bu model tasarım çalışma noktasına yakın ileri akış koşullarında sade ve kullanışlıdır. Ancak \(h(Q)\) tek değerli olduğundan sıfır debi, ters debi veya ters dönme gibi durumları eksiksiz temsil edemez.
İkinci model, sistem kodlarında kullanılan homolog veya dört bölgeli pompa yaklaşımıdır. Debi, dönme hızı ve basma yüksekliği nominal değerlere göre boyutsuzlaştırılmaktadır:
\[ q^*=\frac{Q}{Q_R}, \qquad \omega^*=\frac{\omega}{\omega_R}, \qquad h^*=\frac{h}{H_R} \]
İleri ve ters debi ile ileri ve ters dönme koşulları dört çalışma bölgesine ayrılmaktadır. Her bölgede pompa yüksekliği ve hidrolik tork, parçalı doğrusal tablolarla hesaplanmaktadır. Bu düzen, geleneksel sekiz oktantlı pompa karakteristiklerini benzer normalleştirmeye sahip oktantları birleştirerek dört bölgede temsil etmektedir.
Pompa durması ve ataletle yavaşlama nasıl hesaplanmaktadır?
Pompa hızı kullanıcı tarafından zaman çizelgesiyle verilebilir veya pompa durduktan sonra rotor dinamiğiyle hesaplanabilir. Motor torku kesildiğinde dönme hızı şu dengeyle değişmektedir:
\[ I\frac{d\omega}{dt} = -T_{\mathrm{hy}}\left(q^*,\omega^*\right) -T_f(\omega) \]
Burada \(I\) rotor atalet momentini, \(T_{\mathrm{hy}}\) akışkanın oluşturduğu hidrolik torku ve \(T_f\) mekanik sürtünme torkunu ifade etmektedir.
Dört bölgeli pompa karakteristiğiyle rotor denklemi birlikte çözüldüğünde pompa normal çalışma noktasından düşük debiye, sıfır debiye ve ters akış bölgesine kesintisiz geçebilmektedir. Bu özellik, pompa kaybı ve akış yönünün değiştiği reaktör geçişleri için önemlidir.
Bağlantılarda enerji korunumu nasıl sağlanmıştır?
Birden fazla kolun birleştiği noktada çıkış entalpisi tek bir komşu borunun entalpisine eşitlenemez. Çıkış akışkanı, bağlantıya giren bütün akımların karışımıdır. Çalışmada giriş entalpisi kütlesel debi ağırlıklı ortalamayla hesaplanmıştır:
\[ h_{\mathrm{mix}} = \frac{\sum_i\dot{m}_ih_i} {\sum_i\dot{m}_i} \]
\[ \dot{m}_i = \max\left(\rho_i\alpha_iA_iU_i,0\right) \]
Toplama yalnızca bağlantıya doğru akan kollarda yapılmaktadır. Bir kol bağlantıdan dışarı akıyorsa entalpi sıfır gradyan davranışıyla taşınmaktadır. Akış yönü geçici olay sırasında değiştiğinde karışıma katılan kollar da dinamik olarak yeniden belirlenmektedir.
Python arayüzü devre kurulumunu nasıl kolaylaştırmaktadır?
foamForNuclear Python API’si, bir boyutlu boru parçalarının başlangıç konumu, yönü, uzunluğu, hidrolik çapı, hücre sayısı ve dirsek yarıçapıyla tanımlanmasına olanak verecek biçimde genişletilmiştir.
Yeni bir borunun başlangıcı bir koordinat vektörü yerine mevcut bir boru nesnesi olarak verilebilmektedir. Bu durumda yeni boru önceki borunun çıkışına otomatik olarak bağlanmakta; gerekli dönme, öteleme ve sınır koşulları yazılım tarafından oluşturulmaktadır.
Bağımsız borular için giriş ve çıkış yüzeyleri, bağlı borular için ise branchConnector1D bağlantıları otomatik üretilmektedir. Böylece büyük bir devrenin her bağlantısı için dönüş merkezi ve geometrik dönüşümün elle hesaplanması gerekmemektedir.
Birinci analitik doğrulama: akış bölünmesi ve yeniden birleşme
İlk doğrulama probleminde giriş debisi 570 L·s−1 olan su akışı, çapları, uzunlukları ve sürtünme katsayıları farklı iki paralel kola ayrılmış ve daha sonra yeniden birleşmiştir. İki T-bağlantı arasında 9 m yükseklik farkı bulunmaktadır.
1090 hücreli ağda elde edilen sonuçlar şöyledir:
| Değişken | Analitik çözüm | FFN sonucu | Bağıl hata |
|---|---|---|---|
| Birinci kol debisi | 121,8 L·s−1 | 122,2 L·s−1 | %0,33 |
| İkinci kol debisi | 448,2 L·s−1 | 447,8 L·s−1 | %0,089 |
| A–B basınç farkı | 147 kPa | 147,15 kPa | %0,10 |
Bu sonuçlar, farklı kollara dağılan kütlesel debinin ve dağıtılmış sürtünme kayıplarının analitik çözümle uyumlu biçimde hesaplandığını göstermektedir.
İkinci doğrulama: pompalı ve doğal dolaşımlı kapalı çevrim
İkinci problem, çapları 80, 90, 100 ve 110 mm olan dört borudan meydana gelen kapalı bir çevrimdir. Her boru 40 m uzunluğundadır. İlk düzende akış bir pompayla, ikinci düzende ise alt borunun 350 K’ye ısıtılması ve üst borunun 300 K’ye soğutulmasıyla oluşan kaldırma kuvvetiyle sağlanmıştır.
| Durum | Analitik debi | FFN debisi | Bildirilen hata |
|---|---|---|---|
| Pompalı çevrim | 0,0976 m3·s−1 | 0,0976 m3·s−1 | <%0,05 |
| Doğal dolaşımlı çevrim | 0,7345 m3·s−1 | 0,7341 m3·s−1 | <%0,05 |
Kaynak, analitik ve sayısal sonuçların birbirine çok yakın olduğunu göstermektedir. Bununla birlikte tabloda m3·s−1 olarak verilen debilerin büyüklüğü; belirtilen boru çapları, uzunlukları ve pompa eğrisiyle birlikte değerlendirildiğinde birim veya ölçek açıklaması gerektirmektedir. Örneğin kaynakta verilen \(H_{\mathrm{pump}}=0{,}15-15Q^2\) eşitliğinde \(Q=0{,}0976\) kullanıldığında pompa yüksekliği yalnızca yaklaşık 0,0071 m olmaktadır. Bu nedenle Tablo 3’teki debi biriminin veya ondalık gösteriminin yazarlar tarafından kontrol edilmesi gerekir.
Üçüncü doğrulama: çoklu kollar ve enerji karışımı
Son analitik problem, art arda iki akış bölünmesi, iki yeniden birleşme ve 30 MW’lık yerel ısı kaynağı içeren daha karmaşık bir boru ağıdır. Giriş hızı 1 m·s−1, dirsek kayıp katsayısı 0,9 ve T-bağlantı kayıp katsayısı 1,2 olarak uygulanmıştır.
Çalışmanın 18. sayfasındaki Şekil 5, ısıtılan üst kolun sıcaklığının yükseldiğini ve karışım bağlantılarından sonra sıcaklığın diğer kollara dağıldığını göstermektedir. 19. sayfadaki Şekil 6 ise paralel kollardaki farklı hız büyüklüklerini görselleştirmektedir.
| Değişken | Analitik çözüm | FFN sonucu | Bağıl hata |
|---|---|---|---|
| Alt kol debisi \(Q_b\) | 0,4367 m3·s−1 | 0,4369 m3·s−1 | %0,06 |
| Isıtılan kol debisi \(Q_h\) | 0,1442 m3·s−1 | 0,1444 m3·s−1 | %0,16 |
| Isıtılmayan üst alt-kol debisi \(Q_l\) | 0,2045 m3·s−1 | 0,2047 m3·s−1 | %0,08 |
| Çıkış sıcaklığı | 309,11 K | 309,20 K | %0,03 |
Bu karşılaştırma, yalnızca toplam debinin değil, paralel kollardaki dağılımın ve farklı sıcaklıklardaki akımların birleşiminden oluşan çıkış entalpisinin de analitik sonuçla uyumlu olduğunu göstermektedir.
TRACE karşılaştırması için nasıl bir reaktör devresi kurulmuştur?
Sistem düzeyindeki değerlendirmede, basınçlı su reaktörüne benzeyen basitleştirilmiş iki çevrimli birincil devre kullanılmıştır. Modelde:
- Reaktör çekirdeği,
- Alt ve üst plenumlar,
- İki simetrik sıcak kol,
- İki buhar üreteci temsili,
- İki ana soğutucu pompası,
- İki soğuk kol,
- Plenumlar arasında çapraz bağlantı
bulunmaktadır. Çalışmanın 20. sayfasındaki Şekil 7, TRACE düğümlemesiyle OpenFOAM/FFN boru ağını yan yana göstermektedir.
Dağıtılmış sürtünme Churchill bağıntısıyla, çekirdek ve buhar üreteçlerindeki taşınım ise El-Genk bağıntısıyla hesaplanmıştır. Çekirdek, yakıt içinde bir boyutlu radyal ısı iletimi ve akışkana eşlenik ısı transferi içeren toplulaştırılmış yakıt çubuğu modeliyle temsil edilmiştir.
Ana pompaların nominal koşulları şöyledir:
| Pompa parametresi | Değer |
|---|---|
| Nominal hacimsel debi | 4,65 m3·s−1 |
| Nominal açısal hız | 124,2 rad·s−1 |
| Rotor atalet momenti | 590 kg·m2 |
21. sayfadaki Şekil 8, iki kodda ortak kullanılan homolog pompa yüksekliği ve tork eğrilerini göstermektedir.
Kararlı durumda FFN ile TRACE ne kadar uyuştu?
Her iki pompanın nominal ve yarı hızlarda çalıştığı dört durum karşılaştırılmıştır:
| Sağ pompa hızı | Sol pompa hızı | TRACE çekirdek debisi | FFN çekirdek debisi | Fark |
|---|---|---|---|---|
| %100 | %100 | 4929,2 kg·s−1 | 4817,5 kg·s−1 | −%2,3 |
| %100 | %50 | 3538,5 kg·s−1 | 3459,2 kg·s−1 | −%2,2 |
| %50 | %100 | 3538,5 kg·s−1 | 3459,2 kg·s−1 | −%2,2 |
| %50 | %50 | 2454,5 kg·s−1 | 2399,1 kg·s−1 | −%2,3 |
FFN bütün çalışma noktalarında yaklaşık %2,2–2,3 daha düşük toplam debi öngörmüştür. Araştırmacılar bu düzenli farkı bağlantı yöntemindeki bir kütle korunumu hatasından çok, iki kodun sürtünme ve yerel hidrolik kayıp modellerindeki farklılıklara bağlamaktadır.
Tek pompanın yarı hızda olduğu asimetrik durumlarda çevrim debisi oranları da yakın bulunmuştur. Sağ pompanın %100, sol pompanın %50 hızda olduğu durumda TRACE sağ/sol debi oranını 3,714, FFN ise 3,870 olarak hesaplamıştır.
Sıfır güçlü pompa durması geçişi
İlk geçici durumda çekirdek gücü kullanılmamış ve yalnızca hidrolik davranış incelenmiştir. Sağ pompa durdurulduğunda rotor ataleti nedeniyle hızı kademeli olarak azalmış, sağ çevrim debisi düşmüş ve çalışan sol çevrimin ortak çekirdek üzerinden oluşturduğu basınç farkı sonucunda sağ çevrimde akış tersine dönmüştür.
23. sayfadaki Şekil 9’da TRACE’ın sürekli çizgileriyle FFN’nin kesikli çizgileri, pompa yavaşlaması, sıfır debiden geçiş ve ters akış boyunca büyük ölçüde üst üste gelmektedir. Çalışma geçiş boyunca sayısal bir hata normu veya en yüksek fark değeri vermemiştir; “çok iyi uyum” değerlendirmesi esas olarak eğrilerin görsel karşılaştırmasına dayanmaktadır.
Nominal güçlü pompa durması geçişi
İkinci geçici durumda çekirdek ısı üretimi de modele dâhil edilmiştir. Sağ pompanın durmasıyla soğutucu debisi azalmış, çekirdek çıkış sıcaklığı yükselmiş ve oluşan yoğunluk farkları kaldırma kuvvetini artırmıştır. Böylece hidrolik dirence ve pompa ataletine ek olarak doğal dolaşım etkisi de devre davranışına katılmıştır.
24. sayfadaki Şekil 10 dört değişkeni karşılaştırmaktadır:
- Sağ çevrim kütlesel debisi,
- Sağ pompanın dönme hızı,
- Çekirdek çıkış sıcaklığı,
- Sağ buhar üretecinin ısı transfer katsayısı.
FFN ve TRACE eğrileri; pompa yavaşlaması, debi yönünün değişmesi, çekirdek sıcaklığının yükselmesi ve karışık taşınım rejimine geçiş boyunca birbirine yakın seyretmiştir. Ancak çalışma bu geçiş için integral hata, zaman gecikmesi, tepe sıcaklık farkı veya belirsizlik aralığı raporlamamıştır.
Çalışmanın desteklediği temel sonuçlar
- OpenFOAM içinde farklı kesit ve yönlerdeki bir boyutlu borular, kütlesel debi korunarak bağlanabilmiştir.
- Üç veya daha fazla kol içeren bağlantılarda akış bölünmesi ve birleşmesi analitik çözümlerle %1’in altında hata vermiştir.
- Kütlesel debi ağırlıklı entalpi karışımı, çok kollu ağdaki çıkış sıcaklığını %0,03 hata ile yeniden üretmiştir.
- Basınç sıçraması modeli; dinamik basınç, yerel kayıp, pompa ve hidrostatik etkileri aynı sınır koşulunda temsil edebilmiştir.
- Dört bölgeli pompa modeli, pompa durması sonrasında düşük debi, sıfır debi ve ters akış bölgelerine geçişi modellemiştir.
- Python API, boru ağlarının bağlantı geometrilerini ve sınır koşullarını otomatik oluşturmuştur.
- Basitleştirilmiş reaktör devresinde FFN ile TRACE kararlı durum çekirdek debileri %3 içinde uyuşmuştur.
- Sıfır ve nominal güçlü pompa durması geçişlerinde FFN, TRACE’ın akış ve sıcaklık eğilimlerini yakından izlemiştir.
Çalışma neyi göstermemektedir?
- Çerçeve iki fazlı akış, kaynama, yoğuşma, boşluk oranı veya kritik ısı akısını modellememektedir.
- Bir boru kırılması, soğutucu kaybı kazası veya buhar oluşumlu reaktör geçişi doğrulanmamıştır.
- Sonuçlar deneysel bir termal-hidrolik tesisle karşılaştırılmamıştır.
- TRACE karşılaştırması bağımsız fiziksel validasyon değildir; iki hesaplama kodunun sonuçlarının tutarlılığını göstermektedir.
- Yalnızca bir basitleştirilmiş reaktör devresi ve iki pompa durması senaryosu incelenmiştir.
- Geçici durumlar için nicel hata normları, belirsizlik aralıkları veya duyarlılık analizi verilmemiştir.
- Bir boyutlu ve yüksek doğruluklu CFD bölgelerinin aynı problemde çift yönlü eşlenmesi bu çalışmada gösterilmemiştir.
- Yazılımın düzenleyici nükleer güvenlik analizlerinde kullanılmaya hazır olduğu gösterilmemiştir.
- Duvar süresi, bellek kullanımı, paralel ölçeklenebilirlik ve TRACE’a göre hesaplama maliyeti sayısal olarak raporlanmamıştır.
- Kaynak kodunun bu çalışma için kullanılan sürümü, commit kimliği ve doğrulama giriş dosyaları paylaşılmamıştır.
Çalışmanın güçlü yönleri
Çalışmanın başlıca güçlü yönü, yalnızca yeni bir çözücü tanıtmakla kalmayıp her temel bileşeni ayrı analitik problemlerle sınamasıdır. Akış bölünmesi, doğal dolaşım, pompa basıncı, yerel kayıplar ve sıcaklık karışımı birbirinden ayrılmış doğrulama problemleriyle değerlendirilmiştir.
Bağlantıların kütle, momentum ve enerji için ayrı fiziksel kurallarla ele alınması yöntemsel açıdan önemlidir. Özellikle entalpi karışımının akış yönüne göre dinamik hesaplanması, geçici olaylardaki akış tersine dönüşleriyle uyumludur.
Polinom pompa modeli yanında dört bölgeli homolog modelin bulunması, çerçeveyi yalnızca normal işletme çalışma noktasına değil, pompa durması ve ters akış koşullarına da taşımaktadır.
Python API, karmaşık devrelerin insan tarafından elle oluşturulması sırasında ortaya çıkabilecek yön, bağlantı ve dönüşüm hatalarını azaltabilecek pratik bir yazılım katmanı sunmaktadır.
Başlıca sınırlılıklar ve açıklanması gereken noktalar
En önemli kapsam sınırlaması tek fazlı formülasyondur. Nükleer santral güvenlik analizlerinde kaynama, buharlaşma, yoğuşma, iki fazlı basınç kaybı ve fazlar arası aktarım birçok olayda belirleyicidir. Yazarlar iki fazlı genişletmenin geliştirilmekte olduğunu belirtmekte, ancak bu sürümde iki fazlı denklem veya sonuç sunmamaktadır.
Analitik doğrulamalar düşük hata vermiştir; ancak ayrıntılı ağ ve zaman adımı yakınsama tabloları yalnızca sınırlı biçimde açıklanmıştır. Birinci problemde 1090 hücrenin uzaysal yakınsama için yeterli olduğu belirtilmiş, fakat farklı ağ çözünürlüklerindeki sonuçlar gösterilmemiştir.
TRACE karşılaştırmasında kararlı durum için sayısal farklar sunulurken geçici durumlarda değerlendirme büyük ölçüde grafiklerin üst üste gelmesine dayanmaktadır. Tepe değer farkı, kök ortalama kare hata veya olay zamanlaması hatası gibi ölçüler verilmemiştir.
Tablo 3’te bildirilen kapalı çevrim debilerinin birimi veya ölçeği, aynı bölümdeki geometri ve pompa eğrisiyle birlikte yeniden kontrol edilmelidir. Bu nokta analitik ve FFN değerlerinin birbirine yakın olmasını ortadan kaldırmaz; ancak her iki değerin fiziksel birimi doğru verilmediyse doğrulama sonucunun yeniden yorumlanması gerekir.
Çalışma OpenFOAM tabanlı bir yazılım geliştirmesi olmasına rağmen kullanılan yazılım sürümü, kaynak kodu etiketi, giriş dosyaları ve otomatik test paketi için açık bir erişim bilgisi vermemektedir. Bu eksiklik bağımsız yeniden üretimi sınırlandırmaktadır.
Çalışmanın Yöntemi ve Bulguları
Modelin kapsamı
| Özellik | Çalışmadaki uygulama |
|---|---|
| Temel platform | OpenFOAM tabanlı foamForNuclear |
| Boyut | Bir boyutlu boru ve sistem ağı |
| Akış rejimi | Tek fazlı |
| Korunan büyüklükler | Kütle, momentum ve enerji/entalpi |
| Geometrik bileşenler | Borular, dirsekler, bağlantılar, dallanmalar ve çap değişimleri |
| Hidrolik bileşenler | Yerel kayıplar, pompalar ve zamanla değişen basınç sıçramaları |
| Termal bileşenler | Isı kaynakları, ısı alıcıları ve bağlantılarda ideal karışım |
| Devre oluşturma | foamForNuclear Python API |
Geliştirilen yazılım bileşenleri
| Bileşen | Görevi | Temel sınırlaması |
|---|---|---|
branchConnector1D | Alan bazında çevrimsel veya sıfır gradyan bağlantısı oluşturmak | Çok kollu bağlantıda diğer kolların katkısı yarı örtük hesaplanmaktadır |
| Kütle koruyan hız bağlantısı | Farklı kesitlerde \(\alpha\rho AU\) değerini korumak | Tek fazlı mevcut formülasyona uygulanmıştır |
pressureCoupling1D | Dinamik basınç, kayıp, pompa ve hidrostatik sıçramayı uygulamak | Kayıp katsayılarının doğruluğu kullanıcı girdilerine bağlıdır |
| Polinom pompa modeli | Normal çalışma noktasında \(h(Q)\) eğrisi oluşturmak | Sıfır ve ters akışı temsil edemez |
| Homolog pompa modeli | İleri/ters dönme ve ileri/ters akışı temsil etmek | Doğru homolog baş ve tork tabloları gerektirir |
| Ataletli pompa durması | Motor torku kesildikten sonra rotor hızını hesaplamak | Atalet ve sürtünme modeli kullanıcı tarafından tanımlanmalıdır |
| Entalpi karışımı | Bağlantıya giren akımları kütlesel debiyle ağırlıklandırmak | Bağlantı içinde ideal ve anlık karışım varsaymaktadır |
| Python boru API’si | Ağ, yön ve sınır koşullarını otomatik oluşturmak | Çalışmada tam giriş dosyaları ve sürüm kimliği paylaşılmamıştır |
Analitik doğrulama matrisi
| Doğrulama problemi | Sınanan özellik | Ağ veya koşul | En yüksek bildirilen hata |
|---|---|---|---|
| Akış bölünmesi ve birleşmesi | Bağlantı kütle korunumu ve sürtünme kaybı | 1090 hücre, 570 L·s−1 giriş | %0,33 |
| Pompalı kapalı çevrim | Pompa yüksekliği ve dağıtılmış kayıplar | 320 hücre, dört farklı çaplı boru | <%0,05 |
| Doğal dolaşımlı çevrim | Sıcaklığa bağlı yoğunluk ve kaldırma kuvveti | 300–350 K sıcaklık farkı | <%0,05 |
| Çoklu dallanma ve ısıtma | Akış dağılımı ve entalpi karışımı | 30 MW ısı kaynağı | %0,16 |
TRACE karşılaştırma modelinin ana özellikleri
| Alt sistem | Modelleme yaklaşımı |
|---|---|
| Reaktör devresi | İki çevrimli basitleştirilmiş PWR benzeri birincil sistem |
| Çekirdek | Toplulaştırılmış yakıt çubuğu ve bir boyutlu radyal iletim |
| Isı transferi | Eşlenik yakıt-soğutucu aktarımı ve El-Genk bağıntısı |
| Buhar üreteçleri | Sabit ikincil taraf sıcaklığına sahip ısı alıcıları |
| Sürtünme | Churchill bağıntısı ve bileşene özgü pürüzlülük |
| Pompalar | Dört bölgeli homolog baş ve tork eğrileri |
| Geçici olay | Sağ ana soğutucu pompasının durması |
| Karşılaştırılan senaryolar | Sıfır güç ve nominal güç |
Ana nicel bulgular
| Bulgular | Sonuç | Yorum sınırı |
|---|---|---|
| Analitik hidrolik doğrulama | Bildirilen hatalar %0,33 veya altında | Seçilmiş idealize problemler |
| Enerji karışımı | Çıkış sıcaklığı hatası %0,03 | Sabit özellik ve ideal karışım varsayımı |
| TRACE nominal kararlı durum | FFN çekirdek debisi %2,3 daha düşük | İki kodun bağıntıları tamamen aynı değildir |
| Asimetrik pompa koşulları | Toplam debi farkı %3’ün altında | Yalnızca dört hız kombinasyonu |
| Sıfır güçlü pompa durması | Akış yavaşlaması ve tersine dönüş görsel olarak yakın | Sayısal geçici durum hata metriği verilmemiştir |
| Nominal güçlü geçiş | Debi, pompa hızı, sıcaklık ve ısı transfer eğrileri yakın | Tek geçici olay ve tek reaktör düzeni |
Tekrarlanabilirlik ve sayısal güvenilirlik
Çalışma temel denklemleri, sınır koşullarını, analitik karşılaştırma değerlerini ve reaktör devresinin ana model varsayımlarını ayrıntılı biçimde vermektedir. Bu durum yöntemsel değerlendirmeyi kolaylaştırmaktadır.
Bununla birlikte bağımsız tekrar için aşağıdaki bilgiler eksik veya sınırlıdır:
- Kullanılan OpenFOAM ve foamForNuclear sürümleri,
- Kaynak kodu commit veya yayın etiketi,
- TRACE giriş modeli ve sürümü,
- Bütün boruların çap, uzunluk, pürüzlülük ve kayıp katsayıları,
- Zaman adımı ve doğrusal çözücü toleranslarının tam listesi,
- Geçici durum ağ ve zaman adımı duyarlılığı,
- Duvar süresi, CPU kullanımı ve paralel ölçeklenebilirlik,
- TRACE karşılaştırması için nicel geçici durum hata ölçüleri,
- Analitik testlerin otomatik çalıştırılabileceği giriş dosyaları.
Deterministik hesaplama çalışması olduğu için deneysel tekrar veya klasik istatistiksel anlamlılık testi beklenmemektedir. Buna karşılık sayısal belirsizlik; ağ çözünürlüğü, zaman adımı, iterasyon toleransı, kapanış bağıntıları ve giriş parametrelerinin duyarlılığıyla değerlendirilmelidir. Bu analizler mevcut çalışmada sınırlı düzeyde sunulmuştur.
Kaynak ve Yöntem Notu
Çalışmanın tam özgün adı: Development of a One-Dimensional Thermal-Hydraulics Framework for Nuclear Reactor Analysis in OpenFOAM
Yazarlar ve sıraları: Giovanni Nervi, Alessandro Scolaro, Thomas Guilbaud, Matteo Bettiol, Mathieu Hursin.
Eş katkı veya eş birinci yazarlık: Çalışmada eş katkı veya eş birinci yazarlık açıklaması bulunmamaktadır.
Sorumlu veya iletişim yazarı: Yüklenen çalışmada sorumlu yazarı belirten yıldız veya e-posta adresi yer almamaktadır. SSRN kayıt sayfası Giovanni Nervi’yi “Contact Author” olarak göstermektedir.
Kurumlar: Giovanni Nervi ve Mathieu Hursin için Laboratory for Reactor Physics and Systems Behaviour, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland; Alessandro Scolaro için Laboratory for Reactor Physics and Thermal-Hydraulics, Paul Scherrer Institute, Villigen, Switzerland; Thomas Guilbaud için Transmutex SA, Geneva, Switzerland; Matteo Bettiol için Laboratory for Simulation and Modelling, Paul Scherrer Institute, Villigen, Switzerland belirtilmiştir.
DOI: 10.2139/ssrn.7197741. Bu DOI, SSRN ön baskı kaydına aittir ve hakemli bir dergi makalesi DOI’si değildir.
Dergi veya konferans: İncelenen sürüm için doğrulanmış bir hakemli dergi veya konferans yayını belirtilmemiştir.
Yayın platformu: SSRN.
Özgün yayınevi: Çalışma için doğrulanmış bir nihai dergi yayınevi bulunmamaktadır. SSRN, ön baskının dağıtıldığı platformdur.
Yayın yılı: 2026.
Kaynak türü: Hesaplamalı nükleer mühendislik, tek fazlı sistem termal-hidroliği ve bilimsel yazılım geliştirme araştırma preprinti.
Hakemlik durumu: Çalışma hakem değerlendirmesinden geçmemiştir. Her sayfada “This preprint research paper has not been peer reviewed” ve “Preprint not peer reviewed” uyarıları bulunmaktadır.
Resmî bağlantılar:SSRN resmî kayıt sayfası ve SSRN DOI bağlantısı.
Finansman: Çalışma, Avrupa Birliği tarafından Grant Agreement 101059543 kapsamında finanse edilen ESFR-SIMPLE projesiyle desteklenmiştir.
Yazar katkıları: Giovanni Nervi kavramsallaştırma, veri düzenleme, biçimsel analiz, araştırma, yöntem, yazılım, doğrulama, görselleştirme ve ilk taslak; Alessandro Scolaro kavramsallaştırma, biçimsel analiz, yöntem, proje yönetimi, yazılım, gözetim ve inceleme; Thomas Guilbaud kavramsallaştırma, yazılım ve inceleme; Matteo Bettiol kavramsallaştırma, yazılım, doğrulama ve inceleme; Mathieu Hursin kavramsallaştırma, proje yönetimi ve gözetim görevleriyle belirtilmiştir.
Veri ve kod erişimi: İncelenen çalışmada doğrulama giriş dosyaları, kullanılan yazılım commit kimliği, TRACE modeli veya sonuçların ham sayısal verileri için açık bir veri deposu belirtilmemiştir. Python API’den örnek kod verilmiştir; ancak bu örnek çalışmanın tamamını bağımsız olarak yeniden üretmek için yeterli değildir.
Bu Türkçe makale, yüklenen 27 sayfalık çalışmanın yöntem açıklamaları, 27 temel denklemi, dört ana sonuç tablosu, T-bağlantı ve devre şemaları, sıcaklık ve hız dağılımları, TRACE düğümlemesi, homolog pompa eğrileri ve pompa durması geçiş grafikleri esas alınarak hazırlanmıştır. Çalışmada bulunmayan yeni bir sayısal sonuç veya dış kaynaktan bilimsel performans bulgusu eklenmemiştir. Dış kaynak kullanımı yalnızca SSRN kaydı, DOI, iletişim yazarı ve kurumsal kimlik bilgilerinin bibliyografik doğrulanmasıyla sınırlandırılmıştır.
Çalışmanın başlıca sınırlılıkları; hakem değerlendirmesinden geçmemiş olması, modelin yalnızca tek fazlı akışlarla sınırlı kalması, fiziksel deneysel validasyon bulunmaması, yalnızca bir basitleştirilmiş reaktör devresinin TRACE ile karşılaştırılması, geçici durumlar için sayısal hata ölçülerinin verilmemesi, performans ve ölçeklenebilirlik ölçümlerinin raporlanmaması, yeniden üretilebilir kod ve giriş dosyalarının paylaşılmaması ve kapalı çevrim doğrulamasındaki debi birimi veya ölçeğinin açıklamaya ihtiyaç göstermesidir.
Sonuçlar, geliştirilen çerçevenin seçilmiş tek fazlı analitik testlerde korunum denklemlerini doğru uyguladığını ve basitleştirilmiş pompa durması geçişlerinde TRACE ile benzer eğilimler üretebildiğini göstermektedir. Bulgular, gerçek bir nükleer tesis için lisanslama yeterliliği, iki fazlı kaza analizi, deneysel validasyon veya düzenleyici kullanım onayı anlamına gelmemektedir.
Bir yorum bırakın
E-posta adresiniz yayınlanmayacaktır. Gerekli alanlar * ile işaretlenmiştir