Skip to content

Yee+RZ laser in vacuum propagation - different results for differently polarised laser beams #4436

Open
@SchroederSa

Description

@SchroederSa

When simulating the laser propagation in vacuum using the Yee solver in RZ geometry, I observe:

  • a0 shows different results for x and y polarised beams.
  • a0 for y-polarised beams shows unreasonable fluctuations.
  • Laser envelope for y-polarised beam shows an unexpected high-frequency signal.
  • Higher azimuthal mode resolves the discrepancy in results for differently polarised beams.

By using different a0 retrieval methods, I intended to rule out a diagnostic-related problem to the best of my knowledge.

Best
Sarah

input script

max_step = 300000.0

stop_time = 1.5e-11

warpx.n_rz_azimuthal_modes = 2.0

amr.n_cell = 700    4000
amr.blocking_factor = 32

warpx.numprocs = 2.0    2.0

amr.max_level = 0

warpx.verbose = 1
warpx.do_dive_cleaning = 0
warpx.use_filter = 0
warpx.filter_npass_each_dir = 0.0    1.0
warpx.cfl = 1.0
warpx.do_moving_window = 1
warpx.moving_window_dir = z
warpx.moving_window_v = 1.0
warpx.start_moving_window_step = 1600.0

algo.particle_shape = 3
algo.maxwell_solver = Yee

geometry.dims = RZ
geometry.prob_lo = 0.0    -8e-05
geometry.prob_hi = 0.000115    0.0

boundary.field_lo = none    pec
boundary.field_hi = absorbing_silver_mueller    pec
boundary.particle_lo = absorbing    absorbing
boundary.particle_hi = absorbing    absorbing

lasers.names = laser
laser.do_continuous_injection = 0
laser.profile = Gaussian
laser.direction = 0.0    0.0    1.0
laser.polarization = 1.0    0.0    0.0
laser.position = 0.0    0.0    -1e-06
laser.a0 = 1.36
laser.profile_focal_distance = 0.0022
laser.profile_t_peak = 2e-13
laser.wavelength = 8e-07
laser.profile_waist = 3.3e-05
laser.profile_duration = 3.4e-14

diagnostics.diags_names = diag1
diag1.intervals = 5000
diag1.diag_type = Full
diag1.dump_rz_modes = 1.0
diag1.fields_to_plot = Er    Et    Ez    rho
diag1.file_prefix = diag_vacuum_x_nom/openPMD
diag1.format = openpmd
diag1.openpmd_backend = h5
diag1.electrons.plot_filter_function(t,x,y,z,ux,uy,uz) = "((x*x+y*y)<20e-6)"

VacuumSimulations-compressed.pdf

Metadata

Metadata

Labels

bugSomething isn't workingbug: affects latest releaseBug also exists in latest release versioncomponent: coreCore WarpX functionalitygeometry: RZaxisymmetric 2D and quasi-3D

Type

No type

Projects

No projects

Milestone

No milestone

Relationships

None yet

Development

No branches or pull requests

Issue actions