## PEMODELAN KADAR AIR PADA SIFAT FISIK STABILISASI TANAH GAMBUT
Mohammad Syaiful Pradana 1 , Awawin Mustana Rohmah 2
1,2 Universitas Islam Darul Ulum Lamongan
Abstract: Peat  soil  is  an  organic  soil  with  very  low  carrying  capacity  and  high compressibility.  The  condition  is  less  profitable  for  civil  engineers  in  building  a  civil foundation  foundation.  One  method  of  peat  soil  improvement  can  be  done  with astabilization method that more environment-friendly and cheaper than other methods. Laboratory  based  peat  stabilization  studies  to  increase  carrying  capacity,  reduce compression  and  improve  peat  soil  physical  properties  have  been  conducted  in Indonesia. The results of laboratory studies shown in the graph are still limited by time and content of the mixture. Therefore, further research is needed on the mathematical model toward the physical properties of peat soil stabilization. In this research will be formulated mathematical model of water content on the physical properties of peat soil stabilization. The model is derived from the fluid equation through porous media based on  the  principle  of  continum  and  controlvolume.  The  model  is  then  resolved numerically  by  different  method  until  MacCormack  scheme  with  two  steps  are predictor  step  using  forward  difference  and  correctorstep  using  backward  difference. The  MacCormack  scheme  has  the  advantage  of  solving  fluid  flow  equations  and continuity.  The  model  is  then  simulated  and  validated  by  comparing  the  simulation results  with  the  real  system.  From  the  simulation  results  obtained  the  water  content gradually decreased, the decrease is almost close to zero. In addition, it can be seen the difference  in  decrease  in  moisture  content  at  each  test  point  although  in  small quantities.
Keywords: Moisture Content, MacCormack, Peat Soil Stabilization.

## PENDAHULUAN
Tanah gambut merupakan tanah organik dengan daya dukung sangat rendah dan kompresibilitas tinggi.  Kondisi  tersebut  kurang  menguntungkan  bagi  para civil  engineer dalam  membangun  pondasi  bangunan  sipil  diatasnya  sehingga  menyebabkan  metode perbaikan  tanah  diterapkan  untuk  meningkatkan  daya  dukung  dan  meningkatkan kompressibilitas  (Yulianto  FE,  2016).  Salah  satu  metode  yang  lebih  berwawasan lingkungan dan murah dibandingkan metode lain yaitu metode stabilisasi. (Ma‟ruf MA, 2016).
Metode  stabilisasi  dikembangkan  dengan  menggunakan  bahan  aditif  seperti campuran  kapur-abu  terbang  atau  campuran  kapur-abu  sekam.  Metode  stabilisasi  untuk tanah  gambut  tersebut  dikembangkan  oleh  Yulianto,  F.E  (2012)  menggunakan  bahan campuran  ( admixture  stabilizer )  70%  abu  sekam  padi  ( Rice  Husk  Ash )  dan  30%  kapur ( Lime )  didapatkan  gambut  yang  distabilkan  dengan  10%  bahan  campuran  memberikan kenaikan kapasitas bantalan optimum dan penurunan kompresibilitas sehingga perbaikan menjadi  lebih  menonjol  seiring  bertambahnya  waktu.  Pengembangan  metode  stabilisasi lain  dilakukan  oleh  Ma‟ruf,  M.A  (2015)  menyimpulkan  bertambahnya  umur  stabilisasi maka reaksi bahan stabilisasi akan semakin lama pula.
Penelitian  stabilisasi  tanah  gambut  berbasis  laboratorium  telah  dilakukan  di Indonesia.  Hasil  penelitian  laboratorium  yang  ditampilkan  dalam  grafik  masih  dibatasi oleh waktu dan kadar bahan campuran.Oleh karena itu, dibutuhkan penelitian lebih lanjut mengenai model matematika ke arah sifat fisik stabilisasi tanah gambut, misalnya kadar air  dan  angka  pori  tanah  gambut.  Menurut  Yulianto,  F.E  (2016)  penurunan  kadar  air menyebabkan  angka  pori  menjadi  lebih  kecil  sehingga  kuat  geser  meningkat  dan kompressibilitas menurun.
Berdasarkan pada uraian diatas, pada penelitian ini diturunkan model matematika untuk kadar air pada sifat fisik stabilisasi tanah gambut. Model diturunkan dari persamaan fluida  melalui  media  berpori  berdasarkan  prinsip  kontinum  dan  volume  kendali. Selanjutnya  model  diselesaikan  secara  numerik  dengan  metode  beda  hingga  skema MacCormack dengan dua langkah yaitu langkah prediktor menggunakan beda maju dan langkah  korektor  menggunakan  beda  mundur.  Skema  MacCormack  memiliki  kelebihan untuk  menyelesaikan  persamaan  aliran  fluida  dan  kontinuitas.  Selanjutnya  model disimulasikan dan divalidasi dengan membandingkan hasil simulasi dengan sistem nyata.

