I 2 N N 2/AU2/P2MI-LIPI/04/2012 ( 2 M M N M -GAS Berbahan ......... ( L M )

## MANAJEMEN KONVERSI TERAS  RSG-GAS BERBAHAN BAKAR SILISIDA TINGKAT MUAT TINGGI
Lily Suparlina Pusat Teknologi Reaktor dan Keselamatan Nuklir
Gedung 80 , kawasan Puspiptek, Serpong
lilyrsg@batan.go.id
Diterima editor 26 Agustus 2013 Disetujui untuk publikasi 23 September 2013

## ABSTRAK
MANAJEMEN  KONVERSI  TERAS  RSG-GAS  BERBAHAN  BAKAR  SILISIDA  TINGKAT MUAT  TINGGI .  Penggunaan  bahan  bakar  tingkat  muat  tinggi  dapat  memperpanjang  siklus  operasi reaktor  sampai  40  hari.  Telah  dilakukan  perancangan  konversi  teras  dari  silisida  tingkat  muat  rendah menuju  tingkat  muat  tinggi.  Manajemen  konversi  teras  dari  teras  silisida  tingkat  muat  2,96  gU/cm 3 menuju  teras  silisida  4,8  gU/cm 3 dilakukan  secara  bertahap  dengan  melakukan  perhitungan  manajemen bahan  bakar  dalam  teras  yang  menggunakan  paket  program  perhitungan  2  dimensi  Batan-FUEL. Penggantian  bahan bakar di teras menggunakan pola pergeseran bahan bakar 5/1 yaitu setiap awal siklus terjadi penggantian 5 buah elemen bahan bakar standar dan 1 buah elemen bahan bakar kendali. Dengan mempertahankan konfigurasi teras yang sudah ada, konversi teras dapat dilakukan melalui simulasi teras campuran  2,96  gU/cm 3 -  4,8  gU/cm 3 dengan  memperhatikan  batasan  keselamatan  reaktor  yang dipersyaratkan.  Oleh  karena  itu  tujuan  dari  penelitian  ini  adalah  untuk  merancang  teras  campuran  yang akan menghasilkan parameter neutronik yang aman pada teras penuh pertama silisida tingkat muat tinggi. Dalam  perancangan  ini  digunakan  batang  kendali  pengaman  untuk  menambah  margin  padam  yang berkurang akibat pemuatan bahan bakar tingkat muat tinggi. Hasil analisis menunjukkan bahwa konversi teras silisida tingkat muat 2,96 g U/cm 3 menuju 4,8 g U/cm 3 dapat dilakukan melalui teras campuran tidak langsung dalam 2 tahap yaitu konversi teras silisida 2,96 gU/cm 3 - 3,55 gU/cm 3 dan konversi teras silisida 3,55 gU/cm 3 - 4,8 gU/cm 3 dengan performa yang baik. Keuntungan utama dari penggunaan bahan bakar silisida  tingkat  muat  tinggi  4,8  gU/cm 3 dibanding  teras  silisida  tingkat  muat  rendah  2,96  gU/cm 3 pada teras  RSG-GAS  ialah  bahwa  panjang  siklus  operasi  dapat  lebih  panjang  18  hari  sehingga  dapat menghemat penggunaan bahan bakar.
Kata kunci : silisida, BKP, teras campuran, pola 5/1, Batan-FUEL

## ABSTRACT
MANAGEMENT  OF  RSG-GAS  HIGH  DENSITY  SILICIDE  FUELED  CORE .  The  usage  of  high density  fuel  can  extend  the  reactor  operation  up  to  40  days.  Designing  of  low  to  high  density  silicide fueled  core  conversion  has  been  carried  out.  The  management  of  core  conversion  was  done  gradually using 2 dimensional diffusion code Batan-FUEL. Replacement of unused fuel elements in core using fuel elements  reshuffle  pattern  5/1,  which  at  the  beginning  of  cycle  there  are  5  fuel  elemnts  and  1  control element replaced. By maintaining the existing core configuration, the core conversion can be performed through  mix  core  of    2,96  gU/cm 3 -  4,8  gU/cm 3 silicide  fueled  with  respect  to  reactor  safety  limits requirements.  Therefore,  the  objective  of  this  work  is  to  design  the  mixed  cores  on  the  neutronic performance to achieve safely at first full-silicide core for the reactor with the high uranium meat density. In  the  design,  safety  rods  were  used  to  increase  the  decreasing  minimum  shutdown  margin  due  to  high density  fuel  loading.  The  analyses  results  show  that    silicide  core  conversion  of  2,96  gU/cm 3 to  4,8 gU/cm 3 density can be performed by means of indirect mix core in two steps i.e. silicide of 2,96 g U/cm 3 - 3,55 gU/cm 3 core conversion and silicide of 3,55 gU/cm 3 - 4,8 gU/cm 3 conversion with good performance. The  main  advantage  of  using  high  density  silicide  fuels  of  4.8  gU/cm 3 than  low  density  silicide  fuels  of 2.96 gU/cm 3 on the RSG-GAS core is that, the  operating cycle length of 18 days could be longer in order to save fuel usage.
Keywords : silicide, safety rod, mix core, 5.1 pattern, Batan-FUEL
J " #$ % " & $ '% ) * + " , -% . " V *." / 5 No.3 Oktober 2013, Hal. 137-149 03 '4' 6$ + . ' % - 7 % + $ 89 )'4 9 4:8 7 ; + 9 . <= / 5 > ?@@ , / A // – 240X Nomor : 402/AU2/P2MI-LIPI/04/2012

