Ders 10 / 15İleri~35 dk kaldıDuderstadt & Hamilton (1976)

Heterojen Kor Kafesleri için Hücre Hesaplamaları

Gerçek güç reaktörleri yakıt, yavaşlatıcı, soğutucu ve yapı malzemesinin homojen karışımları değildir; yavaşlatıcı içine daldırılmış ayrık yakıt elemanlarının (çubuklar, plakalar veya demetler) düzenli bir kafesinden oluşur. Bu heterojenlik, homojenleştirilmiş bir kora kıyasla nötron fiziğini derinden değiştirir: termal akı, güçlü soğuran yakıt içinde bastırılır ve yavaşlatıcıda tepe yapar; epitermal rezonans soğurması ise güçlü biçimde öz-perdelenir, böylece homojen bir tahminin öngöreceğinden daha az nötron U-238 yakalamasında kaybedilir. Bölüm 10, tekrarlayan tek bir birim hücre üzerinde ayrıntılı bir transport hesabı yaparak ve sonuçları kaba örgülü tam-kor difüzyon modelinde kullanılabilecek uzaysal ortalamalı, öz-perdelenmiş az-gruplu tesir kesitlerine indirgeyerek bu etkileri çok-gruplu difüzyon çerçevesinde yakalayacak araçları geliştirir. Bölüm, termal heterojen etkileri termal dezavantaj faktörü ve termal kullanım üzerinden inceler; yavaşlatıcı ile yakıt arasındaki akı oranını hesaplamak için çarpışma-olasılığı yöntemlerini (Wigner rasyonel yaklaşımı, ABH yöntemi ve Amouyal–Benoist–Horowitz işlemi) tanıtır. Ardından hızlı heterojen etkileri ele alır: eşdeğerlik teoremi ve Dancoff faktörü aracılığıyla rezonans öz-perdelemesi, yakıtın yumrulaştırılmasıyla rezonans kaçış olasılığının artması ve hızlı fisyon faktörü. Dört-faktör formülü, kafes geometrisinin termal kullanımı rezonans kaçışına karşı nasıl takas ettiğini göstermek için faktör faktör yeniden incelenir ve böylece optimum bir yakıt-yavaşlatıcı oranının varlığı gerekçelendirilir. Bölüm, çok-gruplu difüzyonun nerede başarılı olduğunu ve transport düzeltmelerinin nerede zorunlu olduğunu dürüstçe değerlendiren kapanış notlarıyla biter ve modern kafes-fiziği kodlarına zemin hazırlar. Hücre hesaplamalarına hâkim olmak, mikroskobik tesir-kesiti verisi ile mühendislikte kullanılabilir tam-kor model arasındaki temel köprüdür.

Nükleer Reaktör Analizinde Kafes Etkileri

Bu kısım, heterojen bir kafesin neden basitçe homojen bir karışıma yayılarak ele alınamayacağını gerekçelendirir. Yakıt, yavaşlatıcı, soğutucu ve kılıf boyutları bir nötron ortalama serbest yoluyla karşılaştırılabilir olan farklı uzaysal bölgeler kapladığından, akı tek bir birim hücre içinde güçlü biçimde konuma bağlıdır. Strateji iki adımlı bir homojenleştirmedir: önce tekrarlayan tek bir hücre üzerinde (kafes periyodikliğini zorlayan yansıtıcı/beyaz sınır koşullarıyla) ince ölçekli bir sabit-kaynak veya özdeğer transport problemi çözülür, ayrıntılı spektrum ve uzaysal akı şekli elde edilir; ardından eşdeğer homojenleştirilmiş, öz-perdelenmiş grup sabitlerini üretmek için tesir kesitleri akı-hacim ortalaması alınır. Bu hücre-ortalamalı sabitler doğru reaksiyon hızlarını korur ve kaba örgülü tam-kor çok-gruplu difüzyon hesabına yerleştirilir. Kısım, birim hücreyi, kare veya altıgen hücreyi eşit alanlı halkasal bir hücreyle değiştiren Wigner–Seitz silindirleştirmesini tanımlar ve hücre içinde difüzyon teorisinin neden tek başına yetersiz kaldığını, küresel olarak difüzyon yeterli olsa bile yerel olarak çarpışma-olasılığı veya transport yöntemleri gerektiğini açıklar.

unit cellhomogenizationWigner–Seitz cellcell-averaged cross sectionslattice periodicityreflecting boundary conditionflux-volume weighting

Termal Nötron Fiziğinde Heterojen Etkiler

Termal bölgede yakıt güçlü bir soğurucudur, bu nedenle termal akı yakıt yumrusunda bastırılır ve nötronların yavaşlama yoluyla doğduğu yavaşlatıcıda yükselir. Bu uzaysal çukur, daima birden büyük olan ve yavaşlatıcı-ortalamalı termal akının yakıt-ortalamalı termal akıya oranı olarak tanımlanan termal dezavantaj faktörüyle nicelenir. Soğurma hızı akı ile makroskobik tesir kesitinin çarpımı olduğundan, bastırılmış yakıt akısı yakıtta soğurulan termal nötron oranını azaltır ve dolayısıyla termal kullanım f'yi homojen değerine göre düşürür. Kısım, akı-dezavantaj ağırlıklandırması kullanarak heterojen bir hücre için f'yi türetir ve akı oranını hesaplamak için pratik çarpışma-olasılığı reçeteleri sunar: kaçış olasılıkları için Wigner rasyonel yaklaşımı, silindirik yakıt çubukları için doğruluğu artıran ABH (Amouyal–Benoist–Horowitz) yöntemi ve silindirik-hücre transportunun kullanımı. Bu rekabet kafes tasarımının merkezindedir: yakıtın yumrulaştırılması termal kullanımı düşürür (k için kötü), ancak bir sonraki kısımda gösterildiği gibi aynı anda rezonans kaçışını iyileştirir (k için iyi) ve net etki optimum kafesi belirler.