## METODE PENELITIAN
Parameter  tanah  gambut  yang  penting  untuk  menentukan  sifat  fisiknya  adalah kadar  air,  angka  pori,  spesifik  gravity,  kadar  organik,  kadar  abu,  keasaman,  kadar  gas, rembesan,  susut,  batas  konsistensi  dan  berat  volume.  Sifat  fisik  tersebut  sangat dipengaruhi  oleh  material  organik  dan  proses  terbentuknya  tanah  gambut.  Dalam penelitian ini, sifat fisik yang akan dirumuskan model matematisnya yaitu kadar air. A. Kadar Air
Kadar air tanah gambut asli dalam kondisi jenuh dapat melebihi 600%, namun secara umum  variasi  kadar  airnya  berkisar  antara  500-1500%  seperti  pengujian  yang  telah dilakukan  oleh  Yulianto,  F.E  (2016),  tetapi  kadar  air  tersebut  akan  berubah  secara drastis apabila tercampur dengan bahan anorganik meskipun dalam kadar yang kecil.
B. Stabilisasi Tanah Gambut
Metode merupakan suatu cara perbaikan tanah dengan mencampurkan bahan kimia ke dalam tanah dengan tujuan agar dapat meningkatkan sifat fisik dan tekniknya. Metode stabilisasi  banyak  diaplikasikan  pada  tanah  lempung  yang  memiliki  daya  dukung rendah  atau  tanah  dengan  nilai  plastisitas  tinggi (Expasive  Soil) .  Material  stabilisasi sering menggunakan kapur hidup CaO (Quick Lime) , kapur padam Ca(OH)2 (Calsium
Hydroxide) , abu sekam padi (Rice Husk Ash) dan abu terbang (Fly Ash) .
Model  matematika  dari  percobaan  stabilisasi  tanah  gambut  di  laboratorium diturunkan  berdasarkan  prinsip  kontinum  dan  volume  kendali  menggunakan persamaan kekekalan massa dan persamaan kekekalan momentum pada media berpori berupa  tanah  gambut  yang  telah  mendapatkan  perlakuan  seperti  yang  ditunjukkan Gambar 1.

## Gambar 1. Model Stabilisasi Tanah Gambut
Dari Gambar 1 sketsa model stabilisasi tanah gambut di laboratorium dihasilkan model  matematika untuk mengetahui prilaku pertumbuhan kristal yang dibentuk dari gambut, kapur, abu terbang (fly ash) dan air. Terbentuknya kristal tersebut akibat dari reaksi kimia dari bahan campuran (kapur, fly ash dan air) (Pradana & Widodo, 2014).

