Differences between revisions 2 and 13 (spanning 11 versions)
Revision 2 as of 2010-02-01 09:47:53
Size: 274
Editor: anonymous
Comment:
Revision 13 as of 2010-02-01 10:04:26
Size: 1183
Comment:
Deletions are marked like this. Additions are marked like this.
Line 6: Line 6:
= Einführung in die Grvundlagen der Numerik und das Programmieren unter Unix/Linux = = Einführung in die Grundlagen der Numerik und das Programmieren unter Unix/Linux =

== Iterative Verfahren zur Bestimmung von Nullstellen ==

== Finite Differenzen Methode zur Lösung der Advektionsgleichung ==
=== upstream-Verfahren ===
{{{#!python

from pylab import *

L=100 #raeumlich
M=500 #zeitschritte
#c=0.1 #oder unten als sinus def
dx=1./L
dt=1./10
X=linspace(0,1,L)
u=zeros((M,L))
u[0,:]=exp(-500*(X-0.5)**2) #gaussfunktion
#u[0,45:55]=1. #rechteck

ion()
figure()
line,=plot(X,u[0,:])

for n in range(M-1):
    
    c=0.1*cos((n*dt)/10*pi)
    CFL=c*dt/dx
    for j in range(L-1):
        if CFL>0:
            u[n+1,j]=u[n,j]-CFL*(u[n,j]-u[n,j-1])
        else:
            u[n+1,j]=u[n,j]-CFL*(u[n,j+1]-u[n,j])
                                 
    line.set_ydata(u[n+1,:])
    draw()
}}}
{{attachment:upstream_start.png}}
{{attachment:upstrem.png}}

'''links:''' Startsituation
'''rechts:''' Endsituation

=== Leapfrog-Verfahren ===

Einführung in die Grundlagen der Numerik und das Programmieren unter Unix/Linux

Iterative Verfahren zur Bestimmung von Nullstellen

Finite Differenzen Methode zur Lösung der Advektionsgleichung

upstream-Verfahren

   1 from pylab import *
   2 
   3 L=100 #raeumlich
   4 M=500 #zeitschritte
   5 #c=0.1  #oder unten als sinus def
   6 dx=1./L
   7 dt=1./10
   8 X=linspace(0,1,L)
   9 u=zeros((M,L))
  10 u[0,:]=exp(-500*(X-0.5)**2) #gaussfunktion
  11 #u[0,45:55]=1. #rechteck
  12 
  13 ion()
  14 figure()
  15 line,=plot(X,u[0,:])
  16 
  17 for n in range(M-1):
  18     
  19     c=0.1*cos((n*dt)/10*pi)
  20     CFL=c*dt/dx
  21     for j in range(L-1):
  22         if CFL>0:
  23             u[n+1,j]=u[n,j]-CFL*(u[n,j]-u[n,j-1])
  24         else:
  25             u[n+1,j]=u[n,j]-CFL*(u[n,j+1]-u[n,j])
  26                                  
  27     line.set_ydata(u[n+1,:])
  28     draw()

upstream_start.png upstrem.png

links: Startsituation rechts: Endsituation

Leapfrog-Verfahren

LehreWiki: NumerikProgrammierenEinfuehrung2010 (last edited 2010-02-01 14:47:00 by AlexandraGronholz)