Abstrak —Studi ini bertujuan untuk melihat karakteristik aliran fluida  yang  melalui backward  facing  inclined  step dengan penambahan  paparan  panas  dari high  rise  building pada  sisi upstream .  Peninjauan dilakukan dengan menggunakan perangkat lunak  CFD  ( Computational  Fluid  Dynamic). Pemodelan  ini menggunakan  aliran steady  flow,  incompressible  flow, dan uniformflow pada  sisi  inlet. S imulasi ini  menggunakan model turbulen k-ω SST dengan kecepatan angin pada inlet 3 m/s dan 9 m/s  (Re=2×10⁴  dan  Re=6×10⁴). Meshing menggunakan Hexa- hedral dengan jumlah node 2,5×10⁶. Variasi inclined step 30⁰, 45⁰, dan 90⁰. Kecepatan fluida panas dari gedung yaitu 1,5 m/s dengan temperatur 64⁰ C. Angin lingkungan memiliki temperatur 27⁰ C. Hasil yang  didapatkan  dari  penelitian  ini adalah variasi  sudut inclined  step mempengaruhi  kecepatan fluida,  panjang reattachment, dan  luas  wilayah  yang  mengalami  kenaikan temperatur.  Perbandingan  kontur  kecepatan  memperlihatkan luas area backflow paling sempit berada  pada inclined step 30⁰. Pada pathline tiap  variasi inclined  step diketahui  bahwa reattachment paling  pendek  berada  pada  model  dengan  sudut inclined step  30⁰ . Berdasarkan  kontur pada  posisi y/h=0,5  dan y/h=1  pada  daerah downstream ,  persebaran  temperatur  paling sempit ada pada domain dengan sudut inclined step 30⁰. Selain itu variasi bilangan Reynolds juga berpengaruh terhadap persebaran temperatur.  S emakin  tinggi  nilai  bilangan  Reynolds  maka semakin kecil kenaikan temperatur lingkungan yang terjadi.
Kata  Kunci — high rise  building, inclined  step, paparan  panas , velocity profile, temperature  contour .

## I. PENDAHULUAN
NDONESIA  memiliki  kekayaan budaya dan keindahan ala m yang  menarik  bagi  wisatawan.  Perlu  adanya  dukungan terhadap  daerah  wisata  berupa  akses  ke  tempat  wisata, transportasi, dan  penginapan.  Adanya pembangunan  fasilitas tersebut  mempengaruhi lingkungan  sekitarnya,  terutama pembangunan penginapan. Metode numerik untuk menganalisa hubungan antara pembangunan dan lingkungan merupakan alat untuk engineer , arsitek, perencanaan kota, dan kebijakan  public untuk  menggabungkan  desain perkotaan  alternatif dan  untuk membuat  kebijakan  pedoman. CFD  adalah salah satu alat dari metode  numerik  tersebut, yang  sering  digunakan  pada  iklim perkotaan  dalam  berbagai  skala  [1].  Pada  simulasi  angin lingkungan,  biasanya menggunakan  model  turbulen.  Terdapat lebih dari satu jenis model  turbulen pada CFD.  Dari  11 turbulen model,  kontur  dari  turbulen  model  dengan  grup  yang  sama memiliki  hasil yang mirip  [2].
Variasi  dari  simulasi  dapat berupa bentuk domain . Dima n a pada kenyataanya wilayah  yang memiliki  kondisi  tanah yang rata  hanya sebagian kecil. Domain dapat berupa enlargement yang memiliki back ward facing incined step . Step dari domain pun  memiliki  sudut  yang  bervariasi.  Semakin  kecil  sudut inclined step makan  semakin  pendek  jarak  titik reattachment dengan step [3]. Namun apabila sudut inclined step terlalu kecil makan  panjnag reattachment tidak  dapat  terukur  pada eksperimen  [3].
Hasil  dari  penelitian  aliran  fluida  yang melewati back ward facing  inclined step dengan penambahan  paparan panas  dari gedung pada sisi upstream adalah mendapatkan data pedukung yaitu data kuantitatif  dan data kulitatif.  Data kuantitatif  berupa grafik  distribusi koefisien tekanan (Cp ), sk in friction coefficient (Cf),  dan  jarak reattachment pada  darah downstream .  Data kualitatif  (visualisasi  aliran)  antara  lain  berupa  kontur kecepatan, kontur tekanan, dan pathline dari vektor kecepatan.

## II. PENELITIAN