thermal disadvantage factorthermal utilizationflux depressioncollision probabilityWigner rational approximationABH methodAmouyal–Benoist–Horowitzcylindrical cell

Hızlı Nötron Fiziğinde Heterojen Etkiler

Baskın hızlı/epitermal heterojen etki, U-238 yakalamasının rezonans öz-perdelenmesidir. Homojen bir karışımda her yakıt atomu tam epitermal akıyı görür, bu yüzden rezonans soğurması büyüktür. Yakıt yumrulaştırıldığında, yumru yüzeyindeki güçlü rezonans soğurması çubuğun derinlerinde rezonans enerjilerindeki akıyı bastırır (enerji öz-perdelemesi) ve yumrunun geometrik boyutu, giren nötronlar için mevcut yüzeyi sınırlar (uzaysal öz-perdeleme). Her iki etki de etkin rezonans integralini keskin biçimde azaltır ve dolayısıyla rezonans kaçış olasılığı p'yi yükseltir. Kısım, etkin rezonans integralini yüzey-kütle oranı S/M'nin bir fonksiyonu olarak geliştirir; heterojen bir yumrunun eşdeğer homojen bir sisteme eşlenmesini sağlayan eşdeğerlik teoremini ve nötronların bir yakıt çubuğundan komşusuna yavaşlatılmadan akışını hesaba katan Dancoff (Dancoff–Ginsburg) düzeltmesini ele alır. Ayrıca, bir kafeste artan hızlı fisyon faktörü epsilon'u, yani U-238'in eşiğinin üzerindeki hızlı fisyonundan kaynaklanan küçük ama ihmal edilemez çarpanı işler; çünkü bir çubukta doğan hızlı nötronlar yakıttan çıkmadan hızlı fisyona yol açabilir.

resonance self-shieldingeffective resonance integralequivalence theoremDancoff factorsurface-to-mass rationarrow/wide resonancefast fission factorresonance escape probability

Çok-Gruplu Difüzyon Yöntemi Üzerine Bazı Kapanış Notları

Son kısım, artık hücre-homojenleştirmesi ve öz-perdelemeyle güçlendirilmiş çok-gruplu difüzyon yöntemini bir bütün olarak değerlendirmek için geri çekilir. Tam-kor difüzyon hesabının doğruluğunun, difüzyon çözücüsünün kendisine değil, hücre hesabının ürettiği az-gruplu sabitlerin kalitesine kritik biçimde bağlı olduğunu vurgular. Kısım, yaklaşım hiyerarşisini—sürekli-enerji transporttan çok-gruplu transporta, oradan uzaysal homojenleştirilmiş çok-gruplu difüzyona—gözden geçirir ve her birinin nerede gerekçeli olduğunu belirler: difüzyon, güçlü soğurucuların, vakum sınırlarının, malzeme arayüzlerinin yakınında ve yüksek anizotropik veya boşluklu bölgelerde başarısız olur; buralarda transport (Sn, çarpışma olasılığı veya Monte Carlo) gereklidir. Reaktör analizinin pratik iş akışını tartışır: tesir-kesiti kütüphaneleri, homojenleştirilmiş sabitler üreten hücre/kafes kodları ve tam kor için küresel difüzyon veya nodal kodlar; ve homojenleştirilmiş difüzyonu altta yatan transport çözümüyle uzlaştıran eşdeğerlik teorisine ve süreksizlik faktörlerine işaret eder. Bölüm, bu hücre yöntemlerini sonraki tüm kor tasarımı, tüketim ve kinetik analizlerinin üzerine inşa edildiği temel olarak çerçeveleyerek kapanır.

homogenization accuracyfew-group constantsvalidity of diffusion theorytransport correctionsequivalence theorydiscontinuity factorsreactor analysis workflow

Temel ilişki

k<1 → kritik altı (reaksiyon söner) · k=1 → kritik (kararlı) · k>1 → kritik üstü (güç artar). Reaktörler tam k=1'de çalışır ve ince ayar için kontrol çubuklarını kullanır.

Temel Denklemler (14)

Termal Dezavantaj Faktörü

(10.1)

\zeta = \frac{\bar{\phi}_M}{\bar{\phi}_F} = \frac{\displaystyle\frac{1}{V_M}\int_{V_M}\phi(\mathbf{r})\,dV}{\displaystyle\frac{1}{V_F}\int_{V_F}\phi(\mathbf{r})\,dV}

Bir birim hücre içinde yavaşlatıcıdaki hacim-ortalamalı termal akının yakıttakine oranı. Soğuran yakıt akıyı bastırdığından birden büyüktür; akı çukurunu doğrudan niceler ve heterojen termal kullanıma girdi sağlar.

sembolanlambirim

\zeta

termal dezavantaj faktörü–

\bar{\phi}_M

yavaşlatıcı-ortalamalı termal akın m^-2 s^-1

\bar{\phi}_F

yakıt-ortalamalı termal akın m^-2 s^-1

V_M

hücredeki yavaşlatıcı hacmim^3

