26.07.2013 Views

Kort introduksjon til bølger relatert til tidevann i havet1 1 Bølger i hav ...

Kort introduksjon til bølger relatert til tidevann i havet1 1 Bølger i hav ...

Kort introduksjon til bølger relatert til tidevann i havet1 1 Bølger i hav ...

SHOW MORE
SHOW LESS

Transform your PDFs into Flipbooks and boost your revenue!

Leverage SEO-optimized Flipbooks, powerful backlinks, and multimedia content to professionally showcase your products and significantly increase your reach.

<strong>Kort</strong> <strong>introduksjon</strong> <strong>til</strong> <strong>bølger</strong> <strong>relatert</strong> <strong>til</strong> <strong>tidevann</strong> i <strong>hav</strong>et 1<br />

Helge Drange<br />

Høst 2012<br />

Kommentarer <strong>til</strong> helge.drange@gfi.uib.no<br />

Oppdatert 24. oktober 2012<br />

Observert vannstand (cm) i Bergen mellom 1. juni og 31. august 2012 (svart farge) og modellert<br />

vannstand basert p˚a <strong>tidevann</strong>skomponentene M2, S2, N2 og K2 (rød farge). Observasjoner fra<br />

http://vannstand.no/.<br />

1 <strong>Bølger</strong> i <strong>hav</strong> og atmosfære<br />

<strong>Bølger</strong> betegner fysiske prosesser som transporterer informasjon i tid og rom, som energi, uten<br />

eller med liten adveksjon av masse knyttet <strong>til</strong> transporten. Dette i motsetning <strong>til</strong> for eksempel<br />

geostrofisk balanse og ageostrofisk strøm i atmosfære og <strong>hav</strong> som alltid er knyttet <strong>til</strong> adveksjon<br />

av masse. Følgelig kan bølgesignal forplantes uten eller med liten p˚avirkning av jordrotasjonen,<br />

selv om forplantningshastigheten kan være svært stor. Dette i motsetning <strong>til</strong> adveksjon av masse<br />

som alltid vil være p˚avirket av Corioliseffekten og som skalerer lineært med væskeelementets<br />

hastighet.<br />

En generell gjennomgang av sentrale egenskaper <strong>til</strong> <strong>bølger</strong> er gitt i appendiks.<br />

1 Teorien er i hovedsak hentet fra INTRODUCTION TO GEOPHYSICAL FLUID DYNAMICS, Physical and<br />

Numerical Aspects (2010): Benoit Cushman-Roisin og Jean-Marie Beckers, Academic Press (Publication date:<br />

September 2010), <strong>til</strong>gjengelig fra http://engineering.dartmouth.edu/∼cushman/books/GFD.html.<br />

1


2 Utledning av grunntvannsligningene<br />

2.1 Utgangspunkt og konfigurering<br />

Utgangspunktet for grunntvannsligningene er standard form av horisontal momentumligning<br />

∂u<br />

∂u<br />

+ u∂u + v<br />

∂t ∂x ∂y<br />

∂v ∂v ∂v<br />

+ u + v<br />

∂t ∂x ∂y<br />

og kontinuitetsligningen (hvor u = (u, v))<br />

1 ∂p<br />

− fv = − + Fx<br />

(1)<br />

ρ ∂x<br />

1 ∂p<br />

+ fu = − + Fy<br />

(2)<br />

ρ ∂y<br />

∂ρ<br />

+ ∇ · (ρ u) = 0 (3)<br />

∂t<br />

Vi betrakter en homogen (ρ = ρ0), friksjonsfri (F H=0) og barotrop (∂u/∂z = ∂v/∂z = 0)<br />

væske med fri overflate η(x, y, t) som illustrert i figur 1. Ligningene (1)-(3) kan da skrives p˚a<br />

Figur 1: Illustrasjon p˚a en homogen væske med fri overflate η(x, y, t) og generell batymetri b(x, y).<br />

h(x, y, t) beskriver total væskehyde. zref = H er et referanseniv˚a som beskriver overflaten uten<br />

<strong>bølger</strong>.<br />

formen<br />

∂u<br />

∂u 1 ∂p<br />

+ u∂u + v − fv = −<br />

∂t ∂x ∂y ρ0 ∂x<br />

∂v ∂v ∂v<br />

1 ∂p<br />

+ u + v + fu = −<br />

∂t ∂x ∂y ρ0 ∂y<br />

(4)<br />

(5)<br />

∂u ∂v ∂w<br />

+ + = 0 (6)<br />

∂x ∂y ∂z<br />

2


2.2 Omskrevet form av kontinuitetsligningen<br />

Kontinuitetslignngen er gitt ved (6), der u og v er uavhengig av z (grunnet antagelsen om<br />

barotrop væske), men ∂u/∂x = 0 og ∂v/∂y = 0, slik at vi kan ha divergent strøm.<br />

Vi betrakter den homogene væsken som vist i figur 1 og integrerer (6) fra bunnen z = b(x, y) <strong>til</strong><br />

den frie overflaten z = b(x, y) + h(x, y, t). Dette gir<br />

b+h<br />

∂u ∂v<br />

+ dz + w|<br />

∂x ∂y b<br />

b+h<br />

b = 0 (7)<br />

Her er w(z = b + h) bevegelsen <strong>til</strong> den frie overflaten. Denne kan uttrykkes ved hjelp av den<br />

totalderiverte (som beskriver hvordan overflaten endrer seg med bevegelsen)<br />

w(z = b + h) = Dz<br />

<br />

<br />

=<br />

Dt<br />

D<br />

(b + h) (8)<br />

Dt<br />

P˚a <strong>til</strong>svarende m˚ate er<br />

b+h<br />

= ∂<br />

∂<br />

∂<br />

(b + h) + u (b + h) + v (b + h)<br />

∂t ∂x ∂y<br />

(9)<br />

= ∂h ∂<br />

∂<br />

+ u (b + h) + v (b + h)<br />