## A. Urban Microclimate
Urban  Microclimate semakin  populer  terutama  di  negara- negara  maju sebab  urbanisasi  yang  cepat.  Beberapa dokumentasi  studi Urban  Microclimate membahahas  tentang konsumsi energi gedung, kondisi kesehatan dan kualitas hidup manusia,  tingkat harapan hidup, dan temperatur  yang nyaman untuk  hidup.  Sebagian  besar  penelitian urban  microclimate , baik  eksperimen  maupun  simulasi,  dilakukan  pada  bidang datar.  Hal  ini  dikarenakan  fokus  pada  penelitian  adalah perkotaan  dimana  sebagian  besar  wilayahnya  datar.  Namun sekarang  pembangunan  telah  sampai  ke  berbagai  wilayah. Termasuk  pedesaan  maupun  pinggiran  kota.  Pemandangan yang masih alami  menjadikan  bangunan modern memiliki  nila i jual tersendiri. Selain  itu, pembangunan dilakukan  pada daerah yang lebih  tinggi  dari  sekitarnya.  Hal  ini  akan  membuat  para penghuni mendapatkan pemandangan terbaik apabila bangunan berada  lebih tinggi  dari sekitarnya.  Tidak menutup kemungkinan  bahwa daerah di  bawah bangunan tersebut akan semakin  padat.  Maka  itu  diperlukan  penelitian urban microclimate dan CFD hadir  untuk  membantu  memberika n simulasi  pada permasalahan urban microclimate .
B. Geometri Benda Uji
Penelitian secara numerik ini dilakukan dengan menyimulasikan  model  (gambar  1)  dengan back ward facing

## Studi Numerik Karakteristik Aliran Melalui

## Backward Facing Inclined Step dengan Penambahan  Paparan Panas Deri Gedung pada Sisi Upstream
Franciska Enstinita  Puspita dan Wawan Aries Widodo Departemen Teknik Mesin,  Fakultas Teknologi  Industri,  Institut Teknologi  Sepuluh  Nopember e-mail : wawanaries@me.its.ac.id
I
inclined  step dengan  penambahan  gedung  pada  bagian upstream. Referensi  y/h=0  berada pada sisi inlet yang sejajar dengan wall bawah (gambar  3).
Gambar  1.  Permodelan dom ain  backward  facing  inclined  step (a) tampak samping;(b) tampak atas; dan (c) dimensi kondensor (bawah).
Spesifikasi  dari domain adalah sebagai berikut:
 Tinggi Step              =  h
 Panjang sisi upstream =  10h
 Panjang sisi downstream = 34h
 Tinggi inlet =  9h
 Tinggi outlet =  10h
 Lebar geometri  penampang  =  10h

## Dimensi  Gedung dan kondensor
 Tinggi gedung            =  h
 Lebar gedung            =  0.5h
 Panjang gedung          =  0.5h
 Tinggi area kondensor      =  0.2h
 Lebar area kondensor      =  0.3h
 Panjang area kondensor    =  0.3h
Variasi  yang akan dibuat pada model adalah sudut back ward inclined step (α), yaitu: 30⁰,  45⁰,  dan 90⁰.  Kemudian  dilakukan pula variasi kecepatan inlet , yaitu 3 m/s (Re=2×10⁴)  dan 9 m/s (Re=6×10⁴).
C. Metode Numerik
Penelitian  ini  adalah  dengan menyimulasikan  proses  fluida melewati domain sehingga dari  penelitian  ini  dapat diketahui kecepatan  fluida  dan  kenaikan  temperatur  pada  bagian downstream . Penyimulasian  pada penelitian  ini  memanfaatka n komputasi  fluida  (CFD)  untuk  melakukan  pendekatan secara numerik  dan membuat  geometri  model.  Pada  pemodelannya, model  dibuat sedemikian  rupa hingga menjadi  seperti geometri pada gambar  2 serta mendiskritisasikan  model  ( meshing ) yang rapat pada dinding sehingga model dapat diproses dengan CFD. Secara  umum,  urutan  kerja  yang  dilakukan  pada  penelitian numerik  ada  dua  tahapan,  yaitu  tahap pre - processing dan dilanjutkan dengan tahap post - processing .
Gambar 2. T ahap Pre-processing (a) rapat dekat wall, (b) rapat dekat gedung.
Gambar 3. Boundary Condition.