## C. Persamaan Konservasi Massa
Prinsip fisik pertama yang diterapkan untuk konversi dari sistem ke deskripsi volume kontrol adalah  prinsip  konservasi  massa  yaitu  massa sistem  tetap  konstan (Pritchard, 2011).
𝑑𝑀 𝑑𝑡 𝑠𝑖𝑠𝑡𝑒𝑚 = 0 (1)
Dengan
𝑀 𝑠𝑖𝑠𝑡𝑒𝑚 = 𝑑𝑚 𝑀(𝑠𝑖𝑠𝑡𝑒𝑚 ) = 𝜌𝑑𝑉 𝑉(𝑠𝑖𝑠𝑡𝑒𝑚 )
Formulasi sistem dan kontrol volume ditunjukkan pada (2) berikut
𝑑𝑁 𝑑𝑡 𝑠𝑖𝑠𝑡𝑒𝑚 = 𝜕 𝜕𝑡 𝜂 𝜌 𝑑𝑉 𝐶𝑉 +   𝜂 𝜌 𝑉  . 𝑑𝐴 𝐶𝑆 (2)
Dengan
𝑁 𝑠𝑖𝑠𝑡𝑒𝑚 = 𝜂 𝑑𝑚 𝑀(𝑠𝑖𝑠𝑡𝑒𝑚 ) = 𝜂 𝜌 𝑑𝑉
𝑉(𝑠𝑖𝑠𝑡𝑒𝑚 )
Untuk  mendapatkan  formulasi  kontrol  volume  dari  konservasi  massa,  di  misalkan 𝑁 = 𝑀 dan 𝜂 = . Selanjutnya dengan subtitusikannya ke (2) didapatkan
Gambut tanpa
perlakuan Gambut tanpa perlakuan Area Gambut yang di stabilisasi menggunakan campuran kapur dan abu terbang
𝑑𝑀 𝑑𝑡 𝑠𝑖𝑠𝑡𝑒𝑚 = 𝜕 𝜕𝑡 𝜌 𝑑𝑉 𝐶𝑉 +   𝜌 𝑉   . 𝑑𝐴 𝐶𝑆 (3)
Selanjutnya dengan membandingkan (1) dan (3) didapatkan formulasi volume kontrol dari konservasi massa sebagai berikut.
𝜕 𝜕𝑡 𝜌 𝑑𝑉 𝐶𝑉 +   𝜌 𝑉   . 𝑑𝐴 𝐶𝑆 = 0 (4)

## D. Persamaan Darcy
Persamaan yang menggambarkan mengenai aliran fluida pada media berpori pertama kali  dikembangkan  oleh  Henry  Darcy  pada  tahun  1856.  Persamaan  tersebut merupakan  persamaan  yang  menunjukkan  kecepatan  aliran  fluida  dengan permeabilitas  batuan,  viskositas  fluida  serta  gradient  tekanan  antar  jarak  tempuh aliran. Penyelidikan mengenai aliran fluida dalam media berpori berupa ground water juga  dilakukan  oleh  Syahruddin,  M.H  (2014)  menggunakan  persamaan  Darcy  (5) berikut.
𝑢 = 𝑞 𝐴 = 𝑘 𝑕 1 − 𝑕 2 𝐿 = 𝑘 ∆𝑕 𝐿 (5)
Dengan 𝑢 adalah kecepan aliran (𝐿/𝑆) , 𝑞 adalah debit aliran (𝐿 3 /𝑆) , A adalah luas  penampang (𝐿 2 ) , ∆𝑕 adalah  selisih  permukaan  fluida (𝐿 2 ) , 𝐿 adalah  panjang media (𝐿) dan 𝑘 adalah  konstanta  konduktivitas.  Nilai 𝑘 bergantung  pada  kondisi sifat padatan (𝑠𝑜𝑙𝑖𝑑 𝑚𝑎𝑡𝑟𝑖𝑥) dan sifat dari cairannya yaitu air.
E. Metode Beda Hingga Skema MacCormack
Dalam metode numerik terdapat beberapa metode salah satunya adalah Metode Beda Hingga  ( Finite  Difference  Method )  dengan  berbagai  macam  skema  eksplisit  dan implisit.  Contoh  skema  eksplisit  yaitu  skema  pengepingan  MacCormack. Natakusumah  (2007)  menggunakan  skema  pengepingan  MacCormack  pada persamaan aliran 1 dimensi ( Saint Venant Equation ) . Untuk lebih jelasnya dapat dilihat persamaanorde pertama (6) berikut.
𝜕𝑢 𝜕𝑡 + 𝑎 𝜕𝑢 𝜕𝑥 = 0 (6)
Tahap prediktor :  pada  tahap prediktor menggunakan  skema  maju,  nilai  sementara dari 𝑢 pada watu 𝑛 + 1 didenotasikan oleh 𝑢 𝑖 𝑛+1 sehingga diperoleh (7) berikut.
𝑢 𝑖 𝑛+1 = 𝑢 𝑖 𝑛 − 𝑎 ∆𝑡 ∆𝑥 𝑢 𝑖+1 𝑛 − 𝑢 𝑖 𝑛 (7)
Tahap korektor :  pada  tahap korektor menggunakan  skema  mundur,  nilai prediktor 𝑢 𝑖 𝑛+1 dikoreksi berdasarkan (8) berikut.
𝑢 𝑖 𝑛+1 = 𝑢 𝑖 𝑛+ 1 2 − 𝑎 ∆𝑡 2∆𝑥 𝑢 𝑖 𝑛+1 − 𝑢 𝑖−1 𝑛+1 (8)
Dengan  step  waktu  yang  digunakan  dalam  tahap korektor adalah ∆𝑡 2 dan  mengganti bentuk 𝑢 𝑖 𝑛+ 1 2 dengan rata-rata sementara
𝑢 𝑖 𝑛+ 1 2 = 𝑢 𝑖 𝑛 + 𝑢 𝑖 𝑛+1
2
Sehingga menghasilkan tahap korektor sebagai berikut
𝑢 𝑖 𝑛+1 = 𝑢 𝑖 𝑛 + 𝑢 𝑖 𝑛+1 2 − 𝑎 ∆𝑡 2∆𝑥 𝑢 𝑖 𝑛+1 − 𝑢 𝑖−1 𝑛+1 (9)
Selanjutnya untuk hasil korektor, persamaan (9) dapat ditulis menjadi :
𝑢 𝑖 𝑛+1 = 1 2 𝑢 𝑖 𝑛 +  𝑢 𝑖 𝑛+1 − 𝑎 ∆𝑡 2∆𝑥 𝑢 𝑖 𝑛+1 − 𝑢 𝑖−1 𝑛+1 (10)
Penggunaan  metode  beda  hingga  skema  MacCormack  memiliki  kelebihan  untuk menyelesaikan  persamaan  aliran  fluida,  berguna  dalam  memecahkan  persamaan kontinuitas  dan  persamaan  momentum  dua  dimensi  yang  merupakan  persamaan diferensial parsial non-linier dalam ruang dan waktu (Natakusumah, 2007).