∂t ∂x ∂y<br />

(10)<br />

w(z = b) = Db ∂b ∂b<br />

= u + v<br />

Dt ∂x ∂y<br />

Uttrykkene (10) og (11) innsatt i (7) gir<br />

<br />

∂u ∂v<br />

+ (b + h − b) +<br />

∂x ∂y<br />

∂h ∂<br />

∂<br />

∂b ∂b<br />

+ u (b + h) + v (b + h) − u − v = 0 (12)<br />

∂t ∂x ∂y ∂x ∂y<br />

eller<br />

(11)<br />

∂h ∂ ∂<br />

+ (hu) + (hv) = 0 (13)<br />

∂t ∂x ∂y<br />

som er kontinuitetsligningen uttrykt ved lokal endring av væskesøylen h og divergensen <strong>til</strong><br />

huH.<br />

Alternativt, siden<br />

kan den tidsderiverte i (13) uttrykkes ved η<br />

h(x, y, t) + b(x, y) = H + η(x, y, t) (14)<br />

2.3 Omskrevet form av momentumligningene<br />

De horisontale momentumligningene er gitt ved (4) og (5).<br />

Trykket p(x, y, z, t) er gitt ved hydrostatisk ligning<br />

∂η ∂ ∂<br />

+ (hu) + (hv) = 0 (15)<br />

∂t ∂x ∂y<br />

∂p<br />

= −gρ0<br />

(16)<br />

∂z<br />

3


Trykkvariasjoner som skyldes endringer i overflatehevningen kan da uttrykkes som<br />

(17) kan skrives<br />

eller<br />

P˚a <strong>til</strong>svarende m˚ate, for et vilk˚arlig dyp z1 (se figur 1) gjelder<br />

hvor ∆z = H − z1.<br />

p| b+h<br />

H = −gρ0(b + h − H) (17)<br />

ps − p(H) = −gρ0η (18)<br />

p(H) = ps + gρ0η(x, y, t) (19)<br />

p(z1) = gρ0∆z + ps + gρ0η(x, y, t) (20)<br />

Over tid er det generelt (meget) sm˚a romlige variasjoner i overflatetrykket ps. Vi kan derfor se<br />

bort fra bidrag fra ∇ps. Fra uttrykkene (19) og (20) følger det da at<br />

∂p(H)<br />

∂x<br />

= ∂p(z1)<br />

∂x<br />

∂η<br />

= gρ0<br />

∂x<br />

Tilsvarende sammenheng gjelder for ∂/∂y. Det er derfor kun overflatehevningen η som gir opp<strong>hav</strong><br />

<strong>til</strong> trykk-kraft i en homogen væske. Av denne grunn kalles sammenhengen p = gρ0η for dynamisk<br />

trykk.<br />

2.4 Grunntvannsligningene<br />

Overst˚aende gir grunntvannsligningene uttrykt ved u, v, h og η<br />

For flat bunn har vi at<br />

3 Gravitasjons<strong>bølger</strong><br />

(21)<br />

∂u<br />

∂u<br />

+ u∂u + v − fv<br />

∂t ∂x ∂y<br />

=<br />

∂η<br />

−g<br />

∂x<br />

(22)<br />

∂v ∂v ∂v<br />

+ u + v + fu<br />

∂t ∂x ∂y<br />

=<br />

∂η<br />

−g<br />

∂y<br />

(23)<br />

∂η ∂ ∂<br />

+ (hu) + (hv)<br />

∂t ∂x ∂y<br />

= 0 (24)<br />

h(x, y, t) = H + η(x, y, t) (25)<br />

Gravitasjons<strong>bølger</strong> forekommer p˚a grenselaget mellom <strong>hav</strong> og atmosære, eller mer generelt mellom<br />

to væsker med ulik tetthet. I det følgende betrakter vi overflate<strong>bølger</strong>. Kraften som virker p˚a<br />

disse bølgene er tyngdekraften, følgelig kalles bølgene for (overflate) gravitasjons<strong>bølger</strong>.<br />

Dersom vi antar at bølgenes bølgelengde er stor i forhold <strong>til</strong> vanndypet, f˚ar vi hva som kalles<br />

grunntvanns<strong>bølger</strong>. Dersom bølgene beskriver sm˚a utslag, kan bølgeligningene (22)–(24) lineariseres.<br />

Dette betyr at ethvert produkt av bølgevariable, som for eksempel adveksjonsleddene<br />

i (22) og (23), kan neglisjeres. Videre neglisjeres effekten av jordens rotasjon.<br />

4


Med f = 0 og ved ˚a betrakte en en-dimensjonal bølgebevegelse (v = 0 og ∂/∂y = 0), kan (22)<br />

og (24) skrives som<br />

∂u<br />

∂t<br />

∂η<br />

∂t<br />

∂η<br />

= −g<br />

∂x<br />

∂u<br />

= −H<br />

∂x<br />

u kan elimineres fra uttrykkene over ved ˚a betrakte ∂(26)/∂x og ∂(27)∂t, som gir den klassiske<br />

bølgeligningen<br />

∂2η ∂t2 − gH ∂2η = 0 (28)<br />

∂x2 Bølgeligningen kan løses ved ˚a søke løsning p˚a formen (se avsnitt A.4)<br />

(26)<br />

(27)<br />

η = Re {η0 exp[i(kx − ωt)]} (29)<br />

hvor Re betegner reell del, η0 er amplitude, k er bølgetall i x-retningen og ω er vinkelfrekvens.<br />

Innestting av (29) i (28) gir dispersjonsrelasjonen<br />

−ω 2 + gH k 2 = 0 (30)<br />

Siden en bølges fasehastighet c er gitt ved c = ω/k (se avsnitt A.4), er gravitasjonsbølgenes<br />

fasehastighet c = c0 gitt ved<br />

c0 = ± gH (31)<br />

Siden fasehastighetene c0 er uavhengig av bølgetallet k er gravitasjonsbølgene ikke-dispersive<br />

