Homec4science
Diffusion FENNECS lite (master) Browse

fennecslite/magnet/gt118/B_Ellip_gt118.m
6b3731566193master

Last Change
Raw File

B_Ellip_gt118.m
View Options

function [B] = B_Ellip_gt118(mode,magnet,I,r,z)
%
% Call: [B] = B_Ellip_gt118(mode,magnet,r,z)
% Examples:
% 1. Plot axis magnetic field profile:
%
% >> z = linspace(0,1,101);
% >> r = 0*z + 0.0;
% >> I = [9.85 -15.80 33.59 33.21];
% >> [B] = B_Ellip_gt118('bz','cryogenic',I,r,z);
% >> figure(1); plot(z,B)
%
% 2. Plot field lines (r* aphi = cst)
%
% >> z = linspace(0,1,51);
% >> r = linspace(0,0.2,21);
% >> [Z,R] = meshgrid(z,r);
% >> I = [9.85 -15.80 33.59 33.21];
% >> [Aphi] = B_Ellip_gt118('aphi','cryogenic',I,R,Z);
% >> figure(2); contour(Z,R,R.*Aphi);
%
% Magnetic field associated with GT118 GHz magnet, with currents Igun and Icav,
% as read on power supplies
% The computation is based on the elliptic integral formulae that can be found in:
% Smythe: "Static and Dynamic Electricity", Third edition, McGraw-Hill, p.290.
%
% Version : 1.0 G. Le Bars Mar. 14, 2023
%
%
% mode : 'bz','br','dbzbz','dbzdr','dbrdz','dbrdr', 'd2bzdz2' or 'aphi'
% magnet : 'X3' only so far
% z : a row vector (in [m]) or a 2D matrix containing the z locations where the field is to be evaluated
% r : a row vector (in [m]) or a 2D matrix containing the r locations where the field is to be evaluated
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 10T magnet for DNP gyrotron coils geometry [m]
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
switch lower(magnet)
case 'x3' % Geometry of the gt118 X3 gyrotron magnet used at TCV
zc = 0.435; % center of the cavity and center of the main coils
Ncoil = 4;
Nturn = [3346 7291 21913 21702];
rmin = [0.14 0.14 0.098 0.098];
rmax = [0.1579 0.1579 0.1592 0.1592];
zmin = [-0.327 -0.282 -0.1216 0.0236]+zc;
zmax = [-0.307 -0.238 -0.0236 0.1216]+zc;
na = 8. *[10 10 10 10]; % Number of subcoils (radial direction
nb = 8. *[10 10 10 10]; % Number of subcoils in longitudinal direction
% default current
% I = [9.85 -15.80 33.59 33.21]
Current = [ I(1) I(2) I(3) I(4)]; % Current flowing in each power supply
rcen = (rmin + rmax)/2;
zcen = (zmin + zmax)/2;
dr = rmax - rmin;
dz = zmax - zmin;
coilsurf= dr .* dz;
JTot = Current .* Nturn; % Total current flowing in each coil
J = JTot./(dr.*dz); % Current density in each coil
Ic1 = JTot ; % Total current flowing in each coil [A]
otherwise
error('Unknown magnet: "%s"', magnet)
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Construct array of current loops
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Number of coils
ncoil = length(rcen);
% Number of subcoils
ntot = sum(na.*nb);
% Preallocate space
r2 = zeros(1,ntot);
z2 = zeros(1,ntot);
Ic = zeros(ntot,1);
% isub is the index of subcoil
isub = 0;
for icoil=1:ncoil
for ib=1:nb(icoil)
for ia=1:na(icoil)
isub = isub +1;
r2(isub) = (rcen(icoil) - dr(icoil)/2) + (dr(icoil)/nb(icoil))*(ib-0.5);
z2(isub) = (zcen(icoil) - dz(icoil)/2) + (dz(icoil)/na(icoil))*(ia-0.5);
Ic(isub) = Ic1(icoil) / (na(icoil)*nb(icoil));
end
end
end
% printf('isub =',isub)
% Transform array of evaluation points in a one-dimensional array
%
r1d = reshape(r,1,prod(size(r)));
z1d = reshape(z,1,prod(size(z)));
%
% Suppress locations where Green function diverges (i.e. on the sources)
%
for i=1:length(r2),
[I] = find((r1d.*r1d -r2(i)*r2(i) == 0) & (z1d.*z1d - z2(i)^2)==0);
if ~isempty(I)
r1d(I)=NaN;
z1d(I)=NaN;
end
end
% Call greenem
b = greenem_jph(mode,r1d,z1d,r2,z2);
if(~iscell(mode))
nbmodes=length({mode});
else
nbmodes=length(mode);
end
B=zeros(size(r,1),size(r,2),nbmodes);
% b is a rectangular matrix (length(r1d) x length(r2)) which needs to be multiplied by the
% column vector Ic to get the proper contribution of each source.
for i=1:size(b,3)
temp = b(:,:,i) * Ic;
B(:,:,i) = reshape(temp,size(r,1),size(r,2));
end
B(find((B==Inf)|(B==-Inf)))= 0;
B(find(isnan(B))) = 0;
return
c4science · Help