AIREK-MOD kodunun yeni versiyonu ve TR-2 transiyent analizlerinde uygulanışı

T Ü R KÎYE ATOM ENERJİSİ K U R U M U ÇEKMECE NÜKLEER ARAŞTIRMA VE EĞİTİM MERKEZİ

Teknik Rapor No: 27

AİRE K — MOD KODUN UN YENİ VERSİYON U VE T R —2 TRA N SİYE N T ANALİZLERİNDE U Y G U LAN IŞI

H.İbrahim A R IK A N , U lvi ADALIO ÛLU

NÜKLEER M ÜH EN DİSLİK BÖLÜMÜ

Eylül 1985

TÜRKİYE ATOM ENERJİSİ KURUMU ÇEKMECE NÜKLEER ARAŞTIRMA VE EĞİTİM MERKEZİ

Teknik Rapor No: 27

AIREK—MOD KODUNUN YENİ VERSİYONU VE TR—2 TRANSİYENT ANALİZLERİNDE UYGULANIŞI

H.lbrahim ARIK AN, Ulvi ADALIOĞLU

NÜKLEER MÜHENDİSLİK BÖLÜMÜ

Eylül 1985

ÖZET

AIREK-MOD KODUNUN YENİ VERSİYONU VE TR-2 TRANSİYENT ANALİZLERİNDE UYGULANIŞI

AIREK-MOD kodu,araştırma reaktörlerinin geri beslemeli kinetiğini ele almaktadır.

Termodinamik parametrelerin sıcaklıkla değişimi kod'a altprog- ramlar halinde dahil edilerek TR-2 ve TRIGA araştırma reaktörlerine uygulaması yapıİmiştır.Rapor kodun yeni versiyonumu tanıtmaktadır.

SUMMARY

TRANSIENT ANALYSIS FOR THE TR-2 REACTOR WITH THE NEW VERSION OF AIREK-MOD

The code AIREK-MOD investigates point kinetics of research reactors,in the presence of temperature feedback.

Two sbroutines which consider the change of thermodynamic parameters with temperature were added to the code and the new code is applied to the TR-2 and TRIGA research reactors,fueled with 93$ and 20$ enriched uranium respectively.This report describes the new version of the code.

İÇİNDEKİLER Sayfa 1. GİRİŞ 1 2. HESAP METODU 1 2.1 AIREK-MOD kodu 1 3. BOZUNMA ISISI 3 4. CONKA PROGRAMI 3 4.1. Viskozite 3 4.2. Yoğunluk 4

4.3* Hacımsal genleşme katsayısı 4

4.4. Kondaktivite 5 4.5. Özgül ısı 5 5. TAKON PROGRAMI 7 6. ÖRNEK P R O B L E M : 1,2 8 REFERANSLAR 15 TABLOLAR

Tablo:1- TR-2 ile ilgili büyüklükler 11 ŞEKİLLER

Şekil:1- TR-2 reaktörü güç transiyentleri 12 Şekil:2- TR-2 reaktörünün tabii konveksiyon klapesiyle

yapılan deneyde ölçülen yakıt zarfı sıcaklığı­

nın değişimi 13

Şekil:3- TR-2 reaktörünün tabii konveksiyon klapesiyle yapılan deney için hesap edilen yakıt ve soğu­

Sayfa Ek 1 Ek 2 EKLER - Giriş datası - Program listesi 16 21

1. GİRİŞ

Daha önce yapılan bir çalışmada(l) araştırma reaktörlerinin geri beslemeli kinetiği AIREK—MOD koduyla(2,3) incelenmiştir. Kod topak parametre modeliyle sabit termodinamik sabitleri kul­ lanarak hesapları yapmaktadır.Daha sonra termodinamik parametre­ lerin sıcaklıkla değişimi kod'a altprogramlarla dahil edilmiş, TR-2 ve TRIGA araştırma reaktörlerine tatbikatı yapılmıştır(4).

B u rapor ise kodda yapılan değişiklikleri ve kodun kullanılı­ şını izah etmektedir.

2. HESAP METODU

2 . 1 AIREK-MOD kodu

AIREK-MOD kodu nokta kinetik denklemlerinden