## PENDAHULUAN
Reaktor Serba Guna G.A Siwabessy (RSG-GAS) yang saat ini beroperasi telah mengalami perubahan penggunaan bahan bakar tipe pelat atau MTR (Material Testing Reactor) dari U 3 O 8 -Al menjadi  U 3 Si 2 -Al  dengan  tingkat  muat  sama  yaitu  2,96  gU/cm 3 .  Konversi  teras  RSG-GAS  dari teras berbahan bakar oksida menuju teras berbahan bakar silisida dilakukan melalui teras campuran oksida-silisida dimulai dari teras 36. Pada teras 45 seluruh teras sudah berisi bahan bakar silisida 2,96 gU/cm 3 pada tahun 2002 dan tidak terjadi perubahan parameter neutronik yang siginfikan pada teras silisida 2,96 gU/cm 3 dibandingkan dengan teras oksida 2,96 gU/cm 3 [1-5].
Untuk meningkatkan panjang siklus operasi reaktor, maka telah dilakukan penelitian tentang penggunaan  bahan bakar silisida dengan tingkat  muat 3,55  gU/cm 3 .  Silisida  dengan  tingkat  muat 3,55  gU/cm 3 dipilih  karena  tidak  perlu  mengubah  konfigurasi  dan  parameter  teras  lainnya.  Hasil perhitungan  menunjukkan    bahwa  penggunaan  bahan  bakar  silisida  dengan  tingkat  muat  3,55 gU/cm 3 dapat menambah panjang siklus operasi reaktor sebesar 7 (tujuh) hari dengan daya penuh (210 MWD)[6].
Untuk  lebih  meningkatkan  panjang  siklus  operasi,  maka  dilakukan  penelitian  penggunaan bahan bakar dengan tingkat muat 4,8 gU/cm 3 pada teras RSG-GAS. Pada saat ini beberapa reaktor riset yang menggunakan bahan bakar silisida tingkat muat tinggi, beberapa diantaranya ialah JMTR dan  JRR-3  di  Jepang  serta  Opal  di  Australia  dan  yang  sedang  dalam  taraf  pembangunan diantaranya reaktor Pallas di Belanda, JHR di Perancis dan JRTR di Jordan.  Beberapa penelitian tentang  penggunaan  bahan  bakar  silisida  tingkat  muat  4,8  gU/cm 3 telah  dilakukan  sebelumnya, namun belum diperhitungkan cara mencapai teras setimbang[7-13].  Dalam penelitian ini dilakukan cara  yang  paling  efisien  untuk  mencapai  teras  setimbang  silisida  4,8  gU/cm 3 tanpa  mengubah konfigurasi teras dan parameter teras serta memenuhi batasan operasi reaktor saat ini. Konfigurasi teras RSG-GAS saat ini digunakan dalam konversi teras silisida 2,96 gU/cm 3 ke silisida 4,8 gU/cm 3 tanpa  mengubah  konfigurasi  serta  parameter  teras.  Peningkatan  tingkat  muat  bahan  bakar  maka akan  terjadi  kenaikan  reaktivitas  yang  tinggi  dan  untuk  itu  ditambah  batang  kendali  pengaman untuk  mengkompensasi  margin  reaktivitas  padam  agar  syarat  keselamatan  operasi  reaktor terpenuhi.
Tujuan dari program konversi teras adalah untuk mendapatkan teras setimbang yang optimal yaitu  teras  setimbang  yang  terbaik  dari  segi  operasi  reaktor  dengan  siklus  operasi  yang  lebih panjang  sehingga  dapat  menghemat  bahan  bakar  serta  strategi  untuk  mencapai  teras  setimbang. Program  konversi  dari  teras  silisida  tingkat  muat  2,96  gU/cm 3 menuju  teras  setimbang  silisida tingkat muat 4,8 gU/cm 3 dilakukan melalui teras campuran silisida 2,96 gU/cm 3 dan 3,55 gU/cm 3 serta 4,8 gU/cm 3 .  Berdasarkan hal itu maka sebelum melaksanakan konversi teras, perlu dilakukan perhitungan manajemen konversi teras RSG-GAS.
Perhitungan  konversi  teras  campuran  dilakukan  dalam  2  (dua)  tahap,  pertama  perhitungan sel  untuk  menggenerasi  tampang  lintang  material  teras  dengan  menggunakan  program WIMSD/B5[14]  dan  kedua  adalah  perhitungan  teras  dengan  metode  difusi  neutron  2  dimensi Batan-EQUIL-2D  yang  merupakan  bagian  dari  Batan-FUEL[15].  Kajian  yang  dilakukan  dalam perhitungan  teras  ini  hanya  ditinjau  dari  aspek  neutroniknya  saja.  Program  konversi  teras  RSG- GAS diharapkan dapat menghasilkan desain teras silisida setimbang yang optimal.
P E F GE H K OH QH R -FUEL
Batan-FUEL adalah paket program perhitungan manajemen bahan bakar di teras reaktor yang sudah dikembangkan dan diverifikasi Program Batan-FUEL merupakan paket perhitungan parameter neutronik yang menggunakan metode difusi neutron banyak kelompok energi  dengan model geometri dua dan tiga dimensi. Saat ini Batan-FUEL digunakan dalam manajemen bahan bakar di teras reaktor RSG-GAS
ST TU W XWW YZX [\ U ] ^ ] _ ` X[ 2/AU2/P2MI-LIPI/04/2012 ab c d c e f _ g ch i j h _ f kl m cd l d n k j o _ l g Z[Wp q b c r c s f ^ f r t ] r u f _ d l vf _cd w T x -GAS Berbahan ......... a y l g z T i oc _ glr c )
Paket program Batan-FUEL terdiri dari 3 bagian, yaitu :
1. Batan-EQUIL  digunakan  untuk  pencarian  teras  setimbang  dan  fraksi  bakar  di  awal  dan akhir siklus.
2. Batan-2DIFF  digunakan  untuk  perhitungan  parameter  neutronik  2  dimensi  ke  arah X-Y. Data parameter teras yang dihasilkan dalam perhitungan ini adalah reaktivitas lebih teras, marjin padam, reaktivitas batang kendali, reaktivitas xenon, perubahan reaktivitas dalam satu  siklus  operasi  ( burnup  swing) ,  perubahan  reaktivitas  dari  panas  ke  dingin  ( hot  to cold) , faktor puncak daya radial dan fluks neutron radial.
3. Batan-3DIFF digunakan untuk perhitungan parameter neutronik 3 dimensi dalam arah X- Y-Z seperti reaktivitas sebagai fungsi pemasukan batang kendali, fluks neutron aksial dan faktor puncak daya aksial.

## Manajemen Bahan Bakar Dalam Teras
Manajemen  bahan  bakar  dalam    teras  atau in-core  fuel  management ,  ialah  penentuan  pola penyusunan  bahan  bakar,  penentuan  fraksi  bakar  ( fuel  burn-up )  dan  seleksi  jadwal  pergantian bahan  bakar,  dalam  rangka  untuk  memenuhi  persyaratan  dari  aspek  nuklir  maupun  ekonomi. Manajemen  bahan  bakar  di  dalam  teras  meliputi  kegiatan  pembentukan  konfigurasi  teras, penentuan  fraksi  bakar,  panjang  siklus  operasi,  pola  pemasukan  bahan  bakar  dan  pelaksanannya yang aman.
Untuk mendukung keselamatan dengan menyediakan reaktivitas negatif dari batang kendali, maka reaktor RSG-GAS didesain harus tetap dalam keadaan subkritis, meskipun satu buah batang kendali  yang  memiliki  reaktivitas  terbesar  tidak  dapat  dimasukkan  secara  penuh  ke  dalam  teras. Kondisi  ini  disebut  sebagai  kondisi  satu  batang  kendali  gagal  masuk  atau one  stuck  rod . Pengendalian  reaktivitas  di  reaktor  RSG-GAS  dilakukan  oleh  8  batang  kendali  yang  didesain sebagai  batang  kendali  pemadam  disamping  berfungsi  untuk  mengkompensasi  reaktivitas  lebih teras  selama  satu  siklus  operasi  dan  mengatur  tingkat  daya.  Konfigurasi  teras  setimbang  reaktor RSG-GAS disajikan pada Gambar 1.
Gambar 1. Teras setimbang silisida teras RSG-GAS
~                    No.3 Oktober 2013, Hal. 137-149                               – 240X Nomor : 402/AU2/P2MI-LIPI/04/2012

## Pola Pemuatan Bahan Bakar Teras Setimbang
Tujuan dan persyaratan pemasukan bahan bakar secara periodik ke teras adalah:
a. Menambah  jumlah  reaktivitas  lebih  teras  (dalam  batas  yang  aman)  agar  dapat  dioperasikan pada daya penuh selama panjang siklus operasi yang tertentu
b. Mendapatkan distribusi pembangkitan panas yang merata di teras
c. Menjamin bahwa fluks panas di sistem proteksi reaktor tidak melewati nilai batas
d. Menjamin  bahwa  tingkat  fluks  neutron  pada  posisi  iradiasi  tetap  dijaga  pada  tingkat  yang diinginkan.
Dengan  pola  5/1  terdapat  5  buah  elemen  bakar  standar    (EBS)  dan  sebuah  elemen  bakar kendali  (EBK)  segar  yang  dimasukkan  ke  dalam  teras  di  awal  siklus  sebagai  ganti  dari  5  buah elemen  bakar  standar  dan  sebuah  elemen  kendali  yang  sudah  mencapai  fraksi  bakar  buang maksimumnya.  Keakuratan  pola  pergeseran  ini  telah  dibuktikan  dalam  pelaksanaan  rutin manajemen bahan bakar teras RSG-GAS sejak tahun 1999.  Pola pergeseran bahan bakar di teras RSG-GAS ditunjukkan pada Tabel 1.
Tabel 1. Manajemen bahan bakar di teras RSG-GAS dengan pola 5/1[19]
¡ ¢ £¤£¤ ¥ ¤ ¦§ ¨ © £ ª§ «© £ ¬¨ © ­ £¤ bakar ® ¯ ° ± ² ³ ´ µ ¡¢ £¤£¤ ¥¤ ¦§ ¨ © £ ¶ -9 · -10 ¸ -9 ¹ -8 º -3 » -6 ¸ -5 ¸ -8 ¶ -4 · -9 ¼ -5 ¶ -6 º -10 ¸ -4 » -7 ½ -8 · -3 » -10 ¹ -3 ¼ -7 ¶ -7 · -7 · -4 · -6 » -3 ¶ -8 » -4 º -5 ¶ -5 ¹ -5 ¼ -8 ½ -5 ¼ -9 ¼ -4 ¹ -10 ½ -4 ¼ -6 ½ -9 » -9 º -8 » -8 · -5 · -8 » -5 º -4 ¹ -9 ¸ -6 ½ -7

## Bahan Bakar Silisida
Penelitian  dan  pengembangan  bahan  bakar  berpengayaan  rendah  telah  banyak  mengalami kemajuan dalam program Reduced Enrichment Research and Test Reactot (RERTR) .  Pada saat ini bahan bakar yang bertingkat muat paling tinggi yang telah dikualifikasi untuk reaktor riset adalah bahan  bakar  pengayaan  rendah  silisida  (U3Si2-Al)  dengan  tingkat  muat  uranium  dalam  meat sebesar 4,8 g U/cm 3 .
Satu perangkat bahan bakar standar terdiri dari 21 pelat elemen bahan bakar dengan dan satu perangkat  bahan  kendali  terdiri  dari  15  pelat  elemen  bakar  dengan  2  sisi  terluar  berisi  bahan penyerap AgInCd yang berbentuk garpu.
Data  teknis,  gambar  irisan  penampang  lintang  bahan  bakar  standar  dan  bahan  bakar  pada elemen kendali ditunjukkan padaTabel 2,

## Tabel 2. Data teknis bahan bakar silisida pada teras reaktor RSG-GAS[1]
¾ ¿À ¿ÁÂ Ã Â À ÄÅ Â Á Â Æ Ç ¿È¿À É¿ Ã¿ Ê Â ÈÆ ËÌ ÍÎ Ï Ð ÑÒ Î Ð Ó Ð Ï Ð Ñ ÔÕ Ö Õ × Ø Ö Ð ÑÙ Õ Ó Î ÚÏÏÛ ÜÜÝÞßàÞßá 00 â Ð ã Ð ÔÕ ÓÕ Ñ ä Ð ÓÕ ã ÔÕ Ö Õ × ÚÏ Ï Û ÞÝ å æ Ð ÔÕ × Ö Õ Ñ Õ Ó ä Ð ÑÙ Î Ñ ç Î Ñ (mm) è Ý éé ê ë ÏÓ Õ ì ä Ð Ó Õ ã ä Ð × Ð ÓÐ Ï Ð Ñ Ô Õ Ö Õ × è Þ ê ë ÏÓ Õ ì ä Ð Ó Õ ã ä Ð × Ð ÓÐ Ï Ð Ñ Ö Ð ÑÙ Õ Ó Î Þ é í Õ ã Ð × ÎÕ Ó Ö Ð Ó î Ñ ç Ò î Ñ ç ÔÕ ì Õ Ñ ÔÕ ÖÕ × ï Ó í ç 2 â Ð ã Ð ÔÕ ÓÕ Ñ Ö Ð Ó î Ñ ç Ò î Ñ ç ÔÕ ì Õ Ñ ÔÕ Ö Õ × ÚÏ ÏÛ ð Ý å à ÍÎ Ï Ð ÑÒ Î fuel meat (mm) ð Ý éñ ß á è ÝÜ 5x600 í Õ ã Ð × ÎÕ Ó ò uel meat ó 3Si2Al
ôõ õö ÷ ø÷÷ ùúø ûü ö ý þ ý ÿ N øû 2/AU2/P2MI-LIPI/04/2012 ( ÿ ÿ ÿ úû÷ M þ M ý M   ÿ ÿ õ -GAS Berbahan ......... ( L õ   ÿ M )
P -235  (w/o) !"# T$ % & ' % ) $ '& * + & meat (g/cm 3 ) 2! , - ' % -235  per elemen bakar (g) 2#. - ' % -235 per elemen kendali (g) " / !, -  % )$ + 0 ) 0 A -In-Cd K % 3 + 3 4 0 ) 0 5 && 6 7 ! 7/ -  % )$ + + 8 9 8 0 ) 0 5 & & 6 SS -321 K % 3 + kelongsong penyerap (mm) .! / #

## Batang Kendali Pengaman
Batang  kendali  pengaman  (BKP)  adalah  batang  kendali  yang  dipasang  di  luar  teras  aktif, yang  fungsinya  untuk  meningkatkan  kemampuan  batang  kendali  yang  sudah  ada.  Penambahan batang  kendali  pengaman  diisyaratkan  tidak  merubah  karakteristik  operasi  teras  reaktor  baik  dari sisi neutronik maupun termohidrolik. Untuk memenuhi syarat tersebut pemilihan penambahan BKP berdasarkan fungsi dan batasan sebagai berikut:
- BKP  hanya  berfungsi  sebagai  batang  kendali  pengaman  yang  terintegrasi  dengan  sistem pancung reaktor.
- Posisi  BKP  diusahakan  berada  di  luar  daerah  bahan  bakar  yaitu  di  daerah  reflektor  dan simetris, sehingga tidak merubah karakteristik operasi reaktor
- Penambahan BKP dapat menghasilkan margin reaktivitas padam teras minimum sebesar - 0,5 % D k/k.
Geometri  dari  batang  kendali  pengaman  sama  dengan  geometri  elemen  kendali,  hanya  bagian bahan  bakarnya  diganti  dengan  dummy  aluminium.  Untuk  perhitungan  reaktivitas,  posisi  BKP semua diatas (all up) dan untuk perhitungan margin padam, posisi BKP di bawah (down).

## METODOLOGI

## Generasi Tampang Lintang
Generasi  tampang  lintang  makroskopik  bahan  bakar  fisil  dan  material  teras  lainnya dilakukan  dalam  perhitungan  sel  dengan  paket  program  WIMSD/B5.  Data  masukan  berupa kerapatan  atom  dan  temperatur  bahan  bakar  serta  pembagian  fraksi  bakar  ke  dalam  17 langkahdigunakan  dalam  perhitungan  sel  dalam  WIMSD/B5  dan  akan  dihasilkan  data  keluaran berupa  data  tampang  lintang  makroskopik  bahan  bakar  untuk  kondisi  dingin  tanpa  xenon  dan samarium, dingin tanpa samarium, panas xenon dan samarium setimbang panas tanpa samarium.

## Pemodelan Teras Campuran
Pemodelan teras campuran silisida  2,96 gU/cm 3 – 4,8 gU/cm 3 , 2,96 gU/cm 3 - 3,55 gU/cm 3 , 3,55  gU/cm 3 – 4,8  gU/cm 3 dibuat  berdasarkan  pola  pergeseran  bahan  bakar  5/1  seperti  yang ditunjukkan pada Tabel 3 dan disajikan pada Gambar 2 – 9. Langkah dan strategi pemuatan bahan bakar  ke  dalam  teras  campuran  tersebut  diatur  sedemikian  rupa  sehingga  parameter  teras  seperti reaktivitas  lebih,  margin  padam  minimum  dan  fraksi  bakar  buang  maksimum  di  akhir  siklus memenuhi batasan keselamatan.
J : ;< = : > < ?= @ B C : D E= F : V BF: G H No.3 Oktober 2013, Hal. 137-149 IO ?Q? R< C F ? = E U = C < WX @?Q X QYW U Z C X F [\ G H ] ^__ D G ` GG – 240X Nomor : 402/AU2/P2MI-LIPI/04/2012
d e f g e h i j kl h e m n e f o p h e q r d e f g e h sj k l h e m n e f o p h e q i d e f g e h t j kl h e m n e f o p h e q s d e f g e h uj k l h e m n e f o p h e q t d e f g e h v j kl h e m campuran 5 d e f g e h w j k l h e m n e f o p h e q v
xy yz { |{{ }| ~ z      |~ 2/AU2/P2MI-LIPI/04/2012                           }~{                     y  -GAS Berbahan .........      y      )
¢ £ ¤ ¥ £ ¦ § ¨ ©ª ¦ £ « ¬ £ ¤ ­ ® ¦ £ ¯ ° ± ª ² ª ¦£ ¯³ £ ¯ ´ ¢ £ ¤ ¥ £ ¦ µ¨ ©ª ¦ £ « ¬ £ ¤ ­ ® ¦ £ ¯ §

## Batas Keselamatan Dalam Perhitungan Teras
Perhitungan  teras  campuran  silsida  langsung  maupun  tidak  langsung  perlu  memperhatikan batasan  keselamatan reaktor  RSG-GAS. Batasan  keselamatan  yang  digunakan  dalam  perhitungan ini ialah :
- Reaktivitas lebih teras awal siklus maksimum 10 % D k/k - Marjin reaktivitas padam minimum ( stuck rod condition ) adalah < -0,5 % D k/k
- Faktor puncak daya (FPD) radial maksimum adalah 1,4.
- Fraksi bakar buang ( discharged burn up ) maksimum di akhir siklus adalah  < 70% burnup - Reaktivitas lebih akhir siklus >2,0 % D k/k - Reaktivitas total target < 2,0 % D k/k
Parameter  yang  dihitung  dalam  mencapai  teras  setimbang  silisida  4,8  gU/cm 3 melalui  teras campuran adalah :
- Harga reaktivitas lebih teras pada saat awal siklus ( BOC) dingin bebas xenon
- Margin padam pada awal siklus saat BOC
- Fraksi bakar buang maksimum
- Faktor puncak daya radial maksimum
- Fluks neutron cepat dan termal
¶ · ¸¹ º · » ¹ ¼º ½ ¾ ¿ · À Áº Â · Ã ¾Â· Ä Å No.3 Oktober 2013, Hal. 137-149 ÆÇ ¼È¼ É¹ ¿ Â ¼ º Á Ê º ¿ ¹ ËÌ ½¼È Ì ÈÍË Ê Î ¿ Ì Â ÏÐ Ä Å Ñ ÒÓÓ À Ä Ô ÄÄ – 240X Nomor : 402/AU2/P2MI-LIPI/04/2012

## HASIL DAN PEMBAHASAN
Perhitungan parameter neutronik melalui konversi teras campuran silisida dari 2,96 gU/cm 3 – 4,8 gU/cm 3 dilakukan dengan menggunakan paket program Batan-2DIFF yang merupakan bagian dari  paket  program  Batan-FUEL  dengan  dua  cara  yaitu  konversi  langsung  dan  konversi  tidak langsung.  Kedua  cara  perhitungan  tersebut  dilakukan  dengan  pembatasan  fraksi  bakar  buang maksimum < 70 burnup dan reaktivitas teras di akhir siklus dalam kondisi panas xenon, samarium setimbang > 2,0 % D k/k.

## Konversi Langsung
Konversi langsung dilakukan dengan perhitungan teras dimulai teras setimbang silisida 2,96 gU/cm 3 . Teras setimbang ini  mempunyai panjang siklus 21,80 hari atau setara dengan energi yang dibangkitkan sebesar 654 MWd. Pengaturan panjang siklus ini dilakukan agar dapat menghasilkan reaktivitas  lebih  di  awal  siklus  yang  cukup,  reaktivitas  lebih  di  akhir  siklus,  fraksi  bakar  buang maksimum  dan  margin  padam  yang  memenuhi  batasan  operasi  reaktor.  Fraksi  bakar  teras setimbang silisida 2,96 gU/cm 3 menurut batasan operasi yang diijinkan oleh Bapeten adalah <60%, sedangkan  fraksi  bakar  untuk  tingkat  muat  yang  lebih  tinggi  menurut  batasan  hasil  penelitian elemen  bakar  adalah  <70  %.  Dengan  mengikuti  langkah  pemuatan  bahan  bakar  sesuai  pola reshuffle ,  maka  di  setiap  awal  siklus  dimasukkan  5  elemen  bahan  bakar  standar  dan  satu  bahan bakar kendali  segar silisida 4,8 gU/cm 3 . Lima bahan bakar standar dan satu bahan bakar kendali silisida  2,96  gU/cm 3 yang  sudah  mencapai  fraksi  bakar  buang  maksimum  dikeluarkan  dari  teras reaktor. Langkah ini dilakukan setiap teras campuran mulai dari teras campuran 1 sampai dengan teras campuran  8, sehingga  di teras campuran  8  seluruh  teras sudah  terisi  seluruhnya  oleh bahan bakar silisida  4,8 gU/cm 3 .
Pada  teras  campuran  ini  terdapat  2  model  konfigurasi  teras  yaitu  pada  langkah  ke  2-a  dan langkah 7-a seperti yang ditunjukkan pada Gambar 10 dan 11.
× Ø Ù Ú Ø Û Ü Ý Þ ß à ÛØ á â Ø Ù ã ä ÛØ å æ ç à å è Ø å æ é êë × Ø Ù Ú Ø Û Ü Ü Þ ß à ÛØ á â Ø Ù ã ä ÛØ å ì ç à å è Ø å í é êë
Hasil perhitungan berupa nilai reaktivitas lebih, fraksi bakar maksimum di akhir siklus dan margin padam pada pembentukan teras silisida 4.8  melalui teras campuran silisida  2,96 gU/cm 3 - 4,8  gU/cm 3 sampai  terbentuknya  teras  penuh  silisida  4.,8  gU/cm 3 dengan  konversi  langsung ditunjukkan pada Tabel 3.
îï ïð ñ òññ óôò õö ð ÷ ø ÷ ù ú òõ 2/AU2/P2MI-LIPI/04/2012 ûü ý þ ý ÿ ( ù ý ù ( ýþ þ ù ôõñ ü ý M ý ( ø ( M ÷ M  ( ù þ ( ùýþ ï -GAS Berbahan ......... û L ï  ý ù M ý )
Tabel 3. Hasil perhitungan konversi langsung teras campuran silisida 2,96 g U/cm 3 - 4,8 g U/cm 3 pada daya 30 MW P bahan bakar Si 4.8 gU/cm 3 (EBS s/EBK) R ! " # " ! $ lebih BOC (% D k/k) F % $ " & % buang maksimum (% D k/k) ' ) * ) + siklus (hari) , % +" ) padam minimum (% D k/k) K ! % ) + ) 0 - . / 2 9 3 40 56 3 0 0 . 7 3 80 2 1,30 S : ; S " 2,96 gU/cm 3 7 5< 7 9 3 98 6 7 3 4 5 = 0 > 00 2 0,70 . 70 < . 7 0 > = 4 2 2 2 0,05 T " ? : ) . -a. 70 < . 9 3 56 5 9 3 4 5 . 5 2 4,28 . @ K ' = 7 5< = 7 0 3 = 8 6 2,52 = 5 2 3,47 A .0 < A 7 7 3 A 4 6 6,33 = 9 2 3,36 5 . 5< 5 7 1,30 6 9,12 A 0 2 2,76 6 =0 < 6 7 1,23 6 9 3 9 0 = 8 2 2,07 4 = 5< 4 7 1,42 6 9 3 4 5 2 2 0,02 T " ? : ) 4 -a = 5< 4 7 . 3 6 0 6 9 3 4 9 = 0 2 1.878 A @ K ' 8 A 0 < 8 7 A 3 7 9 6 9 3 4 5 A 0 2 0.056 T " ? : )
Pada  Tabel  3  dapat  dilihat  bahwa  langkah  ke  0  merupakan  teras  setimbang  silisida  2,96 gU/cm 3 .  Pada  langkah  ke  2,  dimana  10  elemen  bakar  dan  2  elemen  bakar  silisida    4,8  gU/cm 3 dengan fraksi bakar rendah berada dalam teras reaktor bersama dengan silisida 2,96 gU/cm 3 dengan fraksi bakar tinggi,  dihasilkan nilai margin padam sebesar -0,05 % D k/k, yang berarti margin padam tersebut  tidak  terpenuhi  karena  melebihi  batasan  operasi  (<-0,5% D k/k)  sehingga  perhitungan reaktivitas  di  akhir  siklus  tidak  dilakukan.  Agar  margin  padam  dapat  terpenuhi,  maka  perlu ditambahkan  2 buah BKP yang diletakkan pada posisi teras G-10 dan B-3 seperti yang ditunjukkan pada  langkah  ke  2-a,  Gambar  10.  Dalam  perhitungan,  penambahan  2  buah  BKP  ini  dapat dipertahankan  sampai  langkah  ke  6,  karena  pada  langkah  ke  7,  margin  padam  yang  dihasilkan sebesar  -0,02  % D k/k  yang  berarti  syarat  batasan  operasi  tidak  terpenuhi  juga,  sehingga  perlu ditambahkan  lagi  2  buah  BKP  yang  diletakkan  pada  posisi  B-10  dan  G-3  dan  dilakukan perhitungan  dengan  langkah  ke  7-a,  Gambar  11.  Pada  langkah  ke  8,  dimana  seluruh  teras  sudah termuati dengan bahan bakar  4,8 gU/cm 3 , margin padam juga tidak memenuhi kriteria, sehingga perhitungan  dihentikan  dengan  anggapan  model  teras  campuran  dengan  konversi  langsung  dari silisida 2,96 gU/cm 3 ke 4,8 gU/cm 3 tidak dapat dilanjutkan, sehingga perlu dilakukan perhitungan dengan model teras campuran dengan konversi tidak langsung.

## Konversi Tidak Langsung
Konversi tidak langsung dilakukan melalui dua tahap perhitungan teras campuran yaitu tahap pertama  adalah  teras  campuran  2,96  gU/cm 3 – 3,55  gU/cm 3 dan  tahap  ke  dua  3,55  gU/cm 3 – 4,8 gU/cm 3 .  Dengan  menggunakan  pola  pergeseran  bahan  bakar  ( reshuffle) 5/1,  pada  tahap  pertama yaitu teras campuran 2,96 gU/cm 3 – 3,55 gU/cm 3 dilakukan sampai seluruh teras berisi bahan bakar silisida  3,55  gU/cm 3 .  Tahap  ke  dua    dilanjutkan  dengan  konversi  teras  dari  teras  silisida  3,55 gU/cm 3 – 4,8 gU/cm 3 dengan cara yang sama dengan tahap pertana.  Langkah ini dilakukan sampai seluruh  teras  berisi  bahan  bakar  silisida  4,8  gU/cm 3 .  Parameter  neutronik  yang  dihasilkan  dari perhitungan  teras campuran tidak langsung ditunjukkan pada Tabel 4 dan Tabel 5. Pada langkah ini  terdapat  dua  teras  campuran  yang  memerlukan  penambahan  BKP  yaitu  pada  teras  campuran seperti yang ditunjukkan pada Gambar 11 dan 12.
J B CE G B H E IG N O Q B U VG W B X OWB Y Z No.3 Oktober 2013, Hal. 137-149 [\ I]I ^E Q W I G V _ G Q E `a NI] a ]b` _ c Q a W de Y Z f ghh U Y i YY – 240X Nomor : 402/AU2/P2MI-LIPI/04/2012
m n o p n q r s t u v qn w x n o y z qn { | } v {~ n { s  
Tabel 4. Perhitungan konversi tidak langsung  I teras campuran silisida 2,96 gU/cm 3 – 3,55 gU/cm 3 pada daya 30 MW                                3                       ¡ ¢ ∆    £                   ¢ ∆ k/k)   ¤            ¥     ¦ §        ¢ ∆ k/k)          ¨ © ª « ¬ ­  ®¨ l  ¨ ¨ ª j  ¯¨ ¬ 1,30       ª  ­ l gU/cm ° j  j ­  l k  ®  ­ ® ª k  j ® ¬ 5,65 ª j ¨  ª ­  ¯ ª  ¯  l ® ª   ¯  ¬ 4,27  j   ­  ¯ ­  ­  ­¯ ª ®   ¬ 4,22 k ª¨  k ­  ­  l j   ¨ ª ®  l ® ¬ 4,25  ª   ­  ­¯ l ª  ¯ ® ª ¯  l ® ¬ 3,95 l  ¨  l j ¨  j ¯ l k  l ® j   ¬ 3,76 ®  ® j ¨  j ¯ l ®   ¯  ª  ¨¨ ¬ 3,34 ¯ k ¨  ¯ ­  ­ ® l l   ¨  ¨  ¨¨ ¬ 3,82         gU/cm °
Pada teras campuran tidak langsung I ini, langkah ke 0 pada Tabel 4 merupakan data teras setimbang  silisida  2,96  gU/cm 3 .  Nilai  reaktivitas  lebih  teras  di  akhir  siklus  dipertahankan  dalam kondisi panas xenon setimbang >2,0 % D k/k dan fraksi bakar di akhir siklus dibatasi < 70 % . Pada Tabel 4 dapat dilihat bahwa untuk mencapai teras penuh berbahan bakar silisida 3,55 gU/cm 3 dapat dilakukan  sampai  langkah  ke  delapan  tanpa  diperlukan  penambahan  BKP,  karena  semua  batasan yang disyaratkan terpenuhi.
Tabel 5. Perhitungan konversi tidak langsung  II teras campuran silisida 3,55 g U/cm 3 – 4,8 g U/cm 3 pada daya 30 MW ± ² ³´ µ ² ¶ · ¸ ¹ º ² » ² ³ ¼ ² ¶² ³ ¼ ² µ ²½ ¾ ¿ À Á Â ´Ã Ä Å ¹ 3 Æ EBS/EBK) Ç ¸ ² µ» ¿ È ¿ » ²É Ê ¸ ¼ ¿ ¶ Ë Ì Í Æ Î D k/k) Ï ½ ² µÉ ¿ ¼ ² µ ²½ ¼ º ² ³´ ¹ ² µÉ ¿ ¹ º¹ Æ Î D k/k) · ² ³ Ð ² ³´ É ¿ µÊ ºÉ Æ ¶²½ ¿ Ñ Ò²½ ´¿ ³ Ó ² Ô ²¹ ¹ ¿ ³ ¿ ¹ º¹ Æ Î D k/k) Õ¸ » ¸ ½ ² ³´ ² ³ Â Ö × Ø ×Ù Ú Ú Ø Û Ü Ý Ü Ø Ü Ü Ö 3.82 ¾ ¸ ¹ º ² ¾ ¿ 3,55 gU/cm 3 × Û Ä Þ × Ø Ú ß Ú × Ø × × ÀÜ Ø Ü Ü Ö 3,90 Þ Ü Þ ÜÄ ß × Ø À× Ú × Ø ÙÝ Ý Ü Ø Ü Ü Ö 3,68 ÞÞ Þ ÛÄ Ý ÞÜ Ø À À Ú × Ø ß À Ý Ü Ø Ü Ü Ö 1,99 Þß ß ÜÄ À Þ Þ Ø ÀÚ Ú × Ø Â × Ý Ý Ø Ý Ý Ö 1,27 Þ Ý ß ÛÄ Û Þ ß Ø Ü Û Ö Ö Ö 0,37 à ¿ Ô ² µ ²¹ ² ³ Þ ÝÖ ² ß ÛÄ Û Þ Þ Ø ß À Ú × Ø × Ù Ý Û Ø Ü Ü ÖßØ ÜÜ ß Ë Õ ·
áâ âã ä åää æçå èé ã ê ë ê ì í åè 2/AU2/P2MI-LIPI/04/2012 îï ð ñ ð ò ó ì ô ðõ ö ÷ õ ì ó øù ú ðñ ù ñ û ø ÷ ü ì ù ô çèäý þ ï ð ÿ ð M ó ë ó ÿ ê ÿ  ó ì ñ ù ó ìðñ â -GAS Berbahan ......... î L ù ô â ö üð ì ôùÿ ð )
3 6 6 3 -  0 3 6 6 - - - 0,22 T - 3 3 6 6 0 3 3 33 -  0 4 0 3 -  0 0 6 3 3 - 30 0 0 - 6 0 6 0 6 -  00 4 0 6 6 3 33 - 0 0 6 6 - 6
Tabel  5  menunjukkan  hasil  perhitungan  teras  campuran  tidak  langsung  tahap  kedua  yang dimulai dengan perhitungan teras berisi bahan bakar silisida 3,55 gU/cm 3 pada langkah ke 8. Nilai reaktivitas lebih awal siklus akan bertambah tinggi sebanding dengan bertambahnya massa bahan bakar  dalam  teras  sehingga  untuk  menurunkan  reaktivitas  lebih  agar  tidak  terlalu  tinggi,  maka panjang siklus dapat diperpanjang dengan mempertahankan fraksi bakar buang maksimum <70%. Margin padam pada langkah 13  kecil  yaitu -0,37 % D k/k  (< -0,5 % D k/k) yang berarti reaktor tidak dapat  subkritis,  dikarenakan  jumlah  bahan  bakar  silisida  4,8  gU/cc  yang  semakin  banyak.  Oleh karena itu diperlukan penambahan BKP yang akan menaikkan kembali nilai margin padam reaktor.
Pada langkah ke 13a seperti yang ditunjukkan pada Gambar 12, setelah penambahan 2 buah BKP, margin padam naik menjadi  -0,8 % D k/k. Penggunaan 2 buah BKP ini dapat dipertahankan sampai langkah ke 14. Pada langkah ke 15 margin padam tidak terpenuhi, sehingga perlu ditambah 2 buah BKP lagi menjadi 4 buah BKP seperti pada langkah 13 a, konfigurasi seperti pada Gambar 11 . Pada langkah ke 15a, nilai margin padam naik menjadi  -2,28 % D k/k yang berarti memenuhi batasan desain. Penggunaan 4 buah BKP ini dapat menghasilkan margin padam yang cukup sampai teras penuh dimana seluruh teras telah terisi bahan bakar silisida 4,8 gU/cm 3 .  Penggunaan 4 buah BKP dianggap kurang ekonomis, sehingga jumlah BKP perlu dikurangi dengan mempertahankan nilai  margin  padam  <-0,5  % D k/k.  Pada  langkah  ke  19  atau  langkah  ke  4  setelah  mencapai  teras penuh  silisida  4,8  gU/cm 3 ,  jumlah  BKP  dapat  dikurangi  menjadi  2  buah  BKP  dan  dapat menghasilkan margin padam sebesar -0,88 % D k/k, sehingga langkah ini dapat dilanjutkan dengan penggunaan 2 buah BKP dan dianggap teras silisida 4,8 gU/cm 3 sudah setimbang. Nilai reaktivitas lebih awal siklus setiap teras campuran serta margin padam minimum disajikan dalam Gambar 13.
Gambar 13. Reaktivitas dan margin padam minimum teras campuran
J ! " # $ % V !% & 5 No.3 Oktober 2013, Hal. 137-149
(' ) * " % $ + " ,. ) . )/, + 2 " . % 78 & 5 9 I:: # & ; && – 240X Nomor : 402/AU2/P2MI-LIPI/04/2012
Dari grafik yang disajikan pada Gambar 13 dapat dilihat bahwa margin padam pada Teras 12 rendah  yaitu  0-1,27  % D k/k  dan  menurut  data  pada  Tabel  13,  pada  langkah  ke  13,    nilai  margin padam pada langkah tersebut tidak memenuhi batasan desain, maka perlu ditambahkan 2 buah BKP sehingga dapat menaikkan margin padam menjadi -2,00 % D k/k (langkah 13a). Langkah yang tidak memenuhi batasan desain berikutnya pada langkah ke 15, kemudian ditambahkan 2 buah BKP lagi sehingga terdapat 4 buah BKP, maka pada langkah 15a, margin padam terpenuhi.
Selain  nilai  reaktivitas,  dari  hasil  perhitungan  juga  dapat  diketahui  besarnya  faktor  puncak daya, fluks neutron cepat dan neutron termal. Nilai fluks neutron cepat, neutron termal dan faktor puncak daya teras campuran disajikan pada Gambar 14,
G? @ A? B C 4. Faktor puncak daya, fluks neutron cepat dan fluks neutron termal teras campuran
Nilai fluks neutron, dan faktor puncak daya yang disajikan pada Gambar 14 merupakan nilai fluks maksimum radial di teras reaktor.Terlihat bahwa semakin besar massa bahan bakar di dalam teras reaktor, maka fluks neutron termal akan mengalami penurunan. Berdasarkan persamaan (2 - 6)  yang  menunjukkan  bila  daya  reaktor  (P)  dan  volume  teras  (V)  konstan,  sedangkan  kerapatan bahan bakar bertambah, maka tampang lintang makroskopik (∑ f ) menjadi lebih besar, maka fluks neutron termal akan menjadi lebih kecil. sementara fluks neutron cepat akan mengalami kenaikan akibat  dari  spektrum  neutron  semakin  keras  pada  tingkat  muat  uranium  yang  lebih  tinggi.  Nilai faktor puncak daya yang tertinggi pada posisi teras tertentu terletak pada Teras ke 9. dimana pada teras  ini  kondisi  teras  belum  merata  dikarenakan  bahan  bakar  dengan  muatan  yang  lebih  tinggi belum banyak yang masuk. Hasil analisis di atas menunjukkan bahwa konversi teras silisida 2,96 gU/cm 3 menuju  teras  setimbang  silisida  4,8  gU/cm 3 dapat  dilakukan  melalui  konversi  teras  tidak langsung.

## KESIMPULAN
Konversi  bahan  bakar    reaktor  RSG-GAS  ke  tingkat  muat  tinggi  4,8  gU/cm 3 yang  dapat dilakukan    adalah  melalui  teras  campuran  tidak  langsung  silisida  2,96  gU/cm 3 -3,55  gU/cm 3 dan teras  campuran  silisida  3,55  gU/cm 3 -  4,8  gU/cm 3 .  Pencapaian  teras  setimbang  melalui  teras campuran  ini  diperlukan  4  buah  batang  kendali  pengaman  (BKP)  sehingga  kriteria  keselamatan dapat terpenuhi. Penggunaan 4 buah BKP pada teras setimbang silisida 4,8 gU/cm 3 dapat dikurangi menjadi  2  buah  BKP  setelah  teras  setimbang  beroperasi  4  siklus  berikutnya.  Keuntungan  utama dari penggunaan bahan bakar silisida tingkat muat tinggi 4,8 gU/cm 3 dibanding bahan bakar silisida
DE F 0.7 0.9 1.1 1.3 1.5 1.7 1.9 0.50 1.00 1.50 2.00 2.50 3.00 H K N O P Q R S U W K H K K K N K O KP K Q K R KS KU KW N H N K XY Z f[ \] ^_ `ab[ f [ \ c d e _ g b ^ l h i j k h m n o p h q r h s t u p h i Fl u k s n e u tro n ( 1 0 1 4 n /c m 2 . s -1 v w x y z { | } ~   x y  x  x
            2/AU2/P2MI-LIPI/04/2012                              ¡  ¢    ¡ £  ¡ ¤     ¥  ¦  § -GAS Berbahan .........  ¨   ©     ¡  )
tingkat  muat  rendah  2,96  gU/cm 3 pada  teras  RSG-GAS  ialah  panjang  siklus  operasi  dapat diperpanjang dari 22 hari menjadi 40 hari sehingga dapat menghemat penggunaan bahan bakar.

## DAFTAR PUSTAKA
1. Anonim. Laporan  Analisis  Keselamatan  RSG-GAS  bab  V-Reaktor , Revisi  10, BATAN. 2008. p: V41-46
2. Surian Pinem, Tukiran Surbakti, Jati Susilo. Pengukuran parameter kinetik teras silisida dengan  kerapatan  2,96  gU/cc  reaktor  RSG-GAS.    Jurnal  Fisika  Himpunan  Fisika Indonesia 2004;A4: 021801 - 06
3. Lily  Suparlina  dan  Jati  Susilo.  Evaluasi  Manajemen  Teras  Reaktor  RSG-GAS Berdasarkan Perhitungan Fraksi Bakar. Prosiding PPI-PDIPTN; 2009 p: 131-136
4. Lily  Suparlina.  Analisis  Peningkatan  Fraksi  Bakar  Buang  Untuk  Efisiensi  Penggunaan Bahan Bakar U 3 Si 2 -Al 2,96 gU/cc di Teras RSG-GAS. Prosiding PPI-PDIPTN. 2010 p: 176-182
5. Liem  P.H,  Bakri  Arbie,  T.M.  Sembiring,  P.  Prayoto,  R.  Nabbi.  Fuel  management strategy  for  the  new  equilibrium  silicide  core  design  of  RSG GAS(MPR-30). ­®¯ ° ± ² ³
Engineering and Design. 1998;
´µ ¶ (3): 207 – 219
6. Jati  Susilo,  Analisis  pola  pemuatan  bahan  bakar  teras  setimbang  RSG-GAS  berbahan bakar U3Si-Al kerapatan 3,55 gU/cc, Jurnal Teknologi Reaktor Nuklir Tri Dasa Mega. 2007; 9(2):75-85
7. Lily  Suparlina,  T.M.  Sembiring.  Application  of  Neutron  Diffusion  Theory  for  High Uranium  Density  Fueled  Reactor  Core  Calculation,  Asian  Physics  Symposium Proceedings.  2005; p: 332-337
8. Lily  Suparlina.  Analysis  of  hot  spot  factor  for  Silicide  4.8  gU/Cc  RSG-GAS  Core. Proceedings  of  The  International  Conference  on  Mathematics  and  Natural  Sciences (ICMNS). 2006; p:1065-68
9. Surian Pinem, Tagor MS, Setiyanto. Analisis transient teras reaktor RSG-GAS pada saat laju  alir  pendingin  turun  dengan  menggunakan  program  MTR-DYN.  Jurnal  Teknologi Reaktor Nuklir Tri Dasa Mega. 2009; 11(3):153-161
10. Setiyanto, T.M. Sembiring, S. Pinem, L. Suparlina. Neutronic analysis of the RSG-GAS silicide  core  with  uranium  density  of  4.8  g/cc . IAEA  Publication  International Conference  on  research  Reactors  Safe  Management  and  Effective  Utilization,  IAEA- CN-156, STI/PUB/1360 ISBN 978-92-0-160108-7 ISSN 1991-2374 Vienna; 2008.
11. Peng Hong Liem and, Tagor Malem Sembiring. Design of transition cores of RSG GAS (MPR-30) with higher loading silicide fuel.  Journal of Nuclear Engineering and Design. 2010 6: 1433-1442
12. Ahmad  Lashkari,  Hossein  Khalafi,  S.  Mohammad  Mirvakili,  Shokufe  Forughi. Neutronic  analysis  for  Tehran  Research  Reactor  mixed-core.  Progress  in  Nuclear Energy. 2012; 60: 31-37.
13. M.  Keyvani,  M.  Arkani,  A.  Hossni  Rokh.    Optimization  of  in-core  fuel  management strategy  of  Tehran  Research  Reactor  (TRR)  using  MCNP-4C.  Annals  of  Nuclear Energy. 2010; 37: 1683 – 1687
14. IAEA. WIMS-D Library Updat e. IAEA, Vienna, 2007
15. Liem  P.H.  Batan-FUEL:  A  general  in-core  fuel  management  code . Journal  Atom Indonesia, July 1996 22;2:67-80