Nükleer Reaktör Korlarının Termal-Hidrolik Analizi
Bir reaktör korunun üretebileceği maksimum güç yalnızca nötron fiziğiyle değil, fisyon ısısını yakıtı eritmeden veya kılıfa zarar vermeden uzaklaştırma yeteneğiyle belirlenir. Bölüm 12, bu sınırları belirleyen termal-hidrolik disiplini geliştirir. Hacimsel fisyon ısı kaynağı q''' (yerel nötron akısı ve fisil yüküyle orantılı) ile başlar; ısıyı UO2 yakıt peletinin içinden radyal iletimle, pelet-kılıf boşluğundan ve kılıftan dışarı izler, ardından akan soğutucuya zorlanmış konveksiyonla aktarır. Her adım bir sıcaklık düşüşü getirir ve düşüşler zinciri, UO2 erime noktasının (~2800°C) altında kalması gereken yakıt merkez ekseni sıcaklığını sabitler. Bölüm, iç ısı üretimli silindirik ısı iletim denklemini formüle eder, parabolik sıcaklık profilini ve yalnızca lineer güç q' ile iletkenliğe bağlı merkez ekseni-yüzey düşüşünü türetir ve güçlü sıcaklığa-bağlı k(T)'yi ele alan iletkenlik-integrali yöntemini tanıtır. Ardından tek-fazlı zorlanmış konveksiyonu Nusselt sayısı ve boyutsuz korelasyonlar (su için Dittus-Boelter, sıvı metaller için Lyon-Martinelli) aracılığıyla; kanal basınç düşüşünü ise Fanning sürtünme faktörü ve Moody diyagramı üzerinden işler. Bölümün kalbi kaynamalı ısı transferidir: kaynama eğrisi, çekirdeksel kaynama, kritik ısı akısı ve çekirdeksel kaynamadan ayrılma (DNB), kuruluk, boşluk oranı ve kayma oranı; bunların hepsi KSR (BWR) ve BSR (PWR) güvenliğini yönetir. Son olarak hidrodinamiği (tek- ve iki-fazlı basınç düşüşü, paralel kanallar arasında akış dağılımı) termal analizle birleştirir, sıcak-kanal faktörlerini ve DNB oranını tanımlar ve termal-hidrolik sınırların reaktör gücünü ve kor tasarımını nasıl kısıtladığını gösterir.
Giriş
Bu kısım, termal-hidroliği nötroniğin mühendislik karşılığı olarak çerçeveler: bir reaktör temelde bir ısı kaynağıdır ve kullanılabilir gücü, ısının ne kadar hızlı ve güvenli uzaklaştırılabileceğiyle sınırlıdır. Bir korun merkezi termal sınırlamalarını belirler—yakıt merkez ekseni sıcaklığı yakıt erime noktasının altında, yüzey ısı akısı kritik ısı akısının (burnout/DNB) altında kalmalıdır—çünkü ikisinden birinin aşılması yakıt erimesine veya kılıf bozulmasına yol açar. Kısım, enerjinin fisyondan başlayarak izlediği yolu çizer: enerji çoğunlukla yakıtta fisyon parçacığı kinetik enerjisi olarak hacimsel ısı kaynağı biçiminde birikir ve pelet, boşluk, kılıf üzerinden soğutucuya geçer. En sıcak yakıt elemanını (sıcak kanal) analiz etme tasarım felsefesini ve ortalama ile maksimum koşulları ilişkilendiren güç-tepe ve sıcak-kanal faktörlerinin kullanımını tanıtır. Problemin bağlaşık doğasını gerekçelendirir: güç dağılımı (nötronik) ısı kaynağını belirlerken, soğutucu sıcaklığı ve yoğunluğu reaktivite katsayıları yoluyla nötroniğe geri besler; bu nedenle termal-hidrolik ve nükleer analizler, koru boyutlandırmak ve maksimum izin verilen gücü belirlemek için nihayetinde birlikte çözülmelidir.
Nükleer Reaktör Korlarında Güç Üretimi
Fisyon olay başına yaklaşık 200 MeV salar, ancak yalnızca yakıtta fisyon parçacığı kinetik enerjisi olarak yerel olarak biriken kısım (kabaca 180-190 MeV; gama ve nötron transportundan başka yerlerde geri kazanılabilen birkaç yüzdeyle) yakıt-içi ısı kaynağı olur. Bu kısım, hacimsel ısı üretim hızı q'''yi (W/cm^3) fisyon hızı yoğunluğuyla orantılı olarak q''' = E_f * Sigma_f * phi şeklinde tanımlar; yüzey ısı akısı q'' (W/cm^2) ve lineer güç q'yü (W/cm) tanımlar ve bunları yakıt geometrisi üzerinden integrasyonla ilişkilendirir. q''' nötron akısını izlediğinden, güç dağılımı kor merkezine yakın tepe yapar (çıplak silindir için eksenel kosinüs, radyal J0 Bessel); bu nedenle bölüm eksenel ve radyal tepe faktörlerini ve nükleer tepeyi imalat ve akış belirsizlikleriyle birleştiren mühendislik sıcak-kanal faktörlerini tanıtır. Kısım ayrıca, kapanıştan sonra süren ve soğutucu-kaybı güvenlik analizine egemen olan fisyon-ürünü bozunmasından gelen bozunma-ısısı katkısını belirtir. Temel mühendislik nicelikleri—toplam kor termal gücü, ortalama ve tepe güç yoğunluğu ve ortalama lineer ısı hızı—önceki bölümlerin çok-gruplu difüzyon yöntemleriyle hesaplanan uzaysal akı şekline doğrudan bağlanır ve böylece güç üretimi nötronikten ısı transferine açık köprüyü oluşturur.
Reaktör Yakıt Elemanlarında Radyal Isı İletimi
Hacimsel kaynak bilindiğinde, bu kısım silindirik yakıt çubuğu boyunca kararlı-durum radyal ısı iletimini çözer. Eksenel iletim ihmal edilerek (çubuk uzun ve ince olduğundan geçerlidir), silindirik koordinatlarda düzgün kaynak q''' içeren iletim denklemi parabolik bir sıcaklık profiline integre olur ve yakıtta merkez ekseni-yüzey sıcaklık düşüşünü verir: T_0 - T_s = q'/(4*pi*k_f); bu, pelet yarıçapına değil yalnızca lineer güce ve iletkenliğe bağlı dikkat çekici bir sonuçtur. UO2 iletkenliği sıcaklıkla güçlü biçimde düştüğünden, kısım iletkenlik integralini (k(T) dT integrali) tanıtır; böylece erime sınırı, mutlak sıcaklık yasasından bağımsız maksimum izin verilen lineer güç olarak ifade edilebilir. Ardından ek sıcaklık düşüşlerini seri olarak zincirler: pelet-kılıf boşluğu boyunca (iletim, ışınım ve teması toplayan boşluk iletkenliği h_g ile yönetilir), kılıf duvarı boyunca (kaynaksız silindirik iletim) ve son olarak soğutucuya film düşüşü. Bunları soğutucudan merkez eksenine toplamak, verilen bir soğutucu sıcaklığı ve ısı akısı için yakıt merkez ekseni sıcaklığını verir. Kısım, merkez eksenini erimenin altında tutan maksimum izin verilen lineer ısı hızını (UO2 için yaklaşık 660 W/cm) elde etmek üzere ilişkiyi tersine çevirerek kapanır; bu, oksit yakıtlı reaktörler için en önemli termal tasarım sınırıdır.
Tek-Fazlı Soğutucularda Zorlanmış Konveksiyonla Isı Transferi
Bu kısım, ısıyı kılıf yüzeyinden akan soğutucuya zorlanmış konveksiyonla aktarır; bu, q'' = h(T_w - T_b) Newton soğuma yasasıyla tanımlanır; burada h taşınımsal ısı transfer katsayısıdır. Türbülanslı akış için h analitik olarak türetilemediğinden, kısım onu boyutsuz Nusselt sayısı Nu = h*D_e/k aracılığıyla ifade etmek için boyut analizini kullanır; bu Reynolds sayısı Re (atalet-viskoz kuvvet oranı) ve Prandtl sayısı Pr (momentum-termal yayınırlık oranı) ile korele edilir. Temel korelasyonlar sunulur: su gibi sıradan akışkanlar (Pr ~ 1) için Dittus-Boelter korelasyonu Nu = 0.023 Re^0.8 Pr^0.4 ve iletimin egemen olduğu sıvı metaller (düşük Pr) için Lyon-Martinelli korelasyonu Nu = 7 + 0.025 Pe^0.8. Çubuk demetinin daire-dışı alt kanallarını ele almak için hidrolik çap D_e tanıtılır. Kısım, kanal boyunca bir enerji dengesinden soğutucu yığın sıcaklığının eksenel artışını hesaplar ve bunu iletim düşüşleriyle birleştirerek sıcak kanal boyunca soğutucu girişinden yakıt merkez eksenine kadar tam sıcaklık dağılımını elde eder. Ayrıca momentum tarafını—Fanning sürtünme faktörü f, Blasius ve Moody korelasyonları ve ortaya çıkan duvar-sürtünme basınç gradyanını—sunarak hidrodinamik kısmında detaylandırılan pompalama-gücü ve basınç-düşüşü kısıtlarını hazırlar.
Nükleer Reaktör Korlarında Kaynamalı Isı Transferi
Kaynama, ısı uzaklaştırmayı çarpıcı biçimde artırır ve KSR'lerde (BWR) kasıtlıdır, BSR'lerde (PWR) yerel olarak mümkündür, ancak kaynama krizi tehlikesini taşır. Bu kısım kaynama eğrisini (ısı akısına karşı duvar aşırı ısınması) izler: doğal konveksiyon, kabarcıkların duvarda oluşup h'yi keskin biçimde yükselttiği çekirdeksel kaynamanın başlangıcı, tam gelişmiş çekirdeksel kaynama, kritik ısı akısı (CHF) denilen tepe ve buhar örtüsünün duvarı yalıttığı ve duvar sıcaklığının sıçradığı, burnout ve kılıf bozulmasına yol açan çekirdeksel kaynamadan ayrılmada (DNB) film kaynamasına felaket geçiş. Aşırı-soğutulmuş/düşük-kuruluk DNB'sini (kabarcık-sınırı kopması, BSR endişesi) yüksek-kuruluk kurumasından (sıvı-film kayboluşu, KSR endişesi) ayırır. İki-fazlı akışın termodinamik ve akış niceliklerini tanımlar: termodinamik denge kuruluğu x, boşluk oranı alpha (buhar hacim oranı) ve daha hızlı buhar hızını sıvıya bağlayan kayma oranı S = v_g/v_l; homojen model S=1 limitidir. Çekirdeksel-kaynama aşırı ısınması için Jens-Lottes korelasyonu ve ampirik CHF korelasyonları (W-3 tipi gibi) tanıtılır. Kısım, kılıfın asla kurumadığını garanti etmek için aşırı-güç geçişlerinde bile yaklaşık 1.3'ün üzerinde tutulması gereken, kritik ısı akısının gerçek yerel akıya oranı olan DNB oranıyla (DNBR) doruğa ulaşır.
Hidrodinamik Kor Analizi
Isı uzaklaştırma, soğutucuyu kor boyunca hareket ettirmeyi gerektirir ve bu kısım pompaların yenmesi gereken basınç düşüşünü ve akışın paralel kanallar arasında nasıl dağıldığını analiz eder. Toplam kanal basınç düşüşü dört katkıya ayrılır: Fanning sürtünme faktörü ve Moody diyagramı kullanan sürtünme (duvar-kayma) düşüşü, ızgaralar, ara parçalar ve alan değişimlerindeki biçim/yerel kayıplar (kayıp katsayıları K ile), kanal boyunca soğutucu yoğunluğu değişiminden gelen ivme düşüşü ve yerçekimsel (yükseklik) yükü. İki-fazlı bölgeler için kısım, iki-fazlı sürtünme çarpanını ve Martinelli-tipi formülasyonu tanıtır, çünkü buhar belirir belirmez sürtünme ve ivme düşüşleri keskin biçimde yükselir. Kritik olarak, bir kordaki tüm paralel kanallar aynı giriş ve çıkış toplama haznelerini paylaşır, bu yüzden aynı toplam basınç düşüşünü yaşar; akış, her kanalın basınç-düşüşü ilişkisi bu ortak sınır koşuluyla eşleşecek şekilde kendiliğinden dağılır. Daha sıcak kanallar daha çok boşluk üretir, iki-fazlı dirençlerini artırır ve dolayısıyla daha az akış alır—bu geri besleme, akış yeniden dağılımına ve kaynamalı sistemlerde yoğunluk-dalgası ve Ledinegg akış kararsızlıklarına yol açabilir. Kısım, termal analizle birlikte korun çalışma noktasını belirleyen pompalama-gücü kısıtını ve kanal akış-denge denklemlerini kurar.
Termal-Hidrolik Kor Analizi
Son kısım, iletim, konveksiyon, kaynama ve hidrodinamiği, maksimum izin verilen reaktör gücünü belirleyen bağlaşık bir kor analizinde birleştirir. Tek-kanal modelini kurar: nötronikten eksenel güç şeklini al, yığın sıcaklık artışını elde etmek için soğutucu enerji dengesini ilerlet, kanal boyunca kılıf-yüzey ve yakıt-merkez ekseni sıcaklıklarını elde etmek için konveksiyon ve iletim düşüşlerini ekle ve kaynamanın olduğu yerde kuruluk ile boşluğu izle. Analiz daha sonra, nominal nükleer tepeyi imalat toleransları, akış düzensizliği ve belirsizlik için paylarla çarpan sıcak-kanal (mühendislik) faktörleriyle belirlenen, en yüksek entalpi artışı ve ısı akısına sahip kanal olan sıcak kanala uygulanır. Aynı anda iki sınır dayatılır: yakıt-merkez ekseni sıcaklığı erimenin altında kalmalı ve minimum DNBR tasarım sınırını (aşırı-güçte ~1.3) aşmalıdır. Hangi sınıra önce ulaşılırsa kor gücünü sınırlar. Kısım, nötronikle çift-yönlü bağlaşımı vurgulayarak kapanır—soğutucu sıcaklığı ve boşluğu, yavaşlatıcı ve Doppler reaktivite katsayıları yoluyla geri besler—böylece öz-tutarlı bir kor tasarımı akı/güç ile sıcaklık/yoğunluk alanları arasında yinelenir; bu, modern bağlaşık nötronik-termal-hidrolik kor benzeticilerinde uygulanan döngünün ta kendisidir.
Temel Denklemler (19)
Bu denklemi çalıştır: Yakıt Çubuğu Radyal Sıcaklık Profili ve Erime Sınırı →Yerel birim hacim başına ısı üretim hızını, fisyon başına geri kazanılabilir enerji ile fisyon reaksiyon hızının çarpımı olarak verir. Nötroniği (akı) termal analize bağlayan kaynak terimidir ve uzaysal güç şeklini izler.
Birim uzunluk başına lineer güç q' ile birim alan başına yüzey ısı akısı q''yü, düzgün üretimli a yarıçaplı silindirik bir yakıt çubuğu için hacimsel kaynağa ilişkilendirir. Bu üç güç ölçüsü termal tasarım boyunca kullanılır.
İç ısı üretimli bir yakıt çubuğu için, eksenel iletimi ihmal eden kararlı-durum radyal iletim denklemi. İntegre edildiğinde yakıt boyunca sıcaklık profilini verir; çubuk-içi tüm sıcaklık düşüşlerinin başlangıç noktasıdır.
Sabit iletkenlik için pelet yüzeyinin üzerindeki maksimum (merkez ekseni) sıcaklık artışı. Dikkat çekici biçimde yalnızca lineer güce ve iletkenliğe bağlıdır, pelet yarıçapına değil; bu da q'yü doğal termal tasarım değişkeni yapar.
İletkenlik sıcaklıkla güçlü biçimde değiştiğinde (UO2'de olduğu gibi) merkez ekseni sıcaklık düşüşünü genelleştirir. k(T)'nin yüzeyden merkez eksenine integrali q'/(4 pi)'ye eşittir; erime sınırı böylece maksimum izin verilen lineer gücü sabitler.
İletim, ışınım ve teması toplayan bir boşluk iletkenliği h_g ile modellenen, pelet ile kılıf arasındaki gazla dolu boşluk boyunca sıcaklık düşüşü. Boşluk, yakıt performansında büyük ve belirsiz bir termal dirençtir.
Kılıf duvarı boyunca kaynaksız silindirik iletim. Soğutucudan merkez eksenine seri zincire kılıf termal direncini ekler; yüksek iletkenlikli ince kılıf (Zircaloy) bu düşüşü küçük tutar.
Kılıf dış yüzeyinden soğutucu yığınına taşınımsal ısı transferi; h taşınımsal katsayıdır. Film sıcaklık düşüşünü tanımlar ve iletimi soğutucuya bağlayan sınır koşuludur.
Su gibi sıradan soğutucularda (Pr ~ 1) Nusselt sayısı (dolayısıyla h) için ampirik türbülanslı-akış korelasyonu; ısıtmada n=0.4. BSR/KSR alt kanallarında taşınımsal katsayıyı hesaplamak için standart araç.
Moleküler iletimin önemli olduğu düşük-Prandtl sıvı-metal soğutucular (örn. LMFBR'lerde sodyum) için ısı-transfer korelasyonu; Peclet sayısıyla korele edilir. Sıvı metallere özgü yüksek h'yi öngörür.
Soğutucu yığın sıcaklığını kanal boyunca ilerleten enerji dengesi; kütle-akışı çarpı ısı-kapasitesi lineer gücü uzaklaştırır. Kosinüs eksenel güç şekliyle integre edildiğinde karakteristik S-biçimli soğutucu sıcaklık profilini verir.
Bir kanal boyunca duvar-sürtünme basınç kaybı; Fanning sürtünme faktörü f, Blasius korelasyonundan veya Moody diyagramından elde edilir. Pompalama basınç düşüşüne tek-fazlı baskın katkı.
Akış rejimini yöneten atalet-viskoz kuvvet oranı; kanallarda türbülanslı akış Re ~ 2100'ün üzerinde başlar. Hem sürtünme faktörüne hem de ısı-transfer korelasyonlarına girer.
Buhar kütle oranını, yerel entalpiyi doymuş-sıvı ve gizli-ısı değerlerine göre tanımlar. Kuruluk, kaynayan bir kanal boyunca artar ve boşluk oranını, iki-fazlı basınç düşüşünü ve DNB davranışını yönlendirir.
Buhar hacim oranını (boşluk) akış kuruluğuna, faz yoğunluk oranına ve buhar ile sıvı hızları arasındaki kayma oranı S'ye ilişkilendirir. S=1 ile homojen modele indirgenir; boşluk oranı hem basınç düşüşü hem de kaynayan korlardaki reaktivite geri beslemesi için temeldir.
Suyun aşırı-soğutulmuş/doymuş çekirdeksel kaynamasında duvar aşırı ısınmasını ısı akısı ve basıncın fonksiyonu olarak veren ampirik ilişki (SI-biçimi, p MPa cinsinden). Kaynama başladığında kılıf-yüzey sıcaklığını tahmin etmek için kullanılır.
Yerel kritik ısı akısının (bir CHF korelasyonundan) gerçek yerel ısı akısına bölünmesiyle tanımlanan güvenlik payı. Sıcak kanal boyunca minimum DNBR, burnout'u önlemek için aşırı-güç geçişlerinde bile tasarım sınırının (~1.3) üzerinde kalmalıdır.
Kanal basınç düşüşünü sürtünme, biçim (ızgaralar/ara parçalar), ivme (yoğunluk değişimi) ve yerçekimi terimlerine ayırır. Tüm paralel kanallar bu toplam düşüşü paylaşır; bu da akış dağılımını ve gerekli pompalama gücünü belirler.
Maksimum soğutucu entalpi artışının (sıcak kanalda) kor-ortalama entalpi artışına oranı. Isı-akısı sıcak-kanal faktörüyle birlikte, ortalama çalışma koşullarını güç sınırını belirlemek için kullanılan sınırlayıcı yerel koşullara ilişkilendirir.
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 lineer gücü ve soğutucu koşullarını ayarlar ve sıcaklığın soğutucu yığınından başlayarak film, kılıf, boşluk üzerinden ve pelet boyunca merkez eksenine kadar ilerleyişini izler. Canlı bir radyal T(r) eğrisi ve çubuğun renkli bir kesiti gerçek zamanlı güncellenir; merkez ekseni UO2 erime noktasına yaklaştığında çubuk kırmızıya döner ve bir erime-sınırı uyarısı tetiklenir, böylece q' tavanı somutlaşır.
Öğrenci kosinüs eksenel güç şekline sahip tek bir ısıtılmış kanalı sürer ve soğutucunun ısınmasını, çekirdeksel kaynamanın başlangıcına ulaşmasını ve z boyunca kuruluk ile boşluk geliştirmesini izler. Yığılmış eksenel grafikler T_b(z), x(z), alpha(z), ısı akısı q''(z) ve kritik ısı akısı zarfını gösterir; minimum DNBR vurgulanır ve 1.3'ün altına düşerse kırmızıya döner. Kütle akışını veya gücü kaydırmak, akış kıtlığının kanalı kaynama krizine ittiğini gösterir.
Öğrenci genel kor gücünü yükseltir ve sıcak kanalda iki termal sınırın aynı anda yarıştığını izler: erimeye karşı yükselen yakıt-merkez ekseni sıcaklığı ve sınırına karşı düşen minimum DNBR. İki-ibreli bir 'sınır göstergesi' hangi kısıtın bağlayıcı olduğunu gösterir ve maksimum izin verilen güç raporlanır. Sıcak-kanal faktörlerini ve giriş akışını değiştirmek, mühendislik paylarının ve pompalamanın elde edilebilir güçle nasıl takas edildiğini gösterir.
Sayısal veriler
Bu bölüm öğrencilere reaktör analizinin termal tarafını kodlamayı öğretir. Bir 1B silindirik ısı-iletim çözücüsü (analitik iletkenlik-integrali artı, sıcaklığa-bağlı k(T) Newton iterasyonuyla ele alınan bir sonlu-fark/sonlu-hacim sürümü) uygularlar, ardından boşluk, kılıf ve taşınımsal film için seri termal dirençleri zincirlerler. Akışkan tarafında boyutsuz-sayı değerlendirmesini (Re, Pr, Pe, Nu) ve Dittus-Boelter ile Lyon-Martinelli korelasyonlarını, ayrıca bir Fanning-sürtünme / Moody basınç-düşüşü rutinini kodlarlar. Kaynama modülü öğrencilerden entalpi, denge kuruluğu ve kayma-düzeltmeli boşluk oranını hesaplamak için bir soğutucu enerji dengesini eksenel olarak ilerletmelerini (açık Euler veya RK4), ardından CHF korelasyonlarını değerlendirmelerini ve DNBR'yi kanal boyunca minimize etmelerini ister. Doruk noktası, erime ve DNBR sınırlarına tabi maksimum izin verilen kor gücünü döndüren, ikiye-bölme kök-bulmalı bağlaşık bir tek-kanal çözücüdür. Öğrenciler üretim kodlarının nasıl düzenlendiğini öğrenir: COBRA-IV / VIPRE / SUBCHANNEL gibi alt-kanal termal-hidrolik kodları, RELAP5 ve TRACE gibi sistem kodları ve bunların tam kor analizi için nötronikle nasıl bağlaştığı. Vurgu, sayısal kararlı ilerletme şemaları, doğrusal-olmayan özellik tabloları altında gürbüz kök-bulma, tablo-tabanlı buhar özellikleri ve ortak bir basınç sınır koşulunu paylaşan paralel-kanal ağı için veri yapıları üzerindedir.