V_F

hücredeki yakıt hacmim^3

Heterojen Termal Kullanım

(10.2)

f = \frac{\Sigma_a^F V_F \bar{\phi}_F}{\Sigma_a^F V_F \bar{\phi}_F + \Sigma_a^M V_M \bar{\phi}_M} = \frac{\Sigma_a^F V_F}{\Sigma_a^F V_F + \Sigma_a^M V_M \,\zeta}

Heterojen bir hücre için yakıtta soğurulan termal nötron oranı. Dezavantaj faktörü zeta yavaşlatıcı soğurma terimini çarpar; yakıt akısı bastırıldığı için f'yi homojen duruma göre düşürür.

sembolanlambirim

f

termal kullanım–

\Sigma_a^F

yakıtın makroskobik soğurma tesir kesitim^-1

\Sigma_a^M

yavaşlatıcının makroskobik soğurma tesir kesitim^-1

\zeta

termal dezavantaj faktörü–

Dört-Faktör Formülü

(10.3)

k_\infty = \eta\, f\, p\, \varepsilon

Sonsuz-ortam çarpan faktörünün üreme faktörü, termal kullanım, rezonans kaçış olasılığı ve hızlı fisyon faktörüne ayrıştırılması. Kafes heterojenliği f, p ve epsilon'u değiştirir; formül bu bölümün merkezindeki tasarım takaslarını açığa çıkarır.

sembolanlambirim

k_\infty

sonsuz çarpan faktörü–

\eta

üreme faktörü (yakıtta termal soğurma başına nötron)–

f

termal kullanım–

p

rezonans kaçış olasılığı–

\varepsilon

hızlı fisyon faktörü–

Rezonans Kaçış Olasılığı

(10.4)

p = \exp\!\left[-\frac{N_F V_F\, I_{\mathrm{eff}}}{\xi\,\Sigma_s\,V_M + \xi\,\Sigma_s\,V_F}\right]

Bir nötronun yakalanmadan U-238 rezonans bölgesini geçerek yavaşlama olasılığı. Etkin rezonans integrali I_eff, sonsuz-seyreltme integralinin yerini alır ve bir yumrudaki öz-perdelemeyle güçlü biçimde azalır; bu da heterojen bir kafeste p'yi yükseltir.

sembolanlambirim

p

rezonans kaçış olasılığı–

N_F

yakıttaki rezonans soğurucusunun sayı yoğunluğum^-3

I_{\mathrm{eff}}

etkin rezonans integralibarn

\xi

çarpışma başına ortalama logaritmik enerji kaybı–

\Sigma_s

makroskobik saçılma tesir kesitim^-1

Etkin Rezonans İntegrali (Tanım)

(10.5)

I_{\mathrm{eff}} = \int \sigma_a(E)\,\frac{\phi(E)}{\phi_\infty(E)}\,\frac{dE}{E}

Gerçek öz-perdelenmiş akının asimptotik 1/E akısına oranıyla ağırlıklandırılmış rezonans integrali. Öz-perdeleme phi(E)/phi_inf(E)'nin rezonans tepelerinde çukurlaşmasına yol açar, bu yüzden I_eff sonsuz-seyreltme integralinden çok daha küçüktür.

sembolanlambirim

\sigma_a(E)

enerjiye karşı mikroskobik soğurma tesir kesitibarn

\phi(E)

gerçek (öz-perdelenmiş) birim enerji başına akın m^-2 s^-1 eV^-1

\phi_\infty(E)

asimptotik 1/E akısın m^-2 s^-1 eV^-1

E

nötron enerjisieV

S/M'ye Karşı Etkin Rezonans İntegrali (Eşdeğerlik)

(10.6)

I_{\mathrm{eff}} = a + b\,\sqrt{\frac{S}{M}}

Ampirik eşdeğerlik bağıntısı: bir yumrunun etkin rezonans integrali, yüzey-kütle oranının kareköküyle doğrusal olarak artar. Daha büyük S/M (daha ince çubuklar) daha az öz-perdeleme ve daha fazla rezonans soğurması demektir; bu da p'yi düşürür.

sembolanlambirim

I_{\mathrm{eff}}

etkin rezonans integralibarn

S

yakıt yumrusunun yüzey alanım^2

M

yakıt yumrusunun kütlesikg

a,b

soğurucu/geometri için ampirik sabitlerbarn, barn (m^2/kg)^-1/2

Dancoff Düzeltme Faktörü

(10.7)

C = \frac{1}{S}\oint_S P_{\mathrm{esc}\to\mathrm{fuel}}\, dS \;,\qquad \Sigma_e' = \Sigma_e\,(1-C)

C, bir yakıt çubuğunun yüzeyinden ayrılan bir nötronun yavaşlatıcıda çarpışma yapmadan başka bir yakıt çubuğuna girme olasılığıdır. Etkin kaçış (yüzey) tesir kesitini (1-C) ile azaltır ve sıkı kafeslerde öz-perdelemeyi artırır; C=0 yalıtılmış çubuk, C->1 yoğun paketleme.

sembolanlambirim

C

Dancoff faktörü–

P_{\mathrm{esc}\to\mathrm{fuel}}

çarpışmadan başka bir çubuğa ulaşma olasılığı–

\Sigma_e

kaçış (yüzey) tesir kesiti, yalıtılmış çubukm^-1

\Sigma_e'

Dancoff-düzeltilmiş kaçış tesir kesitim^-1

Wigner Rasyonel Yaklaşımı (Kaçış Olasılığı)

