
Bir nükleer reaktörün çekirdeğindeki üç boyutlu güç dağılımının doğru bilinmesi; yakıt tükenmesinin, soğutucu sıcaklık alanının ve termal-hidrolik güvenlik sınırlarının izlenmesi açısından temel öneme sahiptir. Bu araştırma, altıgen yakıt demetlerine sahip küçük bir basınçlı su reaktörü modelinde seyrek sayıdaki çekirdek içi dedektörden elde edilen sinyaller kullanılarak üç boyutlu güç dağılımının yeniden oluşturulmasını incelemektedir. Araştırmacılar, yapılandırılmış dikdörtgen ağlara gereksinim duymadan düzensiz geometriler üzerinde çalışabilen thin-plate spline (TPS) yöntemini çekirdek içi dedektörlerin gerçek uzaysal tepki fonksiyonlarıyla birleştirmiştir.
OpenMC v0.14.0 Monte Carlo nötron taşınım koduyla oluşturulan sayısal model 200 MW termal güce, 583 altıgen yakıt demetine ve demet başına 61 yakıt çubuğuna sahiptir. Çekirdek içine dokuz dedektör demeti yerleştirilmiştir. Hesaplamalar 0–30 MWd/kgU arasındaki çeşitli tükenme seviyelerinde gerçekleştirilmiş; gerçek çalışma koşullarındaki olası model-belirsizliklerini temsil etmek amacıyla yakıt zenginleştirmeleri değiştirilmiş, ayrıca maksimum %10 mekânsal pertürbasyon ve ±%1 rastgele gürültü içeren ek testler yapılmıştır.
Dedektör tepki hesapları, bir dedektörün yalnız bulunduğu yakıt demetindeki gücü ölçmediğini göstermektedir. Tepki kuvvetli biçimde yerelleşmiş olsa da komşu demetler de ölçüme katkı vermektedir. Radyal yönde etkinin büyük bölümü yaklaşık 15 cm içindeki ilk üç komşuluk halkasında yoğunlaşmaktadır ve daha uzaktaki katkılar ihmal edilebilir düzeye inmektedir. Dedektör #5 için demet-ortalama tepkinin yaklaşık %44,1'i dedektörün bulunduğu merkez demetten gelirken yakın komşular da ölçülebilir katkılar sağlamaktadır.
Bu fiziksel tepki ihmal edilerek TPS uygulanması, araştırmadaki bütün referans grupları ve tükenme noktalarında hatayı azaltmak yerine büyütmüştür. Dedektör tepki fonksiyonu TPS denklemine dahil edildiğinde ise rekonstrüksiyon doğruluğu belirgin biçimde iyileşmiştir. Orta eksenel katmanda düşük tükenmede maksimum göreli sapma yaklaşık %8'den %3'e, yüksek tükenmede yaklaşık %5'ten %4'e düşürülmüştür. Yapay pertürbasyonlu daha zor testlerde de yeniden yapılandırılmış maksimum göreli hata düşük ve yüksek tükenme koşullarında yaklaşık %4 düzeyinde tutulmuştur.
Türkiye açısından anlamı: Araştırmadaki 200 MW termal küçük reaktör gerçek bir Türkiye reaktörünün mühendislik modeli değildir; yazarlar da çekirdek tasarımının algoritma doğrulaması amacıyla basitleştirilmiş olduğunu açıkça belirtmektedir. Dolayısıyla yakıt düzeni, dedektör yerleşimi, hata değerleri veya güvenlik sonuçları Türkiye'deki herhangi bir reaktöre doğrudan aktarılamaz. Buna karşılık çalışma, gelecekte küçük modüler reaktör veya gelişmiş reaktör teknolojileri için çevrim içi çekirdek izleme sistemleri geliştirilecekse dedektör sinyalinin yalnız yerel güç değeri olarak değil, çevresindeki nötron alanının uzaysal ağırlıklı toplamı olarak modellenmesinin önemini göstermektedir. Gerçek uygulama için reaktöre özgü nötronik model, gerçek dedektör geometrisi, kalibrasyon verisi, ölçüm belirsizliği ve düzenleyici doğrulama ayrıca gereklidir.
Neden küçük ve altıgen çekirdeklerde güç rekonstrüksiyonu daha zor?
Küçük reaktörler, geleneksel büyük güç reaktörlerine kıyasla daha yüksek yüzey/hacim oranına sahip olduklarından çekirdek çevresindeki nötron kaçağının göreli etkisi büyüyebilir. Aynı zamanda altıgen yakıt demetleri, kare yakıt demetleri için geliştirilmiş birçok geleneksel interpolasyon yönteminin dayandığı ortogonal ağ varsayımıyla doğrudan uyuşmaz.
Çekirdeğin küçük olması bir başka problemi de büyütmektedir: dedektör ile komşu yakıt demetleri arasındaki mesafeler azalır. Bu nedenle bir çekirdek içi dedektörün sinyali yalnız içinde bulunduğu yakıt demetinin gücünü temsil etmez; komşu demetlerdeki fisyon nötronları da dedektör sinyaline katkı sağlar.
Araştırmanın temel hipotezi, bu uzaysal dedektör etkisinin matematiksel rekonstrüksiyon sistemine açık biçimde dahil edilmesinin TPS'nin doğruluğunu artıracağıdır.
Thin-plate spline yöntemi ne yapıyor?
TPS yönteminde ölçülen dedektör sinyali ile nötronik hesaplamadan elde edilen sinyal arasındaki oran, çekirdek boyunca değişen bir düzeltme fonksiyonu olarak ele alınmaktadır. Çalışmada sabit nötron enerji spektrumu altında güç, nötron akısı ve dedektör ölçümlerinin birbirleriyle orantılı olduğu varsayılmıştır.
Temel oran ilişkisi:
\[ \frac{P^{M}(\mathbf r)}{P^{C}(\mathbf r)} = \frac{I^{M}(\mathbf r)}{I^{C}(\mathbf r)} = F(\mathbf r) \]
biçimindedir. Burada \(P^{C}\) nötronik hesaplamadan gelen güç dağılımını, \(P^{M}\) gerçek durumu temsil eden güç dağılımını, \(I^{C}\) hesaplanan dedektör sinyalini, \(I^{M}\) ölçülen dedektör sinyalini ve \(F\) TPS ile oluşturulan düzeltme yüzeyini temsil eder.
Yeniden oluşturulan güç:
\[ P^{Rec}(\mathbf r) = P^{C}(\mathbf r)\,F(\mathbf r) \]
yaklaşımıyla belirlenmektedir.
TPS yüzeyi nasıl tanımlanıyor?
Üç boyutlu TPS fonksiyonu bir küresel eğilim bileşeni ile dedektörlerin çevresindeki yerel düzeltmeleri temsil eden çekirdek fonksiyonlarının toplamından oluşturulmaktadır:
\[ F(x,y,z) = T(x,y,z) + \sum_{n=1}^{N_d} b_n K(x,y,z,x_{d,n},y_{d,n},z_{d,n}) \]
Küresel eğilim:
\[ T(x,y,z)=a_0+a_1x+a_2y+a_3z \]
olarak tanımlanmıştır.
Çalışmada kullanılan TPS çekirdek fonksiyonu:
\[ K_n=l_n^2\log l_n \]
şeklindedir ve \(l_n\), yeniden oluşturulacak nokta ile n'inci dedektör arasındaki üç boyutlu Öklid uzaklığıdır:
\[ l_n= \sqrt{ (x-x_{d,n})^2+ (y-y_{d,n})^2+ (z-z_{d,n})^2 } \]
TPS'nin önemli avantajı, yakıt demetlerinin düzenli bir Kartezyen ağ üzerinde bulunmasını zorunlu kılmamasıdır. Bu özellik altıgen demet geometrisini doğal biçimde işleyebilmesini sağlar.
Dedektör aslında ne ölçüyor?
Kaynakta dedektör sinyalinin nötron akısıyla:
\[ S=C\phi \]
şeklinde orantılı olduğu kabul edilmektedir. Ancak bir çekirdek içi dedektörün değeri tek bir geometrik noktadaki güç değildir. Dedektör sinyali bütün çekirdekten gelen katkıların uzaysal tepki fonksiyonuyla ağırlıklandırılmış integralidir:
\[ S= \int_V P(\mathbf r)\, w(\mathbf r)\, d\mathbf r \]
Çekirdek düğümlere ayrıldığında:
\[ S=\sum_i w_iP_i \]
haline gelir. Dolayısıyla aynı dedektör sinyalini oluştururken dedektöre yakın bir yakıt bölgesi ile daha uzak bir yakıt bölgesi aynı ağırlığa sahip değildir.
Radyal ve eksenel dedektör tepki fonksiyonları
Araştırmacılar uzaysal tepkiyi:
\[ w(\mathbf r)=R(x,y)A(z) \]
şeklinde radyal ve eksenel bileşenlere ayırmıştır. \(R(x,y)\) dedektörün yatay düzlemde çevresindeki yakıt demetlerine duyarlılığını, \(A(z)\) ise eksenel yöndeki duyarlılığı temsil etmektedir.
Bu ayrıştırma sayesinde ölçülen sinyal üç boyutlu güç dağılımıyla:
\[ I^M \propto \sum_{x,y,z} R(x,y)A(z)P(x,y,z) \]
şeklinde ilişkilendirilebilmektedir.
Dedektör tepki fonksiyonu nasıl hesaplandı?
Çalışmada ileri nötron taşınım yöntemi kullanılmıştır. Her yakıt konumunun bir nötron kaynağı gibi davranması durumunda dedektörde oluşturduğu fisyon reaksiyon hızı hesaplanmış ve bu değer uzaysal tepki fonksiyonunun temeli olarak alınmıştır.
Hesaplamalar açık kaynak OpenMC v0.14.0 Monte Carlo koduyla gerçekleştirilmiştir. Yazarlar kullandıkları sürümde doğrudan gerekli sabit-kaynak düzenini gerçekleştirmek amacıyla yakıt bölgelerindeki fisyon kanallarını yapay biçimde kapatan bir kesit değiştirme yaklaşımı kullanmıştır. İlgili fisyon kesitleri sanal yakıtlarda yakalama kesitlerine taşınmıştır.
Bu işlem gerçek reaktörde fiziksel olarak kesitlerin değiştirildiği anlamına gelmez; dedektörün her kaynak konumuna olan matematiksel/Monte Carlo tepkisini ayrı ayrı çıkarmak için kullanılan hesaplama tekniğidir.
Verianla Live: Dedektör tepkisinden 3B güç rekonstrüksiyonuna
Süreç, çalışmanın Bölüm 2.1–2.2'de tanımladığı TPS ve dedektör tepki yaklaşımının sıralı özetidir.
| Aşama | İşlem | Kaynak |
|---|---|---|
| 1. Nötronik temel çözüm | OpenMC ile çekirdeğin hesaplanan üç boyutlu güç dağılımı elde edilir. | Bölüm 2.3 |
| 2. Dedektör tepki hesabı | Her kaynak bölgesinin çekirdek içi dedektör sinyaline katkısı ileri taşınım hesabıyla belirlenir. | Bölüm 2.2 ve 3.1 |
| 3. Hesaplanan dedektör sinyali | Güç dağılımı radyal ve eksenel tepki fonksiyonlarıyla ağırlıklandırılır. | Denklem 17 |
| 4. TPS düzeltme yüzeyi | Ölçülen ve hesaplanan dedektör sinyallerinden üç boyutlu düzeltme fonksiyonu oluşturulur. | Denklem 18 |
| 5. 3B güç rekonstrüksiyonu | Hesaplanan güç dağılımı TPS düzeltme fonksiyonuyla güncellenir. | Bölüm 2.1–2.2 |
Küçük reaktör modeli nasıl oluşturuldu?
Çalışmada kullanılan çekirdek, gerçek bir lisanslanmış reaktör tasarımı değil, algoritmanın test edilmesi amacıyla tipik küçük PWR özelliklerini koruyan basitleştirilmiş bir modeldir. Yazarlar bu ayrımı açıkça vurgulamaktadır.
| Çekirdek parametresi | Kaynakta verilen değer |
|---|---|
| Termal güç | 200 MW |
| Çekirdek çapı | 212,5 cm |
| Çekirdek yüksekliği | 180 cm |
| Yakıt demeti sayısı | 583 |
| Demet adımı | 7,5 cm |
| Demet karşılıklı düz yüz mesafesi | 5,8 cm |
| Demet başına yakıt çubuğu | 61 |
| Yakıt çubuğu adımı | 0,75 cm |
| Soğutucu | Su |
| Soğutucu sıcaklığı | 280–320 °C |
| Reflektör | Su |
Kaynak içi not: Tablo 1'de “cladding thickness” 5 cm olarak verilmekte, Tablo 2'de yakıt çubuğu kaplama kalınlığı 0,05 cm olarak tanımlanmaktadır. İlk değerin hangi yapısal kaplamayı temsil ettiği yeterince açıklanmadığından burada yeni bir mühendislik anlamı atanmadı.
Yakıt çubuklarının özellikleri
| Parametre | Değer |
|---|---|
| Yakıt çubuğu uzunluğu | 160 cm |
| Dış çap | 0,6 cm |
| Yakıt | UO2 |
| Yakıt sıcaklığı | 900 K |
| Zenginleştirme | %3,5 ve %4,5 |
| Yakıt yoğunluğu | 10,42 g/cm³ |
| Kaplama | Zirkonyum alaşımı |
| Kaplama sıcaklığı | 600 K |
| Kaplama yoğunluğu | 6,55 g/cm³ |
| Yakıt çubuğu kaplama kalınlığı | 0,05 cm |
| Yanabilir zehir | Gd2O3 |
| Gd2O3 kütle oranı | %8 |
Üç farklı yakıt çubuğu tipi kullanılmıştır: %3,5 zenginleştirilmiş UO2, %4,5 zenginleştirilmiş UO2 ve %8 Gd2O3 içeren %4,5 zenginleştirilmiş UO2. Bunlar beş halkalı, 61 yakıt çubuklu altıgen demetler halinde düzenlenmiştir.
Dokuz çekirdek içi dedektör nasıl yerleştirildi?
Çalışmada çekirdeğin farklı bölgelerini temsil etmek üzere 9 çekirdek içi dedektör demeti kullanılmıştır. Her dedektör ilgili yakıt demetinin merkezine yerleştirilmiştir.
Şekil 4, dedektör konumlarının tüm çekirdek üzerinde seyrek bir örnekleme ağı oluşturduğunu göstermektedir. Bu düzen daha sonraki sonuçların yorumunda önemlidir; çünkü çekirdek kenarında ve özellikle dedektör kapsamasının zayıf olduğu bölgelerde TPS hatalarının sistematik olarak yükseldiği gözlenmiştir.
Tükenme hesabı hangi noktaları kapsıyor?
OpenMC hesapları aşağıdaki tükenme seviyelerinde gerçekleştirilmiştir:
| Hesap noktası | Tükenme (MWd/kgU) |
|---|---|
| 0 | 0 |
| 1 | 0,01 |
| 2 | 0,11 |
| 3 | 1 |
| 4 | 3 |
| 5 | 5 |
| 6 | 8 |
| 7 | 10 |
| 8 | 15 |
| 9 | 20 |
| 10 | 30 |
Şekil 5'te Gd2O3 içeren yakıtın başlangıçtaki güçlü nötron soğurumu nedeniyle düşük tükenmede etkin çoğalma katsayısının bastırıldığı, gadolinyum tükendikçe bu etkinin azaldığı görülmektedir. Kaynak yaklaşık 32 MWd/kgU civarında \(k_{\mathrm{eff}}=1\) elde edildiğini belirttiğinden güç rekonstrüksiyonu için üst sınırı 30 MWd/kgU olarak seçmiştir.
Referans “gerçek” güç dağılımı nasıl oluşturuldu?
Modelin kendi nominal hesaplamasını yine kendisiyle karşılaştırmak yerine üç farklı zenginleştirme pertürbasyonu uygulanmıştır:
- %3,5 zenginleştirilmiş yakıt demetlerinin zenginleştirilmesi kaynak ifadesiyle “%0,2 azaltılmıştır”.
- Yanabilir zehir içeren %4,5 zenginleştirilmiş demetlerde aynı tür azaltma uygulanmıştır.
- Bütün yakıt demetlerinde aynı azaltma uygulanmıştır.
Yorum sınırı: Kaynak “reduced by 0.2% enrichment” ifadesini kullanmaktadır; bunun 0,2 yüzde puanlık mutlak bir değişim mi yoksa nominal değerin göreli %0,2'si mi olduğu açıkça tanımlanmamıştır. Bu nedenle burada matematiksel olarak yeniden yorumlanmamıştır.
Neden ayrıca %10 pertürbasyon ve ±%1 gürültü eklendi?
Nominal model ile zenginleştirilmiş referans çözümleri arasındaki fark bazı durumlarda zaten oldukça küçük ve düzgün dağılımlı olduğundan TPS'nin iyileştirme potansiyelini açık biçimde göstermek zorlaşmaktadır. Araştırmacılar bunun için daha zor bir sentetik test oluşturmuştur.
Referans çözüme x doğrultusunda dış kenarda en fazla %10'a ulaşan sistematik bir uzaysal pertürbasyon ve buna ek ±%1 rastgele gürültü uygulanmıştır.
Bu test gerçek bir ölçülmüş reaktör arızasını temsil etmemektedir; algoritmanın daha büyük ve belirgin uzaysal sapmalar altında nasıl davrandığını sınamak için oluşturulan sayısal stres testidir.
Rekonstrüksiyon doğruluğu nasıl ölçüldü?
Çalışma üç temel hata ölçütü kullanmaktadır: maksimum göreli hata, RMSE ve yüksek güç bölgelerine daha dengeli ağırlık vermek üzere tanımlanan normalize edilmiş RMSE.
Rekonstrüksiyon sonrası göreli hata:
\[ e_i^R= \frac{P_i^{Rec}-P_i^R}{P_i^R} \]
şeklindedir.
RMSE:
\[ \varepsilon^R= \sqrt{ \frac{1}{N} \sum_{i=1}^{N} \left( \frac{P_i^{Rec}-P_i^R}{P_i^R} \right)^2 } \]
ile hesaplanmaktadır.
Çalışmada ayrıca düşük güçlü çevre düğümlerinin göreli hatasının genel metriği aşırı etkilemesini azaltmak amacıyla ortalama referans gücüne göre normalize edilmiş NRMSE kullanılmıştır.
Dedektör #5 çevresindeki gerçek tepki ne kadar yerel?
Şekil 6 ve Şekil 7 araştırmanın en önemli fiziksel sonuçlarından birini göstermektedir. Dedektör #5'in tepki fonksiyonu merkezden uzaklaşıldıkça hızlı biçimde azalmaktadır.
Demet-ortalama radyal grafikte dedektörün bulunduğu merkez demetin katkısı %44,1 olarak gösterilmektedir. En yakın komşu demetlerde katkılar yaklaşık %6 mertebesindeyken daha dıştaki halkalarda yaklaşık %1–2 düzeyine gerilemektedir.
Araştırmacılar toplam pratik etki alanını yaklaşık 15 cm içerisindeki üçüncü komşuluk halkası ile sınırlandırmaktadır. Daha uzaktaki demetlerin katkısı ihmal edilebilir kabul edilmiştir.
Bu sonuç, “dedektör okuması = bulunduğu demetin gücü” varsayımının neden yetersiz olabileceğini doğrudan göstermektedir.
Eksenel tepki nasıl değişiyor?
Şekil 8'de dedektörün eksenel duyarlılığı, kaynak demetinin merkezi, ikinci halka ve üçüncü halkada bulunması için karşılaştırılmıştır.
Kaynak dedektörün bulunduğu merkez demetteyse, aynı eksenel katmanda keskin ve yüksek bir tepki tepesi oluşmaktadır. Kaynak farklı eksenel katmana taşındığında tepki hızla azalmaktadır.
Kaynak radyal olarak ikinci ve üçüncü halkalara uzaklaştıkça eksenel tepki genişlemekte ve düzleşmektedir. Buna rağmen en güçlü sinyal yine kaynakla aynı eksenel seviyenin çevresinde görülmektedir.
Dedektör tepki fonksiyonu kullanılmazsa ne oluyor?
Şekil 9'da TPS, dedektörün uzaysal entegrasyon etkisi hesaba katılmadan uygulanmıştır. Üç referans grubunun tamamında ve incelenen bütün tükenme noktalarında:
- maksimum hata oranı,
- RMSE oranı,
- NRMSE oranı
rekonstrüksiyon sonrasında 1'in üzerinde kalmaktadır.
Başka bir ifadeyle yöntem hatayı düzeltmek yerine büyütmektedir. Bunun nedeni TPS'nin noktasal güç değerini, gerçekte çevredeki birçok düğümden gelen ağırlıklı nötron katkısını içeren dedektör sinyaliyle doğrudan eşleştirmesidir.
Bu sonuç araştırmanın ana mühendislik mesajlarından biridir: ileri düzey bir interpolasyon algoritması kullanmak tek başına yeterli değildir; interpolasyona verilen ölçümün fiziksel oluşum mekanizmasının da doğru modellenmesi gerekir.
Dedektör tepkisi hesaba katılınca ne değişiyor?
Şekil 10'da dedektör tepki fonksiyonları TPS denklemine dahil edildiğinde hata metriklerinin büyük bölümü 1'in altına düşmektedir. Bu, rekonstrüksiyon sonrasındaki farkın başlangıç sapmasından daha küçük olduğu anlamına gelir.
Yüksek tükenme seviyelerinde bazı global hata oranlarının yeniden 1 civarına veya biraz üzerine çıkabilmesinin nedeni kaynak tarafından güç dağılımının giderek daha düz hale gelmesi ve başlangıç sapmasının zaten küçük olmasıyla açıklanmaktadır. Böyle bir durumda TPS'nin düzeltebileceği marj da küçülmektedir.
Düşük ve yüksek tükenmede maksimum hata nasıl değişti?
Orta eksenel katman olan 6. katmanın radyal haritalarında:
- 0 MWd/kgU: maksimum başlangıç sapması yaklaşık %8 → rekonstrüksiyon sonrası yaklaşık %3,
- 30 MWd/kgU: maksimum başlangıç sapması yaklaşık %5 → rekonstrüksiyon sonrası yaklaşık %4.
Dedektör dağılımı neden hâlâ sınırlayıcı?
Şekil 11 ve 12'de özellikle çekirdeğin alt-sol çevresinde diğer bölgelere göre daha büyük yerel hatalar kalmaktadır. Bu bölge, dokuz dedektörün kapsamasının daha seyrek olduğu kısımlardan biridir.
TPS ölçüm noktaları arasındaki yüzeyi interpolasyonla oluşturduğundan, bir noktaya yakın yeterli dedektör bulunmaması doğruluğu sınırlar. Çekirdek kenarındaki yakıt demetlerinin fiziksel olarak daha az komşuya sahip olması da bu sorunu güçlendirmektedir.
Bu nedenle araştırmanın sonucu “dokuz dedektör her durumda yeterlidir” değildir. Tam tersine sonuçlar dedektör yerleşiminin ayrıca optimize edilmesi gerektiğine işaret etmektedir.
Yapay pertürbasyon testinde ne oldu?
Daha zor testte başlangıç dağılımına %10'a kadar sistematik mekânsal pertürbasyon ve ±%1 rastgele gürültü eklenmiştir.
Birinci referans grubu için bütün üç boyutlu çekirdekte maksimum başlangıç sapması düşük tükenmede yaklaşık %18, yüksek tükenmede yaklaşık %14 düzeyindedir. RMSE ve NRMSE başlangıçta yaklaşık %6–8'e çıkmakta ve yüksek tükenmede yaklaşık %5 civarında seyretmektedir.
TPS ve dedektör tepki fonksiyonları birlikte kullanıldığında maksimum hata başlangıç değerinin yaklaşık %30–50'sine, RMSE ve NRMSE ise yaklaşık %20–25'ine indirilebilmiştir.
Şekil 16–17'de yalnız 6. katmanın radyal dağılımına bakıldığında başlangıç maksimum sapması 0 MWd/kgU'da yaklaşık %15, yüksek tükenmede yaklaşık %11'dir. Bu değerlerin önceki %18 ve %14 değerlerinden farklı olması bir çelişki değildir; önceki değerler daha geniş veri kümesindeki maksimumları, Şekil 16–17 ise seçilmiş orta katmandaki radyal dağılımı göstermektedir.
Her iki seçilmiş tükenme durumunda da TPS sonrasında orta katmanın maksimum göreli hatası yaklaşık %4 seviyesine indirilmiştir.
%8 sınırı nasıl yorumlanmalı?
Yazarlar bütün rekonstrüksiyon sonrası göreli hataların National Nuclear Safety Administration'a atfettikleri %8 hata kriterinin içinde kaldığını belirtmektedir.
Ancak çalışmanın bu ifade için verdiği [32] numaralı kaynak doğrudan bir düzenleyici standart veya resmî NNSA dokümanı değildir; 2011 yılında Nuclear Engineering and Design dergisinde yayımlanmış bir güç dağılımı rekonstrüksiyonu çalışmasıdır. Bu nedenle burada %8 değeri bu araştırmanın kullandığı kabul kriteri olarak aktarılmakta, bağımsız olarak doğrulanmış güncel bir mevzuat sınırı olduğu ileri sürülmemektedir.
Çalışmanın desteklediği sonuçlar
- TPS, altıgen yakıt demetleri gibi yapılandırılmamış geometrilerde üç boyutlu güç rekonstrüksiyonu için kullanılabilir.
- Çekirdek içi dedektör sinyali güçlü biçimde yerel olsa da tek bir yakıt demetinin gücünü temsil etmez.
- Dedektör #5 için en büyük radyal katkı dedektörün bulunduğu demetten gelmektedir.
- Dedektörün anlamlı radyal tepki alanı yaklaşık 15 cm içindeki komşu demetlerle sınırlanmaktadır.
- Dedektör uzaysal tepki fonksiyonunun ihmal edilmesi TPS rekonstrüksiyonunda sistematik hata büyümesine neden olmuştur.
- Tepki fonksiyonunun TPS'ye dahil edilmesi farklı tükenme seviyelerinde genel rekonstrüksiyon hatasını azaltmıştır.
- 0 MWd/kgU'daki orta katmanda yaklaşık %8 maksimum sapma yaklaşık %3'e düşürülmüştür.
- 30 MWd/kgU'daki yaklaşık %5 maksimum sapma yaklaşık %4'e indirilmiştir.
- Yapay pertürbasyon testlerinde seçilen orta katmanda maksimum hata yaklaşık %4 seviyesine indirilebilmiştir.
- Rekonstrüksiyon doğruluğu dedektörlerin uzaysal dağılımına güçlü biçimde bağlıdır.
Çalışmanın desteklemediği sonuçlar
- Kullanılan 200 MW termal çekirdek lisanslanmış veya mühendislik optimizasyonu tamamlanmış ticari bir SMR tasarımı değildir.
- Çalışma gerçek işletme halindeki bir reaktörden alınmış çevrim içi dedektör verileriyle doğrulanmamıştır.
- Yapay %10 pertürbasyon gerçek bir reaktör kazasını veya fiziksel arızayı temsil etmemektedir.
- Yöntemin bütün altıgen küçük reaktör tasarımlarında aynı hata düzeyini sağlayacağı gösterilmemiştir.
- Dokuz dedektörün optimal sayı veya optimal yerleşim olduğu kanıtlanmamıştır.
- Geçici rejimlerde ve hızlı güç değişimlerinde yöntem doğrulanmamıştır.
- Ölçüm, model ve nükleer veri belirsizliklerinin tam olasılıksal propagasyonu yapılmamıştır.
- Kaynağın kullandığı %8 kriteri bu çalışma içinde doğrudan güncel düzenleyici dokümanla doğrulanmamıştır.
Çalışmanın Yöntemi ve Bulguları
Araştırma tasarımı
| Bileşen | Uygulama |
|---|---|
| Reaktör türü | Basitleştirilmiş küçük PWR modeli |
| Yakıt demeti geometrisi | Altıgen |
| Termal güç | 200 MW |
| Yakıt demeti sayısı | 583 |
| Yakıt çubuğu / demet | 61 |
| Eksenel ayrıklaştırma | 10 katman |
| Çekirdek içi dedektör demeti | 9 |
| Nötron taşınım kodu | OpenMC v0.14.0 |
| Rekonstrüksiyon yöntemi | Thin-plate spline |
| Dedektör tepki hesabı | İleri nötron taşınımı |
| Tükenme aralığı | 0–30 MWd/kgU |
| Referans senaryosu | 3 zenginleştirme pertürbasyonu |
| Ek stres testi | %10'a kadar uzaysal pertürbasyon + ±%1 rastgele gürültü |
Dedektör tepki bulguları
| Gözlem | Kaynak sonucu |
|---|---|
| Merkez demet katkısı, dedektör #5 | %44,1 |
| Yakın komşular | Demet başına yaklaşık %6 düzeyi |
| Daha uzak halkalar | Yaklaşık %1–2 düzeyi |
| Etkin radyal tepki alanı | Yaklaşık üçüncü komşuluk halkasına / ~15 cm'ye kadar |
| Dördüncü halka ve ötesi | Kaynak tarafından ihmal edilebilir kabul edilmiştir |
| Eksenel davranış | Aynı eksenel düğümde güçlü pik; uzaklık arttıkça hızlı zayıflama |
Dedektör tepki fonksiyonu olmadan TPS
Üç referans grubu ve bütün tükenme noktalarında rekonstrüksiyon sonrası/öncesi:
- maksimum hata oranı > 1,
- RMSE oranı > 1,
- NRMSE oranı > 1
olmuştur. Bu test TPS'nin kendisinin başarısız olduğunu değil, TPS'ye fiziksel olarak yanlış temsil edilen dedektör verisi verildiğinde rekonstrüksiyonun bozulduğunu göstermektedir.
Dedektör tepki fonksiyonu ile TPS
Dedektör hesaplanan sinyalleri çevredeki güç dağılımının \(R(x,y)A(z)\) tepki fonksiyonuyla ağırlıklandırılmış toplamından üretildiğinde rekonstrüksiyon sonrası hata oranlarının büyük kısmı 1'in altına düşmüştür.
| Koşul | Başlangıç maksimum sapması | TPS sonrası maksimum hata |
|---|---|---|
| 0 MWd/kgU, katman 6 | Yaklaşık %8 | Yaklaşık %3 |
| 30 MWd/kgU, katman 6 | Yaklaşık %5 | Yaklaşık %4 |
| Yapay pertürbasyon, 0 MWd/kgU, katman 6 | Yaklaşık %15 | Yaklaşık %4 |
| Yapay pertürbasyon, yüksek tükenme, katman 6 | Yaklaşık %11 | Yaklaşık %4 |
Yöntemin başlıca sınırlılıkları
- Çekirdek mühendislik optimizasyonu tamamlanmış gerçek bir ticari tasarım değildir.
- Rekonstrüksiyon doğrulaması gerçek işletme dedektör ölçümlerinden ziyade Monte Carlo ve sentetik pertürbasyon verilerine dayanmaktadır.
- Dedektör dağılımının seyrek olduğu çevre bölgelerinde hatalar artmaktadır.
- Dedektör yerleşim optimizasyonu yapılmamıştır.
- Belirsizlik propagasyonu kapsamlı biçimde modellenmemiştir.
- Geçici çalışma koşulları bu araştırmanın kapsamı dışındadır.
- Kaynakta bazı geometrik parametrelerin tanımı ve zenginleştirme pertürbasyonunun yüzdesel anlamı belirsizdir.
Kaynak ve Yöntem Notu
Temel kaynak: Kullanıcı tarafından sağlanan 28 sayfalık çalışma.
Tam özgün başlık: Power reconstruction in hexagonal-assembly small reactor using thin-plate spline and in-core detector response
Yazarlar, kaynak sırasıyla: Hao Yang; Tianxiang Wang; Zhuo Li; Shengli Chen.
Sorumlu yazarlar: Zhuo Li — lizhuo@xpu.edu.cn; Shengli Chen — chenshli23@mail.sysu.edu.cn.
Eş katkı/eş birinci yazar: Kaynakta böyle bir beyan bulunmamaktadır.
Kurumlar: Sino-French Institute of Nuclear Engineering and Technology, Sun Yat-sen University, Zhuhai 519082, Guangdong, China; School of Science, Xi’an Polytechnic University, Xi’an, China.
Kaynak türü ve hakemlik durumu: Yüklenen çalışma SSRN üzerinde bulunan ve hakem değerlendirmesinden geçmemiş bir preprinttir.
SSRN Abstract ID: 6933221
DOI: 10.2139/ssrn.6933221
Resmî bağlantı: https://ssrn.com/abstract=6933221
Sürüm notu: Bibliyografik dış kontrolde aynı başlıkla 6933220 numaralı başka bir SSRN indeks sonucu da görünmektedir. Bununla birlikte yüklenen çalışmanın her sayfasında 6933221 belirtilmekte ve DOI kaydı 10.2139/ssrn.6933221 kimliğini doğrulamaktadır. Bu nedenle bilimsel içerik veya bibliyografik bilgiler iki kayıt arasında birleştirilmemiş; yüklenen 6933221 sürümü esas alınmıştır.
Hakemli dergi sürümü: Yapılan tam başlık ve DOI kontrolünde bu çalışma için ayrı bir hakemli dergi yayını doğrulanmamıştır.
Lisans: Yüklenen sürümde şekillerin yeniden kullanımına açık bir Creative Commons lisans beyanı tespit edilmemiştir. Bu nedenle özgün grafikler ve çekirdek haritaları burada bire bir kopyalanmamış; bulgular yeniden yazılmış metin, tablo ve Verianla Live yapılarıyla aktarılmıştır.
Finansman: Yüklenen 28 sayfalık metinde ayrı bir finansman beyanı tespit edilmemiştir.
Çıkar çatışması: Yüklenen sürümde ayrı bir çıkar çatışması beyanı tespit edilmemiştir.
Veri erişilebilirliği: Yüklenen sürümde ayrı bir veri erişilebilirliği beyanı tespit edilmemiştir.
CRediT katkı beyanı: Kaynakta ayrıntılı CRediT yazar katkıları verilmemiştir.
Kaynak içi teknik/editoryal belirsizlikler: Tablo 1'de “cladding thickness” 5 cm, Tablo 2'de yakıt çubuğu kaplama kalınlığı 0,05 cm olarak verilmiştir ve ilk parametrenin fiziksel kapsamı açık değildir. “Assembly flat-to-flat distance = 5,8 cm” Tablo 1'de iki kez tekrarlanmaktadır. Referans çözümlerindeki “0.2% enrichment reduction” ifadesinin göreli yüzde mi yoksa yüzde puan değişimi mi olduğu açıkça belirtilmemiştir.
%8 kriteri hakkında not: Çalışma, rekonstrüksiyon sonrası hataların National Nuclear Safety Administration tarafından öngörülen %8 sınırının içinde kaldığını belirtmektedir. Ancak bu ifade için kullanılan [32] numaralı referans doğrudan resmî düzenleyici doküman değil, Wang ve arkadaşlarının 2011 tarihli Nuclear Engineering and Design makalesidir. Bu nedenle %8 değeri burada çalışmanın kullandığı kabul kriteri olarak aktarılmış, bağımsız biçimde doğrulanmış güncel mevzuat hükmü olarak sunulmamıştır.
Bilimsel yorum sınırı: Çalışma, Monte Carlo nötron taşınımı ve sentetik referans/pertürbasyon senaryoları kullanarak TPS ile dedektör tepki fonksiyonlarının birleştirilmesini doğrulamaktadır. Sonuçlar gerçek bir ticari küçük reaktörün işletme performansını, lisanslanabilirliğini veya güvenlik marjlarını doğrudan kanıtlamamaktadır.
Verianla yöntem notu: Bilimsel açıklamalar yüklenen çalışmanın yöntem, denklem, tablo, şekil ve sayısal sonuçlarına dayanmaktadır. Dış araştırma yalnız bibliyografik kimliğin, SSRN DOI'sinin ve yayın durumunun kontrolü için kullanılmış; başka kaynaklardan yeni nükleer performans sonucu eklenmemiştir.

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