## DESAIN ALAT UJI NOSEL DENGAN MENGGUNAKAN PRINSIP TEROWONGAN ANGIN SUPERSONIK
Bagus H. Jihad *) , Dedi Priadi **) *) Peneliti Bidang Propulsi, Pustekwagan, LAPAN **) Dept. Material dan Metalurgi, FT-Universitas Indonesia e-mail: bagusjihad@y7mail.com

## ABSTRACT
The  accuracy  of  nozzle  design  can  be  improved  by  validation.  Two  methods  of validation is software or hardware application. The software can be either CFD software or  special  software  for  the  nozzle  design,  such  as  NOZZLE  ®  or  Aerospike  ®.  The limitations  of  CFD  software  is  on  unsteady  flow  conditions.  while  the  limitations  of special  software  is  the  obscurity  of  assumptions  used.  Therefore,  the  flow  validation experimental  is important to be implemented.
Lapan  has  three  units  of  wind  tunnels,  namely  subsonic,  transonic  and  supersonic. They  have  the  capability  to  measure  the  external  flow  of  the  body.  While  the  nozzle design  purpose  is  mostly  to  investigate  the  internal  flow  of  the  nozzle.  Therefore,  the nozzle flow testing system which is called a mini wind tunnel has been designed.
This  wind  tunnel  is  designed  to  have  speed  of  Mach  2.4,  but  the  effective  speed  of  2 Mach. By using a calculation we obtain the run time of 46 seconds.  If we assume the pressure drop is 40%, then actually the run time is 28 seconds.
Keyword: Windtunnel, Nozzle, Supersonic, Run-time, Schlieren

## ABSTRAK
Ketelitian  pada  Perancangan  Nosel  dapat  ditingkatkan  dengan  melakukan validasi  terhadap  rancangan  tersebut.  Dua  metode  validasi  yang  dapat  digunakan adalah  menggunakan software atau  menggunakan hardware . Software dapat  berupa software CFD  atau software khusus  untuk  rancangan  nosel  seperti  NOZZLE®  atau Aerospike ®.  Keterbatasan software CFD  dan software khusus  adalah  pada  ketidak jelasan  asumsi  yang  digunakan.  Oleh  karena  itu  validasi  aliran  dengan  eksperimen ( hardware )  penting  untuk  dilaksanakan.  Lapan  mempunyai  tiga  unit  terowongan angin,  subsonik,  transonik,  dan  supersonik.  Ketiga  terowongan  angin  tersebut mempunyai  kemampuan  untuk  pengujian  aliran  disekitar  benda  ( external  flow ). Sedangkan  keperluan  desain  nosel  lebih  kepada  penelitian  aliran  pada  bagian  dalam ( internal  flow ).  Oleh  karena  itu,  telah  dirancang  satu  sistem  pengujian  aliran  dalam nosel yang disebut dengan terowongan angin mini karena bentuknya.
Terowongan  angin  ini  dirancang  hingga  kecepatan  2,4  Mach,  tetapi  efektif  pada kecepatan  2  Mach.  Melalui  perhitungan  kami  memperoleh run  time sebesar  46  detik. Dengan asumsi terjadi drop tekanan sebesar 40%, maka run time sesunguhnya adalah 28 detik.
Kata kunci: Terowongan angin, Nosel, Supersonik, Run-time, Schlieren

## 1  PENDAHULUAN
Terowongan  angin  supersonik telah digunakan untuk kegiatan Research  and  Development  (R&D)  lebih dari  5  dekade.  Aerodinamika,  untuk permulaan propulsi dan akustik merupakan  bidang  ilmu  yang  banyak menggunakannya.  Walaupun  aplikasi CFD  pada  rancangan  wahana  antariksa modern meningkat pesat, tetapi penggunaan  terowongan  angin  juga terus  berkembang  terutama  untuk memvalidasi  hasil  perhitungan  meng- gunakan CFD. Pengujian pada terowongan angin  dapat  memberikan  data  pada medan aliran yang kompleks dan sulit.
Terowongan  angin  supersonik, secara  umum  diklasifikasikan  menjadi empat  [Matsumoto,  2000;  Hodge  dan Koenig,  1995;  Zucrow  dan  Hoffman, 1976], yaitu
 Blowdown  Wind  Tunnel (BWT)  dengan media  penyimpanan  gas  bertekanan tinggi.
 High-enthalpy  Wind  Tunnel dengan sebuah shock tube .
 Atmospheric-entry  Wind  Tunnel dengan sebuah pompa vakum.
 Continuous-flow  Wind  Tunnel dengan sebuah kompresor.
