J. Math. and Its Appl. ISSN: 1829-605X Vol. 3, No. 2, Nov 2006, 67–80

## Pengaruh Amplitudo dan Frekuensi terhadap Fenomena Pemuncakan

## Wuryansari Muharini Kusumawinahyu
Jurusan Matematika FMIPA, Universitas Brawijaya Jl. Veteran Malang 65145, Indonesia.
e-mail: wmuharini@yahoo.com

## Abstrak
Telah diketahui bahwa suatu signal bikromatik yang diberikan seba- gai masukan di pembangkit gelombang di laboratorium hidrodinamika mengalami fenomena pemuncakan ketika ia merambat di sepanjang kolam pengujian. Dengan memanfaatkan besaran yang disebut Maximal Temporal Amplitude (MTA), dikaji pengaruh amplitudo dan frekuensi selubung signal gugus gelombang di pembangkit gelombang terhadap faktor amplifikasi dan posisi terjadinya elevasi maksimum. Data untuk penyelidikan ini dibangkitkan dengan menggunakan suatu simulasi numerik pembangkitan dan perambatan gelombang di sepanjang kolam pengujian dua dimensi. Hasil yang diperoleh dibandingkan dengan hasil yang diperoleh secara analitik menggunakan aproksimasi orde ke tiga persamaan KdV.
Kata Kunci: gelombang bikromatik, fenomena pemuncakan, amplitudo, frekuensi, Maximal Temporal Amplitude , faktor amplifikasi.

## 1. Pendahuluan
Gelombang bikromatik merupakan superposisi dua gelombang monokromatik de- ngan amplitudo dan frekuensi tertentu sehingga dapat dinyatakan sebagai gugus

## Pengaruh Amplitudo dan Frekuensi terhadap Fenomena Pemuncakan
gelombang dengan selubung yang khas. Berbagai penelitian mengenai pembang- kitan gelombang memberikan kesimpulan bahwa suatu signal bikromatik yang diberikan sebagai masukan di pembangkit gelombang di laboratorium hidrodina- mika akan mengalami perubahan bentuk dan elevasinya meninggi ketika ia me- rambat di sepanjang kolam pengujian. Lo dan Mei[11] adalah peneliti pertama yang melaporkan evolusi gugus gelombang bikromatik yang tak stabil dalam jang- ka waktu evolusi yang lama. Hasil pengukuran yang dilakukan oleh G.J. Keller dalam [11] memperlihatkan bahwa gelombang yang semula simetri menjadi tak simetri dan meningkat kecuramannya di bagian depan dan melandai di bagian belakang. Dengan menggunakan skema numerik, Lo dan Mei [11] juga menga- mati terjadinya pemisahan ( splitting ) gugus gelombang bikromatik menjadi dua gugus gelombang dengan kecepatan yang berlainan. Eksperimen serupa dilakukan pula oleh Stansberg pada tahun 1997. Hasil eksperimen Stansberg [13] memper- lihatkan bahwa gelombang permukaan yang pada pembangkit gelombang berupa signal bikromatik, dengan amplitudo q = 0 , 08 m dan frekuensi ω 1 = 3 , 3 rad/dt dan ω 2 = 2 , 99 rad/dt, dalam perambatannya mengalami deformasi selubung, berupa gejala peaking (pemuncakan) dan splitting .
Dalam penelitian ini dikaji pengaruh amplitudo dan frekuensi selubung signal gugus gelombang di pembangkit gelombang terhadap faktor amplifikasi dan posisi terjadinya elevasi maksimum, yang selanjutnya disebut posisi maksimum . Untuk itu diperkenalkan besaran yang dapat digunakan untuk menentukan faktor ampli- fikasi dan posisi maksimum secara global. Besaran yang dimaksud adalah Maximal Temporal Amplitude (MTA) yang merupakan nilai maksimum elevasi signal gelom- bang di setiap posisi. Data yang digunakan dalam penyelidikan ini dibangkitkan dengan menggunakan simulasi numerik pembangkitan dan perambatan gelombang di sepanjang kolam pengujian dua dimensi, yang disebut HUBRIS. Hasil yang di- peroleh secara numerik tersebut dibandingkan dengan hasil yang diperoleh secara analitik dengan menggunakan aproksimasi orde ke tiga persamaan KdV.
Makalah ini disusun mengikuti urutan berikut. Pertama diulas secara rinci mengenai gugus gelombang bikromatik untuk memperlihatkan pentingnya gugus gelombang tersebut. Kemudian dijelaskan mengenai persamaan lengkap dan mod- el numerik HUBRIS. Pada bagian berikutnya diberikan definisi mengenai Maximal Temporal Amplitude (MTA), posisi maksimum dan faktor amplifikasi . Secara anal- itik, perhitungan MTA dapat dilakukan dengan menggunakan hampiran orde ke tiga solusi persamaan KdV. Dengan menggunakan MTA diselidiki pengaruh ampli- tudo dan frekuensi selubung signal bikromatik di pembangkit gelombang terhadap posisi maksimum dan faktor amplifikasi. Hasil perhitungan dan pembahasannya diberikan pada bagian 5. Makalah ini diakhiri dengan kesimpulan pada bagian 6.

