Ke Daftar Isi
PRO~RA~
HECTIC
II UNTUJ: ANALISIS
PADA REA~TOR BERPENDINGIN
PERPINDAHAN
PANAS
GAS ATAU LIQUID
Tri Wulan Tjiptono Bambang Soemarsono Pusat Penelitian Nuklir Yogyakarta ABSTRAK Suatu paket program komputer untuk Analisis Perpindahan Panas Reaktor Berpendingin Cas atau Liquid dengan Program HECTIC II yang aslinya untuk komputer IBM 7090, telah dapat diproses pada komputer IBM PC dengan kompailer WATFOR 77. Latar belakang dari pengalihan program tersebut adalah dalam rangka studi termohidrolik reaktor. Untuk dapat dikompilasi dengan WATFOR 77 pad a komputer IBM PC program asli tersebut harus disesuaikan sedemikian rupa sehingga pernyataan-pernyataan dari program baru tersebut dapat diterima dan dimengerti oleh komputer IBM PC. Setelah program tersebut berhasil dikompilasi, untuk uji cobanya digunakan untuk memproses data masukan program aslinya. Dari output program baru tersebut kemudian dibandingkan dengan output program aslinya dan ternyata sarna. Selainitu, juga telah dicoba untuk perhitungan reaktor PSG dengan hasil yang mendekati dengan kenyataan. ABSTRACT The package programme for Heat Transfer Analysis for Gas or Liquid Cooled Reactor which is sappased to be with IBM 7090 has been used with IBM PC Computer and compiled by Micrc Soft Watfor 77. This program can be used for studying thermohydrolics. The actual program must be modified in order that the statements in the new program can be recived ter and compiled by MIcro Soft WATFOR77. This program can help for studying thermohydrolik. The actual program must be modified in order the statemens in the new program can be received and recognised by IBM PC and WATFOP77 compiler. After occomphithing the compilation the actual input data was processed as a tryout. The output of the new program was then compared to that of the actual program. It proves that the new program and the actual one gave similar reluts. This program is also used for the calculation of PSG and the results are class to reality.
506
Jawaban
a. Spesifikasi
kelongsong
b. Sebenarnya
lebih baik bahan kelongsong
yang
dimuatkan
sarna, agar mudah
diperoleh
dalam teras, perlakuanya.
dari General
terdiri
Atomic
(bahan
dari jenis
bakar) yang
507
I. PENDAHULUAN
Analisis nya
dapat
maupun
dilakukan
pertama dengan
penulis judul"
bangunnya.
REACTOR
mencoba
tersebut
maka sebagai
perpindahan
ANALYSIS
panas
program
COMPUTER
tersebut
panas dan
dalam
LIQUID
untuk menghitung
perpindahan
langkah
suatu paket
OF GAS OR
••. Paket program
laju
keperluan
mentransfer
yang dapat digunakan
aliran,
reaktor,
II an IBM 7090 FORTRAN
TRANSFER
PASSAGES
program
masalah
nuklir umum-
keselamatan
Untuk
dengan komputer,
HECTIC
HEAT
FOR
panas suatu reaktor
dengan analisis
untuk mempelajari
reaktor,
laju
dikaitkan
dalam rancang
biasanya
GRAM
perpindahan
PROCOOLED
adalah rugi
suatu
tekanan,
temperatur
pada
elemen bakar dari jenis reaktor
ber-
pemindah
panas maupun
pendingin
gas atau liquid dengan "single
phase"
(~atu
fase)
"much number" rendah, aliran pendingin axial. Analisis· berdasar suatu geometri minimum "flow tube" simetri dari elemen bakar
reaktor
"plates"), suatu
dimana setiap permukaan
dan
alur
"minimum
diperoleh
Batasan-batasan dalam "pin"
penyiapan (batang)
dibagi sedemikian
sector"
M permukaan
suatu perangkat
yang
diberlakukan
model
"lumping",
konfigurasi
beberapa
bakar).
Panas yang dibangkitkan
bagian
sikan pada kelongsongnya. alur
besar besar,
untuk
kebutuhan,
bakar
akan
pada
program
berbentuk
elemen bakar dibuat
sek-
secara sebarang permukaan
dibagi elemen
oleh elemen bakar disimula-
Program lamp"
dengan komputer
itu statemen
ini digu-
(aliran maupun
dengan 24 "surface
lagi asalkan
bahan
pada
dan simetri,
Sektor tersebut
menjadi
30
rupa sehingga
input untuk elemen 'bakar
digunakan
tor simetri minimal.
maupun
dan N alur.
data
nakan untuk mengurang
panas ("rods"
dimensi
ini dapat digunakan dan
masih
dapat
yang memorinya
untuk lebih
juga lebih
harus disesuaikan
dengan
508
Paket dengan
program
WATFOR77
Kbyte. Program problem"
d1 atas d1mod1f1kas1
pada komputer
faatkan
asli.
terdiri
Sehingga
tentang
bilangan
diameter
Reynold,
diharapkan
laju aliran,
(dalam hal ini ada 5 alur),
(arah
lokal
aksial)
sebagai
Out
program
dan hasil perhi-
laju aliran rerata, parameter
untuk
heat flux lokal
pada tiap permukaan
"sample
program dapat diman-
reaktor.
hidrolik,
640
sesuai dengan output
dari listing data masukan
tungan-perhitungan
temperatur
memori
baru telah dicoba dengan data masukan
untuk studi thermo hidrolik
tersebut
d1komp1las1
dengan
dari program asli dan hasilnya
dari program
alur
IBM PC,
dan
sebagai
tiap
absolut,
fungsi
jarak
dan suhu lokal pending in pada tiap alur
fungsi jarak (arah aksial),
let, panas yang hilang kelingkungan, error suhu pendingin
juga
serta suhu campuran rugi tekanan,
out
perk iraan
dan perk iraan "mach number".
II. TEORI
11.1.
Perhitungan Dalam
menjadi
Aliran Fluida
analisis
N alur (passage),
alur identik. percepatan, analisis
dan •• turbulent
aliran
•• turbulent
sehingga
Efek gesekan dinding,
alur semua ditinjau
Untuk
dengan HECTIC
"boundary
eddy momentum
dengan parameter-parameter ini berarti
antara •• turbulent
shear dengan perkalian
tekanan
Dengan
Baussineg
rerata.
telah membuat
ngan memasukkan
analogi rumusan
"turbulent
"spacer",
exchange"
Persoalan
dibagi
pada semua
drop terhadap
adalah mendapatkan
layers"
total
rugi tekanan
dalam perhitungan.
turbulen
stress"
II aliran
diantara
dasar dalam
hubungan
antara
aliran rerata.
menentukan
hubungan
dari kecepatan
terhadap
dan
"Hukum
Newton",
tentang •• turbulent
shear"de-
eddy diffusivity"
sehingga
509
(1)
Em : turbulent Dalam
program
eddy diffusivity
ini "eddy diffusivity"
disajikan
dalam
yang sama seperti
sering digunakan
dalam perhitungan
parameter"
konduksi
padat.
Stress"
dari
dari arus diantara
banding
dengan
berbanding
selisih
terbalik
permukaan) •
Data
pada tiap alur, ter
hydrolik
dihitung
benda
dua susunan
kecepatan
tabung dianggap
dengan jarak diantara
ini
program didasarkan
permukaan~
alur (tegak
Shear sedan ~urus
"Eddy Diff·usivity·· dihitung
pada bilangan
dari masing-masing
··lamped
"Turbulent
dan luas
cara
Reynold
alur.
Eddy
dan diameDiffusivity
dengan persamaan Re =
~m
20
102
f = faktor gesekan Re = bilangan Adanya
efek
deffusivity universal
Reynold
turbulent
maka
dan dinyatakan
didefinisikan
dengan persamaan
Re = v v
Faktor
20 =
viskositas
gesekan
f
=
CF
f
kinematik
dihitung
(Re)-EFYF
dengan persamaan
pula
turbulent
provil kecepatan
510
Harga
CF,
suaikan
EF dan YF
dengan
adalah tetapan
persoalan
fluida pending in).
yang hanya dapat dise-
yang akan dianalisis
Jika pada alur ada "spacer"
(tergantung atau "pro-
tensions" (pengatur jarak atau benjolan) akan menarobah adanya tahanan (harobatan) gesek yang disebabkan oleh "shear stress" sepanjang
(CDS)
dan
(harganya
alur dan dinyatakan
besarnya antara
tergantung
sebagai drag koefisien
dari
bilangan
Reynold
0.2 hingga 5.0).
Rugi tekanan pada tiap fluida dalaro dari analisis momentum pada
alur ditentukan .,flow tube". Untuk keperluan tersebut maka ditinjau suatu alur k. Atas dasar gambar tsb di atas dibuat neraca gaya yang bekerja pada fluida dalam alur
Wall shear persaroaan
L
force dihubungkan
dinyatakan
faktor
sebagai
gesekan
oleh
=
wk
2 Inter passage
L.1k
dengan
k
gc
shear diarobil dari
= ~turb. ~
~i.k.=
=
.1k
jarak an tara "nominal centroids" dari slur dan k tegak lurus antara keduanya
i
[I
511
Spacer drag force adalah
F sk
Diameter
A f"sk
=
.
hidroulik
dinyatakan
sebagai
4 AOk
bk
Juga didefinisikan
Eddy difusivity
faktor geometri
rerata
campuran
dinyatakan
=
2 Sehingga 1:4
fungsi tekanan dapat dinyatakan
Pk =
(P l.n - Pout)
-tk- v.J J
V
-.
2
~YEDM
Ck
C1.k
2 go +-favVk2 . 4 fk -
Lk AOk Dk
t
+ -
Lk
l
AOk
V
C de
out
l.n
+
2( -
A£sk
sebagai
-
Fav fav
512
II.2 Perhitungan
pembangkit
Panas pembangkit plat
panas
persatuan
atau pin dihitung
panjang pada tiap permukaan
berdasarkan
asumsi bahwa
distribusi
fluks panas aksial adalah identik untuk ti~p permukaan dan distribusi fluks panas melintang tak tergantung dengan X, jarak dari inlet. Distribusi fluks panas aksial dinormalisasi,
dan dinyatakan
dengan
NF.
Integral
pembangkit
panas harus sarna dengan satu.
pembanding
dari
(pada
seluruh
rerata
daya yang dibangkitkan permukaan)
yang dibangkitkan
mukaan),
dimana menyatakan
silkan pada permukaan persatuan adalah
KW
panjang
di NF sepanjang
NF
persatuan
pada suatu titik X persatuan
juga
panjang
dengan
fraksi dari daya total yang diha-
tertentu.Daya
pada permukaan
distribusi
total yang dibangkitkan ke j
pada suatu
fluks aksial harus memenuhi
(untuk tiap permukaan) o
PF~
harus memenuhi
juga
M J-;/
PF j
= 1
day a
(pada semua per-
1
2"
panjang
titik
= menyatakan input daya total dalarn Kilowatt seluruh elemen.
Normalisasi
Sedangkan
sebagai
(untuk tiap persoalan)
X
untuk
51.3
11.3 Perhitungan Proses
perpindahan
perpindahan
ini adalah perpindahan ra permukaan kaan
dan
panas
panas yang dihitung
panas secara konveksi,
bahan bakar ke pendingin,
panas
dalam
radiasi,
konduksi
yang hilang dari permukaan
(dalam hal ini adalah dari calandria Untuk perhitungan
perpindahan
radiasi
an tara permukaan
fisien
perpindahan
program
antar permu-
ke
lingkungan
tube misa1nya).
panas secara konveksi
bahan bakar dengan pendingin,
panas konveksi
panas radiasi masing-masing
dihitung
anta-
dan
faktor
dengan
St( i)= CE*( Re( i ))-EH*YST
dan koe-
perpindahan
:
i =1,2 •...m
di mana
f
St( i)= h( i )/ Cf VLi menyatakan
alur ke L, dari harga h(i) di atas maka perpin-
dahan panas secara dihitung
konveksi
dari permukaan
ke j ke a1ur ke i
dengan
REFj~=
+ -+ P~ *EM~ adalah dengan
faktor 1,
perpindahan
perpindahan
masing-masing panas radiasi
panas radiasi dari
dengan emisivitas antara permukaan
dengan Q'rQd
J~=
V REFJ~(TJ4_T~4)
permukaan
j
EMj dan EM~. Laju
j
dan 1
dihitung
514
Laju dan
perpindahan 1
panas secara konduksi
antara perroukaan
(ingat bahwa dalaro hal ini satu bahan
roenjadi beberapa konduksi
bakar
j
dibagi
jadi yang diroaksud di sini adalah
segroen ,
antar segroen dalaro satu bahan bakar) dihitung
de-
ngan
Kj~
adalah konduktivitas
dapat
panas antara segroen j dan
diarobil saroa dengan nol bila tak ada proses
antara
segroen j dan 1.
akhir
Proses perpindahan
1,
dan
konduksi
panas yang
ter-
dihitung
perroukaan dria
adalah panas yang dipindahkan dari suatu dalaro hal ini terutaroa dari calankelingkungan,
tube
digunakan
roisalnya ke daerah untuk keperluan
GAj adalah konduktivitas adalah
rooderator.
tersebut
Persaroaan
yang
adalah
panas segroen ke j, di roana segroen j
segroen yang berhubungan
dengan
lingkungan
(bukan
perobangkit panas)
11.4 Perhi~ungan
temperatur
permukaan,
pendingin
dan panas
yang hilang ke lingkungan Suhu pada perroukaan j neraca panas pada perroukaan
q.j = q·rad j
q·conv
dihitung j
berdasarkan
persaroaan
yaitu
j~ + q·con j~ + q·rad j + q·am j
dapat ditulis sebagian ~
q·rad j = VREFj1(Pj(Tj2+T~2)(Tj+T~»(tj-t~)
•••••••••••
(1 )
515
Dari
persamaan
persamaan
(1)
untuk m>1 maka
linier dalam tj,
menyelesaikan
maka
akan
diperoleh
m
laju tj ditentukan
set
dengan Suhu
persamaan
tersebut secara simultan. berdarasarkan neraca panas pada dihitung pada tiap alur, dimana neraca panas pada alur i alur yang bersangkutan, adalah dt.1. =
M
L
h.1.b.1.j(tj - t.1.)+
Jzt
dx IV YEDH • CP)Jv
L
*-=-1
2 parameter
jarak nominal
di alur i yang kontak dengan alur k
diantara
alur i dan k,ke 2 permukaan
Seperti pada suhu permukaan,
maka suhu pendingin
dengan menyelesaikan
i
Dalam sebagai dihitung
sejumlah
dua hal tersebut
fungsi jarak (aksial). berdasarkan
dari permukaan pipa terluar)
dinding
terluar
ke lingkungan
Harga integral
dari persamaan
diperoleh simultan.
diatas mesing-masing
dihitung
Panas yang--hilang kepermukaan secara konduksi
(dalam hal ini adalah dinding
dan dinyatakan
GA,j(tj - ta) dx
Simpson.
secara
panas yang dipindahkan
=
q'am1:>
persamaan
tsb.
sebagai
diatas dihitung
:
dengan aturan
516
Hal
terakhir
ngin keluaran
yang dihitung
adalah suhu rerata
dalam hal ini dihitung
pendi-
dengan
tV ~ +-
"'mm
-
1..-::.(
-
III. DISKRIPSI
INPUT/OUTPUT
Data masukan Data-data program
ditulis
utamanya,
program
ini adalah sebagai
dalam berkas masukan berkas tersebut
yang
berikut. terpisah
dari
adalah
1. HECTIC1 .INP,yang berisi data: KTR
dengan FORMAT
12 adalah kontrol
parameter,
yang berhu-
bungan dengan beberapa pilihan PILIHAN/ARTI NILAI o
baca semua data masukan
1
baca saja EMJ GAJ PFJ suhu ternormalisasi fluk masukdan data PFJ,GAJdan saja
EMJ
2 baca
baca
3 4
5 an gas PROB dengan
format F10.6 dapat diisi tanggal
N
dengan
M
menyatakan dengan format 110 bahan bakar maksimum 24
ITO
dengan
format
format
110
110
manyatakan
menyatakan
penggunaan
banyak alur maksimum banyaknya inlet
30
segmen permukaan
temperatur
opsion
nilainya: 1
jika suhu masukan
diberikan
pada data matrik
o
jika suhu masukan
diset sama dengan TIN
TIN I
517
2
jika
suhu masukan
pada persoalan dua persoalan NFQ
dengan
format
diset sama dengan
sebeluronya,
suhu
keluaran
aslkan N nya saroa antara
tsb.
I10
menyatakan
tentang
normalisasi
fluk
yang dipilih o
jika
distribusi
fluk aksial pada
semua
saroa, jadi hanya 21 yang diperlukan 1
NOPT
dengan dicetak
1111110000
format 10I2
untuk menyatakan
dengan
CF EH
format F10.4
dengan
F10.6
molekul
faktor
format f10.6 koefisien
gesekan
0.065
dengan
format F10.6 faktor pangkat bilangan
Stanton
<nominal
dengan dengan
harganya
lb/lb-mol.
format F10.6 konstanta jenis pendingin
bilangan 0.2 faktor
dibilanganReynold
0.2 pada persamaan
bilangan
= 0.0291)
pada persaroaan kapasitas
jenis pendingin
format F10.6 density
persaroaan
0.023xPr-2/S
format F10.6 konstanta
pendingin
pada
stanton
format F10.6 koefisien
DLIQ dengan
akan
dicetak
faktor gesekan besarnya
dengan
tas tergantung
gas
pang kat
dengan
panas tergantung
CAYW
berat
pada persamaan
Bs,B~ dengan
yang
berikut:
Reynold
B1,B2 dengan
TIN
NF
cair)
format
pada persamaan CH
data
agar semua data dan hasil perhitungan
untuk pendingin EF
pada matrik
jika distribusi fluk pada semua permukaan berbeda, maka ada M x 21 nilai atau NF adalah matrik M x 21
Dalaro hal ini dpilih data sebagai GMW
permukaan
reaktor
yang digunakan
pada persamaan reaktor
viskosi-
yang digunakan
lb/ftS <0 untuk
cairan,
gas)
format F10.5 suhu masukaro rerata oC atau format F10.5 total daya yang dibangkitkan
semua permukaan
Kw
W
dengan
format F10.5 aliran pendingin
P
dengan
format F10.5 tekanan masuk,
total
psi
lb/hr
of pada
518
ELF
inc. bakar of inc. format format F10.5 F10.5 suhu koefisien bahan diffusivity gesekan spacer faktor eddy perubah drag untuk dengan yang panas dengan lingkungan panjang gesekan
untuk momentum YST
dengan
format F10.5
YEDH dengan
faktor perubah
format F10.5
bilangan
faktor perubah
eddy
stanton diffusivity
untuk panas ERMAX dengan
format F10.5 maksimum
ENFJK dengan
format 7F7.4 menyatakan
sebagai TINI
dengan
N
dalam
fluk neutron
dibangkitkan
ke j.Yaitu
ke j, j = 1,2,3,
.••••.•. M
data konduktivitas
antara permukaan
nya, dalam setiap M BTU/hr-ftoF, o jika permukaan
atau oC
fraksi dari daya total yang
pada permukaan panjang
untuk
data faktor daya dari
format F7.4 menyatakan
persatuan
berbeda)
data suhu masukan
of
format F7.4 menyatakan
plat/permukaan GAJ dengan
mempunyai
format 5F9.5 menyatakan
alur 1,2,3, PFJ dengan
data fluk termalisasi
jika NFQ = 1, maka datanya matrik
fungsi x/Lh,
Mx21 (tiap permukaan
error yang diijinkan
tidak kontak
dengan
panas
lingkungan-
dalam hal ini harganya
langsung
dengan
lingkung-
annya EMJ
dengan
format F7.4 menyatakan
permukaan 1 ,2,3, untuk pendingin cair AI
dengan
format
1,2,3 ••..•.• M, CAYL
F8.4
M,
menyatakan
emisivitas
format 10F7.5
antara
permukaan
harga sama dengan data
menyatakan
j dan
1 dimana
data ke-2
lansung,
data dalam matrik
dari
luas
0
alur
dalam inc.kwadrat
dengan kontak
••~ •.•
data
MxM, dalam BTU/hr-ftoF.
konduktivitas permukaan
ada
519
AFSl dengan
format F8.5 menyatakan
data luasan frontal dari
spacer pada alur 1,2,3 ....•.• ada K
spacer dengan
pada alur tersebut FJL dengan
menyatakan
permukaan ClK
k.A,
dengan
j
jika pada suatu alur
luas frontal A maka harga
data dari faktor pandang
dan 1, berupa matriks
fraksi dari radiasi
j dan diterima
permukaan
format 5F7.4 menyatakan
faktorn yang merupakan
yang
Data
perimeter
ditulis
datanya
OUTPUT
PROGRAM
data nmising
perbandingan
10F8.5
aliran
masa
geometry
perimeter
dialur
jarak nominal
dari alur i yang kontak dengan input HECTIC3.INP ,berupa matrik
Output dari program Laju
meninggalkan
ini adalah
pendingin ekivalen
dinyatakan
j
dengan BlJ
NxM, dalam inc.
g
alur
1 ,2, .•••• ~••N
termohidrolik
- Bilangan
Reynold
alur 1,2, .•...... N
- Bilangan
Stanton
a1ur 1,2, ••••.... N
alur 1,2, •••••..• N
gesekan Fanning
- Koefisien
dian-
permukaan
- Diameter
- Faktor
i
dan saling tegak lurus
pada berkas
Format
MxM.
1 .
yang kontak dengan alur k terhadap tara keduanya
AFSl
dalam inc.2
format 10F7.4 menyatakan
antara permukaan FJL
N,
alur 1,2, •..••..• N
PP konveksi
alur 1,2, ..••.•.. N
Harga - Panas yang dibangkitkan
permukaan
1,2, ........•• M
- Suhu Permukaan
1,2 .•••• M,
- Suhu pendingin
alur 1,2, ••.••. N,
Selain
tersebut
di atas mas~h
campuran
rerata
besaran
lain yaitu
sebagai
fungsi X/Lh dilengkapi
:
- Suhu keluaran
- Panas yang hilang ke lingkuan - Rugi tekanan Perkiraan Mach number
bahan bakar
besaran
520
- Jumlah nodal yang digunakan - Error untuk perhitungan Dimensi
pada perhitungan
suhu pendingin
dari masing-masing
besaran
di atas,
dimensi dari data masukannya. Dimensi disesuaikan dengan persoalan yang ada.
III.
adalah suatu program but dapat diproses beberapa
mengikuti
data masukan
BAB.I, bahwa program
untuk komputer
dapat
perubahan WATFOR77,
tersendiri
tidak disajikan
Pada dari
dasarnya
program
pad a makalah
tersebut
untuk
komputer
IBM PC 640 KB.
komputer
IBM PC program
Program
baru
program
tersebut
Problem" nuklir
yang Oleh
dengan
karena
(FILE) yaitu berkas programnya,
pacta INPUT
masukan perubahan
baik FORMAT bacamaupaun
telah dicoba
untuk
aslinya,
dan kemudian asli,
samadengan
memproses outputnya
dan ternyata
output program
berpending
in gas nitrogen
batang elemen bakar yang disusun
Output program
prinsipial,
dikompilasi
disajikan
II.
dibanding-
aslinya.
"Sample
bahan
dan terdiri
dalam kisi
pada Lampiran
"sample
output program
terse but adalah data suatu perangkat Reaktor
pada disket
pada FORMAT,
kan dengan output program baru
dengan
pada disket sehingga
tidak pada TAPE. Tentang
dilakuakan terutama FORMAT tulis.
problem"
dapat
ditulis
OUTPUT ditulis pada berkas
dikompi-
ini.
tidak ada perubahan
pada
dilakukan
yang kemidian
baru disajikan
II terse-
terse but dikompilasi
dan program
WATFOR77
dan keluaran,
IBM PC, asalkan
dan penyesuaian,
mikrosoft
HECTIC
IBM 7090. Program
dengan komputer
lasi lagi. Dalam hal ini program
19
akan
HASIL DAN PEMBAHASAN Telah diuraikan'pada
dan
suhu aksial
bakar dari
heksagonal.
521
Selain dalaro hal perangkat harga
itu juga dicoba untuk analisis elemen bakar RSG, ini "minimum sector" elemen bakar adalah 1/8
dengan asumsi bahwa tiap plat sebagai
fraksi daya yang dibangkitkan
revevisi jiakan
6 Gb.4-23, pada
dilihat
demikian
Laropiran III.
adanya
kesesuain
diarobil dari buku
PSAR
data lainnya dan hasilnya
disa-
Dari hasil
dapat
perhitungan,
dengan data
yang
tidak persis saroa. Kita tinjau suhu pemukaan pada
ada,
walaupun
dan
pendingin
ujung elemen bakar dari hasil perhitungan
sesuai'
dengan data yang ada,
dengan perhitungan model
perpindahan
perhitungan
ini saroa ini
dan aeauai pula dengan
hasil
suhu pada permukaan
yang dihitung dengan Dari hasil panas pada elemen bakar plat.
terlihat
bahwa suhu permukaan
dan pada tengah-tengah
IV.
1 permukaan,
yang ada di tengah
plat yang bersangkutan.
KESIMPULAN
Suatu paket prograro dapat di proses pada jenis komputer yang
lain
dengan
asal dilakukan
Micro
untuk analisis lis mencoba kaan
penyesuaian
Soft yang tersedia. perpindahan
untuk membuat
dan dikompilasi
Dari
:lagi
hasilperhitungan
panas pada elemen bakar RSG penusuatu kesimpulan
elemen bakar dan pendingin
bahwa suhu permu~
pad a bagian-bagian
ujungnya
saroa, hal ini sesuai dengan data yang ada. ditengah-tengah menurun
aesuai
permukaan aksial, daan
mempunyai
Plat yang berada dan suhu suhu yang paling tinggi,
dengan letak plat. Dalaro perhitungan
maupun
pendingin
dan dibagi menjadi
dihitung
penyederhanaan
fungsl
21 titik tinjauan.
hasil dengan data yang ada,
adanya
sebagai
persoalan
poslsi
Adanya perbe-
hal tersebut dan dalam
ini suhu
hal
dikarenakan penyusunan
data yang diperlukan. Dilain pihak penulis dengan
program
juga menganalisis
SFRTEMP-R21
dengan hasil peritungan
di mana hasilnya
prograro HECTIC
II.
elemen bakar RSG juga
mendekati
""'--,/ I
/
~ ~ H ~ 6;.
•
t
I
H'
i
01l.b608E-01 )HI RfI ~ or FF S3.2~Y4[-03 •.2f1f7I:-G I,9SI7U~bF. I0 C') [E -- 0'~ 00'1 G 'I~:$1.1 7./1.~4E-0) I. PAS 1.L96'J7E L21261.: 3.29~5E W r0 011.5264E-01 011 1~.55IdjE .3.631 2.1549E 32.t5215E 2.7Y4?E 2.7dGIE ).56211 17.59BOE-O) 1.IIISOt 1.4070E 1.1404[ •••ST94 'I~820E S 2 I616 '2.) 37IIE [-0 (J (,2 02 002 O~ 1 •• 22 .~Y )E .:l.71~OE ~.~473[-C3 11.2fHIE-OI 1.92)6E-Oi 1.1610E 02 II. 04 2.uu06E 1 .)'~ E U 1.60/2[-O~
="'0
o s:;.
S
c-f
'd'd s:: s:; ct-c+ (1)
ti
8'
Ho
Ii
b::10t1
:s: -
IU
oe
-.:J
o
'-DP:: IU CD
•.... •.....
'd Ii o
.; III
PC!)
~
§
,
.
\J1~
Ni N
.30 .80 75 .3'1 .25 .1i5 .65 .70 .60 . .A5 .90 .'10 .95 .1.0 1S 1.00 .')5 10
I
J\)
lId I1~ 'd'l. i.~It. I.n.i 2n. nt!.
C. o. 10 O. 491. <,I 20.~2. IU76. 14 7 16. 19'J6. 651. 4 6 oGL !.J6'J. . '10'). 577. HJ2. 1!:J. 1472. 91,0. 490. 201J6. 2062 2017. II!:! 1614. 1117 1 1290. 11J1J2. 1378. I1232 II. 11,2 1628. 1977 IIJ f'~ t, 400. 1458. I 3S!:!. G. o. IJ111 <,I IJI I<,I 16 7u9. i.309. O. 2032. 2046. 701 2 2062. J124. 134. IJ IJIJO. 9 IJ72. CJ3U. <,15 1,159. e22. 751. 1d0. 3'16. cc. 7'JO. 3':14. 0'1. 10 '13. 1,0• 5U. /j II. YIJ• C 911. /jO. o. o. O. o. c. 77. bII19. fj 14 I. 7 1•'J. .••• •• ••. ••7. 1• 1769. 1272. IH81J. 4'161. HB6. '.19 127;> 1d!.JI,. .15&. H6 .1O.~ 77L 1956. 158. ~1,941J. I.')d. 1,10:,1. 501. 876. 51i• I.,. 1 .1. I• I. • 3 Ie 1442. 9 670. 0 '). 1070. 497. 1070. 260. 781. 746. 3 191 .130. 621:! 60C. 1241 251. 251:J. 121,I. 382. c. 1064. 1C 11.82. 162H. 29':1. 1769. 67H. 1232. 1002. 1:16<,1 eB6. 26'1. IJ09. U60. 47'.1. 'I 22') 5'11. I~72. B61. )11112. 92IJ. 1 1501J 430. 6 456. 822. 76 6 Ie 75. 12. O. U. Ie. .~. 3BO. 2 1427. 3'16. ItJ76. 475. 351. 51. I. 520. 217. 651. 327. 209. 174. IU. 409. 105~. 132. .25H. (1406.
[rJTU!fJi<
FTI
UPJ -
HEAT
GE"IE~ATIUI\,
\Jl \1~
0.50 O. ~) 70 0.C5 0.S5 0.65 £) 3I)07~ ('I~ . 0." 0.90 1.00 iJ.?O aO. O.no .l,O 10 0.&5 0.95
-n6. ~u. '~'J. J'. j IInt>. IIH6. 1.1n. I" I',li'll. en.
<,I 1',1a (3:) 11 353. 1:1 SQ. GO. •I.•.'. I, 11. 1000. 110)( 'I'0 ') I) ~ 91, I. I4I)!J 102 1? ,y12"'6. 03 1<). 13 12[\ 3) u6. 13. I', ')5. 16. lid 1,,15. 12~<). II,:!~ 100' 11H". 7III IIaH. 1,,32. 1376. 12 11"'1. Ie 13'" I" '0'.2. 10') tI.LB. 1.1'1 1', It' tur. U~) 119). let 610C5. In 10 971\. 1', I11(,)1. .~ 127 00<). 121)1, b1033. 1311. ') 1.30" 1',06. Ie l('fJ~. bY(,. 9 "'td. oj!:J" 'J (~II 1I102. I01b. VC 311IJ. O')b~ 31111. I, ? HI. 3'1 'If! ~ H11. 'UO. H<,I. 333. I1,,2. 1011. 152. 5116 71, (0. (5. It III. <)8. Q~. U 00. Oi. U'I. CII. I:EG-r J2tJ t1. n. It ((. j6. I17. III. 3'9. ~ 9. (). l'J. It O. 4. .1• I, ') I. •J. 7. t(7 •1 •••• 11,02. 35', 0') ~:! 3. 10 1011'1. 11115 10111. 12,,3. 1216. 1<' £Yb 12'i2. C I,.~. 11)L 11','1. I13<) I10611. 132 134d. 13". 37'\. 3I,6/,. 3',1> 2. 3/J. •I! SURrflCE I"I,? II, I106'1. ., 10(;'1. 126') 1213. 14 'it!\I, 110',. 10'16. 1.3 112 I', III lIt'l. 12('6. 11<)0. 12H. II, 1.146. IUO-l. C(I,. 10'>2. I? II, IJ/O. H(,? /1.17. U'J'J. 13(( \II IHu. 121'1',:1 O~G. I, .\1~O). 3H<). 3" ,W.,J. 1t'9 11 I', 76. 22. yC II" 6t) a'l. /)!I II 0). 5.11.. IaO. 3" .II). fJ £5b. 1'1. 71, 1375; (I, 15. It! '1 II. I: 2. S. .~ 1. 10U6. tl1OC? 15. 1.\0 1•••• S'I In2. 1011. TEI'I'[RIITLH(;S, 1O'j<). 1/, 5U 'II~• ~ <) b. I1110) 'all,
~N Vl
TSJ
-
N
'.
-
L
'769. tII~~ aj o.. ' ~!.~ I,~ I'I',
fll,,;, 180. 11025. 10117. , ?II b116(. UY5. If (,1II!. 1J~5. 9117. 81,10. I) lIe lJ71. lJGO. In' 1,111 B811. '> 812. U n1122. 115~. l~ 91169. 999. 11 () 2q800. BOO. ')3IC83. "Y. <) Y I.J L 0/1 3Ij10 51. (,1.L j'2. 70. 71 IS. 1h,/11 ')2/1 CEG-r 7I, II'0. 1,II)/I2. ') IIi SO. tJ68. 1020. Ie 1051. 9:; 27. tCO. COUll\NTT [fJ I) 1303. 1'2.96. 1051, 1LOf. 1230. 811220. t"1 826. ;~I U fJ9(;. 105'). "1060. 10lV. 1'2.'2.2. 12.36. 11112• 976. 115 a951 ell 03') I1151, 211 23. 1.3.L 19.3. 1110. 152. I)'2. 16. I) 5~)• UU\ 0L7() 2. 7.~• 3I(,1••I.•• ••••. T L I1253. 1i86. I12.' 12 ,IOP2. 1C 0996. 10"( 2(7!j C. d.1 12.~ 1061. 111 1080. 12011.
K[
S ,
VI VI
Tel
, ~.bl1 £:.05 0.15 1J.2P {\.3~ 6.55 t.OC e i.35 ;.25 :.15 3f! .98 (i.1S C.65 C.7D &.56 &.18 C.% :.80 007.17 ~.15 If I.
.~
, UI~
----1---
IIIm~.35 lU9 L99 c.~~ m.17 5~.72 i! 26~.'" ~ 815.85 .!9 '*HI t6~.3~
!
%
I).DD qIS! D1 t."6HE+05 0.3397£-02 257.5~ 35l!.23 292.10 8;9.75 1291.J.18 350.52 1 DD-i. 4H.37 3%.19 C.Ot 0.7587£-02 O.46I-f£+D5 0.3391£-02 '157.91 C.BE!T 12IJCJE 22-4.57 216.72 <17&.1.3 408 3?'J.(j7 257.54 131.90 391.5ij 326.75 1121. 1307.28 11 .,886.81 of "! 316.57 145.63 651.35 0.00 85D.DF !la{,.c? 7ae. (-$.'31 G.GO 18~ 3~.05 98C.46 G.O& G.OO C.I!~ C.De ~D.f17 qz; t.ert:. !L6~ Ii.DC C.O~ 3C.Or. HiE 3~19.D-1 3<J2.11 :K .52 ~+02D.1288£+00 r, .24 .28 .21 72 tCJ 0.7851£-02 0.3891 C.7850£-02 0.3515£-02 STY H?\ H27.2-1 -t55.% £+05 .66 CJt{J 1010.08 2 REI n.De 1 C.7587F-D2 0.1525£+00 329.67 618.10 711.04 1052.25 G CJ<J~.48 •.28 2~.t7 t.OO G.OD e.DO 6.3316£+02 -IOC."1 25~ -1!:"44 -11i'.61 358.23 0/,'; 1377.61 131~.86 t"'12.42 if! .• 8604.52 172.61 576.15 859.75 9D!i.21 915.62 5~I.t2 e.oo t.on 6H~.fC 860.11 Ii.UD C.3~ LCD 96C.~ 515."5 400.09 302.75 ~.Or. t.or 67.n.oe 538. 8Dc;.B6 603.07 CJ.37 791.22 30i3.l'5 5211.13 26\).01 651.35 t.9D c.or. 863.77 O.D~ 301.51 OD2.26 i-f9.% 5.1 .51) !!t! .81 ••.,., 4'i' n.1526£+00 1 CC1.13 129C.-I6 1371.61 1271.77 1142.44 12411.81 195b.+i Hr.2 188-i.34 WJ2.70 1187.61 ItDE.'? 151)-\.39 671i.+I 956.51 ~1.36 671.66 lege.52 ~1.D7 822.22 1977.27 39~.61 480.-i<3 -415.43 CJ99.36 4<J7.2-4 421i.6a 1612..55 581.13 •t 8 2162.20 2e3t.15 2~i6.17 1232.46 176fj, 13(,3.1 1956 H25.61 121C.81 1'957. 1211.7: 1186.'7 1971.21 1136.76 t613~5!i 1C9!.03 10:•.•.• en .51 2862.20 2031.15 2~6.17 1168. 1157.% Ilel.~~ «rl 1038.09 C.2905E+33 295.25 m.31 tli7 1005.9D -175.16 4Q7.43 7.8.5~ 1.•.• 911.93 591.12 /J.OP 0.01.\ 872.12 156.68 2.42 D 0.3265£-02 C.1293E-02 C.3511£-02 0.3893£+05 0.1219£+63 .5625£+05 1163.59 B6U.ao 127j.77 0.1287£+OC 1170.20 C.1861E+OO U12.51 ,~ry.6b ~9C.23 35~.3e 1816.33 1627.97 1008.10 291;.33 1115.79 1-195.52 1271.77 1070.20 187£.33 1et3.D~ 1627.97 2617.33 1089.&6 1&30.29 1037.26 1031.54 111~.52 C.06 42fJ.88 1232 1768.(11 .48 1766.97 ID3~.2a 1001.27 1'B.73 10~1.52 Hj%.19.IC%.6-! Hlt.57 1003.23 ~n02.~1 l~e1.97 101~.7i 1915.26 1&49.86 GEDEiATIOI, IBYD/DF. iT) 5.26
~ D.2888£+a3 &.2808£+03 0.1121t+&3 0.11211+03 ~059.40 TS.~-1095.1! SURFACE tEftPERATur.r OPJ BEG. 2 3 11.26'9£+83 &.29051+03 0.1698£+03 D.1219E+G3 PAS
AM SESUDAE DIMODIFI~~SI
fAR
"m.43 I
1385.32 1469.85 1383.76 1415.61 1440.'1 1~.37 1376.44 1351.10 1~5.16 1381.63 13e9.61 1372.0B 1310.75 80£ 800.80 .00 899.21 798 820.52 892.66 921.54 894 4)311.3-1 873.59 811.19 823.31) 11 2%8.-41 825.10 8t2.H 913.84 822.51 852.52 931.20 aoc.oo 88&.55 flll..43 869.24 796.36 874.31 80C.DD 880.00 5962.78 851 81D.2! 9%.51 9':16.H trJ5.t8 IBO.'H 820.73 B7i 80<).22 C22.23 31108.74 .36 .13 .88 .HI O.8376IiD I3tC.3~ 1328.2; 981.53 940.59 '53,12 m.6S 91~.@7 9D8.48 88<1.36 919.'" 88&.»8 974.30 9n.10 854 85-1.67 M9 836.46 865.45 .18 1281.98 1280.86 1129.57 1115.40 142~.9V 1111.81 1403.29 1369.20 14~I.OB 1334.88 1361.58 1383.45 1333.48 13D4.85 1383.47 13el.63 1300.72 1233.C7 0.20 916.10· 961.21 84£.39 812.21 M4.12 939.74 CJ.4)3 917.19 m.69 11~1.73 1117.12 111&.29 1139.15 1139.19 1118.15 1118.33 981.10 4)52.19 946 985.82 .92t128.52 1196.84 1210.84 1113.40 '5.1.33 1115.86 113~.IH 866.40 e~.76 636.52 636.75 tt17.H 1058.ao 188£.82 9.18 0.45 0.05 0.3D 0 •.•• fAR1388.55 22~.~7 t231.8~ 0.55 1139.95 1112.03 0.25 0.20 '.90 1126.81 0.85 1115.31 TlilCOO1.AITTmEIlATUlE DIG. 0.70 1440.76 1468.6e 1446.30 144:.91 1117.89 1113.~ 139u.6D 1355.0D 1389.13 0.55 1386.30 1426.47 1410.81 1407.61 1386.11 86.17 02.55 1824.98 1107.36 0.58 0.35 C.ge IlL 1106.66 1811.72 1283.23 1020.29 0.65 1.0e 0.80 0.60 0.15 AlUlun C.76 1214.13 1228.32 IHC).t7 130i.73 1028.43 1116.30 1256.15 1203.59 1268.41 BUT 1222.17 OrrLI! LOSS 1051.115 1235.81 103&.36 1112.4~ 1160.12 1180.77 1131.23 TmYIutm:11~.2850&DEC 1121.71 1238.16 1132.36 1176.32 1191.62 1155.46 1039.21 1079.22 1097.18 1150.&1 1&59.19 nR 1041.86 10~t.87 1099.16 1151.86 1062.46 0.15 30.72 1053.19 0.95 1136.15 1.8D 0.75 1137.69 0.88 1143.39 1D85.~ 1181.66 1141.75 1163.59 1887.97 132.68 1406.78 1112.74 1373.26 1438.36 1371~71 1451.76 1424.42 1419.83 1393.83 0.650••1 1429.88 1460.87 1140.67 1436.69 1112.79 1410.13 1390.26 1388.75 364.89 268.01 366.63 1363.48 1267.73 1355.70 1353.26 1354.76 1269.15 0.30 1335.56 1353.50 0.18 1268.13 1287.76 1203.72 1331.58 124£.28 1330.17 1329.31 1242.03 1321.81 1240.67 1319.81 t2~.I£ 122~.18 1233.61 325.24 1319.51 0.98 1125.96 1134.06 1405.60 13419.92 1374.11 0.25 1156.S6 1280.03 1197.77 11~6.61 1187.38 1187.53 1356.46 1191.66 1191.31 ~. 1008.50 0.75 tlIU 1446.41
rmn
106&.33
1070.33
6
1RESSUIE DROP
0.0169
151 IOftBEr. OF conLAn TErlPIUT1RI TB! ESTIIIUn
eoouIT
? lOBES EftrLOYED
TEnf •ElIROI: IS TBf ISTltlATID tltCE 'VIlBE~ IS 0.09
£ .2t::J[G
50 ;.1.1:
529
0.40 . 0.45 0.50 0.55 0.60 0.65
0.70 0.75 0.80 0.85 0.90 0.95 1.00
1378.43 1451.94 1470.32 1451.94 1378.43 1286.53 111\.66 1010.85 827.06 643.27 459.48 183.79 0.00
1553.13 1635.96 1656.67 1635.96 1553.13 1449.59 1320.16 1138.96 931.88 124.19 511.1\ 207.08 0.00
TSJ---
IlL 0.00
0.05 0.10 0.15 0.20 0.25 0.30
0.35 0.48 0.45
0.50 0.55
0.60 0.65 0.70 0.75 0.80 0.85 0.90 0.95
1.00 IlL 0.00 0.05 O. fG
0.15 0.20 0.25 0.30 0.35
C.40 0.45 C.50 0.55 0.60 0.65 0.70 0.75
0.80 0.85
1
40.00 42.69 47.18 50.98 55.39 60.D2 64.85 69.28 13.65 17.97 81.66 84.84 87.18 89.07 90.40 90.80 90.48 89.71 88.52 85.58 83.27 9
4C.80 42.57 +6.63 49.72 53.09 56.55 59.95 62.83 65.52 68.05 6<1.90 11.25 11.80 12.02 71.81 7C.84 6CJ.37 67.70
SUlFiCI
2
40.00 42.83 41.38 50.96 54.93 59.06 63.19 66.79 70.22 73.51 76.08 78.18 79.24 79.98 80.22 79.58 78.32 76.78 74.95 71.+3 6B.85
TEftPERifDRE
3
40.00 42.92 47.57 51.17 55.13 59.22 63.28 66.78 70.B8 73.22 15.59 77.40 71;.31 18.81 78.83 77.95 76.47 74.12 72.12 69.00 66.32
10
40.00 42.61 46.73 49.86 53.27 56.79 60.23 63.'" 65.86 6B.42 70.28 71.64 72.19 72.40 72.18 71.19 6CJ.69 67.98
5
4
40.00 42.89 47.49 51.02 54.89 58.89 62.83 66.21 69.38 72.39 74.63 76.31 77.11 77.50 71.42 76.46 74.91 73.12 71.08 67.37 64.70
BEC.
FAR
5
40.80 42.90 47.52 51.04 54.90 58.88 62.80 66.1+ 69.27 72.23 14.42 76.05 n.78 77.12 76.98 75.96 74.31 72.53 70.46 66.71 64.D2
6
7 40.BO
8 40.80
12.69
42.52 46.51
51~05
46.4)5 50.20
54.90
53.74
52.85
58.88
57.34)
56.26
60.91
59.61
64 .02 66.86 69.55 71.51
62.15
40.00 42.91 47.53
62.78 66.09 6CJ.20 72.13 .74.2CJ 75.88 76.58 76.88
72.4)6
73.57
49.54
65. to 67.59 69.42
70.15 71.31
73.83
71.54
71.35 70.41
66.29
73.64 72.65 71.13 69.40 61.45 63 •CJ6
63.58
61.-46
59.88
76.70 75.65 71.03
72.16 70.07
68.97 67.34 65.58 62.22
5i.;~ I
12.30 "'3.16 42.31 42.ecJ H.S3 52.50 56.84 II210.45 51.87 55.59 59.13 57.Db 11.05 63.90 59.18 61.42 40.DD 52.89 50.01 50.80 58.1M 58.41 60.2-4 41.53 61.81 56.05 63.46 63.41 40.5" 54.16 60.61 50.08 57.36 49." 83.27 71 40.1)<J 5<)
[email protected] 60.8<3 40.'" H.'O
mum iJ
&.99
in
530
531
DAFTAR PUSTAKA
1.N. Kattehee
at all, HECTIC
PROGRAM
FOR HEAT TRANSFER
REACTOR
PASSAGE,
2.W.C. Reynolds, TRANSFER
4.E.R.G Graw
ANALYSIS
IDO 28595,
1962. PROGRAM FOR HEAT
OF GAS COOLED REACTORS,
NUCLEAR HEAT TRANSPORT,
Eekert,ANALYSIS
COMPUTER
OF GAS OR LIQUID COOLED
HECTIC AN IBM 704 COMPUTER
ANALYSIS
3.M.M. Elwakil, 1971 .
II AN IBM 7090 FORTRAN
AGN-TM 381,1961. Int. Textbook
Co,
OF HEAT AND MASS MASS TRANSFER,
Me
Hill Book Co, 1972.
TANYA JAWAB
1. Kasijanto a. Bagaimana
eara
reaktor memakai b. Apa perbedaan
menentukan
bahan
gas atau eairan
parameter
pendingin
?
fisis untuk pendingin
eairan pada program komputernya
suatu gas atau
?
Jawaban
a. Pendingin eara
reaktor ditentukan
menentukannya
oleh jenis reaktor,
tergantung
pada jenis reaktor
jadi apa
yang akan dibuat. b. Program menulis kita liquid
ini dilengkapi memilih
program
harus memasukkan
daya parameter
kontrol
yang akan dihitung, data BM Gas,
kita harus memasukkan
untuk gas
sedangkan
data rapat masa
untuk untuk liquid,
532
parameter panas,
lain
berbeda
yang
viskositas,
serta
antara
lain
kapasitas
antara
density
gas
dan
liquid. 2. A. Syaukati
a. Bagaimana dialiri? b. Apakah
aliran
tidak
transversal
perlu
adanya
energi untuk menghitung
kanal
antar
persamaan
suhu pendingin
dapat
kesetimbangan ?
c. Bagaiman"a dapat diramal adanya burn out ? Jawaban
a. Dalam
hal ini di eliminir
merupakan
besaran
passage·· sehingga aliran transversal b. Suhu pendingin persamaan permukaan c. Samp~
yang
data harga C~k yang
menyatakan
dengan memasukkan
diameter
"inter
harga C~k = 0 maka
antar kanal dapat tereliminir.
sebetulnya
kesetimbangan elemen bakar.
saat
melalui
sudah ditentukan perpindahan
panas
berdasarkan sepanjang
ini saya belum dapat menganalisis
lebih
lanjut dari hasil perhitungan.
Ke Daftar Isi
