Add 2D circle EB test (BLAST-WarpX#2538)

* Added embedded_circle test

* Add embedded_circle test files

* Removed diag files

* removed PICMI input file

* Update to use default regression analysis

* Added line breaks for spacing

* Added description

* Fixed benchmark file

* Added load balancing to test

* Commented out load_balancing portion of test.
This will be added back in once load balancing is fixed.

Co-authored-by: Axel Huebl <[email protected]>

Author: 2 people

Examples/Tests/embedded_circle/inputs_2d

@@ -0,0 +1,92 @@
+# This is a 2D electrostatic simulation containing a circle embedded boundary at
+# the center of the grid domain, pre -seeded with equal density of argon ions
+# and electrons.
+
+max_step = 11
+warpx.const_dt = 3.99e-13
+warpx.do_electrostatic = labframe
+warpx.self_fields_required_precision = 1e-06
+warpx.eb_implicit_function = -((x-0.00005)**2+(z-0.00005)**2-1e-05**2)
+warpx.eb_potential(x,y,z,t) = -10
+warpx.self_fields_absolute_tolerance = 0.02
+
+# algo.load_balance_intervals = 5
+# algo.load_balance_efficiency_ratio_threshold = 1.001
+# algo.load_balance_with_sfc = 0
+# algo.load_balance_knapsack_factor = 2
+
+amr.n_cell = 128 128
+amr.max_grid_size = 16
+amr.max_level = 0
+geometry.coord_sys = 0
+geometry.prob_lo = 0.0 0.0
+geometry.prob_hi = 0.0001 0.0001
+boundary.field_lo = pec periodic
+boundary.field_hi = pec periodic
+boundary.particle_lo = absorbing periodic
+boundary.particle_hi = absorbing periodic
+algo.particle_shape = 1
+
+particles.species_names = electrons ar_ions
+electrons.mass = m_e
+electrons.charge = -q_e
+electrons.injection_style = nuniformpercell
+electrons.initialize_self_fields = 1
+electrons.num_particles_per_cell_each_dim = 2 2
+electrons.density = 10.2e16
+electrons.momentum_distribution_type = gaussian
+electrons.profile = constant
+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.save_particles_at_zlo = 1
+electrons.save_particles_at_zhi = 1
+electrons.save_particles_at_eb = 1
+ar_ions.mass = 6.633520326333601e-26
+ar_ions.charge = q_e
+ar_ions.injection_style = nuniformpercell
+ar_ions.initialize_self_fields = 1
+ar_ions.num_particles_per_cell_each_dim = 2 2
+ar_ions.profile = constant
+ar_ions.density = 10.2e16
+ar_ions.momentum_distribution_type = gaussian
+ar_ions.ux_m = 0.0
+ar_ions.uy_m = 0.0
+ar_ions.uz_m = 0.0
+ar_ions.ux_th = 2.6285641070031447e-06
+ar_ions.uy_th = 2.6285641070031447e-06
+ar_ions.uz_th = 2.6285641070031447e-06
+ar_ions.save_particles_at_zlo = 1
+ar_ions.save_particles_at_zhi = 1
+ar_ions.save_particles_at_eb = 1
+
+collisions.collision_names = coll_electrons coll_ar_ions
+coll_electrons.type = background_mcc
+coll_electrons.species = electrons
+coll_electrons.background_density = 2.1458888524486985e+23
+coll_electrons.background_temperature = 450.00000000000006
+coll_electrons.background_mass = 6.633520326333601e-26
+coll_electrons.scattering_processes = excitation1 elastic ionization
+coll_electrons.excitation1_cross_section = ../../../../warpx-data/MCC_cross_sections/Ar/excitation_1.dat
+coll_electrons.excitation1_energy = 11.5
+coll_electrons.elastic_cross_section = ../../../../warpx-data/MCC_cross_sections/Ar/electron_scattering.dat
+coll_electrons.ionization_cross_section = ../../../../warpx-data/MCC_cross_sections/Ar/ionization.dat
+coll_electrons.ionization_energy = 15.7596112
+coll_electrons.ionization_species = ar_ions
+coll_ar_ions.type = background_mcc
+coll_ar_ions.species = ar_ions
+coll_ar_ions.background_density = 2.1458888524486985e+23
+coll_ar_ions.background_temperature = 450.00000000000006
+coll_ar_ions.scattering_processes = back elastic charge_exchange
+coll_ar_ions.back_cross_section = ../../../../warpx-data/MCC_cross_sections/Ar/ion_back_scatter.dat
+coll_ar_ions.elastic_cross_section = ../../../../warpx-data/MCC_cross_sections/Ar/ion_scattering.dat
+coll_ar_ions.charge_exchange_cross_section = ../../../../warpx-data/MCC_cross_sections/Ar/charge_exchange.dat
+
+diagnostics.diags_names = diag1
+diag1.diag_type = Full
+diag1.format = plotfile
+diag1.intervals = 1
+diag1.fields_to_plot = phi rho_electrons rho_ar_ions

Regression/Checksum/benchmarks_json/embedded_circle.json

@@ -0,0 +1,27 @@
+{
+  "ar_ions": {
+    "particle_cpu": 31740.0,
+    "particle_id": 3219770578.0,
+    "particle_momentum_x": 2.661332598290843e-18,
+    "particle_momentum_y": 2.6484537197836753e-18,
+    "particle_momentum_z": 2.6446763485694346e-18,
+    "particle_position_x": 3.174118714858889,
+    "particle_position_y": 3.174074331596828,
+    "particle_weight": 988031616.2109375
+  },
+  "electrons": {
+    "particle_cpu": 30737.0,
+    "particle_id": 1041042475.0,
+    "particle_momentum_x": 3.014851468379113e-20,
+    "particle_momentum_y": 3.016132024233842e-20,
+    "particle_momentum_z": 3.0169481587733156e-20,
+    "particle_position_x": 3.0760407489559514,
+    "particle_position_y": 3.072305059445631,
+    "particle_weight": 956794738.7695312
+  },
+  "lev=0": {
+    "phi": 56888.72655917703,
+    "rho_ar_ions": 257.78992245244984,
+    "rho_electrons": 250.2854445348176
+  }
+}

Regression/WarpX-tests.ini

@@ -2245,6 +2245,24 @@ compareParticles = 0
 analysisRoutine = Examples/Tests/ElectrostaticSphere/analysis_electrostatic_sphere.py
 tolerance = 1.e-12
 
+[embedded_circle]
+buildDir = .
+inputFile = Examples/Tests/embedded_circle/inputs_2d
+runtime_params =
+dim = 2
+addToCompileString = USE_EB=TRUE
+restartTest = 0
+useMPI = 1
+numprocs = 2
+useOMP = 1
+numthreads = 1
+compileTest = 0
+doVis = 0
+compareParticles = 1
+particleTypes = electrons ar_ions
+analysisRoutine = Examples/analysis_default_regression.py
+tolerance = 1.e-12
+
 [initial_distribution]
 buildDir = .
 inputFile = Examples/Tests/initial_distribution/inputs