Jurnal Matematika, Sains, dan Pembelajarannya

## ISSN 2622-8904 (online)
Perbandingan Suction Potential dan Water Content pada Saluran
Irigasi Menggunakan Dual Reciprocity
Jahring* 1 , Nasruddin 2 1,2 Universitas Sembilanbelas November Kolaka e-mail: 1 jahring.usn@gmail.com , 2 nash.math@gmail.com

## Abstrak
Penelitian  membahas    penerapan  metode  elemen  batas  dual-reciprocity  (MEBDR)  pada  saluran irigasi  jalur  yang  periodik  pada  tanah  homogen.  MEBDR  merupakan  salah  metode  numerik  yang  dapat menyelesaikan  masalah  infiltrasi  stasioner  pada  saluran  periodik  pada  tanah  yang  homogen.  Beberapa penelitian  sebelumnya  hanya  menhitung  secara  analisis.  Oleh  karena  itu,  MEBDR  ini  sangat  diperlukan. Penelitian  ini  meruapakn  penelitian  studi  literatur  dengan  menggunakan  simulasi  MATLAB.  Dalam peneyelesaiannya dibangun dari Persamaan Richards. Persamaan Richards ini kemudian ditransformasikan menjadi Persamaan Helmholtz termodifikasi. Selanjutnya dengan DRBEM, solusi numerik dari Persamaan Helmholtz termodifikasi diperoleh. Menggunakan solusi numerik yang diperoleh, dapat dihitung nilai numerik dari suction potential dan water content. Dari beberapa saluran periodik, diperoleh saluran non-flat tanpa lapisan memiliki infiltrasi yang baik dibanding saluran flat dan non flat dengan lapisan impermeabel. Dalam hal ini saluran rectangular memiliki suction potential dan water content paling baik jika dibandingkan dengan trapezoidal  dan  semi  circular.  Karena  perbedaan dari  ketiga  saluran  tersebut  tidak  begitu  signifikan.  Jika diterapkan dalam bidang pertanian, maka akan membantu dalam meningkatkan hasil pertanian, khususnya tanaman yang menggunakan irigasi jalur.
Kata kunci: Infiltrasi, saluran periodik, metode elemen batas, dual reciprocity, irigasi jalur, suction potential, water content

## 1. PENDAHULUAN
Air  merupakan  unsur  yang  sangat  diperlukan  oleh  tanaman  dan  segala  hal  yang  berkaitan  dengan kegiatan pertanian. Dengan ketersediaan air yang cukup dapat menjadi faktor pendukung meningkatnya hasil pertanian. Namun tidak dapat dipungkiri, tidak setiap daerah atau lahan pertanian memiliki ketersediaan air yang  cukup. Hal  itu  disebabkan oleh  curah  hujan  yang  rendah  sehingga  akan  sangat  sulit  untuk  memenuhi kebutuhan air tanaman. Berkaitan dengan hal tersebut di atas, sistem irigasi dipandang perlu untuk mengatasi masalah  lahan  pertanian  dengan  intensitas  air  yang  kurang.  Adanya  sistem  irigasi  diharapkan  dapat mempertahankan tingkat kadar air yang cukup tinggi dalam tanah sehingga akar tanaman dapat mengekstrak air dari tanah secara efisien.
Penelitian ini merupakan penelitian berkelanjutan yang dilakukan oleh peneliti dengan menambahkan beberapa saluran irigasi jalur disetiap penelitian. Penelitian Jahring (2016) tentang penerapan metode elemen batas dual-reciprocity pada saluran irigasi jalur dengan saluran yang dilengkapi dengan lapisan impermeabel, yaitu trapezoidal dan rectangular dengan  lapisan  impermeabel,  kemudian  penelitian  Jahring  (2018) dilanjutkan  dengan  menambahkan  saluran trapezoidal dan rectangular tanpa  lapisan  impermeabel.  Namun pada  kenyataannya  bahwa  saluran  yang  ada  di  lapangan  tidak  hanya  saluran trapezoidal dan rectangular dengan dan atau tanpa lapisan impermeabel, tetapi juga terdapat saluran irigasi jalur tanpa lapisan impermeabel yang lain yaitu saluran semi circular dan sebagai pembanding ditambahkan suatu saluran flat . Oleh karena itu, dalam penelitian ini akan mengembangkan penelitian sebelumnya dengan menambahkan saluran irigasi jalur
https://jurnal.unsulbar.ac.id/index.php/saintifik
non-flat semi circular dan saluran flat yang dirangkum dalam satu judul “Penerapan Metode Elemen Batas Dual Reciprocity (MEBDR) pada Saluran Irigasi Jalur pada tanah homogen”.

