Bandung,
Proceedings Seminar Reaktor Nuklir dalum Penelitian Sains dan Teknologi Menuju Era Tinggal Landas
8 -10
Oktober 1991 PPTN - BATAN
SIMULASI REAKTOR SUR-tOOBERDAYANOLDENGAN MENGGUNAKAN PROGRAM TUTSIM DAN FANSIM Sudjatmi K.A, Dudung A.R. Pusat Penelitian Teknik Nuklir - Badan Tenaga Atom Nasional ABSTRAK SIMULASI REAKTOR SUR-IOO BERDAYANOL DENGAN MENGGUNAKAN PROGRAMTUTSIM DAN FANSIM. Karakteristik dinamis dari reaktor SUR-IOOyang beroperasi pada daya rendah telah disimulasikan. Pada kondisi ini, pengaruh temperatur dan xenon terhada p reaktivitas dapat diabaikan. Model matematis sistem reaktor dimasukkan ke dalam program TUTSIM, kemudian data dalam daerah frekuensi dianalisis dengan program FANSIM. ABSTRACT THE SIMULATION OF ZERO-POWERED-SUR-IOOREACTOR WITH TUTSIM AND FANSIM PROGRAM.The dynamic characteristics ofSUR-lOO reactor operating in low power has been simulated. In this condition, the temperature and xenon effect of reactivity is negligible. The mathematical model of reactor system fed in the TUTSIM program. Then the data is analyzed in frequency domain with the FANSIM program. PENDAHULUAN
Simulasi adalah sarana untuk penelitian hubungan besaran fIsis suatu sistem dalam suatu waktu tertentu. Bila model sistem telah tersusun dalam persamaan matematis maka komputer dapat dipergunakan untuk mengamati tingkah laku sistem tersebut tanpa harus membuat lebih dulu sistemnya. Keuntungan simulasi dengan mempergunakan komputer adalah : - Penampilan sistem dapat diamati padasemua keadaan. - Hasil penampilan sistem lapangan dapat diekstrapolasikan dengan membuat peniruan dengan model untuk keperluan peramalan. - Keputusan - keputusan yang berkaitan dengan masa depan pada tingkat konseptual dapat pula ditinjau dan diamati. - Uji coba sistem dapat dilakukan dengan lebih cepat. - Biaya lebih murah dibandingkan dengan membuat sistem itu sendiri. - Penelaahan sistem secara hipotetis dapat dilakukan, yang pada keadaan sesungguhnya tidak dapat dilakukan. - Simulasi dengan komputer adalah cara yang aman dan dapat diterima untuk menganalisis sistem. Dalam mempelajari fungsi pindah suatu reaktor, di sini dibatasi oleh suatu diagram blok yang cukup sederhana, dimana reaktor hanya
22
sebagai media pengali yang berhubungan dengan besaran p, ~, 1, dan n. Pembatasan masalah dalam studi awal simulasi reaktor ini digunakan data-data dari reaktor berdaya sangat rendah yaitu reaktor SUR- 100 (SiemensUnterrichts-:Reactor) yang dapat menghasilkan daya 100 m\v, sehingga pengaruh temperatur dan xenon terhadap reaktivitas dapat diabaikan. Besaran-besaran yang disebutkan di atas dapat diukur lebih dulu dan selanjutnya dengan perhitungan yang berhubungan dengan reaktor dapat ditentukan. yang lain dianggap no!.
Persamaan diferensial dari reaktor kemudian diselesaikan dengan paket program TUTSIM, sedangkan sinyal data hasilnya diproses dengan paket program FANSIM. Hasil simulasi reaktor SUR 100 dengan program TUTSIM dan FANSIM memberikan gambaran respo waktu dan rapat neutron dengan reaktivitas sebagai parameter, dan daerah kestabilan reaktor dapat diketahui dari respon frekuensi. BAHAN DAN TATA KERJA Reaktor SUR -100 (Siemens-U nterrichts-Reactor)
SUR-100 adalah reaktor termis dengan inti yang homogen dan padat, serta menghasilkan daya 100 mW.
Proceedings Seminar Reaktor Nuklir dalam Penelitian Sains dan Tekrwlogi MenuJu Era Tinggal Landas
Data-data reaktor SUR-100 adalah sbb : Tipe reaktor : termis, padat-homogen Bahan bakar : U30S' pengayaan 20% U-235 Moderator : polyethylen (homogen dengan U30S) Inti reaktor : sHinder (Dia. 24 em, tinggi 26 em)
Batang kendali : cadmium 2 buah Sumber neutron : 10 mCi-Radium-Berillyum Daya : 100 mW Kerapatan fluks : 6.106 Neutron /em2 detik neutron termis maks. Reaktivitas : 0,6 $ / buah rata-rata ti80,400 0,880 3,320 32,800 8,940 (detik) ba tang kendali Koeff. temperatur: - 0,034/ 0,332°C Umur prompt. : 4.10-5detik neutron
Bandung,
8 - 10 Oktober
U,'91 PPTN - BATi\N
pembelahan terjadi yaitu sebagai neutron fierempak, melainkan dilahirkan setelah mengalami proses peluruhan beberapa isotop hasil pembelahan tertentu, yaitu sebagai neutron kasip. Bagian dari neutron kasip ini tergantung dELri bahan yang terbelah. Sebagai eontoh : U-235 menghasilkan % 0,7'% neutron kasip dari seluruh neutron yang dihaTabell. Sifat-sifat neutron kasip yang ditimbulkan akibat terbelahnya U-235 ruh dari 0,00125 0,0124 0,00074 0,3010 1,1400 0,00253 0,1110 0,00027 0,00140 0,00021 0,0305 seluruh jumKonstanta pa\Bagian (detik -1) lah neutron r1i Umur neutron 3,0100
Kinetika reaktor
Seeara pendekatan distribusi ruang kerapatan neutron dalam inti reaktor dapat digambarkan melalui persamaan titik kinetika reaktor : 6
I I
dn/dt = n(p -~) / +
AiCi
(1)
L-1
silkan( Tabel1). Kondisi suatu reaktor didefinisikan olf~h faktor pengali k, k=l dapat disebut sebagai reaktor kritis, k <1 di bawah kritis dan k>l di at:lS kritis. Faktor pengali dapat juga didefinisikan sebagai hubunganjumlah neutron generasi tertentu dibagi jumlah neutron generasi sebelumnya. Besaran k diberikan sebagai faktor pengali untuk reaktor yang tak terhingga besarnya. Untuk reaktor terhingga atau untuk reaktor normal digunakan faktor pengali keff (k effektif),
(2) ~ =
I ~i
n = kerapatan neutron (em-3) p = reaktifitas Ci = konsentrasi isotop induk neutron kasip grup i (em-3) ~i = bagian neutron kasip grup i Ai = konstanta paruh neutron kasip grup
(detik-I)
Persamaan (1) dapat diartikan sebagai berikut: dn/dt = perubahan kerapatan neutron/waktu (p/l)n = neutron yang dihasilkan melalui pembelahan - (~/l)n = neutronyangdipertahankan (neutron lambat)
p = (keff-
Teknik pengaturan dilaksanakan berdasarkan pada teknik umpan balik dan analisis sistem seeara linier serta memadukan kOnSepkonsep teorijaringan. Suatu sistem pengaturan adalah suatu hubungan timbal balik antara komponen - komponen yang membentuk suatu konfigurasi sistem yang memberikan hasil yang dikehendaki. Hasil ini dinamakan respon sistim. Sedangkan dasar untuk menganalisis suatu sistem adalah landasan yang diberilum oleh teori sistem linier, yang menganggap ada-
6
I AiCi= neutron i-I
1) / keff
Teknik pengaturan.
kasip yang timbul
Semua neutron, yang timbul karena proses pembelahan, dilepaskan spontan dan berumur 1.Umur neutron I adalah waktu rata-rata antara dua generasi neutron yang berurutan. Tetapi kenyataannya, sebagian keeil dari neutron yang dihasilkan karena proses pembelahan, tidak langsung dihasilkan pada saat proses
23
Bandung,
Proceedings Seminar Reakwr Nuklir dalmn Penelitian Sains dun Teknologi MenuJu Era Tinggal Landas
Okwber 1991 PPTN - BATAN
Setelah terjadi perubahan reaktivitas maka kerapatan neutron dalam reaktor juga berbentuk sinus dengan harga dasar nO.
nya hubungan linier antara sebab dan akibat suatu sistem. Untuk memahami sistem pengaturan yang rumit, maka lebih dulu harus didapatkan model matematisnya yang bersifat kuantitatif. I<~rena itu adalah perlu untuk menggunakan hubungan antara variabel sistem dan model matematisnya. Karena sistem yang ditinjau berada dalam keadaan dinamis, berubah-ubah, persamaan yang melukiskannya adalah persamaan differensial. Untuk menangani persoalanpersoalan tersebut, secara ringkas dapat diikuti cara-cara tersebut di bawah ini: 'rentukan sistem dan bagian-bagiannya. -Rumuskan model matematisnya dan susunlah asumsi-asumsi yang perlu. -'l'uliskan persamaan differensial yang menggambarkan cara kerja sistem. -Selesaikanlah persamaan untuk variabel keluaran yang dikehendaki. -Periksa hasil penyelesaian - penyelesaian itu dan pengandaian yang dipakai, lalu -Analisa kembali. Suatu pendekatan yang sangat praktis llntuk menganalisis suatu sistem ialah dengan cara respon frekuensi. Tanggapan frekuensi suatu sistem didefinisikan sebagai respon keadaan mantap sistem tersebut terhadap sinyal rnasukan sinusoidal. Sinausoidal adalah masukan yang unik, dan sinyal keluaran yang dihasilkannya untuk sistem linier adalah juga sinusoidal dalam keadaan mantap yang berlainan bentuknya dari gelombang yang pertama hanya pada amplituda dan sudut fasenya. Salah satu keuntungan dari cara respon frekuensi ini ialah tersedianya sinyal penguji dalam praktek untuk berbagai frekuensi dan amplituda. ICarena itu penentuan eksperimentasi untuk respon frekuensi dapat dengan mudah diperoleh, dan merupakan cara yang paling dapat diandalkan dan tidak rumit untuk menganalisis s,istem secara eksperimental. Pada percobaan ini sebagai keluaran mula - mula ditentukan bahwa reaktor beroperasi pada daya tertentu dengan reaktivitas awal sarna dengan no!. Dengan menggerakkan batang pengendali maka reaktivitas dirubah sesuai dengan fungsi sinus yang merupakan fungsi waktu. p = a sin wt
8 - 10
n = nO+ b sin (wt +
(7)
cj»
Dari hubungan b/a dan pergeseran fase yang merupakan fungsi frekuensi, dapat diketahui besarnya amplitudo dan sistem fase tersebut, sehingga dapat diketahui kestabilan dari sistern. cj>
,/
," "-
w-2
/
,,, " I "
I
-30°
;0-310-210-1
1
10
102 103
fr~kUVllbi
10'
105 lOb
Jctlk-1
Gambar 1. Perubahan amplituda dan fase pada reaktor tanpa pengaruh umpan balik dengan 1= 10 -5 detik
(6)
24
Dalam teknik kontrol, garis - garis lurus yang digambarkan dalam diagram Bode (Gambar 1) dapat diartikan sbb: Penurunan penguatan terhadap frekuensi pada frekuensi rendah merupakan hubungan integral. Perubahan amplitudo suatu pengatur integral dalam diagram Bode dapat dilihat sebagai garis lurus dengan kenaikkan-1. Pengatur diferensial orde satu dalam diagram Bode digambarkan sebagai garis horisontal dengan gradien + 1. Garis patah-patah pada Gambar1 mula-mula turun dengan naiknya frekuensi. Mula mula besarnya amplitudo konstan dari frekuensi w1 sampai w2, dan di sini garis mulai turun lagi. Untuk menampilkan proses ini dalam diagram Bode, maka garis yang dibentuk oleh pengatur integral di-
Proceedings Seminar Reaktor Nuklir dalum Penelitian Sains dan Tekrwlogi Menuju Era Tinggal Landas
jumlahkan dengan garis lurus yang dibentuk oleh pengatur differensial, dimulai dari frekuensi wl. Untuk menggambarkan penurunan penguatan pada frekuensi tinggi , harus ditambahkan kembali penguatan yang dihasilkan oleh pengatur diferensial. Oleh sebab itu timbul kembali garis patah-patah dengan sudut tertentu pada frekuensi w2. Ditinjau dari segi teori kontrol dapat dilihat bahwa reaktor merupakan sistem yang tidak linier. Bila orang akan melihat dengan bantuan teori kontrol yang linier maka harus dibuat pendekatan sehingga sistem tersebut dapat dianggap sebagai sistem yang linier. Dengan menganggap bahwa perubahan kerapatan neutron dalam bentuk sinus dengan amplitudo yang kecil dibandingkan dengan nilai dasar, maka masalah ini dapat dianggap sebagai linierisasi. Untuk mengetahui besarnya penguatan maka masukan yang berbentuk reaktivitas p diubah menjadi gelombang yang harmonis. Perhitungan akan menjadi lebih mudah bila gelombang ditulis dalam bentuk bilangan kompleks sebagai berikut: p = ~p eiwt
Bandung,
8 -10
Oktober lS'91 PPTN - BATAN
6
Fr(w) = ~n/~p = (no / jw).l/ (1+~ (~i / i-I
(Ai
+ jw))) (13)
~n= amplitudo gelombang harmon is disekitar harga dasar nO ~p= perubahan reaktivitas yang kecil, dimana besarnya amplitudo gelombang harmonis sekitar harga dasar nol. I = umur serentak neutron rata-rata ~i = bagian dari neutron kasip Ai =konstanta pembelahan isotop penghafiil neutron kasip
(8) ..
Karena biasanya reaktor dianggap sebagai fungsi pindah maka oleh gelombang yang harmonis dari suatu gangguan akan mengakibatkan bergelombangnya semua variabellainnya, misalnya kerapatan neutron n atau konsentrasi Ciisotop penghasil neutron kasip dengan frekuensi yang sarna seperti pada gangguan. Untuk kedua besaran itu dapat dituliskan : n = nO + ~n ejwt C.1 = Go 1 + ~C.1 eiwt
fr~I:i."'.t\'it del.ik-1
Gambar 2. Perubahan amplitudo dan fase pada reaktor tanpa pengaruh umpan balik untuk bermacam-macam umur serentak neutron bahan belah U-235.
(9) (10)
Dalam rumus ini pergeseran fase ditulis tidak eksplisit. Secara linier, persamaan diatas dapat ditulis : 6
jw~n = (no/l)~p- jw ~ ~Ci i-I
(11)
~ =jumlah ~i dari ~1 sampai ~6 Gambar 2 menunjukkan harga dan penguatan fase (~n/~p)(~/no)bahan belah U-285 dalam grafik logaritmis bermacam-maeam harga umur serentak neutron rata-rata. Dalam hal ini dapat dilihat, bahwa reaktor yang dilinierkan pada frekuensi yang sangat keeil menunjukkan hubungan integrasi. Harga pada kurva pengatur integral digambarkan sebagai garis lurus dalam grafik logaritmis, yang berjalan dari kiri atas ke kanan bawah. Sedangkan penguatan didasarkan pada kestabilan.
(12) Diagram per!>'amaan dasar (kinetika reaktor).
Dari persamaan ini dapat dituliskan persamaan penguatan hubungan sebab akibat reaktor tanpa pengaruh umpan balik.
Dengan mempergunakan bahasa yang diperlukan program TUTSIM dan FANSIM maka persamaan di atas dapat diubah menjadi diagram sebagai berikut:
25
Bandung,
Proceedings Seminar Reaktor Nuklir cktlam Penelitian Sains dan Teknologi Menuju Era Tinggal Landas
26
8 - 10
Oktober 1991 PPTN - BATAN
UL27 UL DIV MUL
Proceeding8 Seminar Reakwr Nuklir dalwn Penelitian Sains dan Tekrwlogi Menuju Era Tinggal Landas
Program dengan perangkat /unak TUTSIM dan FANSIM.
-54 54 56 731 19 518 30 626 15 10 -52 -56 23 14 55 157 16 8100 313 12 11 97 52 11 50 531951 51 429 17 -7 719 58 -62 -60 1 32 28 24 31 27 62 60 -58 22 25 21 20 23 Model:: Neutron (nfcm3) 40,000E-06
211 CON 15 7INT 23INT ICON 19 INT 54 31 27INT 50 56 53 55 52 57 51 59 58 60 CON 63 61 62
: Umur Reaktivitas :: Neutron lambat neutron (detik) 'peluruhan Sumber neutron (nfcm3) Konstanta peluruhan
27
Bandung,
8 - 10
Okwber 1991 PPTN - BA1'AN
02 104TIM SIN UM IV 101
Bandung,
Proceedings Seminar Reaktor Nuklir dalcun Penelitian Sains clan Teknolcgi Menuju Era Tinggal Landas
104 101 79 892 90 27 120 82 70 120 76 81 77 91 103 120 : -75 relay :relay 0,0376710 76 SUM
8 - 10
Oktober 1991 PPTN-BATAN
102 75 80 79 REL 91CON 120CON :Batang pengendali 82 92 90 EUL GAl 81 :Kp ::Ki Batang pengendali
HASILDAN DISKUSI Dari kurva rapat neutron terhadap waktu dengan reaktivitas sebagai parameter, terlihat bahwa bilamana reaktivitas diperbesar maka I..' iI ,~ ~--1-.L..L-L-L-.-..l-L· •. u,5 konstanta waktunya makin besar. Hal i.iJm:;Ii ini me-U~ .-- I ;lIlIlronIR 1-(a/cul 1 I i,i1i'?m ....• -I ,. i \ -••__ • _0 iil~~ ~i " 1/ i.i!~~ I nunjukkan program komputer sebagai simulai £'Jektu I .' i,l~ U~ U4~ .::: n i.1~11' r tor beJjalan lebih baik untuk reaktivitas yang i.i1~ ( .• -0.1 !.~~i r . V.I lebih besar karena waktu resolusi (sampling) lebih memadai. Dari kurva tanggapan frekuensi terlihat bllhwa untuk frekuensi di bawah 0,5 Hz terdapat perbedaan penguatan antara sistem reaktor loop terbuka dan sistem reaktor loop tertutup. Ini disebabkan karena pada sistem loop tertutup terdapat faktor umpan balik yang mempunyai pengaruh dominan ketika masukan frekuensi di bawah 0,5 Hz. Untuk frekuensi di atas 0,5 Hz kedua kurva hnmpir berimpit.
-----,
rI
,
nnlll
c,
r
1.11 iWi\J
%!'i
I
KESIMPULAN Pada studi awal dari simulasi reaktor dengan daya rendah, perangkat lunak TUTSIM dan FANSIM dapat menyelesaikan persamaan kinetika reaktor titik dengan baik. Hal ini dapat dilihat dari diagram respon waktu dan respon frekuensi yang dihasilkan.
Untuk daya yang lebih tinggi dimana terdapat pengaruh suhu dan racun xenon akan diteliti lebih lanjut pada penelitian yang akan datang.
28
Ban-dung, 8 - 10 Oktober 1991 PPTN - BA1'AN
Proceedings Seminar Reaktor Nuklir dalam Penelitian Sains don Tekrwlogi Menuju Era Tinggal Lan-das
DAFfAR
PUSTAKA
1. Geoffrey Gordon, Systemsimulation, R.OldenbourgVerlag Munchen & Wien (1972). 2. M.A.Schultz, Control of Nuclear Reactors and Power Plans, Mc. Graw-Hill Book Company Inc. (1961). 3. Werner Jentsch, Digitale Simulation Kontinuerlicher Systeme, R. Oldenbourg Verlag Munchen & Wien (1969). DISKUSI
Setyadi : Dart presentasi yang anda sampaikan, bahwa daya reaktor 100 mWatt dan dilakukan di Jerman. Apakah andajuga bisa melakukan hal tersebut di PPTN, dan apakah hambatannya. Sudjatmi : Hulah yang kami harapkan. Motor penggerak batang kendali harus dibuat lebih dulu sehingga batang kendali dapat diperintah oleh komputer. Budiono : 1. Gambar orang yang ditayangkan sistem manual, tetapi yang konvensional adalah sistem kontrol otomatis yang menggerakkan P.I, P.LD. dan lain-lain dengan cara analog. 2. Bagaimana cara menghubungkan komputer dengan batang kendali ? Apakah sistem kontrol yang menggunakan 2 relai tidak on-off? Bagaimana hubungan P.LD. dengan relai ? 3. Respons yang menggunakan kontroler tersebut selalu ada offset, apa benar hasil yang ditayangkan? Saudara menggunakan paket program, kemudian buat percobaan, apakah tidak lebih baik dibuktikan dengan perhitungan secara teori respons tersebut ? Sujatmi : 2. Dengan melalui interface (AID dan D/A), batang kendali dihubungkan dengan komputer. Pada dasarnya P.I.D. berjalan terus tanpa memerlukan ani off atau relai. Tapi disini onl off atau relai digunakan untuk membatasi kerja P.LD. Untuk gangguan yang sangat kecil (pada batas yang diberikan untuk on) relai tidak on, sehingga P.LD. tidak berfungsi. 3. Benar. Kalau anda berminat untuk membuktikan dengan perhitungan, silahkan ! 29