Tipe BWT dipilih untuk rancangan ini. Akan tetapi model ini membutuhkan biaya  konstruksi  yang  lebih  besar daripada  sebuah shock  tunnel untuk mengontrol tekanan pada settling chamber .
Ide  untuk  menghasilkan  semburan supersonik  menggunakan  gas  tekanan tinggi  (atau  tekanan  rendah)  dengan sebuah  nosel  supersonik  dipelopori  oleh Carl G. P. de Laval, pada awal abad XIX. Penemuannya menggunakan turbin satu tingkat  yang  digerakan  oleh  uap  panas melalui  sebuah  nosel  konvergen-divergen [Anderson, 2001].
Tipe  BWT  terdiri  atas  tabung penyimpan  yang  diisi  dengan  udara bertekanan dan sebuah nosel konvergen-divergen. BWT membutuhkan tekanan  yang  sangat  tinggi  untuk bilangan  Mach  yang  lebih  besar. Pengujian dimulai ketika gate valve pada ujung  tangki  mulai  dibuka,  udara bertekanan  dari  tangki  akan  mengalir dan  diekspansikan  ke  dalam plenum chamber .  Di plenum  chamber ,  udara bertekanan  akan  mengalami  penenangan ( settle  down ),  oleh  karenanya  bagian  ini juga  sering  dinyatakan  sebagai settling chamber ,  dimana  udara  mengalami penurunan kecepatan dan dipertahankan pada  tekanan  yang  konstan.  Setelah plenum  chamber adalah  nosel  konvergen- divergen  yang  digunakan  untuk  meng- akselerasikan  aliran  jet  supersonik.  Jet ini  kemudian  diturunkan  kecepatannya melalui  sebuah  diffuser  dan  dibuang  ke atmosfir.
Untuk memvalidasi hasil rancangan nosel, maka Bidang Propulsi Lapan telah melakukan rancangan terowongan angin mini.  BWT  yang  akan  digunakan  untuk mendukung  penelitian  dan  pengem- bangan aliran inlet dan nosel supersonik, serta aerodinamika permukaan. Semburan keluar  dari  nosel  dirancang  pada kecepatan  Mach  1,5  –  2,2,  rancangan dioptimalkan pada kecepatan Mach 2.

## 2 DESAIN DAN KONSTRUKSI
2.1  Sistem Suplai Udara
2.1.1 Kompresor
Suplai  udara  bertekanan  tinggi diperoleh  dari  kompresor  motor  listrik 3-phase  berdaya  3HP  dengan  tekanan maksimal 18 Kg/cm 2 dan tekanan kerja sebesar  16  Kg/cm 2 (Gambar  2-1). Kompresor  yang  digunakan  merk  PUMA model  TK30-160  dengan  dua  silinder. Udara  dipompakan  hingga  10.5  atm,
yang merupakan batas maksimal kekuatan tangki penyimpan.
2.1.2 Filter dan air-dryer
Sebuah air-dryer digunakan untuk menghilangkan  uap  air  dari  udara  yang disuplai.  Sebelum  masuk  ke  dalam air- dryer ,  udara  terlebih  dahulu  disaring (dalam air-filter )  untuk  menghambat partikel  dan  uap  air  yang  mungkin terbawa  oleh  proses  kompresi  dalam kompresor.  Setelah  melewati air-dryer , udara  kembali  disaring  oleh  sebuah air- filter sebelum  masuk  ke  dalam  tangki penyimpan.