(merk at dette gjelder for <strong>bølger</strong> med stor bølgelengde i forhold <strong>til</strong> vanndypet). Dette betyr at<br />

alle gravitasjons<strong>bølger</strong> forplanter seg med en og samme fasehastighet, uavhengig av bølgenes<br />

bølgetall eller -lengde.<br />

4 Sverdrup <strong>bølger</strong><br />

Ogs˚a her betrakter vi lange overflate<strong>bølger</strong>. Men i <strong>til</strong>legg <strong>til</strong> gjennomgangen over, inkluderer<br />

vi rotasjon. Videre antar vi at væsken har uendelig, horisontal utbredelse i to dimensjoner.<br />

De resulterende bølgene er ofte kalt Sverdrup <strong>bølger</strong>. Variasjoner i jordens rotasjon neglisjeres<br />

neglisjeres, s˚a f = konst. De grunnleggende likningene følger da fra (22)–(24)<br />

∂η<br />

+ H<br />

∂t<br />

Dersom vi søker en løsning p˚a formen<br />

∂u<br />

− fv<br />

∂t<br />

=<br />

∂η<br />

−g<br />

∂x<br />

(32)<br />

∂v<br />

+ fu<br />

∂t<br />

<br />

∂u ∂v<br />

+<br />

∂x ∂y<br />

=<br />

=<br />

∂η<br />

−g<br />

∂y<br />

0<br />

(33)<br />

(34)<br />

η ∝ exp[i (kx + ly − ωt)] (35)<br />

5


hvor l = 2π/λy er bølgetallet i y-retningen (λy ier <strong>til</strong>hørende bølgelengde, se avsnitt A.4), gir<br />

dette følgende algebraiske sammenhenger<br />

−iωu − fv = −igkη0<br />

−iωv + fu = −iglη0<br />

(36)<br />

(37)<br />

−iωη0 + iH(ku + lv) = 0 (38)<br />

Ligningene (36)-(38) er tre ligninger med tre ukjente, og kan uttrykkes p˚a vektorform<br />

⎛<br />

−iω −f igk<br />

⎞ ⎛ ⎞<br />

u<br />

⎝ f −iω igl ⎠ ⎝ v ⎠ = 0 (39)<br />

iHk iHl −iω η<br />

Ikke-triviell løsning finnes n˚ar ligningsystemets determinant er lik null. Dette gir følgende sammenheng<br />

ω[ω 2 − f 2 − gH(k 2 + l 2 )] = 0 (40)<br />

eller, for det horisontale bølgetallet k 2 h = k2 + l 2 ,<br />

eller uttrykt med gravitasjonsbølgenes fasehastighet c 2 0 = gH<br />

ω[ω 2 − f 2 − gHk 2 h] = 0 (41)<br />

ω[ω 2 − f 2 − c 2 0k 2 h] = 0 (42)<br />

Uttrykkene (41) og (42) er Sverdrupblgenes dispersjonsrelasjon. Det er generelt ulik fysisk mekanisme<br />

for de ulike løsningene gitt ved en dispersjonsrelasjon. Derfor tolkes gjerne de ulike<br />

løsningene hver for seg.<br />

Siden uttrykket i klammeparantes i (41) og (42) gir at fasefarten c = ω/kh avhenger av bølgetallet<br />

kh, vil resulterende <strong>bølger</strong> med ulik bølgelengde bre seg med ulik fart. Sverdrup<strong>bølger</strong> er derfor<br />

dispersive <strong>bølger</strong> (se avsnitt A.5).<br />

4.1 Løsning ω = 0<br />

Løsningen ω = 0 av (41) er konsistent med ∂/∂t = 0. Fra de grunnleggende ligningene (32) og (33)<br />

ser vi at denne løsningen gir geostrofisk kraftbalanse. Den stasjonære (det vil si tidsinvariante)<br />

løsningen av grunntvannsligningene er alts˚a geostrofisk balanse.<br />

4.2 Løsning ω 2 = f 2 + c 2 0 k 2 h<br />

Denne løsningen viser at |ω| ≥ |f| for alle bølgeløsninger.<br />

Videre er dispersjonsrelasjonen symmetrisk med hensyn p˚a x- og y-retningene. Dette betyr at<br />

verken x- eller y-retningen har en spesiell betydning for bølgefeltet. Vi orienterer derfor koordinatsystemet<br />

slik at x-aksen er rettet i retning av bølgens forplantning. Da vil kh = k og l = 0.<br />

Den resulterende dispspersjonsrelasjonen er da<br />

ω 2 = f 2 + c 2 0 k 2<br />

6<br />

(43)


Bølgens fasehastighet c bestemmes fra sammenhengen c = ω/k. Dette gir<br />

c = ω<br />

k<br />

1 2 2<br />

= f + c0 k<br />

k<br />

21/2 1/2<br />

2 f<br />

= c0 + 1<br />

= c0<br />

k 2 c 2 0<br />

<br />

1 + 1<br />

k 2 L 2 ρ<br />

hvor Lρ = c0/f er Rossby deformasjonsradius. Fra uttrykket over framkommer det at c øker med<br />

avtagende bølgetall eller økende bølgelengde.<br />

Bølgens gruppehastighet er gitt ved cg = dω/dk. Differensieres dispersjonsrelasjonen (43), f˚ar<br />

vi<br />

2ω dω = 2kc 2 0 dk (47)<br />

eller<br />

cg = k<br />

ω c20 = c20 c<br />

1/2<br />

Følgelig gjelder cgc = c 2 0. Gruppehastigheten avtar derfor med økende bølgelengde.<br />

Løsningen av (43) er illustrert i figur 4. For sm˚a bølgetall k (lange <strong>bølger</strong>), vil bølgeegenskapene<br />

Figur 2: Grafisk frams<strong>til</strong>ling av dispersjonsrelasjonen gitt ved (43). I <strong>til</strong>legg er dispersjonsrelasjonen<br />

