Homec4science
Diffusion FENNECS lite (master) Browse

fennecslite/magnet/bi-freq/B_Ellip_bifreq.m
6b3731566193master

Last Change
Raw File

B_Ellip_bifreq.m
View Options

function [B] = B_Ellip_bifreq(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 'normal' % Geometry of the gt118 X3 gyrotron magnet used at TCV
Ncoil = 5;
Nturn = [3567 1659 4307 11398 13384];
rmin = [0.1457 0.1584 0.1452 0.1632 0.1482];
rmax = [0.1582 0.1615 0.1632 0.1946 0.1837];
zcen = [0.2598 0.2598 0.4184 0.4184 0.6022];
deltaz = [0.0540 0.0540 0.1057 0.1057 0.1094];
zmin = zcen-deltaz;
zmax = zcen+deltaz;
na = 5. *[10 10 10 10 10]; % Number of subcoils (radial direction
nb = 5. *[10 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(1) I(1) I(1)]; % 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