magnetic_frill.m

Feeding antennas with proper signals can be difficult. The signal is often described as a voltage, a

% COMSOL Multiphysics Model M-file
% Generated by COMSOL 3.5 snapshot (COMSOL 3.5.0.468, $Date: 2008/08/15 16:04:11 $)

flclear xfem

% COMSOL version
clear vrsn
vrsn.name = 'COMSOL 3.5';
vrsn.ext = ' snapshot';
vrsn.major = 0;
vrsn.build = 468;
vrsn.rcs = '$Name:  $';
vrsn.date = '$Date: 2008/08/15 16:04:11 $';
xfem.version = vrsn;

% Geometry 2
g1=rect2('0.5','0.6','base','corner','pos',{'0','-0.3'},'rot','0');
g2=rect2('0.005*0.3','0.47*0.3','base','corner','pos',{'0','-0.47*0.15'},'rot','0');
g3=geomcomp({g1,g2},'ns',{'g1','g2'},'sf','g1-g2','edge','none');
parr={point2(0.0015,0)};
g4=geomcoerce('point',parr);
parr={point2(0.01,0)};
g5=geomcoerce('point',parr);
g6=revolve(g3,'angles',[0,-1.5707963267948966],'revaxis',[0 0;0 1],'wrkpln',[0 0 0;0 1 0;0 0 1]);
g7=revolve(g4,'angles',[0,-1.5707963267948966],'revaxis',[0 0;0 1],'wrkpln',[0 0 0;0 1 0;0 0 1]);
g8=revolve(g5,'angles',[0,-1.5707963267948966],'revaxis',[0 0;0 1],'wrkpln',[0 0 0;0 1 0;0 0 1]);

% Geometry 1
flclear fem

% Analyzed geometry
clear c s
c.objs={g7,g8};
c.name={'REV2','REV3'};
c.tags={'g7','g8'};

s.objs={g6};
s.name={'REV1'};
s.tags={'g6'};

fem.draw=struct('c',c,'s',s);
fem.geom=geomcsg(fem);

% Initialize mesh for geometry 1
fem.mesh=meshinit(fem, ...
                  'hmax',[0.07], ...
                  'hmaxfact',1.9, ...
                  'hcurve',0.8, ...
                  'hgrad',1.7, ...
                  'hcutoff',0.05, ...
                  'hnarrow',0.3, ...
                  'zscale',0.75, ...
                  'hcutofffac',[4,0.01,5,0.01,7,0.01,8,0.01], ...
                  'hcurvefac',[4,0.2,5,0.2,7,0.2,8,0.2]);
xfem.fem{1}=fem;

% (Default values are not included)

fem=xfem.fem{1};