## 2. Gugus Gelombang Bikromatik
Gelombang bikromatik yang dikaji dalam makalah ini adalah gugus gelombang yang dihasilkan oleh superposisi dua gelombang monokromatik dengan frekuensi yang berbeda, yaitu
η ( x, t ) = q cos( k
1 x − ω 1 t ) + q cos( k 2 x − ω 2 t ) = 2 q cos( κx − νt ) cos( kx − ωt ) ,
dengan
k = k 1 + k 2 2 , κ = k 1 − k 2 2 , ω = ω 1 + ω 2 2 , dan ν = ω
1 − ω 2 2 .
Dari persamaan tersebut terlihat bahwa gugus gelombang bikromatik terdiri atas gelombang pembawa berfrekuensi ω yang dimodulasi oleh selubung gelom- bang berfrekuensi ν . Cepat rambat selubung gelombang tersebut adalah ν/κ dan disebut sebagai kecepatan gugus gelombang bikromatik.
Hasil menarik yang diperoleh Lo, Mei, dan Stansberg mendorong peneliti lain melakukan penyelidikan lebih mendalam terhadap fenomena deformasi gu- gus gelombang bikromatik. Penjelasan analitis terhadap sebagian fenomena yang ditemukan Stansberg diberikan oleh van Groesen dkk. [8]. Mereka mengguna- kan persamaan fase-amplitudo dari persamaan NLS untuk menjelaskan bahwa terjadinya fenomena tersebut disebabkan oleh adanya koreksi terhadap kecepatan fase dan kecepatan gugus yang tidak lagi hanya dipengaruhi oleh bilangan gelom- bang melainkan dipengaruhi pula oleh amplitudo. Westhuis dkk. [15] melaku- kan penyelidikan dengan membangkitkan gelombang bikromatik di laboratorium hidrodinamika menggunakan beberapa amplitudo dan frekuensi selubung yang berbeda-beda. Hasil eksperimen Westhuis dkk. [15] memperlihatkan bahwa defor- masi yang besar dengan peningkatan elevasi gelombang yang tinggi dipengaruhi oleh perbandingan antara amplitudo dan selisih frekuensi kedua gelombang mono- kromatik yang membentuk gelombang bikromatik. Salah satu eksperimen mereka menggunakan parameter yang sama dengan eksperimen Stansberg. Selain itu, dalam [15] dan [14] disajikan pula hasil rekonstruksi eksperimen Stansberg yang dilakukan secara numerik menggunakan HUBRIS.
Usaha untuk memahami fenomena deformasi gelombang bikromatik dilakukan pula oleh Cahyono [4] dengan membangun analytical wave code (AWC) yang mem- berikan koreksi terhadap bilangan gelombang. Penelitian yang dilakukan oleh Cahyono kemudian diperbaiki oleh Andonowati dkk. dalam [12, 3], yang dike- mukakan sekilas pada bagian 4.