(10.8)

P_{\mathrm{esc}}(E) = \frac{\Sigma_e}{\Sigma_t(E) + \Sigma_e}\;,\qquad \Sigma_e = \frac{1}{\bar{\ell}} = \frac{S}{4 V_F}

Dışbükey bir yumru için kapalı-form kaçış olasılığı: bir nötron, kaçış tesir kesitinin toplam tesir kesitine oranına eşit olasılıkla kaçar. Ortalama kiriş uzunluğu l-bar = 4V/S eşdeğer yüzey tesir kesitini verir; bu, eşdeğerlik teorisi ve rasyonel-yaklaşım rezonans işleminin köşe taşıdır.

sembolanlambirim

P_{\mathrm{esc}}

yumrudan ilk-uçuş kaçış olasılığı–

\Sigma_t(E)

yakıtın toplam makroskobik tesir kesitim^-1

\Sigma_e

kaçış (yüzey) tesir kesitim^-1

\bar{\ell}

yumrunun ortalama kiriş uzunluğu (4V/S)m

Çarpışma Olasılığı Karşılıklılığı

(10.9)

V_i\,\Sigma_{t,i}\,P_{i\to j} = V_j\,\Sigma_{t,j}\,P_{j\to i}

Bir hücrenin i ve j bölgeleri arasındaki ilk-uçuş çarpışma olasılıkları için karşılıklılık bağıntısı. Hücre-içi nötron dengesini çözmek ve akı oranını üretmek için kullanılan çarpışma-olasılığı (P-matrisi) yönteminin cebirsel omurgasıdır.

sembolanlambirim

P_{i\to j}

i'de doğan bir nötronun ilk çarpışmayı j'de yapma olasılığı–

V_i

i bölgesinin hacmim^3

\Sigma_{t,i}

i bölgesindeki toplam tesir kesitim^-1

Hızlı Fisyon Faktörü

(10.10)

\varepsilon = \frac{\text{total fissions (fast + thermal)}}{\text{thermal fissions}} \approx 1 + \frac{P_{\mathrm{ff}}\,\nu_{28}\,\sigma_{f,28}}{\dots}

Tüm fisyonlardan üretilen toplam nötronların yalnızca termal fisyonlardan üretilenlere oranı; U-238'in ~1 MeV eşiğinin üzerindeki hızlı fisyonunu hesaba katar. Birden biraz büyüktür (HSR'lerde yaklaşık 1.03) ve hızlı nötronlar yakıt içinde daha uzun kaldığından yumrulaştırmayla artar.

sembolanlambirim

\varepsilon

hızlı fisyon faktörü–

P_{\mathrm{ff}}

bir hızlı nötronun U-238 hızlı fisyonuna yol açma olasılığı–

\nu_{28}

U-238 hızlı fisyonu başına nötron–

\sigma_{f,28}

U-238 hızlı fisyon tesir kesitibarn

Wigner–Seitz Eşdeğer Hücre Yarıçapı

(10.11)

\pi R_{\mathrm{cell}}^2 = A_{\mathrm{cell}} \quad\Rightarrow\quad R_{\mathrm{cell}} = \sqrt{\frac{A_{\mathrm{cell}}}{\pi}} = p_{\mathrm{pitch}}\sqrt{\frac{1}{\pi}}

A_cell alanlı kare (veya altıgen) bir kafes hücresini, dış yarıçapı R_cell olan eşit alanlı halkasal bir silindirle değiştiren silindirleştirme; yavaşlatıcı-yakıt hacim oranını korurken 1-B radyal transport/çarpışma-olasılığı çözümünü mümkün kılar.

sembolanlambirim

R_{\mathrm{cell}}

eşdeğer silindirik hücre dış yarıçapım

A_{\mathrm{cell}}

gerçek kafes hücresinin alanım^2

p_{\mathrm{pitch}}

kafes adımı (çubuk merkez aralığı)m

Akı-Hacim Homojenleştirilmiş Grup Sabiti

(10.12)

\bar{\Sigma}_{x,g} = \frac{\displaystyle\sum_i V_i \,\Sigma_{x,g}^{(i)}\,\phi_{g}^{(i)}}{\displaystyle\sum_i V_i\,\phi_{g}^{(i)}}

x reaksiyonu ve g grubu için, hücre bölgeleri i üzerinde akı-hacim ağırlıklandırmasıyla elde edilen hücre-homojenleştirilmiş makroskobik tesir kesiti. Bu, heterojen hücrenin toplam reaksiyon hızını korur ve tam-kor difüzyon modelinde kullanılan sabitleri üretir.

sembolanlambirim

\bar{\Sigma}_{x,g}

homojenleştirilmiş tesir kesiti, reaksiyon x, grup gm^-1

\Sigma_{x,g}^{(i)}

i bölgesindeki tesir kesitim^-1

\phi_g^{(i)}

i bölgesindeki grup akısın m^-2 s^-1

V_i

i bölgesinin hacmim^3

Dar-Rezonans Öz-Perdelenmiş Akısı

(10.13)

\phi(E) = \frac{\Sigma_p + \Sigma_e}{\Sigma_t(E) + \Sigma_e}\,\frac{1}{E}

Heterojen bir yakıt yumrusu içindeki akı için dar-rezonans (NR) yaklaşımı; kaçış tesir kesiti Sigma_e yavaşlatıcıya sızıntıyı temsil eder. Toplam tesir kesiti bir rezonansta sıçradığında akı çukurlaşır; bu, enerji öz-perdelemesinin özüdür.