## 2.1.3 Tangki penyimpanan
Udara  yang  keluar  dari air  filter kedua,  selanjutnya  masuk  ke  adalam tangki  penyimpanan.  Tangki  mampu menyimpan  udara  kering  sebesar  2  m 3 , pada suhu hingga 150 ° C. Pada tekanan 10.5  Kg/cm 2 atau  sekitar  10,5  atm. Tangki ini dilengkapi dengan water-drain di  bagian  bawah,  dan vent-pressure di bagian  atas  untuk  mencegah  terjadinya over-pressure . Di bagian samping dilengkapi dengan sebuah pressure-gage dan sebuah valve .

## 2.2 Komponen
Komponen-komponen yang diguna- kan  dalam  rancangan  BWT  ini  dapat dilihat  pada  Gambar  2-1.  Pada  gambar [a]  kompresor  2-silinder  digerakan  oleh motor  listrik,  serta  dilengkapi  tabung penyimpan.  Pada  bagian  ini  dipasang- kan  regulator  tekanan  untuk  menjaga agar kompresor dapat berhenti beroperasi pada  tekanan  maksimal  10.5  Kg/cm 2
(terbaca  pada  indikator).  Udara  yang keluar kompresor difilter yang dilengkapi
dengan  regulator  tekanan.  Kapasitas maksimal  tekanan  filter  adalah  16 Kg/cm 2 ,  sedangkan  pada  BWT  ini  diset pada tekanan 10,5 Kg/cm 2 . Bila tekanan berlebih  akan  terbuang  di  bagian  atas tabung,  sehingga  keamanan  dapat terjaga.
Tangki  penyimpan,  gambar  [d] berbentuk dome terdiri  atas;  dua  buah lubang inlet dan outlet (3”);  lubang manhole di bagian bawah depan; indikator tekanan diset pada tekanan 10,5 Kg/cm 2 . Udara  yang  masuk  ke  dalam  tangki penyimpan,  sebelumnya  dikeringkan oleh air-dryer (gambar  [b]).  Udara  yang keluar  dari  tangki  dapat  diatur  dengan gate valve berupa butterfly valve [c] yang dapat diatur bukaanya untuk masuk ke dalam storage  tube .  Di  ujung  dari storage  tube dipasangkan screw  valve tipe butterfly [e],  dimana  setelahnya dipasangkan  sebuah  diffuser  yang berguna  untuk  mengerem  kecepatan udara  yang  masuk  settling  chamber. Udara  yang  masuk  diatur  sedemikian rupa  pada  tekanan settling  chamber [f] yang  diharapkan,  dengan  menggunakan screw  valve .  Kondisi  ini  dapat  dibaca pada pressure transducer yang dipasang- kan  pada settling  chamber ( plenum chamber ) .
Bagian setting atau plenum chamber digunakan  untuk  menenangkan  aliran sebelum  masuk  ke  dalam  model  nosel yang  diuji.  Panjang  bagian  ini  harus sedemikian rupa, sehingga pressure loss yang  terjadi  tidak  berlebihan.  Dengan menggunakan persamaan mekanika fluida  maka  panjang  bagian  ini  dapat diperoleh,  dan  selanjutnya  diverifikasi menggunakan software CFD.
Indikator tekanan dan Kontrol Otomatis (diset pada 10 bar) [a] Air Dryer (untuk mengeringkan udara yang akan menuju tabung) [b] Katup kupu-kupu yang dipasang tepat setelah tabung (  3”) [c] Tabung udara (2 m 3 ) [d] Katup kupu-kupu yang dipasang sebelum memasuki difuser(  3”). Diameter difuser (  5”) [e] Katup , difuser dan settling chamber (tekanan diatur sesuai dengan Number of
Pressure Ratio/NPR) [f]
Gambar 2-1: Komponen pada Terowongan angin mini
Gambar 2-2: Rancangan  Nosel  planar  2-Dimensi  (kiri),  setelah  dibuat  dan  terpasang (kanan)

