diff options
Diffstat (limited to 'Exec/test_particle_injection')
| -rw-r--r-- | Exec/test_particle_injection/ParticleProb.cpp | 109 | ||||
| -rw-r--r-- | Exec/test_particle_injection/inputs | 29 |
2 files changed, 20 insertions, 118 deletions
diff --git a/Exec/test_particle_injection/ParticleProb.cpp b/Exec/test_particle_injection/ParticleProb.cpp deleted file mode 100644 index 7cdee1d46..000000000 --- a/Exec/test_particle_injection/ParticleProb.cpp +++ /dev/null @@ -1,109 +0,0 @@ - -// -// Each problem must have its own version of PhysicalParticleContainer::InitData() -// to initialize the particle data. It must also initialize charge and mass. -// - -#include <cmath> - -#include <AMReX_BLProfiler.H> - -#include <ParticleContainer.H> -#include <WarpXConst.H> - -using namespace amrex; - -void -PhysicalParticleContainer::InitData() -{ - BL_PROFILE("PPC::InitData()"); - - charge = -PhysConst::q_e; - mass = PhysConst::m_e; - - const int lev = 0; - - const Geometry& geom = Geom(lev); - const Real* dx = geom.CellSize(); - - Real weight, ux, uy, uz; - Real particle_xmin, particle_xmax, particle_ymin, particle_ymax, particle_zmin, particle_zmax; - int n_part_per_cell; - { - ParmParse pp("langmuirwave"); - n_part_per_cell = 1; - pp.query("num_particles_per_cell", n_part_per_cell); - weight = 1.e25; - pp.query("n_e", weight); - #if BL_SPACEDIM==3 - weight *= dx[0]*dx[1]*dx[2]/n_part_per_cell; - #elif BL_SPACEDIM==2 - weight *= dx[0]*dx[1]/n_part_per_cell; - #endif - - particle_xmin = particle_ymin = particle_zmin = -2.e-5; - particle_xmax = particle_ymax = particle_zmax = 2.e-5; - pp.query("particle_xmin", particle_xmin); - pp.query("particle_xmax", particle_xmax); - pp.query("particle_ymin", particle_ymin); - pp.query("particle_ymax", particle_ymax); - pp.query("particle_zmin", particle_zmin); - pp.query("particle_zmax", particle_zmax); - - ux = 0.; - uy = 0.; - uz = 0.; - pp.query("ux", ux); - pp.query("uy", uy); - pp.query("uz", uz); - - ux *= PhysConst::c; - uy *= PhysConst::c; - uz *= PhysConst::c; - } - - std::array<Real,PIdx::nattribs> attribs; - attribs.fill(0.0); - attribs[PIdx::w ] = weight; - attribs[PIdx::ux] = ux; - attribs[PIdx::uy] = uy; - attribs[PIdx::uz] = uz; - - for (MFIter mfi = MakeMFIter(lev); mfi.isValid(); ++mfi) - { - const Box& tile_box = mfi.tilebox(); - RealBox tile_real_box { tile_box, dx, geom.ProbLo() }; - - const int grid_id = mfi.index(); - const int tile_id = mfi.LocalTileIndex(); - - const IntVect& boxlo = tile_box.smallEnd(); - for (IntVect iv = boxlo; iv <= tile_box.bigEnd(); tile_box.next(iv)) - { - for (int i_part=0; i_part<n_part_per_cell;i_part++) - { - Real particle_shift = (0.5+i_part)/n_part_per_cell; -#if (BL_SPACEDIM == 3) - Real x = tile_real_box.lo(0) + (iv[0]-boxlo[0] + particle_shift)*dx[0]; - Real y = tile_real_box.lo(1) + (iv[1]-boxlo[1] + particle_shift)*dx[1]; - Real z = tile_real_box.lo(2) + (iv[2]-boxlo[2] + particle_shift)*dx[2]; -#elif (BL_SPACEDIM == 2) - Real x = tile_real_box.lo(0) + (iv[0]-boxlo[0] + particle_shift)*dx[0]; - Real y = 0.0; - Real z = tile_real_box.lo(1) + (iv[1]-boxlo[1] + particle_shift)*dx[1]; -#endif - - if (x >= particle_xmax || x < particle_xmin || - y >= particle_ymax || y < particle_ymin || - z >= particle_zmax || z < particle_zmin ) - { - continue; - } - else - { - AddOneParticle(lev, grid_id, tile_id, x, y, z, attribs); - } - } - } - } -} diff --git a/Exec/test_particle_injection/inputs b/Exec/test_particle_injection/inputs index 26e88f6c6..6681cbb3c 100644 --- a/Exec/test_particle_injection/inputs +++ b/Exec/test_particle_injection/inputs @@ -31,12 +31,23 @@ algo.particle_pusher = 0 # CFL warpx.cfl = 1.0 -langmuirwave.particle_xmin = -20.e-6 -langmuirwave.particle_xmax = 0.e-6 -langmuirwave.particle_ymin = -20.e-6 -langmuirwave.particle_ymax = 20.e-6 -langmuirwave.particle_zmin = -20.e-6 -langmuirwave.particle_zmax = 20.e-6 -langmuirwave.num_particles_per_cell = 10 -langmuirwave.n_e = 1.e25 # number of electrons per m^3 -langmuirwave.ux = 0.01 # ux = gamma*beta_x +particles.nspecies = 1 +particles.species_names = electrons + +electrons.charge = -q_e +electrons.mass = m_e +electrons.injection_style = "NDiagPerCell" + +electrons.xmin = -20.e-6 +electrons.xmax = 0.e-6 +electrons.ymin = -20.e-6 +electrons.ymax = 20.e-6 +electrons.zmin = -20.e-6 +electrons.zmax = 20.e-6 + +electrons.num_particles_per_cell = 10 +electrons.profile = constant +electrons.density = 1.e25 # number of electrons per m^3 + +electrons.momentum_distribution_type = "constant" +electrons.ux = 0.01 # ux = gamma*beta_x |