## 2. METODE PENELITIAN

## 2.1. Formulasi Masalah
Berikut  akan  diberikan  formulasi  masalah  infiltrasi  dari  berbagai  saluran-saluran  irigasi  jalur  yang  berbeda pada tanah yang homogen, yaitu:
a. Saluran flat b. Saluran semi circular c. Saluran r ectangular dengan dan tanpa lapisan impermeabel d. Saluran trapezoidal dengan dan tanpa lapisan impermeabel.
Selanjutnya akan ditentukan domain dari masing-masing saluran irigasi di atas. Dalam menentukan domain saluran irigasi, perlu diberikan asumsi-asumsi tertentu yaitu:
a. Panjang penampang masing-masing saluran irigasi adalah sama yaitu 2L b. Saluran sangat panjang
c. Jarak antar saluran adalah sama yaitu 2D , atau jarak antar titik tengah saluran adalah sama yaitu 2(L+D) d. Saluran selalu dipenuhi air
e. Infiltrasi pada permukaan saluran irigasi konstan yaitu v 0
f. Tidak ada infiltrasi pada permukaan tanah di luar saluran irigasi
g. Tidak ada infiltrasi pada permukaan saluran impermeabel.
Gambar 1. Saluran Irigasi Jalur (Solekhuddin, 2013 dan Lobo, 2008)
https://jurnal.unsulbar.ac.id/index.php/saintifik
Dengan  sifat  simetris,  dapat  didefinisikan  domain  yaitu
D L X +   0 dan 0  Z yang  dinyatakan dengan R dengan batas C , dan tidak ada fluks pada 0 = X dan D L X + = .

## 2.2. Persamaan Fundamental
Menurut  Azis  (2003),  untuk  masalah  infiltrasi  stasioner,  persamaan  yang  sering  digunakan  adalah persamaan  Richards,  yang  merupakan  persamaan  differensial  parsial  non  linear  yang  merepresentasikan perpindahan air dalam tanah tak jenuh dalam dua dimensi adalah:
Z K Z K Z X K X   =           +             (1)
dengan K adalah konduktivitas hidrolik tanah tak jenuh dan Z > 0 mengarah ke bawah secara vertical. Dalam makalah  yang  sama  menjelaskan  potensial  fluks  matrik  Θ  dihubungkan dengan  konduktivitas  hidrolik  oleh Persamaan
( )   − =   dq q K (2)
Hubungan antara konduktivitas hidrolik K dari tanah tak jenuh ( unsaturated soil ) dan konduktivitas hidrolik K 0 dari tanah jenuh ( saturated soil ) didefinisikan sebagai
0 , 0  =   e K K (3)
dengan ψ adalah suction potential dan α adalah suatu konstanta empiris.
Menggunakan Persamaan (2) dan (3) diperoleh bentuk linear dari Persamaan (1) yaitu
Z Z X    =    +     2 2 2 2 (4)
Komponen fluks ke arah horizontal dan vertikal adalah
Z V X U    −  =    − =  dan , (5)
Fluks  normal  pada  permukaan  tanah  dengan  vektor  normal  terhadap  batas  domain  yang  mengarah keluar domain dengan n = ( n 1 , n 2 ) diberikan oleh
2 1 n Z n X F          −  +    − =  (6)
Didefinisikan peubah-peubah tak berdimensi sebagai berikut
F L v f V L v v U L v u L v Z z X x          0 0 0 0 2 ; 2 ; 2 ; ; 2 ; 2 = = =  =  = = (7)
dan mensubstitusikan z e  =  ke Persamaan (4) diperoleh suatu Persamaan berbentuk Persamaan Helmholtz termodifikasi, yaitu
   =   +   2 2 2 2 z x (8)