## HASIL DAN PEMBAHASAN
Pada  stabilisasi  tanah  gambut  menggunakan  bahan  campuran  kapur, fly  ash dan air.  Pencampuran  tersebut  menyebabkan  terjadinya  reaksi  kimia  yang  membuat meningkat  kepadatannya  sehingga  daya  dukung  tanah  gambut  meningkat.  Berikut persamaan reaksi kimia (11) yang dihasilkan.
3 2 2 2 ( ) CaCO H O Ca OH CO    2 2 2 3
SiO H O H SiO
  2 2 3 3
2 ( ) 2 ( ) 2 Ca OH H SiO Ca HSiO H O   
(11)
Dengan 2 ( ) Ca OH adalah  kalsium  hidroksida, 2 CO adalah  karbon  dioksida,
2 3 H SiO adalah asam silicid.
Tabel 1. Reaksi Stoikiometri Senyawa Koefisien Molar 2 ( ) Ca OH 1 74,09 2 3 H SiO 2 78,10 3 2 ( ) Ca HSiO 1 194,26 2 H O 2 18,02 dengan       3 2 2 2 3 ( ) 3500,57 0, 60 ( ) 5786, 43 Ca HSiO H O k Ca OH H SiO   
Laju  reaksi  dipertimbangkan  karena  pentingnya  kemampuan  untuk  meramalkan kecepatan  campuran  reaksi  mendekati  kesetimbangan.  Laju  reaksi  dapat  ditentukan dengan  cara  mengikuti  perubahan  sifat  selama  terjadinya  reaksi.  Dengan  menganalisis campuran  reaksi  dalam  selang  waktu  tertentu,  maka  konsentrasi  reaktan  dan  produk reaksi  dapat  dihitung.  Persamaan  laju  reaksi  ditentukan  berdasarkan  konsentrasi  awal setiap zat, dipangkatkan orde reaksinya. Orde reaksi berpengaruh pada laju reaksi jika
[ ] d A v dt   Persamaan laju reaksinya (12) adalah     x y
v k A B  (12)
x dan y adalah orde reaksinya, untuk setiap konsentrasi yang sama pada produk reaksi dinyatakan :
Untuk orde reaksi orde ke nol :
produk reaksi A    0 v k A k
 