## D.Pengambilan Data
Pengambilan  hasil  setelah  dilakukan  iterasi.  Hasil  tersebut berupa data  kuantitatif  dan kualitatif.  Data  kuantitatif  beru pa distribusi nilai kecepatan, temperatur, Cp dan Cf. Data kualitatif berupa penampilan pathline velocity profile, kontur kecepatan, dan  kontur  temperatur  serta  variasi-variasi  yang  digunakan pada penelitian ini.
Gambar  4.  Posisi  pengambilan  data  (a)  kecepatan  dan  temperatur;  (b) temperatur; (c) Cp, Cf, jarak reattachment, kecepatan, dan temperatur.
Selain  itu, untuk menunjang data kualitatif kontur kecepatan dan  kontur  temperatur  diperlukan  data  kuantitatif.  Data kuantitatif  dipresentasikan  dalam  bentuk  grafik.  Adapun pengambilan  grafik  kecepatan  ditampilkan  pada  Gambar  5c pada posisi x/h=0,75  (melewati  gedung) dan sudut inclined step 30 ⁰  (Re=2×10⁴).  Data  grafik  kecepatan  diambil  dari  posisi y/h=0 sampai y/h=2,2 di titik z/h=-12,
(c) (a) (b) y z x (b)
-15,  -20,  dan  -30.  Pengambilan  data  grafik  distribusi temperatur  ditampilkan  pada  Gambar  4b.  Pada  Cp  dan  Cf diambil  pada wall bawah  sepanjang downstream seperti ditunjukkan  pada  Gambar  4c.  Sedangkan  kontur  kecepatan diambil  pada  posisi  gambar  4a  dan  4c.  Kontur  temperatur diambil  pada semua posisi pada gambar  4. Pathline kecepatan diambil  pada posisi gambar 4c.
E. Grid Independency
Analisa grid  independency menggunakan  model  turbulen dengan menganalisa  kecepatan  rata-rata  ( Vavg) untuk  setiap variasi mesh yang  dibuat  pada  posisi  z/h=-8  dengan  sudut inclined step 90⁰,  sehingga dapat dilihat  pada tabel 1.
T abel 1. Grid independency ditinjau dari kecepatan rata-rata pada cross section z/h=-8. Mesh Jumlah node (juta) Kecepatan rata-rata (m/s) Error (%) Mesh A 1,5 2,9284 - Mesh B 2 2,9381 1,04 Mesh C 2,5 2,9376 0,05 Mesh D 3 2,9370 0,06
Dari  hasil grid independency yang telah dilakukan  diperoleh variasi  kerapatan mesh yang tidak  terdapat perubahan secara signifikan  yaitu  pada  variasi mesh tiga dengan jumlah mesh nodes sebesar 2,5  juta  nodes. Selain  itu,  untuk  memperjela s hasil  dari grid  independency dibuatlah  suatu  grafik  yang menunjukan  hubungan antara jumlah  node dan kecepatan rata- rata.
Gambar 5. Grid independency kecepatan ditinjau pada z/h=-8.