## 3. Persamaan Lengkap dan Model Numerik
Perhatikan gelombang yang merambat pada suatu lapisan fluida dengan kedala- man h di atas dasar rata pada suatu kolam. Misalkan fluida tersebut diasumsikan homogen, tak kental, tak termampatkan dan tak terotasi. Bila diasumsikan gelom- bang merambat tanpa batas hanya dalam arah horisontal searah dengan sumbu x dan tekanan permukaan diabaikan, maka potensial kecepatan fluida φ ( x, z, t ) adalah solusi persamaan Laplace
∂ 2 φ ∂x 2 + ∂ 2 φ ∂z 2 = 0 , − h ≤ z ≤ η ( x, t ) . (1)
Persamaan tersebut dilengkapi dengan syarat-syarat batas berikut:
∂φ ∂z = 0 , z = − h, (2)
∂η ∂t + ∂φ ∂x ∂η ∂x = ∂φ ∂z , z = η ( x, t ) , (3)
gη + ∂φ ∂t + 1 2 "µ ∂φ ∂x ¶ 2 + µ ∂φ ∂z ¶ 2 # = 0 , z = η ( x, t ) , (4)
dengan η ( x, t ) menyatakan elevasi permukaan fluida yang diasumsikan bernilai tunggal. Persamaan (1) hingga (4) dikenal sebagai persamaan lengkap ( full equa- tion ), sebagaimana diungkapkan dalam buku-buku teks baku tentang gelombang seperti [1, 5, 6].
Solusi eksak persamaan lengkap tersebut tidak mudah ditentukan, karena syarat batas kinematik dan dinamik melibatkan potensial kecepatan φ ( x, t, z ) dan elevasi permukaan η ( x, t ) secara simultan dalam bentuk tak linier. Selain itu, persamaan lengkap (1)–(4) belum melibatkan syarat batas pada arah horisontal karena arah horisontal dianggap tak terbatas. Model matematika gerakan air di ko- lam hidrodinamika memerlukan syarat batas pada sisi pembangkit gelombang dan pada sisi lain yang menyerupai pantai dan berfungsi sebagai penyerap agar tidak terjadi gelombang pantul. Selama ini belum diperoleh solusi analitik eksak yang mendeskripsikan dinamika gelombang permukaan air di sepanjang kolam hidrodi- namika karena rumitnya formulasi syarat batas. Berdasarkan persamaan Laplace yang dilengkapi dengan syarat batas di pembangkit gelombang dan sisi penye- rap, Westhuis [14] telah membangun suatu model numerik yang disebut HUBRIS. Model numerik tersebut mensimulasikan pembangkitan dan perambatan dalam satu arah gelombang permukaan bebas tak linier dua dimensi yang deterministik di kolam hidrodinamika. HUBRIS telah diuji dengan percobaan di laboratorium hidrodinamika Maritime Research Institute Netherlands (MARIN), Wageningen, Belanda, dan menunjukkan kesesuaian yang cukup akurat. Bila pada model nu- merik tersebut diberikan masukan berupa signal η ( t ) pada pembangkit gelombang,
maka akan diperoleh keluaran berupa data elevasi gelombang permukaan air η ( x, t ) dan potensial kecepatan air φ ( x, z, t ), dengan ( x, z, t ) menyatakan koordinat ruang dan waktu fisis di laboratorium hidrodinamika. Model numerik HUBRIS yang di- gunakan dalam penelitian ini adalah HUBRIS versi 2.0 dengan masukan berupa signal polikromatik tanpa tambahan fase
η ( t ) = m X j =1 a j e i ( ω j t ) + c.c.,
dengan 1 ≤ m ≤ 10 dan c.c. menyatakan complex conjugate (kompleks sekawan) suku sebelumnya.

## 4. Maximal Temporal Amplitude (MTA)
Untuk mengamati perilaku global elevasi gelombang dalam perambatannya di sepanjang kolam pengujian, digunakan suatu parameter yang disebut Maximal Temporal Amplitude (MTA). Parameter yang menyatakan nilai maksimum elevasi signal gelombang di setiap posisi tersebut telah digunakan dalam [2] untuk menye- lidiki deformasi gelombang optik. Untuk masalah perambatan gelombang air di sepanjang kolam hidrodinamika, MTA yang didefinisikan sebagai
M ( x ) = max t η ( x, t ) (5)
dapat diartikan sebagai tinggi bagian dinding kolam pada sisi panjang yang ter- basahi oleh air selama gelombang merambat, bila diukur dari permukaan air ten- ang z = 0. Untuk mengetahui besarnya perubahan elevasi gelombang secara glo- bal, diamati perbandingan antara nilai MTA tertinggi dan nilai MTA di posisi pembangkit gelombang. Perbandingan itu disebut Faktor Amplifikasi Amplitudo (FAA) yang didefinisikan sebagai
F AA = max M ( x ) M (0) . (6)
Posisi maksimum yang dinotasikan sebagai x max dapat ditentukan dengan menggunakan MTA dengan menentukan posisi di mana pertama kali MTA men- capai maksimum. Nilai elevasi η ( x, t ) dan tinggi selubung A ( x, t ) yang digunakan untuk menghitung MTA dapat diperoleh secara analitik atau numerik.
Secara analitik, Andonowati dan Marwan dalam [3] merumuskan nilai η ( x, t ) untuk gelombang bikromatik dan gelombang trikromatik melalui aproksimasi orde ke tiga solusi persamaan KdV, yang selanjutnya disingkat dengan AOT-KdV. Persamaan KdV merupakan hampiran persamaan lengkap yang memodelkan eleva- si gelombang permukaan yang cukup rendah dan cukup panjang, serta merambat