<strong>til</strong> gravitasjons<strong>bølger</strong> (avsnitt 3) og treghets<strong>bølger</strong> (avsnitt 5) vist. Lρ er Rossby deformasjonsradius,<br />

TBR betegner omr˚adet hvor jordrotasjonen er viktig (treghets<strong>bølger</strong>egime) og GBR<br />

betegner omr˚adet hvor jordrotasjonen spiller liten rolle (gravitasjons<strong>bølger</strong>egime).<br />

nærme seg treghets<strong>bølger</strong> (se avsnitt 5) og virkningen av jordens rotasjon f er viktig. For store<br />

bølgetall (korte <strong>bølger</strong>) er jordroatasjones virkning liten. Vi er da i regimet <strong>til</strong>hørende gravitasjons<strong>bølger</strong><br />

(avsnitt 3).<br />

4.2.1 Resulterende bølgebevegelse<br />

For l = 0 har vi fra (36) og (37) at<br />

−iωu − fv = −igkη0<br />

7<br />

(44)<br />

(45)<br />

(46)<br />

(48)<br />

(49)


(50) gir at<br />

som innsatt i (49) gir<br />

−iωv + fu = 0 (50)<br />

u = i ω<br />

v (51)<br />

f<br />

v = −i f<br />

η<br />

kH<br />

(52)<br />

(52) innsatt i (51) gir<br />

u = ω<br />

η<br />

kH<br />

(53)<br />

Siden |ω| ≥ |f|, følger det at |u| ≥ |v|, eller at bølgen forplanter seg i x-retningen.<br />

Bare reell del av løsningen har en fysisk mening. Dersom vi antar at amplituden η0 er reell, har<br />

vi at<br />

η = Re{η0 exp[i(kx − ωt)]} = η0 cos(kx − ωt) (54)<br />

(54) innsatt i (53) og (52) gir<br />

u = ω<br />

kH η0 cos(kx − ωt) (55)<br />

v = f<br />

kH η0 sin(kx − ωt) (56)<br />

Det følger at det er ingen effekt av jordrotasjonen i x-retningen, mens y-retningen p˚avirkes av<br />

jordrotasjonen.<br />

Bølgens bevegelse i et fast punkt i rommet, for eksempel for x = 0, og p˚a den nordlige halvkule<br />

(f > 0) gir følgende sammenhenger (for enkelhets skyld bruker vi A = η0/kH i det<br />

følgende):<br />

ωt = 0 : u = ωA > 0, v = 0 (57)<br />

ωt = π/2 : u = 0, v = −fA < 0 (58)<br />

ωt = π : u = −ωA < 0, v = 0 (59)<br />

Hastighetskomponentene spenner ut en ellipse som vist i figur 3 med en resulterende bevegelse<br />

rettet med klokken p˚a den nordlige halvkule. Videre er ellipsens hovedakse rettet i x-retningen,<br />

det vil si i retning av bølgens forplantningsretning. Forholdet mellom ellipsens amplituder er gitt<br />

ved forholdet |ω|/|f|.<br />

Bølgens karakteristikk for et fast tidspunkt, for eksempel for t = 0, gir følgende sammenhenger<br />

(A = η0/kH som over, og vi antar vi er p˚a nordlige halvkule, f > 0):<br />

kx = 0 : u = ωA > 0, v = 0 (60)<br />

kx = π/2 : u = 0, v = fA > 0 (61)<br />

kx = π : u = −ωA < 0, v = 0 (62)<br />

kx = 3π/2 : u = 0, v = −Af < 0 (63)<br />

Den <strong>til</strong>hørende bevegelsen er vist med ellipsene nederst i figur (4). Hastighetskomponenten i<br />

bølgens forplantningsretning, u, er vist med røde horisontale piler i samme figur. Det følger at<br />

det er konvergens i venstre halvdel og divergens i høre del av figuren. Dette medfører at overflaten<br />

heves hvor det er konvergens og senkes hvor det er divergens (røde, vertikale piler i figuren). Siden<br />

bølgens form er bevart med bevegelsen, betyr dette at bølgen forplanter seg i positiv x-retning.<br />

Det er alts˚a væskepakkenes oscillasjon i x-retningen som driver bølgen framover.<br />

8<br />

(64)


Figur 3: Illustrasjon av u- og v-komponentenes endring med tiden t p˚a den nordlige halvkule.<br />

Den resulterende bevegelsen er rettet med klokken p˚a den nordlige halvkule.<br />

Figur 4: Illustrasjon av bølge- og partikkelbevegelse <strong>til</strong> en Sverdrupbølge p˚aden nordlige halvkule.<br />

Fra (91) følger det at v og η har motsatt fortegn. Se tekst for beskrivelse av figuren.<br />

9


5 Treghetssvingninger<br />

N˚ar ω → f vil partikkelbanen som beskrevet i avsnittet g˚a over <strong>til</strong> ˚a være sirkulære. Dispersjonsrelasjonen<br />

(41) og (42) gir da at kh → 0, s˚a de resulterende bølgene har ingen romlige variasjoner,<br />

inkludert ingen trykkgradient.<br />

De resulterende ligningene, fra (22) og (23), blir da<br />

For u(t = 0) = u0, v(t = 0) = 0, gir dette<br />

∂u<br />

− fv<br />

∂t<br />

= 0 (65)<br />

∂v<br />

+ fu<br />

∂t<br />

= 0 (66)<br />

u = u0 cos ft (67)<br />

v = −u0 sin ft (68)<br />

hvor u 2 0 = u 2 + v 2 . Uttrykkene over beskriver en sirkelbevegelse i retning med klokken p˚a den<br />

nordlige halvkule. Merk at denne bevegelsen vil vare <strong>til</strong> evig tid n˚ar en ser bort fra friksjon, dette<br />

p˚a tross av at trykkgradienten er null.<br />

Banen <strong>til</strong> en partikkel gitt med posisjonen (xp, yp) følger fra uttrykkene<br />

