BAB II STUDI PUSTAKA. Propagated wave area. Shallow water. Area of study. Gambar II-1. Ilustrasi Tsunami

BAB II STUDI PUSTAKA

II.1 Rambatan Tsunami Gelombang tsunami terbentuk akibat adanya pergesaran vertikal massa air. Pergeseran ini bisa terjadi oleh gempa, letusan gunung berapi, runtuhan gunung es, dan meteor yang jatuh di laut.

Tsunami, Ilustration Zone 1

Dislocation

Deep sea

Has been studied widely

Propagated Shallow wave area water

Zone 2

Flooded area

Area of study

Gambar II-1. Ilustrasi Tsunami

Gelombang tsunami dapat dikategorikan sebagai gelombang di laut dangkal (d/L<0.05, d =

kedalaman laut, L = panjang gelombang) karena panjang

gelombangnya yang sangat panjang. Panjang gelombangnya mencapai ratusan kilometer dengan amplitudo ± 1 meter ketika merambat di laut dalam. Gelombang ini biasanya merambat dengan kecepatan ± 30-1000 km/jam, periode 5-90 menit.

II-1

Pemodelan Aliran Permukaan 2 D Pada Suatu Lahan Akibat Rambatan Tsunami

Gambar II-2. Panjang Gelombang dan Kecepatan Tsunami Terhadap Kedalaman (Hamzah,2006)

Hubungan antara cepat rambat gelombang dan kedalaman adalah: C=

gd ...( II-1 )

Dimana C: Cepat rambat gelombang g: Gravitasi d: Kedalaman Ketika memasuki daerah perairan dangkal, tsunami akan mengalami perlambatan. Akibatnya, gelombang yang berada di depan akan bertumpuk dengan gelombang yang dibelakangnya, sehinggat gelombang yang tadinya hanya memiliki amplitudo 1 meter, akan membesar hingga mencapai ± 30 meter. Rambatan gelombang tsunami di laut dangkal dapat dibagi ke dalam dua zona. Zona pertama adalah zona dimana gelombang yang terjadi diakibatkan oleh energi dari gelombang yang datang dari laut, sedangkan pada zona kedua adalah zona dimana gelombang yang terjadi akibat adanya gravitasi. Pada zona kedua ini, gelombang yang terjadi mirip dengan aliran air pada saat banjir. Pada kedua zona tersebut, gelombang dipengaruhi oleh friksi, kekasaran dasar.

II-2


II.2 Gelombang

Gambar II-3. Parameter Gelombang

Potongan memanjang perambatan gelombang pada kedalaman tetap di elevasi z=0 ditunjukkan pada gambar diatas. Dasar saluran berada pada koordinat z = -d, dengan kedalaman d. Profil gelombang (η) merupakan perubahan koordinat z terhadap jarak dan waktu. Profil gelombang memiliki cepat rambat sebesar C, dengan panjang gelombang L dan tinggi gelombang H. Periode gelombang dinyatakan dalam T, dimana dalam satu periode, gelombang akan menempuh jarak sebesar L. Hal ini berarti bahwa C=:L/T. Beberapa parameter gelombang: k = 2π / L (wave number) ...( II-2 ) σ = 2π / T (wave angular frequency) ...( II-3 ) Persamaan Laplace untuk aliran dua dimensi umum digunakan untuk menggambarkan pergerakan gelombang: ∂ 2 φ ∂ 2φ = 0 ...( II-4 ) + ∂x 2 ∂z 2

Syarat batas ditetapkan untuk dasar dan permukaan. Syarat batas di dasar adalah sebagai berikut:

II-3


w=

∂φ = 0 pada z = –d ...( II-5 ) ∂z

Syarat batas kinematik di permukaan: w=

∂η ∂η pada z = η +u ∂t ∂x

...( II-6 )

Syarat batas dinamik di permukaan: p ∂φ 1 2 u + w 2 + gz + + = 0 ...( II-7 ) ρ ∂t 2

(

)

Pada permukaan dimana tekanan tidak ada, maka syarat batas dinamik menjadi: ∂φ 1 2 u + w 2 + gη + = 0 pada z = η ...( II-8 ) 2 ∂t

(

)

Solusi dari persamaan laplace yang digunakan harus memenuhi syarat-syarat batas yang digunakan. Syarat batas di dasar memiliki bentuk linear, akan tetapi syarat batas kinematik dan dinamik di permukaan tidak. Dengan asumsi bahwa tinggi gelombang relatif kecil dibandingkan dengan kedalaman, maka syarat-syarat batas diatas dapat diterapkan pada level muka air normal. Kinematik: w=

