|
| 1 | +################################# |
| 2 | +####### GENERAL PARAMETERS ###### |
| 3 | +################################# |
| 4 | +max_step = 500 |
| 5 | +amr.n_cell = 8 |
| 6 | +amr.max_level = 0 |
| 7 | +geometry.dims = 1 |
| 8 | + |
| 9 | +warpx.do_electrostatic = labframe |
| 10 | +algo.field_gathering = energy-conserving |
| 11 | +algo.particle_shape = 2 |
| 12 | +warpx.use_filter = 0 |
| 13 | + |
| 14 | +################################# |
| 15 | +############ CONSTANTS ############# |
| 16 | +################################# |
| 17 | +my_constants.n0 = 1.e30 # plasma density, m^-3 |
| 18 | +my_constants.Ti = 100. # ion temperature, eV |
| 19 | +my_constants.Te = 100. # electron temperature, eV |
| 20 | +my_constants.wpe = q_e*sqrt(n0/(m_e*epsilon0)) # electron plasma frequency, radians/s |
| 21 | +my_constants.de0 = clight/wpe # skin depth, m |
| 22 | +my_constants.dt = ( 0.2 )/wpe # time step size, s |
| 23 | +warpx.const_dt = dt |
| 24 | + |
| 25 | +################################# |
| 26 | +####### GENERAL PARAMETERS ###### |
| 27 | +################################# |
| 28 | +geometry.prob_lo = 0.0 |
| 29 | +geometry.prob_hi = 10.*de0 |
| 30 | +warpx.serialize_initial_conditions = 1 |
| 31 | +warpx.verbose = 1 |
| 32 | + |
| 33 | +################################# |
| 34 | +###### BOUNDARY CONDITIONS ###### |
| 35 | +################################# |
| 36 | +boundary.field_lo = periodic |
| 37 | +boundary.field_hi = periodic |
| 38 | +boundary.particle_lo = periodic |
| 39 | +boundary.particle_hi = periodic |
| 40 | + |
| 41 | +################################# |
| 42 | +############ PLASMA ############# |
| 43 | +################################# |
| 44 | +particles.species_names = electrons protons |
| 45 | + |
| 46 | +electrons.species_type = electron |
| 47 | +electrons.injection_style = "NUniformPerCell" |
| 48 | +electrons.num_particles_per_cell_each_dim = 4 |
| 49 | +electrons.profile = constant |
| 50 | +electrons.density = n0 |
| 51 | +electrons.momentum_distribution_type = gaussian |
| 52 | +electrons.ux_th = sqrt(Te*q_e/m_e)/clight |
| 53 | +electrons.uy_th = sqrt(Te*q_e/m_e)/clight |
| 54 | +electrons.uz_th = sqrt(Te*q_e/m_e)/clight |
| 55 | +electrons.xmin = 0 |
| 56 | +electrons.xmax = 10.*de0 |
| 57 | +electrons.zmin = 0 |
| 58 | +electrons.zmax =10.*de0 |
| 59 | + |
| 60 | +protons.species_type = proton |
| 61 | +protons.injection_style = "NUniformPerCell" |
| 62 | +protons.num_particles_per_cell_each_dim = 4 |
| 63 | +protons.profile = constant |
| 64 | +protons.density = n0 |
| 65 | +protons.momentum_distribution_type = gaussian |
| 66 | +protons.ux_th = sqrt(Ti*q_e/m_p)/clight |
| 67 | +protons.uy_th = sqrt(Ti*q_e/m_p)/clight |
| 68 | +protons.uz_th = sqrt(Ti*q_e/m_p)/clight |
| 69 | +protons.xmin = 0 |
| 70 | +protons.xmax = 10.*de0 |
| 71 | +protons.zmin = 0 |
| 72 | +protons.zmax =10.*de0 |
| 73 | + |
| 74 | +diagnostics.diags_names = diag1 |
| 75 | +diag1.intervals = 500 |
| 76 | +diag1.diag_type = Full |
| 77 | + |
| 78 | +warpx.reduced_diags_names = EP EF |
| 79 | +EP.type = ParticleEnergy |
| 80 | +EP.intervals = 100 |
| 81 | +EF.type = FieldEnergy |
| 82 | +EF.intervals =100 |
0 commit comments