(D

le gecikmiş nötron konsantrasyon denklemleri

< 2 ^ 1 = X;N(t )_\-VC£(t) U İ Û (₂)

ve soğuma denklemlerinin birleştirilmesi ile elde edilen

o

İt

2 K-6; k = n

> (3 )

denklem takımını çözmektedir.

Gözönüne a l m a n soğutma denklemlerinde,reaktör kalbi içinde bir tek yakıt çubuğu veya plakası etraflarındaki soğutucu kanal ile soğutmanın fiziki temsili yapılmıştır.Yakıt plakası (veya çubuğu) için yerden bağımsız,zamana bağlı bir T^(t) sıcaklığı,kanaldaki soğutucu için ise aynı şekilde bir T (t) sıcaklığı izafe edilmiştir.

Termodinamiğin birinci kanunu yakıt plakası (çubuğu) için dTf (t) dt = hAf (t)+-j£2iN(t)+ hA. ¥ fcf f V T 0 (t) soğutucu için ise

(4) dTc (t) dt ‘ hA 2v f c c| fV c c c c c c co -f hA 2v A f T^(t)-»— ^ . ( t ) denklemlerini vermektedir.Soğutucu sıcaklığı için

(5)

T

c(t)-f (t )+Tc (t ) (

6

)

bağıntısı kullanılmıştır.Burada

h Konveksiyon katsayısı (J/m2 .San.°C)

Af Yakıt plakası soğuma yüzeyi (m2 )

V f Yakıt plakası hacmi (m3 )

A

c Soğuma kanal alanı (m )

V

c Yakıt kanalında akışkan hızı (m/San.)

Tf Kararlı halde yakıt sıcaklığı

(°c)

T

c Kararlı halde akışkan sıcaklığı (°C) T.

1 Kararlı halde kanal giriş sıcaklığı

(°c)

Yakıt plakasının yoğunluğu (kg/m3 )

»e o Akışkanın yoğunluğu (kg/m3 )

°f s Yakıtın özğül ısısı (j/kg.°C)

c

Pisyon tesir kesiti(bir plaka üzerinde) (1/m)

Ortalama termal nötron hızı (m/San.)

3. BOZUNMA ISISI

Fisyon ürünlerinin bozunumundan ortaya ç ı k a n / v e

Y

enerjisi Way-Wigner formülü ile verilmektedir(5).Fisyon olayından t zaman sonra yayılan toplam radyasyon enerjisi watt cinsinden

(7) dır.

Burada P watt olarak reaktör gücünü.T sabit P gücünde gün olarak

O 0 0

çalışma zamanını ^TgUn olarak reaktör durdurulduktan sonraki herhangi bir zamandır.Reaktör seram yapıldıktan sonra bozunma ısısının

reaktörden atılmasının hesabında seramden hemen sonraki bozunma ısısı sabit ısı kaynağı değeri olarak kullanılmaktadır.Zamanla azalan ısı kaynağı ifadesinin kullanılması yerine bu sabit değerim kullanılması hesapları emniyetli tarafa bırakmaktadır.Hesap edilen toplam ısı kaynağı 173 kw değerindedir.Program modelinin özelliğin­ den dolayı fiktif bir fisyon enerjisi kabul edilerek bozunma ısısı

—15

kaynağı programa dahil edilmiştir.Fisyon enerjisi 0.2653x10 kcal/fis dır.

4, CONKA PROGRAMI

Suyun fiziksel özelliklerinin sıcaklıkla değişimi(6) ile ısı transferi sabitlerinin hesabı CONKA altprogramı içinde yapılmaktadır.

4.1. Viskozite(kg.nT1 .San""1 .)

Kanaldaki suyun sıcaklığı T ;0°c (T ^20°C arasında ise,viskozite

c c

1 3 0 1 iogıoM* 998.333+8.1855(Tc-20)+0.00585(Tc- 2 0 )

2

- 4 . 3 0 2 3 3

(

8

) 20°C^Tc^L00°C arasında ise logl0 r H o 1.3272(20-Tc )-0.001053(Tc- 2 0 ) Tc+105 2 i (9) ifadeleri kullanılmıştır. 4.2. Yoğunluk(kg.m )

Sıcaklık 0°C<Tc<17°C arasında ise yoğunluğun

?=1000-0.007 6(Tc-4)2 +0.49xlO- 5 (Tc-4)4 (10)

şeklinde değiştiği verilmektedir.Buna mukabil 17oC<Tc<,102°C ise

f

ıoc?

(T„-3.982)2(Tn+273)(350-T„)

( u )

466.7(Tc+67)(365-Tc ) şeklinde alınmıştır.

4.3. Hacımsal genleşme katsayısı(°C~^)

Sıcaklık 0°C<Tc<17°C iken hacımsal genleşme katsayısı

/=10“ 6 [o.26+15.94625(Tc-4)-0.243125(Tc-4)2+ 3 .3593xlO“ 3 (t c- 4 )^ (12)

şeklinde bir değişim gösterir. 17°C<Tc<100°C arasında ise

A

t df

T d % A ' - i k i [g(*c+273)(350-Tc )+(Tc-3.982)(77-2Tc ) (Tc- 3 .982) (Tc t-273) (350-Tc ) ( 298-2t g )1 (1 3) (Tc+67)(365-Tc )

3

şeklinde yazılır. 4.4. Kondaktivi t e (k a l .m ”^ .San”1 .°C”1 )

Kanaldaki suyun sıcaklığı 20°C<Tc<(100oC ise, suyun kondaktivitesi

K = 0 .4856+1.5757xlO_

3Tc- 5 .71429xlO-6Tc

(14)

ile verilmektedir.

4.5. Özgül ısıCkcal.kg”1 .0^ 1 )

Su sıcaklığı 0°C^Tc^35°C ise,özgül ısının eğrisi

cp= -0.2117xl0“6 (Tc-35)

3

+

0

.291xl0

"6

(Tc- 35)

2

+

0

.9979!) (

1 5

)

şeklinde bir polinom olmaktadır.

35°C'£Tcd00oC arasında ise,bu durumda

şeklindedir.

Yine CONKA altprogramı içinde kanaldaki ısı konveksiyon katsa­ yısı hesaplanmaktadır.Isı transferi sabitlerinden Nusselt sayısı

Nu=0.023Re° *8 .Pr° *4 R e y n o l d sayısı Re*-*p-ve Prandtl sayısı VrmSüL k

şeklinde alınmıştır.Burada eşdeğer (hidrolik) çaptır. Reaktörde yakıt tarafından üretilen ısı

«Ifis-Vf ^ | v thN(t) (17)

konveksiyon yolu ile yakıttan kanala geçen ısı

^con^k^f^Tf_ T c ^ ( 1 8 )

ve yakıt kanalından çekilen ısı

^ e k ^ c A c ^ p C T o - T i ) (19)

da kontrol için CONKA altprogramında hesaplanmaktadır.

Tabii konveksiyonda soğuma için Nusselt sayısı sabit 6.23 ve vc=0.06 m/San. alınmaktadır.

(₂₀) <lcon*( ^1+ *

2

) ( \ “ *2^,Ac D* c„-(To- T i ) S

İSİ t*

ile hesaplanabilmektedir(5)•

5. TAKON PROGRAMI

(4) ve (5) denklemleri soğutma denklemleriyle elde edilen (3) denklem takımının son üç denklemleriyle karşılaştırılırsa

Ym =*- » Ujh=0 , Pjjj** ^ m*lf•••fS Y8 - _{T v 7 + “T T “}A-ph . 2vcAc c*c , o8=o , p8-0 Ug=»0 y Pg**0 ve ayrıca (_{2 1}) (_{2 2}) (23) (24) Q^=0 _İsi

f

•• Q19s _fCfVf• lı _{Ql6*Qiys*Qiö*Q20*®}

Q24aAf.h_fcOoVc ^21*^22=^23 r* 2vf,An

«2 5- v ^ - ^2 6=Q27!S^28=<52 9®°

n

2Afh „

Q3 0 - - i r ^ Hk

şeklinde olur.Bu denklem katsayıları TAXON altprogramı içinde yazdırılarak hesaplanmaktadır.

6. ÖRNEK PROBLEM:1

TR-2 reaktörü 5 MW güçte kararlı olarak çalıştırıldığı zaman çeşitli reaktivitiler verilerek güç transiyentleri incelenmiştir. Örnek problemde Tablo:1 deki değerler fark denklemleri için kulla­ nılarak; j=0 . 003 lük pozitif bir reaktiviteye karşılık izafi nötron yoğunluğunun,yakıt ve soğutucu sıcaklıkları,reaktivitenin zamana

göre değişimi Şekil:1 de verilmektedir.

İzafi akı önce bir tepe değişimi gösterdikten sonra ilk değerinden daha büyük bir değerde kararlı olmaktadır.Yakıt ve moderatör sıcaklıkları devamlı yükselerek takriben 350 saniyeden sonra sabit kalmaktadır./«O.003 değeri için soğutucu kanalının çıkış sıcaklığı Tq ,kaynama noktasının altında olup 72°C civarın­ da sabitleşmektedir.

ÖRNEK PROBLEM:2

TR-2 reaktörü normal güçde çalışırken primer devre pompa

motoru arızalınıyor ve duruyor.Pompa volan sayesinde yavaşlıyarak dururken lüzumlu soğutmayı kısmen temin ediyor.Pompanın yavaşla­ ması esnasında debi nominal debinin #80 ine erişdiğinde reaktör seram yapıyor.Pompa motorunun arızalanmasından itibaren 180-200 Saniye içinde yani volan hemen hemen bütün enerjisini kaybettikten sonra reaktör taburesinin altındaki tabii dolaşım klapeleri açı­ lıyor ve kalbin soğutulması tabii konveksiyonla oluyor.

Bu olay esnasında yakıt zarf sıcaklığının ölçülen sıcaklığının zamanla değişimi Şekil:2 de verilmektedir.

Aşağıdaki kabullerle olayın matematik modeli kurulmuş ve hesaplar yapılmıştır.

i) Pompa motoru durduktan hemen sonra seram yapıldığı kabul edilmektedir.Esasında debi normal debinin #30'nine erişinceye kadar kalbe fisyon enerjisi doğmaktajseram olduktan sonra artık bozunma ısısı ortaya çıkmaktadır.Olayı basitleştirmek için olay başlangıcında seram olmakta ve bozunma ısısının doğduğu düşünül­ müştür.

ii) Pompa volanının enerji kaybı basit bir yaklaşımla ifade edilmektedir.200 Saniye içinde lineer bir debi aza?ması en basit bir yaklaşım olarak alınmıştır.

iii) 200 Saniye sonunda klapelerin açıldığı ve tabii konvek­ siyonun belli bir zaman içinde (mesalâ 20 saniye içinde) teşekkül ettiği kabul edilmiştir.

iv) Bu 20 saniye içinde soğutma kanalındaki akışkan hızı ile tabii konveksiyonun gene lineer bir değişim göstererek nominal değerlerine erişdikleri kabul edilmiştir.TR-2 reaktörü için nomi­ nal akışkan hızı 0.06 m/San. olarak verilmiştir.Nu sayısının tabii konveksiyondaki değeri

8.23

olarak alınmıştır.

v) Deneysel neticeler yaz şartlarında elde edilmiştir.Hesap­ larda elde mevcut kış şartları datası girilmiştir.Yani havuz suyu sıcaklığının ve başlangıç değer sıcaklıklarının kış şart­ larındaki değerleri alınmıştır.

