Description
Discussed in #5538
Originally posted by moisesae January 6, 2025
I have taken the <u_x> moment using the particle_fields.do_average feature in the input file to get the hydrodynamic flow profile that WarpX interpolates onto the grid which I can then analyze. I have been taking numerical gradients of this profile and have noticed that I consistently get periodic spikes showing up. I thought that if I did a time-average of the velocity profile (~1000 time-steps), I would be able to see that the spikes go away when I apply the gradient. However, they persist even with time-averaging and I am having a difficult time trying to figure out why. They seem to occur every ~1mm but there isn't a length-scale that would correspond to this such as the debye length, gyro-radius, etc, which are sufficiently resolved as well (time-scales too). Pre-gradient, the spikes look small but once you apply the gradient, they get significantly amplified. I can see this in all of the species being simulated, with and without the B-field.
I'm not sure if it is an issue of how the moments are getting interpolated onto the grid, but I would appreciate if anyone could point me to some sanity-checks I can try to pinpoint this down.
(Each curve is normalized by their maximum value so they can be plotted on the same scale and be dimensionless)
The following is an example where I used significant time-averaging but the periodic spikes persist.
