TESIS

σ
m
B
LADO L1
= σ max
⎡ σ ∆t
⎤
⎢
1−
2ε
⎥
Cexp1 = ⎢ ⎥
⎢ σ ∆t
1+
⎥
⎢⎣
2ε
⎥⎦
Chyp1
⎛ x ⎞
⎜ ⎟⎠
⎝ PML
Para L1 σy = σ * y = 0, por lo tanto la condición de
∗
acoplamiento es:
σ σ
= y se calculan las
ε 0 µ 0
componentes dirigidas en dirección “x”: Ezx y Hy x x
σ → Conductividad eléctrica para la zona PML.
Este valor va disminuyendo conforme se acerca al
PEC.
∗ ⎡ σ ⎤
σ = µ 0 ⋅ ⎢ ⎥
σ
⎣ε
0 ⎦
* → Conductividad magnética para la zona PML.
D
m
⎡ ∆t
⎤
⎢ ∆ε
⎥
Cexp2 = ⎢ ⎥
⎢ σ ∆t
1+
⎥
⎢⎣
2ε
⎥⎦
∗ ⎡ σ ∆t
⎤ ⎡ ∆t
⎤
⎢1−
⎥
⎢
2 µ
⎢ ∆ µ ⎥
⎥ Chyp2 = ⎢ ⎥
⎢ σ ∆t
⎥
⎢
1+
⎢ σ ∆t
1+
⎥
⎥
⎣ 2 µ
⎢ ⎥
⎦ ⎣ 2 µ ⎦
= ∗
Cex1=Cexp1
Cex2=Cexp2
Chy1=Chyp1
Chy2=Chyp2
Cálculo de las constantes eléctricas para el
lado L1.
Igualar las matrices para L1 con la “Matriz de
constantes eléctricas 1”
85