Homec4science
Diffusion FENNECS lite (3Dmpimode) Browse

fennecslite/geometries/experiment_both/exp_both_withvessel.m
2e0aa43721aa3Dmpimode

Last Change
Raw File

exp_both_withvessel.m
View Options

%% Create the geometry 3 for the T-REX experiment
% uses the upper ellipse and center column with extruded ellipse.
% bottom axial position
zbottom=305; % mm
radialindent=8; %mm
% Outer vacuum vessel
splid=1;
geomcells(splid).Dval=0; % V
geomcells(splid).order=3;
geomcells(splid).dim=2;
geomcells(splid).name='vacuum vessel';
geomcells(splid).type=0;
geomcells(splid).periodic=0;
rmax=50; %mm
lw=linspace(34,rmax,8);
rwellipse=linspace(zbottom,480,20);
uw=flip(linspace(0,rmax,30));
geomcells(splid).points=flip([zbottom*ones(size(lw')) lw';
rwellipse' rmax*ones(size(rwellipse'));
480*ones(size(uw')) uw']/1e3);
% Central electrode
splid=2;
geomcells(splid).Dval=-20000; % V
geomcells(splid).order=3;
geomcells(splid).dim=2;
%splines(1).epsge=1e-6;
%splines(1).epsce=1e-6;
geomcells(splid).name='central electrode';
geomcells(splid).type=0;
geomcells(splid).periodic=0;
lw=linspace(10-radialindent,10,10);
zwellipse=linspace(zbottom+30.5,414.5,150);
rwellipse=10*ones(size(zwellipse));
uw=flip(linspace(0,10,30));
%tilted ellipse
r_c=10;
z_c=250+125;
cosa=1;
sina=0;
Lz=25;
Lr=2.5;
deltar=(rwellipse-r_c);
deltaz=(zwellipse-z_c);
D=((deltaz)/Lz).^2+((deltar)/Lr).^2;
w=1-D;
ellipsez=zwellipse(w>=0);
deltaz=ellipsez-z_c;
a=(Lr^2*sina^2+Lz^2*cosa^2);
b=2*deltaz*sina*cosa*(Lr^2-Lz^2);
c=-Lr^2*Lz^2 +deltaz.^2*(Lr^2*cosa^2+Lz^2*sina^2);
rwellipse(w>=0)=r_c+(-b-sqrt(b.^2-4*a.*c))/2/a;
%ellipse=
lowerring=linspace(zbottom,zbottom+29.5,40);
geomcells(splid).points=[ lowerring' (10-radialindent)*ones(size(lowerring'));
(zbottom+30)*ones(size(lw')) lw';
zwellipse' rwellipse';
(zwellipse(end)+0.5)*ones(size(uw')) uw']/1e3;
% Bottom plate
splid=3;
geomcells(splid).Dval=0; % V
geomcells(splid).order=3;
geomcells(splid).dim=2;
%splines(1).epsge=1e-6;
%splines(1).epsce=1e-6;
geomcells(splid).name='vacuum vessel';
geomcells(splid).type=0;
geomcells(splid).periodic=0;
lw=linspace(10-radialindent,24+radialindent,50);
geomcells(splid).points=flip([
zbottom*ones(size(lw')) lw';]/1e3);
% Outer electrode
splid=4;
geomcells(splid).Dval=0; % V
geomcells(splid).order=3;
geomcells(splid).dim=2;
%splines(1).epsge=1e-6;
%splines(1).epsce=1e-6;
geomcells(splid).name='outer electrode';
geomcells(splid).type=0;
geomcells(splid).periodic=0;
lbaser=linspace(24+radialindent,24,20);
base=linspace(zbottom+30.5,349,80);
lowerbase=linspace(zbottom,zbottom+29.5,50);
%tiltedwall
alpha=0.1745;
tiltedz=linspace(350,180+250,500);
tiltedr=(tiltedz-tiltedz(1))*tan(alpha)+24;
%tilted ellipse
cosa=cos(alpha);
sina=sin(alpha);
r_c=28;
z_c=375;
Lz=25;
Lr=5;
deltar=(tiltedr-r_c);
deltaz=(tiltedz-z_c);
D=((deltaz*cosa+deltar*sina)/Lz).^2+((deltaz*sina-deltar*cosa)/Lr).^2;
w=1-D;
ellipsez=tiltedz(w>=0);
deltaz=ellipsez-z_c;
a=(Lr^2*sina^2+Lz^2*cosa^2);
b=2*deltaz*sina*cosa*(Lr^2-Lz^2);
c=-Lr^2*Lz^2 +deltaz.^2*(Lr^2*cosa^2+Lz^2*sina^2);
tiltedr(w>=0)=r_c+(-b-sqrt(b.^2-4*a.*c))/2/a;
axialend=linspace(tiltedr(end)+0.5,39.5,15);
upperend=flip(linspace(zbottom,180+250,60));
geomcells(splid).points=flip([lowerbase' (24+radialindent)*ones(size(lowerbase'));
zbottom+30*ones(size(lbaser')) lbaser';
base' 24*ones(size(base))';
tiltedz' tiltedr';
(180+250)*ones(size(axialend')) axialend';
upperend' 40*ones(size(upperend'))]/1e3);
% Dielectric ring
splid=5;
geomcells(splid).Dval=0; % V
geomcells(splid).order=2;
geomcells(splid).dim=2;
%splines(1).epsge=1e-6;
%splines(1).epsce=1e-6;
geomcells(splid).name='dielectric rings';
geomcells(splid).type=2;
geomcells(splid).periodic=0;
ow=flip(linspace(zbottom,zbottom+30,130));
lw=flip(linspace(10.1,23.5,30));
rw=linspace(zbottom,429.5,60);
uw=flip(linspace(0,10,10));
geomcells(splid).points=[flip([upperend', 40*ones(size(upperend'))]);
flip([(180+250)*ones(size(axialend')) axialend']);
flip([tiltedz' tiltedr']);
flip([base' 24*ones(size(base))']);
ow' 24*ones(size(ow'));
zbottom*ones(size(lw')) lw';
zwellipse' rwellipse';
(zwellipse(end)+0.5)*ones(size(uw')) uw']/1e3;
%% Plots
f=figure;
for k=1:length(geomcells)
plothandle=plot(geomcells(k).points(:,1), geomcells(k).points(:,2),'k-x','linewidth',1.5);
hold on
%geomcells(k).points=[geomcells(k).Z; geomcells(k).R];
order=geomcells(k).order;
knots=linspace(0,1,size(geomcells(k).points,1)-(order-2));
knots=augknt(knots, order);
sizec=size(geomcells(k).points,1);
order=length(knots)-sizec(end);
coeffs=geomcells(k).points';
pp=spmak(knots,coeffs);
s=linspace(0,1,100000);
fittedpos=fnval(pp,s);
plot(fittedpos(1,:),fittedpos(2,:),'x-')
end
legend(plothandle,{'Gun geometry'},'location','southwest')
f.PaperUnits='centimeters';
f.PaperSize=[12,8];
xlabel('z [m]')
ylabel('r [m]')
axis equal
savegeomtoh5('exp_both_with_vessel_geom.h5',geomcells,1e-2,true);
% print(f,name,'-dpdf','-fillpage')
% savefig(f,name)
% set(f, 'Color', 'w');
% export_fig(f,name,'-eps')
hold off
c4science · Help