Dette gir<br />

eller<br />

u = dxp<br />

dt<br />

xp − x0 = u0<br />

f<br />

yp − y0 = u0<br />

f<br />

og v = dyp<br />

dt<br />

(xp − x0) 2 + (yp − y0) 2 = u2 0<br />

f 2<br />

(69)<br />

sin ft (70)<br />

cos ft (71)<br />

Det følger da at partikkelen beveger seg i en lukket sikrel rundt punktet (x0, y0). Sirkelens radius<br />

er gitt ved r = u0/f. For u0 = 10 cm s −1 og f = 10 −4 s −1 , som er en representativ verdi for<br />

Coriolisparameteren p˚a midlere breddegrader, er r = 1 km og T = 2π/f ∼ 17 timer.<br />

N˚ar <strong>tidevann</strong>sfrekvensen er <strong>til</strong>nærmet lik frekvensen <strong>til</strong> treghetssvingningene vil det generelt<br />

oppst˚a sterk vekselvirkning (resonnans) mellom de to bølgebevegelsene. De breddegradene dette<br />

skjer for – rundt 30 og 75 grader sør og nord – blir gjerne kalt kritiske breddegrader. For disse<br />

breddegradene kan det forventes spesielt sterk blanding grunnet overgang fra <strong>tidevann</strong>senergi <strong>til</strong><br />

sm˚askala, vertikal blanding.<br />

6 Kelvinbølge mot kyst<br />

6.1 Utgangspunkt<br />

Vi betrakter en overflatebølge som brer seg langs en kyst, for eksempel langs y-aksen som vist i<br />

figur 5. Det kan ikke være hastighet p˚a tvers av kysten, s˚a<br />

10<br />

(72)


Figur 5: Illustrasjon p˚a en overflatebølge som brer seg langs en kyst, i dette <strong>til</strong>fellet langs y-aksen<br />

(x = 0).<br />

u = 0 ved x = 0 (73)<br />

Vi antar videre at u = 0 over alt, slik at bølgen brer seg kun i y-retningen. De lineariserte<br />

grunntvannsligningene (32)–(34) blir da<br />

6.2 Løsningsmetode<br />

−fv = −g ∂η<br />

∂x<br />

∂v<br />

∂t<br />

= −g ∂η<br />

∂y<br />

(74)<br />

(75)<br />

∂η ∂v<br />

+ H = 0 (76)<br />

∂t ∂y<br />

Kelvinproblemet kan løses ved ˚a søke løsning p˚a formen (se avsnitt A.4)<br />

(v, η) = Re {(v0, η0) exp[i(kyy − ωt)]} (77)<br />

Merk at v0, η0 kan ha en x-avhengighet, slik at v0 = v0(x) og η0 = η0(x). Dette betyr at<br />

v = v(x, y, t) og η = η(x, y, t) i (77), der x-avhengigheten kommer fra amplituden v0, η0, mens yog<br />

t-avhengigheten kommer fra bølgeformen exp[i(kyy − ωt)].<br />

Uttrykk (77) innsatt i (74)-(76) gir<br />

−fvo = −g ∂η0<br />

∂x<br />

−iωv0 = −igkyη0<br />

11<br />

(78)<br />

(79)


I der første uttrykket er det brukt at η0 = η0(x).<br />

−iωη0 + iHkyv0 = 0 (80)<br />

De to siste ligningene kan uttrykkes p˚a vektorform<br />

<br />

−iω<br />

iHky<br />

igky<br />

−iω<br />

<br />

v0<br />

<br />

= 0 (81)<br />

Ikke-triviell løsning finnes n˚ar ligningsystemets determinant er lik null. Dette gir<br />

eller<br />

Bølgens fasefart er derfor<br />

c = ω<br />

η0<br />

−ω 2 + gHk 2 y = 0 (82)<br />

<br />

ω = ±ky gH (83)<br />

ky<br />

= ± gH = ±c0<br />

Siden c ikke avhenger av ky er dette en ikke-dispersiv bølge, det vil si at enhver bølgekomponent<br />

brer seg med samme hastighet (se avsnitt A.5). Videre er bølgens gruppehastighet (avsnitt<br />

A.5)<br />

cg = ∂ω<br />

= ±<br />

∂ky<br />

gH = ±c0<br />

(85)<br />

Alts˚a er bølgens fasehastighet og gruppehastighet lik.<br />

Bølgens x-avhengighet følger ved ˚a eliminere v0 fra (78) og (80)<br />

Siden ω/ky = c og gH = c 2 0, gir dette<br />

f ω<br />

η0 = g<br />

Hky<br />

∂η0<br />

∂x<br />

∂η0<br />

∂x<br />

= ± f<br />

|c0| η0<br />

Uttrykket over sier at løsningen kan ha begge fortegn. For positivt fortegn er fasehastigheten<br />

c positiv og bølgen brer seg i positiv y-retning (se 84). For negativt fortegn er fasehastigheten<br />

negativ og bølgen brer seg i negativ y-retning.<br />

Uttrykket over har løsning<br />

η0 = η ∗ 0e ±xf/|c0|<br />

hvor η ∗ 0 = η(x = 0) er bølgens høyde ved kysten. P˚a den nordlige halvkule (f > 0) krever fysisk<br />

løsning at minustegnet velges, ellers ville bølgens amplitude g˚a mot uendelig for økende x (det<br />

vil si n˚ar vi beveger oss bort fra kysten). Følgelig forplanter bølgen seg i negativ y-retning (det<br />

vil si at c = cg < 0), eller med kysten <strong>til</strong> høyre for bevegelsen.<br />

Løsningen kan n˚a skrives p˚a formen<br />

(84)<br />

(86)<br />

(87)<br />

(88)<br />

u = 0 (89)<br />

η = η ∗ 0 cos(kyy − ωt)e −x/Lρ (90)<br />

<br />

g<br />

v = − η (91)<br />

H<br />

12


√<br />