AIREK-MOD'un verdiği neticeler Şekil:3 de özetlenmiştir.Model­ deki sezişe dayanan lineer değişimler deneysel neticelerin gidişi­ ni kalitatif olarak tekrar etirmektedir.Tabii konveksiyon için daha doğru korelasyonlar ile pompa yavaşlamasının doğru ifadesinin kullanılması hesap ile deneyin mutabakatını daha artıracaktır.

Olay başlangıcında bir miktar fisyon enerjisinin doğması ve soğutma debisinin azalması sebebiyle deneyde ortaya çıkan sıcaklık

yükselmesi,olay başlangıcında hemen seram kabulünden dolayı hesaplarda görülmemektedir.

Tablo:1- TR-2 ile ilgili büyüklükler

Parametreler Değerler

Yakıt plakası soğuma yüzeyi,

A f 797.226xl0”4 m 2 Yakıt plakası hacmi,

vf 50.624xl0"6 m 3 -4 2

Soğuma kanal alanı, A 1.3737x10 * m

Soğuma kanal hacmi,

TR-2 küçük kalbindeki toplam

c V

c 83.708x10“ 6 m 3

kanal sayısı, 298

Yakıt kanalında akışkan hızı, V

c 2.39 m/San. Kararlı halde kanal giriş sıcaklığı, T.

1 o

o

o

Kararlı halde kanal çıkış sıcaklığı, T 0

o

o

C

sJ

Kararlı halde akışkan sıcaklığı, T

c 36 °c

Kararlı halde yakıt sıcaklığı, T_r VJ1 H O O

Fisyon başına çıkan enerji,

i

0.765xlü“14 kcal/fis. Ortalama termal nötron hızı,

V th 2481.6 m/San. Ortalama makroskopik termal fisyon

tesir kesiti(bir plaka üzerinde).

t

_-P 30.87 1/m Ortalama termal nötron yoğunluğu,

X

N 1.351xl014 1/m3 TR-2 havuz suyunun kütlesi, M

c 45X104 kg. Soğutucu geri besleme katsayısı, -2xl0"4 1/°C

— I s s s s î e g i s s g g (\)°i — - g g s 3 s flite 0 ( 0 ^ 0 ^ 0 1 0 ^ 9 ^ 0 TT-*-- 'Î ( n n N N N r r O O O Ö Ö Ö Ö Ö Ö Ö O Ö Ö Ö c m in <v E r* İSİ Ş e k i l : 1 -T R -2 r e a k t ö r ü g ü ç t r a n s i y e n t l e r i

13 O 50 10 0 15 0 I 20 0

__

* t (s a n .) 2 5 0 Ş ek i l: 2 -T R -2 rea kt ör ü nü n t a b i i konv ek siy o n k l ap es i y le yapıl an d e ne y d e ö l çü l e n y a k ı t z a r f ı s ıc ak l ı ğ ı n ı n d e ğ i ş i m i ( R e a k t ö r i ş l e t m e b öl üm ünd en a l ı n m ı ş t ı r )

Z a m a n ( S a n .) Ş e k i l : 3 -T R -2 r e a k t ö r ü n ü n t a b i i k o n v e k s i y o n k l a p e s i y l e y a p ı l a n d e n e y i ç i n h e s a p e d i l e n y a k ı t v e s o ğ u t u c u s ı c a k l ı k l a r ı d e ğ i ş i m i

REFERANSLAR:

(1) H.İbrahim ARIKAN,Ulvi ADALIOĞLU,"Araştırma Reaktörlerinin Geri beslemeli Nokta Kinetik Analizi Üzerine Bir Çalışma" TAEK,ÇNAEM Teknik Rapor :18,1984.

(2) R.TUNCEL,U.ADALIOĞLU,"Reaktör Kinetik Denklemlerinin Nümerik Çözümü ve Transfer Fonksiyonu"ÇNAEM-NMB-30,1973.

(3) C.Mongini-Tomagnini,"AIREK-MOD:A Modified Version of the Kine tics Code AIREK ll",EUR.1914-e,1964.

(4) İ .ARIKAN,U.ADALIOGLU,"Araştırma Reaktörlerinin Geri Beslemeli Kinetiği"TAEK,11.Ulusal Nükleer Bilimler Kongeresi,

10-12 Ekim 1984.

(5) S.GLASSTONE,"Principles of Nuclear Reactor Engineering"1961. (6) J.MOTA,J.BRAZAO,"LFEN-N1.75 Notas Internas,Subprogramas

Fortran Para a Determmaçao de Algumas Proprıedades Fisicas da Agua e do Ar, "Sacaue'm,1972.

Ek 1- Giriş datası

Kart No

Format Değişkenler İzahat

1 2 6(16,E 1 2 .5)DATA(1) DATA(2) 6E16.Ö DATA(3) DATA(4) DATA(5) DATA(6) DA T A (7) DATA(ti; DATA(9) DATA(10)

Okunan kart numaraları(K İ ,K 2 ) Çözülecek toplam diferansiyel denklem sayısı(=10)

Gecikmiş nötron öncü konsantras­ yon denklem sayısı(=6>

=0.0 =0.0

*-1.0 RANAL dakı reaktıvıte açılımının birimi % dır. =+1.0 RANAL daki reaktıvıte açılımının birimi pcm dır.

=0.0

0.0 SUBROUTINE INHOR'a giriş yoktur.r(O) veya O) büyuKluic- lerı DATA(i>)'e bağlı olarak veri­ lir. Aynı zamanda W (O) değerierıde giriş datası olarak verilir.

W(Ü) SUBROUTINE INHOR'a giriş var­ dır. W(û) reaktör ters peryoau

verilerek W^(0) değerleri ve r(uj, j\0) büyüklükleride DATA(ö)'e bağ­ lı olarak INHOR dan hesaplanır. Başlangıçta verilen zaman ai'alıgı uzunluğu(=0.01)

Programda müsade edilen minimum mtegrasyon aralığı uzunluğu(=lu Programda müsade edilen maKsırııum

No Kart Format 4E16.8 5-6 6 7

8

6E16.8 3E16.8 3E15.7 Değişkenler D A T A (11) D A T A (12) D A T A (13) D A T A (14) D A T A (15) D A T A (16) D A T A (17) D A T A (18) D A T A (19) D A T A (20),D A T A (25) D A T A (26) D A T A (27) D A T A (28) D A T A (50) DATA(51) DATA(52) İzahat

integrasyon aralığı uzunluğu(10 ) — 2

^ 1 x 1 0 ) Bu iki data nötron yoğun­ luğunun eğimini test etmede kulla­ nılır.

=.1x10 ^ Bu test neticesinde integ­ rasyon intervali büyültülür veya küçültülür.

Gecikmiş nötronların toplam kesri(/) Nötron ömr ü ( l)

Sabit kaynak terimi(Sq )

=A Eğer/^^y//M(l) |^A ve r>B ise bir pike yaklaşılıyor demektir. Program otomatik olarak integ­ rasyon intervalini küçültür. =B Eğer |>A ise pik

geçilmiştir ve kontrol DATA(ll), DATA(12) ile yapılmaktadir.

-1 ise (3) tipinde geri besleme denklemleri kullanılmaktadır. t=s0 anındaki reaktörün nötron yoğunluğu(Nq )

( t ) = | ^ C (t)

Başlangıçtaki yakıt sıcaklığı(T„) " akışkan "

" kanal giriş " Okunan kart numaraları(K İ ,K 2 ) Kararlı halde yakıt sıcaklığı

" " akışkan " " " kanal giriş "

( T c )

Kart No

Format Değişkenler _İzahat

9 _{Okunan kart numaraları(K İ ,K 2 )}

10-11 _{DATA(69),DATA(93) Sıra ile gecikmiş nötron beta} kesri (f - 4 )

12 _{Okunan kart numaraları(Kİ,K2)}/

13 3E16.6 DATA(94),DATA(142) Sıra ile katsayıları.

15 Okunan kart numaraları(K İ ,K 2 )

16 _DATA(143) =0J0

16-17 3E16.8 DATA(144)»DATA(192) Sıra ile katsayıları.

18 _{Okunan kart numaraları(Kİ,K2)}

19-20 3E16.8 DAT A (194),DATA(242) Sıra ile katsayıları.

21 Okunan kart numaraları.

22 DATA(243) =0.ıD

DATA(244) =0.0

DATA(245) DATA(5)*e bağlı olarak r(0) veya f(0) değerlerini gösterir.

22-25 DATA(246),DATA(261) _&2»••«katsayılarıdır.

26 _{Okunan kart numaraları(K İ ,K 2 )}

27 DATA(263),DATA(264) a2 > ...,a sabitleridir.

28 _{Okunan kart numaraları(K İ ,K 2 )}

29 DAT A (311) Stop testinde görülen VV^ nin mertebe sayısı.

DATA(312) Stop testinde görülen VV nin tl mertebe sayısı.

D A T A (313) _=cı c^VV^+c^VV eşitsizliği DATA(314)

~°2,gerçekleştiği taktirde DATA(315) _{~C 3 program durmaktadır.}

DATA(316) Eğer N(t)^DATA(3l6) ise program durur.

_No__ 30

Kart Format Değişkenler

30-35 35 36 37 38 39 40 41 5E15.7 DATA(358) DATA(359) DATA(360) DATA(361) DATA(362) DATA(491) izahat

DATA(317) Eğer | <sÜp/N(t)|>DATA(317) ise program durur.

DATA(318) Eğer t^DATA(3l8) ise program durur.

DATA(319),DATA(348)(3) denklemindeki Q^(i=l,2,...,30) katsayıları.

DATA(349) =1 ise sıcaklığa bağlı katsayılar hesabı yapar.

Okunan kart numarası(Kİ,K2) DATA(350) Yer çekimi ivmesi(=f)

DATA(351) Yakıt plakası kalınlığı(=a) DATA(352) Yakıt plakası eni(=b)

DATA(353) Akışkan kanalının genişliği(= c ) DATA(354) Yakıt plakası yüksekliği(=H)

DATA(355) Akışkanın kanala giriş hızı(v.=0) DATA(356) Fisyondan çıkan enerji(=)^)

DATA(357) Ortalama makroskopik termal fis-yon tesir kesiti(=^J

Ortalama yakıt yoğunluğu(= f ) Yakıtın ısı kapasitesi(«c^)

Ortalama termal nötron hızı(=v ) th Reaktör havuzundaki su kütlesi(=M )

r c

Eğer DATA(362)=0 ise normal soğuma,

Eğer DATA(362)=1 ise bozunma ısı­ sının taşınması.

Okunan kart numarası(Kİ,K2) =0.0

Kart No

Format Değişkenler İzahat

41-42 DATA(492),DATA(497)Bu datalar program içinde sıcak­ lık adresleri olarak kullanılacak

42 DATA(498) =0.0

DATA(499) =0.0

E K 2 - P r o g r a m l i s t e s i

ANA PROGRAM

DIMENSION D A T A <500) »REAC<2) »ANS<20)»T V B <20)»TIME<2)»X <2>»0MEGA<3)» 1CAF’R<80) »ALPHN<20) » C AF’C ( 60 ) »DELTA (20)