## III. HASIL DAN DISKUSI
A. Variasi Sudut Back ward Facing Inclined Step
Distribusi  kecepatan  aliran  fluida  yang  melewati  gedung dengan variasi sudut inclined step memiliki  grafik  yang sama. Gambar  garfik  6  yang diambil  dari  bebrapa titik  pada daerah downstream memperlihatkan  bentuk grafik  yang sama. Gra fik pada  z/h=-12  terlihat  penurunan  kecepatan  akibat  adanya back flow kemudian  kecepatan  berangsur-angsur  normal. Apabila  data diambil  lebih  jauh  dari  step, pada titik  z/h=-15, maka  akan  terlihat  pengaruh back flow tidak sebesar pada titik z/h=-12.
Gambar  6.  Perbandingan  profil  kecepatan  {u/Uinlet)  dengan  variasi  sudut incined step pada posisi x=0,75 (Re=2×10⁴); (a) z/h=-12; (b) z/h=-15.
Gambar  7  memperlihatkan  perubahan  warna  yang  mewakili perubahan kecepatan. Nilai  kecepatan dari  setiap warna  dapat dilihat pada  skala  warna. Warna  biru  pada  kontur mengindikasikan  adanya back flow .
Gambar 7. Visualisasi kontur kecepatan dengan variasi sudut inclined step pada posisi x=0,75 (Re=2×10⁴).
Perbandingan kecepatan setiap variasi step juga dapat dilihat pada gambar 8 dimana terdapat perbedaan luas area back flow .
Gambar  8. (a)Visualisasi  aliran 3d  kontur velocity magnitude dengan sudut inclined step 30⁰ dan (b) Perbandingan velocity magnitude pada z=-12h dan z=- 15h (Re=2x10⁴).
Pada  gambar  9,  perbandingan kontur  kecepatan  pada  variasi bilagan  Reylonds  memiliki  hasil  yang sama .  Domain dengan step α=30⁰  memiliki  area back flow yang lebih sempit dari pada domain dengan sudut step lainnya.
2.925 2.93 2.935 2.94 1,5 2 2,5 3 Kecepatan rata-rata (m/s) jumlah node (juta) 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 0 0.1 0.2 0.3 0.4 0.5 0.6 y/h u/U inlet 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 0 0.1 0.2 0.3 0.4 0.5 0.6 y/h u/U inlet (a) (b) (a) (b) m/s m/s
Gambar 9.  Perbandingan profil kecepatan dengan variasi   Reynolds number pada posisi x=0,75 (α=30⁰).
Variasi  sudut inclined  step juga  berpengaruh  terhadap jangkauan  kenaikan  temparatur  akibat  paparan  panas  dari gedung.  Jangkauan  kenaikan  temperatur  terjauh  ada  pada domain dengan α= 90⁰. Hal tersebut terlihat pada gambar 10 dan 11.  Gambar  10  dan  11  merupakan  kontur  temperatur  pada daerah downstream .  Gambar  12  adalah  grafik  perbandingan temperatur  lokal  dengan  tempratur  inlet  pada  daerah downstream .
Gambar 10. Kontur temperatur pada posisi y/h=0,5 (Re=2×10⁴).
Gambar 11. Kontur temperatur pada posisi y/h=1 (Re=2×10⁴).
Gambar 12. Grafik distribusi temperatur sepnjang z/h pada x=0,75 (Re=2×10⁴); (a) y/h=0,5;(b) y/h=1.
Gambar  kontur  diambil  pada plane y/h=0,5  dan y/h=1 atau setinggi step. Pada Gambar  10 terlihat jelas adanya peningkatan temparatur di sepanjang sumbu z sampai di outline . Pada α=30⁰ terlihat  bahwa  kenaikan  temperatur  bisa  mencapai  1,18⁰  C. Namun  area pada temperatur  ini  lebih  sempit  dari  pada sudut step yang lain.  Jangkauan kenaikan  temperatur  terjauh  berada pada downstream dengan α=90⁰.  Gambar  11,  jika  ketinggian plane  kemudian  dinaikan  menjadi  sejajar  dengan  step maka akan didapatkan gambaran yang serupa. Dimana  area kenaikan temperatur  tertinggi,  yaitu sekitar  1,64⁰  C,  paling  luas  berada pada domain yang  memiliki  sudut inclined step 90⁰.  Dari gambar  12  menunjukan  perbandingan  kenaikan  temperatur local  terhadap temperatur  inlet mencapai  1,0042  (y/h=0.5)  dan 1,0055(y/h=1).
Area  terjauh  yang mengalami  kenaikan  temperatur  berada pada α=90⁰  sedangkan area tersepit berada pada α=30⁰.  Hal  ini karena  fluida  panas  dari  gedung  terjebak  pada wak e akibat adanya englagement . Pada domain α=90⁰  memiliki wak e yang lebih luas areanya sehingga fluida panas dari gedung lebih lama terjebak  pada  daerah downstream. Sedangkan  pada domain α=30⁰ memiliki wake yang lebih  sempit sehingga fluida  panas yang terjebak juga lebih sedikit. Hal ini membuat area kenaikan temperatur tidak begitu luas dari pada domain lainnya.
Pada  profil  temperatur,  perbedaan  kecepatan  aliran  fluida mempengaruhi  transfer panas ke lingkungan sekitarnya. Secara teori,  semakin  kencang laju  aliran  fluida  maka  semakin  dingin temperatur  lingkungan  sebab tidak  terjadi  penumpukan  panas pada daerah downstream .
(a) (b) m/s (⁰C ) (⁰C )
Gambar  13.  Perbandingan  visualisasi  kontur  temperatur  dengan  variasi Reynolds number pada posisi y/h=1 dan α=30⁰.
Gambar 14. Perbandingan visualisasi kontur temperatur dengan variasi Reynolds number pada posisi y/h=1 dan α=45⁰.
Gambar 15. Perbandingan visualisasi kontur temperatur dengan variasi Reynolds number pada posisi y/h=1 dan α=90⁰.
Gambar  diambil  pada plane y/h=0,5. Pada Gambar  13 terlihat adanya peningkatan temparatur di sepanjang sumbu z hingga di outline.  Pada  α=30⁰  terlihat  bahwa  kenaikan  temperatur  bisa mencapai  1,19⁰  C.  Sedangkan pada sudut lainnya,  gambar  14 dan 15,  kenaikan  temperatur  bisa mencapai  1,3⁰C.  Kemudian juga untuk area terjadinya kenaikan  temperatur,  area tersempit berada pada model  dengan sudut α=30⁰.  Semakin  curam  sudut inclined step maka semakin  jauh dan luas area yang mengala mi kenaikan temperatur.