I (90) er ω = −ky gH og Lρ = |c|/f.<br />

Siden ω = −|c0|ky fra (84), kan (90) ogs˚a uttrykkes som<br />

η = η ∗ 0 cos[ky(y + |c0|t)]e −x/Lρ (92)<br />

(90) og (92) sier at bølgens amplitude avtar n˚ar vi fjerner oss fra kysten og at Lρ er lengdeskalaen<br />

for bølgedempingen. Videre sier (92) at bølgen brer seg i negativ y-retning. Dette følger siden<br />

bølgens fase, ky(y + |c0|t), er bevart med bevegelsen. S˚a for økende tid t m˚a y avta for at<br />

y+|c0|t = konst. At bølgen brer seg i negativ y-retning er konsistent med at bølgens fasehastighet<br />

c < 0 (fra 84).<br />

6.3 Egenskaper<br />

• Kelvin<strong>bølger</strong> brer seg langs en kyst med kysten <strong>til</strong> høyre p˚a nordlige halvkule og <strong>til</strong> venstre<br />

p˚a sørlige halvkule.<br />

• Bølgefarten (b˚ade fasehastighet og gruppehastighet) er konstant og absuluttverdien er gitt<br />

ved √ gH. Kelvinbølge langs en kyst er derfor ikke-dispersiv.<br />

• Bølgefart for <strong>hav</strong>dyp H = 100 og 1000 m er p˚a henholdsvis 32 og 100 m s −1 .<br />

• Bølgetopp og bølgebunn er rettet normalt p˚a kysten.<br />

• Bølgetopp og bølgebunn avtar eksponensielt fra kysten med lengdeskala gitt ved Rossby<br />

deformasjonsradius Lρ = √ gH/f.<br />

• P˚a 40 breddegrader og for et <strong>hav</strong>dyp p˚a 100 m, er Lρ ≈ 340 km.<br />

• Kraftbalansen normalt p˚a kysten er geostrofisk (fra 74), som betyr at det er ingen bølgeforplantning<br />

i x-retningen.<br />

• Kraftbalansen langs kysten,<br />

∂v<br />

∂t<br />

= −g ∂η<br />

∂y<br />

(se 75), er drevet av trykkforskjellen generert av overflatehevningen som forklart i figur 6.<br />

• Mens Kelvinbølgen brer seg med konstant fart √ gH i negativ y-retning, beskriver væsken<br />

en sirkulær bevegelse i yz-planet (figur 6).<br />

• Kelvin<strong>bølger</strong> er generert av <strong>tidevann</strong> og av vind nær kyst.<br />

• Nordg˚aende Kelvin<strong>bølger</strong> forklarer hvorfor <strong>tidevann</strong>sutslagene er mye større p˚a fransk relativt<br />

<strong>til</strong> engelsk side av Den engelske kanal.<br />

• Kelvin<strong>bølger</strong> brer seg langs ekvator fra vest mot øst (for begge halvkuler); en kan da<br />

betrakte ekvator som en vegg slik at de ekvatorielle Kelvinbølgene forplanter seg med<br />

ekvator (veggen) p˚a høyre side p˚a den nordlige halvkule og med ekvator <strong>til</strong> venstre p˚a den<br />

sørlige halvkule. Se oppgave 3 fra eksamen i GEOF110, 14. juni 2010.<br />

13


Figur 6: Mekanisme for en Kelvinbølge som brer seg langs en kyst som er rettet langs y-aksen<br />

(p˚a den nordlige halvkule, se figur 5). Fra (91) følger det at v og η har motsatt fortegn. Positiv<br />

bølgeamplitude svarer da <strong>til</strong> negativ v-komponent, og vice versa, vist med svarte piler. Det er<br />

følgelig vekselvis konvergens og divergens mellom bølgetopp og bølgebunn (svarte, vertikalstiplede<br />

linjer). Overflaten m˚a stige der det er konvergens og falle der det er divergens, vist med røde<br />

piler. Dette, sammen med at bølgens form er bevart med bevegelsen (uttrykk 90), medfører at<br />

endring i overflateniv˚aet grunnet konvergens og divergens (røde piler) kan kun skje ved at bølgen<br />

brer seg i negativ y-retning, vist med lang bl˚a pil.<br />

14


A Generelt om <strong>bølger</strong><br />

A.1 Hva er en bølge<br />

En bølge i atmosfæren og <strong>hav</strong>et kan forklares som<br />

fysiske prosesser som transporterer informasjon i tid og rom, som for eksempel energi,<br />

uten eller med liten adveksjon av masse knyttet <strong>til</strong> transporten, og som brer seg med<br />

fart og retning som generelt er ulik den generelle atmosfære- eller <strong>hav</strong>sirkulasjonen<br />

A.2 Størrelser og egenskaper i én romlig dimensjon<br />

Enhver perturbasjon (liten endring) kan uttrykkes som summen av trignometriske (sinus eller<br />

cosinus) <strong>bølger</strong>, hvor hver bølge generelt har ulik amplitude, bølgelengde, periode og fase.<br />

A.2.1 Stasjonær bølgeform<br />

En stasjonær bølgeform i x-retningen, i <strong>til</strong>fellet under uttrykt med <strong>hav</strong>niv˚a η (enhet m), kan<br />

skrives p˚a formen<br />

<br />

η(x) = a cos 2π x<br />

<br />

λx<br />

(93)<br />

Bølgeformen η er karakterisert ved<br />

Amplitude a, slik at η varierer mellom ±a. Enhet som for den avhengige variabel,<br />

i dette <strong>til</strong>fellet m.<br />

Bølgelengde λx. Bølgens form repeteres n˚a x = ±nλx, der n er et heltall. Enhet er<br />

m.<br />

A.2.2 Bølgeform som brer seg i tid<br />

En bølgeform vil generet forplante seg i tid. Bølgens<br />

Fart, alternativt fasefart benevnes cx, se ogs˚a fasefart under. Bølgen kan bre seg<br />