COMMON Ef ATA » RE AC » J 1 »J 2 »J 3 »J 4 »ANS»T V B »T İME»X »OMEGA»SFW»CAPR»J R »CAPC 1 » ALF'HN» JC» J7» JN» IN2» J6» DELTA

WRITE<A»4)

4 FORMAT (1H1/////27H GIRISDATALARI I »DATA<I)///) 1 READ < 5 » 2 > K 1» K 2 READ <5 »13) < lı AT A ( I ) »I=K1»K2) WRITE < 6 »18)K 1» K2 W R I T E (6»3)(I»DATA<I ) »I=K1»K2) I F (K2-500) 1»10»10 10 J 1 = D A T A <1) J2=DATA<2)+1.0 J3=ABS < DAT A < 3)) J4=ABS<DATA<4)) JA=DATA<A) 14 IF<DATA<7)) 15»20»15 15 CALL INHOR 20 DO 25 I=1»J1 ANS <I)=DATA<1 + 18) 25 TVB<I )=ANS<I> ANX=ANS<1) INZ = 0 J7--0 JN =0 TIME <1)=0.O T I M E < 2 ) = O ♦ O X < 1 )=DATA<8) CALL RANAL R E A C d ) = RE AC < 2 ) SFW = 0 . O DO 45 1= 2» J2 45 3FW=SFW+DATA<I+A7)*ANS<I )

OMEGA <1)= D A T A (13)/DATA <14)*< (REAC <1)-1.O) + (SFW + DATA(15))/ANS < 1 ) ) WRITE<A»50) <EiATA< 1+92) »DATA< 1 + 142) » D A T A < 1 + 192) » I = 2 » J 1 )

50 FORMAT<İHO»43H FEEDBACK DENKLEMLERİNİN GIRIS KATSAYILARI /<3E1A.8

1 > )

W R I T E (A » 55) <DATA<I )»1 = 7»12)

55 F ORMAT<İHO»82H ILK.OMEGA ILK.INTERVAL IR'v'AL ALT LIMIT UST LIM I T / <AE1A . 8 ) >

WRITE < A » A O )< DAT A<I)»I = 13»15)» SFW

MIN.INTERVAL MAX.INTE AO F O RMAT<1H0»A1H BETA EFFECTIVE NEUTRON LIFETIME SOURCE TERM S

1UM FI*WI/<4 E 1 A .8)> AN3S=ANS<1)/ANX

WRITE<A» A5)TIME<1 ) » R E A C d )»O M E G A d )»ANSS

AS FORMAT <1H0»A3H TIME IN SEC. REACTIVITY AVG.INVERSE PER. NE 1UTR0N AKIŞI/<4E1A.8))

WRITE<a770) < ANS <I > »1 = 2»J İ )

70 F O R M A T d H O »52H GECİKMİŞ NEUTRON GURUPLARI FEEDBACK DEĞİŞKENLERİ l/(AE1A.8) )

