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

fennecs/matlab/@fennecshdf5/displaypot.m
72a909f1f81bmaster

Last Change
Raw File

displaypot.m
View Options

function f=displaypot(obj,timestep)
% Display the 2D potential at time obj.t2d(timestep)
% if rcoord is true, the potential is evaluated at grid points in r,z coordinates
% if false, the potential is evaluated at grid points in magnetic field line coordinates
% clims are the values limits in eV for the color coding of the
% potential well
if iscell(timestep)
timestep=cell2mat(timestep);
end
f=figure('Name',sprintf('%s phi',obj.name));
ax1=gca;
pot=zeros(obj.pot.nr,obj.pot.nz,length(timestep));
if(any(timestep==0))
pot(:,:,timestep==0)=obj.potxt(:,:,1);
end
pot(:,:,timestep~=0)=obj.pot(:,:,timestep(timestep~=0));
N0=obj.N(:,:,1);
id=find(timestep==0);
timestep(id)=[];
Nend=obj.N(:,:,timestep);
if(~isempty(id))
N0=zeros(obj.N.nr,obj.N.nz);
Nend=cat(3,Nend(:,:,1:id-1),N0,Nend(:,:,id:end));
end
Nend=mean(Nend,3);
Nend=(Nend-0.01*max(Nend(:)))/max(Nend(:));
geomw=obj.geomweight(:,:,1);
%z(isnan(Pot))=[];
%r(isnan(Pot))=[];
%Pot(isnan(Pot))=[];
pot=mean(pot,3);
pot(obj.geomweight(:,:,1)<=0)=NaN;
contourf(obj.zgrid*1e3,obj.rgrid*1e3,pot,50,'edgecolor','none','Displayname','\phi')
xlabel('z [mm]')
ylabel('r [mm]')
xlim([obj.zgrid(1) obj.zgrid(end)]*1e3)
ylim([obj.rgrid(1) obj.rgrid(end)]*1e3)
hold(gca, 'on')
if(any(timestep~=0))
contour(obj.zgrid*1e3,obj.rgrid*1e3,Nend,[0 0],'m--','Displayname','n_e')
end
rdisp=obj.rgrid;
%% Magnetic field lines
% Blines=obj.rAthet;
% levels=linspace(min(Blines(obj.geomweight(:,:,1)>0)),max(Blines(obj.geomweight(:,:,1)>0)),10);
% 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');
geomw=obj.geomweight(:,:,1);
if (~isempty(id) && isempty(timestep)) % the
title(sprintf('Vacuum potential'))
else
title(sprintf('Potential t=%1.2f [ns]',mean(obj.t2d(timestep))*1e9))
end
c=colorbar;
colormap('jet');
% Grey outline showing the metalic walls
[c1,hContour]=contourf(ax1,obj.zgrid*1000,rdisp*1000,-geomw, [0 0]);
drawnow;
xlim(ax1,[obj.zgrid(1)*1000 obj.zgrid(end)*1000])
if(obj.conformgeom)
%ylim([ax1 ],[obj.rgrid(1)*1000 obj.rgrid(rgridend)*1000])
else
% ylim([ax1],[obj.rgrid(1)*1000 obj.rgrid(end)*1000])
end
%ylim(ax1,[0.05*1000 obj.rgrid(end)*1000])
%xlim([obj.zgrid(1) 0.185]*1e3)
xlabel(ax1,'z [mm]')
ylabel(ax1,'r [mm]')
view(ax1,2)
c.Label.String= '\phi [V]';
f.PaperUnits='centimeters';
grid(ax1, '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
% Plot of the geometry and of the different boundaries
% considered for the calculation of the currents
% geomw=obj.geomweight(:,:,1);
% geomw(geomw>1e-6)=NaN;
% geomw(geomw<1e-6)=max(Pot(:)/1e3);
% [c1,hContour]=contourf(ax1,obj.zgrid*1000,obj.rgrid*1000,geomw, max(Pot(:)/1e3)*[1 1]);
% hold on
dirw=obj.dirichletweight(:,:,1);
mask=obj.geomweight(:,:,1)<=1e-6&dirw>=1e-6;
dirw(obj.geomweight(:,:,1)<=1e-6)=max(pot(:)/1e3);
dirw(mask)=min(pot(:)/1e3);
dirw(obj.geomweight(:,:,1)>=1e-6)=NaN;
[c1,hContour2]=contourf(ax1,obj.zgrid*1000,obj.rgrid*1000,dirw,'edgecolor','none');
grid on;
drawnow;
%legend([h1,h2],{'Magnetic field lines','Equipotentials [kV]'},'location','northeast','Orientation','horizontal')
%colormap(ax,([255,165,0;255,255,255;150,150,150])./255)
% We add the coloring for the solid parts of the geometry
% Dielectrics and metallic parts
hFills=hContour2.FacePrims;
[hFills.ColorType] = deal('truecoloralpha'); % default = 'truecolor'
try
hFills(1).ColorData = uint8([255;165;0;255]);
hFills(2).ColorData = uint8([150;150;150;255]);
for idx = 3 : numel(hFills)
hFills(idx).ColorData(4) = 0; % default=255
end
catch
end
drawnow;
axis equal
pause(0.1)
% add central and external metallic walls if we have a coaxial
% configuration
if( obj.walltype >=2 && obj.walltype<=4)
rectangle('Position',[obj.zgrid(1) max(obj.r_a,obj.r_b) obj.zgrid(end)-obj.zgrid(1) 0.001]*1e3,'FaceColor',[150 150 150]/255,'Edgecolor','none')
ylimits=ylim;
ylim([ylimits(1),ylimits(2)+1])
end
if sum(obj.geomweight(:,1,1))==0
rectangle('Position',[obj.zgrid(1) obj.rgrid(1)-0.001 obj.zgrid(end)-obj.zgrid(1) 0.001]*1e3,'FaceColor',[150 150 150]/255,'Edgecolor','none')
ylimits=ylim;
ylim([ylimits(1)-1,ylimits(2)])
end
axis equal
%xlim([-100 200])
papsize=[14 8];
obj.savegraph(f,sprintf('%s/%s_phi%i',obj.folder,obj.name,floor(mean(timestep))),papsize);
end
c4science · Help