∂η pada z = 0 ...( II-9 ) ∂t

Dinamik: gη +

∂φ = 0 pada z = 0 ...( II-10 ) ∂t

Solusi analitik persamaan laplace yang memenuhi syarat batas diatas adalah: Bentuk kecepatan potensial

φ=

gH cosh k (d + z ) sin (kx − σt ) ...( II-11 ) 2σ cosh kd

II-4


Profil permukaan Dengan memasukkan kecepatan potensial dengan syarat batas dinamik linier dan z = 0: H cos(kx − σt ) ...( II-12 ) 2

η= atau

η=

H ⎛x t⎞ cos 2π ⎜ − ⎟ ...( II-13 ) 2 ⎝L T⎠

Cepat Rambat Gelombang Dengan mengkombinasikan syarat batas dinamik dan kinematik: ∂ 2φ ∂φ +g = 0 pada z = 0 ...( II-14 ) 2 ∂z ∂t

dengan memasukkan kecepatan potensial, persamaan diatas dapat ditulis dalam bentuk:

σ 2 = gk tan kd ...( II-15 ) atau C=

σ k

=

g tanh kd ...( II-16 ) k

Dengan menggunakan hubungan C = L / T, persamaan di atas menjadi C=

gT 2πd ...( II-17 ) tanh 2π L

dan L=

gT 2 2πd ...( II-18 ) tanh 2π L

Teori gelombang diatas hanya berlaku jika kedalaman air relatif cukup dalam dan pengaruh dasar laut diabaikan. Persamaan aliran dinamik digunakan untuk kondisi kedalaman dimana efek dari dasar laut tidak dapat lagi diabaikan. II-5


II.2.1 Klasifikasi Gelombang Berdasarkan Kedalaman Relatif Pada saat geombang merambat dari laut dalam hingga ke laut dangkal, panjang gelombangnya akan berkurang. Akan tetapi, kedalaman air akan berubah lebih cepat sehingga, rasio kedalaman terhadap panjang gelombang (d/L) akan terus berkurang. Berdasarkan rasio ini, kedalaman dapat diklasifikan menjadi, laut dalam, laut peralihan dan laut dangkal. 1. d/L > 0.5, laut dalam 2. 0.05< d/L < 0.5, laut peralihan 3. d/L < 0.05, laut dangkal

Apabila kedalaman relatif untuk gelombang lebih besar dari 0.5 maka disebut laut dalam. Persamaannya adalah C0 =

gL0 2π

...( II-19 )

C0 =

gT 2π

...( II-20 )

L0 =

gT 2 2π

...( II-21 )

dan

Apabila kedalaman relatif kurang dari 0.5 maka disebut intermediate range. Persamaannya adalah C L 2πd = = tanh C 0 L0 L

...( II-22 )

Apabila kedalaman relatif kurang dari 0.05 maka disebut shallow water. Persamaannya adalah C = gd ...( II-23 )

II-6


sehingga kecepatan gelombang, dan hubungan antara perioda gelombang dengan panjang gelombang adalah L = CT atau L = gd T

...( II-24 )

II.2.2 Gelombang Pecah Gelombang pecah terjadi ketika kecepatan partikel gelombang mendekati atau sama dengan cepat rambat gelombangnya. Mekanisme keruntuhan gelombang sangat kompleks, mencakup interaksi kestabilan profil gelombang, puncak gelombang dan bentuk asimetris dari gelombang. Pada umumnya gelombang pecah dapat dikategorikan menjadi 4, yaitu: •

Spilling

•

Plunging

•

Collapsing

•

Surging

II-7


Gambar II-4. Profil Gelombang Pecah

Semua tipe gelombang pecah tersebut dapat terjadi di area shallow water. Tapi hanya spilling dan plunging yang muncul di laut dalam. Berdasarkan definisi dari Horikawa (1988), perbandingan tinggi gelombang dengan kedalaman air dapat dijadikan batasan untuk spilling, plunging dan surging. Spilling Æ H/d = 0.8 – 1 Plunging Æ H/d = 1 – 1.2 Collapsing Æ H/d > 1.2

II-8


II.3 Run Up Pada saat rambatan berada di zona dua aliran yang terjadi diakibatkan oleh run up gelombang datang. Tinggi run up maksimum akan bergantung kepada besarnya gelombang datang, kedalaman air normal dan kemiringan pantai. Solusi analitis untuk aproksimasi nonlinear besarnya tinggi Run Up untuk solitary wave (nonbreaking) diberikan oleh Synolakis (1986). 5/ 4

