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

fennecs/magnet/Cryogenic/B_Mat88_Cryogenic_170_test.m
b121cf310431master

Last Change
Raw File

B_Mat88_Cryogenic_170_test.m
View Options

% This script tests the coefficients of the matrix found by
% 'B_Mat88_Cryogenic_170
%
% JHE, July 31st, 2017
%
I= 100*[0.57712 , 0.68596 , 1.07780 , 1.07780] ; % Nominal currents used on MAy 29th, 2017 to apply first SP pulses
% r_cat_min = 0.05419; % largest emitter radius of first prototype 170GHz 1MW
% r_launcher = 0.02068; % launcher radius
%
% z_cat_min = 0.1335;
% z_launcher = 0.78848;
% Include GT170 Geometry
GT170_1MW_Geometry
%
% Taking parameters defined with on-axis fields
%gamma = [0.002933681588584, -0.006273546452915, -0.007626766134967, 0.009858220586519];
% Taking more stringent condition (on axis field at gun, off-axis field at launcher
% gamma = [ 0.00293370270241, -0.00627361773021, -0.00762699351991 , 0.00987214759867];
% Taking parameters defined with potential vector (defined as delta in
% B_Matrix_Cryogenic_170, with another factor 1000 omitted
gamma = [0.081474637274171, -0.174106309671082, -0.203575059328030, 0.260894824884893];
z = linspace(z_launcher-0.01,z_launcher+0.01,41);
r = linspace(0.015, 0.030,16);
[Z,R] = meshgrid(z,r);
[aphi_cat] = B_Ellip_Cryogenic_170('aphi','cryogenic',I,r_cat_min,z_cat_min);
raphi_cat = r_cat_min * aphi_cat;
[Aphi] = B_Ellip_Cryogenic_170('aphi','cryogenic',I,R,Z);
RAphi = R.*Aphi;
figure(1)
contour(Z,R,RAphi,[raphi_cat raphi_cat],'--','LineWidth',2)
hold on
I(4) = I(4) + 10;
[aphi_cat] = B_Ellip_Cryogenic_170('aphi','cryogenic',I,r_cat_min,z_cat_min);
raphi_cat = r_cat_min * aphi_cat;
[Aphi] = B_Ellip_Cryogenic_170('aphi','cryogenic',I,R,Z);
RAphi = R.*Aphi;
figure(1)
contour(Z,R,RAphi,[raphi_cat raphi_cat],'-.','LineWidth',2)
I4lim = - (gamma(1) * I(1) +gamma(2)*I(2) +gamma(3)*I(3)) / gamma(4)
I(4) = I4lim;
[aphi_cat] = B_Ellip_Cryogenic_170('aphi','cryogenic',I,r_cat_min,z_cat_min);
raphi_cat = r_cat_min * aphi_cat;
[Aphi] = B_Ellip_Cryogenic_170('aphi','cryogenic',I,R,Z);
RAphi = R.*Aphi;
figure(1)
contour(Z,R,RAphi,[raphi_cat raphi_cat],'LineWidth',2)
plot(z_launcher,r_launcher,'r*')
c4science · Help