13 F O R M A T (6 E 1 2 •5) 18 FO RM A T <2112) I F ( D A T A ( 3 4 9 ) . N E . 1 . ) G 0 TO 74 D A T A (100)= D A T A ( 3 5 6 ) * D A T A ( 3 5 7 ) * D A T A ( 3 6 0 ) / ( D A T A ( 3 5 8 ) » D A T A ( 3 5 9 ) ) D A T A ( 1 5 0 ) = - D A T A ( 100) I F ( D A T A ( 3 6 2 ) * N E * 1.) GO TO 31 D A T A ( 1 5 0 ) = D A T A ( 1 0 0 ) D A T A (100)=0♦ 31 C O N T I N U E D A T A ( 3 8 0 ) = 2 . » D A T A (353) D A T A ( 3 8 1 ) = 2 . » D A T A ( 352 )» D A T A (354) D A T A ( 3 8 2 ) = D A T A ( 3 5 1 ) » D A T A (352)» D A T A (354) D A T A ( 3 8 3 ) = D A T A ( 3 5 2 ) » D A T A (353) D A T A ( 3 8 4 ) = D A T A ( 3 5 2 ) » D A T A (353)» D A T A (354) 74 CALL INTER A N S S = A N S ( 1 > / A N X I F ( J N ) 9 9 » 2 5 1 » 9 9 251 W R I T E (A f 6 5 ) T I M E ( 1 ) »R E A C ( 1 ) »O M E G A ( 2 ) »ANSS W R I T E ( 6 » 7 0 ) ( A N S ( I ) »1 = 2 » Jl) W R I T E ( 6 » 7 5 ) A N S ( 8 ) »ANS(9) 75 F O R M A T ( İ H O »' YA KIT VE M O D E R A T O R S I C A K L I K L A R I '/ 2 E 1 6 ♦8) CALL TEST I F (J 7 ) 9 9 »7 4 » 99 99 C O N T I N U E STOP END S U B R O U T I N E CNFBT D I M E N S I O N D A T A ( 5 0 0 ) »R E A C ( 2 ) »A N S ( 2 0 ) »T V B ( 2 0 ) » T I M E ( 2 ) » X ( 2 ) » O M E G A ( 3 ) » 1C A P R (8 0 ) »A L P H N ( 2 0 ) » C A P C ( 6 0 ) » D E L T A (20) CO M M O N D A T A » R E A C » J l »J 2 » J 3 » J 4 » A N S » T V B » T I M E » X » O M E G A »S F W »C A P R »J R »CAPC 1 » A L P H N » J C » J 7 » J N » I N Z » J 6 » D E L T A JC = 1 DO 1170 I=1»J1 I F (A B S (A L P H N (I ) ) - l .0) 1 1 5 5 »1 1 6 5 »1160 1155 C D = 10.0 C N = 1.0 DO 1157 L = 1 » 7 C N = C N / C D » A L P H N ( I ) + l *0 1157 C D = C D - 1 »0 C A P C (J C + 2 )= C N / 3 ♦0 C A P C ( J C + 1 ) = ( C A P C ( J C + 2 ) * A L P H N ( I ) + 1 . 0) / 2 * 0 C A P C ( J C ) = C A P C (J C + 1 ) * A L P H N ( I ) + 1 » 0 J C = J C + 3 GO TO 1170 1160 IF(ABS(ALF‘HN (I )) -88 .0 ) 1165 » 1165 »11 61 1161 I F (A L P H N (I )) 1 1 6 2 »1 1 6 2 »1180 1162 C A P C ( J C ) = - 1 . 0 / A L P H N ( I ) GO TO 1166 1165 CAPC(JC) = ( E X P ( A L P H N ( I) ) - l ♦ 0 ) / A L P H N ( I ) 1166 C A P C (J C + 1 ) = ( C A P C ( J C ) - l «0)/ A L P H N (I ) C A P C ( J C + 2 ) = ( 2 * 0 » C A P C ( J C + 1 )- 1.0)/ A L P H N (I )

1170 CONT I NUE GO TO 1171 1180 J N = - 1 1171 RETURN END S U B R O U T I N E I N H O R DI MENSI ON D A T A ( 5 0 0 ) » R E A C < 2 ) » A NS < 2 0 ) » T V B < 2 0 ) » T I M E ( 2 ) » X < 2 ) » OMEGA< 3 ) f 1CAPR( 8 0 > » ALPHN < 2 0 ) » C A F C ( 6 0 ) » D E L T A ( 2 0 ) C O M M O N D A T A »R E A C » J 1 »J 2 »J 3 »J 4 »A N S »T V B »T İ M E »X »O M E G A »S F W »C A P R »J R »C A P C 1 . A L P H N . J C . J 7 * J N . I N Z . J A . D E L T A T E M P 1= 0♦ 0 DO 1000 1 = 2 »32 T E M P 2 = D A T A ( 7 ) - D A T A < 1+192) D A T A < 1 + 1 8 ) = - D A T A ( 1 + 1 9 2 ) » D A T A ( 1 9 ) / T E M P 2 1000 T E M P 1 = T E M P 1 + D A T A < I + A 7 ) / T E M P 2 R E A C ( 1 ) = ( T E M P I + ( D A T A ( 1 4 ) / D A T A ( 1 3 ) > ) » D A T A (7) I F ( D A T A ( 5 ) ) 1005 » 1 0 0 5» 1010 1005 R E A C ( 1 ) = R E A C ( 1 ) * D A T A ( 13) 1010 I F ( D A T A ( 5 ) - 1 .0) 1 0 1 5 » 1 0 1 5 » 1020 1015 D A T A ( 2 4 4 )= R E A C (1> 1020 D A T A ( 3 4 9 ) = R E A C (1) R E T U R N EN D S U B R O U T I N E R N F B T D I M E N S I O N D A T A ( 5 0 0 ) » R E A C ( 2 ) »A N S ( 2 0 ) »T V B ( 2 0 ) » T I M E ( 2 ) » X ( 2 ) » 0 M E G A ( 3 ) » 1C A P R ( 8 0 ) »A L P H N ( 2 0 ) » C A P C ( 6 0 ) » D E L T A (20) C O M M O N D A T A » R E A C » J 1 » J 2 » J 3 » J 4 » A N S » T V B » T I M E » X » 0 M E G A » S F W » C A P R » J R » C A P C 1 » A L P H N » J C » J 7 » J N » I N Z » J 6 » D E L T A S F W = 0 ♦0 DO 1 150 1 = 2 » J2 1150 S F U = S F U + D A T A ( I + A 7 ) * T V B ( I ) C A P R (J R ) = ( ( R E A C ( 2 ) - R E A C ( 1 ) ) * T V B ( 1 ) + S F W + D A T A ( 1 5 ) ) » D A T A ( 1 3 ) / D A T A (14) I F (D A T A (18)) 3 » 3 » 2 2 DO 4 I = 2 » J1 J R = J R + 4 C A P R ( J R ) = D A T A ( 1 + 9 2 ) » T V B ( 1 ) + D A T A ( 1 + 1 4 2 ) » D A T A (19) 4 C O N T I N U E G O T O 6 3 J 1 = J 1-6 DO 1153 1 = 2» J1 J R = J R + 4 1153 C A P R ( J R ) = D A T A ( 1 + 9 2 ) » T V B ( l ) + D A T A ( 1 + 1 4 2 ) * D A T A ( 19) S L B 1 = 0 . 0 DO 8 4 7 A 1 = 1 »5 J 1 A I = J 1 + I + 1 8 4 7 6 S L B 1 = S L B 1 + D A T A ( I + 3 1 8 ) » T V B ( J 1 A I ) D A T A ( 4 9 2 ) =SLB1 S L B 2 = D A T A ( 3 2 4 ) # T V B ( J l + 1 ) DO 8 4 7 7 1 = 1 »4

1170 C O N T I N U E GO TO 1171 1180 J N - - - 1 1171 RETURN END S U B R O U T I N E I N H O R D I M E N S I O N D A T A ( 5 0 0 ) » R E A C ( 2 > *A N S ( 2 0 ) » T V B ( 2 0 > » T I M E ( 2 ) » X < 2 ) » O M E G A (3) » 1C A P R < 80) » A L P H N ( 2 0 ) » C A P C ( 6 0 > » D E L T A (20) C O M M O N D A T A » R E A C » J 1 »J 2 » J 3 » J 4 » A N S » T V B » T I M E » X » O M E G A » S F W » C A P R » J R » C A P C 1 .A L P M N » J C .J 7 .J N »I N 7 . J 6 * D E L T A T E M P 1 = 0 » 0 DO 1 0 0 0 1 = 2» 32 T E M P 2 = D A T A ( 7 > - D A T A < 1+19 2 ) D A T A ( I + 1 8 ) = - D A T A ( 1 + 1 9 2 ) * D A T A ( 1 9 ) / T E M P 2 1 0 0 0 T E M P I = T E M P 1 + D A T A ( 1 + 6 7 ) / T E M P 2 R E A C ( 1 ) = ( T E M P 1 + ( D A T A ( 1 4 > / D A T A ( 1 3 ) ) > * D A T A ( 7 ) I F ( D A T A ( 5 ) ) 1 0 0 5 » 1 0 0 5 » 1010 1005 R E A C ( 1 ) = R E A C ( 1 > * D A T A ( 13> 1 0 1 0 I F ( D A T A ( 5 ) - 1 »0) 1 0 1 5 » 1 0 1 5 » 1 0 2 0 1 0 1 5 D A T A ( 2 4 4 ) = R E A C ( 1 > 1 0 2 0 D A T A ( 3 4 9 > = R E A C ( 1 > R E T U R N E N D S U BROU T I NE RNFBT D I M E N S I O N D A T A ( 5 0 0 ) » R E A C ( 2 ) »A N S ( 2 0 ) »T U B ( 2 0 ) »T I M E ( 2 ) » X ( 2 ) » O M E G A (3 ) » 1C A P R ( 8 0 ) »A L P H N ( 2 0 ) > C A P C ( 6 0 ) » D E L T A (20) C O M M O N D A T A » R E A C » J 1 »J 2 »J 3 »J 4 »A N S »T U B »T I M E »X »O M E G A »S F U »C A P R »J R »C A P C 1 » A L P H N » J C » J 7 » J N » I N Z » J 6 » D E L T A SF W = 0 . 0 DO 1150 I = 2 » J 2 1150 S F U = S F U + D A T A ( I + 6 7 ) * T U B ( I ) C A P R ( J R ) = ( ( R E A C ( 2 ) - R E A C ( 1 ) ) * T V B ( 1 ) +SFW + D A T A ( 1 5 > ) * D A T A ( 1 3 ) / D A T A ( 1 4 ) I F ( D A T A ( 1 8 ) ) 3 » 3 » 2 2 DO 4 I =2 » J 1 J R = J R + 4 C A P R ( J R ) = D A T A ( 1 + 9 2 ) * T U B ( 1 > + D A T A ( 1 + 1 4 2 ) » D A T A ( 1 9 ) 4 C O N T I N U E GOTO 6 3 J 1 = J 1 -- 6 DO 1153 1=2 » J 1 J R = J R + 4 1 153 C A F ' R ( J R ) = D A T A ( I + 9 2 ) » T V B ( 1 ) + D A T A ( I + 1 4 2 ) » D A T A ( 1 9 ) S L B 1 = 0 . 0 DO 8476 1 = 1 » 5 J 1 A I = J 1 +1+1 8476 S L B 1 = S L B 1 + D A T A ( 1 + 3 1 8 ) » T V B < J 1 A I ) D A T A ( 4 9 2 ) = S L B 1 SLB2 = D A T A ( 3 2 4 ) » T U B ( J l + 1 ) DO 3 477 1 = 1 »4

8481 8477 SLB 2=S LB2+DATA(I+324>*TVB<J1BI) DATA(493 )=S LB2 SLB3=0.0 DO 8478 1=1»2 J1CI=J1+I 8478 S L B 3 = S L B 3 + D AT A ( I + 3 2 8 ) * T V B < J 1 C I ) DO 8479 1 = 1,3 JlDI=Jl+I+3 8479 S L B 3 = S L B 3 + D A T A ( 1 + 3 3 0 >*TVB<J1DI) DATA(49 4)=SLB 3 SLB4 = 0 * O DO 8480 1 = 1,3 J1EI=J1+I 8480 S L B 4 = S L B 4 + D A T A ( I + 3 3 3 ) # T V B ( J 1 E I ) DO 8481 1 = 1,2 JlFI=Jl+I+4 SLB4=SLB4+DATA<1+336> * T U B ( J 1 F I > DAT A(495)=SLB4 S L B 5 = 0 .O DO 8482 1=1,4 J1GI=J1+I 8482 SL B5 =S LB 5 + D AT A ( I +3 38 )* T U B (J 1 G I ) S LB5 = S L B5 + D A T A ( 3 4 3 ) * T V B< J l + 6 ) DATA(496)=SLB5 SLB6=0 * O DO 8483 1=1,5 J1HI=J1+I SLB 6=S LB6+DATA<I+343)#TVB(J1HI> DATA(497)=SLB6 J 1 T = J 1 +1 J1 = J 1 + 6 DO 8484 I=J1T,J1 JR=JR+4 LLL=I+13~Jl+484 C A P R ( J R )=DATA <1+92)# T V B < 1> + D A T A < 1 + 142)* D A T A < 19)+DATA < L L L ) RETURN END 8483 8484 6 SUBROUTINE RANAL

DIMENSION D A T A (500),R E A C <2 ),A N S (20),T V B <20),T İ M E <2)>X <2),O M E G A <3> , İCAF’R ( 80 > , ALF'HN ( 20) , CAPO < 60 ), DELTA ( 20 )

COMMON D A T A , R E A C , J 1 ,J 2 ,J 3 ,J 4 ,A N S ,TVB»T I M E , X , O M E G A ,SFU,CAPR,JR,CAPC 1 yALPHN,JC,J 7 » J N , INZ»J6,DELTA

REA C(2 >=D ATA<24 5)+D AT A<2 46 >+DATA<247 )*TIME<2>+DATA<248)*(TIM E<2)* * 12>+DAT A<2 49>* <TIME (2> **3>+ nATA<250>*SIN(DATA(259>*TIME<2> )+DATA(25 2 1 ) * S I N < D A T A < 25 9> *TI ME <2> *2 .)

RE AC (2 )= REA C( 2)+ DA T A < 2 5 2 ) # S I N ( D A T A (2 5 9 > * T I M E ( 2 ) * 3 . >+ D A T A <253>* C O S ( 3 D A T A ( 2 6 0 ) * T I M E < 2 > >+ D A T A <254>* C O S <D A T A <260>* T I M E (2)*2♦)+ D A T A (255)*C 4 0 S ( D A T A ( 2 6 0> * T I M E < 2 > * 3 . )

R E A C ( 2 > = R E A C( 2 > + D A T A < 2 5 6 ) * E X P < D A T A < 2 6 1 >*TIME<2>)+ D A T A <257)* E X P (DAT 1 A < 2 6 1 ) * T I M E( 2 > * 2. 0)+ DA TA< 25 8)*EXP<D ATA(261>*TI M E <2)*3.0>

J J J J = J 2 J J J J DO 5 K = 1 »J J J M I = K + J J J J 5 S U M= S U M+ D A T A ( K + 2 6 2 ) * T V B ( MI > R E A C ( 2 ) = R E A C ( 2 ) + S U M I F ( DAT A C 5 ) ) 1065 » 1 0 7 0 » 1 0 7 0 1065 R E A C ( 2 ) = R E A C ( 2 ) / D A T A ( 1 3 ) 1070 RETURN END S U B R O U T I N E TEST D I M E N S I O N D A T A (500> »R E A C (2> *A N S < 2 0 ) »T U B (20> »T I M E ( 2 ) » X ( 2 ) » O M E G A (3) » 1 C A P R ( 8 0 ) » A L P H N ( 2 0 ) » C A P C ( 6 0 ) » D E L T A ( 2 0 ) C O M M O N D A T A » R E A C » J 1 »J 2 » J 3 » J 4 » A N S » T V B » T I M E » X » O M E G A » S F U » C A P R » J R » C A P C 1 » A L P H N » J C » J 7 » J N » I N Z » J 6 » DELTA I F < A N S < 1 > - D A T A ( 3 1 6 ) ) 3»3»4 3 IF ( T I ME ( 1 > --DATA (318) ) 5»5»4 5 IF < O M E G A < 1 > - D A T A ( 3 1 7 ) ) 6»6»4 6 M A = D A T A (311) M B = D A T A (312) I F ( A N S < M A ) » D A T A ( 3 1 3 ) + A N S ( M B ) * D A T A < 3 1 4 ) - D A T A (315)) 7»7»4 4 J7= 1 7 R E T U R N END SUBROUTI NE INTER-D I M E N S I O N INTER-D A T A ( 5 0 0 ) » R E A C ( 2 ) »A N S ( 2 0 ) »T U B ( 2 0 ) » T I M E ( 2 ) » X ( 2 ) » O M E G A ( 3 ) » 1 C A P R ( 8 0 ) » A L P H N ( 2 0 ) »C A P C ( 6 0 ) » D E L T A (20) C O M M O N D A T A » R E A C » J 1 »J 2 » J 3 » J 4 » A N S » T V B » T I M E » X » O M E G A » S F W » C A P R »J R » CAPC 1 » A L P H N »J C » J 7 » J N » I N Z » J 6 » DELTA I F ( D A T A ( 3 4 9 ) .£0*1.) CALL TACON 75 J R = 1 CALL RANAL CALL RNFBT X ( 2 ) = X ( 1 ) /2 » 0 A L P H N ( 1 > = ( R E A C ( 1 ) - l .0)* X (2)* D A T A (13)/ D A T A (14) DO 95 I ~ 2 » J1 95 A L P H N ( I ) = X ( 2 ) * D A T A ( 1+192) CALL CNFBT I F ( J N ) 96 » 97 » 96 96 CALL DI REC J F ( J N ) 256 » 75 » 256 97 J R = 1 J T = 0 JC = 1 DO 105 I = 1 » J 1

B ELT A ( I ) = ((A L P H N (I )* A N S (I ) ) + C C A P R (J R )* X ( 2 ) ) )*CAPC(JC)

J R=J R+4 J C = J C + 3

I F ( D E L T A ( I ) - ( 2 .0 * * ( - 1 4 . ) )*ANS(I)) 10 5 » 10 0 » 1 0 0

JR = 2 T I M E ( 2 ) = T I M E ( 1 ) + X ( 2 ) CALL RANAL CALL RNFBT JR = 1 JC = 1 DO 125 1 = 1 , J1

D E L T A ( I ) =DELT A< I > + ( X< 2 > * CA PC< J C+ 1 > * < C APR< J R+ 1 ) - CAPR< J R) ) ) JR=JR+4 J C=JC+3 125 T V 8 ( I ) = A N S ( I > + D E L T A < I > JR = 3 CALL RANAL CALL RNFBT DO 150 1 = 1 , J1 150 ALF’HN ( I ) =X (1 > #ALPHN ( I > / X ( 2 ) CALL CNFBT I F ( J N ) 151>152, 151 151 CALL DIREC I F ( J N ) 256 »75 , 256 152 JR = 1 JC=1 DO 160 1 = 1 , J1 D E L T A ( I ) = 2 . 0 * X ( 1 ) * C A P C < J C + 1 > * < C A P R < J R + 2 ) - C A P R < J R ) ) + C A P C ( J C ) * ( X ( 1 >* 1CAPR<JR) +ALPHN<I ) *ANS<I > >

T O B ( I ) = D E L T A ( I ) + A N S ( I ) JC=JC+3 160 JR=JR+4 T I M E < 2 ) = T I M E ( 1 ) + X < 1 ) JR = 4 CALL RANAL CALL RNFBT JR = 1 JC=1 DO 168 1 = 1 , J1 D E L T A ( I ) = D E L T A ( I > + X ( 1 >*<CAPR<JR+3) - 2 . 0*CAPR<J R+ 2 ) + C A P R < J R ) > * < 2 . 0 * C 1APC( JC + 2 ) - C A F ‘C< JC + 1 ) ) JR=JR+4 J C=JC+3 168 T O B ( I ) = A N S ( I > + D E L T A ( I > CALL RANAL CALL RNFBT

0MEGA( 3) =DAT A<1 3 ) / DAT A<14>*< < R E A C < 2 ) - 1 . 0 ) + < SFW+DATA( 1 5 ) >/T VB<1 ) ) I F ( T O B ( 1 ) / A N S ( 1 ) ) 1 8 0 , 1 8 0 , 1 9 0 180 JN=1 CALL DIREC GO TO 256 190 0 ME G A ( 2 ) = ( A L 0 G ( T V B < 1 > / A N S ( 1 > ) ) / X ( l ) I F ( J T ) 200 , 250 »200

200 Q = X ( 1 >*AB3<0MEGA<3>- 2. 0*0MEGA( 2>+0MEGA( 1) > *CAPC<2>/<CAPC< 1 > + 1 . 0 ) I F ( D A T A ( 1 6 ) ) 2 0 5 , 2 3 0 , 2 0 5

205 I F ( D A T A ( 1 6 ) - 0 M E G A ( 3 ) ) 2 3 0 , 2 1 0 , 2 1 0

210 I F ( A B S ( D A T A ( 1 5 ) / T V B < 1 ) + R E A C ( 2 ) ) - D A T A < 1 7 ) ) 230*230, 215 215 I F ( I N Z ) 2 2 5 , 2 2 0 , 2 2 5