% Application mode 1
clear appl
appl.mode.class = 'ElectromagneticWaves';
appl.module = 'RF';
appl.gporder = 4;
appl.cporder = 2;
appl.assignsuffix = '_rfw';
clear prop
prop.solvefor='H';
appl.prop = prop;
clear edg
edg.Im0 = {0,1};
edg.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,1,1,1,1,1,1];
appl.edg = edg;
clear bnd
bnd.srcpnt = {{0;0;0},{0;0;0},{0;0;0.3}};
bnd.type = {'E0','H0','SC'};
bnd.wavetype = {'plane','plane','cyl'};
bnd.ind = [2,2,3,1,1,3,1,1,3];
appl.bnd = bnd;
fem.appl{1} = appl;
fem.frame = {'ref'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;

% Coupling variable elements
clear elemcpl
% Integration coupling variables
clear elem
elem.elem = 'elcplscalar';
elem.g = {'1'};
src = cell(1,1);
clear edg
edg.expr = {{{},'4*(tHx*t1x+tHy*t1y)'}};
edg.ipoints = {{{},'4'}};
edg.frame = {{{},'ref'}};
edg.ind = {{'1','2','3','4','5','6','7','8','9','10','11','12','13', ...
  '15','16','17','18','19','20','21','22'},{'14'}};
src{1} = {{},edg,{},{}};
elem.src = src;
geomdim = cell(1,1);
geomdim{1} = {};
elem.geomdim = geomdim;
elem.var = {'I'};
elem.global = {'1'};
elemcpl{1} = elem;
fem.elemcpl = elemcpl;
xfem.fem{1} = fem;

flclear fem
fem.sdim = {'x','y'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;
xfem.fem{2} = fem;

% Global expressions
xfem.globalexpr = {'Z','1/I'};

% ODE Settings
clear ode
clear units;
units.basesystem = 'SI';
ode.units = units;
xfem.ode=ode;

% Multiphysics
xfem=multiphysics(xfem);

% Generate GMG mesh cases
xfem=meshcaseadd(xfem,'mgauto','shape');

% Extend mesh
xfem.xmesh=meshextend(xfem, ...
                      'geoms',[1]);

% Solve problem
xfem.sol=femstatic(xfem, ...
                   'complexfun','on', ...
                   'solcomp',{'tHxHyHz21','tHxHyHz20','tHxHyHz10'}, ...
                   'outcomp',{'tHxHyHz21','tHxHyHz20','tHxHyHz10'}, ...
                   'blocksize','auto', ...
                   'linsolver','gmres', ...
                   'prefuntype','right', ...
                   'prefun','gmg', ...
                   'prepar',{'presmooth','sorvec','presmoothpar',{'seconditer',2},'postsmooth','soruvec','postsmoothpar',{'seconditer',2},'csolver','pardiso'}, ...
                   'mcase',[0 1]);

% Save current fem structure for restart purposes
fem0=xfem;

% Plot solution
postplot(xfem, ...
         'slicedata',{'Wav_rfw','cont','internal','unit','J/m^3'}, ...
         'slicexspacing',5, ...
         'sliceyspacing',0, ...
         'slicezspacing',0, ...
         'slicemap','jet(1024)', ...
         'title','Slice: Total energy density, time average [J/m^3]', ...
         'grid','on', ...
         'campos',[-2.1945391483183814,-2.9357854429921457,2.318404662437931], ...
         'camtarget',[0.25,0.25,0], ...
         'camup',[0,0,1], ...
         'camva',12.051516421930526);

% Remove generated GMG mesh cases
xfem=meshcasedel(xfem,[1]);

% Plot solution
postplot(xfem, ...
         'slicedata',{'normE_rfw','cont','internal','unit','V/m'}, ...
         'slicedlim',[0 20], ...
         'slicexspacing',[0 0], ...
         'sliceyspacing',0, ...
         'slicezspacing',1, ...
         'slicemap','jet(1024)', ...
         'title','Slice: Electric field, norm [V/m]', ...
         'grid','on', ...
         'campos',[-2.1945391483183814,-2.9357854429921457,2.318404662437931], ...
         'camtarget',[0.25,0.25,0], ...
         'camup',[0,0,1], ...
         'camva',12.051516421930518);

% Geometry 2
fem=xfem.fem{2};

% Geometry objects
clear p s
p.objs={g4,g5};
p.name={'PT1','PT2'};
p.tags={'g4','g5'};

s.objs={g3};
s.name={'CO1'};
s.tags={'g3'};

fem.draw=struct('p',p,'s',s);
xfem.fem{2}=fem;

% COMSOL Multiphysics Model M-file
% Generated by COMSOL 3.5 (COMSOL 3.5.0.494, $Date: 2008/09/19 16:09:48 $)

% COMSOL version
clear vrsn
vrsn.name = 'COMSOL 3.5';
vrsn.ext = '';
vrsn.major = 0;
vrsn.build = 494;
vrsn.rcs = '$Name:  $';
vrsn.date = '$Date: 2008/09/19 16:09:48 $';
xfem.version = vrsn;

% Geometry 2
fem=xfem.fem{2};

% Geometry objects
clear p s
p.objs={g5,g4};
p.name={'PT2','PT1'};
p.tags={'g5','g4'};

s.objs={g3};
s.name={'CO1'};
s.tags={'g3'};

fem.draw=struct('p',p,'s',s);
xfem.fem{2}=fem;