AW: kleine numerische Tricks
Hi
dw/dt + A(x)*dw/dx + C = 0
dw1/dt+q/a21/2*( M1 + M2 ) +C1 = 0
dw2/dt+1/2*(-M1 + M2 ) +C2 = 0
M1=(a11-q)*a21/q*dw1/dx-(a11-q)*dw2/dx
M2=(a11+q)*a21/q*dw1/dx+(a11+q)*dw2/dx
In einer Wellengleichung waeren M1 und M2 jeweils die einlaufenden und Auslaufenden Wellenanteile. (Eine Wellengleichung hat stets zwei Loesungen)
Ueber M1 und M2 lassen sich diese Anteile vorgeben. Bei einer Wellengleichung kann man z.B. eine Ausbreitungsrichtung unterdruecken (Absorbtion am Rand).
Laesst sich das System entkoppeln, so stellt jede Gleichung einen Transportvorgang dar. Ueber eine Koordinatentransformation ueber die man sich in der Transportgeschwindigkeit mitbewegt wird d/dt=0.
Mittels der Charakteristikenmethode lassen sich somit einfache Wellengleichungen auch analytisch loesen. In krummlinigen Koordinaten laesst sich in der Regel das System nicht mehr entkoppeln. Es treten geometrische Quellterme in der Konstanten C auf. Das System bleibt darueber verkoppelt.
Gibt man im iterativen Kreis spezielle Bedingungen fuer M vor, so werden die Variablen selbst komplexwertig und erzeugen einen Kreis. Das ist weniger interessant.
Ge?ndert von richy (01.04.10 um 15:58 Uhr)
|