Ecosyste.ms: Awesome
An open API service indexing awesome lists of open source software.
https://github.com/cchandre/guiding-center
Guiding-center dynamics in plasma physics
https://github.com/cchandre/guiding-center
gyrokinetic hamiltonian numpy plasma-physics python3 scipy sympy
Last synced: about 1 month ago
JSON representation
Guiding-center dynamics in plasma physics
- Host: GitHub
- URL: https://github.com/cchandre/guiding-center
- Owner: cchandre
- License: bsd-2-clause
- Created: 2021-07-05T13:47:11.000Z (over 3 years ago)
- Default Branch: main
- Last Pushed: 2024-08-17T08:07:09.000Z (4 months ago)
- Last Synced: 2024-08-17T09:24:38.977Z (4 months ago)
- Topics: gyrokinetic, hamiltonian, numpy, plasma-physics, python3, scipy, sympy
- Language: Python
- Homepage:
- Size: 6.19 MB
- Stars: 3
- Watchers: 1
- Forks: 0
- Open Issues: 0
-
Metadata Files:
- Readme: README.md
- License: LICENSE
Awesome Lists containing this project
README
# Guiding-Center (GC) dynamics in plasma physics
- [`nondimensionalization.mlx`](https://github.com/cchandre/Guiding-Center/blob/main/nondimensionalization.mlx): Matlab Live Script for the computation of the dimensionless parameters *A*, *ρ* and *η* used in the code and in the analysis
- [`gc2d_dict.py`](https://github.com/cchandre/Guiding-Center/blob/main/gc2d_dict.py): to be edited to change the parameters of the GC computation (see below for a dictionary of parameters)
- [`gc2d.py`](https://github.com/cchandre/Guiding-Center/blob/main/gc2d.py): contains the GC classes and main functions defining the GC dynamics
- [`gc2d_modules.py`](https://github.com/cchandre/Guiding-Center/blob/main/gc2d_modules.py): contains the methods to integrate the GC dynamics
Once [`gc2d_dict.py`](https://github.com/cchandre/Guiding-Center/blob/main/gc2d_dict.py) has been edited with the relevant parameters, run the file as
```sh
python3 gc2d.py
```
or
```sh
nohup python3 -u gc2d.py &>gc2d.out < /dev/null &
```
The list of Python packages and their version are specified in [`requirements.txt`](https://github.com/cchandre/Guiding-Center/blob/main/requirements.txt)
___
## Parameter dictionary- *Potential*: string; 'KMdCN' or 'turbulent'
- *Method*: string
- 'potentials' (only for Potential='turbulent'): plots the electrostatic potential as well as the first and second order guiding-center potentials
- 'diffusion_fo': computes the diffusion coefficient for the full orbits
- 'diffusion_gc': computes the diffusion coefficient for the guiding centers
- 'poincare_fo': plots the full orbits in the plane (*x*, *y*) for every period of the potential (stroboscopic plot)
- 'poincare_gc': plots the guiding-center trajectories in the plane (*x*, *y*) for every period of the potential (stroboscopic plot)
####
- *FLR*: tuple of 2 strings; 'all', 'pade' or 'none'; if 'all', FLR to all orders is taken into account; if 'pade', a Padé approximant is considered for the FLR effects; if 'none', no FLR effects are taken into account
####
- *A*: float or array of floats; amplitude(s) of the electrostatic potential [theory: *A*=εδ/*B*]
- *rho*: float or array of floats; value(s) of the Larmor radius; for full orbits, this value corresponds to the thermal Larmor radius
- *eta*: float or array of floats; value(s) of the coefficient in front of the GC order 2 potential; η>0 for positive charge, η<0 for negative charge [theory: η=1/(2Ω)]
####
- *Ntraj*: integer; number of trajectories to be integrated
- *Tf*: integer; number of periods for the integration of the trajectories
- *threshold*: float; value used to discriminate between trapped and untrapped trajectories (recommended: 4)
- *TwoStepIntegration*: boolean; if True, computes trajectories from 0 to 2π*T*mid, removes the trapped trajectories, and continues integration from 2π*T*mid to 2π*T*f
- *Tmid*: integer; number of periods for the integration of trajectories in the first step (if *TwoStepIntegration*=True)
- *TimeStep*: float; time step used by the integrator (recommended: 5x10-3 for guiding centers and 5x10-4 for full orbits)
- *check_energy*: boolean; if True, the autonomous system is integrated, and the output (`.mat` file) includes the total energy (only if *SaveData*=True)
- *init*: string; 'random' or 'fixed'; method to generate initial conditions
####
- *SaveData*: boolean; if True, the results are saved in a `.mat` file; Poincaré sections and diffusion plots *r*2(*t*) are saved as *fig_extension* files; NB: the diffusion data are saved in a `.txt` file regardless of the value of *SaveData*
- *PlotResults*: boolean; if True, the results are plotted right after the computation
- *Parallelization*: tuple (boolean, int); True for parallelization, int is the number of cores to be used or int='all' to use all available cores
####
- *modulo*: boolean; if True, *x* and *y* are represented modulo 2π (only for Method='poincare' and PlotResults=True)
- *grid*: boolean; if True, show the grid lines on plots
- *darkmode*: boolean; if True, plots are done in dark mode
- *fig_extension*: string; e.g., '.png', '.pdf', '.svg'; format of the figures to be saved
####
- *M*: integer; number of modes (default = 5 for 'KMdCN' and 25 for 'turbulent')
- *N*: integer; number of points on each axis for 'turbulent' (recommended: 212)---
Reference:
- M. Stanzani, F. Arlotti, G. Ciraolo, X. Garbet, C. Chandre, *Transition to super-diffusive transport in turbulent plasmas*, Physical Review E **110**, 025204 (2024); [arXiv:2309.02461](https://arxiv.org/abs/2309.02461)
```bibtex
@article{PhysRevE.110.025204,
title = {Transition to superdiffusive transport in turbulent plasmas},
author = {Stanzani, Matteo and Arlotti, Filippo and Ciraolo, Guido and Garbet, Xavier and Chandre, Cristel},
journal = {Phys. Rev. E},
volume = {110},
issue = {2},
pages = {025204},
numpages = {8},
year = {2024},
month = {Aug},
publisher = {American Physical Society},
doi = {10.1103/PhysRevE.110.025204},
url = {https://link.aps.org/doi/10.1103/PhysRevE.110.025204}
}
```
For more information: