normcoef_dke_yp

   Calculate normalized profile parameters for the 3D electron relativistic drift kinetic solver.

   by Y. Peysson (CEA-DRFC) (yves.peysson@cea.fr) and J. Decker (CEA-DRFC) (joan.decker@cea.fr)

0001 function [xTe_norm,xne_norm,xzni_norm,xzTi_norm,xnloss_norm,xbetath,xlnc_e,xnhu,xrnhuth] = normcoef_dke_yp(mksa,equilDKE,gridDKE)
0002 %
0003 %   Calculate normalized profile parameters for the 3D electron relativistic drift kinetic solver.
0004 %
0005 %   by Y. Peysson (CEA-DRFC) (yves.peysson@cea.fr) and J. Decker (CEA-DRFC) (joan.decker@cea.fr)
0006 %
0007 [qe,me,mp,mn,e0,mu0,re,mc2,clum,alpha] = pc_dke_yp;%Universal physics constants
0008 %
0009 xTe_norm = equilDKE.xTe/mksa.Te_ref;
0010 xzTi_norm = equilDKE.xzTi/mksa.Te_ref;
0011 xne_norm = equilDKE.xne/mksa.ne_ref;
0012 xzni_norm = equilDKE.xzni/mksa.ne_ref;
0013 %xnloss_norm = xnloss/mksa.ne_ref;
0014 xnloss_norm = [];
0015 %
0016 xTen = equilDKE.xTe/mc2;
0017 xzTin = equilDKE.xzTi/mc2;
0018 %
0019 xbetath = sqrt(xTen);%Normalized pth/mc as prescribed by Karney
0020 xlnc_e = 31.3 - 0.5*log(equilDKE.xne) + log(equilDKE.xTe*1000);%Coulomb logarithm (Sauter et al. Phys. Plasmas, 6 (1999) 2834)
0021 xnhu = qe^4*equilDKE.xne.*xlnc_e./(4*pi*e0^2*me^2.*(clum*xbetath).^3);%Electron collision frequency (s-1)
0022 xrnhuth = xnhu./mksa.nhu_ref;
0023 
0024