Limits: Journal of Mathematics and Its Applications E-ISSN: 2579-8936 P-ISSN: 1829-605X Vol. 18, No. 2, Nopember 2021, 95-107 DOI: http://dx.doi.org/10.12962/limits.v18i2.6170

## Simulasi Numerik Persamaan Gelombang Air Dangkal 1D dengan Topografi Tidak Datar Menggunakan Metode Beda Hingga
Ririn Setiyowati 1 , Venna Ade Riestiana 2
1,2 Universitas Sebelas Maret; Jl.Ir. Sutami No.36 A, Puncangsawit, Jebres, Surakarta
1,2 Program Studi Matematika FMIPA UNS Surakarta e-mail: ririnsetiyowati@staff.uns.ac.id, vennaaderiestiana.av@gmail.com
Diajukan: 28 Nopember 2019, Diperbaiki: 22 April 2021, Diterima: 19 Mei 2021

## Abstrak
Tsunami merupakan salah satu contoh gelombang panjang yang dapat dimodelkan  secara matematis menggunakan persamaan gelombang air dangkal.  Persamaan gelombang air dangkal dikonstruksi berdasarkan hukum konservasi massa dan momentum. Persamaan gelombang air dangkal 1D terdiri dari variabel ruang 𝑥𝑥 dan variabel waktu 𝑡𝑡 . Metode beda hingga secara umum dapat digunakan untuk menentukan solusi numerik dari persamaan gelombang air dangkal nonlinear. Pada metode beda hingga terdapat beberapa skema salah satunya adalah skema Lax-Friedrichs.  Dalam penelitian ini dilakukan kontruksi ulang persamaan gelombang air dangkal 1D dengan topografi tidak datar. Kemudian ditentukan solusi numerik dari persamaan gelombang air dangkal 1D menggunakan metode beda hingga skema Lax-Friedrichs.  Simulasi numerik dilakukan dengan tujuan untuk mengetahui hubungan antara nilai awal amplitudo gelombang dan fungsi topografi terhadap tinggi gelombang tsunami yang dihasilkan. Hasil simulasi menunjukkan bahwa semakin besar amplitudo gelombang awal dan semakin curam fungsi topografi yang diberikan maka semakin besar tinggi gelombang yang dihasilkan waktu yang dibutuhkan gelombang untuk sampai ke daratan semakin singkat, dan jangkauan  gelombang masuk ke daratan semakin jauh.
Kata Kunci : persamaan gelombang air dangkal, topografi, metode beda hingga, Lax-Friedrichs

## Abstract
Tsunami is an example of long waves that can be mathematically modeled using shallow water wave equations. Shallow water wave equation is constructed based on the conservation law of mass and momentum. The 1D shallow water wave equation consists of the space variable x and the time variable t. Finite difference method can be used generally to determine the numerical solution of a nonlinear shallow water wave equation. Finite difference method has several schemes, one of them is Lax- Friedrichs scheme.  In this  paper, we reconstruct the 1D shallow water wave equation with non-flat topography. Then, we determine numerical solution using finite different method : Lax-Friedrichs scheme.  Numerical simulations are carried out in order to determine the relationship between the initial value of the wave amplitude and the topographical function of the resulting tsunami wave height. The simulation results show that the greater the initial wave amplitude and the steeper the topographic function given, the greater the resulting wave height, the shorter the time it takes for the waves to reach the land, and the farther the reach of the incoming waves inland.
Keywords : Tsunami, shallow water wave equations, topography, finite difference method, Lax-Friedrichs

## 96 Simulasi Numerik Persamaan Gelombang Air Dangkal 1D dengan Topografi Tidak Datar Menggunakan Metode Beda Hingga

## 1 Pendahuluan
Tsunami adalah gelombang berurutan yang memiliki amplitudo gelombang yang sangat besar. Menurut Synolakis [1], gelombang tsunami merupakan salah satu bentuk gelombang panjang. Perambatan gelombang panjang dapat dimodelkan secara matematis menggunakan persamaan gelombang air dangkal. Flouri et al.[2] menyebutkan bahwa gelombang air dangkal adalah gelombang yang terjadi pada air dangkal yang mana panjang gelombangnya cukup besar dibandingkan kedalamannya.
Pada tahun 2013, Flouri et al. [2], menjelaskan terkait diskritisasi persamaan gelombang air dangkal non linear 2D dengan memperhatikan topograﬁ menggunakan skema tipe Godunov serta disimulasikan pada kasus tsunami Hazard. Roberts  [3]  mengkaji tentang solusi numerik dari hukum konservasi yang diaplikasikan pada persamaan gelombang air dangkal. Crowhurst [4] dan Saiduzzaman et al.[5] mengkaji tentang solusi numerik dari persamaan gelombang air dangkal satu dimensi menggunakan metode beda hingga. Pada tahun 2019, Setiyowati  et al.[6] dalam tulisannya mengkaji tentang penurunan model persamaan gelombang air dangkal 2D dan 1D, serta melakukan simulasi menggunakan metode volume hingga skema Lax-Friedrichs.
Pada penelitian ini disajikan hasil kontruksikan ulang persamaan gelombang air dangkal 1D dari Robert  [3]  dan Setiyowati  et al.[6]. Selanjutnya, ditentukan penyelesaian numerik persamaan gelombang air dangkal 1D menggunakan metode beda hingga skema Lax-Friedrichs. Selain itu, pada penelitian ini dilakukan simulasi  numerik dengan tujuan untuk mengetahui hubungan antara nilai awal amplitudo gelombang dan fungsi topografi terhadap tinggi gelombang tsunami yang dihasilkan. Selain itu, dilakukan interpretasi dari hasil penerapan dan simulasi numerik gelombang air dangkal pada kasus gelombang pantai dengan topograﬁ tidak datar.

