1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
|
# Parameterdatei für Plasmapropagation
allgemein
# runge-Kutta-14
# runge-Kutta-12
# runge-Kutta-10
runge-Kutta-3/8
# runge-Kutta-4
# euler-Vorwärts
ortsschritt 10^-2 * λ
zeitschritt 10^-2 * T
maximalimpuls 5
impulsschritt 2*10^-2
zeit 5 * T
!setze $breite: (5 * λ)
breite $breite
mit Fortschrittsanzeige
Ausgabeverzeichnis aufräumen
allgemeinEnde
ausgaben
prefix /nfs/synologynfs/erich/Plasmapropagation/Daten/
suffix _test.dat
felder AX,AY,dAYDT,EX,BZ,VX1,PX1,PXSqr1,N1,N2,PXRipple1,Rho1,Rho2
phasenraum N1,N2
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
!setze $IonenMassenFaktor: (1836.15267245 + 1838.68366158)
teilchen1
ladung -qe
masse me
maximaldichte 10
minimaldichte 0 # 10^-2
maximales dLnN/dX 10
maximales dLnN/dP 2
impulsbreite 1/20
!setze $profilbreite: (4 * λ)
!setze $randbreite: (0.1 * λ)
verteilung stückweise
0
($breite-$profilbreite)/2 - $randbreite
sin((x - ($breite-$profilbreite)/2 - $randbreite)*π/2/$randbreite)^2
($breite-$profilbreite)/2
1
($breite+$profilbreite)/2
sin((x - ($breite+$profilbreite)/2 + $randbreite)*π/2/$randbreite)^2
($breite+$profilbreite)/2 + $randbreite
0
stückweiseEnde
teilchen1Ende
teilchen2
ladung qe
masse $IonenMassenFaktor * me
maximaldichte nmax1
minimaldichte nmin1
maximales dLnN/dX 10
maximales dLnN/dP 2
impulsbreite 1/20
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
|
