/vpic_reconnection_tutorial

A tutorial on running reconnection simulations using VPIC

Primary LanguageC++MIT LicenseMIT

Running Reconnection Simulations Using VPIC

The tutorial is based on magnetic reconnection simulations. The simulation decks are in decks. We will assume you can access the NERSC's Cori or Perlmutter or TACC's Frontera, but the procedure can be generalized to other machines without too much additional effort.

NOTE: The old version of VPIC uses binary (N-to-N) output (it is still the default for the master branch). The newer versions now support HDF5 file output. We will recommend using the HDF5 version of VPIC.

Compiling and running

Check out how to compile VPIC HDF5 branch and do a test run at HDF5_version_of_VPIC. After that, you should be able to produce some nice pictures of the reconnection layer.

Modifying the decks

If you decide to modify the reconnection deck, see how to do that at modify_reconnection_deck. Make sure to recompile the deck after you change some variables.

Particle tracking

If you hope to tracer some particles in high cadence to plot some particle trajectories as the simulation is running, checkout particle_tracking.

Particle data

Check out how to dump particle data at modify_reconnection_deck. You can plot some phase space diagram using the particle data. Check out a simple script phase_diagram.py to do that.

Contributing

Please submit issues if you find problems or directly contact me for help.