Lärare
Föreläsare och examinator:
Nils E M Svanstedt,
telefon: 772 5346, nilss@chalmers.se
Matematiska Vetenskaper
Datorstudio handledare:
Bobo Feng bobo@chalmers.se
Övningsskrivning.
Skannade lösningar.
p1,p2,p3,p4,
p5,
p6,p7,p8,p9,
p10,p11
Kursbeskrivning
Kursen behandlar matematiska modeller i 1D och 2D
av processer där reaktion och produktion, samt
transportmekanismer som diffusion
och konvektion, ingår. Dessa modeller, som
typiskt baseras på
konservering av massa och värme samt vissa
konstitutiva samband,
utgörs av partiella differentialekvationer.
Huvudmomentet i kursen
består i att lära dig beräkna lösningar till dessa
ekvationer
med Finita Elementmetoden (FEM), vilken bygger på
approximation
med styckvisa polynom. Vi inleder kursen med att
studera styckvisa polynom och hur de kan användas för
att approximera givna funktioner. Vi går sedan
vidare och beskriver hur en approximation av den
okända lösningen till en differentialekvation kan
beräknas. Under kursens gång kommer du att få arbeta
med både teori, implementation av metoden, samt
modellering. Speciellt utvecklar du din egen
FEM-lösare i Matlab. Du kommer aven att lara Dig att
anvanda programvaran Puffin
från Body & Soul.
Kursen omfattar även grunderna om funktionsserier.
Kursform
Föreläsningar (4 tim/v) plus övningar i datorstudio (4 tim/v).
Arbetet med inlämningsuppgifterna spelar en viktig roll under kursen
och ger genomgång av hela kursinnehållet från teori
till praktik.
Schema
Föreläsning:
Måndag 13-15, i KS11. Torsdag 13-15, i KS101.
Datorstudio:
Tisdag 13-15, i KD1. Torsdag 8-10, i KD1.
Kurslitteratur
CDE Computational Differential Equations, Eriksson, Estep, Hansbo and Johnson, Studentlitteratur Publ.
Calculus, Adams.
De avsnitt ur CDE som främst behandlas i kursen är
Kapitel 5.1-5.6 (delar av)
Kapitel 6.2
Kapitel 8.1
Kapitel 9.2 (se även extra material om tidsberoende problem)
Kapitel 14 (Inga bevis)
Kapitel 15-20 Mer information om Poissons ekvation (Session E1), Värmeledningsekvationen (Session E2), Konvektion-Diffusion (Session E3)
Etc.
Kapitel 1-4 ger en kort och trevlig sammanfattning av analys och linjär algebra.
Detaljerat kursprogram
Detaljerad (preliminär) plan:
LV1:
Föreläsning 1, måndag 16 mars 13-15:
Vektorrummet av linjära funktioner på ett intervall. Vektorrummet av styckvis
linjära, kontinuerliga funktioner på ett intervall. Linjär interpolation.
Styckvis linjär, kontinuerlig interpolation. Interpolationsfeluppskattningar.
Bevis av interpolationsfeluppskattningar i max normen.
Kap 5.1-5.4 i CDE valda delar
Föreläsning 2, torsdag 19 mars 13-15:
Introduktion till L2-projektion, kap 5..6 i CDE.
Definition samt härledning av linjärt ekvationssystem.
Formulering och bevis av feluppskattningar av felet i L2-projektionen.
Studioövning 1, tisdag 17 mars 13-15 och torsdag 19 mars 8-10:
Piecewise Polynomial Lab.
Om du vill kan du ladda ner
("Shift-klicka" på länkarna)
de två filer som behövs för
att köra PP Lab till din egen dator:
PP.fig
och PPmod.m.
LV2:
Föreläsning 3, måndag 23 mars 13-15:
Genomgång av kvadraturformler (kap 5.5 i CDE samt extra material från denna hemsida)
med tillämpning på numerisk beräkning av integraler samt beräkning av
ekvationssystem för L2-projektion.
Finita element metoden (FEM). Definition och härledning av linjart ekvationssystem.
Kap 6 samt kap 8 i CDE valda delar.
Studiövning 3, tisdag 24 mars 13-15:
Repetition av kvadraturformler med tillämpning på numerisk beräkning av integraler.
Dessa kommer nu omedelbart att användas i numerisk L2-projektion i dagens övning.
Skapa en mapp Studio3 i din kurskatalog och ladda ned filerna
LoadVector.m och
MassMatrix.m till
denna mapp.
Du guidas sedan genom laborationen genom att klicka på
Kvadratur. L2-projektion.
Studioövning 4,torsdag 26 mars 8-10:
Fortsatt arbete med Kvadratur. L2-projektion..
Föreläsning 4, torsdag 26 mars 13-15:
Fortsättning FEM. Existens och entydighet samt feluppskattning i energinorm.
Kap 6 samt kap 8 i CDE valda delar.
LV3:
Forelasning 5, måndag 30 mars 13-15 Robin randvillkor (extra material från denna hemsida). Problemlösning FEM.
Genomgång av Neumann och Robinvilkor i FEM. Vi visar ocksa att Dirichlet villkor kan approximeras med Robinvillkor.
Sedan löser vi ett stationärt tvåpunkts randvärdesproblem med FEM.
Studiovning 5, tisdag 31 mars 13-15 Börja med att öppna, samt spara hos dig själv, funktionsfilerna
PoissonSolver1D.m och
PoissonAssembler1D.m
Du guidas nu igenom laborationen genom att klicka på
FEM i 1D
Studioovning 6, torsdag 2 april 8-10: FEM för Tidsberoende problem.
Genomgång av vad som krävs för Inlamningsuppgift 1 (dvs att ha en fungerande FEM-lösare för tidsberoende problem som
Du demonstrerar med en körning för Din lärare.)
Föreläsning 6, torsdag 2 april 13-15: Framat och bakat Euler for numerisk losning av ordinara differential ekvationer.
Finit element diskretisering av tidsberoende problem.
Påskuppehåll!
LV4:
Forelasning 7 , måndag 20 april 13-15:
Problemlösning bland Räkneuppgifter lv1-lv3.
Studiovning 7, tisdag 21 april 13-15: Arbete med Tidsberoende problem.
samt Inlämningsuppgift 1.
Studiovning 8, torsdag 23 april 8-10: Slutforande av arbete med Inlamningsuppgift 1.
Forelasning 8 , torsdag 23 april 13-15: Problemlösning. Gemensam räkning samt enskild räkning
bland Räkneuppgifter lv1-lv3.
LV5:
Föreläsning 9, måndag 27 april 13-15:
Introduktion till tvådimensionella problem, triangulering,
basfunktioner. Stationära värmeledningsekvationen i 2D.
Formulering av FEM för Poisson's ekvation. Kap 14.1-3 samt kap 15.1-2.
Ladda hem Robinrandvillkor 2D
Studioövning 9, tisdag 28 april 13-15: Introduktion till
Puffin
Samt arbete med Session E1 i datorsessionerna
i Body & Soul.
LV6:
Föreläsning 10, måndag 4 maj 13-15
Fortsättning på Finita element metoden för tvådimensionella problem. Kap 16.2, 16.4 i CDE
Tidsberoende problem i 2D. Ladda hem
Tidsberoende problem i 2D
Studioövning 10, tisdag 5 maj 13-15
Arbete med Sessionerna E1 och E2 i datorsessionerna
i Body & Soul
Studioövning 11, torsdag 7 maj 8-10:
Fortsatt arbete med Session E2 i datorsessionerna
i Body & Soul
Föreläsning 11, torsdag 7 maj 13-15:
Fortsättning på Finita element metoden för tvådimensionella problem. Kap 16.2, 16.4, 17.1-17.4 i CDE
Tidsberoende problem i 2D.
LV7:
Föreläsning 12, måndag 11 maj 13-15:
Konvektions-diffusionsproblem och reaktions-diffusionsproblem Kap 18 och 19 i CDE valda delar.
Studioövning 12, tisdag 12 maj 13-15 : Arbete med
Sessionerna E3-E5 i datorsessionerna
i Body & Soul.
Studioövning 13, torsdag 14 maj 8-10 : Arbete med
Sessionerna E3-E5 i datorsessionerna och projektet.
Föreläsning 13, tordag 14 maj 13-15:
Serier (ur Adams).
LV8:
Föreläsning 14, måndag 18 maj 13-15:
Serier
Studioövning 14, tisdag 13-15:
Räkneövning på serier
Inlämningsuppgiften
Användande av Puffin for losning av valt problem.
En skriftlig redovisning skall lamnas in. Den skall innehålla
problembeskrivning. Modell (lämplig PDE). Variations- och FEM formulering.
Diskreta ekvationssystemet. Matlab-kod. Numerisk lösning (Visualisering). Kommentar
om lösningen och modellen.
Inspiration kan Du tex få om Du studerar innehället i sessionerna E3-E5 noga.
Under resterande studiopass kommer jag att presentera ett antal
exempelproblem och visa några tidigare rapporter.
Examination
1. Att bli godkänd på två obligatoriska
uppgifter. Uppgift 1 består i att skriva en tidsberoende FEM-lösare i en rumsdimension
och demonstrera en körning på denna vid dator. Uppgift 2 beskrivs ovan. Uppgifterna är individuella
men samarbete i grupp uppmuntras.
2. En skriftlig tentamen. Tentamensproblemen kommer att bygga
till stor del på räkneuppgifterna.
Räkneuppgifter:
Räkneuppgifter vecka 1
Lösningar vecka 1
Räkneuppgifter vecka 2
Lösningar vecka 2
Räkneuppgifter vecka 3
Lösningar vecka 3
Extra materiel:
Extra materiel om Kvadratur (1D)
Extra materiel om Robin randvillkor (1D)
Extra materiel om Tidsberoende problem (1D):
Inga nya räkneuppgifter för vecka 4. Slutfor Inlamningsuppgift 1.
Räkneuppgifter vecka 5
Lösningar vecka 5
Räkneuppgifter med lösningar vecka 6
Räkneuppgifter med lösningar vecka 7
|