J. Math. and Its Appl. ISSN: 1829-605X Vol. 2, No. 2, Nov 2005, 93–101

## Simulasi Model Gelombang Pasang Surut dengan Metode Beda Hingga

## Lukman Hanafi, Danang Indrajaya
Jurusan Matematika FMIPA ITS Kampus ITS Sukolilo Surabaya, 60111 lukman@matematika.its.ac.id

## Abstrak
Beberapa pesisir laut di Indonesia memiliki kondisi pasang surut yang cukup tinggi. Hal ini dapat mempengaruhi pembangunan sektoral dan regional yang dilakukan pemerintah maupun swasta yang berlangsung secara intensif di daerah pesisir. Untuk itu, diperlukan suatu perangkat untuk memperkirakan kondisi pasang surut perairan dalam mendukung perencanaan suatu bangunan laut. Salah satu metode numerik yang digunakan untuk menghitung ketinggian dan kecepatan suatu gelombang pasang surut adalah metode beda hingga. Dalam makalah ini dibahas mengenai pembentukan model matematika gelombang pasang surut dan implementasinya, yaitu perhitungan numerik dengan metode beda hingga dengan simulasi grafiknya.
Kata kunci : beda hingga, gelombang pasang surut, simulasi

## 1. Pendahuluan
Indonesia merupakan negara kepulauan yang dikelilingi oleh lautan. Dua lautan besar yakni Samudera Indonesia dan Samudera Pasifik dan posisinya yang bera-

## Simulasi Model Gelombang Pasang Surut dengan Metode Beda Hingga
da di garis katulistiwa menyebabkan kondisi pasang surut, angin, gelombang, dan arus laut cukup besar. Perubahan pasang surut permukaan air laut dapat mempe- ngaruhi pembangunan sektoral dan regional yang dilakukan pemerintah maupun swasta yang berlangsung secara intensif di daerah pesisir. Mengingat saat ini daerah pantai telah dimanfaatkan dan dikembangkan menjadi daerah pariwisata, transportasi, eksplorasi lepas pantai, industri, pusat perdagangan dan lain-lain. Perubahan pasang surut permukaan air laut juga dapat dimanfaatkan sebagai sum- ber energi. Sumber energi pasang surut bersifat bukan saja dapat diperbaharui tetapi juga tidak menghasilkan polusi sebagaimana energi air. Mengingat sum- berdaya energi hidrokarbon Indonesia akan terus menyurut dan bahkan habis di suatu saat, perlu sejak dini dilakukan eksplorasi dan pengkajian pemanfaatan en- ergi pasang surut mengingat posisi negara Indonesia yang dikelilingi lautan dengan potensi energi pasang surut cukup besar. Untuk itu, diperlukan suatu perangkat untuk membantu proses pengkajian gelombang pasang surut. Dengan simulasi dari model matematika , dapat diprediksikan suatu karakteristik dari gelombang pasang surut.

## 2. Konstruksi Model Matematika Gelombang
Gambar 1 menunjukkan situasi volume kendali dari aliran tidak tunak (non-steady flow) dengan kedalaman air h dan kecepatan rata-rata adalah fungsi pada posisi x dan waktu t. Lebar aliran b adalah fungsi x.
Gambar 1: Volume kendali dari aliran tidak tunak terbuka
Debit Q = R A u dz = b hu Perubahan massa fluida dalam waktu ∆ t berkaitan dengan perubahan ketinggian permukaan fluida adalah:
ρb ( h + δh δt )∆ x − ρbh ∆ x = ρb δh δt ∆ t ∆ x (1)
Perubahan massa fluida dalam waktu ∆ t yang berkaitan dengan nilai keluar masuknya fluida adalah:
ρQ ∆ t − ρ ( Q + δQ δx )∆ t = − ρ δQ δx ∆ t ∆ x (2)
Dengan menyamakan ekspresi Persamaan 1 dan 2 diperoleh persamaan kontinu- itas:
b δh δt + δQ δx = 0 (3)
Persamaan dasar kedalaman rata-rata (depth-average) yang melukiskan fenom- ena gelombang panjang adalah persamaan aliran permukaan bebas (free surface flow) untuk aliran homogen satu dimensi
Gambar 2: Diagram skematik aliran permukaan bebas
Persamaan gerak pada sumbu- x berdasarkan Gambar 2 diatas, dimana gaya- gaya eksternal diabaikan dan ρ adalah konstan [3] adalah :
δQ δt + δ δx µ α Q 2 A ¶ + gA δz s δx + Q ρ τ b − gi b = 0 (4)
Dengan merepresentasikan tekanan alas (bottom-shear stress) sebagai maka persamaan geraknya menjadi:
δQ δt + δ δx µ α Q 2 A ¶ + gA δz s δx + g C 2 A 2 R Q | Q | − gi b = 0 (5)
dimana: Q = Au , A = bh , R = A O dan z s = h + z b b = Lebar aliran (konstan pada waktu).
h = Rata-rata kedalaman air laut
Q = Debit A = Daerah aliran R = Radius hidraulik