Untuk orde reaksi orde satu :
produk reaksi
A 
  v k A 
Untuk orde reaksi orde dua :
2 produk reaksi A 
  2 v k A 
Karena  laju  reaksi  berbanding  lurus  dengan  konsentrasi-konsentrasi  zat  pereaksi,  maka dapat menggunakan laju terintegrasi pada (13) berikut.
    n t d A k A dt       2 t d A kdt A     (13)
Dengan 𝑛 = orde  reaksi, 𝑘 = konstanta  laju  reaksi, 𝑡 = waktu.  Untuk  orde  reaksi  orde dua yang telah diintegralkan diperoleh
  1 t kt c A   (14)
Berthoud  (1912)  dan  Valeton  (1924)  menggambarkan  model  perubahan  massa dengan pertumbuhan dua tahap, yaitu proses difusi, dimana molekul solute berpindah dari bulk fase  liquid  ke  permukaan  solid.  Daya  dorong  terjadinya  kedua  tahap  ini  adalah perbedaan konsentrasi yang dapat ditunjukkan oleh (15) berikut:
( ) eq dm kA c c dt   (15)
Persamaan  (15)  dapat  dibangun  persamaan  massa  padatan  dengan 𝑚 = 𝜌 𝑠 ∀ . Karena gambut merupakan media berpori, besarnya nilai porositas berpengaruh terhadap besarnya  nilai  volume  padatan  dalam  gambut  tersebut,  volume  padatan  tersebut  dapat dituliskan  sebagai ∀= (1 − 𝜙) ,  konstanta  laju  reaksi k dan  Diasumsikan  bahwa  suatu konstan pada waktu 𝑡 diperoleh
( ) (1 ) eq s kA c c d dt      (16)
Persamaan  (16)  diketahui  bahwa  perubahan  massa  padatan  gambut  dipengaruhi oleh  perubahan  laju  pertumbuhan,  perubahan  konsentrasi  serta  perubahan  porositas  dari gambut. Perubahan tersebut dapat didapatkan dengan menggunakan Persamaan kekekalan massa dan momentum berdasarkan volume kendali melalui media berpori.
Menurut  prinsip  kontinum  hukum  kekekalan  massa  disebutkan  bahwa  aliran massa yang masuk dikurangi dengan aliran massa yang keluar sama dengan jumlah massa yang disimpan.

## Gambar 2. Volume Kendali Aliran Massa
Jika  volume  kendali  pada  gambar  2  diperbesar  dan  digambarkan  arah  aliran massanya  dari  permukaan  yang  berjarak x  dari  pusat  elemen  dinyatakan  berturut-turut sebagai ( ) ( ) ( ) x u x A x  dan ( ) ( ) ( ) x x u x x A x x        ,  maka  dari  hukum kekekalan massa diatasdapat dijabarkan menjadi    

 ( ( ) ( ) ( )) ( ( ) ( ) ( )) ( ( ) ( ) ( )) ( ( ) ( ) ( )) ( ( ) ( ) ( )) ( ( ) ( ) ( )) ( ) ( ) e w n s t b x u x A x x x u x x A x x t y v y A y y y v y y A y y t z w z A z z z w z z A z z t m t t m t  
                                  (17) Dengan ( ) ( ) , ( ) ( ) , ( ) ( ) e w n s t b A x A x x y z A y A y y x z A z A z z x y
                 
Selanjutnya  membagi  (17)  dengan ∆𝑡 ,  dan  mengambil  limit ∆𝑡 → 0 pada  ruas  kanan, maka diperoleh
      0 ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) lim t u x u x x y z v y v y y x z m t t m t m w z w z z x y t t
       
                         (18) Karena p m V   , jika p V V   adalah volume pori maka m x y z      . Selanjutnya membagi (18) dengan x y z    diperoleh
  

  ( ) ( )
