TESIS
TESIS
TESIS
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
A continuación se describe la técnica de PML considerando un espacio de dos dimensiones<br />
y las polarizaciones TE y TM.<br />
4.3.1 PML en un espacio de 2D, modo TE<br />
Considere el conjunto de ecuaciones para un espacio de 2D del grupo TE representadas en<br />
la Tabla 3.1. Los términos σ y ∗<br />
σ son las conductividades eléctricas y magnéticas<br />
respectivamente y representan pérdidas en el espacio libre.<br />
Si una onda electromagnética que se propaga en el espacio libre e incide en un espacio con<br />
pérdidas eléctricas y magnéticas, su comportamiento va a ser de reflexión, sin embargo, si<br />
se cumple la relación:<br />
∗<br />
σ σ<br />
= (4.3.1)<br />
ε µ<br />
o<br />
donde ε o es la permitividad eléctrica en el vacío y µ 0 es la permeabilidad magnética en el<br />
vacío, la impedancia en el vacío es igual a la impedancia de un espacio con perdidas<br />
eléctricas y magnéticas, lo que equivale a eliminar la posibilidad de reflexión [4].<br />
La relación (4.3.1) evita la posibilidad de que una onda que viaja en el espacio vacío e<br />
incide normalmente en un medio con pérdidas eléctricas y magnéticas sufra reflexión, pero<br />
aun así existe el problema de las ondas incidentes en forma oblicua, ya que con ellas no se<br />
satisface esta relación.<br />
Berenguer tuvo la idea de separar las componentes de campo eléctrico y magnético de las<br />
celdas que forman el límite del espacio numérico y con ello considerar las ondas incidentes<br />
oblicuas. Para el grupo TE la componente z E es dividida en dos componentes: E zx y E zy y<br />
ello representa el siguiente conjunto de ecuaciones:<br />
o<br />
( E + )<br />
∂H x ∗ ∂<br />
+ σ yH<br />
x = − zx Ezy<br />
µ0 (4.3.2a)<br />
∂t<br />
∂y<br />
( E + E )<br />
∂H y ∗<br />
µ0 + σ xH<br />
y<br />
∂t<br />
∂<br />
=<br />
∂x<br />
zx zy<br />
(4.3.2b)<br />
∂Ezx<br />
ε0 + σ xEzx<br />
∂t<br />
∂H<br />
y<br />
=<br />
∂x<br />
(4.3.2c)<br />
∂Ezy<br />
+ σ yEzy<br />
∂t<br />
∂H<br />
x = −<br />
∂y<br />
ε 0 (4.3.2d)<br />
38