## 2.3 Sistem Imaging
Telah dilakukan pemilihan sistem imaging yang  akan  digunakan  dalam pengujian.  Dari  4  macam  metode  yang ditelaah,  yaitu interferometer , shadow- graph ,  cermin schlieren dan  lensa schlieren menggunakan  analisis  Tekno- Ekonomi,  maka  bobot  terbesar  terletak pada sistem cermin schlieren . Pembobotan diberikan  pada  Tabel  2-1.  Dengan demikian  sistem  cermin schlieren yang
akan  dipergunakan  dalam  penelitian. Telah disiapkan cermin schlieren dengan diameter 6 inci, dengan area yang dapat ditangkap  ±  10  cm 2 .  dengan  demikian maka  model  nosel  yang  dipabrikasi harus  mempunyai  area  lebih  kecil  dari nilai itu.
Guna  melengkapi  sistem  imaging tersebut,  juga  disiapkan high  speed camera yang  dapat  menangkap  obyek hingga 1200 fps.

## Gambar 2-3: Skema instalasi cermin schlieren

## Tabel 2-1: SISTEM PENILAIAN BEDASARKAN BOBOT KEPENTINGAN
No.  Karakteristik Nilai ideal Bobot Perbandingan nilai Shadowgraph  Interferometer Schlieren Lensa  Cermin 1. 2. 3. 4. 5. 6. Man.ability (a) Man.ability (b) Dimension Adaptability Image quality Environment 4 4 4 4 4 4 5 5 8 9 10 7 4 4 4 4 1 4 2 2 2 2 4 2 4 4 4 4 3 4 3 4 3 4 3 4 Total 24 44 21 14 23 21 Bobot x Nilai 176 146 108 166 153 Nilai Teknis 1 0.830 0.614 0.943 0.869 Nilai Ekonomis 1 1 10 1.2 1 Tekno-Ekonomis 1 0.830 0.061 0.786 0.869 [Jihad, Samosir, Soemardi, 2006]

## 3  PREDIKSI  GAYA  DORONG  DAN ESTIMASI RUN TIME

## 3.1 Prediksi Gaya Dorong
Dengan  menggunakan  program Matcad,  diperoleh  bentuk  desain  nosel planar 2-dimensi. Nosel tersebut dirancang untuk  mendapatkan  kecepatan  keluar sebesar  2  Mach  dengan  media  udara. Pada  kecepatan  itu,  rasio  ekspansi diperoleh  sebasar  ε  =  1.6875  [Jihad,  et al., Implikasi ,  2008  dan  Jihad  et  al., Aliran , 2008].
Fluida  yang  akan  digunakan adalah  udara  dengan  nisbah  panas jenis,  γ  =  1,4  dan  tetapan  gas  universal R  =  287  J/(Kg.K).  Untuk  kecepatan desain M d = 2 dengan A e /A t = 1.687, dari tabel  isentropik  untuk  P  =  7.825  atm [Anderson,  2001],  diperoleh:  P e /P 0 = 0.1278; dan T e /T 0 = 0,5556
Untuk kondisi optimum, dimana P 0 = 1 atm, maka diperoleh:
- P e = 0.1278 x 7.825 = 1.000035 atm
- T e = 0.5556 x 300 K = 166.68 K Kecepatan keluar yang dihasilkan nosel:
- T R M V e . .   diperoleh  harga  Ve  =
517.405 m/s
- T R P e e .   diperoleh harga 2.0905 kg/m 3
Dengan  demikian,  maka mass  flow  rate saat keluar nosel:
e e A Ve m . .    =  2.0905  x  517.34  x  3.848
10 -3  4.1616 Kg/s
Dengan  demikian  gaya  dorong yang  dapat  dihasilkan  oleh  nosel  yang telah  dirancang  adalah  (pada  tekanan 7.825 atm):
Ve m F .   =  4.1616  x  517.405  2153.2327 N (215.3233 Kg).
M3 viewing plane Test Section Knife edge M1 M2 F2 F1 light source slit
3.2 Estimasi Run Time
Run  time terowongan  angin dihitung menggunakan persamaan (3-1), Pada  persamaan  ini ,  run  time hanya merupakan  fungsi  dari  kondisi storage tank dan penampang throat nosel .
                n i S f S t i S t S run P P P T A P T V t 1 2 2 * 2 1 . . . . . 03526 . 0 (3-1)