## 2 Metode Penelitian
Dalam melakukan penelitian digunakan metode kajian pustaka. Kajian pustaka dilakukan dengan mengumpulkan bahan referensi berupa buku, dan jurnal mengenai persamaan gelombang air dangkal, konstruksi ulang model, simulasi, dan interpretasi.  Pada tahap konstruksi ulang model diawali dengan penjelasan terkait persamaan gelombang air dangkal 2D menggunakan hukum konservasi massa dan hukum konservasi momentum. Selanjutnya, diasumsikan variabel 𝑦𝑦 diabaikan sehingga diperoleh persamaan gelombang air dangkal 1D. Selanjutnya, model tersebut diterapkan pada kasus tsunami selat sunda dan diselesikan secara numerik menggunakan metode beda hingga skema Lax-Friedrichs. Berikut diberikan uraian tentang persamaan
gelombang air dangkal 2D,  persamaan gelombang air dangkal 1D, dan algoritme beda hingga skema Lax-Friedrichs.

## 2.1 Persamaan Gelombang Air Dangkal 2D
Model persamaan gelombang air dangkal 2D mengasumsikan proﬁl gelombang air sebagai masalah dua dimensi yakni sebagai fungsi dua variabel ruang yaitu x dan y dengan t adalah variabel waktu non-negatif yang mana proses penurunannya mengacu pada Robert [3]  dan Setiyowati [6].
Setiyowati [6] memaparkan bahwa persamaan gelombang air dangkal berlaku untuk ﬂuida homogen yang memiliki massa jenis konstan, bersifat inviscid (tidak kental), incompressible (tidak dapat dimampatkan), dan mengalir secara irrotational. Berikut notasi-notasi dasar yang digunakan pada persamaan gelombang air dangkal 2D yang dapat dilihat pada Gambar 1.
Gambar 1 . Struktur dasar persamaan gelombang air dangkal 2D oleh Setiyowati[6]
Berdasarkan Gambar 1, 𝐷𝐷 menyatakan ketinggian rata-rata permukaan air (diukur dari dasar  topograﬁ  bawah)  dan  bernilai  konstan, 𝑏𝑏 ( 𝑥𝑥 , 𝑦𝑦 ) menyatakan  fungsi  topograﬁ  bawah, 𝐻𝐻 ( 𝑥𝑥 , 𝑦𝑦 )   = 𝐷𝐷 − 𝑏𝑏 ( 𝑥𝑥 , 𝑦𝑦 ) menyatakan kedalaman air di bawah rata-rata permukaan air, 𝜂𝜂 ( 𝑥𝑥 , 𝑦𝑦 , 𝑡𝑡 ) menyatakan fungsi perpindahan air di atas rata-rata permukaan air, ℎ ( 𝑥𝑥 , 𝑦𝑦 , 𝑡𝑡 )  = 𝐻𝐻 ( 𝑥𝑥 , 𝑦𝑦 ) + 𝜂𝜂 ( 𝑥𝑥 , 𝑦𝑦 , 𝑡𝑡 ) menyatakan total kedalaman air, 𝑢𝑢 ( 𝑥𝑥 , 𝑦𝑦 , 𝑡𝑡 ) menyatakan kecepatan pada arah 𝑥𝑥 , dan 𝑣𝑣 ( 𝑥𝑥 , 𝑦𝑦 , 𝑡𝑡 ) menyatakan kecepatan pada arah 𝑦𝑦 . Persamaan gelombang air dangkal diturunkan dari hukum konservasi.

## Hukum Konservasi Massa
Hukum konservasi massa adalah prinsip yang menyatakan bahwa laju perubahan massa ﬂuida di suatu ruang pada selang waktu 𝑡𝑡 harus sama dengan perbedaan laju massa yang masuk dan laju massa yang keluar.

## 98 Simulasi Numerik Persamaan Gelombang Air Dangkal 1D dengan Topografi Tidak Datar Menggunakan Metode Beda Hingga
Dimisalkan 𝜌𝜌 menyatakan  massa  jenis  suatu  ﬂuida,  diambil  permukaan  ﬂuida  yaitu ([ 𝑥𝑥 0 , 𝑥𝑥 1 ], [ 𝑦𝑦 0 , 𝑦𝑦 1 ]) dan dideﬁnisikan suatu massa ﬂuida merupakan perkalian antara massa jenis ﬂuida 𝜌𝜌 dan ketinggian ﬂuida ℎ ( 𝑥𝑥 , 𝑦𝑦 , 𝑡𝑡 ) , sehingga diperoleh total massa pada permukaan ﬂuida ([ 𝑥𝑥 0 , 𝑥𝑥 1 ], [ 𝑦𝑦 0 , 𝑦𝑦 1 ]) pada waktu ke 𝑡𝑡 yaitu
𝑚𝑚 total = � � 𝜌𝜌 𝑦𝑦 1 𝑦𝑦 0 𝑥𝑥 1 𝑥𝑥 0 ℎ ( 𝑥𝑥 , 𝑦𝑦 , 𝑡𝑡 ) 𝑑𝑑𝑦𝑦𝑑𝑑𝑥𝑥 .
Total massa ﬂuida berubah sesuai dengan aliran masuk dan aliran keluar ﬂuida. Laju aliran massa ﬂuida adalah massa ﬂuida dikalikan dengan kecepatan rata-rata ﬂuida u untuk variabel x dan v untuk variabel y .  Menggunakan teorema fundamental kalkulus, perubahan total massa ﬂuida yang berlaku untuk cell yang dipilih dalam ruang dan waktu dan massa jenis ρ bernilai konstan adalah
ℎ 𝑡𝑡 + ( ℎ𝑢𝑢 ) 𝑥𝑥 + ( ℎ𝑣𝑣 ) 𝑦𝑦 = 0. (1)
Persamaan (1) merupakan hukum konservasi massa dari persamaan gelombang air dangkal 2D.