X ( 2 >=DA T A( 9 ) GO TO 245

225 I F ( ABS ( OMEGA ( 3 > > -DATA-< 16 > ) 2 5 0 , 2 3 0 . 2 3 0 230 I F ( Q - D A T A ( 1 1 ) ) 2 3 5 , 2 5 0 , 2 4 5 235 I F ( Q - D A T A (12> ) 2 4 0 , 2 5 0 , 2 5 0 240 X (1 )=2. 0*X< 1 ) GO TO 250 245 CALL DIREC I F ( J N ) 2 5 6 , 7 5 , 2 5 6 250 CONTINUE I F ( X ( 1 ) . LT. 1. ) X(1 )=1. IF* X < 1 > . GT. 5. ) X < 1 ) =5. IF< TIME< 2) . GT. 200. > X ( l ) = l. I F < TI ME < 2 >. GT. 240. ) X<1)=5. TIME(1)=TIME(2 > OMEGA( 1 i - OMEGA<3) REAC( 1 >=REAC( 2 ) DO 255 1 = 1 , J1 255 ANS( I )=TV3< I ) 256 RETURN END S U B R O U T I N E DIREC

D I M E N S I O N D A T A ( 5 0 0 ) , R E A C ( 2 > , ANS<20), TVB<20>, TIME(2), X < 2 ) ,0MEGA<3>, 1CA P R (8 0) .A L P H N (20), CA PC (60), D E L T A (20)

