Skip to main content
Sandbox Physics
Kinetic plasma · A course with experiments

L2 · Chapter 3

Where a fading ripple hides, and how a grid brings it back

The density ripple fades because the distribution folds into stripes across velocity, and those stripes get finer at a steady rate. A computer stores the stripes on a velocity grid. When they become finer than one grid cell, the grid misreads them as coarse stripes, and the ripple reappears at a time set only by the grid.

Derivation

The stripes tighten at a fixed rate

Write the ripple as the complex mode eikxe^{ikx}, keeping in mind that the physical ripple is its real part. Under free streaming, eikxe^{ikx} becomes eik(x−vt)=eikxe−iktve^{ik(x-vt)}=e^{ikx}e^{-iktv}. At a fixed position, the distribution now oscillates in velocity with wave number κ=kt\kappa = kt: the stripes tighten in proportion to time.

Density integrates over velocity. Integrating a smooth Maxwellian against a fast oscillation gives almost zero, and the faster the oscillation, the closer to zero. That is phase mixing. With the field on, the electrons also push on each other, but the stripes keep tightening at very nearly the same rate.

eik(x−vt)=eikx e−iκv,κ(t)=kte^{ik(x-vt)}=e^{ikx}\,e^{-i\kappa v},\qquad \kappa(t)=kt

A grid can only hold stripes so fine

A velocity grid with spacing Δv\Delta v samples e−iκve^{-i\kappa v} only at vj=jΔvv_j=j\Delta v. When κΔv\kappa\Delta v reaches 2π2\pi, every sample equals its value at κ=0\kappa=0. The grid then holds the same numbers as the original smooth ripple, the velocity integral no longer cancels, and the field returns.

The return time is the recurrence time. It depends on the wave number and the velocity spacing, and on nothing physical. Because it arrives on a schedule set by the grid, refining the velocity grid is the clean test: a recurrence moves later, while a physical signal stays put.

k TR Δv=2π⟹TR=2πk Δvk\,T_R\,\Delta v=2\pi\qquad\Longrightarrow\qquad T_R=\frac{2\pi}{k\,\Delta v}

Worked example

Worked example: two grids, one ripple

The phase-space Lab spans velocities from −6-6 to 66. With 256256 cells, Δv=12/256≈0.047\Delta v=12/256\approx0.047, and a k=0.5k=0.5 ripple recurs at TR=2π/(0.5×0.047)≈268T_R=2\pi/(0.5\times0.047)\approx268, far beyond the 70-unit run. At t=40t=40 one stripe spans 2π/(kt)≈0.312\pi/(kt)\approx0.31 in velocity, still about seven cells. The grid resolves it.

With 3232 cells, Δv=0.375\Delta v=0.375 and TR≈33.5T_R\approx33.5. Run the coarse-grid preset and the field climbs back near t=33.5t=33.5, exactly where the dashed recurrence line sits on the history chart. The damping rate measured before the recurrence is still −0.1536-0.1536, the same as on the fine grid, because the early fit window ends before the grid fails.

Three things that can look like a returning signal

A returning field can be physics, a grid recurrence or a reflection from the edge of the velocity box. The Lab’s velocity box is periodic, so a stripe that leaves at the top re-enters at the bottom; that does no harm while the distribution is negligible near the limits, which is why the Maxwellian runs use limits of ±6\pm6 and the beam runs ±9\pm9.

Each suspect has its own test. Doubling the number of velocity cells doubles a recurrence time. Widening the velocity box at fixed spacing changes a boundary effect. A physical signal survives both. The project benchmarks apply both tests to the echo of chapter 6: with 128128, 256256 and 512512 cells, and with limits of 66 and 1212, the echo arrives at the same time.

An experiment to think through

Predict, change, explain

  1. Predict the recurrence time for a k=0.5k=0.5 ripple on a 64-cell velocity grid.
  2. Choose 64 velocity cells in the phase-space Lab and watch for the field to return. Then try 128.
  3. Use the stripe spectrogram to explain the return: the measured stripe power follows the gold line until the stripes reach the grid limit.

The three phase-space Labs solve a one-dimensional electron model in plasma units. Each Lab exports its parameters, model version and field history.

Try explaining it yourself

1. What recurrence time do you predict for k=0.5k=0.5 and 64 velocity cells over [−6,6][-6,6]?

Read an explained answer

Δv=12/64=0.1875\Delta v=12/64=0.1875, so TR=2π/(0.5×0.1875)≈67T_R=2\pi/(0.5\times0.1875)\approx67. That falls just inside the 70-unit run.

2. Why does the recurrence happen at a single time instead of slowly?

Read an explained answer

Every velocity sample returns to its starting phase at the same moment, because kTRvj=2πjkT_R v_j=2\pi j for all grid points. The cancellation that hid the ripple undoes itself everywhere at once.

3. A colleague reports a revived field at t=67t=67 on a 64-cell grid and calls it an echo. What one test settles it?

Read an explained answer

Rerun with 128 cells. A recurrence moves to t≈134t\approx134; a physical revival stays at t=67t=67.

Sources to read alongside

  1. R. Fitzpatrick, Plasma Physics lecture notes, University of Texas at Austin: Landau damping, resonant particles and the plasma dispersion function

    Three linked pages derive the linear response of a Maxwellian plasma, the role of electrons near the wave speed, and the plasma dispersion function used for our dielectric roots.

  2. C. Z. Cheng and G. Knorr, The integration of the Vlasov equation in configuration space, Journal of Computational Physics 22, 330 (1976)

    The usual reference for splitting a Vlasov step into position and velocity shifts, the approach our solver takes. Only the bibliographic record was checked; the solver is tested against exact limits and conservation laws instead.

FeedbackSomething broken, confusing, or missing on this page?Tell us
What kind of feedback?

Sent with this page’s address, its Lab ID, the site language and your device type. Please leave out personal details. Feedback is kept for up to 12 months.Privacy policy