ANALISIS ALIRAN DAN PERPINDAHAN PANAS FLUIDA SISKO DALAM KEADAAN STEDI NURI ANGGI NIRMALASARI 1207 100 017
BAB I PENDAHULUAN LATAR BELAKANG RUMUSAN MASALAH BATASAN MASALAH TUJUAN MANFAAT
LATAR BELAKANG Fluida Sisko digunakan dalam bidang industri dan teknik Fluida Sisko merupakan fluida non-Newtonian Sulit memprediksi perilaku fluida tersebut Menentukan model matematika dari kecepatan aliran dan perpindahan panasnya Menyelesaikan kedua model matematika tersebut secara numerik Menampilkan penyelesaian yang didapat dalam bentuk grafik
Menganalisis bagaimana profil kecepatan aliran dan perpindahan panas berdasarkan grafik
RUMUSAN MASALAH Bagaimana model matematika dari kecepatan aliran dan perpindahan panas fluida sisko dalam pipa
visualisasi profil kecepatan dan perpindahan panas fluida sisko didalam pipa dalam bentuk grafik.
Bagaimana penyelesaian numerik dari model matematika kecepatan aliran dan perpindahan panas fluida sisko dalam pipa
BATASAN MASALAH Tipe aliran fluida sisko yang mengalir dalam pipa adalah seragam stedi. Model matematika dari permasalahan tersebut diselesaikan secara numerik dengan metode beda hingga pusat. Diasumsikan pipa yang digunakan adalah pipa lurus dengan panjang (L). Penampang pipa berupa silinder dengan diameter (D) Luas penampang pipa adalah konstan
TUJUAN Menurunkan model matematika dari kecepatan aliran dan perpindahan panas fluida sisko dalam pipa Menyelesaikan model matematika kecepatan aliran dan perpindahan panas fluida sisko dalam pipa secara numerik visualisasi profil kecepatan dan perpindahan panas fluida sisko didalam pipa dalam bentuk grafik.
MANFAAT BAGI BIDANG TEKNIK DAN INDUSTRI diharapkan dapat memberikan informasi tentang profil kecepatan aliran dan perpindahan panas fluida sisko dalam pipa sehingga dapat digunakan sebagai bahan pengetahuan untuk mengembangkan aplikasinya
BAGI BIDANG MATEMATIKA untuk mengetahui bagaimana peran metode beda hingga dalam menyelesaikan masalah yang terjadi di bidang teknik
BAB II TINJAUAN PUSTAKA FLUIDA Fluida merupakan zat yang berubah bentuk secara kontinu, bila terkena tegangan geser. Fluida terdiri dari 2 macam yaitu fluida cair ( tak mampu-mampat) dan fluida gas (mampu-mampat). Pada fluida cair terdapat viskositas, yaitu sifat dari fluida untuk melawan tegangan geser pada waktu mengalir. Untuk fluida pada umumnya, tegangan geser dan lau regangan geser (gradien kecepatan) dapat dikaitkan dalam suatu hubungan sebagai berikut, (Munson, 2004): •
𝜏 = 𝜇𝛾
dimana : 𝜏 = tegangan geser 𝜇 = kekentalan (viskositas) 𝛾 = laju regangan geser
Karakteristik Fluida Fluida Cair Non-Newtonian
Newtonian
Shear thinning
Shear tickening
Bingham plastic
𝜏 = 𝜇𝛾+a
𝜏 = 𝜇𝛾 𝜏 >, 𝜇<
𝜏 >, 𝜇>
Fluida Sisko Fluida Sisko termasuk dalam karakteristik Bingham Plastic, yang pada beberapa kasus dialirkan dalam pipa annulus, yaitu pipa yang terdiri dari pipa luar dan pipa dalam dengan pusat jari-jari adalah sama. Tensor teganga fluida sisko sbb (M.Khan, 2010):
𝑻 = −𝑝𝑰 + 𝑺
dengan 𝑆 = 𝑎 + 𝑏 𝐴1 = 𝐿 + 𝐿𝑇 ,
1 𝑡𝑟 2
𝑛−1
𝐴12
𝐴1
𝐿 = 𝑔𝑟𝑎𝑑 𝑉
Dimana: 𝑻 : Tensor tegangan 𝑝 ∶ tekanan 𝑺 ∶ tegangan geser pada fluida sisko a, b : Parameter material n : Power index, termasuk sebagai parameter material 𝑉 : Kecepatan, 𝑉=v(r) T : Temperatur, T=T(r)
Koordinat Polar Tempat kedudukan sebuah titik ditunjukkan oleh koordinat-koordinat 𝑟, 𝜃 dan 𝑧. 𝑟 = jarak radial dari sumbu- 𝑧 𝜃 = sudut yang diukur dari garis sejajar sumbu- x 𝑧 = koordinat sepanjang sumbu-z Misal pada kecepatan, komponen-komponennya adalah: 𝑣𝑟 = kecepatan radial 𝑣𝜃 = kecepatan tangensial 𝑣𝑧 = kecepatan aksial Sehingga kecepatan pada sebuah titik, dinyatakan:
𝑽 = 𝒗𝒓 𝒆𝒓 + 𝒗𝜽 𝒆𝜽 + 𝒗𝒛 𝒆𝒛 Dimana: 𝒆𝒓 = vektor arah 𝒓 𝒆𝜽 = vektor arah 𝜽 𝒆𝒛 = vektor arah 𝒛
Persamaan Kontinuitas Hukum kekekalan massa: 𝜕 𝜕𝑡
𝜌 𝑑∀ +
𝜌. 𝐴𝑉. 𝑛 = 0
Dimana: • 𝜌 : kerapatan fluida • 𝑉 : komponen kecepatan fluida yang tegak lurus bidang 𝐴 Pada aliran steady state Dalam keadaan steady 𝜌. 𝐴𝑉. 𝑛 = 0 𝛻. 𝜌. 𝐴𝑉. 𝑛 = 𝛻. 0 𝛻. 𝑉 = 0 ………..Persamaan kontinuitas Dalam koordinat polar silinder: 1 𝜕(𝑟𝑣𝑟 ) 𝑟 𝜕𝑟
+
1 𝜕𝑣𝜃 𝑟 𝜕𝜃
+
𝜕𝑣𝑧 =0 𝜕𝑧
Persamaan Momentum Linier 𝜕 𝜕𝑡
𝜌 𝑉 𝑑∀ +
𝜌. 𝐴𝑉 2 . 𝑛 =
𝑔𝑎𝑦𝑎 − 𝑔𝑎𝑦𝑎 𝑦𝑎𝑛𝑔 𝑏𝑒𝑘𝑒𝑟𝑗𝑎
Berdasarkan pergerakan fluida, dijelask+an oleh persamaan Navier-Stokes, yaitu: 𝑑𝑉 𝜌 = 𝛻. 𝑻 𝑑𝑡 Pada fluida Newtonian yang mengalir dalam pipa pada arah-z (sejajar dinding), persamaan diatas menjadi: 𝜕𝑣𝑧 𝜕𝑣𝑧 𝑣𝜃 𝜕𝑣𝜃 𝜕𝑣𝑧 𝜕𝑝 1 𝜕 𝜕𝑣𝑧 1 𝜕 2 𝑣𝑧 𝜕 2 𝑣𝑧 ρ +𝑣𝑟 + +𝑣𝑧 =− + ρ𝑔 + 𝜇 𝑟 + 2 + 𝜕𝑡 𝜕𝑟 𝑟 𝜕𝜃 𝜕𝑧 𝜕𝑧 𝑟 𝜕𝑟 𝜕𝑟 𝑟 𝜕𝜃 2 𝜕𝑧 2
Aliran fluida dalam pipa annulus Kecepatan aliran dalam pipa diasumsikan sebagai berikut: 𝜕𝑣𝑧 =0 𝜕𝑧 𝜕𝑣 steady 𝑧 =0, 𝜕𝑡
Aliran sejajar dengan dinding sehingga 𝑣𝑟 = 0 dan 𝑣𝜃 = 0 akibatnya (berdasarkan persamaan kontinuitas). Selain itu pada keadaan dengan demikian persamaan Navier-Stokes menjadi: 1𝜕 𝜕𝑝 𝜕𝑣𝑧 0=− +𝜇 𝑟 𝑟 𝜕𝑟 𝜕𝑧 𝜕𝑟 Dengan kondisi batas 𝑣𝑧 = 0 pada 𝑟 = 𝑟0 𝑣𝑧 = 𝑣1 pada 𝑟 = 𝑟1 Dimana : 𝑟0 adalah jari-jari silinder dalam 𝑣1 dan 𝑟1 merupakan kecepatan dan jari-jari silinder luar
Persamaan Distribusi Panas pada aliran fluida Persamaan secara umum (Lienhard, 2005): 𝜕𝑇 𝜌𝐶𝑝 = 𝛻. 𝑘𝛻𝑇 + 𝑞 𝜕𝑡 Sedangkan pada fluida sisko dinyatakan sebagai berikut: 𝜌𝐶𝑝
𝜕𝑇 = 𝑻. 𝑳 − 𝑑𝑖𝑣 𝑞 𝜕𝑡
Dimana: • 𝜌 adalah densitas, • 𝐶𝑝 adalah kapasitas panas pada tekanan konstan, • 𝑞 adalah fluks panas yang persamaannya ditentukan sebagai berikut: 𝑞 = −𝑘 𝑔𝑟𝑎𝑑𝑇
Metode Beda Hingga Deret Taylor: • 𝑢(𝑥0 + ℎ) = 𝑢 𝑥0 + ℎ. 𝑢 𝑥0 +
ℎ2 𝑢 2!
• 𝑢 𝑥0 − ℎ = 𝑢 𝑥0 − ℎ. 𝑢 𝑥0 +
ℎ2 𝑢 2!
𝑥0 + … +
ℎ𝑛 𝑢 𝑛!
𝑥0 − ⋯ +
𝑥0 + 𝑅𝑛
ℎ𝑛 𝑢 𝑛!
𝑥0 + 𝑅𝑛
Didapat pendekatan turunan pertama metode beda hingga pusat: 𝑑𝑢 𝑥0 𝑢 𝑥0 + ℎ − 𝑢 𝑥0 − ℎ ≈ 𝑑𝑥 2ℎ Pendkatan turunan kedua: 𝑑𝑢2 𝑥0 𝑢 𝑥0 + ℎ − 2𝑢 𝑥0 + 𝑢 𝑥0 − ℎ ≈ 𝑑𝑥 2 ℎ2
BAB III PROSEDUR KERJA Studi Literatur
Model kecepatan aliran
Model perpindahan panas
Persamaan diferensial, n=0
Persamaan diferensial, n=1
Persamaan diferensial, n=0
Persamaan diferensial, n=1
Didapat kecepatan secara numerik
Didapat kecepatan secara numerik
Didapat temperatur secara numerik
Didapat temperatur secara numerik
Setiap penyelesaian divisualisasikan dalam bentuk grafik
Analisis profil kecepatan berdasarkan nilai b yang bervariasi
Analisis distribusi temperatur berdasarkan nilai b dan Br yang bervariasi
Membandingkan profil kecepatan dan distribusi temperatur antara Fluida Sisko dan Fluida Newtonian
BAB IV PEMODELAN DAN PENYELESAIAN NUMERIK • Model didapat dengan menurunkan persamaan kontinuitas, momentum linier, dan perpindahan panas yang dipengaruhi oleh gradient kecepatan dan gaya-gaya pada fluida sisko 4.1 Model kecepatan aliran
Persamaan aliran dalam pipa 0=−
𝜕𝑝 𝜕𝑧
+𝜇
1 𝜕 𝑟 𝜕𝑟
𝑟
𝜕𝑣𝑧 𝜕𝑟
𝜕𝑝 1 𝜕 𝜕𝑣𝑧 = 𝑟. 𝜇 𝜕𝑧 𝑟 𝜕𝑟 𝜕𝑟 𝜕𝑣
𝜇 𝑧 = 𝜏𝑟𝑧 (tegangan 𝜕𝑟 geser fluida newtonian)
𝜕𝑝 1 𝑑 = 𝑟𝑺𝑟𝑧 𝜕𝑧 𝑟 𝑑𝑟
𝑺𝑟𝑧 (tegangan geser fluida sisko)
𝜕𝑝 1 𝜕𝑝 𝜕𝑝 𝑑𝑝 𝜕𝑝 𝜕𝑝 𝜕𝑝 𝑑𝑟 + 𝑑𝜃 + 𝑑𝑧 = Dimana, = = 0, sehingga 𝜕𝑟 𝜕𝜃 𝜕𝑟 𝑟 𝜕𝜃 𝜕𝑧 𝑑𝑧 𝜕𝑧 𝑛−1 𝑑𝑝 1 𝑑 dengan 𝑑𝑣 𝑑𝑣 Persamaan Kecepatan Aliran : = 𝑟𝑺 𝑺= 𝑎+𝑏 𝑑𝑧 𝑟 𝑑𝑟 𝑑𝑟 𝑑𝑟 𝑑𝑝 =
Persamaan diferensial Kecepatan Aliran fluida sisko 𝑑𝑝 1 𝑑 𝑑𝑣 = 𝑟 𝑎+𝑏 𝑑𝑧 𝑟 𝑑𝑟 𝑑𝑟
Dimana : 𝑑𝑝 : gradient kecepatan 𝑑𝑧
𝑟 𝑎, 𝑏 𝑣
: jari-jari penampang pipa : parameter material : fungsi kecepatan terhadap jari-jari
𝑛−1
𝑑𝑣 𝑑𝑟
4.2 Model Distribusi panas
Persamaan Distribusi Panas: 𝑑𝑇 𝜌𝐶𝑝 = 𝑻. 𝑳 − 𝑑𝑖𝑣 𝑞 𝑑𝑡 𝑞 = −𝑘 𝑔𝑟𝑎𝑑 𝑇 𝑻 = −𝑝𝑰 + 𝑺 𝐿 = 𝑔𝑟𝑎𝑑 𝑉
𝜌𝐶𝑝
𝑑𝑇 𝑑𝑣 = (−𝑝𝑰 + 𝑺). + 𝑘 𝛻2𝑇 𝑑𝑟 𝑑𝑡 Steady
(−𝑝𝑰 + 𝑺).
𝑑𝑣 𝑑𝑟
+ 𝑘 𝛻 2 𝑇= 0 T=T(r)
Persamaan diferensial Distribusi panas fluida sisko 𝑘 𝑑 𝑑𝑇 𝑑𝑣 𝑟 + 𝑎+𝑏 𝑟 𝑑𝑟 𝑑𝑟 𝑑𝑟
𝑛−1
𝑑𝑣 𝑑𝑟
2
=0
Model Kecepatan aliran dan Perpindahan panas non-dimensional • Variabel-variabel non-dimensional 𝑟 𝑟∗ = , 𝑟0 𝑝∗ = 𝑃𝑟 =
𝑧 𝑧∗ = , 𝑧0
𝑝
𝑎𝑣1 𝑟0
𝑎𝑐𝑝 𝑘
,
𝑣 𝑣∗ = , 𝑣1 𝑇 − 𝑇0 𝑇∗ = , 𝑇1 − 𝑇0
𝐵𝑟 = 𝑃𝑟 𝐸𝑐
,
Model matematika Distribusi panas fluida sisko
Model matematika Kecepatan aliran fluida sisko
𝑑𝑝 𝑑𝑧
=
𝑑 𝑑𝑟
1+𝑏
𝑑𝑣 𝑛−1 𝑑𝑟
𝑑𝑣 𝑑𝑟
+
1 + 𝑟
𝑛−1
𝑏 𝑣1 𝑏∗ = , 𝑎 𝑟0 𝑣1 2 𝐸𝑐 = , 𝑐𝑝 𝑇0 − 𝑇1
1+𝑏
𝑑𝑣 𝑛−1 𝑑𝑟
𝑑𝑣 𝑑𝑟
1𝑑 𝑑𝑇 𝑑𝑣 𝑟 + 𝐵𝑟 1 + 𝑏 𝑟 𝑑𝑟 𝑑𝑟 𝑑𝑟
𝑛−1
𝑑𝑣 𝑑𝑟
2
=0
Penyelesaian Numerik model matematika kecepatan aliran Akan dibandingkan bagaimana profil kecepatan dengan power index n=0 dan n=1 𝑑 2 𝑣 1 𝑑𝑣 𝑏 𝑑𝑝 + + = 𝑑𝑟 2 𝑟 𝑑𝑟 𝑟 𝑑𝑧
n=0
Skema numerik dengan metode beda hingga pusat 𝑣𝑖+1 − 2𝑣𝑖 + 𝑣𝑖−1 1 𝑣𝑖+1 + 𝑣𝑖−1 𝑏 𝑑𝑝 + + = 𝑟𝑖 𝑟𝑖 𝑑𝑧 (∆𝑟)2 2∆𝑟 Dengan : 𝑟𝑖 = 1 + 𝑖∆𝑟, dimana ∆𝑟 = 𝑑 − 1 𝑁 Atau bisa ditulis 1 1 2 1 1 𝑑𝑝 + 𝑣 − 𝑣 + − 𝑣 = − 𝑏. 𝑞𝑖 ∆𝑟 2 2(1 + 𝑖∆𝑟)∆𝑟 𝑖+1 ∆𝑟 2 𝑖 ∆𝑟 2 2 1 + 𝑖∆𝑟 ∆𝑟 𝑖−1 𝑑𝑧
Didapat skema numerik untuk i= 1,2,3,…,N-1 4𝑝2 + 𝑝𝑞𝑖 𝑣𝑖+1 − 8𝑝2 𝑣𝑖 + 4𝑝2 − 𝑝𝑞𝑖 𝑣𝑖−1 =
Dengan: 𝑝=
1 2∆𝑟
dan
𝑞𝑖 =
1 1+𝑖∆𝑟
𝑑𝑝 − 𝑏. 𝑞𝑖 𝑑𝑧
Didapat matriks penyelesaian sebagai berikut: − 8𝑝2 4𝑝2 − 𝑝𝑞2 0 0 ⋱ 0
0 0 4𝑝2 + 𝑝𝑞1 0 2 2 4𝑝 + 𝑝𝑞 0 2 0 − 8𝑝 2 2 0 4𝑝 + 𝑝𝑞3 4𝑝2 − 𝑝𝑞3 − 8𝑝 2 4𝑝 + 𝑝𝑞4 4𝑝2 − 𝑝𝑞4 − 8𝑝2 0 ⋱ ⋱ ⋱ ⋱ 0 0 0 0 𝑑𝑝/𝑑𝑧 𝑞1 𝑑𝑝/𝑑𝑧 𝑞2 𝑞3 𝑑𝑝/𝑑𝑧 −𝑏 𝑞 = 4 𝑑𝑝/𝑑𝑧 ⋮ ⋮ 𝑞𝑛−1 𝑑𝑝/𝑑𝑧 − 4𝑝2 + 𝑝𝑞𝑛−1
… 0 ⋱ 0 0 ⋱ 0 ⋱ ⋱ ⋱ 2 … 4𝑝 − 𝑝𝑞𝑛−1
0 0 0 0 ⋱ − 8𝑝2
𝑣1 𝑣2 𝑣3 𝑣4 ⋮ 𝑣𝑛−1
𝑑2 𝑣 1 𝑑𝑣 1 𝑑𝑝 + = 𝑑𝑟 2 𝑟 𝑑𝑟 1 + 𝑏 𝑑𝑧
n=1
Dengan cara dan definisi yang sama dengan yang dikenakan pada n=0, didapat skema numerik : 4𝑝2 + 𝑝𝑞𝑖 𝑣𝑖+1 − 8𝑝2 𝑣𝑖 + 4𝑝2 − 𝑝𝑞𝑖 𝑣𝑖−1 =
1 𝑑𝑝 . 1 + 𝑏 𝑑𝑧
Dan dinyatakan dalam bentuk matriks sebagai berikut: − 8𝑝2 4𝑝2 − 𝑝𝑞2 0 0 ⋱ 0
•
•
=
1 1+𝑏
0 4𝑝2 + 𝑝𝑞1 2 4𝑝 + 𝑝𝑞2 − 8𝑝2 2 4𝑝2 − 𝑝𝑞3 − 8𝑝 4𝑝2 − 𝑝𝑞4 0 ⋱ ⋱ 0 0
𝑑𝑝/𝑑𝑧 0 𝑑𝑝/𝑑𝑧 0 0 𝑑𝑝/𝑑𝑧 − 0 𝑑𝑝/𝑑𝑧 ⋮ ⋮ 2 + 𝑝𝑞 4𝑝 𝑛−1 𝑑𝑝/𝑑𝑧
0 0 0 0 2 0 4𝑝 + 𝑝𝑞3 2 4𝑝 + 𝑝𝑞4 − 8𝑝2 ⋱ ⋱ 0 0
… 0 ⋱ 0 0 ⋱ 0 ⋱ ⋱ ⋱ 2 … 4𝑝 − 𝑝𝑞𝑛−1
0 0 0 0 ⋱ − 8𝑝2
𝑣1 𝑣2 𝑣3 𝑣4 ⋮ 𝑣𝑛−1
Penyelesaian Numerik model matematika distribusi panas 1𝑑 𝑑𝑣 𝑑𝑇 𝑟 + 𝐵𝑟 +𝑏 𝑟 𝑑𝑟 𝑑𝑟 𝑑𝑟
n=0
𝑑𝑣 =0 𝑑𝑟
Skema numerik dengan metode beda hingga pusat, untuk i=1,2,3,…,N-1 4𝑝2 + 𝑝𝑞𝑖 𝑇𝑖+1 − 8𝑝2 𝑇𝑖 + 4𝑝2 − 𝑝𝑞𝑖 𝑇𝑖−1 = − 𝐵𝑟
𝑣𝑖+1 −𝑣𝑖−1 2 2∆𝑟
+𝑏
𝑣𝑖+1 −𝑣𝑖−1 2∆𝑟
Dengan 𝑣 bergantung pada penyelesaian sebelumnya, selanjutnya penyelesaian dalam bentuk matriks
•
•
− 8𝑝2 4𝑝2 − 𝑝𝑞2 0 0 ⋱ 0
= −𝐵𝑟
0 4𝑝2 + 𝑝𝑞1 2 4𝑝 + 𝑝𝑞2 − 8𝑝2 2 4𝑝2 − 𝑝𝑞3 − 8𝑝 4𝑝2 − 𝑝𝑞4 0 ⋱ ⋱ 0 0
𝑝(𝑣2 ) 2 𝑝(𝑣3 − 𝑣1 ) 2 𝑝(𝑣4 − 𝑣2 ) 2 − 𝐵 . 𝑏 𝑟 𝑝(𝑣5 − 𝑣3 ) 2 ⋮ 𝑝(1 − 𝑣𝑛−2 ) 2
0 0 0 0 2 0 4𝑝 + 𝑝𝑞3 2 4𝑝 + 𝑝𝑞4 − 8𝑝2 ⋱ ⋱ 0 0
𝑝(𝑣2 ) 0 𝑝(𝑣3 − 𝑣1 ) 0 0 𝑝(𝑣4 − 𝑣2 ) − 0 𝑝(𝑣5 − 𝑣3 ) ⋮ ⋮ 2 4𝑝 + 𝑝𝑞𝑛−1 𝑝(1 − 𝑣𝑛−2 )
… 0 ⋱ 0 0 ⋱ 0 ⋱ ⋱ ⋱ 2 … 4𝑝 − 𝑝𝑞𝑛−1
0 0 0 0 ⋱ − 8𝑝2
𝑇1 𝑇2 𝑇3 𝑇4 ⋮ 𝑇𝑛−1
n=1
𝑑𝑇 1 𝑑 𝑟 + 𝐵𝑟 1 + 𝑏 𝑟 𝑑𝑟 𝑑𝑟
𝑑𝑣 𝑑𝑟
2
=0
Skema numerik dengan metode beda hingga pusat, untuk i=1,2,3,…,N-1
4𝑝2 + 𝑝𝑞𝑖 𝑇𝑖+1 − 8𝑝2 𝑇𝑖 + 4𝑝2 − 𝑝𝑞𝑖 𝑇𝑖−1 𝑣𝑖+1 − 𝑣𝑖−1 2 = −𝐵𝑟 1 + 𝑏 2∆𝑟 𝑣 bergantung pada penyelesaian model kecepatan aliran untuk n=1, sehingga didapat matriks sebagai berikut:
•
•
− 8𝑝2 4𝑝2 − 𝑝𝑞2 0 0 ⋱ 0
0 4𝑝2 + 𝑝𝑞1 2 + 𝑝𝑞 2 4𝑝 2 − 8𝑝 2 2 − 8𝑝 4𝑝 − 𝑝𝑞3 4𝑝2 − 𝑝𝑞4 0 ⋱ ⋱ 0 0
0 0 0 0 0 4𝑝2 + 𝑝𝑞3 4𝑝2 + 𝑝𝑞4 − 8𝑝2 ⋱ ⋱ 0 0
𝑝(𝑣2 ) 2 0 𝑝(𝑣3 − 𝑣1 ) 2 0 2 0 = −𝐵𝑟 (1 + 𝑏) 𝑝(𝑣4 − 𝑣2 ) − 2 0 𝑝(𝑣5 − 𝑣3 ) ⋮ ⋮ 2 4𝑝 + 𝑝𝑞𝑛−1 𝑝(1 − 𝑣𝑛−2 ) 2
… 0 ⋱ 0 0 ⋱ 0 ⋱ ⋱ ⋱ 2 … 4𝑝 − 𝑝𝑞𝑛−1
0 0 0 0 ⋱ − 8𝑝2
𝑇1 𝑇2 𝑇3 𝑇4 ⋮ 𝑇𝑛−1
BAB V VISUALISASI DAN PEMBAHASAN Algoritma program: • Mendefinisikan parameter-parameter yang dibutuhkan. • Mendefinisiskan kondisi batas yang telah ditentukan pada bab 4. • Memasukkan kondisi batas ke dalam skema numerik penyelesaian model matematika kecepatan aliran, yaitu persamaan . • Skema numerik yang berupa matrik tridiagonal diselesaikan, sehingga didapat nilai kecepatan pada titik-titik sepanjang jari-jari pipa. • Selanjutnya nilai kecepatan dimasukkan ke dalam skema numerik penyelesaian model matematika perpindahan panas yaitu persamaan • Selanjutnya dihitung nilai error berdasarkan penyelesaian numerik dan penyelesaian eksak.
Grafik kecepatan aliran n=0 𝑑𝑝
Gradien tekanan = −0.4. dengan variasi parameter material 𝑏 = 0 untuk 𝑑𝑧 fluida Newtonian, dan 𝑏 = 0.2, 0.5 dan 0.8 untuk fluida Sisko. Berdasarkan grafik akan dianalisis distribusi kecepatan untuk banyak pendiskritan dan jarijari silinder luar yang berbeda-beda, sebagai berikut:
𝑁 = 50 dan 𝑑 = 10
𝑁 = 150 dan 𝑑 = 10
𝑁 = 150 dan 𝑑 = 20
Pada gambar (5.2) terlihat kecepatan disekitar jari-jari masih belum mendekati satu, sedangkan pada gambar (5.3) kecepatan mendekati 1 pada jari-jari dsekitar 20 semakin besar nilai b maka kecepatan aliran fluida semakin besar
Grafik Kecepatan Aliran
n=1 𝑑𝑝
Gradien tekanan = −0.4. dengan variasi parameter material 𝑏 = 0 untuk 𝑑𝑧 fluida Newtonian, dan 𝑏 = 0.2, 0.5 dan 0.8 untuk fluida Sisko. Berdasarkan grafik akan dianalisis distribusi kecepatan untuk banyak pendiskritan dan jarijari silinder luar yang berbeda-beda, sebagai berikut:
𝑁 = 50 dan 𝑑 = 10
𝑁 = 150 dan 𝑑 = 10
𝑁 = 150 dan 𝑑 = 20
Berbeda dengan distribusi kecepatan dengan power index 𝑛 = 0, dari grafik dapat disimpulkan bahwa untuk power index 𝑛 = 1, kecepatan aliran pada fluida sisko lebih kecil dibandingkan fluida Newtonian, atau dengan kata lain semakin besar nilai b, maka kecepatan aliran semakin kecil
Grafik Distribusi Panas n=0 𝑑𝑝
Gradien tekanan = −0.4 dan bilangan Brinkman 1, dengan variasi 𝑑𝑧 parameter material 𝑏 = 0 untuk fluida Newtonian, dan 𝑏 = 0.2, 0.5 dan 0.8 untuk fluida Sisko. Berdasarkan grafik akan dianalisis distribusi panas untuk banyak pendiskritan dan jari-jari silinder luar yang berbeda-beda, sebagai berikut:
𝑁 = 50 dan 𝑑 = 10
𝑁 = 150 dan 𝑑 = 10
𝑁 = 150 dan 𝑑 = 20
Dengan memberikan nilai parameter material yang berbeda-beda terlihat bahwa semakin besar nilai b yang diberikan maka distribusi temparatur semakin besar, yang berarti pada distribusi panas, temperatur fluida sisko lebih besar dari temperatur fluida Newtonian.
Grafik Distribusi Panas n=1 𝑑𝑝
Gradien tekanan = −0.4 dan bilangan Brinkman 1, dengan variasi 𝑑𝑧 parameter material 𝑏 = 0 untuk fluida Newtonian, dan 𝑏 = 0.2, 0.5 dan 0.8 untuk fluida Sisko. Berdasarkan grafik akan dianalisis distribusi panas untuk banyak pendiskritan dan jari-jari silinder luar yang berbeda-beda, sebagai berikut:
𝑁 = 50 dan 𝑑 = 10
𝑁 = 150 dan 𝑑 = 10
𝑁 = 150 dan 𝑑 = 20
Dari grafik distribusi temperatur diatas terlihat bahwa gafik yang dihasilkan tidak berbeda dengan grafik distribusi temperatur untuk power index 𝑛 = 0, berbeda dengan distribusi kecepatan aliran, untuk distribusi panas dengan power index 𝑛 = 0 dan power index 𝑛 = 1, temperataur semakin tinggi untuk nilai parameter material b yang lebih besar.
Analisis distribusi panas berdasarkan bilangan Brinkman Grafik distribusi panas, dengan variasi bilangan Brinkman:
Bilangan Brinkman merupakan bilangan yang mempengaruhi besarnya temperatur pada fluida. Dari grafik distribusi panas pada gambar diatas terlihat bahwa semakin besar bilangan Brinkman yang diberikan, maka temperatur fluida sisko semakin besar.
BAB VI SIMPULAN Kesimpulan Model Kecepatan Aliran
𝐝 𝐝𝐫
Model Perpindahan Panas Dari grafik disimpulkan bahwa
𝐝𝐰 𝟏+𝐛 𝐝𝐫
𝐧−𝟏
𝐝𝐰 𝟏 𝐝𝐰 + 𝟏+𝐛 𝐝𝐫 𝐫 𝐝𝐫
𝟏 𝒅 𝒅𝑻 𝒅𝒘 𝒓 + 𝑩𝒓 𝟏 + 𝒃 𝒓 𝒅𝒓 𝒅𝒓 𝒅𝒓
a.
b.
c.
𝒏−𝟏
𝐧−𝟏
𝒅𝒘 𝒅𝒓
𝐝𝐰 𝐝𝐩 = 𝐝𝐫 𝐝𝐳
𝟐
=𝟎
Untuk Power index n=0, kecepatan aliran fluida sisko lebih besar daripada fluida Newtonian, atau semakin besar nilai b, kecepatan semakin besar, begitu juga dengan temperatur. Untuk power index n=1, semakin besar nilai parameter b, kecepatan aliran semakin kecil, namun temperatur semakin tinggi. Semakin besar bilangan Brinkman yang diberikan, temparatur semakin tinggi
Saran 1. Pada Tugas Akhir ini analisis yang dilakukan menggunakan asumsi bahwa aliran fluida sisko dalam pipa annulus dalam keadaan steady, selanjutnya dapat dikembangkan penelitian untuk menganalisis profil kecepatan aliran dan perpindahan panas fluida sisko dalam pipa dalam keadaan unsteady. 2. Tugas Akhir ini masih bersifat analitis pada tahap pemodelan dan numerik untuk penyelesaiannya, belum ada data laboraturium yang dipakai sebagai pembanding. Diharapkan kedepannya bisa dilakukan uji laboraturium sehingga model tersebut dapat diterapkan di lapangan.
.
DAFTAR PUSTAKA Abdia, Gunaidi. 2006. Matlab Programming. Bandung: Informatika Alfijar, Julian. Mekanika Fluida II. http://alfijar.files.wordpress.com/2008/01/pertemuaniii-dan-iii.ppt-Mirip. Diakses pada tanggal 1 Maret 2011 pukul 11.00 WIB. Arhami, Muhammad dan Desiani, Anita. 2005. Pemrograman Matlab. Yogyakarta: ANDI. Khan, M. et. al. 2010. Steady Flow and Heat Transfer of a Sisko Fluid In Annular Pipe. Journal of Heat and Mass Transfer. 53: 1290-1297. Departmen of Mathematics, Pakistan. Lienhard IV, John H dan Lienhard V, John H. 2005. A Heat Transfer Textbook. University of Houston. USA. Munson, Bruce R. et. al. 2004. Mekanika Fluida. Edisi Keempat Harinaldi dan Budiarso, penerjemah. Jakarta: Erlangga. Terjemahan dari: Fundamental of Fluid Mechanics. Sajid, M and Hayat, T. 2008. Wire Coating Analysis by Withdrawal From A Bath of Sisko Fluid. Journal of Applied Mathematics and Computation. 199: 13-22. Departmen of Mathematics, Pakistan Saragi, Elfrida. Solusi Analitik dan Numerik Konduksi Panas Pada Pembangkit Energi. http://www.batan.go.id/ppin/lokakarya/LKSTN_10/Elfrida-.pdf. Diakses pada tanggal 2 Maret 2011 pukul 12.00 WIB. Siddiq, A.m. et. al. 2009. On Taylor’s Scraping Problem and Flow of A Sisko Fluid. Journal of Mathe matical Modelling and Analysis. 14: 515-529. Department of Mthematics, York Campus, York, PA 17403, USA.
Streeter, Victor L and Wylie, E Benjamin. 1999. Mekanika Fluida. Edisi Delapan Arko Prijono, penerjemah. Jakarta: Erlangga. Terjemahan dari: Fluid Mechanics Sweet, Erik. 2003. Analytical and Numerical Solutions of Differential Equations Arising In Fluid Flow and Heat Transfer Problems. University of Central Florida Orlando, Florida. Ruwanto, Bambang. 2003. Matematika Untuk Fisika dan Teknik. Yogyakarta: Adicita Karya Nusa.