Persamaan
ini memberikan estimasi  yang  cukup  akurat  untuk terowongan dengan run-time yang cukup lama.  Untuk  terowongan  angin  dengan run  time yang  lebih  pendek,  beberapa faktor  harus  ditambahkan.  Proses  awal dimulai  ketika  udara  bertekanan  me- masuki plenum  chamber sampai  posisi tekanan yang diinginkan. Jumlah udara yang  dapat  dihembuskan  pada  proses ini  diasumsikan  sekitar  2  kali  volume dari plenum chamber [Matsumoto, 2000]. Masa  udara  efektif  dapat  dihitung dengan persamaan,
1 1 . 2 . . . . . V V V t A U t S f S S i S run        (3-2)
Drop tekanan melalui katup dan difuser dinyatakan  oleh  (1  –  χ),  diasumsikan sekitar  40%,  sehingga  persamaan  (2) menjadi:
         
                                         1 S i S f S S 2 t 2 i S i 2 t run V 2 V P P V P T . * A P . T 03526 . 0 t (3-3)
Hodge  dan  Koenig,  1995,  memberikan rumusan  dengan  memasukan  variabel fluida yang digunakan,
 
      1 / 1 2 1 1 2 . * 1 1 / . 2                           oi oi of s RT A P P V t (3-4)
Dengan  menggunakan  kondisi desain, dimana tekanan atmosfir 1 atm. Dengan  kecepatan  keluar  nosel  M e =  2, tekanan  tangki  penyimpan  10,5  atm (154 psi)  pada  suhu 30°C  (540  R).  Pada rasio  ekspansi,  ε=A e /A t ,  1,687  diperoleh rasio tekanan P of /P o = 0.12780.
Suplai  tekanan  untuk  mendapat- kan  kecepatan  M  =  2  akan  berlaku hingga  tekanan  tangki  drop  ke  tekanan dimana timbul shock wave normal pada bagian exit nosel. Reduksi tekanan lebih lanjut  di  dalam  tangki  penyimpan  akan menghasilkan shock normal yang bergerak  ke  belakang  nosel.  Pengujian harus  dihentikan  pada  kondisi  ini karena  aliran  tidak  lagi  supersonik. Untuk  M e =  2,  bilangan  Mach  setelah shock  wave normal  adalah  M 2 =0,5774 dengan  rasio  tekanan  pada shock  wave P 2 /P 1 =4.5. Dengan asumsi bahwa proses berjalan  dalam  kondisi quasi-steady , tekanan tangki pada kondisi ini adalah:
(atm) 739 . 1 7 , 14 5 , 4 1 1278 , 0 1 P P P P P P P P b b 2 2 1 1 o f S   
Menggunakan persamaan (3-4) diperoleh run  time sebesar  46,82  detik.  Dengan mengambil  drop  tekanan  sebesar  40%, maka run  time sesunguhnya  sebesar 28.092 detik.
Gambar 3-1: Instalasi terowongan angin
Keterangan:
1. Kompresor 3-phase, 3HP.
2. Valve.
3. air-filter.
4. Regulator.
5. Oil-filter.
6. Diffuser adaptor ½” - ¾”.
7. Air-dryer.
8. Diffuser-Adaptor ¾” - 1”.
9.   Vent-valve.
10. Tangki penyimpan.
11. Pressure indikator tangki.
12. Water-oil separator.
13. Butterfly valve.
14. Water drain.
15. Butterfly valve (screw type).
16. Diffuser 3”-5” (Plenum chamber).
17. Nosel.