## Simulasi Model Gelombang Pasang Surut dengan Metode Beda Hingga
O = Garis keliling dari dasar aliran
C = Koefisien Chezy u = Rata-rata kecepatan aliran. z b = Ketinggian alas rata-rata diatas bidang horizontal. z s = h + z b = Ketinggian permukaan laut diatas bidang horizontal. α = koefisien ( ≈ 1). i b = Kemiringan alas. x = Koordinat horizontal. g = Percepatan gravitasi.
Dengan subtitusi Q = bhu , persamaan kontinuitas (3) menjadi
δh δt + δ δx ( uh ) = 0 (6)
Pada kasus pada saluran yang luas ( b >> h ) dengan lintang penampang segi- empat dan luas yang konstan, diasumsikan R = h . Ketinggian alas rata-rata ( z b ) diatas bidang horizontal sama dengan nol dengan nilai α = 1., z b = 0, dimana A = bh , Q = bhu , z s = z b + h maka persamaan gerak menjadi :
δu δt + u δu δx + g δh δx + g C 2 h u | u | − gi b = 0 (7)

## 3. Persamaan Kontinuitas dan Gerak Gelombang Pasang Surut
Dasar persamaan kontinuitas dan persamaan gerak gelombang pasang surut adalah persamaan gelombang harmonik progresif [3]. Gelombang harmonik pro- gresif adalah sebuah gelombang berkala yang dilukiskan oleh fungsi sinus dan cosinus. Amplitudo gelombangnya adalah relatif kecil pada perairan dalam se- hingga kecepatannya juga kecil. Percepatan konvektif ( uδu/δx ) pada Persamaan 7 dapat diabaikan dan dengan mengasumsikan konstanta friksi, c f = g | u | C 2 h , (dimana i b = 0), maka Persamaan (7) menjadi:
δu δt + g δh δx + c f u = 0 (8)
Substitusi Persamaan (6) dengan h = D + ζ , dengan ( ζ << h ), dimana D adalah konstan, maka diperoleh persamaan kontinuitas
δζ δt + D δu δx = 0 (9)
Gambar 3: Gelombang Pasang Surut Laut
Substitusi Persamaan (8) dengan h = D + ζ , dengan ( ζ << h ), diperoleh persamaan gerak
δu δt + g δζ δx + c f u = 0 (10)
dimana : ζ ( x, t ) = Ketinggian air diatas referensi acuan. u ( x, t ) = Kecepatan rata-rata diatas penampang yang melintang. D = Kedalaman laut (m). c f = Konstanta friksi = 0 , 0002 s − 1 .
Kondisi awal dan syarat batas dari persamaan gelombang pasang surut diberikan sebagai berikut:
ζ ( x = 0 , t ) = ζ b ( t ) , ζ ( x, t = 0) , u ( x, t = 0) = 0 , u ( x = L, t ) = 0 (11)
Dalam melakukan simulasi Persamaan (10) didiskritkan dengan menggunakan be- da maju terhadap waktu dan beda pusat terhadap posisi sehingga diperoleh
ζ i,k +1 − ζ i,k ∆ t + D u i +1 ,k − u i − 1 ,k 2∆ x = 0 (12)
Persamaan (12) diatas dapat ditulis dalam bentuk:
ζ i,k +1 = ζ i,k − D ∆ t 2∆ x u i +1 ,k + D ∆ t 2∆ x u i − 1 ,k (13)
dengan i = k = 1 , 2 , · · · , n .
Syarat batas dan syarat awal pada Persamaan (11) adalah : ζ ( x = 0 , t ) = ζ b ( t ), ζ ( x, t = 0) = 0, u ( x, t = 0) = 0, u ( x = L, t ) = 0, sehingga ζ 0 ,k = ζ b , ζ i, 0 = 0, u i, 0 = 0, u n,k = 0. Dengan demikian untuk i = 1 → ζ 1 ,k +1 = ζ 1 ,k − D ∆ t 2∆ x u 2 ,k + D ∆ t 2∆ x u 0 ,k .

