-
Notifications
You must be signed in to change notification settings - Fork 200
/
Copy pathinputs_test_rz_psatd_LL2
119 lines (105 loc) · 2.89 KB
/
inputs_test_rz_psatd_LL2
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
# Iterations
max_step = 25
# Domain
amr.n_cell = 32 64
amr.max_level = 0
warpx.numprocs = 1 2
# Geometry
geometry.dims = RZ
geometry.prob_lo = 0 -220e-6
geometry.prob_hi = 200e-6 10e-6
boundary.field_lo = none damped
boundary.field_hi = none damped
# Algorithms
algo.current_deposition = direct
algo.charge_deposition = standard
algo.field_gathering = energy-conserving
algo.maxwell_solver = psatd
algo.particle_pusher = vay
algo.particle_shape = 3
# Numerics
warpx.do_moving_window = 1
warpx.moving_window_dir = z
warpx.moving_window_v = 1.
warpx.n_rz_azimuthal_modes = 1
warpx.use_filter = 1
warpx.verbose = 1
warpx.abort_on_warning_threshold = medium
warpx.cfl = 1.
#warpx.gamma_boost = 1.
#warpx.boost_direction = z
warpx.do_dive_cleaning = 1
warpx.do_divb_cleaning = 1
psatd.do_time_averaging = 1
# PSATD
psatd.update_with_rho = 1
#psatd.v_galilean = 0. 0. -0.9373391857121336
psatd.JRhom = "LL2"
# Particles
particles.species_names = driver driver_back plasma_e plasma_p
particles.use_fdtd_nci_corr = 0
particles.rigid_injected_species = driver
driver.species_type = electron
driver.injection_style = "gaussian_beam"
driver.x_rms = 5e-6
driver.y_rms = 5e-6
driver.z_rms = 20.1e-6
driver.x_m = 0.
driver.y_m = 0.
driver.z_m = -80e-6
driver.npart = 1000000
driver.q_tot = -1e-9
driver.momentum_distribution_type = "gaussian"
driver.ux_m = 0.
driver.uy_m = 0.
driver.uz_m = 2000.
driver.ux_th = 4.
driver.uy_th = 4.
driver.uz_th = 20.
driver.zinject_plane = 2e-3
driver.rigid_advance = true
driver.initialize_self_fields = 0
driver.do_symmetrize = 1
driver.symmetrization_order = 4
driver_back.species_type = positron
driver_back.injection_style = "gaussian_beam"
driver_back.x_rms = 5e-6
driver_back.y_rms = 5e-6
driver_back.z_rms = 20.1e-6
driver_back.x_m = 0.
driver_back.y_m = 0.
driver_back.z_m = -80e-6
driver_back.npart = 1000000
driver_back.q_tot = 1e-9
driver_back.momentum_distribution_type = "at_rest"
driver_back.initialize_self_fields = 0
driver_back.do_symmetrize = 1
driver_back.symmetrization_order = 4
plasma_e.species_type = electron
plasma_e.injection_style = "NUniformPerCell"
plasma_e.zmin = 0.
plasma_e.zmax = 0.5
plasma_e.xmin = -180e-6
plasma_e.xmax = 180e-6
plasma_e.profile = constant
plasma_e.density = 1e23
plasma_e.num_particles_per_cell_each_dim = 2 2 1
plasma_e.momentum_distribution_type = "at_rest"
plasma_e.do_continuous_injection = 1
plasma_p.species_type = hydrogen
plasma_p.injection_style = "NUniformPerCell"
plasma_p.zmin = 0.
plasma_p.zmax = 0.5
plasma_p.xmin = -180e-6
plasma_p.xmax = 180e-6
plasma_p.profile = constant
plasma_p.density = 1e23
plasma_p.num_particles_per_cell_each_dim = 2 2 1
plasma_p.momentum_distribution_type = "at_rest"
plasma_p.do_continuous_injection = 1
# Diagnostics
diagnostics.diags_names = diag1
diag1.intervals = 25
diag1.diag_type = Full
diag1.fields_to_plot = Er Ez Bt jr jz rho rho_driver rho_plasma_e rho_plasma_p
diag1.species = driver plasma_e plasma_p