% 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 versionclear vrsnvrsn.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 2g1=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 1flclear fem% Analyzed geometryclear c sc.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 1fem.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 1clear applappl.mode.class = 'ElectromagneticWaves';appl.module = 'RF';appl.gporder = 4;appl.cporder = 2;appl.assignsuffix = '_rfw';clear propprop.solvefor='H';appl.prop = prop;clear edgedg.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 bndbnd.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 elementsclear elemcpl% Integration coupling variablesclear elemelem.elem = 'elcplscalar';elem.g = {'1'};src = cell(1,1);clear edgedg.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 femfem.sdim = {'x','y'};fem.border = 1;clear units;units.basesystem = 'SI';fem.units = units;xfem.fem{2} = fem;% Global expressionsxfem.globalexpr = {'Z','1/I'};% ODE Settingsclear odeclear units;units.basesystem = 'SI';ode.units = units;xfem.ode=ode;% Multiphysicsxfem=multiphysics(xfem);% Generate GMG mesh casesxfem=meshcaseadd(xfem,'mgauto','shape');% Extend meshxfem.xmesh=meshextend(xfem, ...                      'geoms',[1]);% Solve problemxfem.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 purposesfem0=xfem;% Plot solutionpostplot(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 casesxfem=meshcasedel(xfem,[1]);% Plot solutionpostplot(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 2fem=xfem.fem{2};% Geometry objectsclear p sp.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 versionclear vrsnvrsn.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 2fem=xfem.fem{2};% Geometry objectsclear p sp.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;