sembolanlambirim

\Sigma_p

potansiyel (asimptotik) saçılma tesir kesitim^-1

\Sigma_e

kaçış tesir kesiti S/4Vm^-1

\Sigma_t(E)

E enerjisinde toplam tesir kesitim^-1

E

nötron enerjisieV

Termal Kullanımın ABH Biçimi

(10.14)

\frac{1}{f} = 1 + \frac{\Sigma_a^M V_M}{\Sigma_a^F V_F} + \frac{\Sigma_a^M V_M}{V_F}\,(E-1)\,\frac{1}{\Sigma_a^F V_F / V_F}

Amouyal–Benoist–Horowitz ayrıştırması; 1/f'yi homojen bir terim ile (E-1) ile orantılı heterojen bir düzeltme olarak ifade eder; burada E, yakıt ve yavaşlatıcı bölgelerindeki çarpışma olasılıklarından elde edilen fazla soğurma (akı-bastırma) faktörüdür. Silindirik çubuklar için elle ve kodla yapılan kafes hesaplarında kullanılır.

sembolanlambirim

f

termal kullanım–

E

ABH fazla-soğurma (akı-bastırma) faktörü–

\Sigma_a^F

yakıt soğurma tesir kesitim^-1

\Sigma_a^M

yavaşlatıcı soğurma tesir kesitim^-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.

Etkileşimli · Nötron çoğalmasıCanlı simülasyon

Kontrol çubuklarık = 1.00Kritik

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ı

SİMÜLASYON TASARIMIYÜKSEK SADAKAT

Çubuk-Hücre Termal Akısı ve Dezavantaj Faktörü

Öğrenci yakıt ve yavaşlatıcı özellikleri ile çubuk geometrisini ayarlar, ardından bir Wigner–Seitz silindirik hücresi boyunca 1-B radyal termal akı profilini izler. Akı soğuran yakıt içinde çukurlaşır ve yavaşlatıcıda yükselir; uygulama dezavantaj faktörü zeta'yı ve ortaya çıkan termal kullanım f'yi raporlar ve yakıt-yavaşlatıcı oranı taranırken f'nin nasıl değiştiğini üst üste bindirir.

GİRDİLER

Fuel rod radius a (0.2–0.8 cm) · Cell outer radius R_cell (0.5–1.2 cm) · Fuel absorption Sigma_a^F (5–60 m^-1) · Moderator absorption Sigma_a^M (0.1–5 m^-1) · Moderator scattering Sigma_s^M (50–400 m^-1) · Diffusion coefficients D_F, D_M (auto from Sigma_tr)

ÇIKTILAR

1-D radial flux phi(r) curve (fuel + moderator) · Disadvantage factor zeta readout · Thermal utilization f readout · f vs V_F/V_M sweep plot

YÖNTEMTwo-region 1-D cylindrical diffusion (modified Bessel functions I0/K0) with interface continuity of flux and current, OR a 2-region collision-probability balance; solve the linear system for region-averaged fluxes and integrate for zeta.

KOD KANCASIÖğrenciler kaçış/çarpışma-olasılığı çekirdeğini (veya analitik Bessel katsayılarını) kendi P_escape(geometri, Sigma_t) fonksiyonlarıyla değiştirir; ekrandaki akı profili ve zeta canlı güncellenir, böylece Wigner rasyonel ile tam silindirik sonuçları karşılaştırabilirler.

SİMÜLASYON TASARIMIYÜKSEK SADAKAT

Rezonans Öz-Perdelemesi ve Etkin Rezonans İntegrali

Öğrenci tek bir U-238 rezonansını (Breit–Wigner) inceler ve yakıt yumru boyutunu / sayı yoğunluğunu artırır. Öz-perdelenmiş akı phi(E) rezonans tepesinde keskin bir çukur oluşturur ve hesaplanan etkin rezonans integrali sonsuz-seyreltme değerinden çöker; bu da bir göstergede gösterilen rezonans kaçış olasılığı p'yi yükseltir.

GİRDİLER

Resonance energy E0 (6.67 eV default) · Peak cross section / Gamma widths (Gamma_n, Gamma_gamma) · Fuel number density N_F (10^27–10^28 m^-3) · Surface-to-mass ratio S/M (rod radius slider) · Dancoff factor C (0–0.6) · Moderator scattering Sigma_p (potential)

ÇIKTILAR

Breit–Wigner sigma(E) curve · Self-shielded flux phi(E) with resonance dip · Effective resonance integral I_eff readout · Resonance escape probability p gauge · I_eff vs sqrt(S/M) line

YÖNTEMEvaluate Breit–Wigner cross section on a fine energy mesh; apply NR (narrow-resonance) flux model phi(E)=(Sigma_p+Sigma_e)/(Sigma_t+Sigma_e)/E with Dancoff-corrected Sigma_e=(1-C)S/4V; numerically integrate I_eff = ∫ sigma_a (phi/phi_inf) dE/E (trapezoid/Simpson).

KOD KANCASIÖğrenciler kendi akı modellerini (NR, WR veya ara-rezonans lambda) ve I_eff integratörünü uygular; çukur şekli, I_eff sayısı ve p göstergesi anında tepki verir ve yumrulaştırmanın p'yi neden yükselttiğini gösterir.

SİMÜLASYON TASARIMIORTA SADAKAT

Dört-Faktör Kafes Optimize Edici

