diff options
| author | Erich Eckner <git@eckner.net> | 2015-08-11 14:17:49 +0200 |
|---|---|---|
| committer | Erich Eckner <git@eckner.net> | 2015-08-11 14:17:49 +0200 |
| commit | 5eaf0e1928f4024ac5c1935dfecfe7d29f93417c (patch) | |
| tree | a4d34bd98d8e7152cc822e1b546a7415ff7e1084 /input.plap | |
| parent | bf51ca3a16d648c5da41719b06f0f3a598a37ec3 (diff) | |
| download | Plasmapropagation-5eaf0e1928f4024ac5c1935dfecfe7d29f93417c.tar.xz | |
lauffaehig, es entstehen aber noch "Ripple" im Phasenraum ...
Diffstat (limited to 'input.plap')
| -rw-r--r-- | input.plap | 32 |
1 files changed, 17 insertions, 15 deletions
@@ -7,11 +7,11 @@ allgemein runge-Kutta-3/8 # runge-Kutta-4 # euler-Vorwärts - ortsschritt 5*10^-3 * λ - zeitschritt 5*10^-3 * T - maximalimpuls 10 - impulsschritt 5 * 10^-3 - zeit 40 * T + ortsschritt 10^-2 * λ + zeitschritt 10^-2 * T + maximalimpuls 5 + impulsschritt 10^-2 + zeit 2 * T !setze $breite: (5 * λ) breite $breite mit Fortschrittsanzeige @@ -20,13 +20,13 @@ allgemeinEnde ausgaben prefix /home_raid/erich/Dokumente/Prograemmchen/Plasmapropagation/Daten/ suffix _test.dat - felder AX,AY,dAYDT,EX,BY,BZ,Rho1,Rho2,JX1 + felder AX,AY,dAYDT,EX,BY,BZ,Rho1,Rho2,VX1,PX1,PXSqr1,N1,PXRipple1 ausgabenEnde !setze $tFwhm: (2.5 * T) !setze $tMitte: (1 * T) -licht von links 2 * 2^(-2*((t-$tMitte)/$tFwhm)^2) * ω * cos(ω*t) # Zeitableitung des A-Feldes +# licht von links 2 * 2^(-2*((t-$tMitte)/$tFwhm)^2) * ω * cos(ω*t) # Zeitableitung des A-Feldes !setze $IonenMassenFaktor: (1836.15267245 + 1838.68366158) @@ -36,7 +36,8 @@ teilchen1 maximaldichte 10 minimaldichte 0 # 10^-2 maximales dLnN/dX 10 - maximales dLnN/dP 10 + maximales dLnN/dP 2 + impulsbreite 1/20 !setze $profilbreite: (4 * λ) !setze $randbreite: (0.1 * λ) verteilung stückweise @@ -58,17 +59,18 @@ teilchen2 maximaldichte nmax1 minimaldichte nmin1 maximales dLnN/dX 10 - maximales dLnN/dP 10 + maximales dLnN/dP 2 + impulsbreite 1/20 # verteilung wie teilchen1 verteilung stückweise 0 - ($breite-$profilbreite)/2 - $randbreite - sin((x - ($breite-$profilbreite)/2 - $randbreite)*π/2/$randbreite)^2 - ($breite-$profilbreite)/2 + (0.1 + $breite-$profilbreite)/2 - $randbreite + sin((x - (0.1 + $breite-$profilbreite)/2 - $randbreite)*π/2/$randbreite)^2 + (0.1 + $breite-$profilbreite)/2 1 - ($breite+$profilbreite)/2 - sin((x - ($breite+$profilbreite)/2 + $randbreite)*π/2/$randbreite)^2 - ($breite+$profilbreite)/2 + $randbreite + (0.1 + $breite+$profilbreite)/2 + sin((x - (0.1 + $breite+$profilbreite)/2 + $randbreite)*π/2/$randbreite)^2 + (0.1 + $breite+$profilbreite)/2 + $randbreite 0 stückweiseEnde teilchen2Ende |