1 files changed, 60 insertions, 0 deletions
diff --git a/CSXCAD/matlab/AddLorentzMaterial.m b/CSXCAD/matlab/AddLorentzMaterial.m
new file mode 100644
index 0000000..b0e9629
--- /dev/null
+++ b/CSXCAD/matlab/AddLorentzMaterial.m
@@ -0,0 +1,60 @@
+function CSX = AddLorentzMaterial(CSX, name, varargin)
+% 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>':  n-th order electric plasma frequency (f_eps_plasma)
+%   'MuePlasmaFrequency_<n>':      n-th order magnetic plasma frequency (f_mue_plasma)
+%   'EpsilonRelaxTime_<n>':        n-th order electric plasma relaxation time (losses)
+%   'MueRelaxTime_<n>':            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>':  n-th order electric lorentz pole frequency (f_eps_plasma)
+%   'MueLorPoleFrequency_<n>':      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
+%
+% CSXCAD matlab interface
+% -----------------------
+% author: Thorsten Liebig
+
+CSX = AddProperty(CSX, 'LorentzMaterial', name, varargin{:});