blob: 7bba3b677815d2d14e1bd736047f96ca85bcc5a7 (
plain) (
blame)
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
|
#################################
####### GENERAL PARAMETERS ######
#################################
max_step = 1000
amr.n_cell = 1024 512
amr.max_grid_size = 128
amr.blocking_factor = 32
amr.max_level = 0
geometry.coord_sys = 0 # 0: Cartesian
geometry.is_periodic = 0 0 # Is periodic?
geometry.prob_lo = -100.e-6 0. # physical domain
geometry.prob_hi = 100.e-6 100.e-6
warpx.verbose = 1
warpx.serialize_ics = 1
#################################
############ NUMERICS ###########
#################################
interpolation.nox = 3
interpolation.noy = 3
interpolation.noz = 3
my_constants.zc = 20.e-6
my_constants.zp = 20.05545177444479562e-6
my_constants.lgrad = .08e-6
my_constants.nc = 1.74e27
my_constants.zp2 = 24.e-6
my_constants.zc2 = 24.05545177444479562e-6
warpx.cfl = 1.0
warpx.use_filter = 1
algo.load_balance_intervals = 66
#################################
############ PLASMA #############
#################################
particles.species_names = electrons ions
electrons.charge = -q_e
electrons.mass = m_e
electrons.injection_style = NUniformPerCell
electrons.num_particles_per_cell_each_dim = 2 2
electrons.momentum_distribution_type = "gaussian"
electrons.ux_th = .01
electrons.uz_th = .01
electrons.zmin = "zc-lgrad*log(400)"
electrons.zmax = 25.47931e-6
electrons.profile = parse_density_function
electrons.density_function(x,y,z) = "(z<zp)*nc*exp((z-zc)/lgrad)+(z>=zp)*(z<=zp2)*2.*nc+(z>zp2)*nc*exp(-(z-zc2)/lgrad)"
ions.charge = q_e
ions.mass = m_p
ions.injection_style = NUniformPerCell
ions.num_particles_per_cell_each_dim = 2 2
ions.momentum_distribution_type = "gaussian"
ions.ux_th = 0.
ions.uz_th = 0.
ions.zmin = 19.520e-6
ions.zmax = 25.47931e-6
ions.profile = parse_density_function
ions.density_function(x,y,z) = "(z<zp)*nc*exp((z-zc)/lgrad)+(z>=zp)*(z<=zp2)*2.*nc+(z>zp2)*nc*exp(-(z-zc2)/lgrad)"
#################################
############# LASER #############
#################################
lasers.names = laser1
laser1.position = 0. 0. 5.e-6 # This point is on the laser plane
laser1.direction = 0. 0. 1. # The plane normal direction
laser1.polarization = 1. 0. 0. # The main polarization vector
laser1.e_max = 4.e12 # Maximum amplitude of the laser field (in V/m)
laser1.wavelength = 0.8e-6 # The wavelength of the laser (in meters)
laser1.profile = Gaussian
laser1.profile_waist = 5.e-6 # The waist of the laser (in meters)
laser1.profile_duration = 15.e-15 # The duration of the laser (in seconds)
laser1.profile_t_peak = 25.e-15 # The time at which the laser reaches its peak (in seconds)
laser1.profile_focal_distance = 15.e-6 # Focal distance from the antenna (in meters)
# Diagnostics
diagnostics.diags_names = diag1
diag1.intervals = 10
diag1.diag_type = Full
|
