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

fennecs/matlab/@fennecshdf5/displayDensPsi.m
2eb9a934d936master

Last Change
Raw File

displayDensPsi.m
View Options

function displayDensPsi(obj,tsteps)
%DISPLAYDENSPSI plot the fluid density averaged in time for times t2d(tsteps)
% Displays also the geometry of the metallic boundaries and the
% magnetic field lines
% Show on top the contour of the Psi function for the Davidson
% loading
f=figure('Name', sprintf('%sfluid_dens_psi',obj.name));
ax1=gca;
if nargin<2
tsteps=length(obj.t2d);
end
geomw=obj.geomweight(:,:,1);
dens=mean(obj.N(:,:,tsteps),3);
dens(geomw<=0)=0;
maxdens=max(dens(:));
[h,curve]=contourf(ax1,obj.zgrid*1000,obj.rgrid*1000,dens,30);
hold on
set(curve,'linecolor','none');
%contour(ax1,obj.zgrid,obj.rgrid,obj.geomweight(:,:,1),[0 0],'r-','linewidth',1.5,'Displayname','Boundaries');
if(obj.conformgeom)
ylim(ax1,[obj.rgrid(1) obj.rgrid(sum(obj.nnr(1:2)))]*1000)
else
ylim(ax1,[obj.rgrid(1) obj.rgrid(end)]*1000)
end
xlim(ax1,[obj.zgrid(1) obj.zgrid(end)]*1000)
xlabel(ax1,'z [mm]')
ylabel(ax1,'r [mm]')
%title(ax1,sprintf('Density t=[%1.2g-%1.2g]s n_e=%1.2gm^{-3}',obj.t2d(fieldstart),obj.t2d(fieldend),double(maxdens)))
c = colorbar(ax1);
c.Label.String= 'electron density [m^{-3}]';
view(ax1,2)
colormap(flipud(hot));
% Draws the 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(ax1,obj.zgrid*1000,obj.rgrid*1000,Blines,real(levels),'-.','color','k','linewidth',1.5,'Displayname','Magnetic field lines');
% Draw the metallic boundaries and the geometry itself
[c1,hContour]=contourf(ax1,obj.zgrid*1000,obj.rgrid*1000,-geomw,[0,0],'linewidth',1.5);
drawnow;
xlim(ax1,[obj.zgrid(1)*1000 obj.zgrid(end)*1000])
% Change the color of the metallic boundaries to grey
hFills=hContour.FacePrims;
[hFills.ColorType] = deal('truecoloralpha'); % default = 'truecolor'
try
hFills(end).ColorData = uint8([150;150;150;255]);
for idx = 1 : numel(hFills)-1
hFills(idx).ColorData(4) = 0; % default=255
end
catch
end
if(obj.species(1).R.nt<2)
h0=obj.H0;
p0=obj.P0;
else
h0=mean(H(obj,{1:obj.species(1).VR.nparts,obj.species(1).VR.nt,false}));
p0=mean(P(obj,{1:obj.species(1).VR.nparts,obj.species(1).VR.nt,false}));
end
Mirrorratio=(obj.Rcurv-1)/(obj.Rcurv+1);
locpot=mean(obj.pot(:,:,tsteps),3);
psi=1+obj.qe*locpot/h0-1/(2*obj.me*h0)*(p0./obj.rgrid+obj.qe*0.5*obj.B0.*(obj.rgrid-obj.width/pi*Mirrorratio*cos(2*pi*obj.zgrid'/obj.width).*besseli(1,2*pi*obj.rgrid/obj.width))).^2;
contour(obj.zgrid*1000,obj.rgrid*1000,psi,[0 0],'k--')
grid on
f.PaperOrientation='landscape';
f.PaperUnits='centimeters';
papsize=[16 14];
f.PaperSize=papsize;
set(ax1,'fontsize',14)
%axis equal
obj.savegraph(f,sprintf('%sfluid_dens_psi',obj.name))
end
c4science · Help