i positiv og negativ x-retning. Dette kan betegnes med fasefart ±cx, hvor cx > 0.<br />

Alternativt kan cx ta b˚ade positive og negative verdier. Enhet er m s −1 .<br />

Siden fart = avstand/tid, vil bølgen bre seg en avstand x ′ = ±cxt i løpet av tiden t. (93) kan da<br />

skrives p˚a formen<br />

<br />

2π<br />

η(x, t) = a cos (x ± cxt)<br />

λx<br />

(94)<br />

Bølgen over er karakterisert ved<br />

Fase gitt ved argumentet 2π(x ± cxt)/λx (enhet rad). Punkt med konstant fase<br />

er punkter hvor bølgeformen har samme verdi, for eksempel bølgetopp eller<br />

bølgebunn.<br />

15


A.2.3 Bølgetall<br />

I stedet for ˚a bruke bølgelengde λx, er det vanlig ˚a uttrykke bølgen med bølgens<br />

Bølgetall kx (enhet m−1 ). Sammenhengen mellom bølgetall og bølgelengde er gitt<br />

ved<br />

kx = 2π<br />

(95)<br />

Bølgetallet er antall bølgelengder begrenset av lengden 2π. For eksempel vil en<br />

bølge med bølgelengde 100 km ha et bølgetall 6.3×10 −5 m −1 . Bølgetallet trenger<br />

derfor ikke ˚a være et heltall.<br />

Sies det at en bølge i atmosfæren har bølgetall én, betyr dette at én bølge dekker<br />

hele jordens omkrets (for eksempel rundt jorden p˚a 60 ◦ N). For bølgetall fire vil<br />

det være fire fulle <strong>bølger</strong> rundt jorden.<br />

Ved hjelp av (95) kan (94) skrives p˚a formen<br />

A.2.4 Andre definisjoner<br />

λx<br />

η(x, t) = a cos[kx(x ± cxt)] (96)<br />

Perioden T er tiden det tar for et punkt <strong>til</strong> ˚a repetere seg selv. T m˚a følgelig være<br />

lik tiden bølgen bruker for ˚a bre seg en bølgelengde<br />

T = λx/cx<br />

Vinkelfrekvensen (ogs˚a kalt vinkelhastigheten) ω er er et m˚al p˚a rotasjonshastigheten<br />

ω = 2π<br />

(98)<br />

T<br />

Fasefarten cx er farten <strong>til</strong> bølgen, for eksempel hvor raskt en bølgetopp eller -bunn<br />

brer seg<br />

cx = λx<br />

T<br />

= ω<br />

kx<br />

hvor vi har brukt (95) og (98) i den andre overgangen.<br />

Uttrykt med ω og kx, kan (96) skrives p˚a formen<br />

(97)<br />

(99)<br />

η(x, t) = a cos(kxx ± ωt) (100)<br />

Uttrykkene (94), (96) og (100) uttrykker det samme: For fast tid t = t0 beskriver η en bølgeform<br />

i x-retningen som repeterer seg selv med bølgelengde λx, det vil si bølgens form repeteres for<br />

x = ±nλx, hvor n er et heltall. Tilsvarende, for et fast punkt x = x0 beskriver η en st˚aende bølge<br />

med periode T , det vil si en bølge som repeterer seg selv for t = ±nT .<br />

Endelig,<br />

Fasen <strong>til</strong> en bølge er gitt ved argumentet kxx ± ωt og uttrykker hvor i bølgesyklusen<br />

en befinner seg. Fasen varierer fra 0 <strong>til</strong> 2π.<br />

16


A.3 Størrelser og egenskaper i to romlige dimensjoner<br />

Overst˚aende kan utvides <strong>til</strong> flere dimensjoner. I to dimensjoner gjelder<br />

der<br />

η = a cos(kxx + kyy ± ωt) = a cos(k · x ± ωt) (101)<br />

bølgetallsvektor k = (kx, ky) = kxˆx + ky ˆy, med lengde<br />

og retning<br />

Fasehastigheten er<br />

og bølgelengden er<br />

hvor bølgetallet k er gitt ved (102).<br />

k 2 = k 2 x + k 2 y<br />

ˆk = k<br />

k<br />

c = ω<br />

k<br />

λ = 2π<br />

k<br />

A.4 Størrelser og egenskaper uttrykt med kompleks notasjon<br />

(102)<br />

(103)<br />

(104)<br />

(105)<br />

I stedet for ˚a regne med cosinus- (eller sinus-) <strong>bølger</strong> som beskrevet over, letter det analysen ˚a<br />

uttrykke en bølge p˚a kompleks form<br />

Re {a exp[i(k · x ± ωt)]} (106)<br />

Her betegner Re reell del, a er (kompleks) amplitude, k = kxˆx + ky ˆy er bølgetallvektor i x- og<br />

y-retning, x = xˆx + yˆy er posisjonsvektor og ω er vinkelfrekvens.<br />

Siden<br />

uttrykker (106) standard bølge p˚a formen<br />

exp iψ = cos ψ + i sin ψ (107)<br />

a cos(k · x ± ωt) eller a sin(k · x ± ωt) (108)<br />

avhengig av om amplituden a er reell (som i dette <strong>til</strong>fellet gir en cosinus-bølge) eller kompleks (i<br />

dette <strong>til</strong>fellet en sinus-bølge).<br />

Bølgeformen gitt ved (106) er særdeles hensiktsmessig grunnet eksponensialfunksjonens egenskap<br />

at den deriverte av funksjonen er lik funksjonen selv, korrigert med noen algebraiske koeffisienter.<br />

Dette fører <strong>til</strong> at derivasjonsoperatorene kan erstattes med algebraiske koeffisienter<br />

∂<br />

→ ikx,<br />

∂x<br />

∂<br />

∂y → iky og ∂<br />

→ ±iω (109)<br />

∂t<br />

Dette betyr at for eksempel grunntvannsligningene kan uttrykkes som et sett av algebraiske<br />