## Hukum Konservasi Momentum
Hukum Newton II menyatakan bahwa perubahan momentum terhadap waktu dari suatu sistem sama  dengan  total  gaya yang bekerja pada ﬂuida. Gaya yang diperhatikan disini adalah gaya tekanan.
Dideﬁnisikan total momentum pada permukaan Ω 𝑡𝑡 yaitu
� 𝐩𝐩 d 𝐱𝐱 Ω 𝑡𝑡
dengan Ω 𝑡𝑡 = {( 𝑥𝑥 , 𝑦𝑦 , 𝑡𝑡 ) ∈ ℝ | 𝑥𝑥 0 ≤ 𝑥𝑥 ≤ 𝑥𝑥 1 , 𝑦𝑦 0 ≤ 𝑦𝑦 ≤ 𝑦𝑦 1 , 𝑡𝑡 0 ≤ 𝑡𝑡 ≤ 𝑡𝑡 1 } , 𝐩𝐩 = 𝐮𝐮ℎ , dan 𝐱𝐱 = ( 𝑥𝑥 , 𝑦𝑦 ) . Hukum Newton II dapat ditulis secara matematis sebagai  berikut
𝐹𝐹 = 𝑑𝑑𝐩𝐩 𝑑𝑑𝑡𝑡
dengan 𝐹𝐹 menyatakan gaya (Newton) dan p menyatakan momentum (kg.m/s). Diasumsikan bahwa tekanan  ﬂuida  dalam  domain Ω meningkat  secara  linear  terhadap kedalaman  ﬂuida  (tekanan hidrostatis), dan tekanan pada permukaan bebas ﬂuida adalah konstan sehingga diperoleh
𝐩𝐩 𝑡𝑡 + ∇ ∙ ( 𝐩𝐩 ⊗ 𝐮𝐮 ) + ∇ � 1 2 𝑔𝑔ℎ 2 � = −𝑔𝑔ℎ∇𝑏𝑏 (2)
dengan 𝐩𝐩 ⊗ 𝐮𝐮 adalah tensor product dan 𝑔𝑔ℎ𝛻𝛻𝑏𝑏 adalah gaya pada topograﬁ bawah. Persamaan (2) merupakan persamaan dari hukum konservasi momentum. Oleh karena, 𝐩𝐩 = 𝐮𝐮ℎ dengan 𝐮𝐮 =
( 𝒖𝒖 , 𝒗𝒗 ) dan 𝐱𝐱 = ( 𝑥𝑥 , 𝑦𝑦 ), maka Persamaan (2) dapat diuraikan menjadi sistem dari dua persamaan diferensial yaitu
( ℎ𝑢𝑢 ) 𝑡𝑡 + �ℎ𝑢𝑢 2 + 1 2 𝑔𝑔ℎ 2 � 𝑥𝑥 + ( ℎ𝑢𝑢𝑣𝑣 ) 𝑦𝑦 = 0 ( ℎ𝑣𝑣 ) 𝑡𝑡 + ( ℎ𝑢𝑢𝑣𝑣 ) 𝑥𝑥 + �ℎ𝑣𝑣 2 + 1 2 𝑔𝑔ℎ 2 � 𝑦𝑦 = 0. (3)
Berdasarkan pembahasan di atas, diperoleh persamaan gelombang air dangkal 2D dari Persamaan (1) dan (3), atau dapat dituliskan dalam bentuk
ℎ 𝑡𝑡 + ( ℎ𝑢𝑢 ) 𝑥𝑥 + ( ℎ𝑣𝑣 ) 𝑦𝑦 = 0 ( ℎ𝑢𝑢 ) 𝑡𝑡 + �ℎ𝑢𝑢 2 + 1 2 𝑔𝑔ℎ 2 � 𝑥𝑥 + ( ℎ𝑢𝑢𝑣𝑣 ) 𝑦𝑦 = −𝑔𝑔ℎ𝑏𝑏 𝑥𝑥 ( ℎ𝑣𝑣 ) 𝑡𝑡 + ( ℎ𝑢𝑢𝑣𝑣 ) 𝑥𝑥 + �ℎ𝑣𝑣 2 + 1 2 𝑔𝑔ℎ 2 � 𝑦𝑦 = −𝑔𝑔ℎ𝑏𝑏 𝑦𝑦 .