C O M M O N DATA, REAC, Jl, J2, J3, J4, ANS, TVB, TIME, X, OMEGA. SFW, CAPR, JR, CAPC 1, ALPHN, JC, J7, JN, INZ, J6. D E L T A IF(JN) 1 1 8 4 , 1 1 8 4 , 2 0 3 5 1184 I F (X (2)- D A T A (9 >) 1 1 8 5 , 1 1 9 0 , 1 1 9 0 1185 IF < X ( 1 ) - D A T A (9) ) 2 0 0 0 , 2 0 0 0 , 1 1 9 0 1190 X < 1 )=X(2) T I M E (2)= T I M E ( 1> R E A C ( 2 ) =REAC <1) DO 1195 1= 1 , Jl 1195 T V B (I ) = A N S ( I > J N = 0 G O TO 2 0 3 5 2 0 0 0 IF(JN) 2 0 3 5 , 2 0 1 5 , 2 0 3 5 2 0 1 5 JN=1 2 0 3 5 R E T U R N END SUBROUTINE TACON

DIMENSION DA T A( 5 0 0 ) , REAC( 2 ) , A N S ( 2 0 ) , T V B ( 2 0 ) , TI M£<2>, X<2>,OMEGA*.3), 1CAPR( 8 0 ) , ALPHN( 2 0 ) , CAPC( 6 0 ) , D E L T A ( 2 0 )

COMMON DATA, REAC, J l , J2, J3, J4, ANS, TVB, TIME, X, OMEGA, SFW, CAPR, JP, CAPC 1, ALPHN, JC- J7, JN, INZ, J6, DELTA

T F = T V B ( J l - 2 ) T C = T V B < J l - 1 ) T I = T V B ( J l ) TG=2. * T C - T I

W R I T E (6* 1 0 ) T F, TC.- TO, T I

10 FORMAT< /2X, ' T - Y A K I T » E15. 7, 2X, ' T - S Q G U T * ' , E l 5. 7, 2X, ' T - C I K = ', E15. 7, 12X, ' T - G I R = ' , E l 5. 7)

