diff options
Diffstat (limited to 'Examples/Physics_applications')
3 files changed, 23 insertions, 28 deletions
diff --git a/Examples/Physics_applications/capacitive_discharge/PICMI_inputs_2d.py b/Examples/Physics_applications/capacitive_discharge/PICMI_inputs_2d.py index b9bc0bb3d..171e410b4 100644 --- a/Examples/Physics_applications/capacitive_discharge/PICMI_inputs_2d.py +++ b/Examples/Physics_applications/capacitive_discharge/PICMI_inputs_2d.py @@ -23,7 +23,7 @@ D_CA = 0.067 # m N_INERT = 9.64e20 # m^-3 T_INERT = 300.0 # K -FREQ = 13.56e6 # MHz +FREQ = 13.56e6 # Hz VOLTAGE = 450.0 @@ -292,6 +292,7 @@ mcc_ions = picmi.MCCCollisions( grid = picmi.Cartesian2DGrid( number_of_cells = [nx, ny], + warpx_max_grid_size=128, lower_bound = [xmin, ymin], upper_bound = [xmax, ymax], bc_xmin = 'dirichlet', diff --git a/Examples/Physics_applications/capacitive_discharge/analysis.py b/Examples/Physics_applications/capacitive_discharge/analysis.py index 22d1c725a..1563438bd 100755 --- a/Examples/Physics_applications/capacitive_discharge/analysis.py +++ b/Examples/Physics_applications/capacitive_discharge/analysis.py @@ -14,5 +14,5 @@ import checksumAPI my_check = checksumAPI.evaluate_checksum( 'background_mcc', 'Python_background_mcc_plt00050', - do_particles=True, rtol=7.5e-2 + do_particles=True, rtol=2.5e-3 ) diff --git a/Examples/Physics_applications/capacitive_discharge/inputs_2d b/Examples/Physics_applications/capacitive_discharge/inputs_2d index 5d7194c46..9917d8248 100644 --- a/Examples/Physics_applications/capacitive_discharge/inputs_2d +++ b/Examples/Physics_applications/capacitive_discharge/inputs_2d @@ -1,9 +1,15 @@ # Input file for MCC testing. This runs an eighth of a voltage input period # of the first benchmark case of Turner et al. (2013). +my_constants.Ngas = 9.64e+20 # m^-3 +my_constants.Tgas = 300 # K +my_constants.Nplasma = 2.56e14 # m^-3 +my_constants.freq = 13.56e6 # Hz +my_constants.Mion = 6.67e-27 # kg + max_step = 50 warpx.verbose = 0 -warpx.const_dt = 1.8436578171091445e-10 +warpx.const_dt = 1.0/(400*freq) warpx.do_electrostatic = labframe warpx.self_fields_required_precision = 1e-06 warpx.use_filter = 0 @@ -18,56 +24,44 @@ geometry.prob_hi = 0.067 0.0041875 boundary.field_lo = pec periodic boundary.field_hi = pec periodic -boundary.potential_hi_x = 450.0*sin(2*pi*13.56e6*t) +boundary.potential_hi_x = 450.0*sin(2*pi*freq*t) # Order of particle shape factors algo.particle_shape = 1 particles.species_names = electrons he_ions -electrons.mass = m_e -electrons.charge = -q_e +electrons.species_type = electron electrons.injection_style = nuniformpercell electrons.initialize_self_fields = 0 electrons.num_particles_per_cell_each_dim = 32 16 electrons.profile = constant -electrons.density = 2.56e14 -electrons.momentum_distribution_type = gaussian -electrons.ux_m = 0.0 -electrons.uy_m = 0.0 -electrons.uz_m = 0.0 -electrons.ux_th = 0.0022492462021094224 -electrons.uy_th = 0.0022492462021094224 -electrons.uz_th = 0.0022492462021094224 +electrons.density = Nplasma +electrons.momentum_distribution_type = maxwell_boltzmann +electrons.theta = (kb*30000/(m_e*clight^2)) -he_ions.mass = 6.67e-27 +he_ions.species_type = helium he_ions.charge = q_e he_ions.injection_style = nuniformpercell he_ions.initialize_self_fields = 0 he_ions.num_particles_per_cell_each_dim = 32 16 he_ions.profile = constant -he_ions.density = 2.56e14 -he_ions.momentum_distribution_type = gaussian -he_ions.ux_m = 0.0 -he_ions.uy_m = 0.0 -he_ions.uz_m = 0.0 -he_ions.ux_th = 2.6332116917541573e-06 -he_ions.uy_th = 2.6332116917541573e-06 -he_ions.uz_th = 2.6332116917541573e-06 +he_ions.density = Nplasma +he_ions.momentum_distribution_type = maxwell_boltzmann +he_ions.theta = (kb*Tgas/(Mion*clight^2)) collisions.collision_names = coll_elec coll_ion coll_ion.type = background_mcc coll_ion.species = he_ions -coll_ion.background_density = 9.64e+20 -coll_ion.background_temperature = 300 +coll_ion.background_density = Ngas +coll_ion.background_temperature = Tgas coll_ion.scattering_processes = elastic back coll_ion.elastic_cross_section = ../../../../warpx-data/MCC_cross_sections/He/ion_scattering.dat coll_ion.back_cross_section = ../../../../warpx-data/MCC_cross_sections/He/ion_back_scatter.dat coll_elec.type = background_mcc coll_elec.species = electrons -coll_elec.background_density = 9.64e+20 -coll_elec.background_temperature = 300 -coll_elec.background_mass = 6.67e-27 +coll_elec.background_density = Ngas +coll_elec.background_temperature = Tgas coll_elec.scattering_processes = elastic excitation1 excitation2 ionization coll_elec.elastic_cross_section = ../../../../warpx-data/MCC_cross_sections/He/electron_scattering.dat coll_elec.excitation1_energy = 19.82 |