Öğrenci bir birim hücrenin yakıt-yavaşlatıcı hacim oranını değiştirir ve f, p, epsilon ve k_inf'in canlı hesaplandığını görür. Yavaşlatıcı oranı arttıkça f yükselir ama p düşer; k_inf = eta*f*p*epsilon çarpımı, açık bir optimum (tepe) yavaşlatma noktasına sahip bir eğri çizer ve heterojen kafes tasarım takasını doğrudan gösterir.

GİRDİLER

Moderator-to-fuel volume ratio V_M/V_F (1–20) · Rod radius (couples zeta & S/M) · Fuel enrichment (sets eta) · Use heterogeneous vs homogeneous toggle · Dancoff factor / lattice tightness

ÇIKTILAR

Live readouts of eta, f, p, epsilon · k_inf vs V_M/V_F curve with optimum marker · Bar comparison: homogeneous vs heterogeneous k_inf · Annotated under/over-moderation regions

YÖNTEMCouple the disadvantage-factor model (for f) and the equivalence/Dancoff resonance-integral model (for p) and the empirical epsilon to assemble k_inf across a parameter sweep; locate the maximum by sampling or golden-section search.

KOD KANCASIÖğrenciler f(V_M/V_F) ve p(S/M, C) fonksiyonlarını sağlar; optimize edici k_inf eğrisini yeniden hesaplar ve optimum işaretçisi kayar; böylece fizik seçimlerinin tepe-yavaşlatma tasarım noktasını nasıl hareket ettirdiğini gösterir.

Sayısal veriler

Typical thermal disadvantage factor (LWR pin cell)zeta ≈ 1.05 – 1.20Hafif-su kafesinde yavaşlatıcı termal akısı yakıt akısını %5–20 aşar.

Fast fission factor (LWR)epsilon ≈ 1.02 – 1.04Küçük U-238 hızlı-fisyon katkısı; daha sıkı, daha metalik kafeslerde daha büyük.

Resonance escape probability (LWR)p ≈ 0.70 – 0.90Yakıt yumrulaştırmasıyla homojen karışıma göre güçlü biçimde artar.

U-238 infinite-dilution resonance integralRI_inf ≈ 280 barnÖz-perdelenmiş bir HSR yakıt çubuğunda etkin değer ~9–25 barn'a düşer.

U-238 fast fission threshold≈ 1.0 – 1.4 MeVÜzerinde U-238 fisyonunun anlamlı hale geldiği enerji; epsilon'un kaynağı.

Mean chord length of convex lumpl-bar = 4V/SCauchy bağıntısı; eşdeğer kaçış tesir kesiti Sigma_e = S/4V'yi belirler.

Typical LWR lattice pitch / rod diameterpitch ≈ 1.26 cm, d ≈ 0.95 cmHücre hesaplarında kullanılan tipik BSR yakıt çubuğu geometrisi.

Dancoff factor rangeC = 0 (isolated) → ~0.3–0.5 (tight lattice)Daha yüksek C, çubuktan-çubuğa daha fazla gölgeleme ve daha güçlü öz-perdeleme demektir.

YAZILIM NOTU

Bu bölüm, kafes fiziğinin hesaplama çekirdeğini öğretir: mikroskobik tesir kesitlerini öz-perdelenmiş, homojenleştirilmiş az-gruplu sabitlere nasıl dönüştüreceğinizi. Öğrenciler şunları uygular: (1) hücre-içi akı oranlarını elde etmek için ilk-uçuş kaçış olasılıkları ve karşılıklılık bağıntısını kullanan çarpışma-olasılığı (Pij) çözücüleri; (2) termal dezavantaj faktörü için Wigner rasyonel ve ABH yaklaşımları ile tam iki-bölgeli silindirik difüzyon çözümü (değiştirilmiş Bessel fonksiyonları); (3) dar-/geniş-rezonans akı modelleri, eşdeğerlik teoremi (Sigma_e = S/4V) ve Dancoff düzeltmesi aracılığıyla rezonans öz-perdelemesi; Breit–Wigner tesir kesitlerini integralleyerek etkin rezonans integralleri elde etme; ve (4) grup sabitleri üretmek için akı-hacim homojenleştirmesi. Bunlar tam olarak üretim kafes kodlarının—WIMS, CASMO, APOLLO, HELIOS, SCALE/NEWT, Dragon ve MPACT/OpenMOC'taki karakteristik-yöntemi çözücülerin—içindeki algoritmalardır; bunların çıktısı tam-kor difüzyon/nodal kodları (PARCS, SIMULATE, NESTLE) besler. Platformda, sayısal olarak yoğun çekirdekler (rezonans-integral kuadratürü, Bessel değerlendirmesi, Pij matris kurulumu, altın-oran optimizasyonu) neredeyse yerel hız için WebAssembly'ye derlenen Rust ile yazılır; etkileşimli kontroller, grafikler ve kod-kanca editörü ise TypeScript/React'tır. Öğrenciler küçük sayısal rutinler (enerji-örgü integrasyonu, kök bulma, doğrusal-sistem çözümleri) yazmayı öğrenir ve etkilerini akı profilleri, rezonans integralleri ve k_inf üzerinde anında görür; böylece ders kitabı teorisi ile gerçek reaktör-analizi yazılım mühendisliği arasında köprü kurarlar.

ALGORİTMALAR

