Fisyon Zincirleme Reaksiyonları ve Nükleer Reaktörler—Giriş
Bu bölüm, fisyonun mikroskobik fiziği (Bölüm 2) ile kendi kendini sürdüren nükleer reaktörlerin mühendisliği arasında köprü kurar. Reaktör analizinin en önemli kavramını tanıtır: çoğaltma faktörü k; iki ardışık fisyon neslindeki nötron sayılarının oranı olarak veya daha kullanışlı biçimde nötron üretiminin nötron kaybına (soğurma artı kaçak) oranı olarak tanımlanır. Kritiklik koşulları k<1 (kritik-altı), k=1 (kritik) ve k>1 (kritik-üstü) kurulur ve üstel nüfus artışı N(t)=N0 exp[(k-1)t/l] ile reaktör periyodu T=l/(k-1) gösteren bilinçli olarak basitleştirilmiş bir nokta-kinetik modeli türetilir. Bölüm, ~1e-4 s'lik bir ani (prompt) nötron ömrünün kontrolü neredeyse imkânsız kılacağını ortaya koyarak gecikmeli nötronların kararlılık sağlayan rolüne işaret eder. Ardından 'nötronu izle' yaşam-döngüsü dengesi geliştirilerek dört-faktör formülü k_sonsuz=eta*f*p*epsilon ve altı-faktör formülü k=eta*f*p*epsilon*P_FNL*P_TNL elde edilir; termal kullanım f, üreme faktörü eta, rezonanstan kaçış olasılığı p, hızlı fisyon faktörü epsilon ve hızlı/termal kaçmama olasılıkları tanımlanır. Hızlı nötronların büyük termal fisyon tesir kesitlerinden yararlanmak için yavaşlatılmasında moderasyonun rolü açıklanır. İkinci yarı başlıca güç reaktörü tiplerini—PWR, BWR, HTGR, LMFBR, GCFR—soğutucularını, moderatörlerini, yakıt çevrimlerini ve yakıt sürdürülebilirliğini belirleyen dönüşüm/üretim oranını inceler. Bölüm, nötroniği, ısıl-hidroliği, malzemeleri ve ekonomiyi birbirine bağlayarak reaktör tasarımının çok disiplinli kapsamıyla kapanır.
Çoğaltma Faktörü ve Nükleer Kritiklik
Kararlı bir zincirleme reaksiyon, her fisyondan tam olarak bir nötronun başka bir fisyon başlatmasını, geri kalanların yakalanma veya kaçakla kaybolmasını gerektirir. Bu denge çoğaltma faktörü k ile nicelenir; önce ardışık nesillerdeki nötron oranı, sonra daha pratik biçimde üretim/kayıp oranı k=P(t)/L(t) olarak tanımlanır. Nötron ömrü l=N(t)/L(t)'dir. Kritiklik sınıflaması (k<1 kritik-altı, k=1 kritik, k>1 kritik-üstü) nüfusun azalıp, sabit kalıp yoksa büyüyeceğini belirler. Basit kinetik model dN/dt=(k-1)N/l, reaktör periyodu T=l/(k-1) ile üstel davranış verir. Biçimsel 'nötronu izle' yaşam-döngüsü hesabı, sonsuz ortam için dört-faktör formülü k_sonsuz=eta*f*p*epsilon ve sonlu boyut kaçağı eklendiğinde altı-faktör formülü k=k_sonsuz*P_FNL*P_TNL getirir. Termal kullanım f, üreme faktörü eta, rezonanstan kaçış p ve hızlı fisyon epsilon her biri tesir kesitlerine ve çekirdek bileşimine bağlıyken, kaçmama faktörleri geometriye bağlıdır.
Nükleer Güç Reaktörlerine Giriş
Bir güç reaktörü, fisyon enerjisini bir termodinamik çevrimi süren kullanışlı ısıya dönüştürür. Bölüm baskın ticari reaktör sınıflarını inceler: hafif-su reaktörleri (LWR)—soğutucunun çekirdekte kaynamasına izin veren kaynar sulu reaktör (BWR) ve kütle kaynamasını bastırmak için ~155 bar'da tutulan basınçlı sulu reaktör (PWR)—her ikisi de sıradan suyu hem moderatör hem soğutucu olarak ve düşük zenginleştirilmiş UO2 yakıt kullanır. Yüksek sıcaklıklı gazla soğutulan reaktör (HTGR), yüksek termal verim ve ileri dönüştürücü (Th/U-233) çevrimi için grafit moderasyonu ve helyum soğutucu kullanır. Hızlı üretken reaktörler (sıvı sodyumla soğutulan LMFBR, helyumla soğutulan GCFR), moderatörsüz olarak U-238/Pu-239 çevriminde hızlı nötronlarla çalışır ve tükettiğinden fazla fisil yakıt üretmeyi amaçlar. Dönüşüm oranı CR (>1 olduğunda üretim oranı BR), fisil üretimi tüketime karşı niceler: LWR ~0.6, HTGR ~0.8, LMFBR >1. Üretim için eta>2 gerekir. Her tip; nötron spektrumu, yakıt kullanımı, güvenlik ve soğutucu kimyası arasında farklı ödünleşimleri barındırır.
Nükleer Reaktör Tasarımı
Reaktör tasarımı, yinelemeli ve çok disiplinli bir optimizasyon problemidir. Nötronik hedef—çalışma çevrimi boyunca k_eff=1'i sağlarken güç dağılımını şekillendirmek—ısıl-hidrolik sınırlara (en yüksek yakıt ve kılıf sıcaklıkları, çekirdekçik kaynamasından ayrılma), malzeme kısıtlarına (radyasyon hasarı, korozyon, mekanik bütünlük), yakıt-çevrimi ekonomisine ve güvenlik/lisanslama gereksinimlerine bağlıdır. Çekirdek boyutu ve geometrisi çoğunlukla, kabul edilebilir güç yoğunluğunda belirli bir gücü sağlama gereksinimi ile belirlenir; sonra bileşim, yakıt tüketimini ve fisyon-ürünü zehirlenmesini telafi edecek fazla reaktivite ile kritikliğe ulaşmak üzere ayarlanır. Kontrol, k'yi değiştirerek sağlanır—gücü artırmak için soğurucuları çekip kritik-üstüne geçmek, kritiğe dönmek için geri sokmak—kontrol çubukları, çözünmüş boron veya yanabilir zehirler kullanılarak. Yansıtıcılar kaçağı azaltır ve akıyı düzleştirir. Bölüm, altı-faktör formülü gibi ayrılabilir faktörlerin pedagojik olarak yararlı olduğunu, ancak gerçekçi tasarımın nötron taşınım/difüzyon denklemlerinin ısıl ve tükenme modelleriyle birlikte bağlaşık, yinelemeli sayısal çözümünü gerektirdiğini vurgular.
Temel Denklemler (16)
k'nin yaşam-döngüsü tanımı: iki ardışık fisyon neslindeki nötron nüfuslarının oranı. Kritikliği sınıflandırmak için kavramsal olarak kullanılır.
k'nin pratik, zamana bağlı tanımı: nötron üretim hızının toplam kayıp hızına (soğurma artı kaçak) oranı. Kinetik modellemenin temeli.
Ortalama nötron ömrünü, anlık nötron nüfusunun kayıp hızına oranı olarak tanımlar; zincirleme reaksiyonun zaman ölçeğini belirler.
Temel nötron nüfus dengesi: değişim hızı üretim eksi kayıba eşittir. Nokta kinetiğin başlangıç noktası.
k ve l'nin dengeye yerleştirilmesiyle elde edilen indirgenmiş nokta-kinetik ADD; gecikmeli nötronları ihmal eder. Ani üstel tepkiyi yönetir.
Sabit k ve l için basit kinetik ADD'nin kapalı-form çözümü. Azalma (k<1), kararlı durum (k=1) veya kontrolsüz büyüme (k>1) gösterir.
Üstel güç değişiminin zaman sabiti. k 1'e yaklaştıkça T sonsuza yaklaşır (kararlı durum). Gecikmeli nötronlar ihmal edildiğinde hızlı ani tepkiyi gösterir.
Soğurulan bir nötronun yakıtta soğurulma koşullu olasılığı: yakıtın makroskopik soğurma tesir kesitinin toplama oranı. Yakıt yüklemesiyle artar.
Yakıttaki soğurmanın ışınımsal yakalama yerine fisyon başlatma koşullu olasılığı: yakıt fisyon tesir kesitinin yakıt soğurma tesir kesitine oranı.
Yakıtta soğurulan başına üretilen ortalama fisyon nötronu sayısı. nu çarpı yakıttaki fisyon olasılığına eşittir; temel yakıt-kalitesi parametresi (üretim için eta>2 gerekir).
Hızlı/rezonans düzeltmelerinden önce, çıplak tek-enerjili model için sonsuz-ortam çoğaltma faktörü: üreme faktörü ile termal kullanımın çarpımı.
Enerji bağımlılığını hesaba katan termal reaktör için sonsuz-ortam çoğaltma faktörü: üreme, termal kullanım, rezonanstan kaçış ve hızlı fisyon.
Sonlu bir reaktör için etkin çoğaltma faktörü: dört faktör çarpı hızlı ve termal kaçmama olasılıkları. Kritiklik koşulu k=k_eff=1'dir.
Sonlu-reaktör etkin çoğaltmasını toplam kaçmama olasılığı aracılığıyla sonsuz-ortam değerine bağlar; malzeme özelliklerini geometriden ayırır.
Kritiklikten kesirsel sapma, doğal kontrol değişkeni. rho=0 kritik, rho>0 kritik-üstü, rho<0 kritik-altı; periyot reaktivite ile ters orantılıdır.
Üretken malzemenin yeni fisil yakıta ne kadar verimli dönüştürüldüğünü ölçer. CR>1 (üretim oranı BR), bir reaktörün tükettiğinden fazla yakıt ürettiği anlamına gelir; LWR~0.6, HTGR~0.8, LMFBR>1.
Deneyin: zincirleme reaksiyon
Her fisyon, yeni fisyonları tetikleyebilecek nötronlar salar. Çoğaltma faktörü popülasyonun büyüyeceğini, sabit mi kalacağını yoksa söneceğini mi belirler. Kontrol çubuklarını sokup çıkarmak için kaydırıcıyı sürükleyin ve nötron popülasyonunun tepkisini izleyin.
Diyagramı okuma
- Çubukları çekin (sola kaydırın): daha az nötron soğurulur, k 1'in üstüne çıkar, popülasyon büyür — kritik üstü.
- Çubukları sokun (sağa kaydırın): çubuklar nötronları yutar, k 1'in altına düşer, zincir söner — kritik altı.
- Gerçek santraller k'yı 1.000'e son derece yakın tutar ve kontrolü fiziksel olarak mümkün kılmak için gecikmiş nötronlara dayanır.
Simülasyon tasarımları
Öğrenciler k ve nötron ömrü l'yi ayarlar ve nötron nüfusu N(t)'nin canlı bir grafikte gerçek zamanlı evrimini izler; azalma (k<1), düz çizgi (k=1) veya patlayıcı üstel büyüme (k>1) görür. Bir gösterge reaktör periyodu T=l/(k-1) ve reaktivite rho'yu gösterir. Öğrenciler, yalnızca ani nötronlar modellendiğinde k'deki küçük %0.1'lik bir değişimin neden devasa, neredeyse anlık güç dalgalanmalarına yol açtığını keşfeder ve gecikmeli nötronlara güdülenir.
Öğrenciler altı-faktör formülünün her faktörünü (eta, f, p, epsilon, P_FNL, P_TNL) ayarlar ve k_inf ile k_eff'in anında güncellendiğini görür; bir kritiklik göstergesi kritik-altı/kritik/kritik-üstüne oturur. f ve p'yi bir yakıt-moderatör oranına bağlamak klasik ödünleşimi gösterir: daha fazla yakıt f'yi yükseltir ama p'yi düşürür ve k'yi maksimize eden optimal bir moderasyon oranı üretir.
Öğrenciler bir reaktör sınıfı (PWR, BWR, HTGR, LMFBR, GCFR) seçer ve moderatörünü, soğutucusunu, nötron spektrumunu, yakıtını ve dönüşüm oranını inceler. Bir spektrum grafiği termalden (0.05 eV yakınında tepe) hızlıya (100 keV+ aralığında tepe) dönüşür. Bir yakıt-envanteri çubuğu, bir çalışma çevrimi boyunca fisil üretimi ile tüketimini canlandırarak üretimi (CR>1) ile yakmayı (CR<1) gösterir.
Sayısal veriler
Bu bölüm, öğrencilerin reaktör fiziğini ilk kez çalışan koda dönüştürdüğü yerdir. Temel hesaplama görevleri şunlardır: (1) nokta-kinetik ADD dN/dt=(k-1)N/l'yi integre etmek, analitik üstel çözümü sayısal RK4 ile karşılaştırmak ve gecikmeli-nötron öncül denklemleri eklendiğinde katılığı (stiffness) incelemek; (2) dört- ve altı-faktör formüllerini değerlendirmek, termal kullanım f ile rezonanstan kaçış p'yi makroskopik tesir kesitlerinden ve malzeme parçacık yoğunluklarından hesaplamak, ardından kritik yakıt/moderatör oranı için 1-B kök bulma (ikiye bölme veya Newton) yapmak; ve (3) termal ile hızlı sistemler için temsili nötron enerji spektrumları üretmek. Bunlar gerçek kodlarda kullanılan üretim düzeyindeki yöntemlerin habercisidir: deterministik nötron-difüzyon ve taşınım çözücüler (CITATION, PARCS ve DIF3D'nin kavramsal ataları), az-grup tesir kesitlerini homojenleştiren spektral/kafes kodları (CASMO/SCALE tarzı) ve güç iterasyonuyla referans k_eff için Monte Carlo taşınımı (MCNP, OpenMC, Serpent). Öğrenciler k_eff'in taşınım operatörünün baskın özdeğeri olduğunu ve gerçek kritiklik aramasının ayrılabilir faktörlerin çarpımı değil, bir özdeğer/iterasyon problemi olduğunu öğrenir. Bu platformda hızlı iç döngüler (ADD integrasyonu, özdeğer/güç iterasyonu, tesir-kesiti birleştirme) tarayıcıda neredeyse yerel hızda çalışması için WebAssembly'ye derlenen Rust ile uygulanırken, TypeScript tepkisel arayüzü, çizimi ve öğrenci kod-düzenleme korumalı alanını yönetir.