## Pengaruh Amplitudo dan Frekuensi terhadap Fenomena Pemuncakan
dalam satu arah. Persamaan KdV diperoleh melalui proses uni-directionalisation . Persamaan KdV telah diperbaiki oleh van Groesen [7] menjadi persamaan berikut.
∂η ∂t + i Ω( − i ∂ ∂x ) η = − C 0 η ∂η ∂x , (7)
dengan Ω( k ) = p k tanh( k ) adalah relasi dispersi linier yang mengaitkan frekuensi ω dengan bilangan gelombang k dalam koordinat yang dinormalkan, C 0 suatu kon- stanta yang menyatakan koefisien ketaklinieran, dan η ( x, t ) adalah elevasi gelom- bang dalam koordinat yang dinormalkan. Kaitan antara koordinat yang dinor- malkan dengan koordinat fisis adalah
η lab = hη, x lab = hx, t lab = s h g t, k lab = k h , ω lab = r g h ω. (8)
Solusi η ( x, t ) yang diperoleh melalui AOT-KdV adalah solusi η ( x, t ) yang berbentuk
η ( x, t ) = εη (1) ( x, t ) + ε 2 η (2) ( x, t ) + ε 3 η (3) ( x, t ) + . . . , (9)
dengan ε adalah suatu parameter positif yang kecil yang menyatakan orde ampli- tudo elevasi. Pada orde pertama diambil solusi
η (1) ( x, t ) = m X p =1 a p e i ( k p x − ω p t + ψ p ) + c.c.,
dengan ω p = Ω( k p ). Untuk m = 2 dan m = 3, maka berturut-turut dipero- leh gugus gelombang bikromatik dan gelombang trikromatik. Untuk menghindari resonansi, dilakukan koreksi terhadap bilangan gelombang, dengan menyatakan bi- langan gelombang sebagai k nl p = k p + εk (1) p + ε 2 k (2) p + O ( ε 3 ) , dengan k p = Ω − 1 ( ω p ). Substitusi persamaan (9) ke persamaan KdV (7) menghasilkan polinom dalam ε . Dengan mengumpulkan koefisien polinom tersebut di setiap orde ε diperoleh persa- maan diferensial linier di setiap orde. Solusi persamaan diferensial pada orde yang lebih tinggi diperoleh dengan menggunakan solusi persamaan diferensial pada orde yang lebih rendah. Dengan cara tersebut diperoleh solusi
η ( x, t ) = η (1) ( x, t ) + η (2) ( x, t ) − η (2) f w ( x, t ) + η (3) ( x, t ) − η (3) f w ( x, t ) ,
dengan η ( j ) f w ( x, t ) , j = 2 , 3 menyatakan solusi free wave , yaitu solusi yang perlu ditambahkan untuk mengkompensasi munculnya solusi orde ke dua dan ke tiga di pembangkit gelombang.
Untuk kasus gelombang bikromatik, dalam [3] diambil gelombang bikromatik khusus dengan a 1 = a 2 = q sebagai solusi orde pertama, yaitu
η ( x, t ) (1) = q cos( k nl 1 x − ω 1 t + ψ 1 ) + q cos( k nl 2 x − ω 2 t + ψ 2 )
Dengan menggunakan solusi tersebut dalam definisi MTA (5), Andonowati dan Marwan dalam [3] merumuskan x max untuk gelombang bikromatik tersebut, se- bagai
x max = π ¯ ¯ k − Ω ¯ ¯ = π 2 ¯ ¯ ¯ ν 2 K 00 ( ω ) + q 2 ³ ˜ γ + C 2 0 σ 2 K ( ω ) Ω 0 ( K ( ω )) ´¯ ¯ ¯ = O µ 1 ν 2 , 1 q 2 ¶ , (10)
dengan k = ( k nl 1 + k nl 2 ) / 2, κ = | k nl 1 − k nl 2 | / 2, ω = ( ω 1 + ω 2 ) / 2, ν = | ω 1 − ω 2 | / 2, Ω = ( K ( ω + 3 ν ) + K ( ω − 3 ν )) / 2, ˜ γ = C 2 0 k ( σ 0 + σ 2 ), σ 0 = 1 Ω 0 ( K ( ω )) − Ω 0 (0) ,
σ 2 = K ( ω ) 2 ω − Ω(2 K ( ω )) dan K 00 menyatakan turunan ke dua K . Dari hasil terse- but dapat disimpulkan bahwa nilai x max tidak dipengaruhi oleh tambahan fase ψ 1 maupun ψ 2 melainkan hanya dipengaruhi oleh amplitudo, q , dan frekuensi selubung, ν .