## B. Panjang Reattachment
Setelah  dilakukan  simulasi  akan  didapatkan  gambaran mengenai  pola  aliran  separasi  yang  jarak  titik reattachment dengan  step  dapat  diperkirakan  dengan  melihat  vector  lapis batas.  Gambar  16  merupakan  gambaran resirk ulasi apabila terdapat  penghalang  aliran  berupa  gedung.  Karena  adanya penghalang  di  bagian upstream membuat  aliran  fluida terganggu  dan  menciptakan  ruang  yang  lebih  besar  untuk terjadinya  separasi sehingga titik  separasi menjadi  lebih  dekat dari step.
Gambar 16. Resirkulasi aliran dengan variasi sudut incined step pada posisi x=0,75 (Re=6×10⁴).
Gambar  16  diambil  pada  x=0.75  atau  saat aliran  melewa t i tengah gedung. Selanjutnya panjang reattachment dapat terlihat pada tabel 2.
T abel 2. Perbandingan panjang reattachment pada Re=6x10⁴. Variasi α Variasi z/h X=0 X=0,75 X=3 30⁰ 5,762 4,534 5,1 45⁰ 6,504 5,321 5,667 90⁰ 6,854 6,05 6,206
Perbedaan  panjang  reattachment  dikarenakan  step  α=30⁰ memiliki  ruang kosong yang lebih sempit dari  pada step α=45⁰ dan step α=90⁰ sehingga ruang untuk terjadinya wak e -pun lebih sempit.  Hal  ini juga terjadi  pada aliran  dengan Re=2×10⁴  yang memiliki  pola  yang  mirip.  Perbedaan  hanya  terdapat  pada panjang reattachment yang ditampilkan  pada tabel 3.
T abel 3.
Perbandingan panjang reattachment pada Re=2×10⁴. Variasi α Variasi z/h X=0 X=0,75 X=3 30⁰ 5,186 4,459 10,7 45⁰ 5,67 5,01 10,8 90⁰ 6,157 6,206 10,9
C. Coefficient of Pressure (Cp)
Coefficient of Pressure (Cp) merupakan  salah satu bilangan nondimensional  yang dapat mempresentasikan  seberapa besar fluida menghasilkan tekanan tanpa tergantung kecepatan fluida.
Gambar 17. Perbandingan grafik Cf darah downstream dengan variasi sudut inclined step pada x/h=0,75 (Re=2×10⁴).
Dari Gambar  17 terlihat bentuk grafik yang sama. Pada posisi x/h=0,75,  aliran  fluida  dari  ketiga  model  pada mulanya  akan mengalami  percepatan yang ditandai dengan penurunan grafik akibat favorable  pressure  gradient. Setelah  itu  mengala mi perlambatan  yang  ditandai  dengan  kenaikan  sangat  tajam diakibatkan  oleh  aliran  yang tidak  mampu  melawan adverse pressure gradient pada saat penampang mengembang  dengan tiba-tiba  sehingga aliran  terseparasi. Separasi  akan  berangsur normal  dan akhirnya akan hilang, ini disebabkan karena adanya
-0.3 -0.25 -0.2 -0.15 -0.1 -0.05 0 0.05 0.1 0.15 0.2 0 4 8 12 16 20 24 28 32 36 Cp z/h 30 45
(⁰C) (⁰C)
transfer energi  antara aliran  yang berada di daerah yang tidak terseparasi  menuju  aliran  yang  berada  didaerah  terseparasi sehingga  aliran  yang  terseparasi  akan  berangsur-angsur kembali  normal.  Pengambilan  data  diambil  pada  bagian downstream area dan terletak  perhimpit  pada dinding bawah. Karena  nilai  cp  tidak  bergantung  terhadap  ukuran  maupun kecepatan fluida  maka  hasil pengambilan  data pada Re=6× 10⁴ akan  menghasilkan  grafik  Cp  yang sama  persis dengan garfik Re=2×10⁴.  Hal  tersebut dikarenakan  kedua  variasi  kecepatan menggunakan bentuk domain yang sama.
Dari  grafik terlihat bahwa adanya variasi sudut inclined step berpengaruh terhadap tekanan statis. Dari ketiga  variasi  sudut step  terlihat  bahwa  model  yang  memiliki  α=30⁰  dan  α=45⁰ memiliki  grafik  yang  nyaris  berhimpitan.    Setelah  terjadi percepatan, domain dengan  α=30⁰  dan  α=45⁰  memilik i percepatan yang lebih  tinggi  dari  pada  α=90⁰.  Hal  ini  karena kemiringan  step membuat  lintasan lebih  sempit  sehingga nilai cp lebih rendah dan terjadi  percepatan yang lebih tinggi.

