Skip to content

Solving issues 224, 222 and 221 and several things for the LinearMHDDriftkineticCC model

Byung Kyu Na requested to merge openMPforhybridmodel into devel

Following issues are solved

  1. Solves #224 (closed)

CurrentCoupling5DGradBxB_dg is removed and related kernels are commented.

  1. Solves #222 (closed)

New properties are added at Derham class

grad_bcfree, curl_bcfree and div_bcfree which are Psydac discrete derivatives without boundary operators.

  1. Solves #221 (closed)

Fix a bug in ProjectorPreconditioner in line 471:

Others

  • New callable functions of mhd_equilibrium gradB and curl_unit_b are added for the class CartesianMHDequilibrium and AxisymmMHDequilibrium.

Currently not implemented for GVECequilibrium yet.

  • New velocity sampling for Particles5D

For the parallel velocity, 1 velocity is drawn according to CDF of 1D Maxwellian (same with Particles6D)

For the perpendicular velocity, perpendicular speed v_\perp = \sqrt{v_1^2 + v_2^2} is drawn from the CDF of coordinate transformed (v_1, v_2) \rightarrow (v_\perp, \theta) 2D Maxwellian.

  • Maxwellian5DITPA class is implemented as a kinetic_background of Particles5D

  • CurrentCoupling5DDensity propagator is implemented.

  • openMP decorators are added for the model LinearMHDDriftkineticCC

Edited by Stefan Possanner

Merge request reports