nbia-astro / skeletor

Parallel PIC code written in Python and based on the skeleton codes provided by PICKSC
GNU General Public License v3.0
1 stars 0 forks source link

Shearing wave ideas #170

Open tberlok opened 7 years ago

tberlok commented 7 years ago

This is a summary of our meeting about the shearing wave. We were encouraged by the following figure: shearing-wave.pdf which looks very much like figure 4 in Toby's paper, but there is still an issue with something that looks like a numerical instability.

Our ideas are listed below in no particular order.

  1. Fix the fix-deposit function in sources.py such that it only erases guard layers on the physical boundaries. It will then work with MPI. I fixed this right now as it was easy.
  2. Run the simulation at varying resolution and time step with everything else fixed.
  3. Figure out if we can initialize a warm plasma with dQ=0.2 or if it's too noisy. Increasing to dQ=2 will have the nonlinear damping start sooner.
  4. Extend finite grid theory to drifting plasma. Test this for a cold, drifting ionacoustic wave using our code.
  5. Try solving the equations for the relative velocities to see if the problem persists.
tobson commented 7 years ago

@tberlok Do you have notes on the finite grid linear theory up on GitHub?