Question about this warning

[kmerbcn]

Hello

I have a problem with a simulation. I think I have done something wrong in the set-up. I have the following warning and once the model is computed, the electric fields are not seen.
What's wrong?

[kmerbcn]

Hello

I have solved the problem related with geometric entities. Nevertheless, I still have a problem dealing with the excitation port

The excitation ports are in the xz-plane and the vector direction of electric field is also z-direction.

Do you want the xml file?

Br

Javier

[thliebig]

The excitation is not meshed as it seems. No energy is inserted.
The mesh looks very coarse and not very smooth too...

[kmerbcn]

Now the geometry is finer. The mesh looks as follows

The code I followed is the one coming from the macro developed in Freecad by Lubomir Jagos.
The structure of code was already working in for other microstrip filter I did couple of months ago. First of all, I was doing the structure in FreeCad in 3D Later on, I did a face using LumpedPort directive by taking the needed points in FreeCad and it was working.

I can give to you the structure of the code as well. The extension of the file was changed from matlab to txt
Filtro_Coupled_Clones_openEMS.txt

[kmerbcn]

With the same script structure, i do not have same issues as before. Look at this example.

This one was working with the following script sctructure

close all
clear
clc

%% Change the current folder to the folder of this m-file.
if(~isdeployed)
mfile_name = mfilename('fullpath');
[pathstr,name,ext] = fileparts(mfile_name);
cd(pathstr);
end

%% constants
physical_constants;
unit = 0.001; % Model coordinates and lengths will be specified in mm.
fc_unit = 0.001; % STL files are exported in FreeCAD standard units (mm).

%% switches & options
postprocessing_only = 0;
draw_3d_pattern = 0; % this may take a while...
use_pml = 0; % use pml boundaries instead of mur

