Sillabus
Sillabus
Sillabus
FACULTAT DE FISICA
In the unlikely event that the author did not send a complete manuscript
and there are missing pages, these will be noted. Also, if material had to be removed,
a note will indicate the deletion.
Disscrrlation Püblish<¡ng
UMI U607741
Published by ProQuest LLC 2014. Copyright in the Dissertation held by the Author.
Microform Edition © ProQuest LLC.
All rights reserved. This work is protected against
unauthorized copying underTitle 17, United States Code.
ProQuest LLC
789 East Eisenhower Parkway
P.O. Box 1346
Ann Arbor, MI 48106-1346
UNIVERSITAT DE VALENCIA
BIBLIOTECA CIÉN C IES
N? Registre .
D A T A ____
Por últim o, con entrañable cariño, a mis padres, herm anos y Bea, que
h an tenido que aguantar m is estadios de m al hum or, y tam bién a mis
amigos que siem pre han valorado mi labor.
gag
Capítulo 1: Introducción
1.1. Introducción. 1
1.2. M étodos aplicables a dispositivos pasivos. 2
1.3. Orígenes del m étodo y limitaciones del pasado. 9
1.4. A ntecedentes. 10
1.5. Posibilidades futuras de DFDT. 20
1.6. Organización de esta tesis. 22
2.1. Introducción. 25
2.2. Las ecuaciones de Maxwell. 26
2.3. Diferencias finitas. 28
2.3.1. A proxim ación del operador derivada. 28
2.3.2. Ecuaciones diferenciales en derivadas parciales. 31
2.3.3. Estabilidad y convergencia. 32
2.3.4. El algoritm o de Yee. 33
2.3.5. Estabilidad del algoritmo de Yee. 40
2.3.6. M odelado. 45
2.3.7. D ispersión del algoritm o de Yee. 45
2.4. Excitación. 52
2.4.1. Separación entre zona de campo total y zona de campo
dispersado. 53
2.5. Condiciones de contorno de radiación o condiciones de
contorno absorbentes. 67
2.5.1. Ecuación de onda de u n solo caminó. 68
2.5.2. Im plem entación num érica de las CCA. 78
2.5.3. CCA de M ur de prim er orden. 82
i
pag
2.5.4. Análisis de la reflexión (CCA). 84
2.6. Simplificaciones com putacionales. 89
2.6.1. N um eración de los vectores. 90
2.6.2. Medios con pérdidas. 96
2.6.3. Región del espacio con un nfl finito de m edios distintos. 98
ii
Capítulo 4: Aplicación de DFDT al análisis
de discontinuidades en el plano H
iii
pag
5.3.2. M odo fundam ental (TEM). 197
5.3.3. Modelización de la guía coaxial cuadrada (DFDT). 207
5.3.4. Com portam iento de las CCA con el m odo fundam ental. 208
5.3.5. Estudio de la atenuación de los m odos de orden
superior por debajo del corte. 212
5.3.6. Excitación con señal tem poral sintética con estructura
espacial TEM. Separación cam po total-campo reflejado. 224
5.3.7. M allado variable en la m odelización del coaxial. 229
5.4. Análisis de una guía coaxial rellena de dieléctrico. 231
5.5. Análisis de una unión T simétrica en guía coaxial cuadrada. 235
5.5.1. Modelado. 236
5.5.2. Esquema del program a principal. 242
5.5.3. Resultados num éricos. 244
5.6. Conclusiones. 262
APENDICES. 265
CONCLUSIONES Y PRINCIPALES APORTACIONES. 269
BIBLIOGRAFIA. 275
iv
Capítulo 1: Introducción
1.1. Introducción.
M ucho de lo que v u lg arm en te llam am os m o d ern a tecnología
depende en g ran parte del análisis y síntesis de sistemas de radiofrecuencia.
Esta ha sufrido en los últim os tiem pos un gran desarrollo en campos tan
diversos como pued en ser telecom unicaciones, radar, geofísica, m edicina y
defensa. La posibilidad de com prender, conocer y controlar los efectos de la
radiación electrom agnética se ha revelado de gran im portancia. Todo el
conocimiento adquirido y todo el futuro a generar en este campo se basan en
la construcción de soluciones precisas de las ecuaciones de Maxwell.
2
rep resen tativ o s de los utilizados norm alm ente en la caracterización de
dispositivos pasivos de microondas.
Diferencias Finitas.
3
T am bién recientem ente, se ha d e sarro llad o el m éto d o de los
Elementos de Contorno [7,8] . Este es una com binación del m étodo de la
Ecuación Integral de Contorno y de una discretización sim ilar a la que se
efectúa con elem entos finitos, pero siem pre en dos dim ensiones (solo se
discretiza el contorno).
Con este m étodo [10,11], la región a estudiar es sim ulada con una red
tridim ensional de lineas de transm isión. Con este m étodo se sim ula en el
dom inio del tiem po el fenóm eno de la propagación de ondas. Es m uy
versátil. Es sim ilar al m étodo de Diferencias Finitas en el D om inio del
tiem po con sus m ism as ventajas e inconvenientes, salvo que en lugar de
sim ular cam pos con diferencias divididas, se sim ulan corrientes en lineas
de transm isión [12].
4
Método de los Momentos.
M
M ode m atching".
Este se aplica típicam ente cuando hay una discontinuidad entre dos
sistem as g u ia d o re s d iferen tes. El cam po en am bos lad o s de la
d isco n tin u id ad es desarrollado como com binación lineal de los m odos
propios de las respectivas regiones [15]. Los coeficientes de los desarrollos
son in d eterm in ad o s. Estos se hallan exigiendo la co n tin u id ad de los
campos en la superficie de separación que define la discontinuidad.
Resonancia transversal.
En este m étodo, dos de las tres dim ensiones son discretizadas. ,E1
p ro b lem a se resu e lv e a n alíticam en te en la terc era d im e n sió n y
num éricam ente en las otras dos [18]. Se suele utilizar para el análisis de los
m odos en guias de sección complicada.
6
Método de la matriz de dispersión generalizada.
Las técnicas espectrales son num éricam ente eficientes, pero su rango
de aplicación es lim itado. A continuación se m uestra en la table 1.1. varios
aspectos com parados de los m étodos num éricos más usuales. La evaluación
que se efectúa de estos es cualitativa no habiendo una clara separación entre
'‘m oderado" y "mucho". D entro de cada m étodo hay m uchas variaciones,
siendo m uy im portante la experiencia del investigador para generar trucos,
y acelerar la obtención de resultados. Por esa razón la tabla 1.1. sólo sirve de
u n m odo grosero para tal comparación.
7
TABLA 1.1
Diferencias Finitas G G MB N u lo
Elementos de Contorno M M MB P
M atriz de línea de
transm isión (TLM) MG MG MB P
Ecuación integral PM PM G M
"Mode m atching" M PM G M
Resonancia transversal PM PM Ma M
G = Grande.
M = M oderado.
P = Poco.
MB = M uy Bueno.
B = Bueno.
Ma = Marginal.
8
1.3. Orígenes del Método y limitaciones del pasado.
A m ed iad o s de los años 60, K.S. Yee [21] in tro d u ce u n nuevo
algoritm o p ara resolver las ecuaciones de M axw ell del rotacional, en el
dom inio del tiem po. Con este m étodo el cam po electrom agnético en el
continuo, es m uestreado en distintos puntos de un m allado espacial y en
p u n to s distintos pero equiespaciados de un m allado tem poral. La evolución
del cam po electrom agnético es m odelizada de una form a autoconsistente,
calculando sucesivam ente paso a paso tem poral, siguiendo la analogía
num érica de las ecuaciones de Maxwell en cada nodo de la red. El algoritmo
asi construido sim ula el avance de las ondas m ediante el avance del campo
inicial -m uestreado- en la m em oria del ordenador. En este esquem a, el
increm ento de espacio y tiempo en el m uestreo son escogidos para evitar el
"aliasing" de la distribución continua de cam po, y garantizar la estabilidad
num érica del algoritm o.
9
c) N o se podían describir estructuras m as pequeñas que una celda.
Esto no perm itía m odelizar con suficiente precisión estructuras de alambres
d e lg a d o s, p e q u e ñ a s h e n d id u ra s etc... Esto d ab a lu g a r a g ru esas
m odelizaciones que en caso de contornos curvos tenía aspecto escalonado.
1.4. Antecedentes.
Yee [21] en 1966 propone discretizar las ecuaciones de M axwell del
rotacional para resolverlas por diferencias finitas. Para ello define un cubo
u n id ad en el que en cada nodo hay definida una com ponente del campo. En
su algoritm o, las com ponentes del campo eléctrico y m agnético se calculan
sucesivam ente en intervales de tiem po alternativos. En el m ism o trabajo, a
m odo de ejem plo, m uestra el "scattering" de u n pulso sobre u n cilindro
cuadrado presentando los perfiles del campo dispersado en varios instantes
tem porales.
10
tem poral. F inalm ente el esquem a es aplicado para an alizar el cam po
dispersado por un cilindro dieléctrico de sección circular.
11
ventaja de utilizar el campo total por el ruido que queda al efectuar la resta.
V alida el algoritm o tratando una esfera dieléctrica y una cavidad cilindrica
conductora con una abertura. Finalm ente analiza la penetración del campo
en la carcasa de la cabeza guiadora de un misil. La conclusión es un nivel de
precisión en los resultados de mas menos ldB.
12
Con otros trabajos de Um ashankar y Taflove [30,31] (1982), [32] (1983),
se continúa con la aplicación de DFDT para el cálculo del cam po dispersado
y RCS (sección recta radar) de objetos en dos y tres dim ensiones. En estos se
hace uso de las innovaciones introducidas por M ur, pero con una novedad,
calculan el cam po dispersado en el infinito a p a rtir de las corrientes
equivalentes que se obtienen del campo dispersado próxim o. De esta forma
calculan la RCS y las corrientes inducidas en los blancos com parando con
resultados que obtienen por el m étodo de los M om entos. En el últim o de
estos trabajos hacen cálculos para objetos de hasta nueve longitudes de onda
y dan resultados de tiempo de cálculo en varios ordenadores [32] (1985). La
conclusión de todos es que la precisión que da DFDT es al m enos
com parable a la que da el m étodo de los m om entos o cualquier otro de los
m étodos existentes.
13
problem as de absorción de cam po por tejidos biológicos (SAR: Specific
A bsortion Rate). A nalizan cilindros con p érd id a s en dos dim ensiones y
en cu en tran que FFT-CGM da errores con polarización TE m ientras que
DFDT da excelentes resultados con TE y TM. Su conclusión es que DFDT es
m ejor que FFT-CGM.
Uno de los prim eros trabajos que aborda antenas con DFDT se debe a
Reineix y Jecko [41], (1989). En este se estudian antenas de parche, analizando
ciertos parám etros en función de la frecuencia: corrientes, im pedancia de
entrada y diagram as de radiación.
14
ondas para m odos TElO y T M ll, com parando con resultados analíticos
(precisión de alrededor de un 1%), aplica luego el algoritm o a un caso no
trivial sin solución analítica (aplicador en guía de ondas cargado con
dieléctrico) y compara con resultados obtenidos por elementos finitos.
15
de ordenador, con lo que puede efectuar cálculos con un sim ple ordenador
IBM PC-AT. C oncluye que la excitación m as conveniente es una delta
tem poral y que el m étodo es insensible a los errores de redondeo, con lo que
puede usar núm eros reales de baja precisión ahorrando m as recursos.
16
interesante lanzar la excitación con la estructura del m odo fundam ental
p a ra hacer que este alcance antes una estructura estable, pero que en
cualquier caso esto no afecta a la distribución del cam po. U tilizan las
condiciones absorbentes de M ur modificadas con el tratam iento de Fang [50]
-para m inim izar reflexiones que introducen error, esta es la causa según los
au to res de que el m étodo del dom inio del tiem po sólo se utilize para
obtener datos cualitativos y no parám etros en función de la frecuencia y
d atos de diseño-, y hacen el m ism o tratam iento que De Pourq en las
in te rfases dieléctricas. C alculan d e n sid ad e s de co rrien te, constantes
dieléctricas efectivas, im pedancias características y la relación entre
com ponentes del cam po tranversales y lo n g itu d in ales p ara ver a que
frecuencia desaparece la estructura cuasi-TEM.
A finales de 1990, Bi, Wu y Litva, [55], aplican DFDT para analizar una
disco n tin u id ad en guía de onda rectangular. Se trata de un codo recto con
17
u n truncam iento. U tilizan señales m onocrom áticas y señalan la dificultad
en la obtención de buenos resultados en estructuras conductoras cerradas. La
justificación de esto la atribuyen a su elevado factor de calidad, y en la
m odelización con DFDT como estructuras conductoras perfectas.
Jarem [56], (1991), calcula por el m étodo de los m om entos y con DFDT
la im p ed an cia de en trad a de una sonda de alim entación de una guía
rectangular que ha sido cortocircuitada por un lado. Hace uso del program a
de DFDT, TSARS (Temporal Scattering and Reflectance Software- Lawrence
Liverm ore N ational Laboratory, USA) y com para resultados obtenidos por
ambos m étodos.
18
una m ayor precisión, sin necesidad de efectuar un m allado excesivamente
túpido cuando se desea tener en cuenta todos los detalles de la estructura.
Blaschak, [63], (1988), hace un estudio com parativo de las condiciones de
contorno absorbentes utilizadas en la sim ulación de ondas por diferencias
finitas. En particular las concernientes a la im plem entación num érica de
ecuaciones de onda unidireccional ("one-way w ave equations"). Investiga
la efectividad de las condiciones de segundo y tercer orden obtenidas de
ecuaciones de onda unidireccional m ediante diferentes aproximaciones. La
conclusión es que la aproxim ación de Padé es la que m ejor se comporta.
M oore, Blaschak, TaflQve y K riegsm ann hace una revisión sim ilar [64],
(1988).
19
C h u y C haudhuri [68], (1990), estudian discontinuidades en guía de
ondas com binando el algoritm o de Yee con análisis m odal. Esto perm ite
una presentación m odal (elegante) de resultados num éricos con DFDT.
20
c) Se trabaja en algoritm os de generación de m allado para efectuar el
m odelado, puesto que este en ciertos casos implica el manejo de decenas de
m illones de puntos.
21
1.6. Organización de esta tesis.
En esta tesis se abordan una serie de problem as en principio distintos,
pero todos inherentes al tratam iento m ediante Diferencias Finitas en el
Dominio del Tiempo de dispositivos pasivos de m icroondas.
22
d isp ersió n de varias d iscontinuidades que son contrastados con datos
procedentes de la bibliografía. En particular, se trata un codo recto en guía
rectangular, una unión T en el plano H, un filtro dieléctrico, y por último
un a guía de planos paralelos abierta al espacio libre. En este últim o caso se
trata u n problem a de radiación, pero que se aborda desde el punto de vista
circuital calculando el coeficiente de reflexión. Los cálculos de este capítulo
se llevan a cabo en simple precisión (reales*4).
23
24
Capítulo 2: Fundamentos
Matemáticos
2.1. Introducción.
El m éto d o num érico de D iferencias Finitas en el Dom inio del
Tiempo (DFDT ó FD-TD), cuya base m atem ático-com putacional se describe
en este capítulo, proporciona una solución directa de las ecuaciones de
Maxwell del rotacional. El algoritmo que se describe discretiza las ecuaciones
d el ro tacio n al em pleando diferencias d iv id id a s p a ra aproxim ar las
derivadas. El resultado final es un sistem a com pletam ente explícito de
m odo que no es necesario resolver u n sistema de ecuaciones y la necesidad
de m em oria y tiem po de cálculo es proporcional al tam año eléctrico del
m odelo.
25
d iferen cias fin itas p a ra la resolución de ecuaciones diferenciales en
derivadas parciales. El algoritm o de Yee o m étodo de las Diferencias Finitas
en el Dom inio del Tiempo (DFDT), se explica detalladam ente en el apartado
2.3.4, m o stran d o tam bién las restricciones de estabilidad y dispersión
im puestas p o r el tratam iento num érico. El tipo de excitación a utilizar y el
m odo de in tro d u cirla se analizan en el ap artad o 2.4; se contem pla la
separación entre zona de cam po total y zona de cam po dispersado y su
im p lem en tació n p a ra casos en dos dim ensiones. Las C ondiciones de
C ontorno A bsorbentes que deben de absorber las ondas para sim ular un
m edio indefinido, se tratan en el apartado 2.5. La disponibilidad de recursos
inform áticos es la principal lim itación al m étodo, p or ello en el últim o
ap artad o se abordan una serie de técnicas com putacionales que a la postre
perm itirán el ahorro tanto de m em oria como de tiem po de cálculo (CPU).
(2 .1)
26
siguiente sistem a de ecuaciones escalares es equivalente -en coordenadas
cartesianas (x, y, z)- a las ecuaciones del rotacional:
dE z 3Ex
(2.3.b)
dx dz
dEx dEy
3t 3y ' d x
3Hz 3Hy
= r (-5 n --5
' T --o E x ) (2.4.a)
3t e 3y dz
i
3Hy 3Hx
= r (-5
3t e 3x ^ L-’ “3^
3y “-aEz) (2.4.c)
27
caso TM
3Hx 1 3Ez
(2.5.b)
3t - ’|i. 3y
(25c)
3t “ ji. 3x u ^ c)
caso TE
28
m odo de calcular la derivada en uno de dichos puntos. La derivada se puede
calcular p o r la derecha, por la izquierda y haciendo uso de diferencias
centradas [9].
h2
F(x+h) = F(x) + h F'(x) + -j - F'(x) + ...
de donde
F(xj+1) - F(xí)
F'(xj) = ---- ------------ L + O(h)
que usando la notación x=xj, x-h=xj_i perm ite escribir la siguiente expresión
que nos da la derivada en xj con un error de orden uno en h.
F(xs) - Fíxj.!)
F(xj) = — 1 h ] -- + 0(h)
F ( x í ) - 2 F(xj_!) + F(x; 2)
F"(xj) = ---- 1------------ LÉ_ + Q(h)
p o d e m o s d e s p e ja r F'(x) :
F (xj+ i ) - F ( x j_ i )
F 'U j) = l+12h 1 + 0(h2)
30
F(xi+1)-2 F(xO + F (x ;.|)
F"(xj) = — ^ ^ ^ + 0(h2)
3u 3u c)2u
G(x, y , u , 0x , , 0X0X/ 0x9y / 3y9y )
c)2u c)2u 3u 3u
A S +2B + c H y ty + a + b + c u =
D=AC - B2
31
si D cam bia de signo en el correspondiente dom inio de definición, se dice
que es de tipo Mixto [72], [73].
Según H adam ard, [73], un problem a físico, descrito por una ecuación
d iferen cial con sus corresp o n d ien tes condiciones de contorno, es un
problem a bien planteado ("well posed") si su solución es única y depende
continuam ente de las condiciones de contorno.
32
Los errores que se cometen cuando se soluciona un sistema, ya sea
siguiendo un esquem a u otro, son de dos tipos:
33
la forma integral de las leyes de Faraday y Ampere al nivel de la celda
unidad.
/
v
A
Figura 2.1. Posición de las componentes del campo en la celda unidad de Yee
(DFDT) en tres dimensiones.
Cada punto del espacio se denota como:
F(iAx,jAy,kAz,nAt)= F^ijjk)
donde Ax, Ay, y Az son respectivam ente las dim ensiones del cubo unidad
usado en el m allado en las direcciones x, y, z; At es el paso tem poral
(incremento del tiempo) e i, j, k, n son enteros. Yee usaba diferencias finitas
centradas para las derivadas espaciales y tem porales, que son fácilmente
p rogram adas y tienen una precisión de segundo orden con respecto a los
increm entos espacial y temporal.
9Fn Fn (i+¿ , j , k ) - F M i - £ , j , k )
(2.8.a)
dx At
1 1
(2.8.b)
At
8 Ay Az ]
35
Ey(n+1>(i,j,k) = Ey<n)(i,j,k)+ (2.9.b)
AtrEz(n>(i,j+|,k)-Ez(n)(i,j-|,k) Ey(n)(i,j,k+|)-Ey(n)(i,j,k-|)i
|X l Ay ’ Az 1
Hy *n + 2 (i/j/k)=Hy^n'^(i,j,k)- (2.10.b)
At E x ^ G ^ ^ - E x ^ a j * - ^ Ez<n)(i+ |/j/k)-Ez(n)(i-|/j,k)i
|i Az ’ Ax
Hz(n+2)(i,j,k)=Hz(n'2)(i,j/k)- (2.lÓ.c)
36
El sistema en diferencias finitas (2.9) y (2.10), puede reducirse para los
casos bidimensionales TE y TM de (2.5) y (2.6).
— r ------------------------------------ . ---------------------------------------- i
e L Ax Ay J
1 1 At EzW(i,j+|)-EzW(i,j-|)
H x(n^ ( i j ) = Hx(n- 2 > ( iJ ) --[ -------------- ^ --------------- ] (2.1 l.b )
1 i At Ez(n)(i+lj)-Ez(n\i-4,j)
H y ^ (i,j)= H y (n-2 >(i,j)+- [ ------------^ ] (2.11.c)
At ^Ey(n)(i4,j)-Ey(n)(i-|j) Ex(n)(i,j+|)-Ex(n)(ij-|) ^
jo. Ax " Ay ^
37
H z (" ^ )(i j + i) - H z ( n^ ) ( i j - | )
Ex(n+1>(i,j) = E x^ -D O jH jt--------------- ^ --------------- ] (2.12.b)
Ey
Hy
Ez Ez
? -------- ----- 9 ------— 9
i
i
i
i Figura 2.3. Posición de las componen
i ----- ó Hx
Hx o ----- r del campo en la celda unid
i DFDT para u n modo TE .
i
i
i
Ez o ------ — ó — ------ó
Hy
38
Con el sistem a (2.9) y (2.10), el valor de cada nueva com ponente del campo
en cada p u n to depende únicam ente de su valor anterior y de los valores
anteriores de las otras componentes del campo en puntos adyacentes. De esa
form a, en u n tiempo dado, el cálculo de una com ponente del campo puede
calcularse en un punto a la vez o en p puntos a la vez (si se utilizan p
procesadores trabajando en paralelo). El algoritm o es como sigue:
- Se im pone la excitación.
En este paso a paso hay que asegurar la estabilidad del algoritmo para
que no se a c u m u le n los e rro re s, que se iría n in cre m e n ta n d o
39
p ro g resiv am ente d u ran te toda la sim ulación. Por ello la discretización
espacial y tem poral son elegidas según una relación, para acotar errores en el
proceso de m uestreo y evitar cualquier posible inestabilidad.
(2.13)
o definiendo:
se convierten en:
(2.14)
40
Para analizar la estabilidad del m étodo es conveniente descom poner
el algoritm o en diferencias finitas en dos problem as de valores propios, uno
para la variable espacial y otro para la variable tem poral. Se estudiará la
propagación de ondas planas propias del espacio numérico de datos, es decir,
de la ecuación (2.13) discretizada. Se determ inarán los m odos propios de la
diferenciación espacial y los modos propios de la diferenciación temporal. La
estabilidad num érica del algoritmo será posible si exigimos que el espectro
de valores propios espaciales este contenido en el dom inio de estabilidad de
la red, esto es de los valores propios tem porales, de esta forma evitaremos el
crecimiento sucesivo de los errores en cada iteración [22]. El objetivo de este
análisis es buscar una relación entre las discretizaciones tem poral y
espaciales que perm ita tal estabilidad.
^(n+W(n-^)
(2.16)
41
q= (2.17)
^ (n )
y se desea que I qi | <1 para todos los posibles m odos espaciales en la red
DFDT, para evitar tener cualquier m odo cuya dependencia tem poral sea
d ivergente (esto es, que crezca el valor del cam po sin lím ite). Esta es
esencialm ente la condición de estabilidad. S ustitu y en d o qi en (2.16) y
sacando factor común, obtenemos la siguiente ecuación para q¿,
cuya solución es
A,At ^I J k A t l s\ /n *r\\
qp ~ 2 ~± y i+ C y ) = o (2.19)
Real(A,) = 0 (2.20.a)
42
j ^ Inum érico A ^ (2.21)
< 2 -2 4 )
Real(A,) = O (2.25.a)
43
para todos los posibles valores de kX/ ky, y kz .
Esta ecuación nos da la restricción que debe verificar el paso tem poral
con respecto a los incrementos espaciales, que desnorm alizada es:
1
cAt < (2.27.a)
45
frecuencia para ilustrar los fenómenos sin realidad física asociados a esta
dispersión.
2 1 o1 o1 9 1
sin focoAt) .sin (^kxAx) sin (^kyAy) sin (^kzAz)
9 = \ 9 + 9 + o ) (2.29)
At Ax Ay Az
2 1 2 1 9 1 2 1
sin foxoAt) sin fckxAx) sin fokyAy) sin fckzAz)
S + H S ) (2.30.a)
(cAt) Ax Ay Az
9 1 9 1 9 1
sin CycoAt) .sin fokxAx) sin fokyAy)
S - = ( S + H ) (2.30.b)
(cAt) Ax Ay
46
En contraste con la relación de dispersión num érica, para el caso
trid im en sio n al, la relación de dispersión p a ra un a onda plana en un
continuo libre de pérdidas es:
A prim era vista presenta poco parecido la relación (2.30.a) con el caso
ideal (2.31). Se puede fácilmente dem ostrar que (2.30.a) se reduce a (2.31) en
el lím ite en que At, Ax, Ay y Az tiendan a cero, esto sugiere que la dispersión
n u m érica p u ed e reducirse al caso ideal sim plem ente con u tilizar un
m allado m ás fino.
Ah cosa
(2.34.a)
47
Ah sin a
B= 2----- (2.34.b)
vp 2%
c
= Tr
kfinal
(2.35)
siendo kfjnai el valor final de k obtenido con (2.33) (usualm ente partiendo de
ko=27t bastan dos o tres iteraciones).
48
celdas, dependiendo del ángulo de propagación. Como resultado de esto, el
m odelado m ediante DFDT de pulsos de corta duración y como consecuencia
de ban d a infinita, da lugar a una progresiva distorsión de estos; las altas
frecuencias espaciales, se propagan con m ayor lentitud que las bajas, y las
que tien en m u y altas frecuencias espaciales con longitudes de onda
inferiores a dos o tres celdas no se propagan. La dispersión num érica hace
que los pulsos de duración finita se ensanchen y dejen un residuo numérico
de alta frecuencia espacial en los bordes del pulso. Para tener un control de
la dispersión es conveniente efectuar la transform ada de Fourier espacial del
p u lso d esead o , entonces o la v ista del espectro, se seleccionan las
dim ensiones de la celda, de m odo que las com ponentes principales del
espectro sean resueltas con al m enos diez celdas por longitud de onda. Esta
m ed id a lim itaría la dispersión a m enos de un 1% a pesar de la debida al
ángulo de propagación con la red.
0.99
0.98
U
0.97
>
0.96
0.95
0.94
0* 10 * 20* 30’ 50' 60* 70' 80’ 90'
50
Ideal
1. 0
0. 8
0.6
0.4
0.2
0. 0
0 .IX
51
2.4. Excitación.
Los parám etros que caracterizan un sistem a, se obtienen con DFDT
m ediante dos posibles vías:
52
luego generan sobreoscilaciones artificiales en la señal que se propaga
("overshoot" y "undershoot"), [26].
Trabajos recientes con DFDT, [55], han m ostrado que el num ero de
ciclos que requiere la onda incidente para alcanzar el estado estacionario
depende inversam ente de las pérdidas de la estructura o fenómeno que se
m odela.
53
Zona 1:
Campo total
Estructura
Zona 2:
Campo reflejado
Interfase zona 1-
zona 2
Zona 1: H
Campo total X
- y
E
t o t
•¡O H
y ►
Zona 2: -+■ —
Campo reflejado HS
x
(b )
Figura 2.6. División de la red DFDT en zona de campo total y zona de campo
reflejado: (a) División de zonas, (b) Componentes del campo en
el plano de conexión.
54
La conexión entre ambas zonas, debe asegurar la consistencia de las
derivadas espaciales a través de la interfase de separación.
55
bastante inferior a la de la onda incidente (esto se supone, pero depende del
problem a). Puesto que el cam po incidente sobre los bordes del m allado con
las CCA sería pequeño, bastante mas pequeño sería el residuo de campo no
absorbido, con lo cual se m ejorarían los resultados.
56
podem os m odificar el algoritmo para conseguir consistencia en los cálculos
cuando estos se efectúen en la interfase de separación entre las dos zonas.
Para ello supongam os que la onda avanza en la dirección "y". Consideremos
ahora el pu n to (i,j0-l/2 ), en la zona de separación:
c H x n+1/2inc(i/j0- l / 2) (2.37)
siendo C = At/(sAh)
57
Campo incidente Puntos de tratamiento especial para Ez
(a) H y (b) H x
\ \ £ ¡i A A a j =i í +1
\ A C^) — t— <•>--+ —O— <— -0-- -<•) A j = jj
1
\ : . . . I -
+ • ♦ <í> A
i
ii - ZONA 1 xi Interfase zona 1- zona 2
A ©1 + * ♦ * ♦ k • + © A
1
• 1
t t
1 1
A -+— Qj~ —f - -© Ai j = j0
A A A A A A _ 0 - j0 - i
r— A ' H
_° _° ZONA 2
II H II M
(a)
Puntos de tratamiento especial para
Campo incidente Ez (a) H x (b) H y
e © © © © /7 i = ¿i + l
® Á ♦ • ♦ * ♦, • ♦ • ♦ Á ©
i i
i i
T t
© A- - A - ■*- -A- - 1- - A - - + - - A - - f - - A ©
j =jo
^ e © © © © © ^ j - J0 - i
ZONA 2
II H II II
(b)
58
Las com ponentes Ez situadas en las otras interfases de la figura 2.7.a
son tratadas análogam ente.
O bservando las figuras 2.6.b y 2.7.b, vemos que aparece otro problema
de consistencia para las com ponentes del campo reflejado H x ref y Hy ref ,
que están situ ad as justo m edio paso espacial fuera de la interfase zo n al-
zona2. Para calcular de form a adecuada el campo m agnético reflejado, es
necesario conocer el campo eléctrico reflejado, Ez ref, medio paso espacial a
am bos lados de la com ponente del campo magnético. Pero justo al lado de la
posición del cam po m agnético reflejado, está el campo eléctrico total situado
sobre la zona de separación. Esta inconsistencia, se resuelve de un m odo
análogo a la del caso anterior. Sabemos que en la interfase de separación el
cam po total es sum a del incidente más el reflejado
entonces se m odifican (2.1 l.b) y (2.1 l.c) para calcular H x ref y Hy ref justo un
p u n to fuera de la interfase. Para H x ref queda:
59
H x n+1/2refGJo-l/2) = H x n+1/2ref(i,jo-l/2) + D ( E z nref(i,]0-1) -Ez ntotal(i,j0) )
+
D Ez ninc(Mo) (2.40)
siendo D = At/(|iAh) }
)=)0- l /2 , i= i0 ... 4
Hyn+1/ 2 (i0 - l / 2 ,j) = Hyn+1/ 2(i0 - l / 2 ,j) | yee - D Ezninc(i0 /j) (2.41.c)
60
del cam po incidente para Hx y Hy, en puntos del m allado m edio paso
espacial fuera de la interfase (indicados por triángulos en la figura 2.7.a)
61
z
i nc
i nc
i nc
Plano de fase
constante
k
1
k
o
J = J o
62
oí oio O's
*©■-
°V>
♦ ♦ + ♦ + + t t» ♦ + + + + + t
♦ + + t + + ♦ i
+ + + + + + + f ♦ + + ♦ + + t
+ + + + + ♦ t i t + + + ♦ + !
t + + t ♦ + _ ♦ (ó f ♦ + + + + ♦ + t
©--- ----------------------------- (i). ■*■■.
°o2 * * °03
°1 °4
(a)
\ ©•- °3 o; ... °3
♦ ♦ ♦ + ♦ + + t j + + + + + ♦ t
+ * + + + + + i ♦ + t t + + +
+ * + + ♦ + + \ i + * + ♦ + + +
t + + I ♦ ♦ ♦ t i ♦ + •+ 1 + ♦ +
+ + + ♦ + ♦ ♦ j i + + * + + + +
^ . (7) J7v-►. . _____________ é .
o® - °2
O 'i? 2
(b)
°ío oí °2&-.. °3
-1©-
+ ♦ + * ♦ f t t + + ♦ ♦ ♦ {
♦ + + ♦ + + + + ♦ + + + t i
+ + + + + + \ + + + + ♦ + +
+ + + t + ♦ t t + + + + + t
+ + + + + + + + i + + + + + t
©- ■* á).+ -
°1 °4 °2 °3
(c)
63
- cara de delante:
i= i0 ... i i , j= j0 , k = k 0+ l / 2 ... k r l / 2
- cara de detras:
i= i0 ... i i , j = j i , k = k 0+ l / 2 ... k j - 1 / 2
- cara de abajo:
i= i0 ... h , j=j0+ l / 2 . . . j i - 1 / 2 , k = k 0
- cara de arriba:
64
Eyn+1(i, j, ki) = Eyn+1(i, j, k!) | yee + C H xn+1/2inc(i, j, lq+1/2) (2.45.b)
- cara izquierda:
- cara derecha:
Hxn+1/ 2 (i, jo-1 / 2 , k) = Hxn+1/ 2(i, jo-1 / 2 , k) | yee+D Ezninc(b j0, k) (2.48.a)
Hzn+1/ 2(i, jo-1/2, k) = Hzn+1 / 2(i, jo-1/2, k) | yee-D Exninc(i, j0, k) (2.48.b)
65
- cara de detras (exterior):
Hxn+1/ 2 (i, ), ko-1/2) = Hxn+1/ 2 (i, ), k0- l/2 ) I yee-D Eyninca j, ko) (2.50.a)
66
i=i0 - l / 2 , j=j0 + l / 2 ... ji- 1 / 2 , k=k0...ki
67
reducir suficientem ente como para perm itir que los resultados num éricos
m ediante DFDT sean fiables durante todas las iteraciones tem porales de la
sim ulación.
Figura 2.11. Simulación num érica de la incidencia de una onda plana sobre
el borde izquierdo de un dom inio de cálculo bidimensional.
68
En este ap artado se va a tratar la im plem entación num érica de
ecuaciones de onda de un solo camino, para utilizarlas como CCA en la
resolución m ediante DFDT de las ecuaciones de Maxwell. Las ecuaciones
que se consideran se pueden derivar por dos caminos diferentes [64]:
1
UxX + Uyy - ^2 Utt - 0 (2.54)
(2.55.b)
de m odo que la ecuación de onda (2.54) se puede escribir de una forma mas
compacta:
LU = 0 (2.56)
69
L U = L+ L” U = O (2.57.a)
siendo L” definido
i Y
L = Dx - ~ "y1 - S (2.57.b)
con
CDy
S= ^ ( 2 .5 7 .c)
L" U = 0 (2.58)
aplicado en x=0, produce una CCA exacta que absorbe cualquier onda
procedente del interior de fí. El operador L+ efectúa la misma función en un
contorno a la derecha del dom inio de cálculo, absorbiendo cualquier onda
procedente de Q. que se propague hacia la derecha.
70
que se haga, la reflexión puede m inim izarse sobre un am plio rango de
ángulos de incidencia. Las CCA que utiliza M ur [27], son aproximaciones de
Taylor de segundo orden de (2.57.b). Estas se van a deducir a continuación.
(2.59.a)
(2.59.b)
(2.60.a)
1 c
U x t+ ^ Utt - 2 u yy = 0 x=H (2.60.b)
y=0 (2.60.c)
1 c
U y t+ c U tt' 2 U xx - 0 y=H (2.60.d)
71
Estas CCA se aplican en DFDT sobre las componentes del campo ? Ó H que
descansan sobre los bordes o son tangentes a ellos.
1
Uxx U yy + U z z ~ =0 (2.61.a)
Se construye el operador L:
^
L= 2 i n x n 2 i n
+ Dy^2 + Dz^ " ^ 2 .2 (2.61.b)
que se puede factorizar como en (2.57.a) para dar el operador L , que es una
CCA exacta con la m ism a form a analítica que (2.57.b) pero con otra
expresión para S :
cD v 9 9
S = A / CQp) + ( ^ " ) (2.61.c)
Dt cD y 2 cD z 2
(D x - T * r 5 r + W ) u = 0 e ' 6 2 )
72
l e e
U x t _ ^ U tt + 2 U y y + 2 U z z ~ ^ (2.63.a)
A nálogam ente se d erivan las CCA para los otros bordes del m allado
tridim en sio n al:
Uy t + ^ U tt - f u x x - | u zz = 0 y =H (2.63.d)
Las ecuaciones (2.63) se aplican a las componentes del campo situadas en los
bordes del mallado.
73
. . Pmfe)
r(s ) = ^ 5 j- (2.64)
V 1 - s2 = p 0 + p 2 S2 (2.65)
P 0 + P 2 S2
r(s) = --------- 2 (2.67)
qo+q2s
que da la siguiente CCA de tercer orden,
2 P0
% U x tt + ^ q 2U x y y " ^ ^ t t t " c P 2^ t y y - 0- (2.68)
74
S=0, por esa razón las correspondientes CCA absorben perfectam ente una
onda que incide norm alm ente al contorno donde están definidas.
H igdon [77], [78], presenta una form ulación m as general, que perm ite
el diseño de CCA para ondas planas con un ángulo de incidencia cualquiera.
Para ello desarrolla la siguiente condición de contorno,
f n ^
n
j= 0
V;
(cosaj ^ - c ^ ) |LJ = 0
J
(2.69)
eos a i eos a 2 + 1
Po ” eos a i +cos a 2 (2.70.a)
-1
p
VL
2 = -------------------
eos a i +cos a 2 (2.70.b)
75
caso B: factorización de la relación de dispersión
^ 2 = kx2 + ky 2 = k2 (2.72)
—> c kx c ky
V = (vx , Vy) = ( , - j p ) = (c C0S<|), c sen<j>) (2.73)
siendo <|) el ángulo en el sentido del reloj que form a el vector 1 ? , con 0 < <|> <
2 tu . La ecuación (2.54) perm ite entonces como solución una onda plana
propagándose en todas direcciones.
76
1 - (2.74)
1 ry ky
kx = - k (1 S2) = - k (1 - 2) (2.78)
2 2 k
/
o
kx k= - k2 + ~2~. (2.79)
Para una onda plana como la definida en (2.71) se deducen las siguientes
relaciones:
1 (-jkx )(jco) (0
7 U xt = _ U = kx — = kx k U (2.80.a)
77
(-jky)(-jky) ky
1 1 1
cuxt= ^ utt - 2 uyy (2 .81)
(2 .82)
78
Las derivadas parciales en (2.66) y (2.68) son aproxim adas en los
puntos del contorno del dominio de cálculo. Para calcular las derivadas
norm ales al contorno con diferencias centradas es necesario conocer los
valores del campo en puntos exteriores al mallado. Esto último se evita
expandiendo las ecuaciones en diferencias para puntos ficticios que se sitúan
a medio camino entre el borde del mallado y la primera columna de puntos
del interior. La situación de estos puntos se muestra en la figura 2.12,
U ( o , j +D 1 J ( U +1)
o () y *= ( j + 1 ) Ay
— y = JAy
U (o , j ) U Ü J) U dd)
o.
w S — — y = (j-l)A y
U(o.j-i) U d .j - D
Figura 2.12. Puntos del mallado cerca del borde x=0, que son utilizados en la
derivación de las CCA de Mur.
79
para el caso de que el borde este en x=0. Para el punto (0, j), las derivadas
p ara establecer las CCA se efectúan en el punto (1/2, j). El campo en el punto
(1 / 2 , 0 ) se pone como prom edio del campo en los puntos vecinos,
1 U ( 0 ,j) + U (Í ,j)
U<2 , j) = ----- 2------ (2'83)
aun+1 aun_1
— a / 2' ,
Uxtl ( 1 / 2 , j )~ 2At (2-84a)
La derivada tem poral en (1/2, j) se pone como prom edio de las derivadas
tem porales en (0 , j) y (1 , j):
a2un a2un
ax at2 i 0 / ] ) + at2 (1' ])
Ut tl(l/2 ,j)= 2 (2-84‘b)
i ( u n+ 1( o , j) - 2 u n ( o , j) + u ^ c o , ¡) u n+1q , p - 2 u n g , ¡) + u n-1( i ,
=2l At2 + At2
+ á f e ( un <°'i> + U"»'»)
2
En x=0.
- En x = H:
2
(cAt) Ax / n n „ „
+ nA 2/ „ . x ( U (H')+1) "2U (H,j) + U (H,j-1)+U (H -l,j+l)
2Ay ^(cAt+Ax)
81
En y = 0:
2Ax
+
cAt+A x(Un (i' 0) + U n « '« )
2
-2 Un (i,l)+ U n(i-l,l))
En y = H: _
2Ax
cAt+Ax
( u n (i,H) + Un (i/H - 1 ))
2
82
sólo ángulo de incidencia (n= l en la fórm ula (2.69)). Supongam os la
siguiente aproxim ación:
Ux - - U t = 0 (2.87)
que son las CCA de M ur de prim er orden [27] para una onda que se propaga
en dirección -x. Para poder utilizar diferencias centradas, la ecuación (2.87) es
aplicada a un punto ficticio que dista m edio paso espacial de la frontera y
esta a m edio paso tem poral de dos instantes consecutivos de tiem po (ver
figura 2 .1 2 ).
n+1/2.
Para calcular U v (1/2, j), calcula el prom edio de esta derivada en t=nAt y
x
(n+l)At:
A .n + 1
U“ ' ~(1 , j)-un+1(o, j) una,j)-un(o, j/
Ax + Ax
n+1/2
Para calcular U . (1/2, j), calcula el prom edio de esta derivada en i=0, y
i=l:
n+1 / 2
/ - - n + 1 /2
au aun + \ 1 /2
. A .n + 1
u“^ a j)-unq,j) ull+^(o/ j) - un(Q,j/
At + At
83
- En x = 0:
- En y = 0 : ¡:
84
Considerem os la onda plana que se aleja del dom inio de cálculo de la
figura 2.11. La expresión de esta es:
UCx^t) = exp (j( cot - ~r) = exp (j( cot + k x cos0- kysen0 )) (2.90)
El campo total en el borde del m allado, deberá satisfacer las CCA definidas
en este. E ntonces, considerando la existencia de una onda reflejada,
generada en el contorno, el campo total en x = 0 tendrá la forma,
COS0 - 1
R= coieTT <294>
En el trabajo de Trefethen y H alpern, [76], se desarrollan siete técnicas de
aproximación (siete familias de coeficientes p y q). Estas técnicas ya han sido
85
com entadas anteriorm ente en la form ulación generalizada de obtención de
CCA (apartado 2.5.1). Las tablas I y II m uestran los ángulos de incidencia para
los cuales cada CCA absorbe perfectamente la onda incidente (R=0).
qg = 1 .0 0 0 0 0 en todas
qQ = 1 .0 0 0 0 0 en todas
Padé 1 .0 0 0 0 0 -0.50000 0 .0 0
p u n to s-N e w m an 1 .0 0 0 0 0 -1 . 0 0 0 0 0 0.0,90.0
86
Las figuras 2.13 y 2.14 m uestran el comportamiento de algunas CCA
(las mejores) en el intervalo de ángulos [0, n/2].
- tipo (2,2)
— tipo (2,0)
0 10 20 30 40 50 60 70 80 90
Angulo de incidencia
(a)
- tipo (2,2)
■*" tipo (2,0)
o 10 20 30 40 50 60 70 80 90
Angulo de incidencia
(b)
87
10°
10'1
Coeficiente de reflexión
— Padé
Cheby en sub.
10
7
10
9
0 10 20 30 40 50 60 70 80 90
Angulo de incidencia
88
2.6. Simplificaciones Computacionales.
El desarrollo del m étodo DFDT ha ido m uy ligado a la evolución de
los g ran d es ordenadores. El m odelado ,de problem as electrom agnéticos
com plicados im plica la utilización de cantidades m asivas de m em oria que
deben ser accedidas con rapidez en la ejecución de las operaciones para que
el algoritm o sea atractivo. La sim plicidad de las operaciones que se ejecutan
y el hecho de que el grueso de los cálculos resida en bucles hace que el
algoritm o de Yee sea altam ente vectorizable. Los ordenadores que poseen
procesadores vectoriales perm iten efectuar las operaciones entre vectores
calculando en bloques en lugar de calcular de elem ento en elem ento
(procesadores escalares). La vectorización (com pilación con el procesador
vectorial) es efectiva cuando da lugar a un program a ejecutable mucho mas
rápido que el obtenido con el tradicional com pilado escalar. Esto sólo es
posible cuando el num ero de iteraciones en los bucles vectorizados es el
óptim o ó próxim o a este, y se debe a que coincide con el num ero de
elem entos de los bloques que calcula el procesador vectorial. Por otro lado,
p ara conseguir efectividad es necesario que las m atrices involucradas sean
grandes. A unque el procesado vectorial m aneja bloques de datos, es mas
lento que el procesado escalar, y para conseguir rentabilidad del prim ero
respecto del segundo es conveniente trabajar con m uchos "bloques". La
vectorización del código inform ático le da a la técnica DFDT una gran
rap id ez a la hora de obtener visualización de fenóm enos y datos de diseño
p e ro las consideraciones anteriores indican que la efectividad de la
vecto rizació n no es absoluta, y dep en d e en gran p arte del problem a
a n a liz a d o y de sus co rresp o n d ien tes condiciones de contorno. Las
operaciones se repiten en el interior de las zonas m odeladas siguiendo el
esquem a de Yee (apartado 2.3.4) hasta que se llegan a los contornos, donde
las condiciones de contorno "cortan" la vectorización.
89
cada uno, o a cada grupo de procesadores les p uede corresponder una
porción del dom inio total. Pero la elaboración de un program a para trabajar
sobre una m áquina paralela es todavía un trabajo complejo y delicado ["2 Q
S ym p o siu m on H ig h P erform ance C om puting", M o ntpellier, 7-9 de
Octubre, 1991].
90
E7 (O, 1 , 1 / 2 ) E
bX (11 / 2 . 1 ''0) ( 1 / 2 -1 - 1 / 2 )
___________ E (1/2 ,1 , 1 )
A
E y (O, 1/2.0) ^
(0,1/2,1/2
, 1/ 2 , 0 ) H7 (1/2,1/2,1)
Ev (0 ,1/2 , 1 )
Ey (1/2 ,0 ,1 )
(0 ,0 , 0 ) Ez (0 ,0 ,1/2 ) Ex ( 1 / 2 , 0 , 0 ) Hy (1/2,0,1/2)
E-, (1 ,1 ,1/2 )
Y
X
1/ 2 ,0 ) Hy (1 ,1 / 2 , 1/2
Ey (1,1/2,1)
-P~Z
(1 , 0 , 0 ) E7 ( 1 , 0 , 1 / 2 )
(a)
E7 ( 0 , 1 , 0 ) E (0,1,0) Hy ( 0 , 1 , 0 )
A Ex (0,1.1)
0 , 0 , 1)
Ey (0 ,0 ,1 )
Ez (0,0,0) Ex (0,0,0) Hy (0,0,0)
E2 (1 ,1 ,0 )
Y
X
y (1,0,0) Hy>( 1, 0, 0)
------ Ey(1 ,0 ,1 )
-+Z
Ez (1,0,0)
(b)
91
El resu ltad o de este criterio se m uestra en la figura 2.15.b, y las
ecuaciones (2.9) y (2.10) quedan de la siguiente forma:
A, Hy(n+2)(i,j,k)-Hy(n+2)(i-l,j,k) Hx<"+¡>(ij,k)-Hx<n+¡>(i,j-l,k)
— r ----------------------------------------- 1
£L Ax Ay J
, U , U AtrEz(n)(i/i+ l/k)-Ez(n)(i,i/k)
Hx(n+¡)(i,j,k) = Hx(n-;>Ü,j,k)-— [---------------------- 1— (2.96.a)
Ey(n)(i,j,k+1)-Ey(n)(i,j,k)
Az ]
, \ \ At Ex(n)(i,j,k+1)-Ex(n)(i,j,k)
Hy(n+2)(i,j,k)=Hy(n_J (i,j,k)^ ( 2 . 9 6 . b )
Ez(n)(i+l,j/k)-Ez(n)(i/j<k)
Ax
Ex<n)(i,j+1 ,k)-Ey(n>(i,j,k)
“ ........ Ay -I
La m em oria que se consum e es del orden de 6 N, siendo N el núm ero de
celdas que se utilizan en el modelado.
92
Para im plem entar el sistema (2.95) y (2.96) cuando se hace uso de una
celda unid ad cuadrada (Ax=Ay=Az=Ah) , es conveniente definir y almacenar
los siguientes coeficientes:
D = ¡s¡;
c - ñ , v -v »
Hy(n+^(i/j,k)=Hy(n-^(i/j,k) (2.99.b)
Hz(n+¡>(i/j/k)=Hz(n-J(i,i/k) (2.99.C)
93
Reducción a casos bidim ensionales
+ D [Ey(n)(i,j)-Ey(n)(i-l,j) - Ex(n)(i,i)-Ex(n)i,j-1) ]
= D if (2.102.a)
94
Hy (i,j+1)
Ez(¡,j+1) p ---------------- 9 ----- o Ez (i+1 ,j+1)
Hx (¡,j) (V
r
t y 0, + i )
E y (¡,i)
95
Introduciendo esta nueva variable en las ecuaciones (2.98) a (2.100) se ahorra
al ordenador efectuar una m ultiplicación en cada punto, y en cada iteración.
Sin em bargo para el caso TE, es conveniente el cambio:
^ = D ÍÍ (2 .1 0 2 .b)
H x ^ ^ i^ k ) = Da Hx(n'^G/j,k) (2.104.a)
96
H z(n+2)(i/j/k) = Da Hz(n'^(i,j,k) (2.104.c)
H z(n+y(i,j) = Da H z ^ ü ,] ) (2.106.a)
+ Db [Ey(n)(i,j)-Ey(n)(i-l,j) - Ex<n)(i,j)-Ex(n)(i,j-1) ]
OAt
-
1
28
Ca = ----- — 77 (2.107.a)
o At
1 +
97
At 1
^ ” e Ah a At (2.107.b)
1 + ^ T
p'At
(2.107.C)
At 1
(2.107.d)
ji Ah p' At
1 +
98
Esto p u ed e hacerse definiendo vectores de constantes, con tantos
elem entos como zonas,
Ca = Ca (m) (2.108.a)
Cb = Cb (m) (2.108.a)
Da = Da (m)_____ __ (2.108.a)
en las interfases entre dos zonas se utiliza un valor especial para estas que se
obtiene de las ecuaciones (2.107),haciendo uso del valor m edio de los
parám etros eléctricos. Por ejemplo en los puntos de la interfase entre dos
m edios de constantes el, e2 , jllI , |i 2 se aplica:
el+e 2
8 = 2 (2.109.a)
ul+u 2
= (2.109.b)
99
100
Capítulo 3: Modos TE y TM en guías
de onda. Análisis con DFDT
3.1. Introducción.
Las líneas de transmisión y las guias de ondas son sistemas guiadores
que se utilizan para transportar energía a la frecuencia de m icroondas [79],
[16]. Las características deseables para estos sistemas son principalm ente el
que p erm itan la transm isión de un sólo m odo con una baja atenuación
sobre u n am plio rango de frecuencias. La estructura es llam ada línea de
tra n sm isió n cu an d o el m odo d o m in a n te es u n m o d o T ran v ersal
ElectroM agnético (TEM), y se denom ina guía de ondas cuando el m odo
dom in an te es un Transversal Eléctrico (TE) ó un T ransversal M agnético
(TM). Es im portante esta distinción dado que hay una diferencia conceptual
y m atem ática en el tratam iento de ambos casos. En el caso de una línea de
transm isión la frecuencia de corte es nula para el m odo fundam ental, y la
estructura del campo se obtiene resolviendo la ecuación de Laplace [80], que
es u n a ecuación de tipo elíptico. Los m odos TE y TM si p resen tan
frecuencias de corte y su estudio puede derivarse por dos vías alternativas.
U na posibilidad es la de resolver la ecuación de H elm holtz en la sección
transversal de la guía para las com ponentes longitudinales del campo. La
ecuación que se resuelve es una ecuación de valores propios y es de tipo
elíptico tam bién. Los valores propios están relacionados con las frecuencias
de corte y los vectores propios perm iten obtener la estructura de campo en la
sección transversal de lá guía de los respectivos m odos. Otra posibilidad de
análisis equivalente que se introduce en este capítulo es la de considerar la
sección transversal de la guía como u n resonador bidim ensional. Supuesta
u n a dependencia arm ónica con la variable espacial en la dirección de
propagación (simetría traslacional), sólo hay que resolver la estructura del
cam po en la sección transversal de la guía. H asta aquí el planteam iento es
idéntico al análisis convencional [79],[80], la diferencia está en que en lugar
de resolver la ecuación de H elm holtz bidim ensional, se resuelven las
101
ecuaciones de Maxwell del rotacional (2.5), (2.6), en la sección transversal de
la guía. Previa im posición de una excitación, la evolución tem poral del
cam po está determ inada por las ecuaciones (2.5) y (2.6). Las condiciones de
contorno en las paredes de la guía solamente perm itirán unas determ inadas
estructuras de campo que serán propias de la sección transversal de la guía.
Estas son los m odos y las frecuencias de estos se corresponderán con las
frecuencias de corte de estos.
102
con claridad su utilidad. N avarro, N uñez y M artin hacen una exposición
breve y concisa de su algoritm o, de la excitación que eligen y del uso que
hacen de la T ransform ada Rápida de Fourier (FFT), para determ inar las
frecuencias de resonancia y la geometría de los modos.
(3.1)
descom poniendo los cam pos en sus com ponentes axiales y transversales, y
suponiendo una dependencia en z del tipo exp(-jPz) (onda que viaja hacia la
izquierda) [79] tendremos:
(a) M odo Transversal eléctrico (TE) o m odo H. Las soluciones tienen Ez=0
pero H z* 0. Todas las com ponentes del cam po se p ueden derivar de la
103
com ponente axial Hz del campo m agnético. La ecuación a resolver con las
c o rresp o n d ien tes condiciones de contorno es la ecuación escalar de
H elm holtz para las coordenadas transversales,
2 9 9
siendo V t = — 2 + — 2
dx dy
{Vt + c o V - p 2} Ez =0 (3.2.b)
104
opuesto al de propagación en el espacio libre). Estas ecuaciones, junto con las
condiciones de contorno del sistem a p e rm itirá obtener la evolución
tem poral de las componentes del campo a partir de su distribución inicial en
el origen de tiempo (t=0 ).
105
ezm(x,y/t)= gm(x,y) exp(j2jcfm t) (3.3)
n = 0 ,l,2 ,-, N -l
106
N-l
Fk(i,j)= X e z (tj) exp(-jíúknAt) (3.6)
n=0
siendo cok=27ik/NAt,
\ 1 l-exp(i27cfmNAt)
Fk(i'Í> = 2 j Cmgm(X,y) 1-exp j27r(fr^fky (3'7)
m
siendo fk=k/N A t,
- í cm gm (x,y)
pk M = ^ , 2kAÍ
X
r=o'
firr(fk+At)
(3.8)
r r+1
que es periódica de periodo N, puesto que lk+N+At - ík+ At
V Í Cm gm (x,y) 1
FkW = « ( 3 -9 )
para k=0,l,2...,N/2.
107
Calculo de las frecuencias de corte
Af = 1/NAt (3.10.a)
fm s
fkl Fkl-fk2 Fk2
=— 5-------- (3.10.b)
m Fkl'Fk2
D eterm inación de las líneas de campo
gm(x,y) Fk(x,y)
gm(x0,y0)
E sta relación nos perm ite obtener las líneas de cam po en la sección
transversal de la guía.
108
3.3.2. Numeración de los modos.
El algoritmo DFDT perm ite la obtención de las frecuencias de corte en
o rd en creciente y la distribución de campo en la sección transversal de la
guía p a ra cada m odo, pero no establece ninguna form a de num erar los
m o d o s. La n um eración se p o d ría hacer a po sterio ri, observando las
estructuras del campo y su distribución en torno a ciertas líneas de simetría
d e la guía. De todos m odos no se va a entrar en esta discusión. En los casos
que se tratan, se compara con la bibliografía existente utilizando su misma
notación. En cualquier otro caso sim plem ente se o rdenaran en orden
creciente sin efectuar clasificaciones.
109
3.4.1. Análisis de la excitación. Guía rectangular
con sección b=2a.
El m odelado de la sección transversal se efectúa con u n m allado de
40x20 celdas cuadradas. Las dim ensiones de la celda que se utiliza en el
m allado de la guía son Ax=Ay=Ah=0.5.10‘^m y el paso tem poral -siguiendo
la condición de estabilidad de C ourant [22] - se elige At=Ah/2c , siendo c la
velocidad de la luz en el vacio. Se sim ula la evolución del cam po duran te
8000 pasos tem porales sobre los cuales se efectúa la TDF, esta elección es
heurística y da lugar a una precisión en la TDF de 0.149896 GHz. Los
resultados en función del tipo de excitación se analizan a continuación.
n<50
n>=50 (3.12)
110
f(GHz)
Figura 3.1. Espectro en los puntos (20,10) y (15,15) del mallado.
TABLA 3.1. Frecuencias de corte de los cuatro prim eros m odos TE en guía
rectan g u lar hueca. D esviación respecto de los resu ltad o s
analíticos. Dimensiones de la guía a/b= 2, a=2cm.
Figura 3.2 Líneas de campo m agnético constante para los modos: (a)TEio;
(b)TE20; (c)TEn; (d)TE2i.
énAt-rioAtA
A0 exp- A
n<2n o
ez(i/j)=
lo. 0 n>=2n o
(3 . 13)
112
a) Esta excitación tem poral es distribuida a lo largo de la línea i=20 desde j=5
a j=15. El espectro se m uestra en línea "a" de la figura 3.3. Las líneas de
Hz=constante para cada m odo se m uestran en las figuras 3.4.a-d.
b) Para com parar resultados se ensaya tam bién la excitación en la línea j=10
desde i=10 a i=30. El espectro se m uestra en línea "b" de la figura 3.3. Las
líneas de Hz=constante para cada m odo se m uestran en las figuras 3.5.a-d
113
5 10 15 20 25 30 35 40 45
f(GHz)
Figura 3.3. Espectro para cuatro tipos de excitación en el mismo punto.
114
(a) (b)
(c) (d)
Figura 3.4. Líneas de campo magnético constante para los modos con
excitación a): (a)TEio; (b)TE2o; (c)TEn; (d)TE2 i-
(a) (b)
(c) (d)
Figura 3.5. Líneas de campo magnético constante para los modos con
excitación b): (a)TEio; (b)TE2o; (c)TEn; (d)TE2 i.
(a) (b)
(c) (d)
Figura 3.6. Líneas de campo magnético constante para los modos con
excitación c): (a)TEio; (b)TE2o; (c)TEn; (d)TE2i-
116
f(GHz)
Figura 3.7.
(c) (d)
Figura 3.8. Líneas de campo magnético constante, modos : (a)TEi0; (b)TE2o;
(c)TEn; (d)TE2i.
117
TABLA 3.II. Frecuencias de corte de los m odos TE en guía rectangular hueca.
D esviación respecto de los resultados analíticos y tiem po de
CPU consum ido en el cálculo del espectro. Dim ensiones de la
! guía a /b = 2 , a= 2 cm.
118
c u a d ra d a no se id en tifican apenas con las líneas que se. obtienen
analíticam ente (figuras 3.2.a-d) [87]. La situación de las "sondas" crea fuentes
de cam po estac io n a ria s que d isto rsio n a n la "im agen de Fourier"-
transform ada de Fourier del campo- para las frecuencias de corte.
119
b) Su situación espacial no debe de favorecer a ningún m odo en particular,
por lo que debe de estar situada fuera de líneas de simetría.
c) La determ inación del espectro se efectuará tam bién sobre puntos fuera de
líneas de sim etría, para evitar posibles nulos del campo.
e) C uando el pico del espectro sea dudoso (poco acentuado). Com probar la
validez de las frecuencias de resonancia m ediante la observación de la
geom etría. Esto es la TDF para esa frecuencia sobre toda la sección
transversal de la guía.
f) R epetir el proceso anterior excitando en otro punto para com probar las
frecuencias y las geometrías obtenidas con la excitación anterior. En caso de
no co in cid ir los resu ltad o s, esto p erm ite d escu b rir otros m odos no
generados con la prim era prueba
120
3.4.2. Análisis de la convergencia. Guía rectangular
i
con sección b=2a.
i
121
TABLA 3 .III. Frecuencias de corte de los m odos TE en guía rectangular
hueca obtenidos con DFDT. D esviación resp ecto de los
resultados analíticos para varios m allados. Dim ensiones de la
guía a /b = 2 , a= 2 cm.
122
Influencia del tiempo de simulación
123
TABLA 3.IV. Frecuencias de corte de los m odos TE en guía rectangular
hueca o btenidos con DFDT. D esviación respecto de los
resu lta d o s analíticos para distinto núm ero de iteraciones
tem porales (NAt). D im ensiones de la guía a /b = 2 , a=2cm.
M allado 20x10.
500 - - - -
124
Una nota im portante a tener en cuenta es que el hecho de aum entar
el núm ero de pasos tem porales aum enta la resolución de la TDF en la
determ inación de las frecuencias de resonancia num éricas. Estas frecuencias
de resonancia num éricas van asociadas a la estructura de campo resonante
num érico que se genera en la simulación, por lo tanto una m ayor tiempo de
sim ulación no supondrá un m ayor acercam iento a las frecuencias reales
sino a las num éricas. Este análisis sólo nos perm ite com probar a partir de
que num éro de iteraciones tem porales ya no se consigue m ejorar una
d e te rm in a d a frecuencia. P ara 500At no se aprecia n in g ú n m áxim o
significativo en el espectro. La sensibilidad de la TDF para lOOOOAt es de 0.06
GHz, para lOOOOOAt es de 0.006 GHz, y para lOOOOOOAt es de 0.0006 GHz sin
embargo no se llega a mejorar la desviación de un 0.08% en el prim er modo,
0.31% dél segundo, 0.22% del tercero y 0.21% del cuarto. Para conseguir una
m ayor precisión con respecto a las frecuencias reales, con lo cual una m ayor
aproxim ación a la estructura de cam po real, será necesario efectuar un
m allado m as denso. Paradójicam ente, el error introducido en el cálculo de
la TDF para un num ero bajo de N corrige la desviación inherente al proceso
de sim ulación ("un clavo saca a otro clavo") y para pocas iteraciones
tem porales se consiguen mejores resultados.
125
TABLA 3.V. Frecuencias de córte de los m odos TE en guía rectangular hueca
obtenidos con DFDT. D esviación respecto de los resultados
analíticos para d istin to n ú m ero de iteraciones tem porales
(NAt). Dimensiones de la guía a/b= 2, a=2cm. M allado 40x20.
126
precisión en los resultados pasan por efectuar un m allado mas denso que se
adapte m ejor a los contornos. Se va a analizar una guía circular con radio de
0.5 cm para los m odos Transversales M agnéticos (TM). Se com pararan los
re s u lta d o s o b tenidos con DFDT con los resu ltad o s que se d eriv an
analíticam ente [87] y se observará como varían estos al efectuar m allados
m as densos.
M allado
R2 = (FLOAT(DIM/2))**2
ATEl = DIM - 1
ATI2 = (D IM /2) - 1
RR = FLOAT(ATI3)
DO 12 J12=0,ATI1/1
AURl=SQRT(R2-(FLOAT(J12)-RR)**2)
IF(DIM.EQ.20)THEN
127
IF((J12.EQ.O).OR.(J12.EQ.19))THEN
AUI3=3
END IF
END IF
TOIl(J12)= ATI3-AUI3
TOI2(J12)= ATI3+AUI3
12 CONTINUE
C onvergencia
128
Figura 3.9 Mallado de la guía circular hueca.
TM 21 49.0078 49.6157 1 .2 0
130
20 25 30 35 40 45 50 55 20 25 30 35 40 45 50 55
(a) (b)
20 25 30 35 40 45 50 oo 20 25 30 35 40 45 50 55
(c) (d )
2 0 2 5 3 0 3 5 4 0 4 5 5 0 5 5 20 25 30 35 40 45 50 55
(e) (f)
Figura 3.11. Espectro para cada discretización,r = 0.5 cm, frecuencia en GHz:
(a) r/10, (b) r/20, (c) r/30, (d) r/40, (e) r/50, (f) r/60.
131
TABLA 3.VII. Frecuencias de corte de los cuatro prim eros m odos TM en
guía circular hueca obtenidos con DFDT. Desviación respecto
de los re s u lta d o s a n a lític o s p a ra v a rio s m a lla d o s.
Dimensiones de la guía r = 0.5 cm.
r/A h TM 1 0 % TM ll % TM 21 % TM 0 2 %
132
3.5. Guía coaxial de sección rectangular. Resultados
numéricos.
Es sabido que una guía de ondas con sección cuadrada o circular,
p u e d e so p o rtar dos ondas in d ep en d ien tes con polarización vertical y
horizontal, respectivam ente. Lo m ism o es aplicable a las guías coaxiales
circulares y rectangulares. El m odo usual de operación de estas últim as suele
ser el TEM . Esto sucede en las m icrocintas blindadas consideradas como
coaxiales rectangulares. Pero se dan situaciones en que los m odos de tipo TE
y TM deben tenerse en cuenta [8 8 ]. Las guías coaxiales circulares han sido las
m as u tilizad as en el pasado, pero se observan aplicaciones en que es
preferible el uso de coaxiales rectangulares. Estas últim as presentan ventajas
en su m ecanizado por presentar paredes planas, asi como por presentar un
p lano de polarización bien definido que perm ita discrim inar con m as
facilidad la polarización cruzada.
Asi como los m odos de orden superior en una guía circular coaxial
son bien conocidos [89], [8 8 ] , este hecho contrasta fuertem ente con los pocos
estudios realizados sobre la guía coaxial cuadrada. La bibliografía referente a
este tipo de guía es m uy incom pleta e inadecuada para m uchos propósitos.
En u n trabajo de Bezlyudova de 1963, [90], se describe un m étodo para la
determ inación de los m odos propios m as bajos de líneas de transm isión con
conductor interior y exterior de sección rectangular (TEl0 ,TE 0 l) pero no se da
inform ación explícita aplicable a las cuadradas. Brackelmann en un trabajo
de 1967, [91], analiza líneas rectangulares con conductor interior y exterior
con centros no coincidentes, tratan d o el caso del coaxial como u n caso
particular, unas pocas frecuencias de corte de unos pocos m odos del coaxial
p u e d en deducirse de su trabajo. T ourneur, [92], obtiene el espectro del
coaxial cuadrado haciendo uso de elem entos finitos y de un procedim iento
v aria cio n al de R ayleigh-R itz, su trabajo m u estra las curvas de las
frecuencias de corte de los m odos TElO, T E ll, TE2 0 y T M ll para razones de
b / a d e sd e 0 a 0.3 (dim ensiones de c o n d u cto r in te rio r y exterior
respectivam ente). Finalmente, G runer [8 8 ], [93], [94], [95] describe de forma
in d e p e n d ie n te una técnica sim ilar a la u tiliza d a p or B ezlyudova y
B rackelm ann pero con una diferente im plem entación aplicable a líneas
coaxiales rectangulares. G runer obtiene el espectro de una guía coaxial
133
c u ad rad a con conductor interior y exterior de dim ensiones arbitrarias
extrapolando resultados desde b / a = 0 hasta b /a = l.
134
TABLA 3.VIII. L ongitudes de onda de corte n orm alizadas en coaxial
cuadrado para los m odos TElO/TEOl, T E ll para varias
relaciones de b /a . Resultados por field m atching de Gruner,
y desviación de los resultados por DFDT con respecto a los
de "field matching".
b /a Xc T E io /a ( f ie ld m a t c h in g ) Xc T E lO /a (FD-TD) d e s v i a c io n ( % )
135
TABLA 3.IX. L ongitudes de onda de corte norm alizadas X c/a en coaxial
c u a d rad o p ara los ocho p rim ero s m odos TE en coaxial
cuadrado para varias relaciones de b /a .
136
Las figuras 3.12.a-h m uestran para el caso b /a = 0.5, la geometría de los
ocho prim eros m odos TE (líneas de Hz constante) y las 3.13.a-h los ocho
prim eros TM (líneas de Ez constante). Dado que el prim er modo de orden
su p e rio r aparece degenerado [88], (TElO con el TEOl), la figura 3.12.a
representa la superposición de los dos modos.
(g) (h )
137
(d)
(e) (f)
(g) (h)
138
3.5.2. Guía coaxial de sección rectangular.
La geometría de la guía con sección cuadrada es simétrica respecto de
cuatro ejes, lo cual hace lógico el hecho de que aparezcan degeneraciones.
Una form a de rom per la degeneración es analizar un caso en que la sección
no es cu ad rad a. A continuación se efectúa un breve análisis del caso
rectan g u lar. En particular el m odo TElO aparece en solitario. La guía
analizada presenta unas relaciones entre dimensiones:
a /a ' = 2
b /b ' = 2
b/a=0.5
139
v
(a) (b)
63
(c) (d)
Io »o :
(a) (b)
(c) (d)
140
3.5.3. Análisis de guías con dos conductores
rectangulares con centros no coincidentes.
La im plem entación en el ordenador del algoritm o num érico DFDT es
lo suficientem ente rápida como para perm itir la obtención de resultados de
form a iterativa. Esta posibilidad ha sido explotada analizando una guía con
dos conductores rectangulares con centros no coincidentes. Se calcula la
distribución de campo así como las frecuencias de corte de varios de los
prim eros m odos de tipo TE. Se evalúan estas características en el segundo
caso, variando la posición del conductor interior.
a '/a = 0.5
b'/b = 0.5
b /a = 0.25
Siendo a y a ' las dim ensiones horizontal y vertical del conductor exterior, y
b y b' las dim ensiones horizontal y vertical del conductor interior. El
conductor interior esta situado a una altura X en el centro del conductor
exterior, esto se ve en la figura 3.16.
141
TABLA 3.XI. L ongitudes de onda de corte norm alizadas (A.c/a) en guía
rectangular con u n conductor rectangular insertado (figura
3.16), para los cuatro prim eros m odos TE para varias alturas X,
a'/a=0.5, b'/b=0.5, y b/a=0.25.
Las figuras 3.17.a-d m uestran las líneas H z=constante, para los cuatro
prim eros m odos TE siendo X/a=0. Si X/a=0.094 la distribución de campo es
la que m uestran las figuras 3.18.a-d.
ái
142
(a) (b)
(c) (d )
Figura 3.17. Modos TE, líneas de Hz constante, X/a=0: (a) l 2 TE; (b) 22 TE; (c)
32 TE; (d) 42 TE.
a »
(a) (b)
44
(c) (d)
Figura 3.18. M odos TE, líneas de H z constante, X/a=0.094: (a) l 2 TE; (b) 22 TE;
(c) 32 TE; (d) 42 TE.
143
3.6. Conclusiones.
Se ha hecho una extensa revisión de la labor efectuada p o r otros
autores en la aplicación del m étodo DFDT a problem as de valores propios.
144
consigue u n ajuste perfecto con la sección de la guía. En estos casos se ha
visto que no tiene dem asiada influencia la variación del tam año de la celda.
Por últim o y como conclusión final cabe decir que: El reducido costo
en m em oria y el poco consumo de tiempo de CPU unido a la flexibilidad del
m étodo DFDT y su precisión hacen que en nuestra opinión sea un serio
rival al m étodo de Elementos Finitos en el tratam iento de estos problem as.
En prim er lugar por el tratam iento de las puntas y demas singularidades del
campo que m ientras que con FE requieren un mallado especial, en el caso de
DFDT pasan inadvertidas, incluso con m allados groseros. En segundo lugar
145
p orque los Elem entos Finitos dan soluciones falsas (m odos espúreos [84])
que p roceden exclusivam ente del m étodo num érico em pleado, esto no se
ha en co n trado en n in g ú n m om ento aplicando DFDT. Por últim o cabe
a ñ a d ir que resu lta m ucho m as fácil confeccionar u n p ro g ram a de
D iferencias Finitas que uno de Elem entos Finitos. A dem ás consideram os
que el m étodo DFDT no está lo suficientem ente explotado, y con esta nueva
ap licació n p rete n d em o s m o stra r u n a p o te n c ia lid a d h asta ahora no
m anifestada por otros autores, y es la de utilizar DFDT para el análisis de los
m odos TM y TE en guías de onda de sección arbitraria.
146
Capítulo 4: Aplicación de DFDT al
análisis de discontinuidades en el
plano H
4.1. Introducción.
A lo largo de los últim os veinte años el análisis de cicuitos se
extendió a las frecuencias de m icroondas. Esto se hizo m ediante la
in c o rp o ra c ió n de m o d elo s que exp licasen el fu n cio n a m ie n to de
d eterm in ad o s elem entos pasivos del circuito. A lgunos de los m odelos
em pleados en el diseño se obtenían de la teoría. Otros eran combinaciones
híbridas entre teoría y m edidas experimentales. Sin embargo pocos circuitos
eran com pletam ente m odelizados haciendo uso exclusivo de las ecuaciones
de M axwell. Esto ha hecho que casi toda la tecnología desarrollada en
radio frecu en cia haya sido conseguida con aproxim aciones estáticas o
quasiestáticas. Por otro lado según el diseño se ha ido complicando, ya no
resulta rentable la técnica de fabricar y probar para obtener un prototipo que
funcione bien.
147
a) Como se basan en el principio de superposición, no pueden aplicarse a
p ro b le m a s n o -lin e a le s , p ro b le m a s en q u e la s c a ra c te rís tic a s
electromagnéticas dependen de la am plitud del campo aplicado.
148
obtener buenos resultados, acusan de ello al elevado factor de calidad de la
estructura m odelada.
4.2. Formulación.
Las d isc o n tin u id ad e s en guía de onda en el plano H , son
discontinuidades en las que se rom pe la sim etría traslacional sin cambiar el
plano de propagación del campo eléctrico (que contiene al eje longitudinal).
Este plano continua siendo el mismo, solo que cambia la dirección de la guía
(eje longitudinal). El modo con el que se trabaja es con frecuencia el m odo
fundam ental. La figura 4.1 m uestra una discontinuidad con tres puertas en
el plano H , aunque en general pueden tener cualquier núm ero de puertas.
Las paredes de la guía son conductores perfectos y el dieléctrico de relleno
tiene contantes eléctricas £, p. Las ecuaciones de M axwell del rotacional
gobiernan la evolución del campo en estas estructuras, y puesto que la banda
149
de o p eración suele ser la del m odo fu n d am en tal (TElO en una guía
rectangular), es suficiente con efectuar un análisis para el caso de que se
propague un m odo tipo TE en dos dim ensiones (ecuaciones (2.5), el criterio
en las guías de onda esta cam biado con respecto a la propagación en el
espacio libre).
plano de excitación
150
4.2.1. Modelado.
El m allado de la discontinuidad se efectúa utilizando la celda de la
figura 2.16, (TM en el espacio libre, TE en guía), y las ecuaciones en
diferencias a resolver (DFDT) son las (2.11). Se utilizarán celdas cuadradas de
lado Ah=W /20, siendo W la anchura de la guía. Esta discretización perm itirá
sim ular con suficiente resolución la evolución del m odo fundam ental,
(apartado 2.3.7). La relación de estabilidad que relaciona el paso tem poral con
el espacial (ecuación 2.27.b) se fija:
Para sim ular la adaptación a la entrada de cada una de las puertas que
confluyen a la discontinuidad se hace uso de las Condiciones de Contorno
Absorbentes, (CCA) (apartado 2.5). Estas se definen en el m allado a la salida
de cada p u e rta , y d u ran te la sim ulación DFDT absorben las ondas
procedentes del interior del dom inio de cálculo, que en este caso es la
m odelización de la discontinuidad.
4.3. Excitación.
Como se ha com entado con anterioridad la banda usual de operación
de las guías de onda es la del m odo fundam ental. Por esa razón, para evitar
reflexiones indeseadas y sim ular el com portam iento real del dispositivo, se
introduce la excitación tem poral con la estructura transversal espacial del
m odo TElO. Considerando el eje x como el eje a lo largo del cual se prolonga
la guía, esto es:
151
La dependencia tem poral, como se ha explicado en 2.4 p uede ser de
dos tipos, arm ónica (CW en la bibliografía anglosajona) ó m ediante el uso
de pulsos tem porales.
Se trabaja en el dom inio del tiem po con una sola frecuencia tem poral
f, (CW en la bibliografía). Si la excitación se sitúa en x = íq A 1 \, el cam po Ez
incidente en esta línea es:
Medio paso espacial hacia atrás, el campo H y incidente (ver figura 2.16 con la
celda u n id ad rem unerada) es:
Hyn+1/2inc(i0-l,j) = (4.4)
-YTEloeTElo(jAh)V0sen(27tf(n+l/2)At+PTEioAh/2)
siendo
YTE1 0 - P tE 1( / (2 n í ^
(4.5.a)
152
De este m odo, la introdución de la señal en x=igAh se realiza según la
(c=^ hx
eAh
Hyn+1 / 2 (i0 -l,j) = Hyn+1/ 2 (i0 -l,j) Iyee - D Ezninc(i0 ,j) (4.7)
donde [V]r y [V] 1 son los vectores de tensión incidente y dispersada, [16].
Cada com ponente de estos vectores, es la tensión incidente/reflejada en cada
puerta, y [S] es la m atriz de dispersión del dispositivo.
153
cam po transm itido Vj, y el parám etro S^(f) se obtiene por sim ple división
de las am plitudes transm itida/incidente:
(4.9.a)
Sij V i w
con lo cual
154
A0 exP ' “a t< 2 t0
f(t) = v y
L0.0 t>=2 t0
(4.10)
155
cam po en pun to s suficientem ente alejados de la discontinuidad, p or las
m ismas razones que se comentan en 4.3.1 para el caso de señales armónicas.
N -l
Vj(fk)= N eTFl2na-Ah) X Ezn(ij'íj) exp(-j27tfk nAt) (4.12)
TE1° J n= 0
fk = k /N A t, k=0,l,...,N /2.
156
velocidad de fase diferente. Según nos acercamos al corte, esta velocidad de
fase d ism in u y e h a sta anularse en el corte. A la vista de esto nunca
podrem os esperar a que el campo se vaya del dispositivo a través de las
puertas. Siem pre quedara campo con una frecuencia próxim a al corte pues
no hay perdidas en el m odelado. De aqui se deduce que al efectuar la TF del
cam po en las p u ertas de salida no lo estam os haciendo sobre una señal
tem poral finita, sino sobre un conjunto infinito de arm ónicos cada uno de
ellos saliendo del dispositivo con una velocidad. Adem ás si el pulso excede
en frecuencia la frecuencia de corte del seg u n d o m odo, se estarán
p ro p ag an d o otras frecuencias que no son de interés para el análisis. El
últim o p ro b lem a que se ha visto que presenta el uso de pulsos, es el
fenóm eno de Gibbs que nos obliga a hacer uso de señales tem porales suaves
de duración finita .
157
4.3.3. Señales sintéticas de banda limitada.
Excitación
m
f(nAt) = ^ Vj sen(27Ef¡nAt) (4.14)
1=1
m
Ez ninc(Í0'j) = eTEio^Ah') X V! se n (2TtflnAt) (4.15)
1=1
Medio paso espacial hacia atrás, el campo H y incidente (ver figura 2.16 con la
celda un id ad rem unerada) es:
H yn+1/2inc(i0-l,j)= (4.16)
m
”^%TEio eTEic$^^ V(f])sen(27if](n+l/2)At+ Piypg^Ah/2)
1=1
siendo,
158
Y1/TE10 “ P l T E l f / ^ r á ^ (4.17.a)
(4.7). En ambos casos con las nuevas expresiones de los campos (4.15) y (4.16).
N -l
Vj(fi)
(4.19)
159
4.4. Condiciones de Contorno Absorbentes.
La ad ap tació n en las guías se sim ula m ediante la aplicación de
C ondiciones de C ontorno A bsorbentes. Estas se definen en el borde del
m allado co rrespondiente al final de la guía para sim ular que esta se
prolonga hasta el infinito. La elección de las CCA, ante tantas posibilidades,
parece com plicada, sin em bargo teniendo en cuenta la relación prestaciones-
econom ía (m em oria y tiem po de cálculo), la decisión se decanta p o r la
condición de M ur de prim er orden y las CCA de segundo orden. Adem ás
son las usadas por la m ayoría de los autores, [75], [27], [43], [51]. Las CCA de
prim er y segundo orden absorben perfectam ente una onda plana que incide
norm alm ente sobre la pared en la que están definidas, pero en la sim ulación
d el avance p o r u n a guía de o n d a s no tenem os u n a o n d a p lan a
propagándose a lo largo de su eje. La propagación del m odo fundam ental en
una guía de ondas rectangular se puede descom poner como combinación de
dos ondas p lan as viajando cada u n a con u n ángulo. Por esa razón,
previam ente a la aplicación de las CCA en la sim ulación de circuitos mas
com plicados, se va a estudiar su com portam iento en el caso sencillo de la
sim ulación de un pedazo de guía.
160
norm alizadas), la am plitud se tom a la m ism a para todas (Vj(fp=100.0). La
excitación se sitúa según la ecuación (4.6) en i=200, y para elim inar la señal
incidente a la izquierda de i = 2 0 0 (zona de cam po reflejado) se usa la
ecuación (4.7). En ambos casos con las expresiones de los cam pos (4.15) y
(4.16) para las frecuencias de trabajo com entadas antes.
161
TABLA 4.1. Máximos y m ínim os del campo entre la excitación y las CCA de
M ur de prim er orden. Relación de ondas estacionarias y m ódulo
del coeficiente de reflexión y fase en función de la frecuencia de
la onda.
K qW / tc m ín im o m áxim o ROE P 0o
1.1 63.110 149.592 2.370 0.407 16.647
1.2 69.962 126.018 1.801 0.286 24.000
1.3 75.600 115.754 1.531 0.209 30.000
1.4 80.474 114.470 1.422 0.174 22.086
1.5 84.396 110.913 1.314 0.136 4.022
1.6 89.275 113.559 1.272 0.120 45.236
1.7 88.794 109.027 1.228 0.102 24.828
1.8 91.802 109.655 1.194 0.089 37.457
1.9 91.990 106.620 1.159 0.074 58.477
120
100
80
60
40
20
162
4.4.2. Condiciones de Contorno Absorbentes de
Mur de segundo orden.
Estas CCA son tam bién llam adas de tipo Padé (2,0), (ver apartado 2.5).
M odelado
Los resultados obtenidos son m ejores que con las CCA de prim er
orden. A p esar de ello se observa todavía una apreciable reflexión para
frecuencias m u y bajas.
163
excitación y las CCA de la derecha, la relación de ondas estacionarias y el
coeficiente de reflexión en las CCA.
TABLA 4.II. Máximos y m ínim os del campo entre la excitación y las CCA de
M ur de segundo orden. Relación de ondas estacionarias y
m ódulo del coeficiente de reflexión y fase en función de la
frecuencia de la onda.
K qW / ti m ín im o m áxim o ROE P 0o
1.1 83.968 117.444 1.399 0.166 16.647
1.2 91.659 107.766 1.176 0.081 24.134
1.3 95.008 104.006 1.095 0.045 30.000
1.4 96.867 102.835 1.062 0.030 26.667
1.5 97.865 101.703 1.039 0.019 4.022
1.6 98.637 101.314 1.027 0.013 45.000
1.7 99.107 100.819 1.017 0.009 20.571
1.8 99.363 100.601 1.012 0.006 46.316
1.9 99.545 100.372 1.008 0.004 58.776
120
100
80
60
40
20
164
4.4.3. Discusión.
Las CCA son necesarias en la m odelización de u n circuito, y su
introducción debe sim ular en todo lo posible la realidad física. Esto no es
posible, las C ondiciones de C ontorno A bsorbentes no son perfectas y
generan una cierta reflexión. Se refleja m ucho cam po en las de prim er
orden, en las de segundo orden la reflexión es sensiblem ente inferior, y en
ambos casos se observa que esta aum enta hacia la parte baja de la banda de
operación. En este sentido cabe decir que las de m ayor orden se com portan
m uy bien en la zona alta de la banda de operación TElO.
Desde el punto de vista circuital, las CCA se com portan como una
im p ed an cia. U na form a de co rreg ir su efecto sería d e te rm in a r su
im pedancia equivalente, y tenerla en consideración cuando se calculan los
parám etros de dispersión de un dispositivo.
165
Este últim o da un núm ero de pasos tem porales m uy elevado, excesivo
(3.98x1015At), (por el escaso espaciam iento en frecuencia de la señal). Ante
esto se ha o ptado po r tom ar u n num ero de pasos de sim ulación lo
suficientem ente grande como para contener varios periodos de cada señal
m onocrom ática. A dem ás como se efectúa una com paración frecuencia a
frecuencia, esto no tiene ninguna influencia en el análisis de la reflexión.
f = oo
90 °
(4 .20 )
Figura 4.4
166
Esto indica que dada una frecuencia de trabajo dentro de la banda TElO en
una guía de ondas rectangular, se conoce el ángulo con que las dos ondas
planas incidirían sobre una hipotética pared transversal a la guía. Y esto nos
perm itirá diseñar unas CCA exactas para cada frecuencia. Nos remitimos al
apartado 2.5.1, en particular a la form ulación general de H igdon [77], [78],
p ara d esarrollar unas CCA para ondas planas con ángulo de incidencia
cualquiera.
fd C d \
(4.21)
- En x = 0:
cAt/cosy-Ax
cAt/cosy+Ax
- En x = H:
U n+1(H, j) = Un (H -l, j)+Cf t^C° Sr^X ( u n+1(H -l, j) - Un (H, j)) (4.22.b)
cAt/cosy+Ax
Siendo cosy
R esultados num éricos
K qW / ti m ín im o m áxim o ROE P
1.1 99.191 100.387 1.0121 0.0060
1.2 99.595 100.155 1.0056 0.0028
1.3 99.775 100.110 1.0034 0.0017
1.4 99.853 100.086 1.0023 0.0012
1.5 99.854 100.115 1.0026 0.0013
1.6 99.846 100.131 1.0028 0.0014
1.7 99.804 100.164 1.0036 0.0018
1.8 99.787 100.201 1.0042 0.0021
1.9 99.725 100.273 1.0055 0.0027
120
100
80
60
40
20
0
0 50 100 150 200 250 300 350 400
x/Ah
169
4.5. Análisis de discontinuidades en el plano H.
Resultados numéricos.
El m étodo DFDT se aplica a la sim ulación y análisis de varias
discontinuidades en guía de ondas.
170
El cálculo de las 6000 iteraciones consume un m inuto y diez segundos
de tiem po de CPU en un IBM 6000/520, este tiempo se incrementa a treinta
m inutos en un Apollo 400t.
0.8
<N
0.6
<M
m
0.4
0.2
0
1 1.2 1.4 1.6 1.8 2
Ko W/ti
Figura 4.7. Coeficiente de transm isión de potencia del codo recto.
• DFDT; ------- EF; □ Análisis modal
171
4.5.2. Unión T en el plano H.
Se utiliza un m allado 400x200 donde cada brazo (uno perpendicular al
otro) viene dado por un m allado 200x20. En este caso por tener el sistema
tres puertas y por lo tanto confluir las reflexiones de tres paredes con CCA,
se re c u rre a las CCA de M u r de se g u n d o o rd en p o r su m ejor
com portam iento. Las características del m odelado son idénticas a las del caso
anterior.
172
o
0.8
0.6
0 .4
0.2
0
1 1.2 1 .4 1.6 1.8 2
K oW /n
173
350
300
250
100
50
174
4.5.3. Filtro.
Se m uestra en la figura 4.11. El dieléctrico de relleno tiene er=2.62. Se
utiliza un m allado 400x40, con Ah=W /40, esto se hace así para dar una
m ejor resolución espacial a la estructura.
0.25W
0.3W
0.35W
0.3W
Figura 4.11.
K qW / tc -20 log 1Sl 2 1 (DFDT) -20 log 1Sl2 1 (EF) d % 1Sil 12+ 1Sl2 12(DFDD
1 .1 48.457 - - 0.887
1 .2 40.751 40.0 1.90 0.873
1.3 33.545 33.5 0.13 0.930
1.4 25.548 25.5 0.19 0.953
1.5 11.610 11.5 1 .0 0 0.890
1.547 0.876 0 .0 - 0.977
1.553 3.025 3.0 0.80 1.069
1 .6 11.903 11.5 3.50 0.966
1.7 16.739 16.0 4.60 1 .0 0 0
176
*00
(b)
Figura 4.13. Distribución de Ez en el filtro, K qW / 7 C = 1.547: (a) t=200At;
(b)t=400At.
177
4.5.4. Guía de planos paralelos abierta.
La guía de planos paralelos abierta al espacio libre ("flanged parallel-
plate w aveguide" en la literatura anglosajona) se m uestra en la figura 4.13.
El problem a de la radiación de esta estructura no es nueva, y ha sido tratada
usando varios m étodos [113]. Esto ha generado cierta discusión entorno a
cual de los m étodos em pleados da resultados mas precisos [114]. En este
apartado se analiza esta discontinuidad. Se concluye que el algoritmo DFDT
es bastante adecuado para el tratam iento de discontinuidades en el plano H
y que incluso puede abordar óptim am ente problem as de radiación con una
aceptable precisión. Por otro lado se confirm a que los resultados m as
precisos son los obtenidos por Lee y G run [114] que coinciden con los de
Leong, Koi y Chandra [113] por el m étodo de los Momentos.
178
Esta estructura en el plano H es idéntica a una guía rectangular que se
abre al espacio abierto, y el m odo fundam ental el TEl tiene la m isma
d istrib u ció n que el TElO en la guía. Ante esta sim ilitud se efectúa el
m odelado como el de una guía que se abre en el plano H con un mallado
400x220, siendo Ah=W /20. Tam bién se efectúa un análisis frecuencia a
frecuencia. Se usan las CCA de prim er orden modificadas del apartado 4.4.4
p ara el interior de la guía y las de M ur de segundo orden para el espacio
abierto exterior. En las esquinas del m allado del espacio abierto exterior, se
usan unas CCA de prim er orden. En estas últim as se supone que la onda
incide sobre las esquinas siguiendo la diagonal de las celdas, y se asigna en
este punto para cada instante el campo que había dos pasos temporales antes
en la línea de la diagonal un paso espacial antes [63], estas se detallan en el
apéndice 4.A.
179
o
180
Figura 4.14.b. Distribución de campo Ez en la guía de planos paralelos
abierta, instante 800At. Vista en perspectiva y líneas de Ez constante.
TABLA 4.V. Coeficiente de relexión para cada frecuencia en la guía de
planos paralelos abierta al espacio libre (figura 4.13). Resultados
con DFDT com parados con los obtenidos por Leong, Koi y
Chandra, [113], usando el m étodo de los Momentos.
182
4.6. Conclusiones.
Se ha presentado un algoritm o alternativo y eficiente para el análisis
de disco n tin u id ades en el plano H. Esto ha sido posible con el uso de
D ifere n cia s F in itas en el D o m in io del T iem po, el tra ta m ie n to
b id im en sio n al, y la utilización de señales sintéticas obtenidas como
com binación lineal de señales m onocrom áticas. Las C ondiciones de
C ontorno A bsorbentes han sido analizadas con detalle en este tipo de
aplicaciones, y se han desarrollado unas nuevas CCA que funcionan
perfectam ente, pero para una sóla frecuencia. Se observa que con este
tratam iento y la separación cam po total-cam po dispersado, la excitación
puede situarse todo lo cerca que se desee de la discontinuidad. La variación
de su situación no influye sobre los resultados. Se observa que el uso de
señales sintéticas es ventajoso por perm itir la obtención de parám etros para
varias frecuencias con una sola sim ulación, pero hay que hacer las puertas
de la de la discontinuidad con una longitud de varias longitudes de onda de
la frecuencia m as baja para discernir la estructura de ondas estacionarias, a la
que luego hay que aplicar la corrección de las CCA. Además hay que hacer la
TDF sobre m uchos puntos. Por otro lado con la aplicación de las CCA de
prim er orden m odificadas son necesarias varias simulaciones pero la TDF se
calcula en u n solo pun to en cada puerta; las p u e rta s no deben ser
necesariam ente tan largas, y basta una longitud de onda de la frecuencia mas
baja. Por últim o añadir que estas CCA no precisan de corrección alguna.
Por otro lado, con los últim os casos abordados, se dem uestra que los
cálculos son siem pre susceptibles de m ejora, pero a costa de un m ayor
consum o de recursos de ordenador, por lo tanto hay que ajustarse a un
com prom iso recursos-precisión. Cuando se presenta una resonancia hay que
183
esperar m ucho tiem po de sim ulación al estado estacionario antes de aplicar
la TDF. Y cuando se presenta u n problem a de radiación, hay que procurar
que las CCA del espacio exterior se com porten bien, lo cual es difícil dado
que el frente de ondas que usualm ente incidirá sobre éstas no es plano
(a p a rta d o 2.5.4); o tra solu ció n consistiría en h acer u n m allad o lo
suficientem ente grande de la zona exterior para que se pueda calcular el
efecto de la discontinuidad antes de que llegue la onda reflejada en las CCA
exteriores.
184
Capítulo 5: Extensión de DFDT al
análisis de discontinuidades en tres
dim ensiones
5.1. Introducción.
Todo lo com entado en el apartado 4.1 del capítulo anterior se podría
repetir en este apartado. H abría que añadir que todos los análisis circuitales
en co ntrad o s h asta el m om ento en la bibliografía han sido efectuados
m ediante u n tratam iento tridim ensional que es m ás general. Este tipo de
tratam iento es absolutam ente necesario cuando no hay n in g ú n tipo de
sim etría que p u ed a sim plificar el problem a y que p u d iera perm itir un
an álisis en dos dim ensiones. Esta red u cció n es posible en alg u n as
discontinuidades, como es el caso de las discontinuidades en el plano H (ver
capítulo anterior), pero en general es casi siem pre necesario un análisis en
tres dim ensiones.
185
Laplace num éricam ente u sando diferencias finitas con u n algoritm o de
relajación. Se extienden los algoritm os DFDT aplicados con anterioridad al
estudio de discontinuidades en el plano H . Las CCA, excitación, y dem ás
algoritm os se aplican en tres dim ensiones (3D), pero con una salvedad,
ahora se trata de abordar u n problem a en 3D. La disponibilidad de recursos
en este caso no es trivial, y se hace necesaria una optim ización exahustiva
de todos los algoritmos. Se estudia la utilización de un m allado variable que
p erm ita ah o rrar recursos inform áticos. Se sim ula un caso con resultados
analíticos conocidos para chequear el program a, y finalm ente se aplica la
m etodología, con toda la inform ación generada en apartados anteriores, al
estudio de una unión T en línea coaxial cuadrada.
186
im pedancia de dicha línea de transm isión en función de los parám etros que
determ inan la sección transversal. La transform ación que efectúa Riblet da
como u n caso particular la situación de conductores cuadrados, ya resuelta
p or Bowm an [106].
187
La u n ió n T es un elem ento básico en la realización de m uchos
elem entos pasivos como p o r ejemplo acopladores direccionales. Se h a n
efectuado intentos de analizar esta discontinuidad haciendo uso del m étodo
de la resonancia transversal [103]. El m étodo de la resonancia transversal
perm ite obtener datos de diseño a partir de las longitudes resonantes de la
línea que contienen la d isco n tin u id ad [107]. Sorrentino [17] calcula las
longitudes resonantes m ediante la aplicación de análisis m odal [1 1 2 ], [16], al
conjunto de cavidades paralelepipédicas que constituyen la cavidad total que
encierra la discontinuidad. Sin embargo el procedim iento es tedioso. En este
capítulo se aborda el mismo problem a con otra técnica, DFDT. El m étodo es
sencillo de im plem entar, rápido en la obtención de datos, y m ás general,
p u e sto q u e el análisis en el dom inio del tiem po p erm ite reso lv er
situaciones con dieléctricos anisótropos, no lineales y con estru ctu ras
arbitrarias que no son fáciles de resolver con m étodos espectrales [107],
adem ás posee el atractivo adicional de la sim ulación tem poral, de perm itir
la visualización del fenóm eno de propagación en el dispositivo, tal y como
sucede en la realidad.
188
es im portante con vista a posibles am pliaciones futuras. La idea es que el
p ro g ra m a conste de u n conju n to de su b p ro g ra m a s en p rin cip io
independientes, que perm itan abordar cualquier problem a electromagnético
en tres dimensiones. Los bloques principales del program a son los de cálculo
del campo eléctrico y magnético. En estos dos bloques es donde se consumirá
la prácticam ente totalidad de recursos de m em oria y tiem po de CPU. Otros
bloques con m enos consum o de recursos, pero no m enos im portantes son
los de las CCA, los de introducción de la excitación, los de extracción de la
excitación (para separación cam po-total/cam po-dispersado) y los de cálculo
de la Transform ada Discreta de Fourier (TDF).
189
solam ente se considerarán dos posibles condiciones de contorno, que son las
m as frecuentes, estas son pared eléctrica (cam po eléctrico tangencial y
m agnético norm al nulos) y pared m agnética (campo eléctrico norm al y
m agnético tangencial nulos). Para diferenciar estas dos posibilidades en cada
cara (seis caras), se define un vector con seis componentes, cada una de ellas
referida a una de las caras, btipo(l:6), ver figura 5.1, en la que se ha
considerado la propagación en el eje x. Las posibilidades que se barajan para
la cara "i” (figura 5.1) son sólo dos:
Previam ente al cálculo del campo, el vector btipo habrá fijado que
condiciones de contorno deberán aplicarse a cada una de las caras del
paralelepípedo. En caso de que las condiciones de contorno en una cara no
sean específicam ente las de paredes eléctrica y magnética, el vector btipo
h ab rá dejado sin calcular las com ponentes del cam po en dicha cara,
perm itiendo establecer en esta cualquier otro tipo de condición.
2
Figura 5.1. Paralelepípedo de cálculo.
190
Límites en el cálculo del campo eléctrico
191
si btipo( 6 )=0 , se calcula desde k= l
192
si btipo(5)=l, se calcula hasta k=L3
193
Conexión entre bloques
194
accesibles al cálculo por ser exteriores a la discontinuidad que se modela. La
ventaja que sin em bargo presenta el uso de una m atriz global está en un
m en o r núm ero de cálculos y en consecuencia una m ayor rap id ez en la
sim ulación. Adem ás como los cálculos que se efectúan con las m atrices son
e n te ram en te v ectorizables, cuanto m ay o r es la m atriz m ayor es el
rendim iento que se obtiene del procesador vectorial.
195
persiste el m odo fundam ental y donde se determ ina la am plitud del campo
transm itido.
196
5.3.2. Modo Fundamental (TEM).
Com o se com enta en la in tro d u cció n del capítulo 3, hay una
diferencia epistem ológica entre una línea de transm isión y una guía de
o ndas. Las p rim eras se caracterizan p o r ser el m odo TEM el m odo
dom inante. Estos m odos no presentan com ponente axial de los campos, de
m odo que recurriendo a las ecuaciones de Maxwell, se llega a que el campo
eléctrico tran sv ersal p u e d e exp resarse como el g rad ien te (gradiente
transversal) de una función potencial [16], que contiene la dependencia del
cam po en las coordenadas transversales. La dependencia en la variable axial
es de la forma exp(jpx), siendo P=coVpS/ y "x" la variable axial,
1* = - Vt O (5.1)
, _ “7^ d ~r> 9
siendo Vf = i — + k — .
1 J dy dz
Vt <D= 0 (5.2)
197
Solución de la ecuación de Laplace en el coaxial
3 2 o a 2 í>
— + =0 (5.4)
1 en el conductor interior
0 .0 en el conductor exterior
(5.5)
198
p a r e d e s c o n d u c to r a s
paredes de simetría
Tom ando los mismos increm entos espaciales, Ay=Az=Ah y con la notación
discreta (y+Ay,z)=(j+l,k), queda,
199
a) para u n punto interior entre los dos conductores,
<D(j,k) = 0 (5.7.b)
j=0,l,...,La-l; k=Ha
j=La; k=0,l,...,Ha
0(j,k) = 1 (5.7.c)
j=0,l,...,Lb-l; k=Hb
j=Lb; k=0,l,...,Hb
j=0; k=Hb,...,Ha
j=Lb,..vL a; k=0
200
A lgoritm o para resolver el sistema en diferencias finitas
Bloque de entradas:
- Entradas
-Cota de error, er
Bloque de cálculo:
- error=negativo
- 0 a= 0 (j,k)
- Se calcula e= I Oa-<J>(j,k) I
- Se compara e con er
201
- Si e > er, entonces
- Final repetición
í>(j+l,k) - í>(j-l,k)
Ey(j,k) = -------- 2Áh (5.8.a)
P(j,k+1) - <£(j,k-l)
Ez(j,k) — 2 Ah (5.8.b)
(59b)
j=0,l,...,La-l; k=Ha
202
b. 2 ) pared vertical
3®(yO-4<D(j-l<k) + ®(j-2,k)
Ey(],k) = ----------------2Ah---------------- (5-9x)
j=La; k=0,l,...,Ha-l
j=Lb; k=0,l,—,Hb-l
Q(j/k+l) - flXjJc-l)
Ez(j,k) — 2Ah (5.11.b)
j=0; k=Hb,...,Ha
203
(5.12.a)
j^b^.JL a; k=0
(b)
205
C onsiderando que la energía eléctrica alm acenada p or un id ad de
longitud (We ) se calcula haciendo uso de la siguiente expresión,
We = | J e IVtO|2 ds
S
La im pedancia es,
C= 68.107 p F /m
Z= 48.976 Q
lo cual arroja u n error de u n 2 % respecto del valor exacto para este tipo de
línea (50 Q. ) .
206
5.3.3. Modelización de la guía coaxial cuadrada
(DFDT).
La modelización de un pedazo de guía coaxial cuadrada se efectúa con
u n a sola m atriz para cada com ponente del cam po en el m allado. Se
aprovecha uno de los dos planos de sim etría, el horizontal, y se establecen
tres bloques, que se m uestran en la figura 5.5. El plano inferior por ser de
sim etría, se m odela con una pared m agnética. La discretización se efectúa
con celdas cuadradas Ax=Ay=Az=Ah, y el paso tem poral de sim ulación se
elige según la condición de C ourant (apartado 2.3.5 (2.27.a))
At=0.5Ah/c (5.13)
z
Figura 5.5. Paralelepípedos de cálculo para la guía coaxial.
207
5.3.4. Comportamiento de las CCA con el modo
fundamental.
U n paso previo al análisis en el dom inio del tiem po debe ser
necesariam ente el estudio del com portam iento de las C ondiciones de
C o n to rn o A bsorbentes. Sin estas es im p o sib le la sim ulación de la
propagación, y se consideran adecuadas en la m edida en que estas sean
capaces de sim ular una línea infinita. Si las elegidas no funcionan bien este
paso las delatará, y nos perm itirá tom ar otras m as adecuadas. En la elección
se tien en en cuenta dos im p o rtan te s factores, su uso de recursos
informáticos y su funcionam iento como CCA y en consecuencia la precisión
de los resultados que querem os obtener.
M odelización
(5.14)
208
l i ( x Q , y , z , f ) = F(f) E^(y,z) (5.15)
siendo
f W /2 f
exp - = exp -3 « 5% (5.17)
lOAx
A =^=j-*12A t (5.18)
Con esta elección la máxima frecuencia que se puede extraer del espectro es
[47],
209
Introducción de la excitación y CCA
210
Resultados numéricos
600At
(a)
700At
x /A h
0.1 (b)
Q.
0.01
0.001
0 0.2 0.4 0.6 0.8 1
f/f.m a x
211
5.3.5. Estudio de la atenuación de los modos de
orden superior por debajo del corte.
En este apartado se va a proceder al análisis de la evolución de los
m odos TE y TM cuando estos se generan por debajo de la frecuencia de corte,
en la banda de operación TEM. Asimismo se analizará la reconstrucción del
m odo fu n d am en tal y la reflexión cuando la estru ctu ra espacial de la
excitación no se corresponde con la de este. Se excita la guía con una
estru ctu ra espacial del cam po que es constante en toda la sección, en
principio arbitraria. Por com odidad, para obtener resultados, se opta por una
excitación tem poral sintética de arm ónicos. Esta se com pondrá de varias
frecuencias, estando com prendidas en el rango TEM en una banda de interés.,
que se justificará m as adelante.
212
TABLA 5.1. Factor de propagación de los dos prim eros m odos TE para las
seis frecuencias del análisis. Valor del cam po norm alizado a su
v alo r m áxim o (valor en la excitación), a una distancia de una
longitud de onda del punto de introducción de la excitación para la
frecuencia mas baja (A,l=2.99793.108 (m /s )/lG H z ), y a una distancia
de lOOAh de su valor máximo (Ah=a/50), (100Ah=0.07A.l=0.4^6).
f (GHz) e x p -P ^ o A l exp-P'j^100-07Xl
Pteio
213
Los resultados de la tabla anterior son claram ente indicativos de la
rápida atenuación de que son objeto los m odos TE y TM. Pero no hay que
olvidar que la atenuación de los m odos de orden superior depende en gran
m edida de la estructura espacial del campo inicial utilizado en la excitación.
Por esa razón aunque la tabla 5.1 da una orientación, esta no es suficiente, y
p ara d eterm inar con seg u rid ad donde los m odos de orden superior son
insignificantes habrá que analizar la evolución de las com ponentes del
cam po bajo d istin ta s d istrib u c io n e s espaciales de la excitación. A
continuación se efectúa u n análisis de la evolución de las com ponentes del
campo m ediante la sim ulación DFDT.
R esu ltad o s num éricos con DFDT p a ra la evolución del cam po con
excitación transversal
214
Para valorar la atenuación de los m odos TE, se observa la evolución
de la com ponente axial del campo m agnético Hx, a lo largo de la guía. La
figura 5.7 m uestra H x en función de la posición en los puntos (i,0,45), para
cada u n a de las frecuencias de la sim ulación (1 a 6 GHz). La figura 5.8
m uestra Hx(i,0,45)/Hx(0,0,45), y la 5.9 Hx(i,20,25)/Hx(0,20,25).
1 GHz
2 GHz
3 GHz
o
4 GHz
0.1
5 GHz, 6 GHz
0.001
100
x/Ah
0.1
0.01
0.001
1 0"4
1 0 '5
0
216
0.1
x
£
$
o '
¿
CN
'~x
0.01
£
0.001
o 20 40 60 100
x/A h
0.01
— e — 1 GHz
□ 2 GHz
LO
CN - O - 3 GHz
o" 0.001
CN
H— 4 G H z
X.
K
■A— 5 G H z
w"
10
,-4 ♦ — 6 GHz
LO
CN
o"
CN
X
w
0 20 40 60 80 100
x /A h
217
100
0.001
- Q 1 GHz
— B — 2 GHz
- O - 3 GHz
1— 4 G H z
- A— 5 G H z
— ♦ — 6 GHz
x"
w
LO
o"
x'
N
W
0.1
0 20 40 60 80 100
x/A h
218
R esu ltad o s n u m éricos con DFDT p a ra la evolución del cam po con
excitación longitudinal
m
Exninc(i(>j'k) = sen(7tk/La)sen( 7rj/La) V \ sen(27cfinAt)
1=1
(5.21)
siendo La=a/2Ah=25, m = 6 , fpl.0.1 GHz, V pl.O , y
219
Figura 5.13. A m p litu d de la com ponente axial del cam po m agnético
n o rm alizad a, H x(i,20,25)/H x (0 ,20,25). R esultados de 1 a 6
GHz.(i=x/Ah). Excitación longitudinal.
LO
CN
0.1
o"
CN
0.01
w
£
0.001
£
><
w* 10"4
220
100
LO
CN
ov
CN
10
w
LO
CN
o"
CN
><
W>
W
0.1 ■ ■ I 1__ I
20 40 60 80 100
x/Ah
10
LO
CN
ov
CN
1 r
w
1 GHz
LO
N' B — 2 GHz
o
- O - 3 GHz
0.1
w i— 4 G H z
A— 5 G H z
♦— 6 GHz
0.01
20 40 60 80 100
x/Ah
221
Análisis de resultados
222
todo el rango de frecuencias del estudio. Esta se puede justificar en las altas
frecuencias de corte de los modos con componente axial y la escasa variación
en frec u en c ia. Esto hace qu e los facto res de p ro p a g a c ió n sean
aproxim adam ente iguales. Por otro lado la longitud de la guía es m uy
reducida, con lo cual tam bién se contribuye con este factor a una m enor
diferenciación en este com portam iento.
223
5.3.6. Excitación con señal temporal sintética con
estructura espacial TEM. Separación campo total-
campo reflejado.
A continuación se desriben los campos incidentes TEM y el m odo de
introducir la excitación separando las zonas de cálculo de campo total y de
campo reflejado. Para ello se hace uso del campo transversal TEM deducido
del esquem a en diferencias finitas del apartado anterior 5.3.2.
(5.23)
En las expresiones para H ^ 1^2 y E.JJ hay que tener en cuenta que el
224
Haciendo uso de las ecuaciones (2.46.a) y (2.46.b) la excitación se introduce en
io:
1
E2n+1 ( v j , k) = E2n+1 ( v j, k)| yee-C H n+ 2 inc(i0'l'Í'k) (5.24.b)
At
siendo C =
£ Ah
2_ J_
H zn+ 2 ( V E ) , k) = H zn+ 2 (i,, j, k) | yee + D e £ in c(i0, j, k) (5.25.b)
At
siendo D = — — ,
JLLAh
225
Resultados numéricos
©— 1 GHz
B— 2 GHz
0.001
O - 3 GHz
4 GHz
ñ — 5 GHz
♦— 6 GHz
in
CN
o'
CN
100
m
CN
- O - - 2 GHz
o' - 0 - 3 GHz
CN 0.1
W -A— 5 GHz
> — 6 GHz
CN
o' 0.01
CN
tu
0.001
-100 -60 -20 20 60 100
x/ Ah
226
-© — 1 G H z
10 i 1 i 1-- 1-- 1-- 1-- 1-- 1-- 1-- r
-B — 2 G H z
O - 3 GHz
H— 4 GHz
LO
CN
O -A — 5 G H z
CN
>< -♦ — 6 GHz
w” -m m-
LO
o"
x
N
W
0.1 _ i_ L
LO
CN
V
1 * -#
o
CN
cT - 0 - 1 GHz
0.1
W — S — 2 GHz
LO - * - 3 GHz
o" 0.01 1— 4 G H z
CN
>< -A — 5 G H z
w 0.001 5 ^ -- -♦ — 6 GHz
-4 J I I I I I L
10
227
Discusión de los resultados
Los resu ltad o s que se obtienen son aceptables. Las com ponentes
axiales que se generan son m ero "ruido numérico". Los resultados son m uy
b u en o s sobre todo p a ra el cam po reflejado. Pero esto es u n a m era
com probación de algo que ha sido chequeado con a n terio rid ad , en el
ap artad o 5.3.4. Por ello esta será la excitación que se elegirá para analizar
discontinuidades en guía coaxial. A dem ás, el buen com portam iento de las
CCA de M ur de 1er orden (reflexión inferior al 0.1%) junto con la sencillez
de su aplicación, hace que estas sean las elegidas. Su com portam iento es
perfecto con la extructura de cam po TEM. La razón de ello es que la
estru ctu ra de cam po TEM avanza axialm ente y se p uede asim ilar a un
conjunto de ondas planas que viajan en la dirección de la guía. Si las CCA se
definen transversalm ente a la guía, y sobre estas incide un a estructura TEM
perfecta, entonces la absorbe perfectam ente porque serán ondas planas que
inciden norm alm ente, y en este caso las CCA de M ur de 1er orden funcionan
a la perfección.
b) V ariando el paso tem poral At junto con Ax, de m odo que se conserve la
relación de estabilidad.
229
TABLA 5.II. Coeficiente de reflexión por cambio de mallado en %.
f (GHz) 1 2 3 4 5
230
5.4. A nálisis de un a guía coaxial rellena de
dieléctrico.
Una forma de com probar la validez del m odelado del coaxial y del
algoritm o DFDT en tres dimensiones, es la de abordar un problem a simple
con solución analítica. La geom etría del caso que se trata es sim ple, se
considera un coaxial cuadrado de los m odelados con anterioridad, en el que
se introduce una sección de long itud finita con constante dieléctrica
distinta. La figura 5.21 m uestra la geometría del problema.
54301167763276304365^75473
^ ►
a
Figura 5.21. Problema de prueba.
231
en estas lám inas de separación la m edia de la de ambos lados (ver apartado
2.6.3 del capítulo 2).
2 (5.26)
La tabla 5.III.a m uestra los resultados num éricos com parados con los
analíticos para la sim ulación durante un periodo de 1 GHz. La tabla 5.III.b
m uestra los resultados calculando la TDF durante cuatro periodos.
TABLA 5.III.a. Resultados para TDF durante un periodo de f=lGHz.
233
TABLA 5.III.b. Resultados para TDF durante cuatro periodos de f=lGHz.
234
5.5. Análisis de una unión T sim étrica en guía
coaxial cuadrada.
La figura 5.22 m uestra la discontinuidad. El brazo principal (puerta 1 a
bi
p uerta 2) es coaxial de relación — = 0'4 con a = 10 mm. La im pedancia
puerta 2
puerta 3
236
C on esto, los lím ites de cada zona son los siguientes, con sus
correspondientes códigos de paredes:
C DIMENSIONAMIENTO
BTIACEX=a/Ah
BTIHCEX=a/2Ah
BTIACI1 =bl / Ah
BTIACI2=b2/Ah
BTIHCI2=b2/2Ah
BTILOPl=Ll
BTILOP2=L2
C ZONA A
BTIIXAA=0
BTISXAA=BTILOPl
BTIIYAA=0
BTISYAA=BTIHCEX
BTIIZAA=0
237
BTISZAA=(BTIACEX-BTIACIl)/2
C ZONA B
BTIIXBB=0
BTISXBB=BTILOPl
BTIIYBB=BTIHCI2
BTISYBB=BTIHCEX
BTIIZBB=((BTIACEX-BTIACIl)/2)
BTISZBB=((BTIACEX+BTIACI1)/2)
C ZONA C
BTIIXCC=0
BTISXCC=(BTILOPl-BTIACI2)/2
BTIIYCC=0
BTISYCC=BTIHCEX
BTIIZCC=(BTIACEX+BTIACI2)/2
BTISZCC=BTIACEX
C ZONAD
238
c
BTIIXDD=(BTILOPl +BTIACI2) /2
BTISXDD=BTILOPl
BTIIYDD=0
BTISYDD=BTIHCEX
BTIIZDD=(BTIACEX+BTIACI2)/2
BTISZDD=BTIACEX
C ZONAH
BTIIXHH=(BTILOPl-BTIACI2)/2
BTISXHH=(BTILOPl+BTLACI2) / 2
BTIIYHH=BTIHCI1
BTISYHH=BTIHCEX
BTIIZHH=(BTIACEX+BTIACI2)/2
BTISZHH=BTIACEX
C ZONA E
BTIIXEE=0
BTISXEE=BTILOP2
239
BTIIYEE=0
BTISYEE=BTIHCEX
BTIIZEE=(BTIACEX+BTIACI1)/2
BTISZEE=BTIACEX
C ZONA F
BTIIXFF=0
BTISXFF=BTILOP2
BTIIYFF=BTIHCI1
BTISYFF=BTIHCEX
BTIIZFF=(BTIACEX-BTIACIl)/2
BTISZFF=(BTIACEX+BTIACIl)/2
C ZONA G
BTIIXGG=0
BTISXGG=BTILOP2
BTÜYGG=0
BTISYGG=BTIHCEX
BTÜZGG=0
240
BTISZGG=(BTIACEX-BTI ACI1) / 2
y o- matriz 1
matriz 2 x linea B
lineas de campo eléctrico tangencial linea A
lineas de campo magnético tangencial
241
5.5.2. Esquema del programa principal.
El program a principal fluye con los siguientes bloques:
Entradas
Bloque DF
Bloque DFDT
- btipo = ( 0 0 0 0 0 0 ).
- btipo = ( 0 0 0 0 1 1 ).
242
- btipo=(l 0 1 0 0 0 ).
- btipo = ( 0 0 0 0 0 0 ).
- btipo = ( 0 0 0 0 1 1 ).
- btipo=(l 0 1 0 0 0 ).
- btipo = ( 0 0 0 0 0 0 ).
- btipo = ( 0 0 0 0 1 1 ).
- btipo=(l 0 1 0 0 0 ).
243
- CCA sobre Hx en puerta 1
- btipo = ( 0 0 0 0 0 0 ).
- btipo = ( 0 0 0 0 1 1 ).
- btipo=(l 0 1 0 0 0 ).
244
tarda mucho tiempo en salir de la estructura. Esto es debido a la ausencia de
p érdidas en el m odelado. Las reflexiones sucesivas en el conductor tardan
m ucho en atenuarse ya que la única forma de que se pierda potencia es a
través de las puertas de salida donde están definidas las condiciones de
contorno absorbentes (figuras 5.25.e-f).
245
(d)
246
Figura 5.25.e) t = 800At; f) t = 1800At.
247
Parámetros de dispersión de la T
Para la sim ulación se usan seis matrices, tres para las puertas 1-2 de
270x25x50 de reales de 4 bytes, y otras tres de 100x25x50 para la puerta 3.
248
Figura 5.26.a) t = 1200At, b) t = 1300At.
249
Figura 5.26.c) t = 1500At, d) t = 1600At.
Figura 5.27.a) t = HOOAt, b) t = 1200At.
251
W&zM
252
Los parám etros de dispersión se calculan por sim ple división de las
am plitudes del cam po calculadas con la TDF para cada frecuencia, como en
el a p a rtad o 4.3.3. Los resu ltad o s se com paran con m edidas obtenidas
experim entalm ente en la em presa Radiación y M icroondas S.A., [116].
253
TABLA 5.IV.a. Resultados para TDF durante un periodo de f=lGHz.
Excitación por la puerta 1.
254
TABLA 5.IV.b. Resultados para TDF durante un periodo de f=lGHz.
Excitación por la puerta 3.
255
TABLA 5.V.a. Resultados para TDF durante tres periodos de f=lGHz.
Excitación por la puerta 1.
256
TABLA 5.V.b. Resultados para TDF durante tres periodos de f=lGHz.
Exdtadón por la puerta 3.
257
TABLA 5.VI.a. Resultados para TDF durante quince periodos de f=lGHz.
Excitación por la puerta 1.
desv % 1 1 .0 0 1 0 .0 0 1 0 .0 0 9.00 8 .0 0 8 .0 0
258
TABLA 5.VI.b. Resultados para TDF durante quince periodos de f=lGHz.
Excitación por la puerta 3.
Lo p rim ero que se p u ed e com entar es que los resu lta d o s son
coherentes, pues la potencia se conserva con un pequeño error. A la vista de
las tablas 5.IV.a-b, 5.V.a-b y 5.VI.a-b cabe decir que los resultados óptimos se
consiguen para la sim ulación m as corta, en el sentido de que se obtienen
m as ráp id am en te y no presentan una g ran variación con respecto a los
obtenidos con una sim ulación m as prolongada. Es arriesgado efectúar la
259
TDF prácticam ente después de llegada la señal, sin esperar al estacionario,
(el cual se alcanza relativam ente pronto), por eso se ha cálculado tam bién
durante tres y quince periodos, para com parar resultados. Por otro lado una
sim ulación m ucho m as prolongada introduce m ucho error en el cálculo de
la TDF, p o r eso los resultados para quince periodos tienen u n error en la
conservación de la potencia de alrededor del 5%. H ay que añadir que los
puntos de m uestreo para las puertas 1 y 2 no están sim étricam ente situados
respecto de la discontinuidad, por esa razón el defase de la señal en estos
puntos con excitación en la puerta 3 es ligeram ente distinto.
estos parám etros com parados con los obtenidos en las tablas 5.V.a-b,
arrojan u n error de u n 0.8% para IS n I, u n 0.2% para IS1 2 I, un 0.3% para
IS 1 3 I y u n 1.4% para IS3 3 I. Cabría añadir tam bién como conclusión, que
efectivam ente, como afirm a el autor, estas m edidas experim entales no son
dem asiado afortunadas [116].
260
em pleado en el cálculo de la TDF, la cual se efectuaba para cada una de las
frecuencias de las excitación para una com ponente del campo a lo largo de
las líneas que unían cada puerta con la discontinuidad. Los dibujos han sido
efectuados en u n Apollo 400t y u n HP 9000/375 haciendo uso de las librerías
gráficas Disspla.
261
5.6. Conclusiones.
El m étodo de Diferencias Finitas en el Dom inio del Tiempo se ha
aplicado al análisis de discontinuidades en tres dimensiones.
262
Por últim o se ha m odelizado una discontinuidad en guía coaxial
cuadrada. Se han calculado sus parám etros de dispersión en seis frecuencias
de la b an d a C, que se han com parado con los obtenidos experim ental y
teó ricam en te p a ra el caso cuasiestático. Los resu lta d o s tam b ién son
aceptables.
263
264
APENDICES
Apéndice 4.A
Condiciones de contorno absorbentes en esquinas al espacio libre
Ah
A
Ak
Ah
Figura 4.A.1
265
Para calcular el campo en el punto A, se aproxim a po r u n desarrollo
de Taylor de prim er orden,
,2 - ^ 2 ”
siendo Ak= ( ^ ).
266
Apéndice 5.A
Interfase de cambio m allado fino-mallado grueso
En los puntos comunes a las dos zonas (fina y gruesa) hay que llevar
cuidado en el uso de las diferencias centradas. H ay que m odificar el
algoritm o p ara que sea coherente a am bos lados. Suponiendo que estos
puntos están en xo=ioAh, que a la izquierda tenemos el m allado fino (Ax=Ah)
y a la derecha el grueso (Ax=aAh), podrem os aproximar,
df(xp) f(xQ-KxAh)-f(xQ-Ah)
031 Al f 1 + ° S
~
Ah( ^ )
1+oc
- En la superfice limítrofe entra ambas zonas, A x= A h(~^
A Ah
V3
267
tom ándose,
a Ah
cAt = -To
cAt' <
I 1 + 1~^ 1
\ (ocAh) Ah^+Ah"
que queda,
a V3 Ah
cAt' <
V20C+1 ^
a ^3 A
cAt = 1' At
\ 2 a +1
268
CONCLUSIONES Y PRINCIPALES
APORTACIONES
269
En este trabajo se estudia la form a de resolver estas dificultades. En
p articular y con m ayor pro fu n d id ad , se ha estudiado el caso de una guía
coaxial cuadrada, que luego se ha dirigido al análisis de una discontinuidad,
una T sim étrica construida en coaxial cuadrado. El motivo del estudio radica
tanto en su utilidad para la im plem entación de redes form adoras de haces,
como por la inexistencia de m odelos para la misma. A continuación y como
resum en del trabajo realizado, se enum eran las principales aportaciones y
conclusiones del mismo:
1.- Se ha realizado una pro fu n d a revisión de las bases m atem áticas del
algoritm o DFDT, que ha p erm itido deducir el factor de estabilidad del
algoritm o, establecer el m odo de m odelar la región problem a. D escribir la
dispersión de tipo num érico generada por el m étodo. Obtener el m odo de
separar la región problem a en dos zonas, una de campo total y otra de
cam po reflejado; y por últim o, diseñar e im plem entar las Condiciones de
Contorno A bsorbentes.
2.- Se ha desarrollado la teoría que justifica el análisis tem poral com binado
con la T ranform ada D iscreta de F ourier (TDF), para d e te rm in a r las
frecuencias de corte de los m odos en u n sistema guiador cerrado. Este se ha
tratado como un problem a resonancia transversal, perm itiendo u n análisis
en dos dim ensiones.
3.- Con el fin de desarrollar y verificar la m etodología que perm ita aplicar
DFDT al análisis de m odos, se ha aplicado a una guía rectangular y a una
guía circular, cuyas soluciones analíticas son perfectam ente conocidas. Se ha
analizado la convergencia tanto en función del m allado como del tiem po de
sim ulación, lo cual ha puesto de m anifiesto una escasa influencia del
tam año de la celda cuando el m allado se ajusta bien a la geom etría.
A u m en tan d o el tiem po de sim ulación se ha p o d id o discernir el error
asociado a discretización espacio tem poral de DFDT del error introducido
por la aplicación de la TDF. Por últim o, se ha m ostrado que los contornos
curvos no son una lim itación.
4.- Se han calculado las longitudes de onda de corte de los ocho prim eros
m odos TE y los dnco prim eros TM en una guía coaxial de sección cuadrada,
para distintas relaciones entre las dim ensiones de sus conductores interior y
270
exterior. Tam bién se ha estu d iad o un a guía coaxial con conductores
rectangulares. En ambos casos se ha m ostrado la distribución del campo
longitudinal para varios de estos modos.
5.- Se ha analizado una guía rectangular con una inserción metálica. Se han
calculado las frecuencias de corte de los cuatro prim eros m odos para varias
alturas del conductor interior, m ostrando su estructura en dos casos.
8.- Se calcula el coeficiente de transm isión de potencia en una codo recto. Las
CCA de prim er orden se m uestran adecuadas para una discontinuidad de
dos puertas. Los resultados son buenos con un m allado W /20, siendo W la
anchura de la guía.
11.- Con el fin de m odelizar situaciones en las que se pierde potencia por
radiación, se ha estudiado una guía de planos paralelos abierta. Este es un
problem a de radiación que se aborda desde el punto de vista circuital. Se
calcula el coeficiente de reflexión, se analiza la influencia de la distancia a las
271
CCA que describen el espacio libre y como interviene el tiem po de
sim u lació n .
272
Para la consecución de todos los puntos anteriores ha sido necesaria la
elaboración de una extensa librería de subrutinas en fortran. Estas librerías
p erm iten la m odelización de cualquier estructura de m icroondas en dos o
tre s d im e n sio n e s p a ra el a n á lisis e sp a c io -te m p o ra l del cam po
electrom agnético. La construcción de los program as ha sido enfocada a la
optim ización en el u so de recursos de ord en ad o r así como al uso de
procesadores vectoriales, contem plando en su estructuración la posibilidad
de crecim ien to asi com o de fu tu ra s m odificaciones. El g rad o de
vectorización alcanzado en el caso de la unión T en coaxial, conseguía que el
80% de tiem po de CPU consum ido fuese de procesador vectorial, porcentaje
que aum entaba al prescindir del cálculo de la TDF. En el análisis de modos
TE y TM la vectorización de los procesos de cálculo se ha acercado al 100%,
pero en geometrías rectangulares.
273
274
BIBLIOGRAFIA
[2]. DALY, P. "H ybrid-m ode analysis of m icrostrip by finite elem ent
m ethod," IEEE Trans.,1971,MTT-19, pp. 19-25.
[10]. HOEFER, W. J. R. and ROS, A. "Fin line param eters calculated w ith
the TLM -method," IEEE MTT-S Int. M icrowave Symp. Dig., pp. 341-
343, Apr.-May, 1979.
[11]. A K H TA R ZA D , S. a n d JO H N S, P. B., "T h re e -d im en sio n a l
transm ission-line m atrix Computer analysis of m icrostrip resonators,"
IEEE Trans.,1985, MTT-33, pp. 882-893.
[18]. SCHULZ, U. and PREGLA, R., "A new technique for the analysis of the
dispersión characteristics of p lanar w aveguides and its application to
m icrostrips w ith tuning septums," Radio Sci., 1981,16, pp. 1173-1178.
[19]. SHIH, Y. C., ITOH, T. and BUI, L. Q., "C om puter-aided design of
m illim eter-wave E-plane filters," IEEE Trans., 1983., MTT-31, pp. 135-
142.
[20]. ITOH, T., "Analysis of m icrostrip resonators," IEEE Trans., 1974, MTT-
22, pp. 946-952.
[21]. YEE, K.S., "N um erical solution of initial b o u n d ary valué problem s
.involving M axw ell's equations is isotropic m edia", IEEE Trans., 1966,
AP-14,pp.302-307.
276
[22]. TAFLOVE, A. and BRODWIN, M.E.:"Numerical solution of steady-state
electro m ag n etic sc ate rin g p ro b lem s u sin g the tim e d e p e n d e n t
Maxwell"s equations", IEEE Trans.,1975, MTT-23,pp.623-630.
[23]. TAYLOR, C.D., LAM, D.H. and SHUMPERT, T.H, "electrom agnetic
pulse scattering in tim e-varging in hom ogeneus m edia, IEEE Trans.
1969, AP-17.
[26]. TAFLOVE, A. :"Aplication of the finite difference tim e dom ain m ethod
to sinusoidal steady-state electrom agnetic penetration problems", IEEE
Trans.,1980, EMC-22, pp.191-202.
[27]. MUR, G. :"Absorbing boun d ary conditions for the finite difference
approxim ation of the tim e dom ain electrom agnetic field equations",
IEEE Trans.,1981, EMC-23, pp.377-382.
277
[30]. TAFLOVE, A. a n d UM ASHANKAR, K.R. : "A H y b rid M om ent
M ethod/Finite-D ifference Time Dom ain A pproach to Electrom agnetic
C oupling and A perture Penetration into Complex Geom etries", IEEE
Trans.,1982, AP-30, pp.617-627.
[33]. BRITT, C.L., : " Solution of electrom agnetic scattering problem s using
tim e dom ain techniques", IEEE Trans.,1989, AP-37, pp.1181-1192.
[35]. FURSE C.M., MATHUR, S.P. and GHANDI, O.P.,: "Im provem ents to
the Finite-Difference Tim e-D om ain m ethod for calculating the rad a r
cross section of a perfectly conducting target", IEEE Trans.,1990, MTT-38,
pp.919-927.
[36]. BORUP, T.D., SULLIVAN, D.M. and GHANDI, O.P.: "Com paration of
the FFT C onjúgate G radient m ethod and the Finite-Difference Time-
D om ain m ethod for the 2-D absortion problem ", IEEE T rans.,1987,
MTT-35, pp.383-395.
[38]. SULLIVAN, D., "T hree-dim ensional Computer- sim ulation in deep
reg io n a l h y p e rth e rm ia u sin g the fin ite-d ifferen ce T im e-D om ain
method", IEEE Trans., 1990, MTT-38, pp.204-211.
278
[39]. SULLIVAN, D., "M athem atical m ethods for treatm entn planing in
deep regional hypertherm ia'', IEEE Trans.,1991, MTT-39, pp. 864-874.
[45]. ZHAN G, X., FANG, J., MEI, K.K. and LIU, Y. :"Calculations of the
d isp ersiv e characteristics of m icrostrip by the tim e-dom ain finite
difference m ethod", IEEE Trans.,1988, MTT-36, pp.263-267.
[47]. ZHANG, X., and MEI, K.K. :"Time dom ain finite difference approach to
the calculation of the frequency dependent characteristics of m icrostrip
discontinuities", IEEE Trans.,1988, MTT-36, pp.1775-1787.
279
[48]. OLIVIER, J.C. a n d M cNAM ARA, D .A .,:"Finite-D ifference T im e-
D o m ain (FD-TD) a n aly sis of d isc o n tin u itie s in h o m o g e n e o u s
dispersive w aveguides", Electronics Letters, Vol.25, 1989, pp.1006-1007.
[49]. LIANG, G.C., LIU, Y.W. and MEI, K.K.,:,,Full wave analysis of coplanar
w av eg u id e and slotline u sin g the Tim e-D om ain Finite-D ifference
m ethod"JEEE Trans.,1989, MTT-37, pp.1949-1957.
[53]. NAVARRO, A., N U Ñ EZ, M.J., and MARTIN, E.: "Finite difference
tim e d o m a in FFT m e th o d a p p lie d to a x ia lly s y m m e tric a l
electromagnetic resonant devices", IEE Proc. H,1990,137,(3),ppl93-196.
[54]. NAVARRO, A., NU Ñ EZ, M.J., and MARTIN, E.: "Study of TEo and
TMo M odes in dielectric resonators by a finite difference tim e-dom ain
m eth o d coupled w ith the discrete fourier transform ", IEEE Trans.,
1991,MTT-39,(l),pp.l4-17.
[55]. BI, Z.-Q., WU, K.-L., and LITVA, J., :"Appücatíon of the FD-TD m ethod
to the A nalysis of H -plane W aveguide D iscontinuities", Electronics
Letters,1990, vol.26, pp.1897-1898.
280
[57]. KOIKE, S., YOSHIDA, N. and FUKAI, I.,: "T ransient analysis of
m icrostrip gap in tree-dim ensional space". IEEE Trans., 1985, MTT-33,
pp.726-730.
[58]. KOIKE, S., YOSHIDA, N. and FUKAI, I.,: "Transient analysis of coupling
betw een Crossing lines in three-dim ensional space". IEEE Trans., 1987,
MTT-35, pp.67-71.
[59]. KOIKE, S., YOSHIDA, N . and FUKAI, I.,: "T ransient analysis of
m icrostrip line on anisotropic substrate in three-dim ensional space",
IEEE Trans., 1988, MTT-36, pp.34-43.
[63]. BLASCHAK, J.G. and KRIEGSMANN, G.A. : "A Com parative Study of
Absorbing Boundary Conditions", J. Comput. Phys, 1988, 77, pp.109-139.
[6 6 ]. KIM, I.S. and HOEFER, W.J.R., :"A local m esh refinem ent algorithm for
the Time D om ain-Finite D ifference m eth o d u sin g M axw elLs curl
equations", IEEE Trans.,1990, MTT-38, pp.812-815.
281
[67].BERGER, M.J. an d OLIGER, J.R.,: "A daptative m esh refinem ent for
hyperbolic p ardal differentíal equation", J. Com put. Phys., 1984, vol.53,
pp.484-512.
[69]. RILEY, D.J. and TURNER, D., :"H ibrid thin-slot algorithm for the
an aly sis of n a rro w a p e rtu re s in Finite-D ifference T im e-D om ain
calculations", IEEE Trans.,1990, AP-38, pp.1943-1950.
[70]. ZIVANOVIC, S.S., YEE, K.S. and MEI, K.K., :"A subgridding m ethod for
th e T im e-D om ain F inite-D ifference m e th o d to solve M axw ell"s
equations", IEEE Trans.,1991, MTT-39, pp.471-479.
[74]. RALSTON, A., Introducción al Análisis Numérico , Lim usa-W iley S.A.,
México, 1970.
282
[79]. COLLIN, R.E., Field theory of Guided Waves : M cGraw-Hill 1966.
[85]. BRIGHAM, E.O., The Fast Fourier Transform and its Applications,
Prentice-Hall, N ew Jersey, 1988.
[8 8 ]. GRUNER, L., "H igher order m odes in square coaxial lines", IEEE
Trans.,1983, MTT-31,pp.770-772.
283
[91]. BRACKELM ANN, W. et al, "Die G renzfrequenzen von h o h eren
Eigenwellen in Streifenleitungen", Arch. Elek. U bertragung, 21, pp.112-
120.
[94]. GRUNER, L., "E stim ating rectangular coax cutoff", M icrow ave J.,
1979,22,pp.88-92.
[97]. PIN CHUP, A.R., CROWLEY, C.W. and SILVESTER, P.P., "Spectrally
correct finite elem ent operators for electrom agnetic vector fields",
J.Appl. Phys.,1988, 63, pp.3025-3027.
[98]. JOHNS, P.B.,: "On the relationship betw een TLM and Finite-Difference
m ethods for MaxwelTs equations", IEEE Trans.,1987, MTT-35, p p .60-61.
[99]. GWAREK, W.K.,: Com m ents on "On the relationship betw en TLM and
Finite-Difference m ethods for M axwell equations", IEEE Trans., 1987,
MTT-35, pp.872-873.
284
[102]. KOSHIBA, M. an d SUZUKI, M.: "A plication of boundary-elem ent
m eth o d to w av eg u id e d iscontinuities", IEEE T ran s.,1986, MTT-
34,pp.301-307.
[104]. METCALF, W.S., "Why not use rectangular coax?", Micro waves, 1968,
April, pp.52-56.
[105]. RIBLET, J.R., "The exact dim ensions of a family of rectangular coaxial
lines w ith given impedance", IEEE Trans.,1972, MTT-20, pp.538-541.
[106]. BOWMAN, F., "Notes on tw o-dim ensional electric field problem s",
Proc. London M ath. Soc.,1935, 39, pp.211-214.
[107]. ITOH, T., N um erical Techniques For M icrow ave A n d M illim eter-
Wave Passive Structures , John W iley & Sons, USA, 1989.
[110]. RIBLET, J.R., "Upper limit on the error of an im proved approxim ation
for the characteristic im pedance of rectangular coaxial line", IEEE
Trans.,1980, MTT-28, p p .666-667.
[113]. LEONG, M.S., KOOI, P.S. and CHANDRA, :"Radiatión from a flanged
parallel-plate w aveguide: solution by m om ent m ethod w ith inclusión
of edge condition", Proc IEE, Pt. H, 1988,135, pp.249-255.
285
[114]. LEE, S.-W., an d GRUN, L., :"Radiation from a flanged parallel-plate
w aveguide: C om parison of solutíons",IEEE Trans., 1982, AP-30, pp.147-
148.
286
|C/.'
ÜNIVERS1TÁT DE VALtNClA FACULTAD DE -&NCISS Rí|
V a le n c ia a 'Y a / d e 19.