SEMINAR NASIONAL MATEMATIKA DAN PENDIDIKAN MATEMATIKA UNY 2015 T – 17
Solusi Numerik Persamaan Gelombang Air Dangkal Linear Menggunakan FEM Nikenasih Binatari Universitas Negeri Yogyakarta
[email protected]
Abstrak— Pada tulisan ini akan disajikan solusi numerik persamaan gelombang air dangkal linear menggunakan finite element method (FEM). Solusi numerik diperoleh dengan mengaproksimasi solusinya menggunakan kombinasi linear fungsi basis yang kemudian disebut tent function. Hasil dari simulasi menunjukkan bahwa untuk bentuk akhir pergerakan yang sama, profile gelombang awal bervariasi bergantung pada lama pergerakan dasar. Semakin cepat pergerakan dasar, profile gelombang awal akan semakin mirip dengan bentuk akhir pergerakan dasar. Kata kunci: Pergerakan Dasar, Persamaan Gelombang Air Dangkal Linear, Finite Element Method
I.
PENDAHULUAN
Sejak terjadinya di Aceh, 26 Desember 2004, Gelombang Tsunami menjadi topik yang terkenal di bidang penelitian fluida karena perbedaan karakteristiknya dengan gelombang-gelombang yang lain. Gempa bumi selalu mengawali terjadinya gelombang tsunami. Pada tsunami aceh, seismograf mencatat besaran gempa bumi adalah 8.8 SR akan tetapi kemudian penelitian lebih lanjut menemukan bahwa getaran yang terjadi di Aceh tersebut mencapai 9.1 – 9.3 SR dan kira-kira terjadi selama 10 menit gerakan patahan. Hanya dalam hitungan jam, lebih dari 230.000 orang terbunuh dan lebih dari 5000 orang dinyatakan hilang. Oleh karena itu, tsunami sering juga disebut dengan gelombang pembunuh. Banyak peneliti mempelajari tentang propagasi gelombang tsunami karena tsunami berbeda dalam model pembangkitnya dan dalam karakteristiknya yaitu periode, panjang gelombang dan kecepatan. Tsunami dapat memiliki periode 10 sampai 120 menit , panjang gelombang 100-200 km dan kecepatan 200m/s. Tsunami terjadi karena gangguan impulsif dari air laut yang disebabkan oleh perubahan bentuk pada dasar laut secara tiba-tiba laut sebagai contoh tanah longsor, ledakan gunung berapi, gempa bumi yang terjadi di bawah dan lain-lain. Dari beberapa penyebab Tsunami, gempa merupakan penyebab utamanya . Karakteristik gelombang tsunami ditentukan oleh penyebabnya. Pada generasi pasif, gelombang permukaan awal yang dihasilkan merupakan translasi dari perubahan bentuk akhir dasar laut. Penelitian lebih lanjut oleh Dutykh dan Diaz (2007 ) menganalisa tentang gelombang di permukaan dan gerakan bawah secara analitis menggunakan transformasi Fourier. Hasil penelitian menunjukkan bahwa gelombang awal di permukaan dipengaruhi oleh profil bawah dan durasi gerakan. Diasumsikan bahwa dasar laut mula-mula adalah datar, untuk kasus sesaat yaitu di mana gerakan terjadi dengan cepat , profil gelombang awal hanya translasi sama seperti dalam Passive Generation. Akan tetapi untuk kasus yang lain, semakin lama pergerakan dasar laut maka gelombang yang dihasilkan juga akan semakin berbeda. Teori lanjutan dari generasi pasif biasa disebut dengan generasi aktif. Penelitian ini dibatasi pada fluida yang memenui tiga sifat tidak termampatkan (incompressible), tidak kental (inviscid) dan tidak berrotasi (irrational). Salah satu model pergerakan fluida yang sangat terkenal adalah persamaan air dangkal. Dimisalkan x, t adalah elevasi permukaan laut, x, t potential free surface dan hx, t kedalamaan pada posisi x dan saat t. Mekipun persamaan air dangkal mewakili kondisi fluida secara lengkap, akan tetapi formulanya sangat kompleks. Oleh karena itu, masalah ini kemudian dibatasi pada gelombang dengan amplitude sangat kecil. Pada penelitian ini diambil perubahan dasaru laut tertinggi adalah 1 m. Jadi di sini kita menempatkan tertinggi tidur laut deformasi adalah 1 m . Tujuan dari batas ini adalah untuk membentuk masalah ke dalam bentuk linear sehingga solusi akan mudah dianalisa .
299
ISBN 978-602-73403-0-5
Di sini kita hanya memperhatikan pada elevasi di permukaan. Dengan demikian, dengan menggunakan asumsi pada batas kita mendapatkan dua persamaan yang dikenal sebagai Linear Air Dangkal Persamaan sebagai berikut t x h x t h
1 2 3
t g x,0 0 dan x,0 0
Persamaan (1) adalah persamaan di permukaan, menunjukkan ketidakmungkinan air mengalir melalui permukaan serta sementara (2) adalah persamaan Bernoulli, yang menyatakan bahwa tekanan di permukaan air harus lenyap . Nilai awal (3) menunjukkan bahwa ketinggian datar sebelum gerakan dan juga fluida tidak memiliki potensi di permukaan bebas . Masalah diferensial ini adalah masalah nilai awal. Untuk menganalisis masalah ini, kita dekati masalah dengan menggunakan metode pengali. Metode pengali adalah strategi untuk mengekstraksi informasi dengan mengalikan persamaan diferensial dengan fungsi yang sesuai dan kemudian mengintegrasikan. Metode ini biasanya disebut sebagai formulasi variational atau formulasi lemah. Formulasi ini memungkinkan penggunaan Metode Galerkin, yaitu cara alami memproyeksikan persamaan ke ruang bagian aproksimasi terbatas dimensi . Pada tulisan ini, akan dianalisis elevasi gelombang menggunakan skema numerik. Persamaan Linear Air Dangkal Persamaan akan didiskritisasi secara spasial. Oleh karena itu, ketinggian gelombang dapat digambarkan sebagai kombinasi linear dari fungsi dasar. Di sini kita menggunakan fungsi hat linear sebagai fungsi dasar. Kemudian, integrasi waktu akan dilakukan dengan menggunakan Runge Kutta atau ODESOLVER pada MATLAB .
II.
IMPLEMENTASI NUMERIK
A. Formulasi Variasional Seperti yang disebut sebelumnya bahwa pada masalah nilai awal ini digunakan metode galerkin. Untuk sebarang fungsi v yang terdefinisi pada , kalikan kedua ruas (1a) dan (1b) dengan v sehingga menghasilkan
v d t
x
h x v d t h v d
v d g v d
(5)
t
(4)
Dengan menggunakan integrasi parsial, (4) dapat ditulis ulang sebagai berikut
v d h v t
x
h v d h v d x
x
t
6
Pada bahasan ini, fungsi v disebut dengan fungsi tes. Untuk menyederhanakan penghitungan, dipilih fungi test yang memenuhi h x v 0. Dari sini diperoleh
300
SEMINAR NASIONAL MATEMATIKA DAN PENDIDIKAN MATEMATIKA UNY 2015
v d h v d h v d t
x
x
7
t
B. Diskritisasi Diskritisasi dilakukan pada domain ruang . Dimisalkan domain ruang dipartisi menjadi N sub domain, oleh karena itu terdapat N+1 node x1 , x2 ,, x N , x N 1 . Disini, subdomain sering disebut dengan elemen. Tanpa mengurangi sifat keumuman, disini dipilih elemen seragam dimana semua elemennya mempunyai lebar yang sama yaitu x . Jadi, disini diperoleh x xi 1 xi , i 1,2,, N . Elevasi gelombang digambarkan sebagai kombinasi linear dari fungsi dasar. Pada penelitian ini, fungsi dasar yang digunakan adalah fungsi hat linear sebagai berikut :
x xk 1 x x , if x xk 1, xk k 1 k xk 1 x vk , if x xk , xk 1 xk 1 xk 0 otherwise Ilustrasi dari fungsi tersebut dapat dilihat pada Gambar 1 sebagai berikut :
v1 ( x) v k (x)
v 2 ( x) v k 1 ( x)
x1 x2
v N 1 ( x )
v k 1 ( x)
v N (x)
v3 ( x )
x3
xk 1
x N x N 1
xk xk 1
Gambar 1. Fungsi Tent
Dengan menggunakan fungsi hat linear, atau sering juga disebut dengan fungsi tent, fungsi elevasi gelombang dan fungsi potensial di permukaan dapat diaproksimasi sebagai berikut :
x, t x, t
N 1
ˆ t v x
8
ˆ t v x
9
k
k
k 1 N 1
k
k
k 1
Oleh karena itu, elevasi gelombang dan potensial di permukaan tidak lagi bergantung atas ruagn akan tetapi hanya bergantung atas waktu saja. Substitusikan (8) dan (9) pada (7) dan (5). Kemudian, dipilih N+1 fungsi tent sebagai fungsi test maka akan diperoleh N+1 persamaan sebagai berikut :
N 1 N 1 t ˆk t vk x vs d h x ˆk t vk x x vs d t h vs d . k 1 k 1
301
ISBN 978-602-73403-0-5
N 1 N 1 t ˆk t v k x v s d g ˆ k t v k x v s d k 1 k 1
untuk s 1,2,...,N 1 .
Persamaan tersebut dapat dinyatakan dalam operasi matriks sebagai berikut :
ˆ A t B ˆ c A tˆ gA ˆ
10 11
dimana
A a ks v k v s d , B bks hv' k v' s d, c cs t h vs x d
Misalkan dasar laut berbentuk datar sebelum terjadi gerakan laut. Oleh karena kedalamannya konstan, h0 . Selanjutnya, dimisalkan perubahan kedalaman dapat dinyatakan sebagai hx, t h0 f x S t dimana f adalah perubahan bentuk dasar laut dan S adalah dependensi waktu fungsi . Misalkan T adalah durasi lama pergerakan dasar laut. Pada penelitian ini, beberapa simulasi dilakukan dengan bermacam-macam nilai T. sebagai ilustrasi, berikut merupakan hasil simulasi untuk T = 10 detik dan T = 50 detik dengan fungsi waktu sebagai berikut :
1 t 1 cos ,0 t T S t 2 T 1 , t T
14
Ilustrasi fungsi 14 dapat dilihat pada Gambar 2 berikut :
Gambar 2. Time Dependencies for T = 50.
Karena f adalah fungsi atas domain ruang, oleh karenanya f dapat diaproksimasi sebagai kombinasi linear atas fungsi dasar sebagai berikut :
302
SEMINAR NASIONAL MATEMATIKA DAN PENDIDIKAN MATEMATIKA UNY 2015
N 1
fˆ v x ,
f x
l l
l 1
Darisini diperoleh
N 1
fˆ v x
hx, t h0 S t
12 .
l l
l 1
Substitusikan (12) pada (10) diperoleh
A tˆ h0 B1 S B2 ˆ d S A fˆ 14 A tˆ gA ˆ
13
dengan A, B1 dan B2 adalah matriks tridiagonal. 1 3 1 6 0 A .. .. .. 0
1 6 2 3 1 6
..
0 1 6 2 3 ..
..
..
..
1 6 .. ..
..
.. .. 1 6
..
fˆ1 fˆ2 fˆ1 fˆ2 0 1 B2 .. 2x .. .. 0
..
.. 2 3 1 6
0 .. .. .. 0 1 1 0 0 0 1 2 1 0 1 2 1 .. .. 1 B1 .. .. .. .. .. .. x .. .. .. .. .. .. .. 1 2 1 1 .. .. .. .. 1 1 ( N 1)( N 1) 6 0 1 3 ( N 1)( N 1)
fˆ1 fˆ2 fˆ1 2 fˆ2 fˆ3 fˆ2 fˆ3
0 .. .. fˆ2 fˆ3 fˆ2 2 fˆ3 fˆ4 .. .. ..
..
.. .. ..
.. fˆN 2 2 fˆN 1 fˆN .. .. fˆN 1 fˆN
..
Berdasarkan nilail awal
x,0 0 and x,0 0
Aproksimasi (3) menggunakan (8) dan (9) diperoleh
303
.. fˆN 1 fˆN fˆN 1 fˆN 0 0 .. ..
ISBN 978-602-73403-0-5
N 1
ˆk 0vk x 0 and
k 1
N 1
ˆ 0v x 0 k
k
k 1
Karena vk x bebas linear maka diperoleh 2N +2 nilai awal yaitu ˆ k 0 0 dan ˆk 0 0 untuk k = 1, 2 , … , N+1. Darisini dapat dilihat bahwa (13) dan (14) adalah sistem persamaan diferensial linear order satu. Jadi, terdapat 2N+2 persamaan dengan 2N+2 fungsi waktu yang tidak diketahui, ˆk t dan ˆ t . k
Solusi dari sistem persamaan ini dapat dicari dengan menggunaan software MATLAB, yaitu ODESOLVER. Terdapat bermacam-macam tipe ODESOLVER. Pada penelitian ini digunakan ODE45, yaitu metode untuk menyelesaikan sistem persamaan diferensial berdasarkan rumus pada metod Runge Kutta. Nilai-nilai parameter yang digunakan pada simulasi program ini dapat dilihat pada Tabel 1 berikut. PARAMETERS Parameter
Simbol
Nilai 4
Domain ruang
x
-8x10 sampai 8 x 104 m
Lebar elemen
dx
250 m
Domain waktu
t
0 s sampai 200 s
Panjang interval waktu
dt
0.5 s
Percepatan gravitasi
g
9.81 m2/s
Kedalaman mula-mula
h0
4000 m
Pada penelitian ini digunakan contoh dimana bentuk dasar laut setelah terjadi gempa adalah persegi atau diilustrasikan pada Gambar 3 berikut. Jika nilai dari S sama dengan 1 maka kedalamannya menjadi fungsi atas ruang saja.
Gambar 3. Deformasi Persegi
Hasil simulasi pada elevasi gelombang menggunakan prosedur ODE45 untuk T = 10 detik dapat dilihat pada gambar 4 berikut :
Gambar 4. Elevasi di permukaan saat t = 1s, t = 3s, t = 7s, t = 10s, t = 25s, t = 50s and t = 100s.
304
SEMINAR NASIONAL MATEMATIKA DAN PENDIDIKAN MATEMATIKA UNY 2015
Durasi pergerakan dasar laut 10 detik dan seperti yang dapat dilihat pada Gambar 4 gelombang akan mulai merambat setelah 10 detik sejak gerakan pertama. Amplitudo tertinggi elevasi berkisar 0,85 m, diperoleh pada t = 10 detik. Setelah 10 detik , gelombang merambat ke kanan dan ke kiri dengan amplitudo yang sama , sekitar 0,46 m . Dibandingkan dengan kasus pertama, berikut ini adalah hasil simulasi elevasi permukaan untuk T = 50 detik. Seperti yang dapat dilihat bahwa gelombang akan merambat setelah dasar laut berhenti bergerak dan mencapai amplitudo tertinggi, sekitar setengah dari profil awal 0,41 m, segera setelah itu. Kemudian setelah 50 detik, gelombang merambat ke kanan dan ke kiri dengan sama amplitudo 0,37 m .
Gambar 4. Elevasi di permukaan saat t = 1s, t = 10s, t = 30s, t = 50s, t = 75s, t = 100s and t = 150s.
III.
SIMPULAN DAN SARAN
Dari simulasi dapat disimpulkan bahwa untuk gerakan yang sama, profil gelombang awal akan berbeda tergantung pada berapa lama dasar laut bergerak. Semakin cepat pergerakan laut, profil elevasi permukaan akan terlihat sama. Pada awalnya, selama gerakan , permukaan akan bergerak ke atas tetapi segera setelah dasar laut berhenti bergerak, gelombang akan merambat ke kanan dan kiri. Amplitudo tertinggi dicapai pada saat akhir gerakan keatas .
DAFTAR PUSTAKA [1]
D. Dutykh, F. Dias, “Water wave generated by a moving bottom” in Tsunami and nonlinear waves, springer berlin heidelberg pp. 65-95. 2007. (references)
[2]
I.M Navon. “Finite element simulation of the shallow-water equations model on a limited-area domain”. Appl. Math. Modelling, 1979, Vol 3, October. pp 337-348.
[3]
F. Imamura, A.C Yalciner, G. Ozyurt, “Tsunami modelling manual” 2006. Unpublished.
[4]
M. Asadzadeh. “An introduction to the finite element method for differential equations”. 2010. Unpublished.
[5]
N. Binatari. “Gelombang yang dibangkitkan oleh pergerakan bawah laut”. Prosiding of Seminar Matematika dan Pendidikan Matematika. Jurusan Pendidikan Matematika FMIPA UNY. pp. MT-259, November 2013.
[6]
A. Sayma. “Computational fluid dynamics, 1st Edition”. 2009.
[7]
H. Zhang, et al. “Modeling and visualization of tsunami”. Pure and Applied Geophysics. 165 (2008), pp 475–496.
305
ISBN 978-602-73403-0-5
[8]
M.S Gockenbach. “Understanding and implementing the finite element method”. Society for Industrial and Applied Mathematics. 2006.
306
