test_microstrip_line.m

% test_microstrip_line.m


% test_horn_antenna.m
warning off;

clc
clear
close all

pkg load geometry
pkg load miscellaneous

modelerPath = './Modeler_2D';
mesherPath = '../mesh2d';
meshWrapperPath = './Mesh_2D';
solver2DPath = './Solver_2D';
solver1DPath = './Solver_1D';
miscPath = './misc';

addpath(modelerPath);
addpath(genpath(mesherPath));
addpath(meshWrapperPath);
addpath(solver2DPath);
addpath(solver1DPath);
addpath(miscPath);

% Simulation properties
simProps = cem2D_createSimPropsStruct(...
                'fMin',0.1,...
                'fMax',5,...
                'polarizationType','TE');

                % Global mesh properties
meshProps = mesh2D_createMeshPropsStruct(...
              'relWLmeshMax',0.005, ...
              'boundingBoxAddSpace',1, ...
              'algorithmType','delfront',...
              'maxRadiusEdgeRatio',1.5);

% Unless overridden, f_sim will be set to the average between fMin and fMax
% f_sim = 3;

fc_target_GHz = 1.5;

% Line Properties
Line_W = 0.35;
Sub_Thick = 0.4;
Sub_W = Line_W*8;
Metal_Thick = 0.052;

warning off;

linePanel = mod2D_createRectangleStruct([-0.5*Line_W Metal_Thick+Sub_Thick],[0.5*Line_W Metal_Thick*2+Sub_Thick]);
gndPanel = mod2D_createRectangleStruct([-0.5*Sub_W 0],[0.5*Sub_W Metal_Thick]);

substrate = mod2D_createRectangleStruct([-0.5*Sub_W Metal_Thick],[0.5*Sub_W Metal_Thick+Sub_Thick]);

bbox = mod2D_createRectangleStruct([-0.5*Sub_W Metal_Thick+Sub_Thick],[0.5*Sub_W Metal_Thick + Sub_Thick*5]);

% Create the materials for this simulation
material_default = cem2D_createMaterialDefs;
material_PEC = cem2D_createMaterialDefs('name','PEC','type','PEC');
material_FR4 = cem2D_createMaterialDefs('name','FR4','type','normal','er',4.3);
% Create the material list with only the default material
materialList = cem2D_addMaterialToList(material_default);
materialList = cem2D_addMaterialToList(material_PEC,materialList);
materialList = cem2D_addMaterialToList(material_FR4,materialList);

% Create polygon list
polList = {substrate,linePanel,gndPanel};

% Subtract all shapes from bounding box
for polIdx = 1:numel(polList)
    bbox = mod2D_booleanOperation(bbox,polList{polIdx},'subtract');
end

% Add bounding box to polygon list
polList = [{bbox} polList];

% Assign materials per-shape
materialAssignement = {'default','FR4','PEC','PEC'};


%%%%%%%%%%%%%%%%%%%%%
% Plot the geometry %
%%%%%%%%%%%%%%%%%%%%%
figHdl = figure;
axHdl = axes;
mod2D_showPolygon(axHdl,bbox,[1 1 1],[0 0 0]);
mod2D_showPolygon(axHdl,substrate,[0 0.7 0],[0 0 0]);
mod2D_showPolygon(axHdl,linePanel,[255 165 0]/255,[0 0 0]);
mod2D_showPolygon(axHdl,gndPanel,[255 165 0]/255,[0 0 0]);

set(gca,'fontsize',14);
xlabel('\xi [mm]','fontsize',16);
ylabel('\eta [mm]','fontsize',16);
%%%%%%%%%%%%%%%%%%%%%
% Plot the geometry %
%%%%%%%%%%%%%%%%%%%%%

lineList = {};

% Test mesh of bounding box
meshData = mesh2D_generateInitialMesh(...
            polList,...             % Entire polygon list. This is mainly to calculate the maximum\minimum edge length
            lineList,...            % List of lines. Used for porst, mainly.
            materialAssignement,... % Material assignments for initial LFS assignment
            materialList,...        % Corresponding material properties
            meshProps,...           % Darrens Mesher properties
            simProps);              % Simulation properties. This is to assign spatial meshing rules, most

% Smooth the mesh
meshData = mesh2D_smoothMesh(meshData,meshProps);

if exist('f_sim','var')
    if isempty(f_sim)
        f_sim = 0.5*(simProps.fMin + simProps.fMax);
    end
else
    f_sim = 0.5*(simProps.fMin + simProps.fMax);
end
[lambda,Et,Ez,edgeTriplets] = cem2D_calcPortModes(meshData,meshProps,materialList,materialAssignement,simProps,f_sim);

k0 = 2*pi*2.55e9/3e8;
lambda_target = -(k0^2)*0.5*(4.3 + 1)

[~,sortIdxs] = sort(abs(lambda_target - lambda));
% [fc_TE,fc_TM,Et,Ez,edgeTriplets] = cem2D_calcModesByCutoff(meshData,meshProps,materialList,materialAssignement,simProps);

% Search for solutions with lowest cutoff.

c0 = physical_constant('speed of light in vacuum');

hold(axHdl,'on');

for tIdx = unique(meshData.tnum).'
    patch('faces',meshData.tria(meshData.tnum == tIdx,1:3),'vertices',meshData.vert, ...
        'facecolor','none', ...
        'edgecolor',[0,0,0]) ;
end

[Xm,Ym] = meshgrid(linspace(-0.5*Sub_W,0.5*Sub_W,25),linspace(0,Metal_Thick + Sub_Thick*5,15));

EI = cem2D_vectorElementInterp(...
    meshData.vert,...
    meshData.tria,...
    edgeTriplets,...
    Et(:,sortIdxs(1)));

Exy = EI(Xm,Ym);

%Eabs = sqrt(sum(abs(Exy).^2,2));



qHdl = quiver(Xm(:),Ym(:),abs(Exy(:,1)),abs(Exy(:,2)));
set(qHdl,'color',[1 0 0]);
%scatter(edgeCent(:,1),edgeCent(:,2),40,e0,'filled');

hold(axHdl,'off');

warning on;