Rs ⎛H⎞ = 2.831 cot β ⎜ ⎟ ho ⎝ ho ⎠

...( II-25 )

Dengan Rs

= maksimum run up

h0

= kedalaman air normal

H

= tinggi gelombang

β

= sudut kemiringan shore line

Salah satu permasalahan dalam pemodelan run up adalah batasan wet/dry. Titiktitik grid di darat yang semula kering (h,u,v = 0) akan berubah basah akibat run up yang terjadi. Pemodelan batasan wet/dry dengan menggunakan perbandingan elevasi dasar di titik kering dengan elevasi air di titik sebelumnya diberikan oleh Synolakis (1986). Ying Li (2002) memodelkan batasan ini dengan menetapkan posisi shoreline sebagai fungsi waktu. Lynett (2002) mengambil nilai ekstrapolasi (h,u,v) dari titik titik di sebelumnya. Batasan wet/dry dengan menetapkan nilai batas kedalaman minimum diterapkan oleh Tawatchai Tingsanchal (1999) untuk pemodelan dambreak. Semua metode diatas telah dikomparasi dengan data eksperimen dan menunjukkan hasil yang baik.

II.4 Persamaan Gerak Aliran 2D Persamaan gerak aliran yang umum digunakan untuk memodelkan gelombang adalah persamaan St.Venant dan persamaan Boussinesq. Persamaan St.Venant:

II-9


∂h ∂ (uh ) ∂ (vh ) + + = S ...( II-26 ) ∂t ∂x ∂y ∂Uh ∂ 2 ∂ ∂ (h + z ) + U h + UVh + gh = − ghS fx ...( II-27 ) ∂t ∂x ∂y ∂x ∂Vh ∂ ∂ 2 ∂ (h + z ) + = −ghS fy ...( II-28 ) UVh + V h + gh ∂t ∂x ∂y ∂x

dimana : u dan v

= kecepatan arah x dan y

h

= kedalaman air

Sfx, Sfy

= kemiringan energi arah x dan y

Persamaan St. Venant tersebut pada umumnya berlaku di laut dangkal. Berdasarkan definisi dari Robert M. Sorensen (1993), tsunami merupakan gelombang panjang dengan d/L<0.05 (laut dangkal). Salah satu contoh penerapan persamaan St.Venant pada pemodelan rambatan gelombang tsunami adalah pemodelan oleh Vasily V Titov dan Costas Emmanuel Synolakis (1997) untuk kasus The Hokkaido-Nansei-Oki Tsunami dengan menggunakan persamaan pengatur St. Venant tanpa adanya friksi dasar (Sf=0) dan tidak membedakan zona rambatan. Pemodelan mencakup dari mulai terbentuknya gelombang hingga merambat di darat. Hasil pemodelan menunjukkan komparasi hasil yang cukup baik dalam skala besar. Persamaan gerak dalam bentuk lain yang juga dapat digunakan untuk memodelkan gelombang adalah persamaan Boussinesq. Perbedaan mendasar antara persamaan ini dengan persamaan St.Venant adalah adanya suku dispersi pada persamaan Boussinesq. Selain itu, persamaan boussinesq tidak memasukkan adanya pengaruh kemiringan dan kekasaran dasar. Aliran yang terjadi sangat dipengaruhi oleh kedalaman. Persamaan Boussinesq sendiri memiliki beberapa bentuk. Pada umumnya persamaan Boussinesq yang digunakan untuk pemodelan diberikan oleh Madsen Sorensen (1992) dan Nwogu (1993) Persamaan Boussinesq oleh Nwogu (1993)

II-10


⎧⎪⎛ zα 2 h2 ⎞ ⎫⎪ h⎞ ⎛ − ⎟ h∇ ( ∇.u ) + ⎜ zα + ⎟ h∇ ( ∇.(hu ) ) ⎬ = 0 ...( II-29 ) 6⎠ 2⎠ ⎝ ⎪⎩⎝ 2 ⎪⎭

ηt + ∇. ⎣⎡( h + η ) u ⎦⎤ + ∇ ⎨⎜

⎧z ⎫ ut + g∇η + ( u.∇ ) u + zα ⎨ α ∇ ( ∇.ut ) + ∇ ( ∇.(hut ) ) ⎬ = 0 ...( II-30 ) ⎩2 ⎭

Dimana: u = kecepatan (u,v) η = elevasi muka air h = kedalaman