## 2.2 Persamaan Gelombang Air Dangkal 1D
Persamaan gelombang air dangkal 1D hanya bergantung pada variabel ruang x dan variabel waktu 𝑡𝑡 yang berarti notasi untuk fungsi topograﬁ, perpindahan air di atas permukaan air, dan total kedalaman air masing-masing adalah 𝑏𝑏 ( 𝑥𝑥 ), 𝜂𝜂 ( 𝑥𝑥 , 𝑡𝑡 ), dan ℎ ( 𝑥𝑥 , 𝑡𝑡 ) .  Total  kedalaman air pada posisi 𝑥𝑥 dan pada waktu ke 𝑡𝑡 dinyatakan dalam persamaan berikut
ℎ ( 𝑥𝑥 , 𝑡𝑡 )   = 𝐷𝐷 − 𝑏𝑏 ( 𝑥𝑥 ) + 𝜂𝜂 ( 𝑥𝑥 , 𝑡𝑡 ).
Persamaan gelombang air dangkal 1D sama dengan persamaan gelombang air dangkal 2D namun ditambah asumsi bahwa  komponen arah 𝑦𝑦 diabaikan, sehingga diperoleh  model persamaan gelombang air dangkal satu dimensi (1D) sebagai berikut
ℎ 𝑡𝑡 + ( ℎ𝑢𝑢 ) 𝑥𝑥 = 0 ( ℎ𝑢𝑢 ) 𝑡𝑡 + �ℎ𝑢𝑢 2 + 1 2 𝑔𝑔ℎ 2 � 𝑥𝑥 = −𝑔𝑔ℎ𝑏𝑏 𝑥𝑥 (4)
Dengan mensubstitusikan ℎ = 𝜂𝜂 + 𝐷𝐷 − 𝑏𝑏 ke Persamaan (4) diperoleh Persamaan (5) yaitu
𝜂𝜂 𝑡𝑡 + ( 𝑢𝑢 ( 𝜂𝜂 + 𝐷𝐷 − 𝑏𝑏 )) 𝑥𝑥 = 0 𝑢𝑢 𝑡𝑡 + � 1 2 𝑢𝑢 2 + 𝑔𝑔𝜂𝜂� 𝑥𝑥 = 0. (5)
100 Simulasi Numerik Persamaan Gelombang Air Dangkal 1D dengan Topografi Tidak Datar Menggunakan Metode Beda Hingga

## 2.3 Algoritme Beda Hingga Skema Lax-Friedrichs
Dimisalkan 𝑥𝑥 ∈ [ 𝑥𝑥 𝑚𝑚𝑚𝑚𝑚𝑚 , 𝑥𝑥 𝑚𝑚𝑚𝑚𝑥𝑥 ] dan 𝑡𝑡 ∈ [0, 𝑇𝑇 ] . Pada metode beda hingga, diskritisasi grid untuk variabel 𝑥𝑥 dibagi 𝑁𝑁 interval dengan jarak yang sama sehingga terdapat 𝑁𝑁 + 1 titik dengan 𝑥𝑥 = 𝑖𝑖∆𝑥𝑥 , ∆𝑥𝑥 = 𝑥𝑥 𝑚𝑚𝑚𝑚𝑚𝑚 𝑁𝑁 , dan 𝑖𝑖 = 0,1, … , 𝑁𝑁 . Kemudian untuk variabel 𝑡𝑡 dibagi menjadi 𝑀𝑀 interval sehingga terdapat 𝑀𝑀 + 1 titik dengan 𝑡𝑡 = 𝑗𝑗∆𝑡𝑡 , ∆𝑡𝑡 = 𝑇𝑇 𝑀𝑀 , dan 𝑗𝑗 = 0,1, … , 𝑀𝑀 . Nilai ∆𝑡𝑡 juga dapat ditentukan menggunakan formula
∆𝑡𝑡 = 𝜆𝜆∆𝑥𝑥 max ( 𝑐𝑐 1 , 𝑐𝑐 2 )
dengan 𝜆𝜆 merupakan nilai Courant dan
𝑐𝑐 1 , 2 = 𝑢𝑢 0 ± �𝑔𝑔ℎ 0 , ℎ 0 = 𝐷𝐷 + 𝐴𝐴
dengan 𝑔𝑔 = 9.81 m/s 2 , serta 𝑢𝑢 0 merupakan kecepatan awal gelombang pada saat 𝑡𝑡 = 0 .
Skema Lax-Friedrichs merupakan bentuk pengembangan dari skema eksplisit sehingga untuk menyelesaikan  persamaan diferensial parsial digunakan pendekatan beda maju  untuk turunan terhadap variabel waktu 𝑡𝑡 , dan pendekatan beda pusat untuk turunan terhadap variabel ruang 𝑥𝑥 . Pembagian partisi digunakan untuk memudahkan perhitungan secara numerik pada skema Lax- Friedrichs dapat dilihat pada Gambar 2.

## Gambar 2 . Skema Lax-Friedrichs
Dasar dari skema Lax-Friedrichs adalah nilai dari ( 𝑖𝑖 , 𝑗𝑗 ) pada skema eksplisit diganti dengan rata- rata dari ( 𝑖𝑖 +  1, 𝑗𝑗 ) dan ( 𝑖𝑖 − 1, 𝑗𝑗 ) sehingga  diperoleh bentuk diskrit  dari 𝜂𝜂 ( 𝑥𝑥 , 𝑡𝑡 ), 𝑢𝑢 ( 𝑥𝑥 , 𝑡𝑡 ) , dan 𝑏𝑏 ( 𝑥𝑥 ) adalah
𝜂𝜂 𝑚𝑚 𝑗𝑗 = 𝜂𝜂 𝑚𝑚+1 𝑗𝑗 + 𝜂𝜂 𝑚𝑚−1 𝑗𝑗 2 , 𝑢𝑢 𝑚𝑚 𝑗𝑗 = 𝑢𝑢 𝑚𝑚+1 𝑗𝑗 + 𝑢𝑢 𝑚𝑚−1 𝑗𝑗 2 , 𝑏𝑏 𝑚𝑚 = 𝑏𝑏 𝑚𝑚+1 + 𝑏𝑏 𝑚𝑚−1 2 , (6)

