HASIL DAN PEMBAHASAN Distribusi Suhu dan Kelembaban Udara pada Kandang Sapi Perah Analisis distribusi suhu dan kelembaban udara dilakukan pada saat kandang tidak diisi sapi (kandang kosong). Karakteristik kandang sapi perah yang meliputi dimensi kandang, karakteristik bahan penyusun kandang, kemiringan atap kandang, sifat fisik udara di dalam dan luar kandang dijadikan input pada analisis distribusi suhu dan kelembaban udara dalam kandang sapi perah menggunakan CFD. Output CFD mengenai distribusi udara di dalam kandang ditampilkan dalam bentuk irisan kontur melintang dan vektor kecepatan aliran dan pembentukan suhu udara di dalam kandang dimana suhu udara ditampilkan dalam total temperature, sedangkan kelembaban udara diperoleh dari hasil perhitungan menggunakan persamaan (21) sampai dengan persamaan (24). Data yang digunakan untuk analisis distribusi suhu dan kelembaban udara adalah hasil pengukuran tanggal 16 Juni 2007 yang dipilih mewakili cuaca cerah pada musim kemarau. Pemilihan waktu untuk analisis distribusi suhu dan RH serta simulasi dilakukan berdasarkan kondisi kecepatan angin dan suhu udara lingkungan yang relatif stabil pada waktu tertentu (± 30 detik) sehingga diperoleh aliran udara yang laminer dalam kandang (Lampiran 1). Aliran udara dalam kandang dengan bukaan yang sangat lebar dianggap sebagai aliran udara di atas bidang datar dan laminer apabila nilai bilangan Reynolds kurang dari 500.000 (Cengel, 2003). Berdasarkan pemilihan waktu tersebut, maka diperoleh waktu yang dapat mewakili pagi hari (09:20 WIB), siang hari (13:00 WIB) dan sore hari (15:20 WIB). Radiasi matahari pada pukul 09:20, 13:00 dan 15:20 WIB masing-masing sebesar 396,04 Watt/m2, 506, 57 Watt/m2 dan 317,32 Watt/m2. Kecepatan dan arah angin pada pukul 09:20, 13:00 dan 15:20 WIB masing-masing sebesar 1,0 m/det dari arah depan kandang, 1,0 m/det dari arah kiri kandang dan 1,0 m/det dari arah depan kandang. Suhu udara lingkungan pada pukul 09:20, 13:00 dan 15:20 WIB masing-masing sebesar 28,8; 32,52 dan 31,8oC dan RH lingkungan pada pukul 09:20, 13:00 dan 15:20 WIB masing-masing sebesar 74,5%, 57,6%
dan 56%. Radiasi matahari pada tanggal 16 Juni 2007 dapat dilihat pada Gambar 8, sedangkan suhu udara dan RH lingkungan disajikan pada Gambar 9. Dari Gambar 8 terlihat bahwa radiasi matahari meningkat dari pagi sampai siang hari dan turun pada sore harinya. Naik turunnya radiasi sangat mempengaruhi suhu udara dan RH lingkungan. Suhu udara lingkungan meningkat dari pagi hari sampai siang hari dan mencapai puncaknya pada pukul 13:00 WIB, sedangkan RH tinggi pada pagi hari kemudian turun pada siang hari dan meningkat kembali pada sore hari.
Radiasi matahari (W/m 2)
600 500 400 300 200 100
6: 00 7: 00 8: 00 9: 00 10 :0 0 11 :0 0 12 :0 0 13 :0 0 14 :0 0 15 :0 0 16 :0 0 17 :0 0 18 :0 0
0
Pukul (WIB)
Gambar 8 Radiasi matahari (Watt/m2) pada tanggal 16 Juni 2007
o
Suhu udara ( C) dan RH (%)
100 90
Suhu udara
RH
80 70 60 50 40 30 20 10
9: 00 10 :0 0 11 :0 0 12 :0 0 13 :0 0 14 :0 0 15 :0 0 16 :0 0 17 :0 0 18 :0 0
8: 00
7: 00
6: 00
0
Pukul (WIB)
Gambar 9 Suhu udara dan RH lingkungan pada tanggal 16 Juni 2007 Dari kondisi iklim mikro lingkungan di atas, analisis distribusi suhu dan RH dalam kandang dilakukan pada tiga waktu berbeda yaitu pada pagi hari (9:20 WIB); siang hari (13:00 WIB) dan sore hari (15:20 WIB). Data input untuk boundary condition pada fluent 6.2 untuk analisis ini dapat dilihat pada Tabel 7 dan Tabel 8.
Tabel 7 Nilai massa jenis, panas jenis dan konduktivitas bahan penyusun kandang Satuan (kg/m3) (kJ/kgoC) (W/moK)
Properties of material Massa jenis (ρ) Panas jenis (Cp) Konduktivitas panas (K)
Concrete 2.310 879 1,2
Asbestos 1922 1,00 4,0
Sumber : Cengel, 2003
Tabel 8 Data input boundary condition untuk fluent 6.2 pada tanggal 16 Juni 2007 Uraian Lingkungan Suhu udara Kecepatan angin Arah angin Atap kanan Tebal Heat fluks Heat transfer coeficient Suhu Atap kiri Tebal Heat fluks Heat transfer coeficient Suhu Tembok kanan Tebal Suhu Free streem temperature Heat transfer coeficient Tembok kiri Tebal Suhu Free streem temperature Heat transfer coeficient Lantai Tebal Suhu Free streem velocity Free streem temperature Heat transfer coeficient Bak air Tebal Suhu Free streem temperature Heat transfer coeficient Depan atas Tebal Heat fluks Heat transfer coeficient Suhu
Satuan 9:20 WIB
Pukul 13:00 WIB
15:20 WIB
C m/det
28,80 1,00 depan
32,52 1,00 kiri
31,80 1,00 depan
m W/m2 W/m2.oC o C
0,005 396,04 5,24 42,00
0,005 506,57 5,38 56,80
0,005 317,32 4,82 48,60
m W/m2 W/m2.oC o C
0,005 396,04 5,14 46,30
0,005 506,57 4,13 41,10
0,005 317,32 3,90 39,80
m C o C W/m2.oC
0,155 27,40 27,40 0,0
0,155 35,00 35,00 2,43
0,155 42,90 42,90 5,92
m C o C W/m2.oC
0,155 29,40 29,40 3,08
0,155 31,90 31,90 0,0
0,155 32,70 32,70 1,22
m C m/det o C W/m2.oC
0,2 28,50 0,00 28,50 1,45
0,2 32,50 0,00 32,50 0,0
0,2 35,90 0,00 35,90 1,22
m C o C W/m2.oC
0,155 29,40 29,40 6,03
0,155 33,90 33,90 3,41
0,155 34,70 34,70 3,73
m W/m2 W/m2.oC o C
0,005 396,04 8,65 34,90
0,005 506,57 5,24 38,20
0,005 317,32 4,68 36,60
o
o
o
o
o
Analisis distribusi suhu dan kelembaban udara dalam kandang dilakukan pada ketinggian 0,6 (posisi sapi berbaring); 1,2 dan 1,6 m (posisi sapi berdiri) dari lantai kandang. Di dalam simulasi menggunakan CFD, material penyusun kandang seperti atap, dinding, lantai, penutup atas dianggap sebagai wall. Inlet merupakan bukaan ventilasi kandang yang tergantung dari arah angin. Pada saat arah angin (inlet) berasal dari depan bangunan kandang, maka outlet-nya adalah bagian bukaan ventilasi yang berada di sebelah kiri, kanan, belakang dan atas (atap). Pada saat angin berasal dari kanan kandang (inlet), outlet berada pada bagian bukaan ventilasi sebelah kiri, depan, belakang dan atas kandang. Suhu udara dalam kandang berasal dari suhu udara lingkungan yang naik pada pagi sampai siang hari dan menurun kembali pada sore hari. Pada pukul 09:20 WIB, suhu udara dalam kandang memiliki kecenderungan meningkat dari posisi dekat lantai menuju posisi dekat atap karena panas matahari yang diterima atap dihantarkan ke dalam kandang sehingga semakin dekat dengan atap suhu udara semakin tinggi. Berbeda dengan kelembaban udara, semakin tinggi suhu udara dalam kandang pada kondisi tekanan uap tetap dan kelembaban mutlak tetap dimana di dalam kandang terjadi proses pemanasan yang dianggap tidak terjadi penambahan uap air,
apabila suhu udara meningkat maka terjadi
penurunan kelembaban udara. Pada atap, suhu udara lebih tinggi karena radiasi matahari yang langsung mengenai atap, dimana suhu dalam kandang masih rendah sehingga panas dari radiasi matahari yang diterima atap dipindahkan secara konveksi ke dalam kandang. Tingginya suhu udara di bagian atap menyebabkan tekanan udara di sekitar atap meningkat dan dengan nilai koefisien tekanan negatif pada bukaan atas, udara terdorong ke luar melalui bukaan atas membawa udara panas dari sekitar atap dan dalam kandang (Gambar 10). Pada pukul 09:20 WIB suhu udara dalam kandang di tiga ketinggian (z=0,6; 1,2 & 1,6 m) lebih rendah dari pada suhu udara lingkungan, karena radiasi matahari yang diterima atap dan konveksi panas dari material penyusun yang dihantarkan masih rendah (Tabel 9). Dari Tabel 9 dapat dilihat bahwa dengan suhu udara lingkungan yang masuk sebesar 28,8oC, suhu udara dalam kandang pada ketinggian 0,6, 1,2 dan 1,6 m paling tinggi sebesar 28,7oC. Sebaliknya, kelembaban udara dalam kandang lebih tinggi dari kelembaban udara lingkungan
karena proses pemanasan dalam kandang masih rendah sehingga uap air dalam kandang belum banyak yang terbuang karena efek panas dan angin lingkungan. Dengan kecepatan angin sebesar 1,0 m/detik dari depan (inlet) kandang pada pukul 09:20 WIB, daerah yang lebih rendah temperaturnya berada di dekat inlet, semakin jauh dengan inlet temperaturnya semakin tinggi. Pada ketinggian 0,6 m udara yang masuk ke kandang terhalang oleh bak penampung air di kanan dan kiri inlet. Dengan tinggi dinding pada bukaan kanan dan kiri sebesar 1,05 m, udara yang masuk tidak dapat keluar melalui outlet sebelah kanan dan kiri sehingga outlet belakang memiliki temperature yang paling tinggi. Pada ketinggian 1,2 dan 1,6 m, udara lingkungan dapat masuk melalui inlet depan sebesar 0,8 (bukaan), outlet kanan dan kiri juga dapat berperan sebagai inlet karena perbedaan suhu udara di dalam dan luar kandang menyebabkan arah gerakan angin sehingga mendorong angin yang berada di sekitar outlet kanan dan kiri masuk ke dalam kandang menuju outlet belakang. Hal ini sejalan dengan yang dikemukakan oleh Takakura (1979), dimana perbedaan tekanan udara dan perbedaan temperatur lingkungan menyebabkan terjadinya pergerakan udara dengan laju yang sangat dipengaruhi oleh kecepatan angin. Masuknya udara lingkungan dalam kandang menyebabkan RH udara dalam kandang juga berubah tergantung dari besarnya suhu udara di dalam dan luar kandang. Secara lebih jelas kontur dan vektor suhu udara dalam kandang pada ketinggian 0,6; 1,2 dan 1,6 m hasil simulasi CFD pada pukul 09:20 WIB dapat dilihat pada Gambar 11 dan Gambar 12 .
Atap Bukaan atas Atap Outlet
Bukaan kanan
Inlet
Dinding kanan (1,05 m) Penampung air
Gambar 10 Sebaran suhu udara dalam kandang pada pukul 09:20 WIB (16 Juni 2007)
Belakang Outlet
Belakang Outlet
Bukaan
Bukaan
Bukaan
Bukaan
kiri
kanan
kiri
kanan
Depan
Depan
Inlet
Z = 0,6 m
Z = 1,2 m
Belakang Outlet
Bukaan
Bukaan
kiri
kanan
Depan
Inlet
Z = 1,6 m Gambar 11 Kontur suhu udara di dalam kandang sapi perah pada pukul 09:20 WIB (16 Juni 2007)
Inlet
Belakang Outlet
Belakang Outlet
Bukaan
Bukaan
Bukaan
Bukaan
kiri
kanan
kiri
kanan
Depan
Inlet
Depan
Z = 0,6 m
Z = 1,2 m
Belakang Outlet
Bukaan
Bukaan
kiri
kanan
Depan
Inlet
Z = 1,6 m Gambar 12 Vektor kecepatan aliran suhu udara di dalam kandang sapi perah pada pukul 09:20 WIB (16 Juni 2007)
Inlet
Tabel 9 Suhu udara dalam kandang sapi perah FH hasil analisis CFD Ketinggian z (m) 0,6
1,2
1,6
Suhu udara (oC)
Nilai
9:20 WIB
13:00 WIB
15:20 WIB
Minimum
28,65
32,37
32,20
Maksimum
28,70
32,65
32,44
Rata-rata
28,69
32,57
32,37
Coefficient of variance (%)
0,0244
0,0649
0,1382
Minimum
28,68
32,52
32,33
Maksimum
28,70
32,65
32,40
Rata-rata
28,69
32,61
32,38
Coefficient of variance (%)
0,0105
0,0000
0,0675
Minimum Maksimum
28,68 28,70
32,55 32,65
32,36 32,39
Rata-rata
28,69
32,63
32,38
Coefficient of variance (%)
0,0070
0,0124
0,0460
Pada siang hari (pukul 13:00 WIB) dan sore hari (pukul 15:20 WIB) ketika komponen penyusun kandang (atap, lantai, dinding dan rangka) telah menyimpan dan menghantarkan panas, suhu udara dalam kandang lebih tinggi dari suhu lingkungan sehingga kelembaban udara (RH) dalam kandang menurun pada pukul 13:00 WIB dan meningkat kembali pada sore hari.
Naik dan turunnya
kelembaban udara dalam kandang dipengaruhi langsung oleh suhu udara lingkungan. Nilai kelembaban udara dalam kandang pada pukul 09:20; 13:00 dan 15:20 WIB (16 Juni 2007) di ketinggian 0,6, 1,2 dan 1,6 m berkisar antara 64,1380,90% (Tabel 10) Tabel 10 RH udara dalam kandang sapi perah FH hasil analisis CFD Ketinggian z (m) 0,6
1,2
1,6
Nilai
RH udara (%) 9:20 WIB
13:00 WIB
15:20 WIB
Minimum
80,70
64,13
68,00
Maksimum
80,90
65,00
68,80
Rata-rata
80,74
64,75
68,23
Minimum
80,70
64,59
68,14
Maksimum
80,78
65,00
68,37
Rata-rata
80,74
64,87
68,20
Minimum
80,70
64,69
68,17
Maksimum
80,78
65,00
68,27
Rata-rata
80,74
64,94
68,20
Pada pukul 13:00 WIB, dengan kecepatan angin 1,00 m/detik dan arah (inlet) dari kiri kandang, suhu udara dalam kandang terdistribusi hampir merata di bagian inlet dan outlet. Suhu udara terendah berada di dekat dinding kanan (bawah outlet) pada katinggian kurang dari 1,05 m. Rendahnya suhu udara di bawah dinding kanan disebabkan udara lingkungan yang masuk melalui inlet (bukaan kiri) langsung menuju outlet (bukaan sebelah kanan) dan akibat terhalangi tembok (1,05 m) yang berada di bukaan kanan, udara dibelokkan kembali ke tengah, arah inlet (bukaan kiri) dan atas untuk diteruskan ke outlet. Hal ini merupakan sifat dari udara yang akan membelokkan pola alirannya apabila mengenai suatu halangan yang tidak dapat dilewatinya. Pada pukul 13:00 WIB, dimana suhu udara dalam kandang lebih tinggi dari suhu udara lingkungan, dengan heat transfer coeficient (h) pada atap kiri yang lebih rendah (0,99 W/m2.oC) dari pukul 09:20 WIB, suhu di dekat (bawah atap) relatif lebih rendah dari suhu di inlet dan outlet. Kondisi tersebut menunjukkan bahwa pada siang hari (pukul 13:00 WIB) dimana suhu kandang lebih tinggi dari suhu udara lingkungan, suhu atap tidak berpengaruh secara signifikan terhadap suhu udara di bawahnya karena berfungsinya bukaan (outlet) atas sehingga proses pemanasan di bawah atap tereduksi oleh udara yang keluar melalui bukaan atas. Rendahnya proses pemanasan dalam kandang akibat panas yang dipindahkan secara konveksi oleh atap (pukul 13:00 WIB) menyebabkan nilai RH yang terdistribusi dalam kandang (64,85%) mendekati RH lingkungan (65,25). Sebaran udara dalam kandang pada pukul 13:00 WIB (16 Juni 2007) disajikan pada Gambar 13.
Atap
Bukaan atas Atap Inlet
Outlet
Dinding kanan (1,05 m) Penampung air
Gambar 13 Sebaran suhu udara dalam kandang pada pukul 13:00 WIB (16 Juni 2007)
Pada pukul 13:00 WIB, di ketinggian (z=0,6 m) dengan kecepatan angin sebesar 1,0 m/detik dari arah kiri kandang (inlet), bukaan outlet di depan dan belakang kandang berfungsi dengan baik, sedangkan pada bukaan sebelah kanan karena terhalangi tembok (1,05 m) udara berubah arah ke bukaan depan dan belakang serta berbalik ke arah tembok kiri. Kondisi ini menyebabkan suhu udara tertinggi berada pada daerah dekat tembok kiri (inlet), sedangkan RH tertinggi berada pada daerah tembok kanan. Tingginya suhu udara di daerah dekat inlet dipengaruhi oleh panas yang dipancarkan atap dan tembok sebalah kanan yang dibawa oleh gerakan angin yang berputar menuju inlet. Kondisi tersebut membuktikan bahwa luas bukaan ventilasi menjadi faktor yang cukup penting pada perhitungan distribusi udara dan RH di dalam kandang selain faktor kecepatan angin dan tekanan udara (Takakura, 1979). Pada z=1,2 dan 1,6 m, bukaan ventilasi yang berperan sebagai outlet adalah bukaan bagian kanan, karena udara dapat bergerak tanpa halangan (tembok). Pada bukaan ventilasi sebelah kiri, tekanan udara di luar kandang lebih tinggi dari tekanan udara di dalam kandang sehingga terjadi aliran udara ke dalam kandang dan keluar melalui inlet di bukaan sebalah kanan. Daerah bukaan sebelah kanan berfungsi sebagai outlet karena tekanan udara di dalam kandang lebih tinggi dari tekanan udara di luar kandang (sebelah kanan), bukaan ventilasi sebelah kanan merupakan daerah di atas bidang tekanan netral. Sementara pada bukaan ventilasi sebelah depan dan belakang memiliki tekanan udara yang sama dengan tekanan udara luar kandang sehingga dapat berfungsi sebagai bidang tekanan netral (Brockett & Albright, 1987). Hal ini terbukti dengan tidak keluarnya aliran udara pada bukaan sebelah depan dan belakang kandang. Pada kondisi seperti ini distribusi suhu udara dalam kandang lebih dominan dipengaruhi oleh efek termal daripada efek angin. Secara lebih jelas kontur dan vektor aliran udara dalam kadang hasil simulasi dengan CFD pada pukul 13:00 WIB di ketinggian 0,6; 1,2 dan 1,6 m dapat dilihat pada Gambar 14 dan Gambar 15.
Bukaan belakang
Bukaan belakang
Kiri
Kanan
Kiri
Kanan
Inlet
Outlet
Inlet
Outlet
Bukaan depan
Bukaan depan
Z = 0,6 m
Z = 1,2 m
Bukaan belakang
Kiri
Kanan
Inlet
Outlet
Bukaan depan
Gambar 14
Z = 1,6 m Kontur suhu udara di dalam kandang sapi perah pada pukul 13:00 WIB (16 Juni 2007)
Bukaan belakang
Bukaan belakang
Kiri
Kanan
Kiri
Kanan
Inlet
Outlet
Inlet
Outlet
Bukaan depan
Bukaan depan
Z = 0,6 m
Z = 1,2 m
Bukaan belakang
Kiri
Kanan
Inlet
Outlet
Bukaan depan
Z = 1,6 m Gambar 15 Vektor kecepatan aliran suhu udara di dalam kandang sapi perah pada pukul 13:00 WIB (16 Juni 2007)
Pada pukul 15:20 WIB, dengan kecepatan angin 1,00 m/detik dan arah (inlet) dari depan kandang, suhu udara dalam kandang terdistribusi merata di bagian inlet dan outlet. Suhu udara terendah berada di dekat atap (penutup) sebelah depan karena tidak terkena radiasi matahari. Pada pukul 15:20 WIB, suhu material penyusun kandang (atap, tembok, lantai, penampung air) lebih tinggi dari suhu udara lingkungan. Tingginya suhu material penyusun atap menyebabkan kandang menjadi panas akibat panas yang dikonveksikan oleh material penyusun bahan kandang ke dalam kandang sehingga suhu di dalam kandang tersebar secara merata. Tingginya suhu di dalam kandang menyebakan tekanan dalam kandang meningkat sehingga udara terdorong ke luar kandang melalui outlet yang tersebar di bukaan belakang, kanan dan kiri (Bockett and Albright, 1987). Meningkatnya suhu udara dalam kandang akibat panas yang dikonveksikan material penyusun kandang menyebakan RH dalam kandang (68,21%) lebih rendah dari RH lingkungan (68,85%). Pada bagian atap, suhunya hampir sama dengan suhu pada semua bagian di dalam kandang yang menunjukkan bahwa proses pemanasan yang relatif kecil dalam kandang terjadi secara merata pada daerah sekitar material penyusun kandang (lantai, dinding kanan/kiri, penampung air dan atap). Sebaran udara dalam kandang pada pukul 15:20 WIB (16 Juni 2007) disajikan pada Gambar 16
Atap Bukaan atas Atap
Bukaan kanan Outlet
Dinding kanan (1,05 m) Inlet
Penampung air
Gambar 16 Sebaran suhu udara dalam kandang pada pukul 15:20 WIB (16 Juni 2007) Pada pukul 15:20 WIB, dengan kecepatan angin sebesar 1,0 m/detik dan arah angin (inlet) dari depan kandang, bukaan outlet di belakang kandang berfungsi dengan baik (pada z=0,6 m), sedangkan pada bukaan sebelah kanan dan kiri udara yang dibawa angin terhalangi tembok (1,05 m). Suhu udara terendah berada di sebelah kanan
dan kiri bukaan, tertinggi berada di tengah kandang. Pada z=1,2 dan 1,6 m, bukaan ventilasi yang berperan sebagai outlet adalah bukaan bagian kanan, kiri dan belakang. Pada bukaan ventilasi sebelah depan, tekanan udara di luar kandang lebih tinggi dari tekanan udara di dalam kandang sehingga terjadi aliran udara ke dalam kandang dan keluar melalui inlet di bukaan sebalah kanan, kiri dan belakang. Distribusi suhu udara pada z=1,2 dan z=1,6 m tersebar merata di seluruh bidang pada kandang yang sangat dipengaruhi oleh efek termal yang ditimbulkan oleh radiasi matahari dan material bahan penyusun kandang yang mengeluarkan panas. Kondisi ini dapat dilihat dari berperannya bukaan ventilasi sebelah kanan, kiri dan belakang sebagai outlet yang menunjukkan bahwa tekanan udara dalam kandang (pada arah kanan, kiri dan belakang) lebih tinggi dari tekanan udara luar kandang. Secara lebih jelas kontur dan vektor aliran udara dalam kadang hasil simulasi dengan CFD pada pukul 15:20 WIB di ketinggian 0,6; 1,2 dan 1,6 m dapat dilihat pada Gambar 17 dan Gambar 18. Dari hasil analisis di atas, distribusi suhu dan kelembaban udara dalam kandang sapi perah FH (tanggal 16 Juni 2007) pada pukul 09:20 dengan suhu dan kelembaban udara maksimum pada ketinggian 0,6, 1,2 dan 1,6 m sebesar 28,7oC dan 90%, pada pukul 13:00 dengan suhu
dan kelembaban udara pada ketiga ketinggian sebesar
o
32,65 C dan 65% menyebabkan terjadi stress sedang pada sapi perah. Demikian juga dengan suhu dan kelembaban udara untuk ketiga ketinggian pada pukul 15:20 (masingmasing sebesar 32,44oC dan 68,8%) menyebabkan sapi perah dalam kondisi stress sedang (Wierema, 1990). Stres tersebut akan menurunkan produktivitas sapi perah FH yang diindikasikan dengan: 1) penurunan nafsu makan; 2) peningkatan konsumsi minum; 3) penurunan metabolisme dan peningkatan katabolisme; 4) peningkatan pelepasan panas melalui penguapan; 5) penurunan konsentrasi hormon; 6) peningkatan temperatur tubuh, respirasi dan denyut jantung (McDowell, 1972); 7) perubahan tingkah laku (Ingram & Dauncey, 1985); 8) meningkatkan intensitas berteduh sapi (Combs, 1996). Untuk mengurangi tingkat stres pada sapi perah FH dapat dilakukan melalui modifikasi disain kandang dengan cara merubah tinggi dan lebar kandang dan memperluas bukaan ventilasi kandang agar suhu dalam kandang lebih rendah.
Belakang Outlet
Belakang Outlet
Bukaan
Bukaan
Bukaan
Bukaan
kiri
kanan
kiri
kanan
Depan
Inlet
Depan
Z = 0,6 m
Z = 1,2 m
Belakang Outlet
Bukaan
Bukaan
kiri
kanan
Depan
Inlet
Z = 1,6 m Gambar 17 Kontur suhu udara di dalam kandang sapi perah pada pukul 15:20 WIB (16 Juni 2007)
Inlet
Belakang Outlet
Belakang Outlet
Bukaan
Bukaan
Bukaan
Bukaan
kiri
kanan
kiri
kanan
Depan
Inlet
Depan
Z = 0,6 m
Z = 1,2 m
Belakang Outlet
Bukaan
Bukaan
kiri
kanan
Depan
Inlet
Z = 1,6 m Gambar 18 Vektor kecepatan aliran suhu udara di dalam kandang sapi perah pada pukul 15:20 (16 Juni 2007)
Inlet
Validasi Suhu dan Kelembaban Udara pada Kandang Sapi Perah Validasi distribusi suhu udara dilakukan dengan cara membandingkan data suhu udara hasil pengukuran dengan data suhu udara hasil simulasi menggunakan CFD di 24 titik dalam kandang. Validasi dilakukan pada kondisi kandang tidak diisi sapi (kandang kosong) sebanyak 3 kali pada tanggal 16 Juni 2007 yaitu pada pagi hari (pukul 9:20 WIB), siang hari (pukul 13:00 WIB) dan sore hari (pukul 15:20 WIB). Hasil validasi distribusi suhu udara dalam kandang sapi perah menunjukkan kecenderungan hasil simulasi CFD mendekati hasil pengukuran dengan nilai standar deviasi dan error yang rendah (Lampiran 2). Nilai minimum, maksimum dan rata-rata standar deviasi hasil validasi pada pukul 09:20; 13:00 dan 15:20 WIB (16 Juni 2007) dapat dilihat pada Tabel 11. Nilai rata-rata error pada pukul 09:20; 13:00 dan 15:20 WIB (16 Juni 2007) masing-masing sebesar 1,90; 1,40 dan 1,28%. Pada beberapa titik terjadi perbedaan yang cukup mencolok karena terkait dengan penentuan jarak grid yang sedikit berbeda antara pengukuran dan simulasi, tetapi masih dalam batasan yang rendah (standar deviasi dan error < 5oC). Nilai standar deviasi dan error yang rendah tersebut menunjukkan bahwa simulasi menggunakan CFD memiliki akurasi yang cukup tinggi sehingga dapat dijadikan acuan untuk perancangan kandang sapi perah FH dalam perspektif distribusi suhu. Selanjutnya, data input dalam solver untuk keperluan simulasi desain kandang diambil dari data pengukuran pada siang hari (pukul 13:00 WIB), karena pada siang hari radiasi matahari mencapai puncaknya, demikian juga dengan suhu udara dalam kandang. Validasi suhu udara hasil pengukuran dan hasil simulasi CFD di 24 titik pada pukul 09:20, 13:00 dan 15:20 WIB (16 Juni 2007) dapat dilihat pada Gambar 19 sampai dengan Gambar 21. Tabel 11 Hasil validasi suhu udara pengukuran dengan suhu udara hasil CFD dalam kandang Nilai standar deviasi (oC)
Nilai
9:20 WIB
13:00 WIB
15:20 WIB
Minimum
0,00
0,02
0,03
Maksimum
1,13
0,79
0,74
Rata-rata
0,39
0,33
0,30
35.00
SD = 0,39 0C;
Error = 1,90%
0
Suhu ( C)
33.00 31.00 29.00 27.00 25.00 1
2
3
4
5
6
7
8
9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
Titik pengukuran T Simulas i
T Ukur
Gambar 19 Validasi suhu hasil simulasi CFD terhadap suhu pengukuran di 24 titik dalam kandang (pukul 09:20 WIB, tanggal 16 Juni 2007)
35.00
SD = 0,330C;
Error = 1,40%
0
Suhu ( C)
33.00 31.00 29.00 27.00 25.00 1
2
3
4
5
6
7
8
9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
Titik pengukuran T Simulasi
T Ukur
Gambar 20 Validasi suhu hasil simulasi CFD terhadap suhu pengukuran di 24 titik dalam kandang (pukul 13:00 WIB, tanggal 16 Juni 2007)
35.00
SD = 0,300C;
Error = 1,28%
0
Suhu ( C)
33.00 31.00 29.00 27.00 25.00 1
2
3
4
5
6
7
8
9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
Titik pengukuran T Simulas i
T Ukur
Gambar 21 Validasi suhu hasil simulasi CFD terhadap suhu pengukuran di 24 titik dalam kandang (pukul 15:20 WIB, tanggal 16 Juni 2007) Validasi kelembaban udara (RH) dalam kandang dilakukan dengan membandingkan RH ukur dengan RH hitung (didasarkan pada suhu hasil simulasi menggunakan CFD) di 3 titik pada tiga waktu (pukul 09:20, 13:00 dan 15:20 WIB, tanggal 16 Juni 2007). Validasi RH dilakukan pada saat kandang tidak diisi sapi (kandang kosong). Nilai RH ukur dan RH hitung dan validasinya dapat dilihat pada Lampiran 3. Secara umum terdapat kecenderungan yang sama antara RH ukur dengan RH hasil perhitungan menggunakan CFD. Perbedaan secara umum dinyatakan dalam standar deviasi sebesar 2,44% dengan standar deviasi rata-rata pada pukul 09:20 WIB sebesar 3,85% (2,61 – 5,55%), pada pukul 13:00 WIB sebesar 2,12% (0,06 – 3,53%) dan pada pukul 15:20 WIB sebesar 1,37% (0,13 – 2,96%). Nilai error rata-rata antara RH hasil pengukuran dengan RH simulasi CFD pada pukul 09:20; 13:00 dan 15:20 WIB adalah sebesar 4,73%. Pada beberapa titik terdapat perbedaan yang mencolok akibat penentuan jarak grid yang sedikit berbeda antara pengukuran dengan simulasi, tetapi secara umum masih dalam batasan standar deviasi dan error yang rendah (standar deviasi dan error < 5%). Rendahnya nilai standar deviasi menunjukkan bahwa validasi RH memiliki akurasi yang tinggi sehingga dapat dijadikan acuan untuk perancangan kandang sapi perah FH dalam perspektif distribusi RH pada saat kandang tidak diisi sapi. Validasi RH ukur dengan RH hitung dapat dilihat pada Gambar 22.
90
SD = 2,44%;
80
Error = 4,73%
70
RH(%)
60 50 40 30 20 10 0 1
2
3
4
5
6
7
8
9
10
11
12
Titik pengukuran RH CFD
RH ukur
Gambar 22 Validasi RH hasil simulasi CFD terhadap RH pengukuran di 4 titik dalam kandang pada pukul 09:20 (titik 1-4), 13:00 (titik 5-8) dan 15:20 WIB (titik 9-12) pada tanggal 16 Juni 2007 Simulasi Disain Kandang Sapi Perah Simulasi dilakukan dengan melibatkan 20 ekor sapi perah yang ditempatkan dalam kandang. Rata-rata bobot badan sapi perah adalah 350 kg dengan rataan luas kulit sebesar 3,47 m2 (Lampiran 4) yang diletakkan secara merata di dalam kandang (Lampiran 5). Peletakan kulit sapi (radiator) dalam simulasi menggunakan CFD dimodelkan dengan hamparan kulit berbentuk persegi panjang pada arah x (193 cm) dan arah y (900 cm) pada dua ketinggian (z) dengan jarak 20 cm dari tembok kiri dan kanan, 62 dan 125 cm dari lantai (Lampiran 6). Nilai koefisien pindah panas konveksi pada kulit sapi perah tergantung dari feed intake (Lampiran 7).
Kondisi awal kandang sebelum
dilakukan simulasi berupa kandang sapi perah FH dengan tinggi 5,75 m, lebar 6,3 m, tinggi dinding kanan dan kiri 1,05 m, sebelah depan dan belakang terdapat bak penampung air dengan tinggi 1,05 m akan disimulasikan dengan diisi sapi (Gambar 23) dengan kecepatan angin 0,7 m dari arah kiri/kanan dan depan/belakang. Data distribusi suhu yang diperoleh akan dibandingkan dengan hasil simulasi (merubah tinggi, lebar dan tinggi dinding kandang serta penempatan bak penampung air). Bentuk geometri kandang simulasi dapat dilihat pada Gambar 24. Data input untuk fluent 6.2 pada simulasi dapat dilihat pada Tabel 12.
Outlet Inlet
Outlet Belakang Kiri
Tembok kiri
1
Kanan
1
1
1
Tembok kanan
Depan Inlet
Penampung air Kulit sapi
Gambar 23 Peletakan kulit sapi perah pada geometri kandang awal
Tembok kiri (0,4 m) Tembok kanan (0,4 m)
Gambar 24 Peletakan kulit sapi perah pada geometri kandang simulasi Simulasi dilakukan pada 9 disain kandang dengan dimensi: tinggi kandang (T1=5,25 m; T2=5,75 m & T3=6,25 m), lebar kandang (L1=6,3 m; L2=7,3 m & L3=8,3 m), tinggi dinding (0,4 m) dan posisi bak penampung air dipindah dari letak awal. Arah angin berasal dari depan/belakang kandang dan kanan/kiri kandang dengan kecepatan yang sama yaitu sebesar 0,7 m/detik.
Tabel 12 Data input boundary condition untuk fluent 6.2 untuk simulasi Uraian Lingkungan Suhu udara Kecepatan angin Arah angin Atap kanan Tebal Heat fluks Heat transfer coeficient Suhu Atap kiri Tebal Heat fluks Heat transfer coeficient Suhu Tembok kanan Tebal Suhu Free streem velocity Free streem temperature Heat transfer coeficient Tembok kiri Tebal Suhu Free streem velocity Free streem temperature Heat transfer coeficient Lantai Tebal Suhu Free streem velocity Free streem temperature Heat transfer coeficient Bak air Tebal Suhu Free streem velocity Free streem temperature Heat transfer coeficient Depan atas Tebal Heat fluks Heat transfer coeficient Suhu Kulit sapi Heat transfer coeficient Suhu
Satuan o
C m/det
Nilai 32,25 0,70 depan/belakang dan kanan/kiri
m W/m2 W/m2.oC o C
0,005 597,2 1,16 55,1
m W/m2 W/m2.oC o C
0,005 597,2 3,23 35,9
m 0 C m/det o C W/m2.oC
0,155 32,7 0,00 32,7 3,68
m o C m/det o C W/m2.oC
0,155 27,3 0,00 27,3 0,0
m o C m/det o C W/m2.oC
0,2 26,9 0,00 26,9 0,0
m o C m/det o C W/m2.oC
0,155 31,1 0,00 31,1 0,0
M W/m2 W/m2.oC o C
0,005 597,2 2,31 36,9
W/m2.oC o C
55,41 36,21
Hasil Simulasi Kondisi Awal Arah Angin (Inlet) dari Bukaan Kanan Hasil simulasi distribusi suhu udara dalam kandang pada kondisi awal dengan arah angin (inlet) dari kanan/kiri menunjukkan bahwa suhu udara dalam kandang terdistribusi secara tidak merata pada tiap ketinggian (Gambar 25). Suhu tertinggi berada di dekat kulit sapi (radiator) yang tersebar dari arah inlet sampai outlet. Pada arah outlet (bukaan kiri) sebaran suhu udara tinggi lebih luas dari inlet (bukaan kanan) karena berfungsinya bukaan kiri sebagai outlet yang membuang panas dalam kandang. Tingginya suhu udara dalam kandang di sekitar radiator menunjukkan bahwa meningkatnya suhu dalam kandang didominasi oleh pengaruh kulit sapi yang memancarkan panas secara terus menerus dengan suhu kulit 36,21oC dan heat transfer coefficient yang cukup tinggi (55,41 W/m2.oC). Panas yang terkonveksi dari material penyusun kandang (atap, dinding kanan dan kiri) tidak signifikan
meningkatkan suhu dalam kandang, selain disebabkan
luasnya bukaan ventilasi, nilai heat transfer coefficient dari material penyusun kandang juga lebih rendah dibandingkan dengan kulit sapi perah FH sebagai radiator (Tabel 12). Tingginya suhu udara dalam kandang menyebabkan tekanan udara dalam kandang meningkat sehingga lebih tinggi dari tekanan udara di luar kandang menyebabkan udara panas dalam kandang keluar melalui bukaan kiri sebagai outlet (Bockett & Albright, 1987).
Atap Bukaan atas Atap
Outlet Inlet
Bukaan depan
Penampung air
Dinding kanan (1,05 m) Bukaan kanan
Gambar 25 Sebaran suhu udara dalam kandang hasil simulasi pada kondisi awal (arah angin (inlet) dari kanan) Suhu udara pada ketinggian 0,6 m lebih tinggi dari suhu udara pada ketinggian 1,2 dan 1,6 m (Tabel 13). Kondisi ini menunjukkan bahwa kulit sapi sebagai radiator yang ditempatkan pada ketinggian 0,62 dan 1,25 m dengan suhu
sebesar 36,21oC cukup dominan dalam memberikan panas ke kandang walaupun suhu atap cukup tinggi (atap kanan 55,1oC dan atap kiri 35,9oC). Pada ketinggian 0,6 m, suhu udara tertinggi berada di tengah kandang, sedangkan suhu udara terendah berada di dekat inlet (bukaan sebalah kanan). Tingginya suhu udara di tengah kandang disebabkan oleh panas yang dihasilkan radiator menuju ke tengah kandang yang selanjutnya menuju bukaan di atas radiator (bukaan kanan, depan, belakang, atas) sebagai outlet (tekanan udara dalam kandang lebih besar dari tekanan udara di luar kandang). Arah distribusi udara menuju ke tengah dan bukaan atas disebabkan oleh kurangnya bukaan ventilasi di bagian kanan/kiri dan depan/belakang akibat dinding setinggi 1,05 m. Martalerz (1977) mengemukakan bahwa laju pertukaran udara dipengaruhi oleh total luas dan arah bukaan, kecepatan angin dan perbedaan temperatur di dalam dan di luar kandang. Secara lebih jelas kontur suhu udara dalam kandang hasil simulasi kondisi awal pada arah angin (inlet) dari bukaan kanan dapat dilihat pada Gambar 26. Tabel 13 Suhu udara (oC) dalam kandang sapi perah FH hasil analisis CFD pada kondisi awal dengan inlet dari kanan/kiri dan depan/belakang z (m)
Inlet Nilai
0,6
1.2
1.6
Kanan/kiri
Depan/belakang
Minimum
34,030
32,250
Maksimum
34,520
35,618
Rata-rata
34,240
34,035
Coefficient of variance (%)
0,339
2,942
Minimum
33,380
32,250
Maksimum
34,420
34,988
Rata-rata
33,970
33,069
Coefficient of variance (%)
1,048
1,925
Minimum
32,830
32,250
Maksimum
34,030
35,596
Rata-rata
33,710
33,786
Coefficient of variance (%)
0,953
2,682
Bukaan belakang
Bukaan belakang
Kiri
Kanan
Kiri
Kanan
Outlet
Inlet
Outlet
Inlet
Bukaan depan
Bukaan depan
Z = 0,6 m
Z = 1,2 m
Bukaan belakang
Kiri
Kanan
Outlet
Inlet
Bukaan depan
Z = 1,6 m Gambar 26 Kontur suhu udara di dalam kandang sapi perah pada kondisi awal simulasi (arah angin (inlet) dari kanan)
Pada ketinggian 1,2 m (diantara radiator 1 dan 2), suhu udara terendah berada di dekat inlet dan suhu udara tertinggi berada di dekat outlet. Udara lingkungan yang masuk ke kandang melalui bukaan kiri membawa panas yang dihasilkan kulit sapi dan pancaran panas dari material penyusun bahan ke outlet. Pada bagian kanan kandang, dimana tekanan udara lebih rendah dari tekanan udara dalam kandang berfungsi sebagai outlet. Sementara pada bukaan bagian depan dan belakang tidak terdapat perbedaan yang signifikan antara tekanan udara di dalam dan luar kandang sehingga tidak berfungsi sebagai outlet (Brockett & Albright, 1987). Pada ketinggian 1,6 m (di atas radiator), suhu udara dalam kandang memiliki nilai rata-rata terendah dibandingkan dengan ketinggian 0,6 dan 1,2 m. Rendahnya distribusi udara pada ketinggian ini diakibatkan udara panas yang dihembuskan kulit sapi tereduksi oleh udara lingkungan yang masuk dan keluar melalui bukaan inlet dan outlet. Besar kecilnya bukaan ventilasi berpengaruh terhadap laju aliran udara yang keluar dan masuk ke kandang (Martalerz, 1977). Rata-rata suhu udara dalam kandang pada ketinggian 1,6 (33,71oC) lebih tinggi dari suhu udara lingkungan yang masuk (32,25oC) sehingga sistem ventilasi belum berjalan dengan baik pada ketinggian 1,6 m. Gardjito (2002) menyatakan bahwa sistem ventilasi alamiah yang baik adalah sistem yang sanggup menurunkan suhu di dalam ruangan sampai sama dengan suhu udara luar yang sangat bergantung pada faktor iklim setempat dan faktor rancangan bangunan dengan sistem ventilasinya. Arah Angin (Inlet) dari Bukaan Depan Suhu udara dalam kandang simulasi (kondisi awal) dengan arah angin (inlet) dari depan terdistribusi secara tidak merata pada tiap sisi dan ketinggian kandang (Gambar 27). Suhu tertinggi berada di dekat sapi (radiator) karena sapi mengeluarkan panas dengan suhu kulit sebesar 36,21oC dan heat transfer coefficient yang cukup besar (55,41 W/m2.oC). Selain itu, adanya dinding penampung air di bukaan depan dan dinding bukaan kanan dan kiri menyebabkan udara panas dalam kandang terjebak karena udara tidak dapat melewati dinding dan berbelok kea rah tengah dan bukaan kanan. Panas yang terkonveksi dari material penyusun kandang (atap, dinding kanan dan kiri) tidak signifikan meningkatkan suhu dalam kandang, selain disebabkan luasnya bukaan ventilasi, nilai heat transfer coefficient dari material penyusun kandang juga
lebih rendah dibandingkan dengan kulit sapi perah FH sebagai radiator (Tabel 12). Meningkatnya suhu udara dalam kandang menyebabkan tekanan udara dalam kandang meningkat menyebabkan udara panas dalam kandang keluar melalui bukaan kanan dan bukaan belakang sebagai outlet.
Atap Bukaan atas Atap Outlet
Dinding kanan (1,05 m) Bukaan kanan Inlet
Penampung air
Gambar 27 Sebaran suhu udara dalam kandang hasil simulasi pada kondisi awal (arah angin (inlet) dari depan) Suhu udara pada ketinggian 0,6 m lebih tinggi dari suhu udara pada ketinggian 1,2 dan 1,6 m (Tabel 13) yang menunjukkan bahwa kulit sapi sebagai radiator (0,62 dan 1,25 m) cukup dominan memberikan panas ke kandang. Pada ketinggian 0,6; 1,2 dan 1,6 m, suhu udara tertinggi berada di sekitar sapi sampai ke arah dinding kanan dan dinding kiri, sedangkan suhu udara terendah berada di tengah kandang mulai dari inlet (bukaan depan) sampai outlet (bukaan belakang). Rendahnya suhu udara di tengah kandang disebabkan oleh panas yang dihasilkan kulit sapi (radiator) yang menuju ke tengah kandang dapat segera dibuang oleh udara lingkungan yang masuk melalui outlet (bukaan belakang). Arah distribusi udara menuju ke tengah disebabkan oleh kurangnya bukaan ventilasi di bagian bukaan kanan dan kiri akibat terhalangi oleh sapi (radiator). Distribusi suhu udara pada kandang simulasi (kondisi awal) dengan arah angin (inlet) dari depan mengahasilkan nilai yang lebih rendah daripada pada saat inlet berasal dari bukaan kanan dengan selisih rata-rata di tiga ketinggian (z=0,6; 1,2 dan 1,6 m) sebesar 0,17oC. Secara lebih jelas kontur suhu udara dalam kandang hasil simulasi kondisi awal pada arah angin (inlet) dari bukaan depan dapat dilihat pada Gambar 28.
Belakang Outlet
Belakang Outlet
Bukaan
Bukaan
Bukaan
Bukaan
kiri
kanan
kiri
kanan
Depan
Depan
Inlet
Z = 0,6 m
Inlet
Z = 1,2 m
Belakang Outlet
Bukaan
Bukaan
kiri
kanan
Depan
Inlet
Z = 1,6 m Gambar 28 Kontur suhu udara di dalam kandang sapi perah pada kondisi awal simulasi (arah angin (inlet) dari depan)
Hasil Simulasi dengan Arah Angin dari Depan/belakang dan Kanan/kiri Simulasi dilakukan pada 9 disain kandang dengan luas kulit sapi sebagai radiator yang sama. Disain kandang simulasi memiliki dimensi: tinggi kandang (T1=5,25 m; T2=5,75 m & T3=6,25 m), lebar kandang (L1=6,3 m; L2=7,3 m & L3=8,3 m), tinggi dinding (0,4 m) dan posisi bak penampung air dipindah dari letak awal. Arah angin berasal dari depan/belakang kandang dan kanan/kiri kandang dengan kecepatan yang sama yaitu sebesar 0,7 m/det. Hasil simulasi disain kandang pada ketinggian (z=0,6; 1,2 & 1,6 m) dengan arah angin dari depan/belakang disajikan pada Tabel 14, sedangkan hasil simulasi dengan arah angin dari kanan/kiri disajikan pada Tabel 15. Dari Tabel 14 dan 15 dapat dilihat bahwa distribusi suhu udara dalam kandang hasil simulasi pada ketinggian 0,6; 1,2 dan 1,6 m sangat dipengaruhi oleh luas bukaan ventilasi. Pada tinggi atap kandang 5,25 m, distrubusi suhu udara dalam kandang akan menurun dengan bertambah lebar bukaan inlet dan outlet dari depan/belakang kandang. Suhu udara tertinggi berada pada ketinggian 1,2 m karena berada pada dua radiator (kulit sapi) yang memancarkan panas sebesar 36,21oC. Kondisi yang sama terjadi pada tinggi atap kandang 5,75 m dan 6,25 m. Semakin tinggi atap kandang maka bukaan ventilasi juga semakin luas sehingga suhu udara dalam kandang akan menurun. Bukaan ventilasi yang semakin besar menyebabkan pertukaran udara di dalam dan luar kandang semakin tinggi sehingga suhu udara dalam kandang akan lebih cepat turun sebanding dengan bertambahnya bukaan ventilasi (Mastalerz, 1977). Bukaan ventilasi pada simulasi diperluas dengan cara menurunkan dinding kanan dan kiri kandang dari 1,05 m menjadi 0,4 m, memindahkan tempat penampung air yang berada di depan dan belakang kandang (posisi awal setinggi 1,05 m menjadi 0 m). Udara lingkungan masuk ke kandang akibat perbedaan temperatur antara di luar (32,25oC) dan di dalam kandang (> 33oC) sehingga tekanan udara di dalam dan di luar kandang juga berbeda yang menyebabkan terjadinya aliran udara masuk ke kandang melalui bukaan inlet searah dengan arah angin. Laju masuknya udara lingkungan ke dalam kandang dipengaruhi oleh luas bukaan ventilasi, kecepatan angin, arah bukaan, perbedaan temperatur di dalam dan luar kandang.
Tabel 14 Distribusi suhu udara hasil simulasi CFD pada beberapa disain kandang dengan arah angin dari depan/belakang Z (m) 0.6
Nilai
Minimum (0C) Maksimum (0C) Rata-rata (0C) Standard deviasi Coefficient of variance (%) 1.2 Minimum (0C) Maksimum (0C) 0 Rata-rata ( C) Standard deviasi Coefficient of variance (%) 1.6 Minimum (0C) 0 Maksimum ( C) 0 Rata-rata ( C) Standard deviasi Coefficient of variance (%) Rata-rata pada z=0.6, 1.2 dan 1.6 m
Tinggi Atap 5.25 m L1 L2 L3 32.250 32.243 32.244 34.212 34.426 34.775 33.017 32.897 32.835 0.565 0.586 0.613 1.711 1.781 1.868 32.250 32.250 32.250 35.601 35.600 35.629 33.756 33.720 33.378 1.123 1.124 1.177 3.326 3.334 3.525 32.250 32.250 32.250 35.698 35.634 35.593 33.688 33.455 33.232 1.119 1.130 1.097 3.321 3.378 3.302 33.487 33.357 33.148
Tinggi Atap 5.75 m L1 L2 32.231 32.247 34.051 34.426 32.944 32.890 0.641 0.590 1.956 1.794 32.250 32.250 35.323 35.621 33.767 33.261 1.154 1.140 3.469 3.376 32.250 32.250 35.433 35.634 33.350 33.455 1.193 1.130 3.578 3.378 33.354 33.202
L3 32.174 34.159 32.776 0.731 2.220 32.244 35.281 33.448 1.175 3.512 32.250 35.015 32.820 0.917 2.793 33.015
Tinggi Atap 6.25 m L1 L2 32.245 32.250 34.143 34.212 32.842 33.017 0.665 0.565 2.024 1.712 32.250 32.250 35.353 35.601 33. 767 33. 366 1.202 1.123 3.604 3.325 32.250 32.250 35.164 35.567 32.943 33. 812 0.992 1.161 3.010 3.461 33.184 33.398
L3 32.250 33.937 32.857 0.658 2.003 32.250 35.075 33.349 1.166 3.496 32.250 35.033 32.546 0.844 2.572 32.917
Kondisi awal 32.250 35.618 34.035 1.001 2.942 32.250 34.988 33.069 0.635 1.921 32.250 35.596 33.786 0.906 2.682 33.630
Tabel 15 Distribusi suhu udara hasil simulasi CFD pada beberapa disain kandang dengan arah angin dari kanan/kiri Z (m) 0.6
Nilai
Minimum (0C) Maksimum (0C) Rata-rata (0C) Standard deviasi Coefficient of variance (%) 1.2 Minimum (0C) Maksimum (0C) 0 Rata-rata ( C) Standard deviasi Coefficient of variance (%) 1.6 Minimum (0C) 0 Maksimum ( C) 0 Rata-rata ( C) Standard deviasi Coefficient of variance (%) Rata-rata pada z=0.6, 1.2 dan 1.6 m Rata-rata 2 arah angin
Tinggi Atap 5.25 m L1 L2 32.251 32.250 35.855 35.809 34.234 34.164 1.144 1.102 3.342 3.226 32.250 32.250 35.708 35.660 33.721 33.729 0.923 0.997 2.738 2.957 32.250 32.250 34.706 35.033 33.342 33.316 0.636 0.741 1.906 2.224 33.766 33.736 33.627 33.547
L3 32.242 35.788 34.135 1.088 3.186 32.250 35.664 33.663 0.954 2.833 32.250 34.652 33.262 0.698 2.099 33.687 33.418
Tinggi Atap 5.75 m L1 L2 32.223 32.250 35.890 35.800 34.242 34.193 1.116 1.117 3.263 3.263 32.250 32.250 35.167 35.565 33.708 33.571 0.808 0.899 2.405 2.666 32.250 32.250 34.183 34.533 33.303 33.256 0.558 0.617 1.677 1.853 33.751 33.673 33.703 33.438
L3 32.226 35.852 34.170 1.071 3.134 32.250 35.574 33.697 0.928 2.753 32.250 35.079 33.251 0.707 2.127 33.706 33.361
Tinggi Atap 6.25 m L1 L2 33.533 32.251 35.729 35.855 34.823 34.234 0.715 1.144 2.052 3.341 32.684 32.250 34.655 35.708 33.795 33.721 0.515 0.923 1.525 2.738 32.635 32.250 33.954 34.706 33.412 33.343 0.294 0.635 0.879 1.905 34.010 33.766 33.597 33.582
L3 32.250 35.775 34.264 1.080 3.152 32.250 35.427 33.693 0.859 2.549 32.250 34.720 33.250 0.625 1.879 33.736 33.327
Kondisi awal
34.03 34.52 34.24 0.12 0.34 33.38 34.42 33.97 0.35 1.05 32.83 34.03 33.71 0.32 0.95 33.97 33.800
Distribusi suhu udara dalam kandang pada 9 disain kandang simulasi selain dipengaruhi oleh bukaan ventilasi, kecepatan angin, juga dipengaruhi oleh efek termal yang terjadi di dalam kandang. Efek termal ini muncul karena panas yang dipancarkan oleh sapi melalui kulitnya menyebabkan suhu udara dalam kandang meningkat dan lebih tinggi dari suhu udara lingkungan. Kulit sapi perah memberikan kontribusi cukup tinggi terhadap
panas dalam kandang. Pada
simulasi ini suhu kulit sapi sebesar 36,2oC dengan heat transfer coeficient sebesar 55,41 W/m2.oC. . Dengan kecepatan angin yang masuk sebesar 0,7 m/detik, maka efek termal dalam kandang tidak dapat diabaikan. Hal ini sejalan dengan yang dikemukakan oleh Papadakids et.al. (1996) yang menyatakan bahwa kecepatan angin yang melebihi 1,8 m/detik efek termal terhadap laju ventilasi dapat diabaikan,
tetapi bila kecepatan angin lebih rendah dari 1,8 m/detik maka efek
termal tidak dapat diabaikan. Distribusi suhu udara dalam kandang di tiga ketinggian (z=0,6; 1,2 dan 1,6 m) pada 9 disain simulasi yang dipilih adalah distribusi suhu udara dalam kandang yang memiliki nilai mendekati suhu udara lingkungan. Berdasarkan Tabel 14 dan 15, distribusi suhu udara dalam kandang yang dipilih adalah disain kandang dengan tinggi atap kandang 6,25 m, lebar 8,3 m, tinggi dinding kanan dan kiri 0,4 m dengan posisi bak penampung air dipindahkan. Dipilihnya disain ini karena bukaan ventilasi yang dibuat telah mampu mereduksi panas dalam kandang dengan suhu udara rata-rata dalam kandang pada tiga ketinggian (z=0,6; 1,2 dan 1,6 m) sebesar 33,327oC (rata-rata pada dua arah angin). Distribusi suhu udara dalam kandang hasil simulasi pada arah angin dari depan/belakang dan kanan/kiri lebih tinggi dari suhu udara lingkungan (32,25oC) yang membuktikan bahwa efek termal dalam kandang tidak dapat diabaikan. Distribusi suhu udara hasil simulasi pada disain kandang terpilih memiliki nilai paling rendah jika dibandingkan dengan disain kandang lainnya. Disain terpilih hasil simulasi dengan tinggi kandang sebesar 6,25 m merupakan disain dengan tinggi atap kandang yang cukup besar jika dibandingkan dengan penelitian terdahulu (tanpa CFD) yaitu 3,6 – 4,2 m (Hahn, 1985), 2-3 m (McDowell, 1972), 3,5 m untuk atap yang terbuat dari seng (Basyarah, 1995).
Disain kandang hasil simulasi memiliki distribusi suhu udara dalam kandang rata-rata sebesar 33,327oC (33,561oC pada z=0,6; 33,521oC pada z=1,2 m
& 32,898oC pada z=1,6 m), lebih rendah
dari disain awal (Tabel 16).
o
Perbedaan terbesar pada ketinggian 1,6 m (0,85 C) akibat diturunkannya tembok kanan dan kiri dari 1,05 m menjadi 0,4 m dan dipindahkannya posisi bak penampung air yang menyebabkan bukaan ventilasi kandang lebih luas dan udara lingkungan dapat masuk ke kandang mulai dari ketinggian 0 m (arah angin dari depan/belakang) atau 0,4 m (arah angin dari kanan/kiri).
Perbedaan suhu udara
pada kondisi awal sebelum simulasi pada arah angin (inlet) dari kanan/kiri (33,97oC) dengan hasil simulasi (33,74oC) sebesar 0,23oC dan arah angin (inlet) depan/belakang (33,630oC
pada kondisi awal dan 32,917oC hasil simulasi)
sebesar 0,713oC dengan rataan pada dua arah angin (inlet) sebesar 0,474oC akan meningkatkan dry matter intake (DMI) sebesar 0,403 kg per hari per ekor (West, et.al., 2003). Distribusi suhu udara rata-rata dalam kandang (pada z=0,6; 1,2 & 1,6 m) sebesar 33,327oC menyebabkan sapi perah akan mengalami stress sedang (Wierema, 1990), kecuali kelembaban udara dalam kandang dapat dipertahankan di bawah 45%. Tabel 16 Distribusi suhu udara pada disain kandang terpilih dan kandang kondisi awal dengan inlet di kanan/kiri dan depan/belakang
Z (m) Uraian 0.6 Minimum (oC) 0.6 Maksimum (oC) 0.6 Rata-rata (oC) 0.6 Standard deviasi 0.6 Coefficient of variance (%) 1.2 Minimum (oC) 1.2 Maksimum (oC) 1.2 Rata-rata (oC) 1.2 Standard deviasi 1.2 Coefficient of variance (%) 1.6 Minimum (oC) 1.6 Maksimum (oC) 1.6 Rata-rata (oC) 1.6 Standard deviasi 1.6 Coefficient of variance (%) Rata-rata pada z=0,6; 1,2 dan 1,6 m
Inlet di kanan/kiri Hasil Disain simulasi awal 34.03 32.250 34.52 35.775 34.24 34.264 0.12 1.080 0.34 3.152 33.38 32.250 34.42 35.427 33.97 33.693 0.35 0.859 1.05 2.549 32.83 32.250 34.03 34.720 33.71 33.250 0.32 0.625 0.95 1.879 33.73 33.97
Inlet di depan/belakang Hasil Disain simulasi awal 32.250 33.937 32.857 0.658 2.003 32.250 35.075 33.349 1.166 3.496 32.250 35.033 32.546 0.844 2.572 32.917
32.250 35.618 34.035 1.001 2.942 32.250 34.988 33.069 0.635 1.921 32.250 35.596 33.786 0.906 2.682 33.630
Upaya yang dapat dilakukan untuk mengurangi stres pada sapi perah adalah pemberian air minum dingin, penyemprotan angin dingin ke seluruh tubuh ternak (Shibata, 1996), modifikasi disain kandang, pemilihan bahan atap kandang, penambahan kecepatan angin melalui kipas dan pemberian shelter di sekitar kandang. Pemberian air minum dingin dapat meningkatkan produksi susu sapi Holstein sebesar 10,86% dari 22,1 kg pada air minum 28oC menjadi 24,5 kg pada air minum 10oC (Milam, et. al., 1986). Wilks et. al., (1990) melaporkan bahwa terjadi kenaikan produksi susu sapi Holstein sebesar 4,85% dari 24,7 kg pada air minum 27oC menjadi 25,9 kg pada 10,6oC. Qisthon (1999) melaporkan bahwa pemberian air minum pada suhu 10oC dapat memperbaiki produktivitas sapi dara FH melalui pertambahan bobot badan dan efisiensi pakan, meningkatkan kecernaan bahan kering dan bahan organik pakan dibandingkan dengan pemberian air minum pada suhu 16, 22 dan 28oC. Penyemprotan air dingin ke seluruh tubuh ternak dapat menurunkan suhu udara dalam kandang. Penurunan suhu udara sekitar kandang sebesar 5oC dapat meningkatkan produksi susu sapi FH sebesar 10 kg/hari yaitu dari 35 kg/hari menjadi 45 kg/hari (Berman, 2005). Pemberian air minum dingin dan penyemprotan air dingin ke seluruh tubuh ternak akan meningkatkan RH dalam kandang sehingga diperlukan perhitungan dalam pemberiannya agar naiknya RH dalam kandang tetap memiliki nilai temperature humidity index (THI) yang lebih rendah sehingga tidak menyebabkan stres pada sapi. Pemilihan bahan atap kandang erat kaitannya dengan efektivitas bahan atap menghantarkan panas radiasi yang diterima ke dalam kandang agar radiasi yang sampai ke ternak rendah. Suhu udara dalam kandang yang atapnya terbuat dari asbes, seng dan rumbia (konduktivitas bahan 0,0001 kal/detoC) berturut-turut 26,5; 27,0 dan 26,4oC (Gatenby & Martawijaya, 1986). Respons fisiologis sapi perah sangat baik terhadap bahan atap kandang rumbia dibandingkan dengan genteng dan seng. Respons fisiologis ini dapat dilihat dari suhu tubuh, suhu rectal, suhu kulit, denyut jantung dan frekuensi nafas yang lebih rendah pada sapi FH yang diberi atap rumbia dibandingkan dengan yang diberi atap seng atau genteng (Soemarto, 1995).
Pemberian kecepatan angin tertentu dapat mereduksi panas dalam kandang. Pemberian kecepatan angin melalui terowongan angin (wind tunnel) yang dibuat dalam kandang dapat menurunkan suhu (4,2oC), menurunkan THI (6,0) dan meningkatkan RH (26%) dalam kandang (Smith et.al., 2005). Pemberian kecepatan angin sebesar 1,125 m/det pada siang hari (pukul 11:00 – 13:00 WIB) menyebabkan perubahan suhu rektal, suhu kulit, suhu tubuh dan frekuensi pernafasan pada sapi (Hadi, 1995).
Pemberian shelter di sekitar kandang dapat
menurunkan suhu udara lingkungan di sekitar kandang, mengurangi kecepatan angin yang masuk ke kandang, mengurangi semprotan air hujan dari sekitar kandang (Esmay, 1986)
