TESIS
TESIS
TESIS
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
3.1 Introducción<br />
CAPÍTULO III<br />
ALGORITMO DE YEE<br />
Las ecuaciones de Maxwell son ecuaciones diferenciales parciales de segundo orden que<br />
describen la presencia de ondas electromagnéticas. La necesidad de encontrar su solución<br />
computacional ha crecido en el transcurso del tiempo y existen en la actualidad diferentes<br />
métodos numéricos que pretenden dar solución a este problema. El método de Diferencias<br />
Finitas en el Dominio del Tiempo (FDTD) es una técnica numérica que resuelve las<br />
ecuaciones de Maxwell en forma diferencial en el dominio del tiempo y espacio.<br />
En este capítulo se presenta la construcción del algoritmo de FDTD comenzando con el<br />
desarrollo de la técnica de Diferencias Finitas, la cual se utiliza para solucionar en forma<br />
discreta las ecuaciones diferenciales. Se definen las ecuaciones de Maxwell en forma<br />
discreta para un espacio rectangular de tres y dos dimensiones.<br />
Para verificar que las soluciones del método numérico converjan a la solución real se<br />
realiza un análisis de estabilidad numérica definido por Courant, Friedrich y Levy (CFL) y<br />
Von Neumann [2]. Se presenta una comparación de la velocidad de propagación de la onda<br />
dentro del método numérico con el real y por último se definen el conjunto de fuentes<br />
utilizadas en este método numérico.<br />
3.2 Diferencias Finitas<br />
El método de diferencias finitas es un método numérico que se utiliza para resolver<br />
ecuaciones diferenciales parciales en forma discreta. Esta técnica consiste en reemplazar las<br />
derivadas parciales por una ecuación definida como “diferencias finitas” aproximada que si<br />
bien no cumple exactamente con la ecuación diferencial, desde el punto de vista práctico se<br />
toma como tal.<br />
Las fórmulas de diferencias finitas son obtenidas por medio de la expansión de las series de<br />
Taylor. Considerando la derivada parcial<br />
∂ F(<br />
x,<br />
t)<br />
, fijando el valor de x y realizando la<br />
∂t<br />
1 1<br />
aproximación en dos puntos t + ∆t<br />
y t − ∆t<br />
se tiene:<br />
2 2<br />
.<br />
y<br />
2<br />
3<br />
⎛ 1 ⎞<br />
∆t<br />
∆t<br />
1 ∆t<br />
1<br />
F ⎜ x,<br />
t + ∆t<br />
⎟ = F(<br />
x,<br />
t)<br />
+ F′<br />
( x,<br />
t)<br />
+ F ′′ ( x,<br />
t)<br />
⋅ + F ′′ ′ ( x,<br />
t)<br />
⋅ + ... (3.2.1a)<br />
⎝ 2 ⎠<br />
2 4 2!<br />
8 3!<br />
18