## 101
Pada persamaan gelombang air dangkal 1D (6) terdapat lima turunan yaitu turunan pertama dari 𝜂𝜂 dan 𝑢𝑢 terhadap 𝑡𝑡 serta turunan pertama dari 𝜂𝜂 , 𝑢𝑢 , dan 𝑏𝑏 terhadap 𝑥𝑥 . Menggunakan pendekatan beda maju untuk turunan terhadap variabel waktu 𝑡𝑡 dan pendekatan beda pusat untuk turunan terhadap variabel ruang 𝑥𝑥 , diperoleh
𝜂𝜂 𝑡𝑡 = 𝜂𝜂 𝑚𝑚 𝑗𝑗+1 − 𝜂𝜂 𝑚𝑚 𝑗𝑗 ∆𝑡𝑡 , 𝑢𝑢 𝑡𝑡 = 𝑢𝑢 𝑚𝑚 𝑗𝑗+1 − 𝑢𝑢 𝑚𝑚 𝑗𝑗 ∆𝑡𝑡 , 𝜂𝜂 𝑥𝑥 = 𝜂𝜂 𝑚𝑚+1 𝑗𝑗 − 𝜂𝜂 𝑚𝑚−1 𝑗𝑗 2 ∆𝑥𝑥 , 𝑢𝑢 𝑥𝑥 = 𝑢𝑢 𝑚𝑚+1 𝑗𝑗 − 𝑢𝑢 𝑚𝑚−1 𝑗𝑗 2 ∆𝑥𝑥 , 𝑏𝑏 𝑥𝑥 = 𝑏𝑏 𝑚𝑚+1 − 𝑏𝑏 𝑚𝑚−1 2 ∆𝑥𝑥 . (7)
Persamaan (6) dan (7) disubstitusikan ke dalam Sistem Persamaan (5) sehingga diperoleh penyelesaian numerik persamaan gelombang air dangkal 1D menggunakan metode beda hingga skema Lax-Friedrichs
𝜂𝜂 𝑚𝑚 𝑗𝑗+1 = � 𝜂𝜂 𝑚𝑚+1 𝑗𝑗 + 𝜂𝜂 𝑚𝑚−1 𝑗𝑗 2 � − ∆𝑡𝑡 2 ∆𝑥𝑥 (( 𝜂𝜂 𝑚𝑚+1 𝑗𝑗 + 𝐷𝐷 − 𝑏𝑏 𝑚𝑚+1 ) 𝑢𝑢 𝑚𝑚+1 𝑗𝑗 − ( 𝜂𝜂 𝑚𝑚−1 𝑗𝑗 + 𝐷𝐷 − 𝑏𝑏 𝑚𝑚−1 ) 𝑢𝑢 𝑚𝑚−1 𝑗𝑗 )
𝑢𝑢 𝑚𝑚 𝑗𝑗+1 = � 𝑢𝑢 𝑚𝑚+1 𝑗𝑗 + 𝑢𝑢 𝑚𝑚−1 𝑗𝑗 2 � − 𝑔𝑔∆𝑡𝑡 2 ∆𝑥𝑥 ( 𝜂𝜂 𝑚𝑚+1 𝑗𝑗 − 𝜂𝜂 𝑚𝑚−1 𝑗𝑗 ) − ∆𝑡𝑡 2 ∆𝑥𝑥 � �𝑢𝑢 𝑚𝑚+1 𝑗𝑗 � 2 + �𝑢𝑢 𝑚𝑚+1 𝑗𝑗 � 2 2 � . (8)