## D.Coefficient of Friction (Cf)
Sama  seperti coefficient  of  pressure  (Cp) , coefficient  of friction  (Cf) merupakan  parameter  nondimensional.  Dengan kata lain  besarnya Cf  tidak  terpengaruh oleh kecepatan fluida. Besarnya sk in friction  coefficient (Cf)  erat  kaitannya  dengan besarnya  tegangan  geser  dinding  ( 𝜏𝜔 ).  Geometri back ward facing step identik dengan perubahan penampang yang melebar sehingga aliran  yang mengalir  didalamnya  akan  memperole h pengaruh dari adverse pressure gradient yang besar.
Gambar 18. Perbandingan grafik Cf daerah downstream dengan variasi sudut inclined step pada x/h=0,75 (Re=2×10⁴).
Turunnya  nilai  Cf  dimulai  dari  nol  mengindikasika n terjadinya  separasi aliran  yang dipengaruhi adverse  pressure
gradient sehingga terjadi  penebalan boundary layer, akibatnya aliran  didekat diverging  wall akan  mengalami defisit momentum  sehingga  terjadi  separasi.  Kemudian  berangsur- angsur kecepatan aliran  kembali  normal  diikuti  oleh kenaikan nilai  Cf.  Gambar  18  menunjukan  grafik  Cp  pada  daerah downstream yang diukur dari step hingga ke outlet. Dari  grafik Cf  didapatkan  bahwa  nilai  panjang reattachment  yang paling kecil  berada di plane x=0  dengan nilainya adalah 4 pada sudut α=30⁰  dan nilai  reattachment yang paling besar yaitu pada jarak sudut α=90⁰ dengan nilainya  adalah 6,5.

## IV. KESIMPULAN/ RINGKASAN
Berdasarkan  berbagai  analisa,  kecepatan,  temperatur, Reattachment, Cp, dan Cf,  setiap variasi inclined step memilik i pengaruh  yang  berbeda.  Grafik dan  kontur  kecepatan memperlihatkan  semakin  sempit  sudut inclined step makan semakin  sempit  pula  area  terjadi back flow pada  adaerah downstream .  Grafik  dan  kontur  temperatur  memperlihatka n , baik  Re=2×10⁴  maupun  Re=6×10⁴,  pada  step  α=90⁰  terjadi kenaikan temperatur yang spesifik hingga mencapai outlet. Hal ini  karena area wak e yang lebih luas dan lambatnya  kecepatan fluida  sehingga  penumpukan  fluida  panas  yang  juga  lebih banyak  pada  daerah downstream. Mempertimbangkan  hasil simulasi,  pembangunan  pada  daerah downstream sebaiknya dilakukan  pada  derah  yang  memiliki  sudut inclined  step kecil.atau  landai.  Pembangunan  juga  bisa  dilakukan  tidak berada  tepat  dibawah  bangunan  yang  mengeluarkan  flu ida dengan tempsaratur tinggi.  Selain  itu,  simulasi  dengan variasi jarak  gedung pada sisi upstream dengan step atau jarak  antar gedung diperlukan  untuk mendapat  gambaran  fenomena  pada kondisi yang lain.

## DAFTAR PUSTAKA
[1] Y. T oparlar and et al, Renewable and Sustainable Energy Reviews . Netherlands: Eindhoven University of Technology, 2017.
[2] L. Biao, “Evaluation of CFD Simulation Using Various Turbulence Models  for  Wind  Pressure on  Buildings  Based  on Wind  T unnel Experiments,” China, 2015.
[3] P.  Louda,  “Numerical  Simulation  of  Flows  Over  2D  and  3D Backward-Facing Inclined Steps,” Czech Republic, 2013.
-0.0015 0.0005 0.0025 0 4 8 12 16 20 24 28 32 36 Cf z/h 3 0