CALL CQNKA( TF, TC, T I , TO, DAT A(381 ) , DAT A(383>, DATA<380>, VC, RO, VI SC, CP 1 , HCONV)

AFCON=DATA<381> »HCONV ROCP=RO*CP*DATA( 3 8 4 )

D A T A (200)— A F C O N / (D A T A < 358 > » D A T A (359 ) » D A T A (382))

D A T A < 2 0 1 > - - < AFCON/RQCP+2. *VC»DATA<3 8 3 ) / DAT A( 3 8 4 ) ) D A T A <202>= - 2 »AFC0N*298. / ( D A T A <361) *CP) I F ( DATA< 362) . NE. 1. ) GO TO 15 D A T A <202) =0. 15 D A T A ( 3 3 7 / - - D A T A ( 2 0 0 ) DATA < 3 4 2 ) =AFCON/ROCP DA T A ( 3 4 3 ) = 2 . * VC » D A T A < 3 8 3 ) /DATA < 384) DATA < 3 4 8 ) » - D A T A ( 202 > RETURN END

SUBROUTINE CONKA(TF, TC, T I , TO, AF, AC, D, VC, RO, VI SC, CP, HCONV)

DIMENSION DATA<5 0 0 ) , REAC( 2 ) , ANS<20) , T VB<20) , TI ME<2>, X<2>, OMEGA(3>> 1CAPR< 80 > , ALPHN( 2 0 ) , CAPC < 6 0 ) , DELTA < 20)

COMMON D A T A , R E A C , J l , J 2 , J 3 . J 4 » A N S , T V B , T I M E , X , OMEGA, SFW, CAPR, JR, CAPC 1, ALPHN, JC, J 7 , JN, INZ, J 6 , DELTA

I F C T C - 2 0 . ) 1 , 2 , 2 1 AVI SC = 1301. / ( 9 9 8 . 333+8. 1855*<TC- 20. > +0. 00585* < T C - 2 0 . ) * * 2 ) - 3 . 30233 VI SC=10. **<AV I S C — 1. > GO TO 5 2 AVISC=<1. 3272*<20. - T O - 0 . 001 0 5 3 * ( T C - 2 0 . > * * 2 ) / ( T C + 105. ) VI SC=. 01002*10. » » ( A V I S C - l . ) WRITE <6, 3 ) A V I S C , VI SC

3 FORMAT(2X, ' A V I S C = ' , E15. 7, 2X, ' V I S C = ' , E15. 7) 5 I F ( T C —17. ) 10, 20, 20

10 R0=1000. - . 0 0 7 7 6 * ( T C —4. )**2+. 4 9 E - 5 * ( T C - 4 . )**4

BETA=. 26+15. 9 4 6 3 * ( T C - 4 . ) - . 24 3 1 3 * ( T C - 4 . )**2+. 00 3 3 6 * ( T C - 4 . >**3 BET A=BET A*1. E - 6

GO TO 30 20 R0=1000. - ( T C - 3 . 982) **2* ( TC+273. ) * ( 3 5 0 . - T C ) / ( 4 6 6 . 7*( T C+67. ) * ( 3 6 5 . ~T 1C ) ) DR0=2. * < TC+273. ) * ( 3 5 0 . - T C ) + ( T C - 3 . 982>*( 77. - 2 . *TC ) DRO=DRO— ( T C - 3 . 982 )*< TC+273. >*<350. - T C > * ( 2 9 8 . - 2 . * T C > / ( ( T C + 6 7 . >*( 365 1. - T O ) DRQ= ( T C - 3 . 982 >*DRQ/ (466. 7 * ( TC+67. >*(365. - T O > 3E i A=DRO/RO 30 GRQN=9.' 80665*BETA*R0**2*D**3* ( T F —TC ) /VI SC**2 CCND=0. 4856+1. ' 3757E- 3*T C- 5. 71 429E- 6*T C**2 C 0 N D=C 0 N D / 3 6 Ö 0. 19 ( T Ç —35. > 50. aG, 60 50 CP=- 0. 2117E- 6*! . T C - 35. >**3+0. 2 9 1 E - 6 * ( T C - 3 5 . >**2+0. 99795 GO TO 70 60 CP=2. 1 3 4 9 E - 6 * ( T C - 3 5 . )**2+. 99795 70 PRN=CP*VISC/COND PA=-GRON»PPN

I F ( T Î . LT. 17. ) GO TO 75 R 1 =1000. — ( T I - 3 . 982>**2*<TH- 273. >*<350. - T I ) / < 4 6 6 . 7*<TI-*-67. >*<365. - T I I ) ) GO TO 80 75 Rî =1000. —. 0 0 7 7 6 * ( T l —4. >**2+. 4 9 E - 5 * C T I - 4 . >**4 80 CONTINUE I F ( TO. LT. 17. > GO TO 85 R2=1000. - C T 0 - 3 . 982>* * 2 * <TO+273. >*<350. - T 0 > / < 4 6 6 . 7*<T0+67. >»<365. - T 10) > GO TO 90 85 R2=1000. - . 0 0 7 7 6 * ( T O - 4 . >**2+. 4 9 E - 5 * < T O - 4 . ) **4 90 CONTINUE I F ( D A T A ( 362 >-1. ) 1 1 , 1 2 , 1 1 12 CONTINUE I F ( T I M E ( 2 > . GT. 200. > GO TO 114 QD=<4 4 0 - 2 2 * T I M E ( 2 ) > * 0 . 80 W R I T E <6,121> QD 121 FORMAT(10X, ' Q D = ' , E 1 5 . 7> VC=QD/<298. » D A T A <383>*3600. > GO TO 13 114 I F i TIME<2> . GT. 220. > GO TO 115 VC=0. 003*< T I M E <2>—200. > C V C = ( R1—R 2>» D A T A<3 5 0)* D* * 2 /(96. *VISC> I F C V C . GT . 0 , 0 6 ) VC=0. 06 GO TO 13 115 VC=Q. 06 GO TO 13 11 VC=2. 39 13 CONTINUE

I F (DATA <362) EG. 1. .AND. T I ME( 2> . GT. 2 0 0 . ) GO TO 34 R E=VC *RO*D/ VISC

GR=RE*PRN*DATA( 3 8 0 ) /DATA < 354 > I F (RE. LE. 2000 > GO TO 33 UNU= 0 2 3 * R E » » 8 » P R N * * . 4 GO TO 35 33 UNU=2. *GR**0. 333 GO TO 35 34 CONTINUE UNU=8. 23 IF T I M E ( 2 > . LT. 2 2 0 . ) UNU=0 2 4 5 * ( T I M E ( 2 > - 2 0 0 >+3. 32 35 HCONV=UNU*COND/D AR1 2 = ( R1+R2) » < R1- R2 > AVE=192 *VISC CPI 2=CP* < TO—T I > QTAS=VC*DATA( 383 > *RO*CP12

QTTT=AR 12*DATA < 383) * < D**2 > *CP 12'*DATA < 350 > /AVE

OF IS--DATA ( 382 > » DATA < 356 > »DATA < 357 > »DATA < 360 > » DATA ( 19 > QCONv»HCONV»AF*( T F - T C >

WR I T E ( 6, 14 > VI S C , RO, BETA, COND, CP

14 FORMAT</2X, ' VI SC= ' , E l 5. 7, 2X\ 'R0= ', E l 5. 7, 2X. /BETA= ' , E15. 7, 2X, 'COND= 1 E15. 7, 2X, 'CP= ', E15. 7)

WRITE!( 6, 1 5 ) GFON, PRN, RA, UNU

15 FORMAT< /2X, "GRAS= ', E 15 7, 2X, 'PRANDL = ', E15. 7, 2X, ' R A Y L E I G = ' , E15. 7, 2X 1, ' N U S S E L T - ' , E15. 7/>

WRI T E(6, 1 6 ) QCONV,OTAS, GFI S, QT T T

16 FORMATCIOX, 'GCONV» E15. 7, 10X, ' G T A S I M A » E15. 7, 10X, ' G F I S Y » '. El 110X, 'QTAB= ', E l 5. 7)

WRI T E(6, 1 7 ) V C , R E - GR» HCONV

17 FORMAT < /1 OX» HI Z - ', El 5. 7, 5X, ' REYNOLD»' » E15. 7, 5X, ' G R A T Z » ', E15. 7 1 'HCONV» £15. 7/ )

RETURN END