Бул изилдөө, ядролук реакторlerin түтүктөр, насостор, байланыштар, жылуулук булактары ve гидравликалык жоготууларdan oluşan bütün контурlerini OpenFOAM ortamında бир өлчөмдүү olarak modellemek için yeni bir термалдык-гидравликалык çerçeve иштеп чыгат. Yöntem, foamForNuclear platformunun mevcut бир фазалуу масса, импульс ve энергия denklemlerini değiştirmek yerine; бир өлчөмдүү bileşenleri масса, импульс ve энергияyi сактап birbirine bağlayan yeni чек шарттары, көп бутактуу шилтеме modelleri, басым секириги hesapları, төрт аймактуу насос модели ve автоматтык тор түзүү araçları eklemektedir. Analitik doğrulama problemlerinde агым, басым ve температура жыйынтыктарı %1’in altında салыштырма катаlar vermiş; basitleştirilmiş iki çevrimli bir реактор контуруnde TRACE ile yapılan karşılaştırmada туруктуу абал өзөк чыгымдары yaklaşık %3 içinde uyuşmuştur. Bununla birlikte изилдөө yalnızca бир фазалуу агымları kapsamaktadır ve тажрыйбаsel bir реактор sistemiyle değil, аналитикалык чечимler ve başka bir sistem koduyla karşılaştırılmıştır.
Yeni çerçevenin temel amacı, düşük hesaplama maliyetli бир өлчөмдүү sistem modelleriyle yüksek çözünürlüklü CFD veya көзөнөктүү чөйрө modellerini aynı OpenFOAM altyapжылуулукnda birleştirebilecek ortak bir temel oluşturmaktır. Geliştirilen branchConnector1D чек шарты, bağlı түтүктөрın kesitleri farklı olduğunda ылдамдыкları eşitlemek yerine массалык чыгымyi korumaktadır. pressureCoupling1D ise dinamik басым, жергиликтүү гидравликалык жоготуу, насос basıncı ve hidrostatik farkları tek bir басым секириги denkleminde birleştirmektedir. Dört bölgeli гомологдук насос modeli, насос токтошу sırasında чыгымnin sıfıra yaklaşmasını, tersine dönmesini ve насосnın инерцияle yavaşlamasını temsil eчыгымlmektedir.
Түркия жагынан баалоо: Изилдөө; Türkiye’de nükleer энергия mühendisliği, реактор güvenlik analizi, araştırma реакторleri, küçük modüler реактор изилдөөları, энергия sistemleri eğitimi ve açık булакlı mühendislik yazılımı geliştirme faaliyetleri açжылуулукndan ыкмаsel маани taşımaktadır. Yaklaşımın Türkiye’deki araştırma ve tasarım projelerinde kullanılabilmesi için hedef реакторün gerçek geometri, материал, насос, vana, жылуулук değiştirici ve kontrol verileriyle yeniden kurulması; тажрыйбаsel контурler 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.
Эмне үчүн бир өлчөмдүү реактор системалык коддору керек?
Nükleer santrallerin термалдык-гидравликалык davranışı; реактор çekirdeği, sıcak ve soğuk kollar, насостор, буу üreteçleri, басым kapları, шилтеме түтүктөрı ve güvenlik sistemlerinin birlikte изилдөөsı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ı, түтүк ve ekipmanları eksen boyunca sırаянтmış kontrol көлөмleri olarak temsil eder. Böylece басым, чыгым, температура ve энергия dağılımı bütün реактор контуруnde görece düşük hesaplama maliyetiyle izlenebilir. TRACE, RELAP5 ve ATHLET gibi araçlar bu yaklaşımın yerleşik örnekleridir.
Ancak yeni реактор tasarımlarında табигый айлануу, büyük havuzlar, termal tabakalaşma ve karmaşık аралашма bölgeleri gibi üç boyutlu olaylar da önem kazanmıştır. Bu nedenle tek bir çözünürlük düzeyi yerine, bazı bölgelerin бир өлчөмдүү sistem modeliyle, bazı bölgelerin көзөнөктүү чөйрө yaklaşımıyla, критикалык bölgelerin ise CFD ile çözüldüğü çok doğruluklu ыкмаler geliştirilmektedir.
foamForNuclear ичиндеги кемчилик эмне эле?
foamForNuclear, OpenFOAM tabanlı çok fizikli ve çok doğruluklu bir nükleer analiz platformudur. Platform, GeN-Foam ile OFFBEAT чечүүчүlerinin geliştirilmesi ve birleştirilmesiyle; агым, жылуулук transferi, nötronik ve yагымt davranışı gibi farklı fiziklerin ortak bir ortamda çalıştırılmasını amaçlamaktadır.
Изилдөөya göre foamForNuclear, CFD ve көзөнөктүү чөйрө tabanlı термалдык-гидравликалык bölgeleri aynı simülasyon içinde çalıştırabilmesine rağmen bütün santral контурsini kurmak için gereken жергиликтүү bir бир өлчөмдүү sistem altyapжылуулукna sahip değildi. Pompa, vana, dirsek, T-шилтеме, çap değişimi ve жергиликтүү жоготуу gibi bileşenler için haricî bir sistem koduna veya OpenModelica шилтемеsına ihtiyaç duyuluyordu.
Изилдөөчүлөр bu bağımlılığı kaldırmak için OpenFOAM içinde çalışan жергиликтүү bir çerçeve geliştirmiştir. Изилдөөnın “OpenFOAM içine doğrudan gömülmüş ilk uygulama” olduğu yönündeki ifade изилдөөчүлөрın kendi öncelik iddiasıdır; изилдөө içinde kapsamlı ve bağımsız bir yazılım önceliği incelemesi sunulmamıştır.
Иштелип чыккан алкак эмнени өзгөртөт?
Изилдөө, foamForNuclear’ın mevcut бир фазалуу көзөнөктүү чөйрө термалдык-гидравликалык denklemlerini yeniden yazmamaktadır. Kütle, импульс ve энтальпия сакталууunu çözen аянт denklemleri aynı kalmaktadır. Yeni изилдөө esas olarak şu altyapıyı eklemektedir:
- Bir boyutlu түтүктөрı farklı yön ve kesitlerde birbirine bağlayan чек шарттары,
- Kesit değişimlerinde массалык чыгымyi koruyan ылдамдык eşleştirmesi,
- Üç veya daha fazla kolun birleştiği байланыштар,
- Dinamik басым, жергиликтүү жоготуу, насос ve hidrostatik etkileri içeren басым секириги,
- İleri ve geri агымı kapsayan насос karakteristikleri,
- Pompa durması ve инерцияle yavaşlama modeli,
- Bağlantılarda массалык чыгым ağırlıklı энтальпия аралашмаı,
- Python üzerinden otomatik бир өлчөмдүү түтүк ağı oluşturma araçları.
OpenFOAMдын учурдагы бир өлчөмдүү ыкмасы эмне үчүн жетиштүү эмес эле?
OpenFOAM, ince bir üç boyutlu ağın bazı yüzeylerinde empty чек шартыnu kullanarak iki veya бир өлчөмдүү çö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 түтүк hatlarında sorun oluşturmaktadır. Örneğin 90° dirsek içeren bir контур, tek doğrultulu düz bir ağ gibi temsil edilemez.
Bir çözüm olarak түтүктөрın dış yüzeylerinde kaymalı duvar, байланыштарda ise cyclic veya uyumsuz ağlar için cyclicAMI kullanılabilir. Bu ыкма түтүктөрı farklı açılarla bağlayabilse de изилдөө üç temel sorun belirlemiştir:
- Kütle сакталууu: Kesit аянтları farklı iki түтүкnun sınırında aynı ылдамдыкın uygulanması, aynı массалык чыгымnin geçtiğini garanti etmez.
- Bütün аянтlara aynı sınır tipinin uygulanması: Basınç, ылдамдык ve энтальпия bağlı olmalıdır; ancak hidrolik çap gibi yardımcı bir аянт sınır boyunca үзгүлтүксүз olmak zorunda değildir.
- Karmaşık ağ oluşturma: Her шилтеме için öteleme vektörü, dönme açжылуулук ve dönme merkezinin elle tanımlanması büyük контурlerde hata riskini artırmaktadır.
branchConnector1D эмне кылат?
Изилдөөчүлөр farklı аянт değişkenlerine farklı davranışlar atayabilmek için branchConnector1D adlı bir шилтеме чек шарты geliştirmiştir. Basınç, ылдамдык, энтальпия ve seçilen türbülans değişkenleri gibi fiziksel olarak bağlanması gereken аянтlar çevrimsel шилтемеyla eşleştirilirken, hidrolik çap gibi шилтеме boyunca farklı kalabilen аянтlarda sıfır gradyan uygulanmaktadır.
Bu ayrım, yalnızca korunması veya үзгүлтүксүз olması gereken değişkenlerin karşı түтүкyla güçlü biçimde bağlanmasını sağlamaktadır. Böylece bütün аянтların zorunlu olarak aynı мааниe getirildiği yapay bir шилтеме önlenmektedir.
Ар түрдүү диаметрдеги түтүктөрдө масса кантип сакталат?
İki түтүкnun kesitleri, тыгыздыкları veya etkin агым аянтları farklı olduğunda ылдамдыкların eşitlenmesi yerine массалык чыгымlerin eşitlenmesi gerekir. Изилдөөnı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ığı катышы, \(\rho\) yoğunluğu, \(A\) kesit аянтını ve \(U\) ылдамдыкı көрсөтөт. \(P\) ve \(N\), шилтемеnın iki tarafındaki kontrol көлөмlerini ifade etmektedir.
Komşu taraftan gelen ылдамдык katkжылуулук kesit, тыгыздык ve faz катышы kullanılarak ölçeklendirilmektedir. Böylece küçük kesitli түтүкda ылдамдык yükselirken büyük kesitli түтүкda düşebilmekte, fakat шилтемеdan geçen toplam массалык чыгым aynı kalmaktadır.
Көп бутактуу байланыштар кантип моделденген?
Standart çevrimsel OpenFOAM шилтемеsı iki yüzey arasındaki ikili ilişki için tasarlanmıştır. T-шилтемеda ise bir ana түтүк iki kola ayrılabilir veya iki kol tek түтүкda birleşebilir. Изилдөөnın 7. sayfasındaki Şekil 1, ana түтүкnun ucunda iki ayrı шилтеме yüzeyi oluşturularak üst ve alt kollara bağlandığı düzeni көрсөтөт.
Üç kollu шилтемеdaki масса сакталууu şu сүрөтde 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 шилтемеya doğru olduğunda pozitif, шилтемеdan uzaklaştığında negatif kabul edilmektedir. Bir sınır yüzeyinin etkili komşu ылдамдыкı, yalnızca doğrudan eşleştirildiği yüzeyden değil, шилтемеdaki diğer bütün kollardan gelen массалык чыгымlerden de hesaplanmaktadır.
Diğer kollardan gelen katkı açık bir булак terimi olarak eklendiği için formülasyon tamamen değil, yarı örtüktür. Dengesiz чыгым dağılımlarında yагымnsamayı iyileştirmek amacıyla bu açık terime isteğe bağlı gevşetme uygulanabilmektedir.
Bir hücrenin ucunda birden fazla шилтеме yüzeyi oluşturulması geometrik hücre hacmini bozabileceğinden, изилдөөчүлөр шилтемеya komşu hücrelerin көлөмlerini ve басым-ылдамдык bağlaşımındaki yüzey katkılarını yeniden hesaplamıştır.
Басым секириги кайсы таасирлерди камтыйт?
Bağlantıdan geçen массалык чыгымnin korunması, шилтемеnın iki tarafındaki басымların mutlaka eşit olması gerektiği anlamına gelmez. Kesit değişimi, dirsek, vana, насос ve yükseklik farkı басымta fiziksel bir sıçrama meydana getirebilir.
Bu amaçla pressureCoupling1D adlı bir чек шарты geliştirilmiştir. Kaldırma kuvveti içeren OpenFOAM чечүүчүlerinde taşınan değiştirilmiş басым:
\[ p_{rgh}=p-\rho gz \]
olarak tanımlanmaktadır. İki taraf arasındaki басым секириги изилдөөda ş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 басым farkını, ikinci terim жергиликтүү hidrolik kaybı, \(S_p\) насосnın oluşturduğu басым artışını ve son terim тыгыздык ile yükseklikten булакlanan hidrostatik düzeltmeyi temsil etmektedir.
Yerel kayıp katsayжылуулук \(K\) ve насос kaynağı zamanla değişen fonksiyonlar olarak verilebildiğinden vana кыймылi veya насос токтошу gibi өтмө olaylar modellenebilmektedir. Akış yön değiştirdiğinde kayıp terimi yeni giriş tarafına göre yeniden yönlendirilmektedir.
Насос үчүн кайсы эки модель иштелип чыккан?
İlk model, насос yüksekliğini көлөмsel чыгымnin polinomu olarak tanımlamaktadır:
\[ h(Q)=\sum_i c_iQ^i \]
Pompa ылдамдыкı değiştiğinde katsayılar benzerlik yasalarına göre ölçeklendirilmektedir. Bu model tasarım изилдөө noktasına yагымn ileri агым koşullarında sade ve kullanışlıdır. Ancak \(h(Q)\) tek мааниli olduğundan sıfır чыгым, ters чыгым veya ters dönme gibi durumları eksiksiz temsil edemez.
İkinci model, sistem kodlarında kullanılan homolog veya төрт аймактуу насос yaklaşımıdır. Debi, dönme ылдамдыкı ve кысуу yüksekliği nominal мааниlere 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 чыгым ile ileri ve ters dönme koşulları dört изилдөө bölgesine ayrılmaktadır. Her bölgede насос yüksekliği ve hidrolik tork, parçalı doğrusal таблицаlarla hesaplanmaktadır. Bu düzen, geleneksel sekiz oktantlı насос karakteristiklerini benzer normalleştirmeye sahip oktantları birleştirerek dört bölgede temsil etmektedir.
Насостун токтошу жана инерция менен жайлашы кантип эсептелет?
Pompa ылдамдыкı kullanıcı tarafından zaman çizelgesiyle verilebilir veya насос durduktan sonra ротор dinamiğiyle hesaplanabilir. Motor torku kesildiğinde dönme ылдамдыкı şu dengeyle değişmektedir:
\[ I\frac{d\omega}{dt} = -T_{\mathrm{hy}}\left(q^*,\omega^*\right) -T_f(\omega) \]
Burada \(I\) ротор инерция momentini, \(T_{\mathrm{hy}}\) суюктукın oluşturduğu hidrolik torku ve \(T_f\) механикалык сүрүлүү torkunu ifade etmektedir.
Dört bölgeli насос karakteristiğiyle ротор denklemi birlikte çözüldüğünde насос normal изилдөө noktasından düşük чыгымye, sıfır чыгымye ve тескери агым bölgesine kesintisiz geçebilmektedir. Bu касиет, насос kaybı ve агым yönünün değiştiği реактор geçişleri için önemlidir.
Байланыштарда энергиянын сакталуусу кантип камсыздалган?
Birden fazla kolun birleştiği noktada çıkış энтальпияsi tek bir komşu түтүкnun энтальпияsine eşitlenemez. Çıkış суюктукı, шилтемеya giren bütün агымmların аралашмаıdır. Изилдөөda giriş энтальпияsi массалык чыгым ağırlıklı орточоyla 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 шилтемеya doğru akan kollarda yapılmaktadır. Bir kol шилтемеdan dışarı агымyorsa энтальпия sıfır gradyan davranışıyla tжешилүүktadır. Akış yönü өтмө olay sırasında değiştiğinde аралашмаa katılan kollar da dinamik olarak yeniden belirlenmektedir.
Python интерфейси контур курууну кантип жеңилдетет?
foamForNuclear Python API’si, бир өлчөмдүү түтүк parçalarının başlangıç konumu, yönü, uzunluğu, hidrolik çapı, hücre sayжылуулук ve dirsek yarıçapıyla tanımlanmasına olanak verecek biçimde genişletilmiştir.
Yeni bir түтүкnun başlangıcı bir koordinat vektörü yerine mevcut bir түтүк nesnesi olarak verilebilmektedir. Bu durumda yeni түтүк önceki түтүкnun çıkışına otomatik olarak bağlanmakta; gerekli dönme, öteleme ve чек шарттары yazılım tarafından oluşturulmaktadır.
Bağımsız түтүктөр için giriş ve çıkış yüzeyleri, bağlı түтүктөр için ise branchConnector1D байланыштарı otomatik üretilmektedir. Böylece büyük bir контурnin her шилтемеsı için dönüş merkezi ve geometrik dönüşümün elle hesaplanması gerekmemektedir.
Биринчи аналитикалык текшерүү: агымдын бөлүнүшү жана кайра биригиши
İlk doğrulama probleminde giriş чыгымsi 570 L·s−1 olan su агымı, çapları, uzunlukları ve сүрүлүү katsayıları farklı iki paralel kola ayrılmış ve daha sonra yeniden birleşmiştir. İki T-шилтеме arasında 9 m yükseklik farkı bulunmaktadır.
1090 hücreli ağda elde edilen жыйынтыктар şöyledir:
| Değişken | Analitik çözüm | FFN sonucu | Bağıl hata |
|---|---|---|---|
| Birinci kol чыгымsi | 121,8 L·s−1 | 122,2 L·s−1 | %0,33 |
| İkinci kol чыгымsi | 448,2 L·s−1 | 447,8 L·s−1 | %0,089 |
| A–B басым farkı | 147 kPa | 147,15 kPa | %0,10 |
Bu жыйынтыктар, farklı kollara dağılan массалык чыгымnin ve dağıtılmış сүрүлүү kayıplarının аналитикалык чечимle ылайык biçimde hesaplandığını көрсөтөт.
Экинчи текшерүү: насостуу жана табигый айланма жабык контур
İkinci problem, çapları 80, 90, 100 ve 110 mm olan dört түтүкdan meydana gelen kapalı bir çevrimdir. Her түтүк 40 m uzunluğundadır. İlk düzende агым bir насосyla, ikinci düzende ise alt түтүкnun 350 K’ye жылуулукtılması ve üst түтүкnun 300 K’ye soğutulmasıyla oluşan kaldırma kuvvetiyle sağlanmıştır.
| Durum | Analitik чыгым | FFN чыгымsi | Bildirilen hata |
|---|---|---|---|
| Pompalı çevrim | 0,0976 m3·s−1 | 0,0976 m3·s−1 | <%0,05 |
| Doğal айланууlı çevrim | 0,7345 m3·s−1 | 0,7341 m3·s−1 | <%0,05 |
Kaynak, analitik ve sayısal жыйынтыктарın birbirine çok yагымn olduğunu көрсөтөт. Bununla birlikte таблицаda m3·s−1 olarak verilen чыгымlerin büyüklüğü; belirtilen түтүк çapları, uzunlukları ve насос eğrisiyle birlikte мааниlendirildiğinde birim veya ölçek açıklaması gerektirmektedir. Örneğin булакta verilen \(H_{\mathrm{pump}}=0{,}15-15Q^2\) eşitliğinde \(Q=0{,}0976\) kullanıldığında насос yüksekliği yalnızca yaklaşık 0,0071 m olmaktadır. Bu nedenle Tablo 3’teki чыгым biriminin veya ondalık gösteriminin yazarlar tarafından kontrol edilmesi gerekir.
Үчүнчү текшерүү: көп бутактар жана энергия аралашуусу
Son analitik problem, art arda iki агым bölünmesi, iki рекомбинация ve 30 MW’lık жергиликтүү жылуулук kaynağı içeren daha karmaşık bir түтүк ağıdır. Giriş ылдамдыкı 1 m·s−1, dirsek kayıp katsayжылуулук 0,9 ve T-шилтеме kayıp katsayжылуулук 1,2 olarak колдонулган.
Изилдөөnın 18. sayfasındaki Şekil 5, жылуулукtılan üst kolun sıcaklığının yükseldiğini ve аралашма байланыштарından sonra sıcaklığın diğer kollara dağıldığını көрсөтөт. 19. sayfadaki Şekil 6 ise paralel kollardaki farklı ылдамдык büyüklüklerini görselleştirmektedir.
| Değişken | Analitik çözüm | FFN sonucu | Bağıl hata |
|---|---|---|---|
| Alt kol чыгымsi \(Q_b\) | 0,4367 m3·s−1 | 0,4369 m3·s−1 | %0,06 |
| Isıtılan kol чыгымsi \(Q_h\) | 0,1442 m3·s−1 | 0,1444 m3·s−1 | %0,16 |
| Isıtılmayan üst alt-kol чыгымsi \(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 чыгымnin değil, paralel kollardaki dağılımın ve farklı температураlardaki агымmların birleşiminden oluşan çıkış энтальпияsinin de analitik жыйынтыкla ылайык olduğunu көрсөтөт.
TRACE салыштыруусу үчүн кандай реактор контуру түзүлгөн?
Sistem düzeyindeki баалооde, басымlı su реакторüne benzeyen basitleştirilmiş iki çevrimli birincil контур kullanılmıştır. Modelde:
- Reaktör çekirdeği,
- Alt ve üst plenumlar,
- İki simetrik sıcak kol,
- İki буу üreteci temsili,
- İki ana soğutucu насосsı,
- İki soğuk kol,
- Plenumlar arasında çapraz шилтеме
bulunmaktadır. Изилдөөnın 20. sayfasındaki Şekil 7, TRACE düğümlemesiyle OpenFOAM/FFN түтүк ağını yan yana көрсөтөт.
Dağıtılmış сүрүлүү Churchill bağıntжылуулукyla, çekirdek ve буу üreteçlerindeki taşınım ise El-Genk bağıntжылуулукyla hesaplanmıştır. Çekirdek, yагымt içinde бир өлчөмдүү radyal жылуулук iletimi ve суюктукa eşlenik жылуулук transferi içeren toplulaştırılmış yагымt çubuğu modeliyle temsil edilmiştir.
Ana насосторın nominal koşulları şöyledir:
| Pompa parametresi | Değer |
|---|---|
| Nominal көлөмsel чыгым | 4,65 m3·s−1 |
| Nominal açısal ылдамдык | 124,2 rad·s−1 |
| Rotor инерция momenti | 590 kg·m2 |
21. sayfadaki Şekil 8, iki kodda ortak kullanılan гомологдук насос yüksekliği ve tork eğrilerini көрсөтөт.
Туруктуу абалда FFN менен TRACE канчалык шайкеш келди?
Her iki насосnın nominal ve yarı ылдамдыкlarda çalıştığı dört durum karşılaştırılmıştır:
| Sağ насос ылдамдыкı | Sol насос ылдамдыкı | TRACE çekirdek чыгымsi | FFN çekirdek чыгымsi | 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 изилдөө noktalarında yaklaşık %2,2–2,3 daha düşük toplam чыгым öngörmüştür. Изилдөөчүлөр bu düzenli farkı шилтеме ыкмаindeki bir масса сакталууu hatasından çok, iki kodun сүрүлүү ve жергиликтүү гидравликалык жоготуу modellerindeki farklılıklara bağlamaktadır.
Tek насосnın yarı ылдамдыкda olduğu asimetrik durumlarda çevrim чыгымsi oranları da yагымn bulunmuştur. Sağ насосnın %100, sol насосnın %50 ылдамдыкda olduğu durumda TRACE sağ/sol чыгым катышыnı 3,714, FFN ise 3,870 olarak hesaplamıştır.
Нөл кубаттагы насос токтошу өтмөсү
İlk өтмө durumda çekirdek gücü kullanılmamış ve yalnızca hidrolik davranış изилденген. Sağ насос durdurulduğunda ротор инерцияi nedeniyle ылдамдыкı kademeli olarak azalmış, sağ çevrim чыгымsi düşmüş ve çalışan sol çevrimin ortak çekirdek üzerinden oluşturduğu басым farkı sonucunda sağ çevrimde агым tersine dönmüştür.
23. sayfadaki Şekil 9’da TRACE’ın үзгүлтүксүз çizgileriyle FFN’nin kesikli çizgileri, насос yavaşlaması, sıfır чыгымden geçiş ve тескери агым boyunca büyük ölçüde üst üste gelmektedir. Изилдөө geçiş boyunca sayısal bir hata normu veya en yüksek fark мааниi vermemiştir; “çok iyi uyum” баалооsi esas olarak eğrilerin görsel karşılaştırmasına dayanmaktadır.
Номиналдык кубаттагы насос токтошу өтмөсү
İkinci өтмө durumda çekirdek жылуулук üretimi de modele dâhil edilmiştir. Sağ насосnın durmasıyla soğutucu чыгымsi azalmış, çekirdek çıkış sıcaklığı yükselmiş ve oluşan тыгыздык farkları kaldırma kuvvetini artırmıştır. Böylece hidrolik dirence ve насос инерцияine ek olarak табигый айлануу etkisi de контур davranışına katılmıştır.
24. sayfadaki Şekil 10 dört değişkeni karşılaştırmaktadır:
- Sağ çevrim массалык чыгымsi,
- Sağ насосnın dönme ылдамдыкı,
- Çekirdek çıkış sıcaklığı,
- Sağ буу üretecinin жылуулук transfer katsayжылуулук.
FFN ve TRACE eğrileri; насос yavaşlaması, чыгым 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 yагымn seyretmiştir. Ancak изилдөө bu geçiş için integral hata, zaman gecikmesi, tepe температура farkı veya belirsizlik aralığı raporlamamıştır.
Изилдөө колдогон негизги жыйынтыктар
- OpenFOAM içinde farklı kesit ve yönlerdeki бир өлчөмдүү түтүктөр, массалык чыгым korunarak bağlanabilmiştir.
- Üç veya daha fazla kol içeren байланыштарda агым bölünmesi ve birleşmesi аналитикалык чечимlerle %1’in altında hata vermiştir.
- Kütlesel чыгым ağırlıklı энтальпия аралашмаı, көп бутактуу ağdaki çıkış sıcaklığını %0,03 hata ile yeniden üretmiştir.
- Basınç sıçraması modeli; dinamik басым, жергиликтүү жоготуу, насос ve hidrostatik etkileri aynı чек шартыnda temsil eчыгымlmiştir.
- Dört bölgeli насос модели, насос токтошу sonrasında düşük чыгым, sıfır чыгым ve тескери агым bölgelerine geçişi modellemiştir.
- Python API, түтүк ağlarının шилтеме geometrilerini ve чек шарттарыnı otomatik oluşturmuştur.
- Basitleştirilmiş реактор контуруnde FFN ile TRACE туруктуу абал өзөк чыгымдары %3 içinde uyuşmuştur.
- Sıfır ve nominal güçlü насос токтошу geçişlerinde FFN, TRACE’ın агым ve температура eğilimlerini yагымndan izlemiştir.
Изилдөө эмнени көрсөтпөйт?
- Çerçeve iki fazlı агым, kaynama, yoğuşma, boşluk катышы veya критикалык жылуулук агымsını modellememektedir.
- Bir түтүк сынууsı, soğutucu kaybı kazası veya буу oluşumlu реактор geçişi doğrulanmamıştır.
- Sonuçlar тажрыйбаsel bir термалдык-гидравликалык tesisle karşılaştırılmamıştır.
- TRACE karşılaştırması bağımsız fiziksel validasyon değildir; iki hesaplama kodunun жыйынтыктарının tutarlılığını көрсөтөт.
- Yalnızca bir basitleştirilmiş реактор контуру ve iki насос токтошу senaryosu изилденген.
- 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 бул изилдөөda gösterilmemiştir.
- Yazılımın düzenleyici nükleer güvenlik analizlerinde kullanılmaya hazır olduğu gösterilmemiştir.
- Duvar мөөнөтsi, bellek kullanımı, paralel ölçeklenebilirlik ve TRACE’a göre hesaplama maliyeti sayısal olarak raporlanmamıştır.
- Kaynak kodunun бул изилдөө için kullanılan sürümü, commit kimliği ve doğrulama giriş dosyaları paylaşılmamıştır.
Изилдөөнүн күчтүү жактары
Изилдөөnın başlıca güçlü yönü, yalnızca yeni bir чечүүчү tanıtmakla kalmayıp her temel bileşeni ayrı analitik problemlerle sınamasıdır. Akış bölünmesi, табигый айлануу, насос basıncı, жергиликтүү жоготууlar ve температура аралашмаı birbirinden ayrılmış doğrulama problemleriyle мааниlendirilmiştir.
Bağlantıların масса, импульс ve энергия için ayrı fiziksel kurallarla ele alınması ыкмаsel açıdan önemlidir. Özellikle энтальпия аралашмаının агым yönüne göre dinamik hesaplanması, өтмө olaylardaki агым tersine dönüşleriyle ылайыкdur.
Polinom насос модели yanında төрт аймактуу homolog modelin bulunması, çerçeveyi yalnızca normal işletme изилдөө noktasına değil, насос токтошу ve тескери агым koşullarına da taşımaktadır.
Python API, karmaşık контурlerin insan tarafından elle oluşturulması sırasında ortaya çıkabilecek yön, шилтеме ve dönüşüm hatalarını azaltabilecek pratik bir yazılım katmanı sunmaktadır.
Негизги чектөөлөр жана түшүндүрүлүшү керек болгон жагдайлар
En önemli kapsam sınırlaması бир фазалуу formülasyondur. Nükleer santral güvenlik analizlerinde kaynama, буулануу, yoğuşma, iki fazlı басым 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 жыйынтык sunmamaktadır.
Analitik doğrulamalar düşük hata vermiştir; ancak ayrıntılı ağ ve zaman adımı yагымnsama таблицаları yalnızca sınırlı biçimde açıklanmıştır. Birinci problemde 1090 hücrenin uzaysal yагымnsama için yeterli olduğu belirtilmiş, fakat farklı ağ çözünürlüklerindeki жыйынтыктар gösterilmemiştir.
TRACE karşılaştırmasında туруктуу абал için sayısal farklar sunulurken өтмө durumlarda баалоо büyük ölçüde grafiklerin üst üste gelmesine dayanmaktadır. Tepe маани farkı, kök орточо kare hata veya olay zamanlaması hatası gibi ölçüler verilmemiştir.
Tablo 3’te bildirilen kapalı çevrim чыгымlerinin birimi veya ölçeği, aynı bölümdeki geometri ve насос eğrisiyle birlikte yeniden kontrol edilmelidir. Bu nokta analitik ve FFN мааниlerinin birbirine yагымn olmasını ortadan kaldırmaz; ancak her iki мааниin fiziksel birimi doğru verilmediyse doğrulama sonucunun yeniden yorumlanması gerekir.
Изилдөө OpenFOAM tabanlı bir yazılım geliştirmesi olmasına rağmen kullanılan yazılım sürümü, булак 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.
Изилдөөнүн ыкмасы жана табылгалары
Моделдин камтышы
| Özellik | Изилдөөdaki uygulama |
|---|---|
| Temel platform | OpenFOAM tabanlı foamForNuclear |
| Boyut | Bir boyutlu түтүк ve sistem ağı |
| Akış rejimi | Tek fazlı |
| Korunan büyüklükler | Kütle, импульс ve энергия/энтальпия |
| Geometrik bileşenler | Borular, dirsekler, байланыштар, dallanmalar ve çap değişimleri |
| Hidrolik bileşenler | Yerel kayıplar, насостор ve zamanla değişen басым sıçramaları |
| Termal bileşenler | Isı булакları, жылуулук alıcıları ve байланыштарda ideal аралашма |
| Devre oluşturma | foamForNuclear Python API |
Иштелип чыккан программалык компоненттер
| Bileşen | Görevi | Temel sınırlaması |
|---|---|---|
branchConnector1D | Alan bazında çevrimsel veya sıfır gradyan шилтемеsı oluşturmak | Çok kollu шилтемеda diğer kolların katkжылуулук yarı örtük hesaplanmaktadır |
| Kütle koruyan ылдамдык шилтемеsı | Farklı kesitlerde \(\alpha\rho AU\) мааниini korumak | Tek fazlı mevcut formülasyona колдонулган |
pressureCoupling1D | Dinamik басым, kayıp, насос ve hidrostatik sıçramayı uygulamak | Kayıp katsayılarının doğruluğu kullanıcı girdilerine bağlıdır |
| Polinom насос модели | Normal изилдөө noktasında \(h(Q)\) eğrisi oluşturmak | Sıfır ve тескери агымı temsil edemez |
| Homolog насос модели | İleri/ters dönme ve ileri/тескери агымı temsil etmek | Doğru homolog baş ve tork таблицаları gerektirir |
| Ataletli насос токтошу | Motor torku kesildikten sonra ротор ылдамдыкını hesaplamak | Atalet ve сүрүлүү modeli kullanıcı tarafından tanımlanmalıdır |
| Entalpi аралашмаı | Bağlantıya giren агымmları массалык чыгымyle ağırlıklandırmak | Bağlantı içinde ideal ve anlık аралашма varsaymaktadır |
| Python түтүк API’si | Ağ, yön ve чек шарттарыnı otomatik oluşturmak | Изилдөөda tam giriş dosyaları ve sürüm kimliği paylaşılmamıştır |
Аналитикалык текшерүү матрицасы
| Doğrulama problemi | Sınanan касиет | Ağ veya koşul | En yüksek bildirilen hata |
|---|---|---|---|
| Akış bölünmesi ve birleşmesi | Bağlantı масса сакталууu ve сүрүлүү 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ı түтүк | <%0,05 |
| Doğal айланууlı çevrim | Sıcaklığa bağlı тыгыздык ve kaldırma kuvveti | 300–350 K температура farkı | <%0,05 |
| Çoklu dallanma ve ысытуу | Akış dağılımı ve энтальпия аралашмаı | 30 MW жылуулук kaynağı | %0,16 |
TRACE салыштыруу моделинин негизги өзгөчөлүктөрү
| Alt sistem | Modelleme yaklaşımı |
|---|---|
| Reaktör контурsi | İki çevrimli basitleştirilmiş PWR benzeri birincil sistem |
| Çekirdek | Toplulaştırılmış yагымt çubuğu ve бир өлчөмдүү radyal iletim |
| Isı transferi | Eşlenik yагымt-soğutucu aktarımı ve El-Genk bağıntжылуулук |
| Buhar üreteçleri | Sabit ikincil taraf sıcaklığına sahip жылуулук alıcıları |
| Sürtünme | Churchill bağıntжылуулук 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 насосsının durması |
| Karşılaştırılan senaryolar | Sıfır güç ve nominal güç |
Ana nicel табылгалар
| Bulgular | Sonuç | Yorum sınırı |
|---|---|---|
| Analitik hidrolik doğrulama | Bildirilen hatalar %0,33 veya altında | Seçilmiş idealize problemler |
| Enerji аралашмаı | Çıkış sıcaklığı hatası %0,03 | Sabit касиет ve ideal аралашма божомолı |
| TRACE nominal туруктуу абал | FFN çekirdek чыгымsi %2,3 daha düşük | İki kodun bağıntıları tamamen aynı değildir |
| Asimetrik насос koşulları | Toplam чыгым farkı %3’ün altında | Yalnızca dört ылдамдык kombinasyonu |
| Sıfır güçlü насос токтошу | Akış yavaşlaması ve tersine dönüş görsel olarak yагымn | Sayısal өтмө durum hata metriği verilmemiştir |
| Nominal güçlü geçiş | Debi, насос ылдамдыкı, температура ve жылуулук transfer eğrileri yагымn | Tek өтмө olay ve tek реактор düzeni |
Кайталануучулук жана сандык ишенимдүүлүк
Изилдөө temel denklemleri, чек шарттарыnı, analitik karşılaştırma мааниlerini ve реактор контуруnin ana model божомолlarını ayrıntılı biçimde vermektedir. Bu durum ыкмаsel баалооyi 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 түтүктөрın çap, uzunluk, pürüzlülük ve kayıp katsayıları,
- Zaman adımı ve doğrusal чечүүчү toleranslarının tam listesi,
- Geçici durum ağ ve zaman adımı duyarlılığı,
- Duvar мөөнөтsi, CPU kullanımı ve paralel ölçeklenebilirlik,
- TRACE karşılaştırması için nicel өтмө durum hata ölçüleri,
- Analitik testlerin otomatik çalıştırılabileceği giriş dosyaları.
Deterministik hesaplama изилдөөsı olduğu için тажрыйбаsel 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 мааниlendirilmelidir. Bu analizler mevcut изилдөөda sınırlı düzeyde sunulmuştur.
Булак жана метод эскертүүсү
Изилдөөnı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: Изилдөөda eş katkı veya eş birinci yazarlık açıklaması bulunmamaktadır.
Sorumlu veya iletişim yazarı: Yüklenen изилдөөda sorumlu yazarı belirten yıldız veya e-posta adresi yer almamaktadır. SSRN kayıt sayfası Giovanni Nervi’yi “Contact Author” olarak көрсөтөт.
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 билдирилген.
DOI: 10.2139/ssrn.7197741. Bu DOI, SSRN ön baskı kaydına aittir ve рецензентli bir dergi makalesi DOI’si değildir.
Dergi veya konferans: İncelenen sürüm için doğrulanmış bir рецензентli dergi veya konferans yayını belirtilmemiştir.
Yayın platformu: SSRN.
Özgün yayınevi: Изилдөө 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, бир фазалуу sistem termal-hidroliği ve bilimsel yazılım geliştirme araştırma препринтi.
Hakemlik durumu: Изилдөө рецензиядан өткөн эмес. Her sayfada “This препринт research paper has not been peer reviewed” ve “Preprint not peer reviewed” uyarıları bulunmaktadır.
Resmî байланыштар:SSRN расмий kayıt sayfası ve SSRN DOI шилтемеsı.
Finansman: Изилдөө, 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, ыкма, yazılım, doğrulama, görselleştirme ve ilk taslak; Alessandro Scolaro kavramsallaştırma, biçimsel analiz, ыкма, 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 билдирилген.
Veri ve kod erişimi: İncelenen изилдөөda doğrulama giriş dosyaları, kullanılan yazılım commit kimliği, TRACE modeli veya жыйынтыктарın ham sayısal verileri için açık bir veri deposu belirtilmemiştir. Python API’den örnek kod verilmiştir; ancak bu örnek изилдөөnın tamamını bağımsız olarak yeniden üretmek için yeterli değildir.
Bu Türkçe makale, yüklenen 27 sayfalık изилдөөnın ыкма açıklamaları, 27 temel denklemi, dört ana жыйынтык таблицаsu, T-шилтеме ve контур şemaları, температура ve ылдамдык dağılımları, TRACE düğümlemesi, гомологдук насос eğrileri ve насос токтошу geçiş grafikleri esas alınarak hazırlanmıştır. Изилдөөda bulunmayan yeni bir sayısal жыйынтык veya dış булакtan bilimsel performans bulgusu eklenmemiştir. Dış булак 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.
Изилдөөnın başlıca чектөөлөрı; рецензент баалооsinden geçmemiş olması, modelin yalnızca бир фазалуу агымlarla sınırlı kalması, fiziksel тажрыйбаsel validasyon bulunmaması, yalnızca bir basitleştirilmiş реактор контуруnin TRACE ile karşılaştırılması, өтмө 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 чыгым birimi veya ölçeğinin açıklamaya ihtiyaç göstermesidir.
Sonuçlar, geliştirilen çerçevenin seçilmiş бир фазалуу analitik testlerde сакталуу denklemlerini doğru uyguladığını ve basitleştirilmiş насос токтошу geçişlerinde TRACE ile benzer eğilimler üretebildiğini көрсөтөт. Bulgular, gerçek bir nükleer tesis için lisanslama yeterliliği, iki fazlı kaza analizi, тажрыйбаsel validasyon veya düzenleyici kullanım onayı anlamına gelmemektedir.
Пикир калтырыңыз
E-mail дарегиңиз жарыяланбайт. Милдеттүү талаалар * менен белгиленген