diff options
Diffstat (limited to 'input.plap')
| -rw-r--r-- | input.plap | 33 |
1 files changed, 23 insertions, 10 deletions
@@ -7,10 +7,10 @@ allgemein runge-Kutta-3/8 # runge-Kutta-4 # euler-Vorwärts - ortsschritt 10^-3 * λ - zeitschritt 10^-3 * T - minZeitschritt 10^-10 * T - zeit 20 * T + ortsschritt 2*10^-3 * λ + zeitschritt 2*10^-3 * T + minZeitschritt 2*10^-3 * T + zeit 40 * T diffusionsterm 2 # max. p/Lp !setze $breite: (5 * λ) breite $breite @@ -26,12 +26,14 @@ ausgabenEnde !setze $tFwhm: (2.5 * T) !setze $tMitte: (1 * T) -licht von links 0.1 * ω * 2^(-2*((t-$tMitte)/$tFwhm)^2) * cos(ω*t) # Zeitableitung des A-Feldes +# licht von links 0.1 * 2^(-2*((t-$tMitte)/$tFwhm)^2) * ω * cos(ω*t) # Zeitableitung des A-Feldes + +!setze $IonenMassenFaktor: (1836.15267245 + 1838.68366158) teilchen1 spezifische Ladung -q/me maximaldichte 10 -# minimaldichte 1 + minimaldichte 0 # 10^-2 !setze $profilbreite: (4 * λ) !setze $randbreite: (0.1 * λ) verteilung stückweise @@ -48,10 +50,21 @@ teilchen1 teilchen1Ende teilchen2 - spezifische Ladung (q/me)/2 - maximaldichte nmax1*2 -# minimaldichte nmin1*2 - verteilung wie teilchen1 + spezifische Ladung q/$IonenMassenFaktor/me + maximaldichte nmax1*$IonenMassenFaktor + minimaldichte nmin1*$IonenMassenFaktor +# verteilung wie teilchen1 + verteilung stückweise + 0 + (0.1*λ + $breite-$profilbreite)/2 - $randbreite + sin((x - (0.1*λ + $breite-$profilbreite)/2 - $randbreite)*π/2/$randbreite)^2 + (0.1*λ + $breite-$profilbreite)/2 + 1 + (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 Dateiende |