Homec4science
Diffusion FENNECS: Finite Element Non-Neutral Electron Cloud Simulator (stable) Browse

fennecs/matlab/@fennecshdf5/displayEffTrap.m
0dc1ab28c2d5stable

Last Change
Raw File

displayEffTrap.m
View Options

function displayEffTrap(obj,time2d,clims,mag)
% Plot the larmor radius for created particles with low energy
% the larmor radius is calculated by considering that the
% initial perpendicular velocity \approx the ExB velocity
% potential well is shown with isolines on top
if iscell(time2d)
time2d=cell2mat(time2d);
end
if nargin<3
clims=inf*[-1 1];
end
if nargin<4
mag=[];
end
[zgrid,rgrid]=ndgrid(obj.zgrid,obj.rgrid);
geomw=obj.geomweight(:,:,1);
inside=griddedInterpolant(zgrid,rgrid,obj.geomweight(:,:,1)');
% Compute the Larmor radius assuming v_perp=v_exb
if time2d>0
Er=obj.Er(:,:,time2d);
Ez=obj.Ez(:,:,time2d);
else
Er=obj.Erxt(:,:,1);
Ez=obj.Ezxt(:,:,1);
end
rl=abs(obj.me/obj.qe*(-Er.*obj.Bz'+Ez.*obj.Br')./(obj.B.^3)');
% Remove the points outside the geometry
rl(obj.geomweight(:,:,1)<0)=0;
% calculate the normal vector to the magnetic field lines at each position
perpB=zeros(size(obj.sinthet,1),size(obj.sinthet,2),2);
perpB(:,:,1)=sign(-Er).*obj.sinthet;
perpB(:,:,2)=-sign(-Er).*obj.costhet;
pos=zeros(size(obj.sinthet,1),size(obj.sinthet,2),2);
pos(:,:,1)=rgrid';
pos(:,:,2)=zgrid';
posu=pos+2*rl.*perpB;
isinside=inside(posu(:,:,2),posu(:,:,1))>0;
model=obj.potentialwellmodel(time2d,true,mag);
z=model.z;
r=model.r;
Pot=model.pot;
Potdepth=scatteredInterpolant(z',r',Pot,'natural');
effectdepth=0.5*(Potdepth(posu(:,:,2),posu(:,:,1))+Potdepth(pos(:,:,2),pos(:,:,1)));
effectdepth(geomw<0)=NaN;
effectdepth(~isinside)=NaN;
effectdepth(effectdepth<0)=NaN;
f=figure;
contourf(obj.zgrid*1000,obj.rgrid*1000,effectdepth,40,'edgecolor','none')
hold on
contour(obj.zgrid*1000,obj.rgrid*1000,geomw,[0 0],'r-')
%% Magnetic field lines
Blines=obj.rAthet;
levels=linspace(min(Blines(obj.geomweight(:,:,1)>0)),max(Blines(obj.geomweight(:,:,1)>0)),20);
Blines(obj.geomweight(:,:,1)<0)=NaN;
[~,h1]=contour(obj.zgrid*1000,obj.rgrid*1000,Blines,real(levels),'k-.','linewidth',1.2,'Displayname','Magnetic field lines');
if time2d>0
N=obj.N(:,:,time2d);
N=N/max(N(:));
contour(obj.zgrid*1000,obj.rgrid*1000,N,'k--')
end
colormap(jet)
xlabel('z [mm]')
ylabel('r [mm]')
c=colorbar;
c.Label.String='Effective well [eV]';
caxis(clims)
% Grey outline showing the metalic walls
geomw(obj.geomweight(:,:,1)>0)=-1;
geomw(obj.geomweight(:,:,1)<=0)=1;
[c1,hContour]=contourf(obj.zgrid*1000,obj.rgrid*1000,geomw, [0 0]);
drawnow;
grid on;
hFills=hContour.FacePrims;
[hFills.ColorType] = deal('truecoloralpha'); % default = 'truecolor'
try
drawnow
hFills(1).ColorData = uint8([150;150;150;255]);
for idx = 2 : numel(hFills)
hFills(idx).ColorData(4) = 0; % default=255
end
catch
end
papsize=[14 8];
obj.savegraph(f,sprintf('%s/%s_effwellr_rl_%i',obj.folder,obj.name,floor(mean(time2d))),papsize);
end
c4science · Help