ligninger som kan løses direkte. Løsningen gir alle mulige kombinasjoner av bølgeparametre og<br />

ligningsparametre som <strong>til</strong>fredss<strong>til</strong>ler de kontinuerlige ligningene. Dette, sammen med en fysisk<br />

tolkning av bølgeløsningen, gir en fullverdig beskrivelse av bølgene.<br />

17


A.5 Dispersjonsrelasjon<br />

Det algebraiske forholdet mellom ω og k, uttrykt som ω = f(k) (det vil si at ω er en funksjon<br />

av k), kalles bølgens dispersjonsrelasjon.<br />

A.5.1 Ikke-dispersive <strong>bølger</strong><br />

Dersom ω har en lineær avhengighet <strong>til</strong> k, det vil si at ω ∝ k, kalles bølgen ikke-dispersiv. I dette<br />

<strong>til</strong>fellet forflytter <strong>bølger</strong> seg med samme fasehastighet, se (99) og (104), uavhengig av bølgens<br />

bølgelengde.<br />

A.5.2 Dispersive <strong>bølger</strong><br />

Dersom ω ikke har en lineær avhengighet <strong>til</strong> k vil enhver bølge med ulik bølgelengde forflytte seg<br />

med ulik fasehastighet. I dette <strong>til</strong>fellet er bølgen dispersiv.<br />

A.6 Gruppefart<br />

For dispersive <strong>bølger</strong> forplanter ikke energien seg med én bølge, men med “totalbidraget” fra<br />

de ulike bølgene. Det er dette “totalbidraget” vi ser som en bølge, eksempelvis som en gravitasjonsbølge<br />

i et vannbasseng, ikke de ulike bølgekomponentene. Følgelig er gruppefarten,<br />

som representerer farten <strong>til</strong> alle bølgene sett som en helhet, gjerne mer viktig enn de ulike<br />

bølgekomponentenes fasefart.<br />

Gruppefart er gitt av uttrykket<br />

cg = ∂ω<br />

(110)<br />

∂k<br />

En god illustrasjon p˚a forholdet mellom fasefart og gruppefart er gitt p˚a siden<br />

http://www.isvr.soton.ac.uk/spcg/tutorial/tutorial/Tutorial_files/Web-further-dispersive.htm.<br />

A.6.1 Utledning<br />

Dette avsnittet er <strong>til</strong> informasjon og er ikke pensum.<br />

(110) kan vises ved ˚a betrakte to <strong>bølger</strong> med nesten like bølgetall og bølgefrekvenser, for eksempel<br />

Fra identiteten<br />

følger det at<br />

Med<br />

og<br />

η = η0 cos(k1x − ω1t) + η0 cos(k2x − ω2t) (111)<br />

cos(a ± b) = cos a cos b ∓ sin a sin b (112)<br />

cos(a + b) + cos(a − b) = 2 cos a cos b (113)<br />

a = 1<br />

1<br />

(k1 + k2)x − (ω1 + ω2)t (114)<br />

2 2<br />

b = 1<br />

1<br />

(k2 − k1)x − (ω2 − ω1)t (115)<br />

2 2<br />

18


kan (111) skrives p˚a formen<br />

η = 2η0 cos<br />

1<br />

2<br />

<br />

1<br />

1<br />

(k1 + k2)x − (ω1 + ω2)t cos<br />

2 2<br />

<br />

1<br />

(k2 − k1)x − (ω2 − ω1)t<br />

2<br />

Siden bølgene antas ˚a være <strong>til</strong>nærmet like, m˚a k1 ≈ k2 og ω1 ≈ ω2, og vi kan definere<br />

k =<br />

Innsatt i (116) gir dette<br />

k1 + k2<br />

2<br />

, ω =<br />

ω1 + ω2<br />

2<br />

En mulig variant av (118) er illustrert i figur (7).<br />

(116)<br />

, ∆k = k2 − k1 , ∆ω = ω2 − ω1 (117)<br />

<br />

1 1<br />

η = 2η0 cos ∆k x − ∆ω t cos(kx − ωt) (118)<br />

2 2<br />

Figur 7: Illustrasjon av en lineær kombinasjon av to <strong>bølger</strong> med nesten lik bølgelengde og<br />

bølgefrekvens. Bølgeformen brer seg mot høyre med fasehastighet c = ω/k og bølgelengde λ = 2π/k<br />

(heltrukken linje), mens amplituden for bølgepakken er 2η0 cos(∆k x/2 − ∆ω t/2), bølgelengden<br />

λ ′ = 4π/δk, perioden 4π/δω og gruppehastigheten cg = ∂ω/∂k (stiplede linjer).<br />

(118) er en bølge som dels brer seg som en standard bølge med formen cos(kx − ωt) med fasefart<br />

c = ω/k, men hvor amplituden 2η0 er modulert med en sakte varierende bølge cos(∆k x/2 − ∆ω t/2)<br />

med (lang) bølgelengde 4π/∆k og (lang) periode 4π/∆ω. Sistnevnte bølge brer seg med fart gitt ved<br />

bølgelengde/periode, det vil si<br />

∆ω<br />

(119)<br />

∆k<br />

eller, for sm˚a ∆ω og ∆k,<br />

∂ω<br />

(120)<br />

∂k<br />

Det er sistnevnte størrelse som angir bølgens gruppefart, som kan sees p˚a som bølgens “totalbidrag”.<br />

Forplantning av energi knyttet <strong>til</strong> bølgen er derfor assosiert med gruppefarten, ikke de individuelle<br />

bølgenes eller totalbølgens fasefart.<br />

19


B Tidal analysis<br />

coming...<br />

20


Figur 8: Computed (from top) M2, M2+S2, M2+S2+K2 and M2+S2+K2+N2 tides in Bergen.<br />

21


Figur 9: Observed (black) and computed M2+S2+K2+N2 (red) tides in Bergen.<br />

22

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!