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 ...
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