Quasilinear modelling of 2D3V Weibel turbulence

📌 Key Features

• Solves the Vlasov-Maxwell equations in 3D velocity-space and 2D wavevector-space grids.
• Simulates the quasilinear interaction between modes of Weibel turbulence in a weakly collisional, magnetized, anisotropic plasma
• Tracks spectral dynamics and field energy transfer
• Implements Stormer-Verlet (Leapfrog) numerical scheme
• MPI-parallelized for efficient computation of mode dynamics across distributed systems (e.g., KIAM K100 supercomputer).

Compilation and Execution

This code was developed on the KIAM K100 supercomputer (Moscow, Russia).

To compile the program, run:

make

After compilation, run the program with:

mpirun -np <number_of_processes> -maxtime <time_of_execution_in_minutes> build/exe9_c

Configuration

Simulation parameters are set in config.txt. Key parameters include:

• start_level_modes — initial perturbation level for each Weibel mode

• Ngarmonik_r, Ngarmonik_phi — number of filamentation modes in r and phi in polar coordinate

• Beta_perp — thermal velocities of background and beam plasma

• A — relative beam density

• MAGNIT — external uniform magnetic field

• type_of_collisions — 0 - without collisions; 1 - anomolous collisions proportional to b²/(8π)

• setkaBB — velocity space grid size

• ratio_Vstepx_multiple_setkaBB_to_Beta_perp — X-step scaling - Determines velocity grid spacing in x-direction: Δv_x = (β_bkg × 8)/setkaBB

• Tmax maximum simulation time (in plasma frequency)

• dt time step (in plasma frequency)

• kmin_r — Minimum k_r - Lower bound of wavevector range in r-direction

• kstep_r — k_r step - Spacing between wavevectors in r-direction

License

Apache License 2.0 - See LICENSE for details.

❓ Support

For questions and bug reports:

• Open an issue
• Email:
  • kuznetsov.alexey@ipfran.ru
  • a.kuznetsov7677@gmail.com