## 3  Hasil dan Pembahasan
Pada penelitian ini dilakukan simulasi numerik persamaan gelombang air dangkal 1D pada persamaan  (5) terhadap kasus tsunami berdasarkan amplitudo gelombang dan memisalkan fungsi topograﬁ dasar laut 𝑏𝑏 ( 𝑥𝑥 ) yang mendekati bentuk topografi Selat Sunda pada penelitian Giachetti et al. [7]. Selain itu nilai amplitudo gelombang awal juga diambil dari Giachetti et al. [7]. Tujuan dari simulasi ini dilakukan untuk mengetahui hubungan antara nilai awal amplitudo gelombang dan fungsi topografi terhadap tinggi gelombang tsunami yang dihasilkan.
Diambil 𝑥𝑥 ∈ [0, 100000] dengan garis pantai Banten terletak pada 𝑥𝑥 = 50000 meter, dan posisi longsoran gunung Anak Krakatau terletak pada 𝑥𝑥 = 100000 meter. Topografi pada Giachetti et al. [7] diambil untuk fungsi topografi pada jarak 9 kilometer dari posisi longsoran mengarah ke garis pantai, dan untuk fungsi topografi setelahnya dimisalkan berupa fungsi linear, sehingga diperoleh fungsi topografi dasar laut seperti pada Gambar 3.
Gambar 3 . Topografi dasar laut Selat Sunda oleh Giachetti et al. [7]
Berdasarkan Gambar 3,  dikembangkan topografi dasar laut selat sunda didekati dengan empat fungsi topograﬁ b ( x ) yaitu
𝑏𝑏 1 ( 𝑥𝑥 ) = ⎩ ⎪ ⎪ ⎨ ⎪ ⎪ ⎧ 200,                                      0 ≤ 𝑥𝑥 ≤ 50000;
−𝑥𝑥 125 + 600,                          50000 < 𝑥𝑥 ≤ 75000;
0   ,                                      75000 < 𝑥𝑥 ≤ 100000;
− ( 𝑥𝑥 − 105000) 2 400000 + 200, 100000 < 𝑥𝑥 ≤ 105000;
200,                                       105000 < 𝑥𝑥 ≤ 110000;
𝑏𝑏 2 ( 𝑥𝑥 ) = ⎩ ⎪ ⎪ ⎨ ⎪ ⎪ ⎧ 200,                                        0 ≤ 𝑥𝑥 ≤ 50000; − ( 𝑥𝑥 − 50000) 3125000 + 200,        50000 < 𝑥𝑥 ≤ 75000; 0   ,                                         75000 < 𝑥𝑥 ≤ 100000; − ( 𝑥𝑥 − 105000) 2 125000 + 200, 100000 < 𝑥𝑥 ≤ 105000;
200,                                        105000 < 𝑥𝑥 ≤ 110000; 𝑏𝑏 3 ( 𝑥𝑥 ) = ⎩ ⎪ ⎪ ⎨ ⎪ ⎪ ⎧ 200,                                        0 ≤ 𝑥𝑥 ≤ 50000; − ( 𝑥𝑥 − 75000)
3125000 ,                  50000 < 𝑥𝑥 ≤ 75000;
0   ,                                         75000 < 𝑥𝑥 ≤ 100000; − ( 𝑥𝑥 − 105000) 2 125000 + 200, 100000 < 𝑥𝑥 ≤ 105000;
200,                                        105000 < 𝑥𝑥 ≤ 110000; (9)
𝑏𝑏 4 ( 𝑥𝑥 ) = ⎩ ⎪ ⎪ ⎪ ⎨ ⎪ ⎪ ⎪ ⎧ 200,                                                  0 ≤ 𝑥𝑥 ≤ 50000;
−𝑥𝑥 100 + 700,                            50000 < 𝑥𝑥 ≤ 60000; 100 ,                                      60000 < 𝑥𝑥 ≤ 700000;
− ( 𝑥𝑥 − 70000) 2 250000 + 200, 70000 < 𝑥𝑥 ≤ 75000;
0,                                              75000 < 𝑥𝑥 ≤ 100000;
− ( 𝑥𝑥 − 105000) 2 125000 + 200, 105000 < 𝑥𝑥 ≤ 11000;
200,                         105000 < 𝑥𝑥 ≤ 11000;
dengan fungsi topografi 𝑏𝑏 1 ( 𝑥𝑥 ) merupakan fungsi topografi yang berbentuk sloping bed , 𝑏𝑏 2 ( 𝑥𝑥 ) merupakan topografi landai, 𝑏𝑏 3 ( 𝑥𝑥 ) merupakan topografi curam, dan 𝑏𝑏 4 ( 𝑥𝑥 ) merupakan fungsi topografi modifikasi antara 𝑏𝑏 1 ( 𝑥𝑥 ) dan 𝑏𝑏 2 ( 𝑥𝑥 ) namun lebih curam dibandingkan fungsi topografi 𝑏𝑏 2 ( 𝑥𝑥 ) dan lebih landai dibandingkan 𝑏𝑏 1 ( 𝑥𝑥 ) . Berikut diberikan sketsa fungsi topografi dasar laut seperti pada Gambar 4.
Gambar 4 . Fungsi topografi dasar laut 𝑏𝑏 ( 𝑥𝑥 )
Kemudian, diberikan syarat awal gelombang yang berupa gelombang soliter dan kecepatan gelombang pada t = 0 yang diambil dari Setiyowati et al. (2019), yaitu
𝜂𝜂 ( 𝑥𝑥 , 0) = 𝐴𝐴 sech �� 3 4 𝐷𝐷 2 ( 𝑥𝑥 − 97500) � , (10)
dan 𝑢𝑢 ( 𝑥𝑥 , 0) = 0 . Selanjutnya, diambil beberapa nilai awal amplitudo gelombang A yaitu 15, 20, 25, dan 30 meter dengan posisi gelombang awal terletak pada 𝑥𝑥 = 97500 meter. diambil amplitudo gelombang awal 𝐴𝐴 dari Giachetti et al. [7] yaitu 45 meter dengan posisi gelombang awal terletak pada 𝑥𝑥 = 97500 meter.
Selanjutnya, simulasi numerik dilakukan pada domain variabel ruang 𝑥𝑥 ∈ [0, 100000] dan variabel waktu 𝑡𝑡 ∈ [0, 1000] . Diskritisasi dilakukan pada grid terstruktur dengan partisi ruang adalah cell ∆𝑥𝑥 = 100000 𝑁𝑁 dan partisi waktu ∆𝑡𝑡 diperoleh dari kriteria CFL. Diambil banyaknya grid 𝑁𝑁 adalah 2000, dan nilai Courant 𝜆𝜆 dari Setiyowati et al. [6] adalah 0.68 sehingga diperoleh
∆𝑥𝑥 = 50 , ∆𝑡𝑡 = (0.68) ∆𝑥𝑥 max ( 𝑐𝑐 1 , 𝑐𝑐 2 )
dengan 𝑐𝑐 1 , 2 = 𝑢𝑢 0 ± �𝑔𝑔ℎ 0 , ℎ 0 = 𝐷𝐷 + 𝐴𝐴 , 𝑔𝑔 = 9.81 m/s 2 , dan 𝑢𝑢 0 = 0 m/s .
Hasil simulasi  numerik  dievaluasi  pada  posisi  garis  pantai ( x = 50000  meter) dengan nilai awal amplitudo gelombang 𝐴𝐴 yaitu 15, 20, 25, dan 30 meter, fungsi topografi 𝑏𝑏 1 ( 𝑥𝑥 ), 𝑏𝑏 2 ( 𝑥𝑥 ), 𝑏𝑏 3 ( 𝑥𝑥 ) dan 𝑏𝑏 4 ( 𝑥𝑥 ) , syarat awal (10), dan 𝑢𝑢 ( 𝑥𝑥 , 0) = 0 terhadap waktu tempuh ( 𝑡𝑡 ) dan maksimum tinggi gelombang ( max( 𝜂𝜂 ( 𝑥𝑥 , 𝑡𝑡 )) ) dapat dilihat pada Tabel 1.
Tabel 1. Hasil simulasi berdasarkan 𝐴𝐴 dan 𝑏𝑏 ( 𝑥𝑥 ) terhadap waktu tempuh dan tinggi gelombang
𝐴𝐴 𝑏𝑏 1 ( 𝑥𝑥 ) 𝑏𝑏 2 ( 𝑥𝑥 ) 𝑏𝑏 3 ( 𝑥𝑥 ) 𝑏𝑏 4 ( 𝑥𝑥 )
𝑡𝑡 max( 𝜂𝜂 ( 𝑥𝑥 , 𝑡𝑡 )) 𝑡𝑡 max( 𝜂𝜂 ( 𝑥𝑥 , 𝑡𝑡 )) 𝑡𝑡 max( 𝜂𝜂 ( 𝑥𝑥 , 𝑡𝑡 )) 𝑡𝑡 max( 𝜂𝜂 ( 𝑥𝑥 , 𝑡𝑡 )) 15  859.52 2.4359 889.14 2.1331 835.83 2.7416 866.93 2.3834 20  857.02 3.1487 886.29 2.7724 833.59 3.5238 865.07 3.0839 25  855.40 3.8225 884.35 3.3822 831.52 4.2554 862.64 3.7474 30  853.21 4.4612 882.56 3.9643 829.59 4.9437 860.37 4.3766
Berdasarkan Tabel 1 dapat dilihat bahwa  untuk fungsi topografi 𝑏𝑏 3 ( 𝑥𝑥 ) dihasilkan tinggi gelombang yang paling besar dengan kisaran antara 2.7 sampai 4.9 meter, sedangkan untuk kasus fungsi topografi 𝑏𝑏 2 ( 𝑥𝑥 ) dihasilkan tinggi gelombang yang paling kecil dengan kisaran antara 2.1 sampai 3.9 meter. Selain itu untuk fungsi topografi 𝑏𝑏 1 ( 𝑥𝑥 ) dan 𝑏𝑏 4 ( 𝑥𝑥 ) menghasilkan tinggi gelombang dengan kisaran 2.4 sampai 4.4 meter. Hal tersebut berlaku untuk semua nilai