## 5. Hasil dan Pembahasan

## 5.1. Pengaruh Amplitudo
Nilai η ( x, t ) diperoleh dari keluaran program paket HUBRIS yang mensimulasikan pembangkitan dan perambatan gelombang di kolam pengujian. Oleh sebab itu koordinat yang digunakan adalah koordinat fisis, dengan kedalaman kolam h = 5 m dan percepatan gravitasi g = 9 , 81 m/dt 2 . Semua nilai diberikan dalam koordinat fisis dengan sistem satuan ( m, kg, dt ). Masukan yang diberikan di pembangkit gelombang adalah signal
η (0 , t ) = q cos( ω 1 t ) + q cos( ω 2 t ) = 2 q cos( νt ) cos( ωt ) .
Kebergantungan η ( x, t ) terhadap rasio antara q dan ν telah dikaji dalam [15, 4, 9], namun belum diketahui kebergantungan perilaku MTA, M ( x ), terhadap kedua parameter tersebut. Untuk menyelidiki pengaruh amplitudo dan frekuen- si selubung gelombang bikromatik terhadap F AA dan x max , digunakan beberapa nilai amplitudo q dan frekuensi selubung ν signal bikromatik di pembangkit gelom- bang, dengan ω = 3 , 145 / dt.
Pertama-tama, untuk memusatkan perhatian pada pengaruh amplitudo, fre- kuensi selubung diambil konstan, yaitu ν = 0 , 155, sedangkan amplitudo diambil bervariasi, yaitu 0 , 01 ≤ q ≤ 0 , 09. Untuk mengilustrasikan perubahan x max dan F AA terhadap q , dalam Gambar 1 diperlihatkan kurva MTA evolusi gelom- bang air di permukaan, yang di posisi x = 0 berupa signal bikromatik dengan 0 , 01 ≤ q ≤ 0 , 03. Dalam [8] gelombang bikromatik tersebut digolongkan sebagai gelombang rendah. Gelombang-gelombang tersebut belum memperlihatkan gejala