( ) ( )
( ) ( ) u x u x x v y v y y x y w z w z z z t
 
 
                    (19) Ambil limit 0, 0, 0 x y z
      pada ruas kiri (19), maka diperoleh
Δ x Δ y Δ z 𝜌(𝑥 + ∆𝑥)𝑢(𝑥 + ∆𝑥)𝐴(𝑥 + ∆𝑥) 𝜌(𝑥)𝑢(𝑥)𝐴(𝑥) 𝜌(𝑧)𝑤(𝑧)𝐴(𝑧) 𝜌(𝑧 + ∆𝑧)𝑤(𝑧 + ∆𝑧)𝐴(𝑧 + ∆𝑧) 𝜌(𝑦)𝑣(𝑦)𝐴(𝑦) 𝜌(𝑦 + ∆𝑦)𝑣(𝑦 + ∆𝑦)𝐴(𝑦 + ∆𝑦)
( ) ( ) ( ) ( ) u v w x y z t                 (20)
Dengan  adalah  massa  jenis  fluida, u adalah  laju  aliran  fluida  arah  longitudinal, 𝑣 adalah  laju  aliran  fluida  arah  lateral, w adalah  laju  aliran  fluida  arah  vertikal,  dan  adalah  porositas.  Dengan  menggunakan  turunan  parsial  pada  (20)  dan  karena  dilihat dalam arah longitudinal maka diperoleh
u u x y t t                         (21)
Perubahan  porositas  terhadap  waktu  didapatkan  dengan  menggunakan  hukum kekekalan massa yaitu perubahan rata-rata massa dalam kontrol volume ditambah dengan arus massa yang melalui kontrol volume sama dengan banyak massa yang diciptakan dari sumbernya atau hilang. Menurut teorema pengangkutan Reynold
Rate of change of mass in cell + net outward mass flux = 0
( ) u. 0
faces d A dt  
   
(22)
Dengan menggunakan aliran massa, maka (22) dapat dijabarkan menjadi
( ) ( )
( ) 0 u v
u x u y z v y v t x t y w x z w z w x y t z                                                                         (
23)
Persamaan  (23)  dibagi  dengan x y z    dan  diasumsikan  aliran incompressible maka  konstan  dan  ditambahkan  dengan  komponen  transfer  massa  dari  hasil  reaksi kimia sehingga diperoleh
0 u v w t x y z 
  
                     (24) Dengan ( ) eq kA c c    .  Selanjutnya  menggunakan  turunan  parsial  dan  diambil  dalam arah longitudinal diperoleh
( ) eq kA c c u t x    
             (25)
Setelah  diperoleh  model  pada  (14),  (16),  (21)  dan  (25)  dilanjutkan  dengan  membentuk persamaan menjadi persamaan non dimensional untuk memudahkan dalam simulasi tanpa bergantung  pada  satuan  dan  dimensi.  Selanjutnya  penyelesaian  secara  numerik
menggunakan  metode  beda  hingga  skema  MacCormack.  Dalam  hal  ini  diberikan  hasil untuk tahap prediktor dan korektor dari (21) sebagai berikut.
Tahap prediktor
1 1
1 1 n n n n
n n n n n n n n i i i i i i i i i i i i u u u t t x x                                      (26)
Tahap korektor
1 1 1 1
1 1 1 1 1 1 2 1 2 n n n n n n n n n n n n n i i i i i i i i i i i i i n i u u t u t x x                     
                                                           ( 27)
Berikut ini merupakan hasil simulasi numerik menggunakan metode beda hingga skema  MacCormack  untuk  laju  reaksi  kimia,  massa  padatan  gambut  dan    kadar  air gambut yang ditunjukkan pada masing-masing gambar 3, gambar 4 dan gambar 5 berikut.

## Gambar 3. Laju Reaksi
Pada  plot  gambar  3  diketahui  bahwa  konsentrasi  laju  reaksi  pada  reaktan mengalami  penurunan  setiap  waktu.  Penurunan  tersebut  berlangsung  hingga  mendekati nol yang artinya pada saat waktu tertentu proses reaksi kimia akan berhenti.

