7.3 Luftmotstand og førsteordens separabel ODE
Drag-modeller, terminalfart, og den separable førsteordens ligningen `v(t) = v_T(1 − e^{−gt/v_T})` — med det som kobler rett videre til Python-sporet.
Kapitlet hører til sjanger G — svingninger og differensiallikninger, i den førsteordens varianten. Den numeriske halen hører til sjanger H — Python-kode og numerisk integrasjon, som er med i 7 av 7 sett (100 %).
Slik ligger poengene, ut fra sensorveiledningene i H2023, V2024 og V2025:
| Steg | Uttelling |
|---|---|
| Frilegemediagram og oppsett av likningen | bærende |
| Klassifisering: første orden, separabel | selvstendig poenggivende |
| Terminalfarten fra | egne poeng — og krever ikke at du løser likningen |
| Gjennomført separasjon og integrasjon med initialbetingelse | kjernen i oppgaven |
| Tolkning av grensetilfellene og | ofte en egen deloppgave |
Felle #4 i bokas feilregister bor her: å skrive kvadratisk luftmotstand som når farten kan skifte fortegn. Det koster typisk ett poeng i utledningen og ett til i koden. Vi bruker en hel løkke på hvorfor absoluttverdien må være der.
Høyeste prioritet — dette må sitte.
Øktforslag: kapitlet er på 65 minutter. Løkke 1–2 (modellene og terminalfarten) er en naturlig første økt, løkke 3–4 (den separable løsningen og grensene) den neste.
Forkunnskaper — sist du var her
Dette kapitlet bygger på kap. 7.1 og kap. 3.2. Det du trenger derfra, ferdig oppfrisket:
Bevegelseslikningen er Newtons 2. lov med akselerasjonen skrevet som derivert. Når ingen kraft avhenger av posisjonen, kan vi arbeide direkte med farten, og ordenen faller til én:
Likevektsløsningen finner du ved å sette den deriverte lik null. For et fallende legeme kalles den terminalfarten, og den er en gratis poenghøster: den krever ikke at likningen løses.
Separabel betyr at alt med den ukjente kan samles på den ene siden og alt med tiden på den andre. Alle likningene i dette kapitlet er separable.
Fra matematikken trenger du to ting: Variabelskifte (substitusjon) og Integrasjon av 1/x og eksponentialfunksjoner. Det ene integralet som gjør hele jobben, er
Fra kap. 1.2 trenger du dessuten for små — den brukes til å vise at løsningen går over i fritt fall like etter slipp.
Notasjonsavtale for kapitlet. Vi regner positiv retning nedover i fallproblemene, slik at tyngden blir positiv og farten positiv når legemet synker. Farten skrives som kursiv skalar med fortegn. Luftmotstandens koeffisienter er (lineær modell, enhet ) og (kvadratisk modell, enhet ) — de er ulike størrelser med ulike enheter og kan aldri byttes om. Terminalfarten er ; en periode ville vært , og snordrag , men ingen av dem opptrer her. Tyngden er , aldri . Vi bruker .
Løkke 1 — de to dragmodellene (~14 min)
Slipp et kaffefilter fra skulderhøyde. Det når full fart nesten med det samme og daler så jevnt ned. Slipp en stein fra samme høyde, og den akselererer hele veien.
Forskjellen er luftmotstanden, og den har to ansikter i dette faget. Hvilket du får, avhenger av hvor raskt og hvor stort legemet er — men i en eksamensoppgave står det alltid i teksten hvilken modell som gjelder. Din jobb er å kjenne dem igjen, sette dem opp med riktig fortegn, og vite hva de gjør med likningen.
En kraft fra lufta eller væsken på et legeme som beveger seg gjennom den, alltid rettet mot bevegelsen.
Kraften har ingen egen kilde i form av et fast legeme — den kommer av at mediet må skyves til side og av friksjon mot overflaten. Størrelsen vokser med farten, og det er nettopp den fartsavhengigheten som gjør bevegelseslikningen til en ekte differensiallikning.
Merk retningsregelen: «mot bevegelsen» er en betingelse som gjelder til enhver tid, også når legemet snur. Det er derfor modellen må skrives slik at kraften bytter fortegn samtidig med farten.
Koeffisienten har enhet . Minustegnet gjør at kraften peker mot bevegelsen: er positiv, er negativ, og omvendt. Modellen håndterer retningsskifte helt av seg selv, siden inngår i første potens.
Lineær motstand gjelder for små, langsomme legemer i et seigt medium: en liten kule i olje, et støvkorn i luft, et lodd i glyserol. Den er dessuten den ene av de to modellene som gir en enkel lukket løsning, og derfor den som brukes når oppgaven vil ha et analytisk svar.
Koeffisienten har enhet og avhenger av tverrsnittet, formen og luftas tetthet. Absoluttverdien er ikke pynt. Uten den ville leddet pekt samme vei uansett fartsretning, og et legeme på vei oppover ville blitt drevet videre oppover i stedet for bremset.
Kvadratisk motstand gjelder for større og raskere legemer i luft: en fallskjermhopper, en ball, en sykkel, en bil. Det er den realistiske modellen for det meste av dagliglivets bevegelser, og derfor den som dominerer i de numeriske oppgavene.
- På vei opp er (med positiv retning oppover). Motstanden skal peke nedover. Med blir den negativ ✓.
- På vei ned er . Motstanden skal nå peke oppover. Med blir leddet ✓, altså oppover.
- Med ville leddet vært negativt i begge tilfeller, og på vei ned ville lufta ha dyttet ballen fortere nedover. Modellen ville brutt sammen fysisk.
Intuisjon: gir størrelsen, og den siste -en gir retningen. Produktet har derfor riktig størrelse og riktig fortegn i ett og samme uttrykk.
Dette er felle #4 i bokas feilregister, og den følger deg videre inn i Del 8: i koden skrives leddet -D*abs(v[i])*v[i], aldri -D*v[i]**2. Skriver du kvadratet, taper du poeng både i utledningen og i koden — to steder for samme feil.
Når kan du forenkle? Bare når du kan garantere at farten ikke skifter fortegn i det tidsrommet du regner på. En fallskjermhopper som faller nedover hele veien, kan trygt beskrives med — men du må skrive ut betingelsen. Å forenkle uten å si hvorfor er å håpe at sensor antar det beste.
Forklar med egne ord, i to–tre setninger hver:
a) hvorfor luftmotstand gjør bevegelseslikningen til en differensiallikning som ikke kan løses med konstant--formlene,
b) hvorfor kvadratisk luftmotstand må skrives og ikke ,
c) hvorfor koeffisientene og ikke kan ha samme enhet.
Løkke 2 — terminalfarten (~12 min)
Det første spørsmålet i en dragoppgave er nesten alltid det samme: hvor fort går det til slutt? Og det gode ved spørsmålet er at du kan svare uten å løse differensiallikningen i det hele tatt.
Intuisjon: når farten er så stor at motstanden akkurat balanserer tyngden, er resultantkraften null. Da er akselerasjonen null, og farten endrer seg ikke mer. Dette er nøyaktig likevektsløsningen fra kap. 7.1, bare med et fysikknavn på.
Terminalfarten er en likevektsløsning av bevegelseslikningen: den konstante funksjonen oppfyller likningen for alle .
Legemet når aldri i endelig tid — det nærmer seg verdien eksponentielt. I praksis er farten innenfor av etter et par tidskonstanter, og da er forskjellen uten betydning.
Terminalfarten er en gratis poenghøster. Den krever bare at du setter akselerasjonen lik null, og den er ofte en egen deloppgave som kan besvares selv om du står fast på resten.
Kvadratisk motstand. gir
Legg merke til skaleringen — den er verdt en kommentar i en besvarelse. Med lineær motstand er terminalfarten proporsjonal med massen; med kvadratisk motstand vokser den bare som kvadratrota av massen. En fallskjermhopper som veier dobbelt så mye, faller derfor bare ganger så fort — ikke dobbelt så fort.
Benevningskontroll for den kvadratiske: , og kvadratrota gir ✓.
En tredje sammenheng er verdt å merke seg, fordi den forenkler alle regnestykkene i løkke 3: for lineær motstand er
altså kan tidskonstanten leses rett av terminalfarten.
En fallskjermhopper med utstyr har samlet masse . I fritt fall før utløsning er den kvadratiske dragkoeffisienten ; med åpen skjerm er den .
a) Finn terminalfarten i begge fasene.
b) Hvor stor er akselerasjonen i det øyeblikket skjermen åpner seg, hvis hopperen da har nådd terminalfart i første fase?
c) Kommenter tallet i b).
a) Fritt fall før utløsning:
nedover, altså .
Med åpen skjerm:
nedover, altså — omtrent farten du får ved å hoppe ned fra halvannen meters høyde, som er nettopp poenget med en fallskjerm.
b) I åpningsøyeblikket er farten fortsatt , men dragkoeffisienten er blitt . Med positiv retning nedover:
altså om lag rettet oppover (bremsing).
c) Det er ganger tyngdeakselerasjonen — helt urealistisk for et menneske, som tåler i størrelsesorden i kort tid. Konklusjonen er at modellen med en brå overgang mellom to konstante -verdier er for grov: i virkeligheten folder skjermen seg ut over et par sekunder, og vokser gradvis. En bedre modell lar variere med tiden, og den kan bare løses numerisk.
Sensorblikk: en deloppgave som ber om en kommentar til et tall, gir full uttelling for et presist resonnement uten regning. Å regne ut og la det stå ukommentert er å la poeng ligge — det er nettopp urimeligheten som er svaret.
Et kaffefilter med masse faller gjennom lufta. Ved de lave fartene det er snakk om, er motstanden godt beskrevet av den lineære modellen med .
a) Finn terminalfarten.
b) Hva er akselerasjonen i slippøyeblikket?
c) Hva er akselerasjonen når farten er ?
d) To like filtre stables oppi hverandre, slik at massen dobles mens formen og dermed er uendret. Hva blir terminalfarten da?
En hagleskål og en tennisball har samme form og samme dragkoeffisient , men ulik masse. Motstanden er kvadratisk.
a) Finn terminalfarten for en ball med masse .
b) Hvor mange ganger så tung må en ball være for å få dobbelt så stor terminalfart?
c) En medstudent sier: «Da faller en dobbelt så tung ball dobbelt så fort.» Hva er galt?
d) Ved hvilken fart er luftmotstanden akkurat halvparten av tyngden for ballen i a)?
Løkke 3 — den separable løsningen (~20 min)
Nå skal vi finne hele forløpet, ikke bare sluttfarten. Vi tar den lineære modellen, som er den ene som gir et pent svar, og vi gjør det steg for steg — for det er utregningen sensor følger.
Bevegelseslikningen for et legeme som slippes fra ro, med positiv retning nedover:
Klassifiseringen først, som alltid: første orden, lineær, inhomogen — og separabel. Det er separabiliteten som gir metoden.
og siden :
Steg 2 — separer. Alt med til venstre, alt med til høyre:
Steg 3 — integrer begge sider fra start til tidspunktet , med :
Steg 4 — løs for . Ta eksponentialfunksjonen på begge sider:
og siden , kan det også skrives
Intuisjon steg for steg: minustegnet i steg 3 kommer av at står med negativt fortegn i nevneren — det er stedet flest mister et fortegn. Og legg merke til at integrasjonsgrensene bærer initialbetingelsen: skriver du dem inn med det samme, slipper du en ubestemt konstant helt.
Kontroll ved innsetting. Deriver løsningen:
siden . Likningen er oppfylt for alle , og ✓.
Etter én tidskonstant er altså av terminalfarten nådd, etter to , etter tre og etter fem over . Tommelfingerregelen «tre tidskonstanter, så er du praktisk talt framme» er verdt å ha.
For den lineære modellen gjelder dessuten : kjenner du terminalfarten, kjenner du tidsskalaen. En tung, langsom gjenstand har lang ; et lite, lett legeme når terminalfarten nesten umiddelbart.
En liten isoporkule med masse slippes fra ro. Luftmotstanden er lineær med .
a) Finn terminalfarten og tidskonstanten.
b) Skriv .
c) Hvor lang tid tar det før kula har nådd av terminalfarten?
d) Hvor langt har kula falt etter , og hvor langt ville den falt i fritt fall?
a)
Kontroll: ✓ — de to uttrykkene må gi samme tall.
b)
c) Vi krever , altså :
altså om lag tre tidskonstanter, som tommelfingerregelen sier.
d) Strekningen er integralet av farten:
I fritt fall ville kula falt på samme tid. Luftmotstanden har altså mer enn halvert fallstrekningen på ett sekund.
Benevningskontroll: har enhet ✓, som seg hør og bør for en strekning.
Rimelighetskontroll: en isoporkule som daler ned med under etter et kvart sekund, stemmer godt med det du ser når du slipper en. En stein av samme størrelse ville hatt hundre ganger større masse og dermed hundre ganger større terminalfart — og fortsatt akselerert.
Et lodd med masse slippes fra ro ned i en beholder med olje. Motstanden er lineær med .
a) Finn terminalfarten og tidskonstanten.
b) Skriv og finn farten etter .
c) Når er farten halvparten av terminalfarten?
d) Vis at løsningen din oppfyller differensiallikningen.
En liten kule slippes fra ro i en væske. Farten måles til etter og til etter . Motstanden antas lineær, og kula slippes fra ro.
a) Sett opp de to likningene målingene gir.
b) Bestem terminalfarten og tidskonstanten.
c) Hva er kulas masse hvis ?
d) Hvorfor er det ikke nok med én måling for å bestemme både og ?
Løkke 4 — grensene, og der den lukkede løsningen ryker (~15 min)
En løsning er ikke ferdig kontrollert før du har sett på grensene. For er det to av dem, og begge er faste deloppgaver på eksamen.
Grensen . Da går , og
Løsningen nærmer seg terminalfarten, som den skal — og det bekrefter at likevektsløsningen vi fant i løkke 2, faktisk er det bevegelsen ender i.
Grensen . Her bruker vi Taylor-utviklingen fra kap. 1.2, :
Det første leddet er nøyaktig fritt fall: like etter slipp er farten null, motstanden derfor null, og legemet faller som om lufta ikke fantes. Det andre leddet er den første korreksjonen, og den er negativ — motstanden begynner å bite.
Intuisjon: modellen må gi fritt fall i starten og terminalfart til slutt. Gjør den ikke det, er noe galt. Dette er den billigste kvalitetskontrollen som finnes på et slikt uttrykk, og på eksamen er den ofte en egen deloppgave med egne poeng.
For et legeme som slippes fra ro og faller nedover hele veien, er , og likningen
er fortsatt separabel. Den kan faktisk løses i lukket form:
Du kan verifisere den ved innsetting akkurat som før. Så påstanden «kvadratisk motstand kan aldri løses analytisk» er feil, og den skal du ikke skrive.
Det som derimot er sant, er hvor smalt gyldighetsområdet er. Løsningen over holder bare så lenge
1. farten ikke skifter fortegn — så snart legemet snur, må skrives ut, og likningen skifter form ved vendepunktet;
2. bevegelsen er i én dimensjon — i to dimensjoner kobles - og -likningene gjennom farten , og systemet har ingen løsning i elementære funksjoner;
3. koeffisientene er konstante — varierer med høyden, som i tynnere luft, ryker den også.
Et skrått kast med luftmotstand bryter både punkt 1 og punkt 2, og det er nettopp det som gjør Del 8 nødvendig. Kastet er den vanligste innpakningen av den numeriske deloppgaven i eksamenssettene, og formelarket hjelper deg ikke der: den ferdigheten må sitte i fingrene.
Å bytte om og . De har ulike enheter ( mot ) og gir helt ulike terminalfarter. En benevningskontroll fanger byttet på fem sekunder.
Å miste fortegnet i separasjonen. Integralet gir , ikke . Mistet minustegn her gir en løsning som vokser uten grense i stedet for å flate ut — og en løsning som vokser mot uendelig er alltid et varsel om at et fortegn er borte.
Å glemme initialbetingelsen. Uten den står du igjen med en ubestemt konstant, og svaret er ikke ferdig. Legg initialbetingelsen inn som integrasjonsgrenser med det samme, så slipper du hele problemet.
Felle #3 — å «integrere» likningen som om var konstant. Gir null poeng; se kap. 7.1.
Å påstå at kvadratisk motstand i én dimensjon ikke kan løses analytisk. Den kan det, med en . Det som ikke går, er 2D-tilfellet og tilfeller der farten skifter fortegn — vær presis på hvilket av dem du snakker om.
Å oppgi farten uten retning. Terminalfarten er en vektorstørrelse: skriv « nedover», ikke bare tallet. Manglende retning på en vektorstørrelse koster typisk 2 poeng.
Ta utgangspunkt i for et legeme som slippes fra ro med lineær luftmotstand.
a) Vis at uttrykket gir fritt fall for små , ved å Taylor-utvikle til og med andreordens ledd.
b) Vis at når , og forklar fysisk hva det betyr.
c) For hvilken fart er akselerasjonen halvparten av ?
d) En medstudent har regnet ut . Hvordan ser du at det er feil uten å regne på nytt?
a) Sett opp bevegelseslikningen for hele bevegelsen, både opp og ned, og forklar hvorfor absoluttverdien er nødvendig.
b) Skriv likningen slik den gjelder på vei opp, og slik den gjelder på vei ned.
c) Kommer kula opp med samme fart som den kom ned igjen med, eller ikke? Begrunn fysisk uten å løse likningen.
d) Kan bevegelsen beskrives med ett enkelt lukket uttrykk for hele veien? Begrunn, og si hva neste steg ville vært.
a) To fallskjermhoppere med samme utstyr og samme kroppsstilling, men ulik masse, hopper samtidig. Hvem når bakken først, og hvorfor?
b) En medstudent skriver: «Luftmotstand gir alltid en analytisk uløselig likning, så alt må gjøres numerisk.» Rett opp utsagnet så presist du kan.
c) Hvorfor er terminalfarten en gratis poenghøster i en eksamensoppgave, og hva må du likevel huske å skrive?
Begrepsbank
Begrepsbanken er flashcard- og repetisjonsstoff — den gjentar det du nettopp har lest. Hopp trygt over ved førstegangslesing; tidsanslaget for kapitlet gjelder kjernestoffet.
og hver side integreres for seg. Metoden virker fordi og kan behandles som differensialer — det er samme regnegrep som ved substitusjon i et integral.
Alle dragmodellene i dette kapitlet er separable. Legger du integrasjonsgrensene inn med det samme (fra til på venstre side, fra til på høyre), slipper du både integrasjonskonstanten og et vanlig fortegnsrot.
Egenskapen at når : i det første øyeblikket etter slipp er farten null, motstanden derfor null, og legemet faller som i vakuum.
Grensen er en obligatorisk kontroll av ethvert dragresultat. Gir uttrykket ditt noe annet enn for små , er det galt — og kontrollen tar ti sekunder med Taylor-utviklingen .
Samme matematiske form dukker opp i radioaktiv nedbrytning, i avkjøling av en kopp kaffe og i utladning av en kondensator. Det er ikke tilfeldig: alle tre er systemer der endringsraten er proporsjonal med avstanden til likevekten.
Grovregelen er at lineær motstand gjelder for små, langsomme legemer i et seigt medium, og kvadratisk motstand for større, raskere legemer i luft eller vann.
Skillet skyldes at strømningen rundt legemet skifter karakter: ved lave farter siger mediet jevnt forbi, ved høyere farter dannes virvler bak legemet. Overgangen beskrives av et dimensjonsløst forholdstall som ikke er pensum i dette emnet.
På eksamen står det alltid i oppgaveteksten hvilken modell som gjelder. Din jobb er å bruke den riktig — ikke å avgjøre hvilken som er fysisk mest treffende.
vokser med tverrsnittet og med luftas tetthet, og synker med en strømlinjeformet form. En fallskjerm har en som er over hundre ganger større enn kroppens egen.
Enhetsforskjellen er den raskeste kontrollen på at du har brukt riktig modell: kommer du ut med en terminalfart som har feil enhet, har du blandet dem.
Funksjonen starter med stigningstall 1 i origo og flater ut mot 1, akkurat som eksponentialløsningen — kurvene ligner, men er ikke like.
Gyldigheten er smal: løsningen krever at farten ikke skifter fortegn og at bevegelsen er endimensjonal. Bryter du en av forutsetningene, må likningen løses numerisk.
Hver komponentlikning inneholder altså begge fartskomponentene, og de kan ikke løses hver for seg. Det er denne koblingen som gjør at et skrått kast med luftmotstand ikke har noen løsning i elementære funksjoner — og som gjør Del 8 til en nødvendig del av faget, ikke et vedheng.
Fordi kraften alltid peker mot bevegelsen, er skalarproduktet negativt i hvert eneste tidssteg.
Konsekvensen er at mekanisk energi ikke er bevart når luftmotstand er med, og at en kule som kastes opp alltid kommer ned med lavere fart enn den ble kastet opp med. Se kap. 5.1 for arbeidsbegrepet.
Dette er også grunnen til at luftmotstand aldri kan «hjelpe» et legeme framover: den er en ren tapskraft i mekanisk forstand.
For store nærmer dette seg : en rett linje med stigningstall , forskjøvet nedover i forhold til den rene terminalfart-linja. Forskyvningen er «etterslepet» fra oppstartsfasen, altså den strekningen legemet tapte mens det ennå ikke hadde full fart.
Den faste siste deloppgaven i en dragoppgave: «skriv et program, eller pseudokode, som løser likningen numerisk». Den er med i 7 av 7 sett (100 %).
Poengene gis for riktige initialbetingelser, et riktig akselerasjonsuttrykk og en riktig oppdateringsløkke — ikke for syntaks. Metoden emnet bygger på heter Euler–Cromer, og den utvikles fra grunnen i Del 8.
Det du kan gjøre allerede nå, er å skrive bevegelseslikningen på den formen koden trenger: akselerasjonen isolert på venstre side, motstandsleddet skrevet med absoluttverdi.
Den er uavhengig av selv og bestemmes bare av tidskonstanten. Det gjør den til et praktisk mål i laboratorieforsøk: mål tiden til halv sluttfart, så har du og dermed forholdet uten å måtte kjenne noen av dem hver for seg.
Betingelsen som må være oppfylt for at kan forenkles til : farten skifter ikke fortegn i det tidsrommet du regner på.
En fallskjermhopper som synker hele veien, oppfyller den. En ball som kastes opp og faller ned igjen, gjør det ikke — der er toppunktet et skille mellom to ulike likninger.
Forenklingen er lov, men betingelsen skal skrives ut. En besvarelse som bruker uten et ord om hvorfor, ser ut som felle #4 selv om regningen er riktig.
Legemet starter altså alltid som i fritt fall, og kurven har stigningstall i origo. Dette gjelder både for lineær og kvadratisk motstand.
Det er en av de raskeste kontrollene på et fartsuttrykk: deriver det og sett . Får du noe annet enn , er svaret galt.
De to er bare like hvis fallhøyden er stor nok til at terminalfarten rekker å bli nådd, altså minst noen tidskonstanter. Faller et lodd i ett tiendedels sekund, er sluttfarten langt under .
Å svare med på et spørsmål om farten ved bakken er en klassisk feil, og den er lett å unngå: sjekk om falltiden er stor sammenliknet med .
Den vanligste innpakningen av den numeriske deloppgaven: et legeme kastes skrått opp, og motstanden er kvadratisk.
Uten motstand er - og -bevegelsen uavhengige, og oppgaven løses med kastformlene fra kap. 2.2. Med motstand kobles de gjennom farten , banen blir usymmetrisk, og rekkevidden blir kortere enn den klassiske formelen sier.
Det finnes ingen lukket løsning, og oppgaven ender derfor alltid i en numerisk deloppgave.
Å bestemme modellens ukjente størrelser — typisk og — ved å sette måledata inn i den kjente løsningsformen.
Regelen er at du trenger like mange uavhengige målinger som du har ukjente parametre. Én måling låser ikke to parametre, uansett hvor nøyaktig den er.
Et nyttig grep er å dele to måleuttrykk på hverandre: da faller den ene parameteren bort, og du står igjen med én likning i én ukjent.
Minustegnet kommer av kjerneregelen: den indre funksjonen har derivert . Dette er stedet flest mister et fortegn i hele Del 7.
Kontrollen er enkel: en løsning som vokser uten grense i stedet for å flate ut mot , har mistet nettopp dette minustegnet.
Dette kapitlet er skrevet av Anthropics toppmodeller (Claude Opus og Claude Fable) og er foreløpig ikke manuelt gjennomgått — kvalitetskontrollen gjøres av uavhengige KI-agenter, og innmeldte feil rettes fortløpende. Funnet en feil? Meld fra, så retter vi den. Les mer om hvordan innholdet lages.
Skolesaga er en uavhengig læringsressurs og er ikke tilknyttet eller godkjent av den aktuelle utdanningsinstitusjonen. Dette er ikke offisielt studiemateriell. Les mer.