## Simulasi Model Gelombang Pasang Surut dengan Metode Beda Hingga
Karena δu δx = 0 pada saat x = 0, maka u 0 ,k = u 1 ,k . Sehingga pada saat i = n − 1 menjadi:
ζ n − 1 ,k +1 = ζ n − 1 ,k + D ∆ t 2∆ x u n − 2 ,k
Persamaan (10) yang memuat konsep gerak linier dengan koefisien friksi, harus dijabarkan dalam bentuk diskrit sesuai dengan titik-titik tinjauan pada grid. Ben- tuk diskritisasi persamaan gerak 10 diatas, yaitu:
u i,k +1 − u i,k ∆ t + g ζ i +1 ,k − ζ i − 1 ,k 2∆ x + c f u i,k = 0 (14)
Persamaan (14) diatas dapat ditulis dalam bentuk:
u i,k +1 = (1 − ∆ tc f ) u i,k − g ∆ t 2∆ x ζ i +1 ,k + g ∆ t 2∆ x ζ i − 1 ,k (15)
dimana i = k = 0 , 1 , 2 , · · · , n . Seperti pada persamaan kontinuitas, pada persa- maan gerak juga terdapat beberapa syarat batas dan kondisi awal, yaitu pada Persamaan (11), dimana ζ ( x = 0 , t ) = ζ b ( t ), ζ ( x, t = 0) = 0, u ( x, t = 0) = 0, u ( x = L, t ) = 0, sehingga ζ (0 , k ) = ζ b , ζ ( i, 0) = 0, u n,k = 0. Karena δζ δx = 0 pada saat x = n , maka ζ n,k = ζ n − 1 ,k . Sehingga pada saat persamaannya menjadi:
u n − 1 ,k +1 = (1 − ∆ tc f ) u n − 1 ,k − g ∆ t 2∆ x ζ n − 1 ,k + g ∆ t 2∆ x ζ n − 2 ,k
Dari pendiskritan diatas, kemudian dibentuk menjadi suatu matrik untuk menye- lesaikan suatu sistem persamaan x k +1 = Ax k + B dengan memisalkan a = 1 − ∆ tc f , b = g ∆ t 2∆ x , c = D ∆ t 2∆ x .
Dalam meyelesaikan sistem persamaan x k +1 = Ax k + B diatas, supaya penyele- saian ini dapat stabil maka diberikan batasan sebagai berikut:
1. 1 − ∆ tc f ≤ 1
2. g ∆ t 2∆ x ≤ 1
3. D ∆ t 2∆ x ≤ 1