currDir = strrep(pwd(), '', '\');
display(currDir);

% --no-simulation : dry run to view geometry, validate settings, no FDTD computations
% --debug-PEC : generated PEC skeleton (use ParaView to inspect)
openEMS_opts = '--disable-dumps';

%% prepare simulation folder
Sim_Path = 'simulation_output';
Sim_CSX = 'Filtro_Calculado_Stepped_Distributed_tx_v2.xml';
[status, message, messageid] = rmdir( Sim_Path, 's' ); % clear previous directory
[status, message, messageid] = mkdir( Sim_Path ); % create empty simulation folder

%% setup FDTD parameter & excitation function
max_timesteps = 3000000;
min_decrement = 1e-04; % 10*log10(min_decrement) dB (i.e. 1E-5 means -50 dB)
FDTD = InitFDTD( 'NrTS', max_timesteps, 'EndCriteria', min_decrement, 'OversSampling', 100);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% BOUNDARY CONDITIONS
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
BC = {"MUR","MUR","MUR","MUR","MUR","MUR"};
FDTD = SetBoundaryCond( FDTD, BC );

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% COORDINATE SYSTEM
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
CSX = InitCSX('CoordSystem', 0); % Cartesian coordinate system.
mesh.x = []; % mesh variable initialization (Note: x y z implies type Cartesian).
mesh.y = [];
mesh.z = [];
CSX = DefineRectGrid(CSX, unit, mesh); % First call with empty mesh to set deltaUnit attribute.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% EXCITATION gauss
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
f0 = 25.01000000000.0;
fc = 20.01000000000.0;
FDTD = SetGaussExcite( FDTD, f0, fc );
max_res = c0 / (f0 + fc) / 20;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% MATERIALS AND GEOMETRY
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
CSX = AddMetal( CSX, 'PEC' );

%% MATERIAL - Aire
CSX = AddMaterial(CSX, 'Aire');
CSX = SetMaterialProperty(CSX, 'Aire', 'Epsilon', 1, 'Mue', 1);
CSX = ImportSTL(CSX, 'Aire', 9600, [currDir '/Aire_gen_model.stl'], 'Transform', {'Scale', fc_unit/unit});

%% MATERIAL - PEC
CSX = AddMetal(CSX, 'PEC');
CSX = ImportSTL(CSX, 'PEC', 9900, [currDir '/Pista_Salida002_gen_model.stl'], 'Transform', {'Scale', fc_unit/unit});
CSX = ImportSTL(CSX, 'PEC', 9800, [currDir '/Masa002_gen_model.stl'], 'Transform', {'Scale', fc_unit/unit});
CSX = ImportSTL(CSX, 'PEC', 10000, [currDir '/Pista_Entrada002_gen_model.stl'], 'Transform', {'Scale', fc_unit/unit});

%% MATERIAL - ROGERS_5880
CSX = AddMaterial(CSX, 'ROGERS_5880');
CSX = SetMaterialProperty(CSX, 'ROGERS_5880', 'Epsilon', 2.2, 'Mue', 1);
CSX = ImportSTL(CSX, 'ROGERS_5880', 9700, [currDir '/Sustrato002_gen_model.stl'], 'Transform', {'Scale', fc_unit/unit});

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% GRID LINES
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%% GRID - xyz_aire - Aire (Fixed Distance)
mesh.x(mesh.x >= -10 & mesh.x <= 40) = [];
mesh.x = [ mesh.x (-10:1:40) ];
mesh.y(mesh.y >= -24 & mesh.y <= 26) = [];
mesh.y = [ mesh.y (-24:1:26) ];
mesh.z(mesh.z >= -5 & mesh.z <= 5) = [];
mesh.z = [ mesh.z (-5:1:5) ];
CSX = DefineRectGrid(CSX, unit, mesh);

%% GRID - xy_masa_sustrato - Masa002 (Fixed Distance)
mesh.x(mesh.x >= 0 & mesh.x <= 27.1) = [];
mesh.x = [ mesh.x (0:0.3:27.1) ];
mesh.y(mesh.y >= -5 & mesh.y <= 5) = [];
mesh.y = [ mesh.y (-5:0.3:5) ];
CSX = DefineRectGrid(CSX, unit, mesh);

%% GRID - xy_masa_sustrato - Sustrato002 (Fixed Distance)
mesh.x(mesh.x >= 0 & mesh.x <= 27.1) = [];
mesh.x = [ mesh.x (0:0.3:27.1) ];
mesh.y(mesh.y >= -5 & mesh.y <= 5) = [];
mesh.y = [ mesh.y (-5:0.3:5) ];
CSX = DefineRectGrid(CSX, unit, mesh);

%% GRID - z_sustrato - Sustrato002 (Fixed Count)
mesh.z(mesh.z >= -0.28 & mesh.z <= 0) = [];
mesh.z = [ mesh.z linspace(-0.28,0,4) ];
CSX = DefineRectGrid(CSX, unit, mesh);

%% GRID - z_pistas - Pista_Salida002 (Fixed Count)
mesh.z(mesh.z >= 0 & mesh.z <= 0.035) = [];
mesh.z = [ mesh.z 0.0175 ];
CSX = DefineRectGrid(CSX, unit, mesh);

%% GRID - z_pistas - Pista_Entrada002 (Fixed Count)
mesh.z(mesh.z >= 0 & mesh.z <= 0.035) = [];
mesh.z = [ mesh.z 0.0175 ];
CSX = DefineRectGrid(CSX, unit, mesh);

%% GRID - xy_pistas - Pista_Salida002 (Fixed Distance)
mesh.x(mesh.x >= 0 & mesh.x <= 26.038) = [];
mesh.x = [ mesh.x (0:0.1:26.038) ];
mesh.y(mesh.y >= -1.125 & mesh.y <= 0.875) = [];
mesh.y = [ mesh.y (-1.125:0.1:0.875) ];
CSX = DefineRectGrid(CSX, unit, mesh);

%% GRID - xy_pistas - Pista_Entrada002 (Fixed Distance)
mesh.x(mesh.x >= 12.838 & mesh.x <= 27.1) = [];
mesh.x = [ mesh.x (12.838:0.1:27.1) ];
mesh.y(mesh.y >= -1.511 & mesh.y <= -0.691) = [];
mesh.y = [ mesh.y (-1.511:0.1:-0.691) ];
CSX = DefineRectGrid(CSX, unit, mesh);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% PORTS
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

portNamesAndNumbersList = containers.Map();

d = 0.28;
V = 1;
E = V/d;

%% PORT - portin - portin
portStart = [0, -0.25, -0.28];
portStop = [0, 0, 0.035];
portR = 50;
portUnits = 1;
portExcitationAmplitude = E;
portDirection = [0 0 1]portExcitationAmplitude;
[CSX port{1}] = AddLumpedPort(CSX, 10000, 1, portRportUnits, portStart, portStop, portDirection, true);
portNamesAndNumbersList("portin") = 1;

portStart = [27.1, -1.511, -0.28];
portStop = [27.1, -0.691, 0.035];
portR = 50;
portUnits = 1;
portExcitationAmplitude = E;
portDirection = [0 0 1]portExcitationAmplitude;
[CSX port{2}] = AddLumpedPort(CSX, 11000, 2, portRportUnits, portStart, portStop, portDirection, false);
portNamesAndNumbersList("portin") = 2;

%% PROBE - efield - efield
CSX = AddDump(CSX, 'efield_efield', 'DumpType', 0, 'DumpMode', 2);
dumpboxStart = [ 0, -5, -0.05 ];
dumpboxStop = [ 27.1, 5, -0.05 ];
CSX = AddBox(CSX, 'efield_efield', 0, dumpboxStart, dumpboxStop );

%% PROBE - hfield - hfield
CSX = AddDump(CSX, 'hfield_hfield', 'DumpType', 1, 'DumpMode', 2);
dumpboxStart = [ 0, -5, -0.08 ];
dumpboxStop = [ 27.1, 5, -0.08 ];
CSX = AddBox(CSX, 'hfield_hfield', 0, dumpboxStart, dumpboxStop );

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% RUN
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
WriteOpenEMS( [Sim_Path '/' Sim_CSX], FDTD, CSX );
CSXGeomPlot( [Sim_Path '/' Sim_CSX] );

if (postprocessing_only==0)
%% run openEMS
RunOpenEMS( Sim_Path, Sim_CSX, openEMS_opts );
end

[kmerbcn]

Hello guys!

Now it seems to be solved. I checked the xy coordinates in the grid and the ones referring the port and it seems that in the y-coordinates, there was a small decimal number not included. Thanks and sorry for my lack of experience!