## Pengaruh Amplitudo dan Frekuensi terhadap Fenomena Pemuncakan
Gambar 1: Grafik nilai MTA, M ( x ), gelombang bikromatik dengan ν = 0 , 155 konstan dan untuk q = 0 , 01, q = 0 , 02, dan q = 0 , 03.
Gambar 2: Grafik nilai MTA, M ( x ), gelombang bikromatik dengan ν = 0 , 155 konstan dan untuk q = 0 , 04, q = 0 , 05, dan q = 0 , 06.
deformasi spatial yang mudah teramati. Oleh karena itu grafik MTA ketiganya menunjukkan pola yang hampir sama.
Gejala deformasi spatial gelombang bikromatik mulai terlihat untuk 0 , 04 ≤ q ≤ 0 , 06. Gelombang-gelombang bikromatik dengan 0 , 04 ≤ q ≤ 0 , 06 digolongkan sebagai gelombang moderat [8]. Grafik MTA ketiga gelombang moderat tersebut memperlihatkan F AA yang meningkat, seperti diperlihatkan dalam Gambar 2. Perubahan F AA dan x max yang signifikan terjadi untuk gelombang bikromatik dengan 0 , 07 ≤ q ≤ 0 , 09 seperti diperlihatkan dalam Gambar 3. Menurut [8], gelombang yang dapat mencapai elevasi setinggi 5% kedalaman air digolongkan sebagai gelombang ekstrim. Dalam gambar tersebut terlihat bahwa tinggi MTA untuk q = 0 , 07 mencapai 0 , 25 m, yaitu 5% dari h = 5 m. Dengan demikian gelom- bang bikromatik dengan q ≥ 0 , 07 dapat digolongkan sebagai gelombang ekstrim. Amplitudo sebesar q = 0 , 08 adalah amplitudo yang digunakan oleh Stansberg
Gambar 3: Grafik nilai MTA, M ( x ), gelombang bikromatik dengan ν = 0 , 155 konstan dan untuk q = 0 , 07, q = 0 , 08, dan q = 0 , 09.
dalam eksperimennya.
Gejala deformasi yang tidak signifikan pada kasus gelombang rendah dan gelombang moderat menimbulkan masalah dalam penentuan x max . Dalam Gam- bar 1 dan 2 terlihat bahwa nilai MTA mencapai nilai maksimum yang hampir sama di beberapa posisi, sehingga tidak mudah menentukan posisi mana yang meru- pakan posisi maksimum yang sebenarnya. Untuk mengatasi masalah tersebut, dilakukan proses averaging (perata-rataan) dengan terlebih dahulu menentukan nilai maksimum dan minimum lokal M ( x ) di setiap posisi. Kemudian nilai rata- rata diperoleh dengan mengambil rata-rata nilai maksimum dan minimum lokal di setiap posisi. Nilai rata-rata yang diperoleh di sini disebut sebagai hasil averaging I. Bila hasil averaging I belum dapat memberikan informasi yang meyakinkan, ma- ka dilakukan proses averaging sekali lagi terhadap hasil averaging I atau diambil nilai rata-rata dua posisi maksimum yang berdekatan yang diperoleh dari hasil averaging I.
Hasil perhitungan x max yang telah melalui proses averaging dan hasil perhi- tungan F AA , dibandingkan dengan hasil yang diperoleh secara analitik dengan menggunakan AOT-KdV, disajikan dalam Tabel 1.
Pada kolom ke dua dan ke tiga Tabel tersebut diperlihatkan perubahan x max bila q berubah, berturut-turut dengan menggunakan data elevasi η ( x, t ) yang di- peroleh dari simulasi numerik HUBRIS dan yang diperoleh dari AOT-KdV. Terli- hat bahwa nilai x max mengecil bila q bertambah. Untuk melihat pola kebergantun- gan x max terhadap q , grafik x max – q yang diperoleh secara numerik dibandingkan dengan grafik x max – q yang diperoleh secara analitik dari AOT-KdV, seperti dis- ajikan dalam Gambar 4 (a). Karena diperoleh pola yang hampir sama, dapat dikatakan bahwa x max = O ¡ 1 /q 2 ¢ .
Dari kolom ke empat dan ke lima Tabel 1 dapat dilihat bahwa untuk q ≥ 0 , 03, F AA meningkat bila nilai q membesar. Untuk melihat pola peningkatan F AA

