HeLaZ (Hermite-Laguerre Z-pinch solver, 2020)

How to run it :

0. Required librairies : hdf5, fftw3-mpi, mpi Fortran compiler (Intel?)

1. Be sure to have correct paths in local/dirs.inc for the different libraries
2. You can compile from HeLaZ/ using make and launch from HeLaZ/wk using ./../bin/HeLaZ
3. To have a better interface, open a script HeLaZ/wk/parameters*.m and run it to set up a wanted simulation.
4. You can obtain various plots and gifs using HeLaZ/wk/analysis_2D.m once the simulation is done. To select the correct output file, run parameters*.m with the corresponding simulation parameters and then run analysis_2D.m (everything with matlab from wk/)

Logbook

(Current versions : 2.3.1.1)

0. Write MOLI matlab solver in Fortran using Monli1D as starting point

0.0 go from 1D space to 2D fourier and from Hermite basis to Hermite-Laguerre basis

0.1 implement linear Poisson equation in fourier space

0.2 implement moment hierarchy linear terms

0.3 RK4 time solver

0.4 Benchmark with MOLI matlab results for Z-pinch (cf. kz_linear script)

		Note : benchmark_*.m compares MOLI and HeLaZ linear results

0.5 Load COSOlver matrices

0.6 Benchmarks now include Dougherty, Lenard-Bernstein and Full Coulomb collision operators

	    Note : for full Coulomb, one must store a precomputed matrix from COSOlver in the iCa folder

1. Implementation of the non linear Poisson brackets term

1.0 FFTW3 has been used to treat the convolution as a product and discrete fourier transform

1.1 Methods in fourier_mod.f90 have been validated by tests on Hasegawa Wakatani system

1.1 Qualitative test : find similar turbulences as Hasegawa Wakatani system with few moments

1.2 Zonal flows are observed in a similar way to Ricci Rogers 2006 with GS2

1.3 Linear analysis showed that a certain amount of PJ are recquired to trigger mode

		1.3.1 The \eta_B = 0.5 case is easier since it converged better in linear analysis than \eta_B = 1.0

		1.3.2 Collisionality helps

1.4 Quantitative study with stationary average particle flux \Gamma_\infty

		1.4.1 Convergence study of \Gamma_\infty w.r.t. P and J

		1.4.2 Direct comparison with GS2 results of Ricci,Rogers 2006

		1.4.3 Code too expensive in sequential to reach PJ convergence

2. MPI parallel version (branch MPI)

2.1 First compilable parallel version (1D parallel along kr)

		2.1.1 Benchmarks, profiling and portability of the code, max speed up is arround 10 for 12 proc. when using large enough resolution or polynomials degree. Should be increased with other dimension distribution.

2.2 Allow restart with different P,J values

2.3 GK Dougherty operator

		2.3.1 Convergence study of steady radial particle transport

			2.3.1.1 Detailed analysis of numerical stability for P,J = 10,5.

2.4 Data distribution along P (under consideration)

2.5 GK Coulomb operator

3. GK 3D version, kr,kz,kpar for linear device

3. DK 3D version, kr,kz,kpar for linear device

4. DK+GK 3D version, kr,kz,kpar for linear device

4. 3D version with curvature