## 4  HASIL DAN PEMBAHASAN
Setelah  tangki  diisi  penuh,  maka peralatan  pendukung  disiapkan,  antara lain:  sistem  schlieren, high  speed camera , pressure transducer beserta data akuisisinya.
Gambar 4-1: Peralatan  disiapkan  untuk pengujian
Pada percobaan pertama peralatan belum dipasang, hal ini untuk mencegah terjadinya hal-hal yang  tidak diharapkan. Semburan  nosel  dibiarkan  ke  udara terbuka  ( free  jet ).  Percobaan  diulangi
hingga  beberapa  kali  agar  komponen- komponen  yang  tersisa  akibat  proses pabrikasi  telah  benar-benar  keluar. Setelah dirasakan cukup, maka terowongan dipasang untuk mengarahkan semburan  keluar  ruangan  uji.  Pada pemasangan  pertama,  terowongan  jebol karena  terhisap  oleh  semburan.  Setelah dilakukan perbaikan menggunakan multipleks, dilakukan pengujian berikut- nya,  kali  ini  multipleks  juga  jebol  dan menyumbat aliran (Gambar 4-2). Besarnya gaya  sedot  dari  semburan  yang  keluar nosel, menyebabkan dinding terowongan harus dimodifikasi untuk ketiga kalinya.
Setelah  semua  hal  dipersiapkan, maka  dilakukan  percobaan  dengan bukaan  katup  setengah  dan full . Indikator  tekanan  tangki  menunjukkan 9.5  bar  pada  saat  udara  dihembuskan. Gambar  4-3  dan  4-4,  memberikan  hasil pengujian  yang  ditangkap  oleh high speed  camera untuk  dua  buah  nosel dengan  konfigurasi  jejari throat yang berbeda.
Hasil yang diperoleh untuk kedua model  tersebut  berbeda,  hal  ini  sesuai dengan yang telah disimulasikan menggunakan CFD [ Jihad et al., 2008].
Pada  pengujian,  tekanan  yang terukur  hanya  pada  bagian settling chamber ,  sehingga  perlu  diupayakan agar tekanan sepanjang dinding divergen nosel dapat terukur. Dengan terukurnya tekanan sepanjang dinding/sumbu nosel,  maka  pola  aliran  dalam  nosel dapat  digambarkan,  hal  ini  penting
untuk  mengetahui  apakah  desain  nosel yang telah dilakukan optimal atau tidak.
Selama  proses  pengujian,  timbul uap  air/embun  pada  dinding  acrilyc nosel  yang  mengganggu  proses  peng- ambilan  gambar.  Hal  ini  diakibatkan oleh penurunan temperatur pada bagian nosel.  Untuk  itu  perlu  diupayakan pemasangan heater ,  untuk  menaikkan suhu udara pada tangki menjadi sekitar 400  K.  Hal  ini  penting  agar  pembacaan tekanan  pada pressure  transducer tidak bias akibat aliran dua-phasa ini.
Gambar 4-2: Terowongan yang rusak akibat sedotan udara yang keluar dari nosel (kiri: seng, kanan: multipleks t = 1 cm).
Gambar 4-3:Hasil imaging pola aliran dengan schlieren pada nosel planar 2D

## Gambar 4-4:Hasil imaging pola aliran dengan schlieren pada nosel N_16
Gambar 4-5: Terjadinya  embun  pada  bagian  acrilyc  nosel  menyebabkan  terjadinya bayangan hitam pada imej yang ditangkap kamera

## 5  KESIMPULAN DAN SARAN
Dari  rancangan  hingga  pengujian yang telah dilakukan dapat disimpulkan beberapa hal:
 Terowongan angin mini yang dirancang untuk mendapatkan pola aliran dalam nosel telah berjalan dengan baik.
 Tekanan  yang  terukur  hanya  pada bagian settling  chamber atau inle t nosel,  sehingga  perimeter  nosel  belum dapat ditentukan dengan baik.
Diperlukan  paling  tidak  tiga  titik tekanan, yaitu pada bagian inlet, throat
dan outlet.
 Terjadinya  embun  pada  acrilyc  nosel menyebabkan  pengambilan  gambar menjadi  kurang  sempurna,  perlu diupayakan  menaikkan  temperatur udara.

