Bagaimana cara mengecek agar hasil penggambaran grafis yang kita buat telah mendekati solusi persamaan kontinuitas? Di tulisan ini saya akan mendemonstrasikan bagaimana menggunakan Excel untuk membuat garis-garis ekuipotensial berdasarkan persamaan kontinuitas!!
Me-review sedikit, pada tulisan yang lalu saya telah membahas mengenai formulasi beda hingga untuk mendapatkan nilai aproksimasi dari persamaan kontinuitas dari tinggi potensial hidrolik, yang mana persamaan tersebut didasarkan pada Hukum Darcy. Di tulisan kali ini, saya akan mencoba menyelesaikan variasi problem rembesan/aliran air pada struktur geoteknik menggunakan pendekatan beda hingga.
Rangkuman rumus beda hingga
Disini saya berikan rangkuman rumus dari formulasi beda hingga yang telah saya jabarkan di posting yang lalu
Nodal interior:
Nodal boundary tak terdefinisi (asumsi no-flow):
(1a). Vertikal kiri:
(1b). Vertikal kanan:
(2a). Horizontal bawah:
(2b). Horizontal atas:
(3a). Sudut kiri bawah:
(3b). Sudut kanan bawah:
(3c). Sudut kanan atas:
(3d). Sudut kiri atas:
Tinggi hidrolik total (hydraulic head/piezometric head)
Agar bisa menghitung kondisi batas pada problem rembesan (seepage), maka kita harus bisa menghitung besarnya tinggi hidrolik 
Bila kita ingat persamaan Bernoulli:
Penjumlahan ketiga terms ini didefinisikan sebagai tinggi hidrolik total. Terms pertama menyatakan velocity head. Ingat bahwa Hukum Darcy berlaku pada aliran Stokes (laminer), sehingga terms ini dapat diabaikan. Sedangkan penjumlahan terms kedua (pressure head) dan ketiga (elevation head) menghasilkan apa yang dikenal sebagai hydraulic head atau piezometric head.
Tinggi hidrolik total atau tinggi tekanan piezometrik sesungguhnya merepresentasikan kontribusi tekanan statik yang diukur pada titik lokal tertentu. Misalkan kita memasang piezometrik pada beberapa 2 posisi tertentu didalam tanah seperti pada gambar dibawah ini:
Tinggi tekanan piezometric
Pada piezometrik pertama (dikiri) pada gambar diatas, diukur tinggi tekanan statik sebesar 8 m, sedangkan pada piezometrik satunya, diukur tinggi tekanan statik sebesar 7 m. Ini adalah tinggi tekanan yang dinamai pressure head.
Sedangkan elevasi dimana dilakukan pengukuran pada kedua piezometrik diatas masing-masing adalah -13 m dan -10 m. Cat: Letak datum bisa dimana saja, asalkan konsisten!!
Sehingga nilai hydraulic head untuk masing-masing piezometrik diatas adalah
Contoh kasus 1: Sheetpile dan ekskavasi
Sebuah pekerjaan konstruksi akan dilakukan ditengah badan aliran, oleh karena itu disepanjang bagian tengah bagian konstruksi tersebut akan ditanam sheetpile. Agar pekerjaan dapat dilakukan, tinggi permukaan air dijaga selalu nol. Untuk lebih jelasnya perhatikan gambar dibawah ini:
Contoh 1: Sheetpile dan ekskavasi
Datum saya ambil di elevasi terendah dari permukaan tanah. Ingat bahwa pemilihan datum tidak akan mempengaruhi hasil akhir!! Jadi silahkan saja pilih datum di elevasi lain. 
Tinggi hidrolik total di bagian kiri pada permukaan pasir:
Pada bagian tengah:
Sedangkan pada bagian kanan:
Menggunakan Excel, dimana setiap sel merepresentasikan panjang atau lebar 0.5 m, saya akan mencoba menyelesaikan problem diatas dengan pendekatan beda hingga.
Untuk mencari solusi beda hingga saya akan memanfaatkan kemampuan Excel untuk menyelesaikan circular reference, sehingga pertama-tama masuk ke File>Options lalu pilih tab Formula, beri tanda centang pada pilihan Iteration untuk mengaktivasi, terakhir jangan lupa masukkan angka yang cukup besar pada Maximum Iterations, misalnya 10000.
Pertama-tama, definisikan seluruh kondisi batas yang ada seperti pada gambar dibawah ini. Warna orange tua di bagian pinggir hanya digunakan untuk membantu saya menghitung koordinat. Warna putih dan warna biru, yang masing-masing merepresentasikan daerah kosong dan air, merupakan bagian diluar sistem problem boundary yang sesungguhnya tidak diperlukan untuk didefinisikan, namun disini saya masukkan nilainya agar hasil tampilan di Excel bisa lebih baik.
Bagian berwarna hijau diambil berdasarkan perhitungan hydraulic head diatas.
Bagian berwarna ungu diisikan formula beda hingga sesuai formula beda hingga untuk nodal boundary.
Kemudian isikan bagian interior (berwarna orange muda) dengan formula beda hingga untuk nodal interior. Excel akan melakukan iterasi, dan kecepatan mencapai konvergensi solusi akan tergantung dari kecepatan performa komputer. Hasilnya adalah sbb:
Bila kita plot hasil kalkulasi diatas dengan Excel dengan menggunakan chart tipe surface, maka kita akan peroleh gambar seperti dibawah ini yang mana garis kontur yang membatasi gradasi warna merupakan garis-garis ekuipotensial!!
Pada gambar diatas saya telah menambahkan garis yang menggambarkan permukaan tanah dan permukaan air, serta garis tebal yang menggambarkan sheetpile. Legenda dibagian bawah gambar menginformasikan tinggi hydraulic head. Garis-garis kontur ekuipotensial sengaja saya buat dengan gradien hidrolik
Sekarang bagaimana dengan garis alirannya?? 😦 Bila kita ingat formula garis aliran yang tergantung dari besar gradien hidrolik, maka kita tahu bahwa garis aliran harus tegak lurus garis ekuipotensial.
Namun penggambaran garis aliran tidak dapat dengan mudah dibuat dengan Excel (bisa tapi tidak mudah/cepat), karena garis aliran adalah vektor sedangkan tinggi hidrolik adalah skalar. Oleh karena itu saya akan sambung cara pembuatannya di next post
Plus saya akan berikan beberapa contoh lainnya juga. To be continued… 😎
James's Files 
Terima kasih, blognya sangat membantu. Kalau boleh bertanya, bagaimana cara menambahkan warna pada bagan surfacenya sekaligus mengubah rentangnya dari 5 menjadi 0.5 ya? Soalnya secara default hanya ada 3 warna saja dengan rentang 5 (0-5, 5-10, dan 10-15). Terima kasih banyak.