Collision probability (Pij) methodFirst-flight escape probabilityWigner rational approximationABH (Amouyal–Benoist–Horowitz) methodNarrow/wide resonance self-shieldingEquivalence theorem & Dancoff correctionBreit–Wigner resonance integral quadratureTwo-region cylindrical diffusion (Bessel functions)Flux-volume homogenizationGolden-section optimization

DİLLER

Rust/WASMTypeScript

İlgili simülasyon

Öğrenciler basamak, rampa, sinüzoidal veya darbe reaktivitesi ekler ve reaktör gücü ile altı öncül konsantrasyonunun gerçek zamanlı evrimini izler. Prompt sıçramayı, kararlı asimptotik periyodu ve ρ bir dolara (hızlı kritik) yaklaştıkça periyodun dramatik çöküşünü gözlemlerler. Canlı bir gösterge inhour denkleminden periyodu gösterir ve hızlı-kritik koşullarını işaretler.

İlgili simülasyonu aç

← Önceki

Termal Spektrum Hesaplamaları ve Termal Grup Sabitleri

Nükleer Reaktör Kalbi Tasarımının Genel Yönleri

Ders 10 / 15İleri~35 dk kaldıDuderstadt & Hamilton (1976)

Heterojen Kor Kafesleri için Hücre Hesaplamaları

Nükleer Reaktör Analizinde Kafes Etkileri

unit cellhomogenizationWigner–Seitz cellcell-averaged cross sectionslattice periodicityreflecting boundary conditionflux-volume weighting

Termal Nötron Fiziğinde Heterojen Etkiler

thermal disadvantage factorthermal utilizationflux depressioncollision probabilityWigner rational approximationABH methodAmouyal–Benoist–Horowitzcylindrical cell

Hızlı Nötron Fiziğinde Heterojen Etkiler

resonance self-shieldingeffective resonance integralequivalence theoremDancoff factorsurface-to-mass rationarrow/wide resonancefast fission factorresonance escape probability

Çok-Gruplu Difüzyon Yöntemi Üzerine Bazı Kapanış Notları

homogenization accuracyfew-group constantsvalidity of diffusion theorytransport correctionsequivalence theorydiscontinuity factorsreactor analysis workflow

Temel ilişki

k<1 → kritik altı (reaksiyon söner) · k=1 → kritik (kararlı) · k>1 → kritik üstü (güç artar). Reaktörler tam k=1'de çalışır ve ince ayar için kontrol çubuklarını kullanır.

Temel Denklemler (14)

Termal Dezavantaj Faktörü

(10.1)

\zeta = \frac{\bar{\phi}_M}{\bar{\phi}_F} = \frac{\displaystyle\frac{1}{V_M}\int_{V_M}\phi(\mathbf{r})\,dV}{\displaystyle\frac{1}{V_F}\int_{V_F}\phi(\mathbf{r})\,dV}

sembolanlambirim

\zeta

termal dezavantaj faktörü–

\bar{\phi}_M

yavaşlatıcı-ortalamalı termal akın m^-2 s^-1

\bar{\phi}_F

yakıt-ortalamalı termal akın m^-2 s^-1

V_M

hücredeki yavaşlatıcı hacmim^3

V_F

hücredeki yakıt hacmim^3

Heterojen Termal Kullanım

(10.2)

f = \frac{\Sigma_a^F V_F \bar{\phi}_F}{\Sigma_a^F V_F \bar{\phi}_F + \Sigma_a^M V_M \bar{\phi}_M} = \frac{\Sigma_a^F V_F}{\Sigma_a^F V_F + \Sigma_a^M V_M \,\zeta}

sembolanlambirim

f

termal kullanım–

\Sigma_a^F

yakıtın makroskobik soğurma tesir kesitim^-1

\Sigma_a^M

yavaşlatıcının makroskobik soğurma tesir kesitim^-1

\zeta

termal dezavantaj faktörü–

Dört-Faktör Formülü

(10.3)

k_\infty = \eta\, f\, p\, \varepsilon

sembolanlambirim

k_\infty

sonsuz çarpan faktörü–

\eta

üreme faktörü (yakıtta termal soğurma başına nötron)–

f

termal kullanım–

p

rezonans kaçış olasılığı–

\varepsilon

hızlı fisyon faktörü–

Rezonans Kaçış Olasılığı

(10.4)

p = \exp\!\left[-\frac{N_F V_F\, I_{\mathrm{eff}}}{\xi\,\Sigma_s\,V_M + \xi\,\Sigma_s\,V_F}\right]

sembolanlambirim

p

rezonans kaçış olasılığı–

N_F

yakıttaki rezonans soğurucusunun sayı yoğunluğum^-3

I_{\mathrm{eff}}

etkin rezonans integralibarn

\xi

çarpışma başına ortalama logaritmik enerji kaybı–

\Sigma_s

makroskobik saçılma tesir kesitim^-1

Etkin Rezonans İntegrali (Tanım)

(10.5)

I_{\mathrm{eff}} = \int \sigma_a(E)\,\frac{\phi(E)}{\phi_\infty(E)}\,\frac{dE}{E}

sembolanlambirim

\sigma_a(E)

enerjiye karşı mikroskobik soğurma tesir kesitibarn

\phi(E)

gerçek (öz-perdelenmiş) birim enerji başına akın m^-2 s^-1 eV^-1

\phi_\infty(E)

asimptotik 1/E akısın m^-2 s^-1 eV^-1

E

nötron enerjisieV

S/M'ye Karşı Etkin Rezonans İntegrali (Eşdeğerlik)

(10.6)

I_{\mathrm{eff}} = a + b\,\sqrt{\frac{S}{M}}