## 105
amplitudo gelombang awal yang diberikan. Tabel 1 menunjukkan juga bahwa fungsi topografi yang menghasilkan gelombang dengan waktu tempuh tersingkat untuk sampai kedaratan adalah fungsi topografi 𝑏𝑏 3 ( 𝑥𝑥 ) , dan yang terlama adalah fungsi topografi 𝑏𝑏 2 ( 𝑥𝑥 ) . Dengan demikian, dapat disimpulkan bahwa semakin landai topografinya, maka  tinggi gelombang yang dihasilkan semakin kecil dan waktu tempuh yang dibutuhkan gelombang untuk sampai ke daratan semakin lama.
Selanjutnya, simulasi numerik dievaluasi pada waktu 𝑡𝑡 = 1000 detik  dengan tujuan untuk melihat seberapa jauh gelombang masuk ke daratan sehingga diperoleh hasil seperti pada Tabel 2.
Tabel 2. Hasil simulasi terhadap posisi dan tinggi gelombang maksimum pada waktu 𝑡𝑡 = 1000 detik
𝐴𝐴 𝑏𝑏 1 ( 𝑥𝑥 ) 𝑏𝑏 2 ( 𝑥𝑥 ) 𝑏𝑏 3 ( 𝑥𝑥 ) 𝑏𝑏 4 ( 𝑥𝑥 )
𝑥𝑥 max( 𝜂𝜂 ( 𝑥𝑥 , 𝑡𝑡 )) 𝑥𝑥 max( 𝜂𝜂 ( 𝑥𝑥 , 𝑡𝑡 )) 𝑥𝑥 max( 𝜂𝜂 ( 𝑥𝑥 , 𝑡𝑡 )) 𝑥𝑥 max( 𝜂𝜂 ( 𝑥𝑥 , 𝑡𝑡 )) 15 44350 2.0294 45150 1.8980 42750 2.1279 44150 2.0145 20 43750 2.6310 45050 2.4677 42650 2.7513 44050 2.6120 25 43600 3.1983 44900 3.0076 42550 3.3375 43950 3.1772 30 43500 4.7393 44860 3.5257 42450 3.8951 43850 3.7163
Berdasarkan hasil pada Tabel 2 dapat dilihat bahwa semakin besar amplitudo gelombang awal maka semakin besar tinggi gelombang yang dihasilkan dan semakin jauh gelombang masuk ke daratan. Pada fungsi topografi 𝑏𝑏 1 ( 𝑥𝑥 ) dapat dilihat bahwa dengan amplitudo awal 30 meter gelombang mencapai posisi terjauh yaitu 43500 meter atau 6.5 kilometer dari garis pantai dengan ketinggian gelombang maksimum adalah 3.7393 meter. Pada fungsi topografi 𝑏𝑏 2 ( 𝑥𝑥 ) dengan amplitudo awal 30 meter gelombang mencapai posisi terjauh yaitu 44850 meter atau 5.15 kilometer dari garis pantai dengan ketinggian gelombang maksimum adalah 3.5257 meter. Pada fungsi topografi 𝑏𝑏 3 ( 𝑥𝑥 ) dengan amplitudo awal 30 meter gelombang mencapai posisi terjauh yaitu 42450 meter atau 7.55 kilometer dari garis pantai dengan ketinggian gelombang maksimum adalah 3.8951 meter. Dan pada fungsi topografi 𝑏𝑏 4 ( 𝑥𝑥 ) dengan amplitudo awal 30 meter gelombang mencapai posisi terjauh yaitu 43850 meter atau 6.15 kilometer dari garis pantai dengan ketinggian gelombang maksimum adalah 3.7163 meter..
Hasil pada Tabel 2 menunjukkan juga bahwa topografi 𝑏𝑏 3 ( 𝑥𝑥 ) yang merupakan topografi curam menghasilkan tinggi gelombang yang paling  besar dan jarak tempuh yang paling jauh dibandingkan tiga topografi lainnya, sehingga dapat disimpulkan bahwa semakin besar amplitudo
106 Simulasi Numerik Persamaan Gelombang Air Dangkal 1D dengan Topografi Tidak Datar Menggunakan Metode Beda Hingga
gelombang awal yang diberikan dan semakin curam topografi dasar lautnya maka semakin besar tinggi gelombang yang dihasilkan dan semakin jauh jangkauan gelombang yang  masuk ke daratan.