Fluks tak berdimensi diperoleh
z fe n n − − =   2   (9)
Fluks pada permukaan saluran adalah v 0 dan besar fluks tak berdimensi adalah 2 𝜋 / 𝛼 L . Sehingga syarat batas pada permukaan saluran permeabel adalah
2 2 . z e n n L     −  = +  (10)
Metode elemen batas Dual Reciprocity (MEBDR) digunakan untuk mencari solusi numerik Persamaan (8).  Pada  tahap  selanjutnya,  mencari  solusi  numerik  dengan  MEBDR  dibutuhkan  suatu  Persamaan  integral yang disebut Persamaan integral batas ( Boundary integral equation ) yaitu
https://jurnal.unsulbar.ac.id/index.php/saintifik
( ) ( ) ( ) ( ) ( ) ( ) ( )
, ; , , , , , , ; , , C x z x z x z x z ds x z n n                 = −         ( ) ( ) , ; , , ,
R x z x z dxdz    +   (11)
dengan solusi fundamental persamaan Laplace dua dimensi
( ) ( ) ( ) 2 2 1 , ; , ln 2 x z x z         = − + −   (12)
dan
( ) ( ) ( ) ( ) 0, Jika , 1 , , Jika , berada pada bagian kurva 2
1, Jika , . R C smooth C R              =    
(13)
Persamaan  integral  (11)  dapat  diselesaikan  secara  numerik  dengan  mendiskritisasi  batas  domain menggunakan ruas-ruas garis yang berhubungan satu sama lain dengan ujung-ujungnya terletak pada kurva C dan mendiskritisasi domain dengan titik kolokasi interior. Misalkan ruas-ruas garis tersebut C (k) , k = 1,2,...,N sehingga jumlah ruas garis sebanyak N ruas garis. Pada setiap C (k) , k = 1,2,...,N , dipilih titik tengahnya sebagai titik kolokasi ( collocation points ) yang dituliskan dengan (a (1) ,b (1) ), (a (2) ,b (2) ),..., (a (N) ,b (N) ) dan titik-titk kolokasi interior ( interior points collocation ) sebanyak L titik yang dapat dituliskan dengan (a (N+1) ,b (N+1) ),...,(a (N+L) ,b (N+L) ) .  Misalkan  ( ξ,η)  =  (x (n) ,y (n) ) untuk n=1,2,...,N+L ,  maka  aproksimasi  dapat dituliskan sebagai
( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) 2 1 1 1 , , , n N L N n n nj k k n n k k n n j k x y F x y p F x y     + = = = + −   (14)
dengan ( ) ( ) ( ) ( ) ( ) , 1, 2,..., n n n x y n N L  = = + dan ( ) ( ) ( ) ( ) ( ) ( ) ( ) 1 , ; , , n n m N L nj mj m m x y x y   + = =   (15)
( ) ( ) ( ) ( ) ( ) 1 , , ; , , , k k C F x z ds x z     =  
(16)
( ) ( ) ( ) ( ) ( ) 2 , ; , , , . k k C x z F ds x z n        =      
(17)
dan
( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) 1 , ,
, ; , , ; ,
, ; , , , ; ,
k k N k k x y k x y x y x y a b
x y a b a b x y a b n n x y         = =    
 = + +        (18)
( ) ( ) ( ) ( ) ( ) ( ) ( ) 1 2 1 , , ; , , . N k k k k k F x y a b F      = −  (19)
( ) ( ) ( ) ( ) ( ) ( ) 1 , ; , mj j j m m x y x y   −   =   (20)
dengan
( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) 2 3 2 4 5 2 2 , ; , 1 , ; , , ; , ,
1
1 1 , ; , , ; , , ; , , ; , , 4 16 25
, ; , .
m m m m m m m m m m m m m m m m m m x z a b r x z a b r x z a b
x z a b r x z a b r x z a b r x z a b r x z a b x a z b
  = + +