sembolanlambirim

I_{\mathrm{eff}}

etkin rezonans integralibarn

S

yakıt yumrusunun yüzey alanım^2

M

yakıt yumrusunun kütlesikg

a,b

soğurucu/geometri için ampirik sabitlerbarn, barn (m^2/kg)^-1/2

Dancoff Düzeltme Faktörü

(10.7)

C = \frac{1}{S}\oint_S P_{\mathrm{esc}\to\mathrm{fuel}}\, dS \;,\qquad \Sigma_e' = \Sigma_e\,(1-C)

sembolanlambirim

C

Dancoff faktörü–

P_{\mathrm{esc}\to\mathrm{fuel}}

çarpışmadan başka bir çubuğa ulaşma olasılığı–

\Sigma_e

kaçış (yüzey) tesir kesiti, yalıtılmış çubukm^-1

\Sigma_e'

Dancoff-düzeltilmiş kaçış tesir kesitim^-1

Wigner Rasyonel Yaklaşımı (Kaçış Olasılığı)

(10.8)

P_{\mathrm{esc}}(E) = \frac{\Sigma_e}{\Sigma_t(E) + \Sigma_e}\;,\qquad \Sigma_e = \frac{1}{\bar{\ell}} = \frac{S}{4 V_F}

sembolanlambirim

P_{\mathrm{esc}}

yumrudan ilk-uçuş kaçış olasılığı–

\Sigma_t(E)

yakıtın toplam makroskobik tesir kesitim^-1

\Sigma_e

kaçış (yüzey) tesir kesitim^-1

\bar{\ell}

yumrunun ortalama kiriş uzunluğu (4V/S)m

Çarpışma Olasılığı Karşılıklılığı

(10.9)

V_i\,\Sigma_{t,i}\,P_{i\to j} = V_j\,\Sigma_{t,j}\,P_{j\to i}

sembolanlambirim

P_{i\to j}

i'de doğan bir nötronun ilk çarpışmayı j'de yapma olasılığı–

V_i

i bölgesinin hacmim^3

\Sigma_{t,i}

i bölgesindeki toplam tesir kesitim^-1

Hızlı Fisyon Faktörü

(10.10)

\varepsilon = \frac{\text{total fissions (fast + thermal)}}{\text{thermal fissions}} \approx 1 + \frac{P_{\mathrm{ff}}\,\nu_{28}\,\sigma_{f,28}}{\dots}

sembolanlambirim

\varepsilon

hızlı fisyon faktörü–

P_{\mathrm{ff}}

bir hızlı nötronun U-238 hızlı fisyonuna yol açma olasılığı–

\nu_{28}

U-238 hızlı fisyonu başına nötron–

\sigma_{f,28}

U-238 hızlı fisyon tesir kesitibarn

Wigner–Seitz Eşdeğer Hücre Yarıçapı

(10.11)

\pi R_{\mathrm{cell}}^2 = A_{\mathrm{cell}} \quad\Rightarrow\quad R_{\mathrm{cell}} = \sqrt{\frac{A_{\mathrm{cell}}}{\pi}} = p_{\mathrm{pitch}}\sqrt{\frac{1}{\pi}}

sembolanlambirim

R_{\mathrm{cell}}

eşdeğer silindirik hücre dış yarıçapım

A_{\mathrm{cell}}

gerçek kafes hücresinin alanım^2

p_{\mathrm{pitch}}

kafes adımı (çubuk merkez aralığı)m

Akı-Hacim Homojenleştirilmiş Grup Sabiti

(10.12)

\bar{\Sigma}_{x,g} = \frac{\displaystyle\sum_i V_i \,\Sigma_{x,g}^{(i)}\,\phi_{g}^{(i)}}{\displaystyle\sum_i V_i\,\phi_{g}^{(i)}}

sembolanlambirim

\bar{\Sigma}_{x,g}

homojenleştirilmiş tesir kesiti, reaksiyon x, grup gm^-1

\Sigma_{x,g}^{(i)}

i bölgesindeki tesir kesitim^-1

\phi_g^{(i)}

i bölgesindeki grup akısın m^-2 s^-1

V_i

i bölgesinin hacmim^3

Dar-Rezonans Öz-Perdelenmiş Akısı

(10.13)

\phi(E) = \frac{\Sigma_p + \Sigma_e}{\Sigma_t(E) + \Sigma_e}\,\frac{1}{E}

sembolanlambirim

\Sigma_p

potansiyel (asimptotik) saçılma tesir kesitim^-1

\Sigma_e

kaçış tesir kesiti S/4Vm^-1

\Sigma_t(E)

E enerjisinde toplam tesir kesitim^-1

E

nötron enerjisieV

Termal Kullanımın ABH Biçimi

(10.14)

\frac{1}{f} = 1 + \frac{\Sigma_a^M V_M}{\Sigma_a^F V_F} + \frac{\Sigma_a^M V_M}{V_F}\,(E-1)\,\frac{1}{\Sigma_a^F V_F / V_F}

sembolanlambirim

f

termal kullanım–

E

ABH fazla-soğurma (akı-bastırma) faktörü–

\Sigma_a^F

yakıt soğurma tesir kesitim^-1

\Sigma_a^M

yavaşlatıcı soğurma tesir kesitim^-1

Deneyin: zincirleme reaksiyon

Etkileşimli · Nötron çoğalmasıCanlı simülasyon

Kontrol çubuklarık = 1.00Kritik

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.