## DAFTAR RUJUKAN
Anderson,  John  D.  Jr.,  2003. Modern Compressible Flow:  With Historical Perspective ,  3 rd ed,  McGraw-Hill Companies, Inc., New York, pp.183-186.
Gamble, E., Terrell, P.E., D., DeFrancesco, R.,  2004. Nozzles  Selection  and Design Criteria , AIAA-2004-3923.
Gokhale,  S.  S.,  dan  Suresh,  R.,    1997.
Numerical Computations of Internal Flows  for  Axisymmetric  and  Two- Dimensional  Nozzles ,  Int.    Journal for  Numerical  Methods  in  Fluids, Vol. 25,  hal 599-610.
Jihad  Bagus  H.,  Priadi,  Dedi,  Soemardi, Tresna  P.,  Siradj,  Eddy  S.,  2008. Desain  Divergen  Nosel  Supersonik 2-Dimensi  Menggunakan  MOC , Seminar  Nasional  Mesin  dan Industri  (SNMI4),  Jakarta  28
Agustus 2008.
Jihad  Bagus  H.,  Priadi,  Dedi,  Soemardi, Tresna  P.,  Siradj,  Eddy  S.,  2008.
Aliran Turbulen pada Nosel Konvergen  Divergen  Tipe  Bel , Seminar Nasional Teknologi Simulasi  IV  (TEKNOSIM  2008), UGM 16 Oktober 2008.
Jihad  Bagus  H.,  Priadi,  Dedi,  Soemardi,
Tresna  P.,  Siradj,  Eddy  S.,  2008.
Aliran Inviscid Melalui Nosel Supersonik Konvergen-Divergen , Seminar Nasional IPTEK Dirgantara XII Tahun 2008, Jakarta, 5 November 2008.
Jihad  Bagus  H.,  Priadi,  Dedi,  Soemardi, Tresna  P.,  Siradj,  Eddy  S.,  2008.
Implikasi  Desain  Kurvatur  Throat Terhadap Kinerja Nosel Konvergen-
Divergen , Seminar
Nasional
Tahunan  Teknik  Mesin  (SNTTM- VII), Manado, 4-6 November 2008. Jihad  Bagus  H.,  Priadi,  Dedi,  Soemardi, Tresna  P.,  Siradj,  Eddy  S.,  2008. Pengaruh  Dimensi  Kontur  Throat Nosel Konvergen-Divergen Terhadap  Gaya  Dorong ,  Seminar Nasional  Tahunan  Teknik  Mesin (SNTTM-VII), Manado, 4-6
November 2008.
Jihad, Bagus H., Samosir, G., Soemardi, Tresna  P.,  2006. Studi  Pemilihan Sistem  Visualisasi  Aliran  Untuk Model  Nosel ,  Seminar  Nasional Tahunan Teknik Mesin V, Departemen Teknik Mesin,
Universitas  Indonesia,  Depok,  21- 23 November 2006.
Matsumoto, Joji, 2000. Design and Testing of  a  Subscale  Supersonic  Aero- propulsion  Wind  Tunnel ,  Thesis, Dept.
Aerospace Engineering, Universitas Texas.
Zucrow, Maurice J. and Hoffman, Joe D., 1976. Gas Dynamics : Vol. I dan II,
John Wiley & Sons, New York.

## Daftar Notasi
A Luas penampang M Bilangan Mach n Eksponen polytopic p Tekanan R Bilangan Reynold t Waktu t run Run time Δt Waktu delai T Temperatur U Kecepatan udara V Volume  Rasio panas spesifik χ Faktor drop tekanan  Densitas udara subscript av Nilai waktu rata-rata S Tangki penyimpan t Waktu total 0 Nilai stagnasi 1 Plenum chamber 2 Seksi uji superscript e Akhir jendela uji i Kondisi awal uji j run ke-j f Kondisi final uji * throat nosel