## Gambar 4. Massa Padatan Gambut
Selanjutnya  pada  plot  gambar  4  diketahui  bahwa  terjadi  peningkatan  massa padatan  gambut  akibat  dari  proses  rekasi  kimia.  Selain  itu,  diketahui  bahwa  luas penampang  yang  lebih  besar  pada  media  gambut  memiliki  nilai  massa  padatan  lebih tinggi dibandingkan dengan luas penampang lebih kecil.

## Gambar 5. Kadar Air Gambut
Pada  plot  gambar  5  diketahui  bahwa  selama  proses  reaksi  terjadi  penurunan  kadar  air pada gambut. Hal ini terjadi akibat air yang diserap oleh gambut selama proses stabilisasi berlangsung.

## SIMPULAN DAN SARAN
Berdasarkan  hasil  dan  pembahasan  disimpulkan  bahwa  konsentrasi  laju  reaksi pada  reaktan  mengalami  penurunan  setiap  waktu,  adanya  peningkatan  massa  padatan
gambut stabilisasi akibat dari proses reaksi kimia, dan penurunan  kadar air pada gambut stabilisasi.  Hal  ini  terjadi  akibat  air  yang  diserap  oleh  gambut  selama  proses  stabilisasi berlangsung.
Saran  untuk  penelitian  lebih  lanjut  dapat  menggunakan  sifat  fisik  gambut  yang lain  misalnya  angka  pori  gambut  dan  diselesaikan  dengan  menggunakan  metode  lain untuk menyelesaikan permasalahan yang lebih kompleks lagi.

## DAFTAR PUSTAKA
Ma‟ruf,  MA.,  Syauqiah,  I,  2015, Efek  Waktu  Pemeraman  Terhadap  Karakteristik  Sifat Fisiktanah  Gambut  Dengan  Bahan  Stabilisasi  Serbuk  Kayu ,  Prosiding  Seminar Nasional Teknik Sipil Unlam, ISBN 978-602-648-300-3, pg 333-344.
Ma‟ruf,  MA.,  Yulianto,  FE,  2016, Tanah  Gambut  Berserat  :  Solusi  dan Permasalahannya  Dalam  Pembangunan  Infrastruktur  Yang  Berwawasan Lingkungan ,  Prosiding  Seminar  Nasional  Geoteknik,  ISBN  978-602-6483-02-7, pg 279-292.
Natakusumah,  DK,  2007, Pemodelan  Hubungan  Hujan  dan  Aliran  Permukaan  Pada Suatu  DAS  dengan  Metoda  Beda  Hingga ,  Prosiding  ITB  Sains  dan  Teknologi, Vo.39 No.1&2, pg 97-123.
Pritchard, PJ., Leylegian, JC, 2011, Introduction to Fluid Mechanic , USA, John wiley & Sons.Inc.
Pradana, M. S., & Widodo, B. (2014). Model Pertumbuhan Kristal Pada Gambut (XVII).
Surabaya: Konfrensi Nasional Matematika XVII. Retrieved from http://knm17.its.ac.id/
Syahruddin, MH, 2014, Persamaan Aliran Air Dalam Media Berpori Sebagai Aliran Air Tanah  (GroundWater) ,  Simposium  Fisika  Nasional  SFN  XXVII,  ISSN  1411- 4771.
Yulianto,  FE,  Mochtar,  NE,  2012, Behavior  of  Fibrous  Peat  Soil  Stabilized  With  Rice Husk  Ash  (RHA)  and  Lime ,  Proceeding  of  8 th International  Symposium  on Lowland Technology, ISBN 978-602-95227-1-6, pg 106-110.
Yulianto,  FE.,  Mochtar,  NE,  2016, The  Effect  Of  Water  Pore  Condition  To  Shear Strength  And  Compression  Behavior  Of  Stabilized  Fibrous  Peat  Mixing  Lime Caco3  And  Rice  Husk  Ash ,  International  Journal  of  Applied  Engineering Research, ISSN 0973-4562, Vol.11 No.15, pg 8578-8582.