Homec4science
Diffusion FENNECS lite (master) Browse

fennecslite/magnet/Cryogenic/Draw_B_Ellip_170_geom_refurb.m
6b3731566193master

Last Change
Raw File

Draw_B_Ellip_170_geom_refurb.m
View Options

% creates the magnetic field h5 file necessary for fennecs
% This uses the geometry of the refurbished 170GHz coaxial gyrotron gun
% The magnetic field is the one created by the asg magnet
% Both the magnet and geometry are definde in the Final report of the
% Development of the european gyrotron CCCGDS6
%% Import and calculate the magnetic field
magnet=load('asg_magnet_red.mat');
%% calculate magnetic field and magnetic vector
idr=1:10:length(magnet.r);
idz=1:10:length(magnet.z);
% Define the coils currents
%I = [ 9.7 7.7 96.2499 91.0124 87.4723];
% I = [0 0 4.5 1.967 88.2799];
% name='phiBprofile_refurbasg_4_5_red';
% I = [0 0 4.8 1.967 88.2799];
% name='phiBprofile_refurbasg_well_red';
% I = [0 0 5 1.967 88.2799];
% name='phiBprofile_refurbasg_5_red';
% I = [0 0 5.15 1.967 88.2799];
% name='phiBprofile_refurbasg_limdwn_red';
%I = [0 0 5.35 1.967 88.2799];
%name='phiBprofile_refurbasg_limup_red';
% I = [0 0 5.5 1.967 88.2799];
% name='phiBprofile_refurbasg_5_5_red';
%I = [0 0 5.6 1.967 88.2799];
%name='phiBprofile_refurbasg_5_6_red';
%I = [0 0 5.7 1.967 88.2799];
%name='phiBprofile_refurbasg_5_7_red';
% I = [0 0 5.75 1.967 88.2799];
% name='phiBprofile_refurbasg_5_75_red';
% I = [0 0 5.8 1.967 88.2799];
% name='phiBprofile_refurbasg_5_8_red';
% I = [0 0 5.85 1.967 88.2799];
% name='phiBprofile_refurbasg_5_85_red';
% I = [0 0 5.9 1.967 88.2799];
% name='phiBprofile_refurbasg_5_9_red';
% I = [0 0 6 1.967 88.2799];
% name='phiBprofile_refurbasg_6_red';
% I = [0 0 6.4 1.967 88.2799];
% name='phiBprofile_refurbasg_6_4_red';
% I = [0 0 6.79541 1.967 88.2799];
% name='phiBprofile_refurbasg_nominal_red';
% for Ib=linspace(4.6,6.8,12)
%Ib=4.4;
%I = [0 0 Ib 1.967 88.2799];
%name=sprintf('phiBprofile_refurbasg_%2d_red',floor(Ib*10));
%
%I = 0.6*[0 0 6.79541 1.967 88.2799];
%name='phiBprofile_refurbasg_nominal_0_6';
I = 0.6*[8.037 7.030 8.154 3.485 86.706];
name='phiBprofile_refurbasg_nominal2_0_6';
Icoil=[I(1) I(2) -I(3)-I(5) -I(3)-I(5) I(4)+I(5) I(4)+I(5) I(5) I(5)];
rmin = [0.13958 0.13958 0.13710 0.16952 0.13670 0.19101 0.13458 0.17197 ];
rmax = [0.142694 0.142694 0.16733 0.20005 0.18941 0.23666 0.16972 0.20412 ];
zmin = [0.07854 0.13854 0.18200 0.18200 0.31778 0.31748 0.50685 0.50685 ];
zmax = [0.09846 0.15846 0.28350 0.28350 0.48977 0.49017 0.71128 0.71128 ];
% Nturn = [ 218 218 2528.5 2626 7689 11517.5 6110 9656.5 ];
% na = 2. *[ 10 10 10 10 10 10 10 10 ]; % Number of subcoils (radial direction
% nb = 2. *[ 10 10 10 10 10 10 10 10 ]; % Number of subcoils in longitudinal direction
% Icoil=Icoil.*Nturn./(na.*nb);
[Z,R] = meshgrid(magnet.z(idz),magnet.r(idr));
Aphi=zeros(length(idr),length(idz));
Bz=zeros(length(idr),length(idz));
Br=zeros(length(idr),length(idz));
for i=1:size(magnet.subcoils,1)
Aphi=Aphi+Icoil(i)*magnet.subcoils{i,1}(idr,idz);
Bz=Bz+Icoil(i)*magnet.subcoils{i,2}(idr,idz);
Br=Br+Icoil(i)*magnet.subcoils{i,3}(idr,idz);
end
%% Plots
f=figure;
%%
%axis equal
raphi=R.*Aphi;
lvls=logspace(-4,log10(max(raphi(:))),50);
[~,cont1]=contour(Z*1e3,R*1e3,raphi,lvls,'b:','linewidth',1.5);
%xlim([min(magnet.z) max(magnet.z)]*1e3)
%ylim([min(magnet.r) max(magnet.r)]*1e3)
%[~,cont2]=contour(Zphi,Rphi,Phi,20,'b');
rectangle('Position',[-0.1, 0.045, 0.292, 0.092-0.045]*1e3,'EdgeColor','black','Linestyle','--','linewidth',2)
%legend([plothandle,cont1],{'Gun geometry', 'Magnetic field lines'},'location','southwest')
%f.PaperUnits='centimeters';
%f.PaperSize=[12,8];
xlabel('z [mm]')
ylabel('r [mm]')
axis equal
grid on
ylim([0,inf])
xlim([-300 750])
for i=1:length(rmin)
rectangle('Position',[zmin(i) rmin(i) zmax(i)-zmin(i) rmax(i)-rmin(i)]*1e3,'edgecolor','r')
text(zmin(i)*1e3-10,(rmin(i)+rmax(i))/2*1e3,sprintf('S%i',i),'fontsize',10,'color','r')
end
print(f,name,'-dpdf','-fillpage')
savefig(f,name)
set(f, 'Color', 'w');
export_fig(f,name,'-eps')
hold off
%% Save magnetic field and geometry to disk
save=true;
overwrite=true;
if save
idr=1:1:length(magnet.r);
idz=1:1:length(magnet.z);
Aphi=zeros(length(idr),length(idz));
Bz=zeros(length(idr),length(idz));
Br=zeros(length(idr),length(idz));
for i=1:size(magnet.subcoils,1)
Aphi=Aphi+Icoil(i)*magnet.subcoils{i,1}(idr,idz);
Bz=Bz+Icoil(i)*magnet.subcoils{i,2}(idr,idz);
Br=Br+Icoil(i)*magnet.subcoils{i,3}(idr,idz);
end
savemagtoh5([name,'.h5'],magnet.r,magnet.z,Aphi,Br,Bz,overwrite);
%savegeomtoh5('geom_ceramics.h5',geomcells,1e-2,overwrite);
end
% end
c4science · Help