## Pengaruh Amplitudo dan Frekuensi terhadap Fenomena Pemuncakan
Tabel 1: Pengaruh beda frekuensi dan amplitudo selubung signal di pembangk- it gelombang terhadap Faktor Amplifikasi Amplitudo ( F AA ) dan posisi maksimum MTA gelombang bikromatik.
ν = 0 , 155 q = 0 , 07 x max F AA x max F AA q HBR AOT HBR AOT ν HBR AOT HBR AOT 0 , 010 162 , 6 156 , 5 1 , 238 1 , 025 0 , 6200 10 , 2 9 , 2 1 , 062 1 , 149 0 , 020 158 , 5 154 , 1 1 , 191 1 , 050 0 , 4650 18 , 2 16 , 3 1 , 063 1 , 157 0 , 030 149 , 9 150 , 3 1 , 190 1 , 083 0 , 3100 40 , 6 36 , 4 1 , 124 1 , 165 0 , 040 141 , 6 145 , 3 1 , 240 1 , 121 0 , 2325 74 , 2 62 , 5 1 , 204 1 , 188 0 , 050 139 , 3 139 , 3 1 , 333 1 , 161 0 , 1938 101 , 7 86 , 4 1 , 327 1 , 212 0 , 055 137 , 2 136 , 0 1 , 379 1 , 187 0 , 1550 130 , 2 125 , 5 1 , 591 1 , 269 0 , 060 135 , 0 132 , 6 1 , 452 1 , 214 0 , 1162 184 , 8 193 , 1 2 , 005 1 , 361 0 , 065 133 , 4 129 , 1 1 , 527 1 , 242 0 , 0969 205 , 4 244 , 3 2 , 408 1 , 462 0 , 070 130 , 2 125 , 5 1 , 591 1 , 269 0 , 0775 223 , 6 312 , 2 2 , 863 1 , 646 0 , 075 127 , 6 121 , 9 1 , 694 1 , 293 0 , 080 126 , 0 118 , 2 1 , 806 1 , 315 0 , 085 123 , 0 114 , 6 1 , 879 1 , 344 0 , 090 111 , 0 110 , 9 1 , 975 1 , 379
terhadap q , dalam Gambar 4 (b) diperlihatkan kurva F AA terhadap q yang di- peroleh secara numerik dengan menggunakan program paket HUBRIS dan yang diperoleh secara analitik dengan menggunakan AOT-KdV. F AA yang diperoleh melalui AOT-KdV lebih rendah daripada F AA yang diperoleh secara numerik, se- bab solusi yang diperoleh dari AOT-KdV menghasilkan interaksi tak linier antar mode yang terbatas hanya sampai orde ke tiga, sedangkan solusi yang diperoleh dari program paket HUBRIS menghasilkan interaksi tak linier antar mode hing- ga orde tinggi. Dalam [8, 15, 4, 9] diungkapkan bahwa sifat ketaklinieran, yang menyebabkan meningkatnya elevasi, meningkat bersama meningkatnya q . Oleh karena itu dapat dipahami bahwa semakin tinggi q maka semakin cepat gelom- bang meninggi, sehingga posisi maksimum semakin dekat dan faktor amplifikasi juga semakin besar.

## 5.2. Pengaruh Frekuensi
Untuk melihat pengaruh frekuensi selubung gelombang bikromatik, ν , terhadap x max dan F AA , amplitudo q diambil konstan sedangkan ν diambil bervariasi dari 0 , 0775 sampai dengan 0 , 62. Dari Tabel 1 dapat diperiksa bahwa x max berkurang sekitar 4 kali lipat bila nilai ν bertambah menjadi dua kali lipat, kecuali untuk
(a)
(b)
Gambar 4: (a) Grafik nilai x max MTA terhadap q dan (b) grafik F AA terhadap q yang diperoleh dari program paket HUBRIS (hitam) dibandingkan dengan hasil yang diperoleh dari AOT - KdV (merah), untuk gelombang bikromatik dengan nilai ν = 0 , 155 dan ω = 3 , 145 konstan.
ν ≤ 0 , 155. Dengan demikian, untuk kasus ini, dapat dikatakan bahwa posisi mak- simum MTA mendekati O (1 /ν 2 ). Hasil ini sesuai dengan posisi maksimum yang diperoleh secara analitik dengan menggunakan AOT - KdV [3], yang diberikan pada persamaan (10). Perbandingan hasil perhitungan x max secara numerik dan analitik tersebut disajikan secara grafis dalam Gambar 5 (a). Dalam Tabel 1 dan Gambar 5 (b) diperlihatkan bahwa nilai F AA berkurang bila ν bertambah. De- ngan demikian untuk memperoleh faktor amplifikasi yang tinggi sebaiknya diambil nilai ν yang kecil. Namun, nilai ν kecil yang menghasilkan amplifikasi tinggi ju- ga mengakibatkan peningkatan x max , sehingga bila diambil ν yang terlalu kecil gelombang akan mencapai maksimum di posisi yang jaraknya dari pembangkit gelombang dapat melampaui panjang kolam pengujian. Nilai ν = 0 , 155 mem- berikan nilai F AA yang cukup tinggi dan posisi maksimum yang tidak terlalu jauh.

## 6. Kesimpulan
Telah dilakukan penyelidikan mengenai pengaruh amplitudo dan frekuensi selu- bung gugus gelombang bikromatik terhadap faktor amplifikasi dan posisi maksi- mum . Dari pengamatan nilai MTA dapat disimpulkan hal-hal berikut.
Pengaruh amplitudo dan frekuensi selubung terhadap posisi maksimum yang diperoleh melalui penyelidikan yang dilakukan secara numerik mendukung hasil yang diperoleh secara analitik dengan menggunakan aproksimasi orde ke tiga persamaan KdV. Posisi maksimum makin dekat ke posisi pembangkit gelombang bila amplitudo maupun frekuensi selubung membesar. Dapat disimpulkan, bahwa