## 4. Simulasi Model Gelombang Pasang Surut
Pada permasalahan mengenai gelombang pasang surut air laut telah diperoleh model dengan menggunakan asumsi-asumsi untuk mempermudah penghitungan dan telah dibuat bentuk pendiskritannya pada Persamaan (12) untuk persamaan
kontinuitas dan Persamaan (14) untuk persamaan gerak. Pada model ini kedala- man yang diperhitungkan adalah 10 meter karena dibawah itu dianggap aliran- nya terlalu deras. Pada beberapa eksperimen mengenai gelombang pasang su- rut air laut, nilai-nilai parameter diasumsikan sebagai berikut: Kedalaman laut ( D = 10 m ), konstanta friksi ( c f = 0 . 0002), jumlah grid ( N = 40 grid, jarak terjauh ( L = 60 m ) dan perubahan jarak (∆ x = L/ ( N + 1) = 60 / 41 m . Agar penyelesaian dari permasalahan diatas dapat stabil, maka diberi batasan sebagai berikut:
1. 1 − ∆ tc f ≤ 1 Sehingga nilai dari ∆ t dapat diperoleh sebagai berikut: ∆ tc f ≤ 0 dimana c f = 0 . 0002, maka ∆ t ≥ 0.
2. g ∆ t 2∆ x ≤ 1
Setelah dimasukkan nilai g = 10 m/dt 2 dan ∆ x = L/ ( N + 1) = 60 / 41 m = 1 , 46 m , diperoleh: 10∆ t 2(1 , 46) ≤ 1, maka ∆ t ≤ 0 , 292 detik.
3. D ∆ t 2∆ x ≤ 1 Setelah dimasukkan nilai D = 10 m dan ∆ x = L/ ( N +1) = 60 / 41 m = 1 , 46 m , diperoleh: ∆ t ≤ 0 , 292 detik.
Maka perubahan waktu untuk tiap iterasi adalah 0 ≤ ∆ t ≤ 0 , 292 detik.
Gambar 4: Plot nilai ketinggian dan kecepatan gelombang dengan dengan perubahan waktu 0,02 detik
Gambar 4 diatas menunjukkan nilai ketinggian dan kecepatan gelombang pasang surut air laut berdasarkan pembentukan matrik dengan pendekatan nu- merik beda hingga. Dengan ketinggian awal sebesar 0 , 5 meter, terlihat ketinggian

## Simulasi Model Gelombang Pasang Surut dengan Metode Beda Hingga
gelombang naik turun di sekitar nilai ketinggian awalnya. Ketinggian gelombang mencapai 1 , 283 meter dengan kecepatan 0 , 206 m/s setelah berjalan selama 200 satuan waktu (4 detik) dengan penjalaran gelombang ke arah sumbu x dari pusat gelombang sejauh 20 satuan (29,2 meter).
Gambar 5: Plot nilai ketinggian dan kecepatan gelombang dengan dengan perubahan waktu 0,02 detik
Pada Gambar 5 diatas terlihat nilai ketinggian dan kecepatan gelombang dengan ketinggian awal 1m. Ketinggian gelombang terlihat naik turun di se- kitar ketinggian 1 meter. Ketinggian gelombang mencapai 2,565 meter dengan kecepatan 0,415 m/det setelah berjalan selama 200 satuan waktu (4 detik) dengan penjalaran gelombang ke arah sumbu x dari pusat gelombang sejauh 20 satuan (29,2 meter).

## 5. Kesimpulan
Berdasarkan analisa dan pembahasan diatas dapat diambil beberapa kesimpulan sebagai berikut:
1. Untuk memprediksi ketinggian dan kecepatan gelombang pasang surut air laut, dapat menggunakan pendekatan numerik dengan metode beda hingga, yaitu menggunakan beda maju terhadap perubahan waktu dan beda pusat terhadap perubahan jarak.
2. Semakin tinggi nilai perubahan waktu (waktu tempuh gelombang), maka semakin tinggi pula ketinggian dan kecepatan gelombang pasang surut air laut tersebut.

## 3. Ketinggian gelombang bergerak naik turun sepanjang nilai ketinggian awal- nya.

## Pustaka
[1] Irmawanti,Wahyu, (2003), Estimasi Filter Kovariansi Akar Kuadrat Dengan Rang Tereduksi Pada Penyebaran Limbah Cair di Laut, Tugas Akhir, Juru- san Matematika FMIPA-ITS, Surabaya.
[2] Smith, G.D., (1985), Numerical Solution of Partial Differential Equation: Finite Difference Methods, Oxford University Press, New York.
[3] Van Rijn, Leo C., (1994), Principles of Fluid Flow and Surface Waves in Rivers, Estuaries, Seas and Oceans, Aqua Publications, Amsterdam.
[4] Verlaan, Martinus, (1998), Efficient Kalman Filtering Algirithms for Hydro- dynamic Models, Technische Universeteit Delft, Delft.