AddLorentzMaterial

AddLorentzMaterial(CSX, name, varargin)

Full definition:

function CSX = AddLorentzMaterial(CSX, name, varargin)

Add a Drude/Lorentz type dispersive material model. Note: openEMS currently only supports a drude material type.

The drude type frequency dependent material:

eps_r(f) = eps_r* ( 1 - f_eps_plasma^2/(f*(f-j/t_eps_r)) )
mue_r(f) = mue_r* ( 1 - f_mue_plasma^2/(f*(f-j/t_mue_r)) )

with

  • f_eps_plasma: the respective electric angular plasma frequency

  • f_mue_plasma: the respective magnetic angular plasma frequency

  • t_eps_r: the respective electric relaxation time

  • t_mue_r: the respective magnetic relaxation time

Use SetMaterialProperty to define the material constants:

  • ‘EpsilonPlasmaFrequency’: electric plasma frequency (f_eps_plasma)

  • ‘MuePlasmaFrequency’: magnetic plasma frequency (f_mue_plasma)

  • ‘EpsilonRelaxTime’: electric plasma relaxation time (losses)

  • ‘MueRelaxTime’: magnetic plasma relaxation time (losses)

Note: all properties must be positive values

Higher order Drude type:

  • ‘EpsilonPlasmaFrequency_’: n-th order electric plasma frequency (f_eps_plasma)

  • ‘MuePlasmaFrequency_’: n-th order magnetic plasma frequency (f_mue_plasma)

  • ‘EpsilonRelaxTime_’: n-th order electric plasma relaxation time (losses)

  • ‘MueRelaxTime_’: n-th order magnetic plasma relaxation time (losses)

The Lorentz type frequency dependent material:

eps_r(f) = eps_r* ( 1 - f_eps_plasma^2/(f^2-f_eps_Lor_Pole^2-jf^2*/t_eps_r)) )
mue_r(f) = mue_r* ( 1 - f_mue_plasma^2/(f^2-f_mue_Lor_Pole^2-jf^2*/t_mue_r)) )

with the additional parameter (see above)

  • f_eps_Lor_Pole: the respective electric angular lorentz pole frequency

  • f_mue_Lor_Pole: the respective magnetic angular lorentz pole frequency

Use SetMaterialProperty to define the material constants:

  • ‘EpsilonLorPoleFrequency’: electric lorentz pole frequency (f_eps_Lor_Pole)

  • ‘MueLorPoleFrequency’: magnetic lorentz pole frequency (f_mue_Lor_Pole)

Note: all properties must be positive values

Higher order Drude type:

  • ‘EpsilonLorPoleFrequency_’: n-th order electric lorentz pole frequency (f_eps_plasma)

  • ‘MueLorPoleFrequency_’: n-th order magnetic lorentz pole frequency (f_mue_plasma)

example:

CSX = AddLorentzMaterial(CSX,'drude');
CSX = SetMaterialProperty(CSX,'drude','Epsilon',5,'EpsilonPlasmaFrequency',5e9,'EpsilonRelaxTime',1e-9);
CSX = SetMaterialProperty(CSX,'drude','Mue',5,'MuePlasmaFrequency',5e9,'MueRelaxTime',1e-9);
[..]
CSX = AddBox(CSX,'drude', 10 ,start,stop);

See also AddBox, AddMaterial, SetMaterialProperty

author: Thorsten Liebig