Termal Spektrum Hesaplamaları ve Termal Grup Sabitleri
Bu bölüm, termal enerji aralığında (E < ~1 eV) nötron enerji spektrumunu belirlemek ve sürekli enerjili tesir kesitlerini çok-gruplu difüzyon teorisinin gerektirdiği az-gruplu sabitlere indirgemek için gereken fiziği ve hesaplama altyapısını geliştirir. Hedef çekirdeklerin durağan kabul edildiği ve saçılmanın yalnızca aşağı-saçılma olduğu yavaşlama bölgesinin aksine, termal nötronlar termal hareket içindeki ve moleküllerde ile kristal kafeslerde kimyasal olarak bağlı çekirdeklerle enerji alışverişi yapar. Sonuç olarak hem yukarı hem aşağı saçılma meydana gelir, spektrum bir Maxwell-Boltzmann dağılımına yaklaşır ancak hiçbir zaman tam olarak ona eşit olmaz ve saçılma çekirdeği ayrıntılı denge (detailed balance) koşulunu sağlamalıdır. Bölüm önce ölçülen termal spektrumların genel özelliklerini açıklar: sızıntı ve 1/v soğurma tarafından 'sertleştirilen' (daha yüksek enerjiye kaydırılan) Maxwell'e yakın bir tepe ve yavaşlama bölgesine bağlanan 1/E birleştirme kuyruğu. Ardından sonsuz-ortam termal spektrum denklemini türetir, serbest gaz (proton-gazı ve ağır-gaz) saçılma modellerini tanıtır ve iki temel yaklaşık modeli sunar: ağır-gaz (diferansiyel) yaklaşımı ve Wigner-Wilkins proton-gazı modeli ile basit etkin-nötron-sıcaklığı reçeteleri. Son olarak, gerçekçi spektrumların tablolanmış saçılma çekirdekleri S(α,β) kullanılarak integral taşınım-denge denkleminin enerjide ayrıklaştırılmasıyla sayısal olarak nasıl hesaplandığını ve elde edilen akının akı-ağırlıklı ortalama ile termal grup sabitleri üretmek için nasıl kullanıldığını açıklar. Bu termal grup sabitleri doğrudan Bölüm 3'ün az-gruplu difüzyon hesaplarına girer ve mikroskobik nötron-çekirdek fiziği ile tam-çekirdek reaktör analizi arasındaki döngüyü kapatır.
Termal Nötron Spektrumlarının Genel Özellikleri
Yaklaşık 1 eV'nin altında nötronlar saçılma çarpışmalarında enerji kaybetmenin yanı sıra enerji kazanabilir de, çünkü hedef çekirdekler termal hareket içindedir ve moleküllerde veya kafeslerde bağlıdır. Varsayımsal sonsuz, soğurmasız bir ortamda nötron spektrumu, ayrıntılı dengenin gerektirdiği denge olan moderatör sıcaklığı T'deki Maxwell-Boltzmann dağılımı M(E,T)'ye tam olarak gevşerdi. Gerçek bir reaktörde iki etki bu ideali bozar: (i) 1/v soğurma tercihen yavaş nötronları uzaklaştırır ve (ii) sürekli yavaşlama kaynağı yukarıdan nötron besler. Her ikisi de spektral tepeyi sıcaklığın ima ettiği değerin üzerine iter; bu olaya spektral sertleşme denir ve genellikle etkin nötron sıcaklığı T_n > T_moderatör ile modellenir. Tam termal spektrum böylece, birkaç kT civarındaki bir geçiş (birleştirme) bölgesi aracılığıyla 1/E yavaşlama kuyruğuna düzgünce birleştirilen, sertleşmiş, Maxwell'e yakın bir termal tepeden oluşur. Bu özellikleri anlamak, bölümün geri kalanında geliştirilen termalleşme modellerini ve grup-sabiti reçetelerini motive eder.
Nötron Termalleşmesinin Yaklaşık Modelleri
Tam saçılma çekirdekleri karmaşık olduğundan, bölüm işlenebilir yaklaşımlar geliştirir. Serbest gaz modeli, moderatörü Maxwell-Boltzmann hız dağılımına sahip çekirdeklerden oluşan ideal tek atomlu bir gaz olarak ele alır; saçılma çekirdeği kapalı biçimde yazılabilir ve iki sınırlayıcı modelin temelini oluşturur. Ağır-gaz modeli, her çarpışmanın az enerji aktardığı şekilde moderatör kütlesinin A >> 1 olduğunu varsayar; saçılma integralini 1/A'da en düşük mertebeye açmak, integral denge denklemini, denge çözümü tam olarak Maxwellian olan ikinci mertebeden bir diferansiyel denkleme (ağır-gaz veya Fokker-Planck termalleşme denklemi) dönüştürür. Wigner-Wilkins modeli, bir proton gazının (A=1, hidrojenli moderatörler) zıt uç durumunu ele alır ve spektrumu elde etmek için integre edilebilen birinci mertebeden bir diferansiyel denklem verir. En basit mühendislik reçeteleri, moderatör sıcaklığı ile soğurma-saçılma oranına bağlı bir etkin sıcaklık T_n'de sertleşmiş bir Maxwellian varsayar ve 1/E kuyruğuna geçiş yapan bir birleştirme fonksiyonu kullanır. Bu modeller titizliği analitik içgörü için takas eder ve hâlâ hızlı tahminler ve sayısal kodların kıyaslanması için kullanılır.
Termal Nötron Spektrumlarının Genel Hesaplamaları
Kimyasal bağlanma, moleküler titreşim/dönme ve kristal kafes dinamiğinin önemli olduğu gerçekçi moderatörler için saçılma çekirdeği serbest gazla yaklaştırılamaz. Bunun yerine, malzemenin dinamik yapı faktöründen türetilen boyutsuz momentum transferi α ve enerji transferi β'nin bir fonksiyonu olan saçılma yasası S(α,β)'dan hesaplanır; bu, termal-saçılma ENDF/B dosyalarında tablolanan ve LEAPR/NJOY gibi kodlar tarafından su, grafit, berilyum vb. için çekirdeklere işlenen veridir. Sonsuz-ortam termal denge denklemi daha sonra enerjiyi ince bir termal grup yapısına ayrıklaştırarak, çekirdekten yoğun saçılma transfer matrisini oluşturarak (ki bu otomatik olarak ayrıntılı dengeyi sağlar) ve ortaya çıkan doğrusal sistemi ayrık akı için çözerek sayısal olarak çözülür. GATHER, SOFOCATE, THERMOS gibi kodlar ve WIMS/NJOY'un termal modülleri bunu uygular. Yakınsamış termal akı φ(E) nihayetinde, çok-gruplu difüzyon denklemlerinde kullanılan Σ_a, Σ_s, νΣ_f ve D değerlerini üreterek, akı-ağırlıklı ortalama ile noktasal tesir kesitlerini geniş-grup termal grup sabitlerine indirgemek için kullanılır.
Temel Denklemler (13)
Bu denklemi çalıştır: Sonsuz-Ortam Termal Spektrum Çözücü →Soğurmasız sonsuz bir ortamın T sıcaklığında gevşediği birim-enerji başına akı spektrumudur. Termalleşme denkleminin denge çözümü ve spektral sertleşmenin ölçüldüğü referanstır.
Termal dengedeki serbest moderatör çekirdeklerinin hız dağılımıdır; serbest-gaz saçılma çekirdeğini oluşturmak ve termal olarak ortalanmış tesir kesitlerini hesaplamak için kullanılır.
Sızıntısız, sonsuz homojen bir ortamda enerjiye bağlı nötron dengesidir. Yavaşlamanın aksine, saçılma integrali tüm E' üzerinden (hem yukarı hem aşağı saçılma) çalışır. S(E), termal aralığı yukarıdan besleyen yavaşlama kaynağıdır. Bu, tüm termalleşme modellerinin çözdüğü ana denklemdir.
Herhangi bir fiziksel termal saçılma çekirdeğinin sağlaması gereken temel kısıttır. Soğurma ve kaynak yokluğunda Maxwellian M(E,T)'nin termalleşme denkleminin denge çözümü olmasını garanti eder ve yukarı-saçılma ile aşağı-saçılma hızlarını bağlar.
Maxwellian hareket içindeki çekirdeklerden (kütle numarası A) oluşan ideal tek atomlu bir gazdan saçılma için kapalı biçimli çekirdektir. Hata-fonksiyonu biçimi (burada şematik) hem aşağı-saçılmayı (E<E') hem yukarı-saçılmayı (E>E') ele alır ve E' >> k_B T için durağan-çekirdek sonucuna indirgenir. Ağır-gaz ve Wigner-Wilkins limitlerinin başlangıç noktasıdır.
İntegral termalleşme denklemi, ağır-kütle limitinde (A >> 1, çarpışma başına küçük enerji transferi) bu ikinci mertebeden Fokker-Planck tipi adi diferansiyel denkleme indirgenir. Kaynaksız, soğurmasız çözümü Maxwellian'dır; soğurma ve kaynak terimleri sertleşmiş spektrumu üretir. Basit sonlu farklarla çözülebilir.
Wigner ve Wilkins tarafından serbest proton gazı (A=1) için elde edilen tam diferansiyel denklemdir. Boyutsuz enerji x = E/k_BT ve yoğunluk ψ = φ/M cinsinden, hidrojenli moderatörlerde (su) termal spektrumu vermek için sayısal olarak integre edilebilir. SOFOCATE tipi kodlarda ve kıyaslama olarak yaygın kullanılır.
Westcott/Coveyou tipi ampirik reçete: termal tepe, moderatör sıcaklığı T_m'yi soğurma-moderasyon oranıyla orantılı olarak aşan bir etkin sıcaklık T_n'deki bir Maxwellian ile yaklaştırılır. C sabiti 1 mertebesindedir (hidrojen baskın sistemler için yaygın olarak ~1.4). Hızlı spektrum tahminleri verir.
Çift-diferansiyel bağlı-saçılma tesir kesiti, moderatörün dinamiğini (bağlanma, titreşim, kafes) kapsayan simetrik saçılma yasası S(α,β) aracılığıyla ifade edilir. α boyutsuz momentum transferi, β boyutsuz enerji transferidir. S(α,β) ENDF termal dosyalarında tablolanır ve titiz termal çekirdek üretiminin temel girdisidir.
Saçılma yasasının bağlı olduğu boyutsuz değişkenlerin tanımları. μ, laboratuvar çerçevesindeki saçılma kosinüsüdür. Bu değişkenler serbest-gaz limitinde S(α,β)'yı neredeyse malzemeden bağımsız yapar ve termal saçılma verilerini tablolamak ve enterpole etmek için doğal koordinatlardır.
Herhangi bir noktasal reaksiyon tesir kesiti Σ_x(E)'nin, termal grup g'deki x reaksiyonu için bir grup sabitine akı-ağırlıklı indirgemesi. Hesaplanan termal spektrum φ(E) ile uygulanarak çok-gruplu difüzyon denklemleri için Σ_a,g, Σ_s,g, νΣ_f,g ve (akım-ağırlıklı biçim aracılığıyla) D_g üretilir.
1/(3Σ_tr)'nin akı-ortalaması ile elde edilen termal grup difüzyon katsayısı; burada Σ_tr, ortalama saçılma kosinüsü ile düzeltilmiş taşınım tesir kesitidir. Termal aralıkta hem anizotropi hem kimyasal bağlanma Σ_tr'yi etkiler, dolayısıyla D_g termal spektruma bağlıdır.
Tam düşük-enerji spektrumunun mühendislik gösterimi: etkin sıcaklık T_n'de bir Maxwellian termal tepe artı, birkaç kT'nin üzerinde kuyruğu açan bir birleştirme fonksiyonu Δ(E) ile harmanlanan 1/E epitermal (yavaşlama) kuyruğu. Westcott g- ve s-faktör tesir kesiti konvansiyonlarının temelini oluşturur.
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ı
Öğrenci moderatör türünü, sıcaklığı ve soğurma gücünü ayarlar, ardından denge denklemi çözülürken kararlı-durum termal akı φ(E)'nin ortaya çıkışını izler. Soğurma arttıkça Maxwellian tepenin daha yüksek enerjiye kaydığını (spektral sertleşme) ve yavaşlama kaynağı yükseldikçe 1/E kuyruğunun belirdiğini görür. Canlı bir bindirme, hesaplanan spektrumu moderatör sıcaklığındaki ve uydurulmuş etkin sıcaklıktaki saf Maxwellian ile karşılaştırır.
Hidrojenli moderatörlere odaklanan bu simülasyon, Wigner-Wilkins adi diferansiyel denklemini boyutsuz enerjide integre eder ve indirgenmiş akı ψ(x)'in 1/v soğurma gücünün bir fonksiyonu olarak birimden (saf Maxwellian) nasıl saptığını gösterir. Öğrenciler tam Wigner-Wilkins sonucunu ağır-gaz yaklaşımı ve basit etkin-sıcaklık Maxwellian'ı ile karşılaştırarak her modelin nerede başarılı veya başarısız olduğunu görür.
Çözücüden gelen spektrumu kullanarak öğrenci, akı-ağırlıklı ortalama ile noktasal tesir kesitlerini bir- veya az-gruplu termal sabitlere indirger. Grup sınırlarını seçer, bir ağırlıklandırma spektrumu seçer (gerçek φ vs. düz Maxwellian vs. sabit) ve Σ_a, νΣ_f ve D'nin nasıl değiştiğini görür; bu, ağırlık olarak ders kitabı Maxwellian'ının değil, öz-tutarlı termal akının neden kullanılması gerektiğini çarpıcı biçimde gösterir.
Sayısal veriler
Bu bölüm, öğrencilerin nötron fiziğini üretim reaktör-fiziği kodlarına ilk kez bağladıkları yerdir. Temel algoritmik beceriler şunlardır: (1) ayrıklaştırılmış sonsuz-ortam termalleşme denge denklemini kurma ve çözme; bu, LU ayrıştırması veya Gauss-Seidel/SOR iterasyonu ile çözülen yoğun (serbest-gaz) veya üç-köşegenli (ağır-gaz) bir doğrusal sisteme dönüşür; (2) ayrıntılı dengeye uyan termal saçılma kernelleri oluşturma; bu, ya analitik serbest-gaz formüllerinden ya da tablolanmış S(α,β) saçılma-yasası verilerini okuyup enterpole ederek yapılır; ve (3) sayısal kareleme yoluyla akı-ağırlıklı grup-sabiti üretimi. Bunlar gerçek kodları yansıtır: NJOY'un THERMR/LEAPR modülleri S(α,β)'dan termal kerneller üretir; GATHER, SOFOCATE ve THERMOS termal spektrumu çözer (THERMOS, ayrık kernelli integral taşınımı kullanır); ve kafes kodları WIMS, CASMO ve SCALE/XSDRN elde edilen spektrumları çok-gruplu kütüphaneleri indirgemek için kullanır. Öğrenciler kernel-matris montajını, yakınsama izlemeli iteratif bir doğrusal çözücüyü, Wigner-Wilkins modeli için ODE integratörlerini (RK4/atış) ve ortalama için Gauss/Simpson karelemesini kodlamayı öğrenir. Ayrıca sayısal kararlılığın neden ayrıntılı dengenin tam olarak uygulanmasını gerektirdiğini (aksi takdirde soğurmasız limit Maxwellian'dan kayar) ve enerji-ızgara çözünürlüğünün yukarı-saçılma yakalama doğruluğunu nasıl kontrol ettiğini öğrenir. Rust/WASM'daki uygulama matris çözümü için neredeyse yerel hız verirken, TypeScript katmanı interaktif çizim ve parametre taramalarını yönetir.