Readme of available codes (http://hsxie.me/codes/) by Hua-sheng XIE:
- amc
AMC (Alfven Mode Code) is an eigenvalue code mainly (but not limited) aimed to study the Alfven physics (continuum spectrums and eigenmodes) in tokamak.
Cite:
[Xie2015] H.S. Xie & Y. Xiao, Parallel equilibrium current effect on existence of reversed shear Alfven eigenmodes, Phys. Plasmas, 22, 022518 (2015). http://scitation.aip.org/content/aip/journal/pop/22/2/10.1063/1.4913487
AMC has been used in [Chen2014a] W. Chen et al, EPL 107, 25001 (2014).
----
Also 'awcon', which calculates only the continuums, and can also calculates acoustic coupling.
Awcon has been used in [Chen2014b] W. Chen et al, Nucl. Fusion 54, 104002 (2014).
Hua-sheng XIE (IFTS-ZJU, huashengxie@gmail.com)
20:56 2015/3/4
===========================================================
Also
1. GFLDR
Numerical Solving of General Fishbone-Like Dispersion Relation in the Accumulating Point
This code is used and updated in [Ma2014].
cite:
[Ma2014] Ruirui Ma, Ilija Chavdarovski, Gaoxiang Ye, and Xin Wang, Linear dispersion relation of beta-induced Alfv¨¦n eigenmodes in presence of anisotropic energetic ions, Physics of Plasmas 21, 062120 (2014).
2. tearing
Linear (both intial value and eigenvalue solvers) tearing mode solution in 2D reduced MHD equation
a. Slab 2D
b. Cylinder
The cylinder version has been used with update in [Liu2015].
Cite:
[Liu2015] Dongjian Liu, Wenlu Zhang, Joseph McClenaghan, Jiaqi Wang, and Zhihong Lin, Verification of gyrokinetic particle simulation of current-driven instability in fusion plasmas. II. Resistive tearing mode, Physics of Plasmas 21, 122520 (2014).
20:58 2015/3/20
- gpdf
Cite:
[Xie2013] H. S. Xie, Generalized plasma dispersion function: One-solve-all treatment, visualizations, and application to Landau damping, Physics of Plasmas, 2013, 20, 092125, http://scitation.aip.org/content/aip/journal/pop/20/9/10.1063/1.4822332
% Hua-sheng XIE, huashengxie@gmail.com, IFTS-ZJU, 2013-05-26 23:37
% Calculate GPDF, Generalized Plasma Dispersion Function, see [Xie2013]
% Modify the parameters in the code for your own usage
Abstract
A unified, fast, and effective approach is developed for numerical
calculation of the well-known plasma dispersion function with
extensions from Maxwellian distribution to almost arbitrary
distribution functions, such as the $\delta$, flat top, triangular,
$\kappa$ or Lorentzian, slowing down, and incomplete Maxwellian
distributions. The singularity and analytic continuation problems
are also solved generally. Given that the usual conclusion
$\gamma\propto\partial f_0/\partial v$ is only a rough approximation
when discussing the distribution function effects on Landau damping,
this approach provides a useful tool for rigorous calculations of
the linear wave and instability properties of plasma for general
distribution functions. The results are also verified via a linear
initial value simulation approach. Intuitive visualizations of the
generalized plasma dispersion function are also provided.
Download: See supplementary material at http://dx.doi.org/10.1063/1.4822332 E-PHPAEN-20-024310 for the MATLAB version numerical routines, where GPDF, root finding, IVP simulation and test plotting codes are included. [ Supplementary Material]
ftp://ftp.aip.org/epaps/phys_plasmas/E-PHPAEN-20-024310
Also: http://arxiv.org/abs/1305.6476, http://arxiv.org/abs/1211.5984, http://arxiv.org/abs/1304.5883
Hua-sheng XIE (huashengxie@gmail.com, IFTS-ZJU)
14:22 2015/3/5
- landaudamping
Landau Damping Simulation Models
Hua-sheng Xie, IFTS-ZJU, huashengxie@gmail.com
Abstract:
Several models for Landau damping simulations are provided: the simplest linear model, particle-in-cell (PIC) simulation and Vlasov continuity simulation as well as the discussions of the numerical aspects of the dispersion relation.
I hope this tutorial can be useful for beginners to understand the basic physics and numerical methods in plasma physics. It can be used for the undergraduate or graduate courses such as "Introduction to Plasma Physics" or "Computational Plasma Physics". And, it can also be the start point of your researches.
This version of the tutorial (notes and codes) was presented at Sichuan University, Chengdu, China.
You can modified/improved them for your own usages. But, you should keep the copyright information. This statement also holds for many of my other codes.
xiehuasheng
2013-10-15 14:48
- orbitm
1. orbit_general
Orbit.m, single particle motion for: space plasma (current sheet, dipole, ...), Tokamak, Laser, ...
A short code, but can treat almost all the common particle motion problems in different EM fields.
2. tokamak
Orbit_tokamak: Using Lorentz force, with gyro. Rewritten from "Orbit.m".
OrbitGC_qfix: Guiding center orbit with q=const.
OrbitGC_qpsip: Guiding center orbit with q=q1+q2*psip+q3*psip^2.
See "OrbitTokamak_doc.doc" for related formula.
Hua-sheng XIE (huashengxie@gmail.com)
10:32 2012-12-28
- pdrf
Cite:
[Xie2014] H. S. Xie, PDRF: A general dispersion relation solver for magnetized multi-fluid plasma, Computer Physics Communications, 2014, 185, 670 - 675, http://www.sciencedirect.com/science/article/pii/S0010465513003408.
% Hua-sheng XIE, huashengxie@gmail.com, IFTS-ZJU, 2013-05-03 22:02
% pdrf.m (Plasma Dispersion Relation - Fluid version), general dispersion
% relation solver for multi-fluid plasma, ver 1.0. Full-matrix method.
Abstract
A general dispersion-relation solver that numerically evaluates the full propagation properties of all the waves in fluid plasmas is presented. The effects of anisotropic pressure, external magnetic fields and beams, relativistic dynamics, as well as local plasma inhomogeneity are included.
Download: http://cpc.cs.qub.ac.uk/summaries/AERF_v1_0.html
Also: http://arxiv.org/abs/1305.1205, http://arxiv.org/abs/1304.5885
Hua-sheng XIE (huashengxie@gmail.com, IFTS-ZJU)
14:00 2015/3/5
- pdrk
Cite:
[Xie2016] H. S. Xie & Y. Xiao, PDRK: A General Kinetic Dispersion Relation Solver for Magnetized Plasma, Plasma Science and Technology, 18, 97 (2016). (also: http://arxiv.org/abs/1410.2678)
% Hua-sheng XIE, huashengxie@gmail.com, IFTS-ZJU, 2014-06-01 17:00
% pdrk_em3d.m, Plasma Dispersion Relastion solver, kinetic, EM3D,
% bi-Maxwellian equilibrium distribution with parallel drift.
% Transform to matrix eigenvalue problem lambda*X=M*X.
% J-pole approximation for Z(zeta)=sum(b_j/(zeta-c_j))
%
% sparse, eigs
Abstract
A general, fast, and effective approach is developed for numerical calculation of kinetic plasma dispersion relations. The plasma dispersion function is approximated by J-pole expansion. Subsequently, the dispersion relation is transformed to a standard matrix eigenvalue problem of an equivalent linear system. The result is accurate for J = 8 except the solutions that are the little interesting heavily damped modes. In contrast to conventional approaches, such as Newton¡¯s iterative method, this approach can give either all the solutions in the system or a few solutions around the initial guess. It is also free from convergent problems. The approach is demonstrated from electrostatic one-dimensional and three-dimensional dispersion relations, to electromagnetic kinetic magnetized plasma dispersion relation for bi-Maxwellian distribution with parallel velocity drift.
Download: https://github.com/hsxie/pdrk
Also: http://arxiv.org/abs/1410.2678, http://arxiv.org/abs/1304.5885
==============================================
From: Hua-sheng XIE [huashengxie@gmail.com]
Sent: Sunday, October 12, 2014 7:37 PM
Subject: FYI: PDRK - General Kinetic Dispersion Relation Solver
Dear, Boys/Girls/Sir/Madam,
I send this to you in case you are interested in it (if not, just ignore it). The manuscript and source code can be found on arXiv (http://arxiv.org/abs/1410.2678). See 'Other formats ' for the source code: [ Download source ] , and rename it to '*.tar' and then decompress it.
This is a successor of the general multi-fluid solver PDRF (arXiv:1305.1205 or Computer Physics Communications, 2014, 185, 670 - 675). The basic idea of them was mentioned in arXiv:1304.5885 .
PDRK includes all the basic kinetic effects (Landau damping, FLR and so on) and can be accurate for the benchmark of PIC or Vlasov solvers. So, I think it would be more useful than PDRF.
Further comments and suggestions are welcome. Thanks.
Best,
Hua-sheng XIE
==============================================
Hua-sheng XIE (huashengxie@gmail.com, IFTS-ZJU)
14:06 2015/3/5
- pplu
PPLU (Plasma Physics Learning Utility) is developed to help beginners to learn and understand the forth state of matter (plasma) more easily and interested with several modules. However, I may not update it anymore.
- shafeq
'Shafeq' is written for GTC to include shifted circular tokamak equilibrium. However, you can modify it for other useages following the document 'shafeq.pdf'.
I haven't tested the part 'Reading data from TRANSP'. If you need this capability, you should check it yourself.
Hua-sheng XIE (IFTS-ZJU, huashengxie@gmail.com)
21:54 2015/3/23
List all codes
Old download page