## 4 Simpulan
Berdasarkan hasil dan pembahasan diperoleh model persamaan gelombang air dangkal 1D dengan topografi tidak flat  dan algoritme  beda hingga skema Lax-Friedrichs dari  persamaan gelombang air dangkal 1D. Pada penelitian ini persamaan gelombang air dangkal dibatasi pada dimensi satu dengan fungsi topografi yang sederhana. Metode numerik yang digunakan untuk mencari penyelesaian persamaan gelombang air dangkal 1D adalah metode beda hingga skema Lax-Friedrichs. Simulasi numerik dilakukan dengan tujuan untuk mengetahui hubungan antara nilai awal amplitudo gelombang dan fungsi topografi terhadap tinggi gelombang tsunami yang dihasilkan.  Berdasarkan hasil simulasi dapat disimpulkan bahwa semakin besar amplitudo gelombang yang diberikan maka tinggi gelombang yang dihasilkan semakin besar, waktu yang dibutuhkan gelombang untuk sampai ke daratan semakin singkat, dan jangkauan gelombang masuk ke daratan semakin jauh. Semakin curam fungsi topografi yang diberikan maka semakin besar tinggi gelombang yang dihasilkan, waktu yang dibutuhkan gelombang untuk sampai ke daratan semakin singkat, dan jangkauan gelombang masuk ke daratan semakin jauh. Pada penelitian selanjutnya dapat dikembangkan persamaan gelombang air dangkal dua dimensi atau lebih, serta diselesaikan secara numerik menggunakan skema yang lain seperti skema Crank- Nicolson. Selain itu, dapat juga dikembangkan untuk fungsi topografi yang lebih bervariasi atau yang mendekati bentuk topografi pada kasus nyata.

## 5  Daftar Pustaka
[1]  C. E. Synolakis, “The runup of solitary waves,” J. Fluid Mech. , vol. 185, no. May, pp. 523–545, 1987, doi: 10.1017/S002211208700329X.
[2]  E. T. Flouri, N. Kalligeris, G. Alexandrakis, N. A. Kampanis, and C. E. Synolakis, “Application of a finite difference computational model to the simulation of earthquake generated tsunamis,” Appl. Numer. Math. , vol. 67, pp. 111–125, 2013, doi: 10.1016/j.apnum.2011.06.003.
[3]  S. G. Roberts, “Numerical solution of conservation laws applied to the Shallow Water Wave Equations,” vol. 2013, p. 199, 2013.
[4]  P. Crowhurst and Z. Li, “Numerical solutions of one-dimensional shallow water equations,” Proc. - UKSim 15th Int. Conf. Comput. Model. Simulation, UKSim 2013 , no. April, pp. 55–60, 2013, doi: 10.1109/UKSim.2013.63.
[5]  M. Saiduzzaman and S. K. Ray, “Comparison of Numerical Schemes for Shallow Water Equation,” Glob. J. Sci. Front. Res. Math. Decis. Sci. , vol. 13, no. 4, 2013.
[6]  R. Setiyowati and Sumardi, “A Simulation of Shallow Water Wave Equation Using Finite Volume Method: Lax-Friedrichs Scheme,” J. Phys. Conf. Ser. , vol. 1306, no. 1, 2019, doi: 10.1088/1742-6596/1306/1/012022.
[7]  T. Giachetti, R. Paris, K. Kelfoun, and B. Ontowirjo, “Tsunami hazard related to a flank collapse of Anak Krakatau Volcano, Sunda Strait, Indonesia,” Geol. Soc. Spec. Publ. , vol. 361, no. 1, pp. 79–90, 2012, doi: 10.1144/SP361.7.