Homec4science
Diffusion FENNECS Lite (3D) Browse

fennecslite/geometries/paschen_curves/paschen_geom.m
6b37315661933D

Last Change
Raw File

paschen_geom.m
View Options

%% Create the geometry 1 for the T-REX experiment
% uses the upper ellipse and constant radius center column.
clear all
close all
zbottom=0; % mm
ztop=50;
rmax=200; %mm
% Top plate
splid=1;
geomcells(splid).Dval=0; % V
geomcells(splid).order=3;
geomcells(splid).dim=2;
geomcells(splid).name='top plate';
geomcells(splid).type=0;
geomcells(splid).periodic=0;
uw=flip(linspace(0,rmax,30));
geomcells(splid).points=flip([ztop*ones(size(uw')) uw']/1e3);
%{
% Outer electrode
splid=2;
geomcells(splid).Dval=0; % V
geomcells(splid).order=3;
geomcells(splid).dim=2;
geomcells(splid).name='vacuum vessel';
geomcells(splid).type=2;
geomcells(splid).periodic=0;
rmax=60; %mm
rw=linspace(zbottom,ztop,20);
geomcells(splid).points=flip([rw' rmax*ones(size(rw'))]/1e3);
% Central electrode
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='central electrode';
geomcells(splid).type=2;
geomcells(splid).periodic=0;
%ellipse=
w=linspace(zbottom,ztop,60);
geomcells(splid).points=[ w' 10*ones(size(w'))]/1e3;
%}
% Bottom plate
splid=2;
geomcells(splid).Dval=0; % V
geomcells(splid).order=3;
geomcells(splid).dim=2;
geomcells(splid).name='bottom plate';
geomcells(splid).type=0;
geomcells(splid).periodic=0;
lw=linspace(0,rmax,30);
geomcells(splid).points=flip([ zbottom*ones(size(lw')) lw']/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('paschen_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