South East Asian Conference on Mathematics and Its Applications Proceedings Institut Teknologi Sepuluh Nopember, Surabaya ISBN 978 – 979 – 96152 – 5 – 1
APLIKASI BEDA HINGGA PADA PROSES STERILISASI MAKANAN KALENG 1
Lukman Hanafi, 2 Mardlijah dan 3 Dedik Ardian 1, 2 , 3
Jurusan Matematika Institut Teknologi Sepuluh Nopember Kampus ITS Keputih Sukolilo Surabaya, Indonesia email :
[email protected]
Abstrak. Proses sterilisasi makanan kaleng dengan pemanasan dilakukan untuk penentuan jaminan keselamatan bahan makanan. Jaminan ini meliputi tepat atau tidaknya tanggal kadaluwarsa yang tercantum pada label, rusak tidaknya kualitas makanan dan jumlah kandungan nutrisi yang ada. Namun demikian, proses sterilisasi makanan yang diberikan tidak semata-mata membunuh mikroba, tetapi juga harus mempertimbangkan mutu akhir dari produk, dimana kerusakan mutu oleh pemanasan harus diminimalkan. Dengan demikian, optimasi proses sterilisasi makanan kaleng diperlukan untuk dapat menentukan kombinasi suhu dan waktu selama pemanasan yang dapat memenuhi kriteria keamanan pangan dan mutu. Penyelesaian dan simulasi model dari permasalahan menggunakan metode beda hingga untuk mendapatkan solusi numeriknya. Metode beda hingga merupakan salah satu bentuk penyelesaian numerik yang digunakan untuk menyelesaikan persamaan diferensial. Kata kunci : sterilisasi, makanan kaleng, beda hingga
1. Pendahuluan Perkembangan ilmu pengetahuan dan teknologi yang ada saat ini memerlukan suatu solusi yang tepat dari permasalahan yang ada, terutama dalam bidang industri.. Dalam kehidupan sehari-hari banyak dijumpai hal-hal yang berkaitan dengan perpindahan panas, terutama dalam bidang industri. Perpindahan panas atau heat transfer adalah ilmu yang meramalkan perpindahan energi yang terjadi karena adanya perbedaan temperatur diantara benda atau material. Ilmu perpindahan panas tidak hanya mencoba menjelaskan bagaimana energi panas itu berpindah dari suatu benda ke benda lainya, tetapi juga dapat meramalkan laju perpindahan yang terjadi pada kondisi-kondisi tertentu. salah satu dari beragam masalah yang dapat diselesaikan dengan ilmu perpindahan panas adalah persamaan panas pada proses sterilisasi makanan kaleng.
2. Proses sterilisasi makanan kaleng Proses sterilisasi merupakan proses utama pada proses produksi pangan steril komersial, khususnya untuk menjamin tercapainya keamanan pangan steril komersial. sterilisasi (processing) pada pengalengan adalah proses pemanasan wadah serta isinya pada suhu dan jangka waktu tertentu untuk menghilangkan atau mengurangi faktor-faktor penyebab kerusakan makanan, tanpa menimbulkan gejala lewat pemasakan (over cooking) pada makanannya. Waktu dan suhu yang diperlukan untuk proses sterilisasi biasanya tergantung pada konsistensi atau ukuran partikel bahan, derajat keasaman isi kaleng, ukuran headspace, besar dan ukuran kaleng, kemurnian uap air (steam) yang digunakan, dan kecepatan perambatan panas. Mikroba mempunyai ketahanan panas yang berbeda - beda,. Tabel 1 memperlihatkan kombinasi suhu dan waktu yang diperlukan untuk membunuh mikroba. Sel vegetatif khamir dan kapang dapat
(60 880 C ) , sedangkan bakteri termofilik dan mesofilik 0 perlu suhu yang lebih tinggi untuk membunuhnya (umumnya pada suhu standar 121 C ). diinaktifkan pada suhu yang lebih rendah
69
70
Aplikasi Beda Hingga Pada Proses Sterilisasi Makanan Kaleng Tabel 1. Kombinasi suhu dan waktu yang diperlukan untuk menurunkan jumlah mikroba pada level yang sama.
Waktu (menit) 10 5 30 - 60
Organisme Sel Vegetatif Khamir Kapang Bakteri Termofilik : Clostridium Thermosacahharilyticum Bacilus
Suhu
(0 C)
80 60 88
3–4
121
4
121
3
121
0,1 – 1
121
1,5
121
0,6 menit
121
Bakteri Mesofilik Clostridium Botulinum Toksin Botulinum A & B Clostridinum Sporogenes Bacillus Subtillis
Apabila suspensi mikroba dipanaskan pada suhu konstan, maka penurunan jumlah mikroba mengikuti pola logaritmik sebagai fungsi dari waktu. Pada suhu tertentu,laju inaktifasi mikroba selama waktu pemanasan pada suhu tertentu dapat dinyatakan sebagai berikut :
dN kN dt
(1)
dimana : N : Jumlah mikroba sisa yang masih hidup setelah pemanasan t : Waktu pemanasan k : laju reaksi Jka persamaan (1) diintegralkan, maka diperoleh persamaan berikut :
ln
N kt N0
(2)
Persamaan (2) menunjukkan plot kurva semilogaritma N terhadap t persamaan tersebut dapat dirubah menjadi lebih sederhana 2.303 log
N N kt kt atau log N0 N0 2.303
Nilai slope 2.303/k sering dinyatakan dengan nilai D, sehingga :
log
N t N0 D
Nilai D = 2.303/k, yaitu waktu penurunan desimal (decimal reduction time). Nilai D dipengaruhi oleh suhu, semakin tinggi suhu pemanasan, maka waktu yang diperlukan untuk menginaktivasi mikroba akan semakin pendek. Gambar 1 memperlihatkan pengaruh suhu terhadap nilai D. Pada suhu pemanasan yang lebih tinggi kurva akan semakin curam yang memberikan slope kurva -1/D semakin besar atau nilai D semakin kecil.
Lukman Hanafi, Mardlijah, Dedik Ardian
71
Gambar 1 Pengaruh suhu terhadap nilai D
Nilai D dari setiap mikroba memiliki sensitivitas yang berbeda terhadap perubahan suhu.. Sensitifitas nilai D terhadap suhu sering dinyatakan dengan nilai Z, yaitu perubahan suhu yang diperlukan untuk merubah nilai D sebesar 90% atau 1 siklus logaritma. Misalnya, Clostridium botulinum memiliki nilai Z 0
0
sebesar 10 C , artinya untuk merubah nilai D mikroba tersebut dari 0.25 menit pada suhu 121 C menjadi 0.025 menit (menurun sebesar 90% atau satu siklus logaritma), suhu pemanasan harus dinaikkan 0 0 sebesar 10 C , yaitu menjadi 131 C .
Gambar 2 Kurva semi-logaritmik hubungan antar nilai D dengan suhu.
Kurva semi logaritmik ini berbentuk linear dengan slope-nya adalah -1/Z secara matematis, nilai DT pada suhu tertentu dapat dinyatakan dengan persamaan berikut :
log
D T Tref DT Z
Dimana : D : waktu yang dibutuhkan untuk mengurangi bakteri sebanyak 90 persen (menit) DT : nilai D pada suhu tertentu (menit) T
Tref Z
0
: suhu pemanasan tertentu ( C /
0
F) 0
: suhu standar yang digunakan untuk nilai D ( C / 0
: sensitifitas nilai D terhadap suhu ( C /
0
F)
0
F)
72
Aplikasi Beda Hingga Pada Proses Sterilisasi Makanan Kaleng
3. Model Matematika dan Metode Beda Hingga Proses sterilisasi makanan kaleng memiliki serangkaian kendala yang harus diperhatikan. Kendala yang ada diantaranya adalah proses perambatan panas yang dimodelkan dengan persamaan difrensial parsial orde dua dan kendala yang lain berupa nilai sterilisasi dengan memperhatikan derajat pertumbuhan mikroorganisme dalam makanan kaleng. Proses pemodelan persamaan panas dimensi tiga disini dimulai dari sistem koordinat Cartesius kemudian ditransformasikan kedalam sistem koordinat tabung. Pemakaian sistem koordianat tabung karena lebih sesuai dengan masalah yang akan diselesaikan, yaitu kaleng berbentuk tabung. Selanjutnya akan dicari model matematis perambatan panas yang terjadi. Namun untuk selanjutnya pemanasan tergantung dari jenis bahannya yang diamati, seperti makanan kaleng tergantung jenis makanannya yang terlihat dari kalor jenis bahan c, konduktifitas suhu bahan k dan masa jenis bahan ρ. Model persamaan panas adalah :
1 T 2T 2T 2T 2 2 2 2T t x y z 1 c dimana adalah konstanta penghamburan panas d k
1 L 2 Z
X
R
Y
1 L 2
Gambar 3 Kaleng berbentuk tabung
Selanjutnya jika T=T(x,y,z,t) ditransformasikan dalam koordinat tabung T =T(r, ,z,t) seperti terlihat pada gambar 3 dengan transformasi sebagai berikut x = r cos , y = r sin dan Z = z, akan diperoleh :
1 T 1 T 2T 1 2T 2T t r r r 2 r 2 2 z 2 Karena kaleng yang diperlihatkan berbentuk tabung yang bersifat simetri sehingga perambatan panas tidak bergantung pada besar sudut .
Lukman Hanafi, Mardlijah, Dedik Ardian
73 Y
R
R
R X
R
Gambar 4 Perambatan panas pada sisi kaleng
Persamaan panas pada kaleng yang berbentuk tabung, dengan T(r,z,t) adalah suhu bahan makanan di lokasi kaleng (r,z) dengan r( 0 r R ) dan z( L / 2 z L / 2 ) saat waktu t . Untuk melengkapi pemodelan kendala yang diberikan pada perambatan panas, ditambahkan syarat awal pada waktu initial time t i suhu awalnya adalah T0
T (r , z, t i ) T0 r [0, R)
L L z [ , ] 2 2
Sedangkan syarat batas suhu di sekeliling, dasar dan atas kaleng diberikan oleh T [r,z,t] = u(t)
L , t ] = u(t) 2 L T [ r , , t ] = u(t) 2 T [ r,
z [
L L , ] 2 2
r [0, R] r [0, R]
sedangkan suhu di sepanjang sumbu z, yaitu r = 0 tidak berubah terhadap r Penyelesaian dan simulasi model dari permasalahan di atas menggunakan metode beda hingga untuk menentukan solusi numeriknya. Konsep yang digunakan pada metode beda hingga adalah mengubah sistem persamaan difrensial parsial ke sistem persamaan diferensial biasa dengan mendiskritkan ruang domainnya. Ruang (r,z,t) didiskritisasi dengan N titik (i=1,2,3,...,N) untuk dimensi r, M titik (j=1,2,3,...,M) untuk dimensi z dan P titik (k=1,2,3,..,P) untuk dimensi t dengan
ri R (i 1)r L z j ( j 1)z 2 t k (k 1)t dengan
R N 1 L z M 1 T t P 1 r
74
Aplikasi Beda Hingga Pada Proses Sterilisasi Makanan Kaleng Persamaan beda hingga digunakan untuk pendekatan turunan parsial dari u(r,z,t) dengan
T (ri r , z j , t k ) u (ri r , z j , t k ) T r, z j , t k r 2 r u (ri r , z j , t k ) 2T (ri , z j , t k ) T (ri r , z j , t k ) 2T r , z , t i j k r 2 2 r
T (ri , z j z, t k ) 2T (ri , z j , t k ) T (ri , z j z , t k ) 2T ri , z j , t k 2 z 2 z u (ri , z j , t k t ) u (ri , z j , t k ) T ri , z j , t k t t sehingga didapat
1 Ti , j ,k 1 Ti , j ,k Ti 1, j ,k Ti 1, j ,k Ti 1, j ,k 2Ti , j ,k Ti 1, j ,k Ti , j 1,k 2Ti , j ,k Ti , j 1,k t 2rr 2 r 2 z untuk i = 2,3,4,...,N-1, j = 2,3,4,...,M-1, k = 1,2,3,...,P Sedangkan untuk syarat batasnya, diberikan oleh Ti , j ,k T (t ) j = 1,2,3,...,M, k = 2,3,...,
Ti , j ,k T (t )
i = 1,2,3,...,N
k = 2,3,...,P
Ti ,M ,k T (t ) i = 1,2,3,...,N k = 2,3,...,P Pada pusat kaleng dan sepanjang tinggi kaleng gradien
T 0 sehingga untuk r = 0, (i = N), r
T ( R r , z j , t k ) = T ( R r , z j , t k ) atau TN 1, j ,k TN 1, j ,k dan kemudian disubstitusikan didapat
2t t TN 1, j ,k TN , j ,k 2 TN , j 1,k 2TN , j ,k TN , j 1,k 2 r z Misalkan diberikan persamaan panas dengan syarat awalnya T0 0 (dapat pula dengan nilai besaran TN , j ,k 1 TN , j ,k
yang lain, misal = 100), konstanta 1,2 10 m / s , jari-jari kaleng R = 0,025 m dan tinggi kaleng L= 0,06 m, maka solusi numerik pada pusat kaleng diberikan oleh gambar 5 berikut 7
2
Gambar 5 Solusi numerik
Lukman Hanafi, Mardlijah, Dedik Ardian
75
4. Pemodelan nilai sterilisasi Fenomena lain dari sterilisasi makanan kaleng yang sangat penting adalah penentuan nilai sterilisasi pada akhir waktu. Nilai sterilisasi ini menentukan derajat kehancuran mikroorganisme didalam bahan makanan dan relatif terhadap konsentrasi mikroorganismenya, nilai sterilisasi disini diberikan dengan besaran F0 jika c(r,z,t) adalah konsentrasi dari mikroorganisme pada lokasi (r,z) saat waktu t, maka c(0,0,t) adalah konsentrasi mikroorganisme pada pusat kaleng saat waktu t. Seperti yang disebutkan di atas, nilai sterilisasi ini lebih dilihat pada pusat kaleng karena tempat ini lebih lambat menerima panas sehingga kemungkinan besar mempunyai konsentrasi mikroorganisme terbesar. Kehancuran mikroorganisme akibat pemanasan diberikan oleh persamaan difrensial berikut :
dc0,0, t kc0,0, t dt c0,0, t i c0
(3)
c 0 adalah konsentrasi awal, parameter k merupakan sebuah fungsi dari suhu. Ketergantungan fungsi ini dapat dilihat pada parameter jumlah waktu yang dibutuhkan untuk mereduksi konsentrasi Dengan
mikroorganisme.
In(10) k
dan suhu yang dihasilkan ref dibutuhkan untuk mereduksi besar dari sebuah suhu acuan
yang berhubungan dengan
Tref oleh sebuah faktor dari 10
T (0,0, t ) Tref ref exp ln 10 ref
dengan T(0, 0, t) adalah suhu pada pusat kaleng. sshingga didapat nilai
k
ln(10)
ref
T (0,0, t ) Tref exp ln 10 ref
k
jika persamaan (3) diselesaikan maka didapat t1
ln ct ln c0 k dt t0
c ln t c0
ln(10) dt t0 ln 10 t1 T 0,0, t Tref c0,0, t exp exp ln(10) c0 ref ref t0 t1
dt
Pada persamaan diatas, nilai sterilisasi pada waktu t didefinisikan dengan
ln(10) t1 F0 t exp T 0,0, t Tref dt t0 ref Jika diberikan derajat sterilisasi F yang diinginkan pada akhir waktu, maka harus dipastikan
F0 t f F
begitu pula untuk suhu akhir harus
T 0,0, t f T f
dan konsentrasi dari mikroorganisme saat akhir waktu harus memenuhi:
76
Aplikasi Beda Hingga Pada Proses Sterilisasi Makanan Kaleng
c
0,0, t c f
0
ln(10) exp F ref
Misalkan diberikan persamaan panas sebagai titik pusat kaleng dengan syarat awalnya
T0 0 (dapat
pula dengan nilai besaran yang lain, misal = 100), konstanta 1,2 10 m / s , jari-jari kaleng R = 0,025 m dan tinggi kaleng L= 0,06 m, dengan grid tinggi = 3, grid jari-jari = 2 dan grid waktu = 10 maka kondisi pada titik pusat kaleng ketika dipanaskan selama 10 menit ditunjukkan oleh gambar 6 7
2
Gambar 6 Kondisi titik pusat kaleng setelah dipanaskan selama 10 menit
5. Kesimpulan 1.
Proses perambatan panas pada proses sterilisasi makanan kaleng yang berbentuk tabung dapat dituliskan dalam bentuk PDP orde dua sebagai berikut. Jika T(r,z,t) adalah suhu bahan makanan di lokasi kaleng (r,z) dengan r (0 r R) dan Z (
L L Z ) saat waktu t, maka persamaan 2 2
1 T 1 T 2T 2T panasnya adalah dengan syarat awal dan syarat batas : t r r r 2 z 2 L L L L L T(r,z, t i ) = T0 r 0, R Z , T(r, ,t) = T(t) r , 2 2 , 2 2 , 2 L T L L T(r, ,t) = T(t) r 0, R, 0 r = 0 Z , , r 2 2 2 7 2 2. Dari simulasi model yang dilakukan untuk beberapa parameter 1.2 x10 m / s R = 0,025m dan L = 0,06m dapat disimpulkan bahwa suhu pada pusat kaleng dengan kondisi syarat awalnya bernilai nol dan syarat batasnya bernilai satu tak melewati suhu syarat batasnya sepanjang pergerakan waktunya.
Lukman Hanafi, Mardlijah, Dedik Ardian
77
Daftar Pustaka 1. Basarudin T, 1994, Metode Beda hingga untuk persamaan difrensial, Elex Media Komputindo, Jakarta 2. Holman J.P., Jasjfi E.1994, Perpindahan kalor. Erlangga, Jakarta 3. Kannan A, Gourisankar P.Ch., Sandaka 2008 . Heat Transfer Analysis of canned food Sterilization in a Still Retort. Journal of Food Engineering 88 213-228. 4. Kreith F, Prijono A 1997, Prinsip-prinsip Perpindahan Panas, Edisi ketiga. Erlangga, Jakarta 5. Kusnandar F ,dkk . 2008. Aspek Mikrobiologi Makanan Kaleng http://www.unhas.ac.id/gdln/dirpan/pengalengan/Topik6/modul/ diakses tanggal 28 Februari 2009. 6. Soehardjo, 1996, Matematika 4, ITS, Surabaya.