= + + = − + − (21)
MEBDR  digunakan  untuk  memperoleh  solusi  numerik (ϕ) ,  dan  menggunakan  Persamaan  berikut untuk memperoleh suction potential (ψ) , (Jahring, 2016).
0 0 1 ln . z v L e K        =     (22)
https://jurnal.unsulbar.ac.id/index.php/saintifik
Untuk memperoleh water content ( 𝜃 ) digunakan Persamaan (Warrick, 2001)
( ) 1 1 m r n
s r        − =   − +     (23)
dengan θ r dan θ s adalah water  content residual  dan  jenuh, α , n dan m  =  1  –  (1/n) adalah  parameter  yang bergantung pada jenis tanah.

## 3. HASIL DAN PEMBAHASAN
MEBDR diterapkan pada masalah infiltrasi stasioner pada beberapa saluran periodik yang berbeda- beda pada tanah homogen dan dipilih jenis tanah Pima clay loam (PCL) dengan nilai α dan K 0 masing-masing 0.014 cm -1 dan 9.9 cm/hari (Jahring, 2016), (Lobo, 2008), (Solekhuddin, 2013), θ r = 0.095 cm -1 , θ s = 0.41 cm - 1 dan n=1,31 cm -1 (Warrick, 2001), dan nilai v 0 adalah 75% dari nilai K 0 (Basha, 1999). Selanjutnya ukuran dari  tiap-tiap  saluran non-flat  impermeabel dapat  dipilih  sedemikian  sehingga  panjang  permukaan  saluran sama yaitu L = 50 cm dan D = 50 cm. Menggunakan nilai L tersebut, maka lebar saluran non-flat adalah 200/ π cm dengan kedalaman (50 – 100/π) cm untuk saluran rectangular , 75/π cm untuk saluran trapezoidal , serta masing-masing 50 cm dan 49.5625 cm untuk saluran rectangular dan trapezoidal dengan lapisan impermeabel. Selanjutnya MEBDR diimplementasikan dalam bahasa pemrograman MATLAB.
Pada  MEBDR banyaknya  segmen garis  ( N )  dan  titik  interior  ( M )  berpengaruh  terhadap keakuratan solusi numerik yang dihasilkan. Berikut diberikan jumlah segmen garis dan titik interior, (Jahring, 2016).
1. Saluran flat diberikan N = 200 dan M = 625 2. Saluran semi circular diberikan N = 201 dan M = 620
3. Saluran rectangular diberikan N = 203 dan M = 625 4. Saluran trapezoidal diberikan N = 202 dan M = 619 5. Saluran rectangular impermeabel diberikan N = 208 dan M = 609 6. Saluran trapezoidal impermeabel diberikan N = 209 dan M = 609 Selanjutnya  variasi  nilai ψ dan θ dengan  kedalaman (Z) pada  jarak  yang  berbeda-beda  dari  pusat saluran (X) secara geometri dapat dilihat pada Gambar 2-6 berikut ini.
Gambar 2. Suction potential dan water content pada X = 15
Dari Gambar 2(a) yaitu pada X = 15 yang terletak di bawah saluran, nilai ψ untuk semua jenis saluran (kecuali saluran dengan lapisan impermeabel) menurun  terhadap  Z.  Hal ini  menunjukkan bahwa  kandungan air  ( water  content )  di bawah saluran pada bagian dangkal lebih besar dibandingkan pada bagian terdalam di dalam  tanah  sehingga  perembesan  air  semakin  ke  dalam  semakin  kecil.  Hal ini  didukung dengan nilai θ yang  semakin  menurun  (lihat  Gambar  2(b))  yang  menunjukkan bahwa  kandungan  air  semakin  ke  bawah semakin kecil. Dari Gambar juga terlihat bahwa khusus untuk saluran dengan lapisan impermeabel, nilai ψ dan
https://jurnal.unsulbar.ac.id/index.php/saintifik
θ meningkat. Ini disebabkan karena pada posisi X  = 15 terletak di bawah saluran yang impermeabel sehingga berdasarkan asumsi awal tidak ada infiltrasi.
Gambar 3. Suction potential dan water content pada X = 30
Gambar 3 di atas dengan X = 30 juga menunjukkan penurunan nilai ψ dan θ terhadap Z karena masih terletak di bawah saluran. Namun ada sedikit perbedaan pada saluran dengan lapisan impermeabel. Terlihat bahwa ψ dan θ meningkat. Hal ini disebabkan karena adanya lapisan impermeabel pada saluran. Pada beberapa saluran non-flat tanpa lapisan impermeabel, tidak ada perbedaan nilai ψ yang signifikan. Namun jika dibandingkan dengan saluran flat dan non flat dengan lapisan impermeabel, terdapat perbedaan nilai ψ yang signifikan khususnya di bawah saluran pada bagian dangkal. Ini disebabkan oleh posisi permukaan saluran.  Pada  saluran flat permukaan  terletak  pada  permukaan  tanah  yaitu  pada Z  =  0 ,  sedangkan  kelima saluran yang lain permukaannya terletak dari Z = 0 sampai kedalaman yang ditentukan.
Gambar 4. Suction potential dan water content pada X = 50
Saluran flat dengan X = 50 pada Gambar 4 terletak pada ujung saluran dengan nilai ψ dan θ yang secara grafik  konstan pada  nilai ψ = –70 dan θ = 36,3 .  Ini dikarenakan  pada X  =  50 terletak  antara  saluran  yang sebelumnya diasumsikan terdapat infiltrasi  dan permukaan tanah yang sudah diasumsikan tidak ada infiltrasi sehingga nilai ψ hampir konstan. Selain itu, kondisi ini menunjukkan bahwa sepanjang titik X = 50 , infiltrasi sama yang menyebabkan kandungan air sama.
Pada saluran dengan lapisan impermeabel (Gambar 4(a), 5(a), dan 6(a)) nilai ψ naik pada kedalaman tertentu  dan  turun pada kedalaman  tertentu  juga  menuju  ke konvergenannya.  Ini menunjukkan  bahwa  pada bagian dangkal sampai kedalaman tertentu di bawah permukaan tanah, kandungan air lebih rendah dari bagian yang  lebih  dalam.  Hal  ini  disebabkan  oleh  panjang  saluran  yang  dialiri  air  cukup  panjang  sehingga menyebabkan adanya aliran air horizontal yang sangat besar pada saluran yang mempengaruhi kandungan air di bawah permukaan tanah di luar saluran. Sehingga grafik nilai ψ dan θ naik dari permukaan tanah sampai pada  kedalaman  yang  dipengaruhi  aliran  air  horizontal.  Selanjutnya  grafik  akan  turun  sampai  menuju  ke
https://jurnal.unsulbar.ac.id/index.php/saintifik
konvergenannya yang dipengaruhi oleh infiltrasi yang ada di sekitar saluran. Hal ini diperjelas pada Gambar 4(b), 5(b), dan 6(b) dengan nilai θ yang naik turun hingga menuju kekonvergenannya.
Gambar 5. Suction potential dan water content pada X = 65
Gambar 6. Suction potential dan water content pada X = 85
Pada X  =  65 dan X  =  85 ,  nilai  ψ  meningkat  terhadap Z (Gambar  5(a)  dan  6(a).  H a l   ini menunjukkan bahwa kandungan air ( water  content ) pada permukaan tanah paling rendah dibandingkan pada semua bagian tanah di dalam domain. hal ini juga ditunjukkan pada Gambar 5(b) dan 6(b) dengan nilai θ yang semakin  meningkat.  Ini  karena  adanya  asumsi awal  bahwa  pada  permukaan tanah  di  luar saluran tidak ada aliran air yang masuk ke dalam tanah.  Pada beberapa saluran non-flat tanpa  lapisan  impermeabel, tidak  ada perbedaan nilai ψ dan θ yang signifikan.  Namun jika  dibandingkan dengan saluran flat dan non-flat dengan lapisan impermeabel, terdapat perbedaan nilai ψ dan θ di luar saluran.
Dari Gambar 2 –.6 dapat dijelaskan bahwa hubungan antara suction potential ( ψ ) dan water  content ( θ )  pada  jenis  tanah Pima  Clay  Loam (PCL)  adalah  berbanding  lurus.  Dengan  kata  lain,  semakin  besar infiltrasi maka  semakin besar pula kandungan airnya. Perhatikan saluran flat pada Gambar 2(a) – 6(a).  Dari hasil yang diberikan, saluran flat memiliki infiltrasi dan kandungan air paling besar.  Namun pada realisasinya, saluran flat tidak mungkin untuk diterapkan dengan asumsi-asumsi yang telah diberikan sebelumnya.  Hal ini dikarenakan sejatinya saluran flat tidak ada, hanya bisa diselesaikan secara analitik saja tanpa bisa diterapkan.
https://jurnal.unsulbar.ac.id/index.php/saintifik
Selanjutnya  perhatikan  pula  saluran non-flat tanpa  lapisan  impermeabel  ( semi  circular , rectangular dan trapezoidal )  pada  Gambar 2(a)  – 6(a).  Pada  Gambar  tersebut  terlihat  bahwa  saluran rectangular memiliki infiltrasi paling baik jika dibandingkan dengan semi circular dan trapezoidal .  Hal ini dikarenakan nilai suction potential di  bawah saluran pada bagian dangkal paling besar dan nilai suction potential di bawah  permukaan tanah di  luar saluran pada bagian dangkal paling kecil.  Namun karena perbedaan nilai suction potential dari ketiga saluran tersebut tidak begitu signifikan, sehingga hanya diberikan distribusi nilai ψ untuk satu saluran saja.
Jika  dipandang  dari  segi  aplikatif,  di  Indonesia  banyak  dijumpai  saluran-saluran  dengan  bentuk trapezoidal baik itu diperoleh dari internet ataupun dari sumber informasi yang lain. Hal ini disebabkan oleh selain cara membuat salurannya sangat mudah, juga memiliki nilai suction potential dan water content paling baik  diantara  saluran  yang  lain. Oleh  karena  itu,  untuk  saluran  periodik non-flat dengan dan  tanpa  lapisan impermeabel hanya akan diberikan distribusi nilai ψ pada saluran trapezoidal ,  dengan anggapan bahwa secara umum infiltrasi dari setiap saluran dengan saluran yang lain dapat dianggap sama.  Secara geometri distribusi nilai ψ disajikan pada Gambar 7 berikut ini.
Gambar 7. Surface plot  distribusi nilai ψ pada saluran trapezoidal
Dari  Gambar 7(a)  di  atas,  terlihat  jelas  distribusi nilai ψ .  Pada bagian  dangkal  di  bawah  permukaan saluran  memiliki  nilai ψ paling  besar  di  antara  bagian  yang  lain  dan  terus  menurun  menuju  ke  titik kekonvergenannya.  Ini  menunjukkan  bahwa  kandungan  air  ( water  content )  di  bawah  permukaan  saluran paling  besar  dibandingkan dengan bagian  lain.  Pada  bagian di  bawah  permukaan tanah di luar saluran yang tidak  jauh  dari  saluran  terlihat  semua  nilai ψ hampir  sama.  Ini  menunjukkan  bahwa  kandungan  air  ( water content )  hampir  sama  pada semua  bagian  kedalaman tanah.  Selanjutnya pada  bagian  di  bawah  permukaan tanah  di  luar  saluran  yang  merupakan  titik  terjauh  dari  pusat  saluran  dan  di  bawah  permukaan  saluran impermeabel memiliki nilai ψ yang paling rendah dibandingkan dengan bagian yang lain dan terus meningkat menuju titik kekonvergenannya pada kedalaman tertentu.  Hal ini menunjukkan bahwa di  bawah permukaan tanah di luar saluran pada posisi jauh dari pusat saluran dan di bawah permukaan saluran impermeabel memiliki kandungan air ( water  content ) paling rendah dibandingkan dengan bagian tanah yang lain dalam domain.
Dari Gambar 7(b) di atas terlihat jelas distribusi nilai ψ .  Pada bagian dangkal  di  bawah  permukaan saluran  memiliki  nilai ψ paling  besar  di  antara  bagian  yang  lain  dan  terus  menurun  menuju  ke  titik
https://jurnal.unsulbar.ac.id/index.php/saintifik
103
kekonvergenannya. Ini menunjukkan bahwa kandungan air ( water content ) di bawah permukaan saluran paling besar dibandingkan dengan  bagian  lain.  Pada  bagian  di  bawah  permukaan tanah di  luar saluran yang  tidak jauh dari saluran terlihat semua nilai ψ hampir sama. Ini menunjukkan bahwa kandungan air ( water  content ) hampir  sama  pada  semua  bagian kedalaman tanah.  Selanjutnya pada  bagian di  bawah  permukaan tanah  di luar saluran yang merupakan titik terjauh dari pusat saluran memiliki nilai ψ yang paling rendah dibandingkan dengan bagian yang lain dan terus meningkat menuju titik ke konvergenannya pada  kedalaman tertentu. Hal ini menunjukkan bahwa di bawah permukaan tanah di luar saluran pada posisi jauh dari pusat saluran memiliki kandungan air ( water  content ) paling rendah dibandingkan dengan bagian tanah yang lain dalam domain.

## 4. KESIMPULAN
Berdasarkan hasil dan pembahasan, maka dapat disimpulkan  bahwa Suction potential pada semua jenis saluran paling besar terletak pada domain di bawah permukaan saluran pada bagian dangkal. Ini menunjukkan bahwa water content pada bagian dangkal di bawah permukaan saluran lebih tinggi dari bagian terdalam pada domain di dalam tanah. Suction potential pada semua jenis saluran paling kecil terletak pada domain di bawah permukaan  tanah  di  luar  saluran  pada  bagian  dangkal.  Ini  menunjukkan  bahwa water  content pada  bagian dangkal di bawah permukaan tanah di luar saluran paling rendah dari pada bagian domain yang lain di dalam tanah.  Terdapat  perbedaan  nilai ψ yang  signifikan  antara  saluran flat dengan  saluran non  flat baik  tanpa maupun dengan lapisan impermeabel. Pada jenis tanah pima clay loam (PCL), suction potential berbanding lurus terhadap water content . Artinya, semakin besar infiltrasi maka semakin besar pula kandungan airnya.

## DAFTAR PUSTAKA
Azis, M.I., Clements, D.L., Lobo, M., 2003, A Boundary Element Method for Steady Infiltration from Periodic Channels , ANZIAM J. 44:C61-C78.
Basha, H.A., 1999, Multidimensional Linearized Nonsteady Infiltration with Prescribed Boundary Conditions at the Soil Surface , Water Resources Research, 35(1):75-83.
Batu, V., 1977, Steady Infiltration from a Ditch ; Theory and Experiment, Soil Sci, Soc, Am, j., 41:677-681.
Jahring,  2016,  The  Application  of  Dual  Reciprocity  Boundary  Element  Method  for  Steady  Infiltration Problems in the Periodic Channel with Impermeable Layer, Journal of Math Sciences (JMS). 1.2;28-34.
Lobo,  M.,  2008, Boundary  Element  Methods  for  A  Class  of  Infiltraton  Problems ,  PhD  Tesis,  Adelaide University.
Solekhudin,  I.,  Ang,  K.C.,  2013, A  Dual-Reciprocity  Boundary  Element  Method  for  Steady  Infiltration Problems , ANZIAM J. 54:171-180.
Warrick, A.W., 2001, Soil Physics Companion , CRC Press LLC, Boca Raton.