Homec4science
Diffusion FENNECS lite (3Dmpimode) Browse

fennecslite/geometries/experiment_arc/exp_arc_withvessel.m
6b37315661933Dmpimode

Last Change
Raw File

exp_arc_withvessel.m
View Options

%% Create the geometry 1 for the T-REX experiment
% uses the upper ellipse and constant radius center column.
C_column_length=400;%mm
C_column_radius=10; %mm
O_electrode_length=400;%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(40,rmax,8);
rw=linspace(250,480,20);
uw=flip(linspace(0,rmax,30));
geomcells(splid).points=flip([250*ones(size(lw')) lw';
rw' rmax*ones(size(rw'));
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='vacuum vessel';
geomcells(splid).type=0;
geomcells(splid).periodic=0;
lw=linspace(-.5,0,10)+C_column_radius;
rw=linspace(270.5,C_column_length-0.5,60);
uw=flip(linspace(0,1,10))*C_column_radius;
geomcells(splid).points=[ 270*ones(size(lw')) lw';
rw' C_column_radius*ones(size(rw'));
C_column_length*ones(size(uw')) uw']/1e3;
% Dielectric ring
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='vacuum vessel';
geomcells(splid).type=2;
geomcells(splid).periodic=0;
lw=flip(linspace(C_column_radius+.1,17.5,10));
rw=linspace(250,C_column_length-0.5,60);
uw=flip(linspace(0,C_column_radius,10));
geomcells(splid).points=[
250*ones(size(lw')) lw';
rw' C_column_radius*ones(size(rw'));
C_column_length*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;
rw=linspace(240,250,10);
lw=linspace(C_column_radius-.5,17.5,40);
uw=flip(linspace(240,250,10));
geomcells(splid).points=flip([
250*ones(size(lw')) lw';]/1e3);
% Outer electrode
splid=5;
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;
base=linspace(250,276,80);
%arc
radius=700; % mm
r_c=716.87399;
z_c=290.26217;
arcz=linspace(277,O_electrode_length);
arcr=r_c-sqrt(radius^2-(arcz-z_c).^2);
%tiltedwall
axialend=linspace(arcr(end)+0.5,39.5,15);
upperend=flip(linspace(250,arcz(end),60));
geomcells(splid).points=flip([base' 17*ones(size(base))';
arcz' arcr';
arcz(end)*ones(size(axialend')) axialend';
upperend' 40*ones(size(upperend'))]/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,1000);
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_arc_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