∇ = (∂ / ∂x, ∂ / ∂y ) g = gravitasi Bentuk standar (Peregrine, 1967) dari persamaan tersebut diperoleh dengan memasukkan nilai Zα/h = -1/2 dan berlaku pada daerah laut dangkal. Nwogu (1993) memberikan Zα/h = 0.531 yang membuat persamaan ini berlaku untuk domain yang lebih luas, yaitu laut dalam hingga laut dangkal. Pemodelan rambatan gelombang dengan menggunakan persamaan ini telah dilakukan oleh Wei dan Kirby (1995) untuk kasus rambatan gelombang pada suatu saluran datar (tidak ada kemiringan dasar) dan memberikan komparasi yang sangat baik dengan data dari model fisik. Persamaan Boussinesq oleh Madsen Sorensen (1992) St + Px + Q y = 0 ...( II-31 )

⎛ P 2 ⎞ ⎛ PQ ⎞ ⎟ + Pt + ⎜ + gdS x + ψ 1 = 0 ...( II-32 ) ⎜ d ⎟ ⎜⎝ d ⎟⎠ y ⎝ ⎠x ⎛ Q 2 ⎞ ⎛ PQ ⎞ ⎟ + Qt + ⎜ + gdS y + ψ 2 = 0 ...( II-33 ) ⎜ d ⎟ ⎜⎝ d ⎟⎠ x ⎠y ⎝

ψ1 ≡

1 3 ⎡⎛ P ⎞ ⎛Q⎞ ⎤ 1 h ⎢⎜ ⎟ + ⎜ ⎟ ⎥ − h 2 Pxxt + Qxyt ...( II-34 ) 6 ⎢⎣⎝ h ⎠ xxt ⎝ h ⎠ xyt ⎥⎦ 2

(

)

II-11


ψ2 ≡

1 3 ⎡⎛ Q ⎞ ⎛P⎞ ⎤ 1 h ⎢⎜ ⎟ + ⎜ ⎟ ⎥ − h 2 Qyyt + Pxyt ...( II-35 ) 6 ⎢⎣⎝ h ⎠ yyt ⎝ h ⎠ xyt ⎥⎦ 2

(

)

dimana indeks x, y, dan t adalah simbol differensial untuk ruang dan waktu, d adalah kedalaman aliran total, h adalah kedalaman aliran, S adalah elevasi permukaan, P dan Q adalah komponen-komponen depth-integrated velocity, dan ψ1 dan ψ2 adalah bentuk-bentuk Boussinesq. Persamaan tersebut telah digunakan

untuk memodelkan transformasi gelombang akibat adanya struktur terendam oleh Nita Yunita (2001) dan menunjukkan komparasi yang baik dengan model serupa.

II.5 Metode Numerik Pemodelan numerik dapat dilakukan dengan pendekatan beda hingga (finite difference) ataupun volume hingga (finite Volume). Perbedaan paling mendasar dari kedua metode tersebut adalah bentuk grid. Pada metode beda hingga, bentuk grid terbatas kotak. Sedangkan pada metode volume hingga, bentuk grid lebih fleksibel (quadrangular). Penerapan metode beda hingga untuk aliran permukaan telah diterapkan oleh Dantje Kardana N,et.al., (2005) dan M. Syahril B.K, et al (2006) dan memberikan hasil yang baik. Pada penyelesaian persamaan diferensial numerik dengan metode beda hingga terdapat 3 macam skema yang dapat diterapkan, yaitu 1. Forward Difference Scheme Pada skema ini, nilai turunan pada suatu titik didekati dengan menggunakan nilai di titik tersebut dan nilai di titik sesudahnya. ∂f ( x) f ( x + Δx) − f ( x) = ...( II-36 ) ∂x Δx

Penyelesaian yang diperoleh dari skema ini memiliki ketelitian orde 1 2. Backward Difference Scheme Pada skema ini, nilai turunan pada suatu titik didekati dengan menggunakan nilai di titik tersebut dan nilai di titik sebelumnya. II-12


∂f ( x) f ( x) − f ( x − Δx) = ...( II-37 ) ∂x Δx

Penyelesaian yang diperoleh dari skema ini memiliki ketelitian orde 1 3. Central Difference Scheme Pada skema ini, nilai turunan pada suatu titik didekati dengan menggunakan nilai di titik tersebut dan nilai di titik sebelumnya. ∂f ( x) f ( x + Δx) − f ( x − Δx) = ...( II-38 ) ∂x Δx

Penyelesaian yang diperoleh dari skema ini memiliki ketelitian orde 2.

Gambar II-5. Pendekatan Numerik dengan Metode Beda Hingga

II-13

BAB II STUDI PUSTAKA. Propagated wave area. Shallow water. Area of study. Gambar II-1. Ilustrasi Tsunami

Recommend Documents