-
Notifications
You must be signed in to change notification settings - Fork 3
/
Copy pathinput-deck
executable file
·171 lines (171 loc) · 8.64 KB
/
input-deck
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
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
#######################################################################
## ##
## ElectroStatic Vlasov-Maxwell (ESVM) code ##
## ##
## INPUT DECK ##
## ##
## Code units : ##
## - elementary charge $e$ for electrical charges, ##
## - electron mass $m_e$ for masses, ##
## - initial average plasma electron density $n_0=Zn_i$ with $Z$ and ##
## $n_i$ the plasma ion atomic number and density for densities, ##
## - inverse of Langmuir plasma electron angular frequency ##
## $1/omega_p = \sqrt{me / 4 pi n_0 e^2}$ for times, ##
## - initial plasma electron thermal velocity (standard deviation) ##
## $v_{T_{e_0}}$ for velocities and ##
## - Debye screening length ##
## $\lambda_\mathrm{Debye}=v_{T_{e_0}}/omega_p$ for space ##
## ##
## Initial commit written by Michaël J TOUATI - Dec. 2015 ##
## ##
#######################################################################
#
#######################################################################
## ##
## Simulation name ##
## ##
## simu = Simulation name ##
## (less than 60 characters without spaces or punctuations ##
## except - or _) ##
## ##
#######################################################################
#
#simu Two-stream-instability
#
#######################################################################
## ##
## Number of OpenMP threads ##
## ##
## N_th = Number of threads used by OPEN_MP ##
## if 0, the maximum possible value is used ##
## ##
#######################################################################
#
#N_th 0
#
#######################################################################
## ##
## Plasma properties - immobile ions ##
## ( Do not affect simulation results! ) ##
## T = electron temperature in eV ##
## ##
## Z = atomic number of fully ionized immobile ions ##
## ##
## ni = ion density in /cm^3 ##
## ##
#######################################################################
#
#T 1000.
#
#Z 1
#
#ni 1.e19
#
#######################################################################
## ##
## Plasma electron phase-space simulation box properties ##
## ##
## x_min = lower spatial simulation box boundary ##
## ##
## x_max = upper spatial simulation box boundary ##
## ##
## d_x = simulation box spatial cell size ##
## ##
## vx_min = lower velocity space simulation box boundary ##
## ##
## vx_max = upper velocity space simulation box boundary ##
## ##
## d_vx = simulation box velocity bin size ##
## ##
#######################################################################
#
#x_min -51.
#
#x_max 51.
#
#d_x 0.5
#
#vx_min -40.
#
#vx_max 40.
#
#d_vx 0.1
#
#######################################################################
## ##
## Simulation times ##
## ##
## cfl = CFL parameter (must be < 1) ##
## ##
## L_t = Simulation time duration ##
## ##
## dt_diag = Time interval between simulation results dumping ##
## ##
#######################################################################
#
#cfl 9.e-1
#
#L_t 20.
#
#dt_diag 0.25
#
#######################################################################
## ##
## Boundary conditions and numerical schemes ##
## ##
## maxwell = 1 for Maxwell-Ampere solver ##
## (Poisson is computed at the 1st iteration only) ##
## 2 for Maxwell-Poisson solver ##
## ##
## b_cond = 1 for absorbing or ##
## 2 for periodic ##
## ##
## scheme = Available numerical advection schemes ##
## - Linear : scheme = 1 (donor_cell), 2 (Lax-Wendroff),##
## 3 (Beam_Warming) or 4 (Fromm) ##
## - Non linear : scheme = 5 (minmod), 6 (superbee), ##
## 7 (Van_Leer), 8 (MUSCL1) or 9 (MUSCL2) ##
## ##
## b_VL = free parameter for Van Leer schemes only (1<=b_VL<=2) ##
## ##
#######################################################################
#
#maxwell 1
#
#b_cond 2
#
#scheme 6
#
#b_VL 2.
#
#######################################################################
## ##
## Plasma Physics academic cases ##
## see esvm.pdf file for the meaning of perturbation parameters ##
## ##
## perturb = 0 for no perturbation ##
## 1 for the electrostatic wakefield test case ##
## 2 for the (non)linear Landau damping test case ##
## 3 for the two stream instability test case ##
## ##
## A_pert = perturbation amplitude ##
## ##
## k_pert = perturbation wave vector ##
## ##
## omega_pert = perturbation angular frequency ##
## ##
## v_d = perturbation/Maxwell-Boltzmann equilibrium ##
## drift velocity(ies) ##
## ##
#######################################################################
#
#perturb 3
#
#A_pert 1.e-1
#
#k_pert 0.06159985595
#
#omega_pert 0.06159985595
#
#v_d 10.
#