@@ -74,22 +74,18 @@ neutron_1.do_not_deposit = 1
74
74
my_constants.background_dens = 1.e26
75
75
my_constants.beam_dens = 1.e20
76
76
77
+ deuterium_2.species_type = deuterium
78
+ deuterium_2.injection_style = "NRandomPerCell"
79
+ deuterium_2.num_particles_per_cell = 8000
80
+ deuterium_2.profile = "parse_density_function"
77
81
## A tenth of the macroparticles in each cell is made of immobile high-density background deuteriums.
78
82
## The other nine tenths are made of fast low-density beam deuteriums.
79
- deuterium_2.species_type = deuterium
80
- deuterium_2.injection_sources = high_density low_density
81
- deuterium_2.profile = constant
82
- deuterium_2.high_density.injection_style = "NRandomPerCell"
83
- deuterium_2.high_density.num_particles_per_cell = 800
84
- deuterium_2.high_density.density = 10.*background_dens
85
- deuterium_2.high_density.momentum_distribution_type = at_rest
86
- deuterium_2.low_density.injection_style = "NRandomPerCell"
87
- deuterium_2.low_density.num_particles_per_cell = 7200
88
- deuterium_2.low_density.density = 10./9.*beam_dens
89
- deuterium_2.low_density.momentum_distribution_type = "parse_momentum_function"
90
- deuterium_2.low_density.momentum_function_ux(x,y,z) = 0.
91
- deuterium_2.low_density.momentum_function_uy(x,y,z) = 0.
92
- deuterium_2.low_density.momentum_function_uz(x,y,z) = sqrt(2*m_deuterium*Energy_step*(floor(z)**2))/(m_deuterium*clight)
83
+ deuterium_2.density_function(x,y,z) = if(y - floor(y) < 0.1, 10.*background_dens, 10./9.*beam_dens)
84
+ deuterium_2.momentum_distribution_type = "parse_momentum_function"
85
+ deuterium_2.momentum_function_ux(x,y,z) = 0.
86
+ deuterium_2.momentum_function_uy(x,y,z) = 0.
87
+ deuterium_2.momentum_function_uz(x,y,z) = "if(y - floor(y) < 0.1,
88
+ 0., sqrt(2*m_deuterium*Energy_step*(floor(z)**2))/(m_deuterium*clight))"
93
89
deuterium_2.do_not_push = 1
94
90
deuterium_2.do_not_deposit = 1
95
91
0 commit comments