## Pengaruh Amplitudo dan Frekuensi terhadap Fenomena Pemuncakan
(a) (b)
Gambar 5: (a) Grafik nilai x max MTA terhadap ν dan (b) grafik F AA terhadap ν yang diperoleh dari program paket HUBRIS (hitam), dibandingkan dengan hasil yang diperoleh dari AOT - KdV (merah), untuk gelombang bikromatik dengan q = 0 , 07 konstan.
posisi maksimum MTA gelombang bikromatik berbanding terbalik dengan kuadrat amplitudo selubung signal di pembangkit gelombang. Di samping itu, posisi mak- simum MTA gelombang bikromatik juga berbanding terbalik dengan kuadrat fre- kuensi selubung signal di pembangkit gelombang, asalkan nilai frekuensi selubung tidak terlalu kecil.
Nilai Faktor Amplifikasi Amplitudo (FAA) evolusi gelombang bikromatik di- pengaruhi oleh besar amplitudo dan frekuensi selubung gugus signal bikromatik di pembangkit gelombang. Nilai FAA meningkat bila nilai amplitudo membesar. Sebaliknya, nilai FAA meningkat bila nilai frekuensi selubung mengecil. Hasil tersebut memberikan kecenderungan yang sama dengan hasil yang diperoleh se- cara analitik yang menggunakan aproksimasi orde ke tiga persamaan KdV dalam [3].

## Ucapan terima kasih
Penulis mengucapkan terima kasih kepada Dr. Jaap-harm Westhuis atas diijinkan- nya penulis menggunakan program paket HUBRIS dalam penelitian ini.

## Pustaka
[1] Akhmediev dan A. Ankiewicz (1997), Solitons, nonlinear pulses and beams , Chapman & Hall.
[2] Andonowati, dan van Groesen, E. (2003), Optical pulse deformation in second order non-linear media, Journal of Nonlinear Optical Physics & Materials , 12 , 221 – 234.
[3] Andonowati, Marwan, dan Kusumawinahyu, W. M. (2005), Phase effect on the Maximal Temporal Amplitude of down-stream running nonlinear water waves, dikirim ke ANZIAM Journal .
[4] Cahyono, E. (2002), Analitical wave codes for predicting surface waves in a laboratory basin , PhD Thesis, University of Twente, the Netherlands.
[5] Debnath, L. (1994), Nonlinear Water Waves , Academic Press, San Diego.
[6] Dingemans, M. W. (1997), Water wave propagation over uneven bottoms, Part 1: Linear Wave Propagation , volume 13 of Advanced series on ocean engineering, World Scientific, Singapore.
[7] van Groesen, E. (1998), Wave groups in uni-directional surface wave models, J. Eng. Math. , 34 , 215– 226.
[8] van Groesen, E., Andonowati, dan Soewono, E. (1999), Non-linear effects in bi-chromatic surface waves, Proc. Estonian Acad. Sci. Phys. Math. , 48 , 206 – 229.
[9] van Groesen, E. W. C. dan Westhuis, J. H. (2002), Modelling and simulation of surface water waves, Mathematics and computers in simulation, 59(4) , 341 – 360.
[10] Korteweg, D. J. dan de Vries, G. (1895), On the change of form of long waves advancing in a rectangular canal, and on a new type of long stationary waves, Phil. Mag. and J. of Science , 39 , 422 – 443.
[11] Lo, E. dan Mei, C. C. (1985), A numerical study of water-wave modulation based on a higher order nonlinear Schr˝odinger equation, J. Fluid Mech. , 150 , 395 – 416.
[12] Marwan dan Andonowati (2003), Wave deformation on the propagation of bi-chromatics signal and its effect to the maximum amplitude, JMS FMIPA ITB , vol. 8 no. 2 , 81–87.

## Pengaruh Amplitudo dan Frekuensi terhadap Fenomena Pemuncakan
[13] Stansberg, C. T. (1997), On the nonlinear behaviour of ocean wave groups, Proc. Waves ’97 , 2 , 1227.
[14] Westhuis, J. H. (2001), The numerical simulation of nonlinear waves in a hydrodynamic model test basin , PhD thesis, University of Twente, the Netherlands.
[15] Westhuis, J. H., van Groesen, E., dan Huijsmans, R. H. M. (2001), Experi- ments and numerics of bichromatic wave groups, J. Waterway, Port, Coastal and Ocean Engineering , 127 , 334 – 342.