PDE 2D

In English
Rapporten som ps-fil.
Rapporten som pdf-fil.

Poissions ekvation

Här har poissons ekvation lösts i en cirkel med radien 1.
-(uxx+uyy) = f
u = 0 på randen
f = 1 då |(x,y)-(.3,.7)|<.1, 0 annars
Det första rutnätet har ungefärlig kantlängd 0.25. Rutnätet har sedan förfinats i fyra steg med villkoret att vinklarna mellan två intilliggande trianglars normaler ej får överstiga ett föreskrivet värde. Den använda matlabkoden som zipfil (kör pdeproj).


Hängmatta

Har har poissions ekvation lösts i en rektangel.
-(uxx+uyy) = f
u = 0 på kortsidan
du/dn = 0 på långsidorna
f = 3 (istället för 1, för att få en snyggare bild) i en rektangel inuti lösningsrektangeln, 0 annars.
Med lite god vilja skulle detta kanske kunna vara en hängmatta. Den använda matlabkoden som zipfil (kör pdeproj).


Poissons ekvation i ett L

Den använda matlabkoden som zipfil (kör pdeproj).



Vågekvationen i en cirkel

Den använda matlabkoden som zipfil (kör pdeproj).
Tryck så rör den sig igen (jag glömde loopkommandot)

Tillbaka till Håkans sida