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

fennecs/matlab/analyse_espic2D.m
3456c8a6442fgeom_testmumps

Last Change
Raw File

analyse_espic2D.m
View Options

file='teststable126.h5';
close all hidden;
if (~exist('M','var'))
M.file=file;
M.type='all';
end
M=load_espic2d(file,M,'all');
tmin=2;
tmax=length(M.ekin);
figure
plot(M.t0d(tmin:tmax),M.ekin(tmin:tmax),'o-',...
M.t0d(tmin:tmax),M.epot(tmin:tmax),'d-',...
M.t0d(tmin:tmax),M.etot(tmin:tmax),'h-',...
M.t0d(tmin:tmax),M.eerr(tmin:tmax),'x--')
legend('ekin', 'epot', 'etot','eerr')
xlabel('Time [s]')
ylabel('Energies [J]')
grid on
figure
semilogy(M.t0d(tmin:tmax),abs(M.eerr(tmin:tmax)/M.etot(2)*100),'h-')
xlabel('t [s]')
ylabel('E_{err} %')
grid on
figure
[Z,R] = meshgrid (M.zgrid,M.rgrid);
subplot(2,1,1)
contourf(Z,R,M.Br')
xlabel('Z')
ylabel('R')
c = colorbar;
c.Label.String= 'B_r [T]';
grid on
subplot(2,1,2)
contourf(Z,R,M.Bz')
xlabel('Z [m]')
ylabel('R [m]')
c = colorbar;
c.Label.String= 'B_z [T]';
grid on
figure
contourf(Z,R,M.B',100,'LineColor','none');
shading flat
colormap('jet')
hold on
quiver(Z,R,M.Bz',M.Br','k')
% h=streamslice(Z,R,M.Bz',M.Br');
% set(h,'Color','w')
% figure
% ax3=subplot(2,1,1);
% contourf(ax3,M.zgrid,M.rgrid,mean(M.Er(:,:,end-2000:end),3)')
% xlabel(ax3,'Z [m]')
% ylabel(ax3,'R [m]')
% c = colorbar(ax3);
% c.Label.String= 'E_r [V/m]';
% title(ax3,'Average Radial Electric field')
% grid on
% ax3.Children(1).ZDataSource='data.Er';
%
% ax4=subplot(2,1,2);
% contourf(ax4,M.zgrid,M.rgrid,mean(M.Ez(:,:,end-2000:end),3)')
% xlabel(ax4,'Z [m]')
% ylabel(ax4,'R [m]')
% c = colorbar(ax4);
% c.Label.String= 'E_z [V/m]';
% title(ax4,'Average Axial Electric field')
% grid on
% ax4.Children(1).ZDataSource='data.Ez';
figure
zpos=33;
tinit=1;
tend=size(M.N,3);
deltat=floor(tend/2);
H0=3.2e-14;
P0=8.78e-25;
Mirrorratio=(M.Rcurv-1)/(M.Rcurv+1);
pot=M.pot(zpos,:,tend);
psi=1+M.qe*pot(:)/H0-1/(2*M.me*H0)*(P0./M.rgrid+M.qe*0.5*M.B0.*(M.rgrid-M.width/pi*Mirrorratio*cos(2*pi*M.zgrid(zpos)/M.width)*besseli(1,2*pi*M.rgrid/M.width))).^2;
[~,I]=max(M.N(:,zpos,tinit));
plot(M.rgrid,M.N(:,zpos,tinit),'bx-','DisplayName',sprintf('t=%1.2f[ns]',M.t2d(tinit)*1e9))
hold on
plot(M.rgrid,mean(M.N(:,zpos,tend-deltat:tend),3),'rx-','DisplayName',sprintf('t=[%1.2f-%1.2f] [ns] averaged',M.t2d(tend-deltat)*1e9,M.t2d(tend)*1e9))
%plot(M.rgrid(I-2:end),1./M.rgrid(I-2:end)*M.N(I,zpos,tinit)*M.rgrid(I),'DisplayName','N=c*1/r')
xlabel('R [m]')
ylabel('N')
J=find(psi>0);
plot(M.rgrid(J),psi(J))
title(sprintf('Number of particles at z=%1.2e[m]',M.zgrid(zpos)))
legend
figure
N=M.N(I+2,zpos,:);
plot(M.t2d,N(:))
xlabel('t [s]')
ylabel('N')
title(sprintf('Evolution of number of particles in cell at z=%1.2e[m] r=%1.2e[m]',M.zgrid(zpos),M.rgrid(I+2)))
fieldstep=1;
data.pot=M.pot(:,:,fieldstep)';
data.pos=[M.Z(:,fieldstep),M.R(:,fieldstep)];
data.Ez=M.Ez(:,:,fieldstep)';
data.Er=M.Er(:,:,fieldstep)';
Rlarmor=(M.me*sqrt(M.VR.^2+M.VTHET.^2)/(M.qe*M.B0));
Rmax=M.R;
Rmin=M.R-Rlarmor;
[Isoutsideuprow,Isoutsideupcol]=find(Rmax-M.rgrid(end)>0);
[Isoutsidedownrow,Isoutsidedowncol]=find(M.rgrid(1)-Rmin>0);
f=uifigure('Name','Grid data');
mf=uipanel(f,'Position',[5 50 f.Position(3)-10 f.Position(4)-55]);
mf.AutoResizeChildren='off';
m=uipanel(f,'Position',[5 5 f.Position(3)-10 40]);
sgtitle(mf,sprintf('step=%d t=%0.5e s',fieldstep*M.it1,M.t2d(fieldstep)))
data.pos=[M.Z(:,fieldstep),M.R(:,fieldstep)];
data.pot=M.pot(:,:,fieldstep)';
data.Ez=M.Ez(:,:,fieldstep)';
data.Er=M.Er(:,:,fieldstep)';
sld = uislider(m,'Position',[10 30 0.6*m.Position(3) 3]);
sld.Value=fieldstep;
sld.Limits=[1 size(M.pot,3)];
edt = uieditfield(m,'numeric','Limits',[1 size(M.pot,3)],'Value',1);
edt.Position=[sld.Position(1)+sld.Position(3)+25 5 40 20];
edt.RoundFractionalValues='on';
Printbt=uibutton(m,'Position',[edt.Position(1)+edt.Position(3)+10 5 40 20],'Text', 'Save');
%Playbt=uibutton(m,'Position',[Printbt.Position(1)+Printbt.Position(3)+10 5 40 20],'Text', 'Play/Pause');
sld.ValueChangingFcn={@updatefigdata,edt,mf,M};
%sld.ValueChangedFcn={@updatefigdata,edt,mf,M,data};
edt.ValueChangedFcn={@updatefigdata,sld,mf,M};
Printbt.ButtonPushedFcn={@plotGridButtonPushed,M,sld};
%Playbt.ButtonPushedFcn={@playButtonPushed,M,data,sld,edt};
PlotEspic2dgriddata(mf,M,fieldstep);
BrillouinRatio=2*M.omepe^2/M.omece^2
c4science · Help