PFC 1955
PFC 1955
PFC 1955
Resumen
En este proyecto fin de carrera se pretende estudiar la técnica de las imágenes espaciales
para el cálculo de funciones de Green en cavidades cilíndricas, optimizarla y extender el método
a cavidades de geometría arbitraría. El cálculo de las funciones de Green y sus potenciales
asociados son fundamentales para el análisis electromagnético de múltiples estructuras y
circuitos mediante la técnica de la ecuación integral.
-II-
-III-
A Alberto
-IV-
-V-
AGRADECIMIENTOS
Finalmente, un recuerdo para todos mis amigos, pidiendo perdón por el tiempo
que he faltado y no he podido estar con vosotros. Aun así habéis estado siempre ahí,
muchas gracias. En especial agradecerte Mónica tu cariño y apoyo (y dotes matemáticas
:P) durante todo este tiempo.
-VI-
-VII-
Índice de Contenidos
1. INTRODUCCIÓN ........................................................................................... 1
1.1 Descripción del proyecto......................................................................................................................1
1.1.1 Planteamiento inicial del proyecto..................................................................................................1
1.1.2 Objetivos.........................................................................................................................................2
1.1.3 Estructura del proyecto ...................................................................................................................2
1.1.4 Aportaciones novedosas .................................................................................................................3
1.1.5 Software desarrollado .....................................................................................................................4
-VIII-
3.2.3.2 Resultados obtenidos ............................................................................................................65
3.2.4. Método de la variación dinámica del radio imagen .....................................................................67
3.2.4.1 Desarrollo teórico..................................................................................................................67
3.2.4.1 Resultados obtenidos ............................................................................................................69
3.2.5 Optimización por descenso de gradiente de las cargas .................................................................77
3.2.5.1 Desarrollo teórico..................................................................................................................77
3.2.5.2 Resultados obtenidos ............................................................................................................83
3.2.6 Optimización por descenso de gradiente de las posiciones ..........................................................88
3.2.6.1 Desarrollo teórico..................................................................................................................88
3.2.6.2 Resultados obtenidos ............................................................................................................95
3.2.7 Empleo de métodos de forma conjunta y comparación ................................................................99
3.2.8 Conclusiones...............................................................................................................................104
-IX-
5.2 Análisis de cavidades eléctricamente cortas...................................................................................176
5.2.1 Transformada de Sommerfeld ....................................................................................................176
5.2.2 Descripción de la cavidad ...........................................................................................................177
5.2.3 Análisis de cavidades circulares .................................................................................................179
5.2.3 Análisis de cavidades arbitrarias.................................................................................................185
-X-
-XI-
Índice de Figuras
-XII-
yy
Figura 2.37 Representación potencial vector G A . Fuente en (0.7 Ȝ, 0.7 Ȝ) y 20 dipolos imágenes. ............. 34
Figura 2.38 Situación una carga imagen respecto al cilindro......................................................................... 36
Figura 2.39 Dipolo fuente en el eje ‘x’ y sus dos dipolos imagen asociados a un punto tangente .................. 38
Figura 2.40 Potencial escalar magnético normalizado en el eje x=0.............................................................. 41
xx
Figura 2.41 Potencial vector eléctrico G F normalizado en el eje x=0 .......................................................... 42
xx
Figura 2.42 Potencial vector eléctrico G F normalizado en el eje y=0 ......................................................... 42
yx
Figura 2.43 Potencial vector eléctrico G F normalizado evaluado con 20 imágenes .................................... 42
Figura 3.1 Representación 15 cargas imágenes alrededor del cilindro, fuente en origen coordenadas ......... 50
Figura 3.2 Representación del potencial escalar en contorno del cilindro en función del ángulo.
Empleo de 15 imágenes, fuente en origen de coordenadas .............................................................................. 50
Figura 3.3 Representación de situación y coste del potencial escalar en un cilindro, empleo de 15
imágenes y fuente situada en (0.7 Ȝ,0.7 Ȝ) ........................................................................................................ 51
Figura 3.4 Distancia entre puntos tangentes en la aproximación del trapecio................................................ 55
Figura 3.5 Aproximaciones integrales cerradas: solape de información ........................................................ 58
Figura 3.6 Características de la aproximación de Simpson’s ......................................................................... 58
Figura 3.7 Aproximaciones integrales semi-abiertas: no se solapa la información........................................ 59
Figura 3.8 Situación de los puntos intermedios en una aproximación semi-abierta ....................................... 60
Figura 3.9 Detalle de los puntos intermedios tomados en un arco del cilindro al realizar una
aproximación por Trapecio general................................................................................................................. 61
Figura 3.10 Representación de situación y coste del potencial escalar en un cilindro, empleo de 14
imágenes y fuente situada en (0.7 Ȝ,0.7 Ȝ) ........................................................................................................ 63
Figura 3.11 Representación de situación y coste del potencial escalar en un cilindro tras ubicar un
punto tangente en un pico de potencial, empleo de 15 imágenes y fuente situada en (0.7 Ȝ,0.7 Ȝ) .................. 64
Figura 3.12 Detalle de la ubicación de un nuevo punto de condición de potencial nulo donde se
encontraba anteriormente un pico de potencial. .............................................................................................. 64
Figura 3.13 Detalle de la distribución del potencial, 25 cargas imágenes habiendo evolucionado su
colocación en picos de potencial desde 19 cargas. Fuente en (0.7 Ȝ, 0.7Ȝ) ..................................................... 66
Figura 3.14 Distribución del potencial a lo largo del contorno del cilindro con 15 imágenes.
a) Situando las cargas en picos de potenciales y evolucionando desde una situación inicial de 10
imágenes. b) Empleando Point-Matching. ....................................................................................................... 66
Figura 3.15 Situación 20 cargas imagen alrededor cilindro, fuente en (0.5Ȝ ,0.5Ȝ)........................................ 67
Figura 3.16 Situación 20 cargas imagen , fuente ubicada en (0.5Ȝ, 0.5Ȝ) ...................................................... 68
Figura 3.17 Situación 20 cargas imagen ubicación circular. Fuente en (0Ȝ,0Ȝ) ............................................. 69
Figura 3.18 Coste en función del radio (0-8000). 20 imágenes, fuente en (0Ȝ,0Ȝ). ......................................... 69
Figura 3.19 Coste en función del radio (0-11). 20 imágenes, fuente en (0Ȝ,0Ȝ). ............................................. 70
Figura 3.20 Coste en función del radio (2-7). 20 imágenes, fuente en (0Ȝ,0Ȝ). ............................................... 70
Figura 3.21 Coste en función del radio, paso grueso y fino. 20 imágenes y fuente en (0Ȝ,0Ȝ) ........................ 71
Figura 3.22 Variación del coste con el radio variable (0-80). 10 cargas imágenes, fuente en (0Ȝ,0Ȝ) ........... 72
Figura 3.23 Comparativa costes método original y variación de radio .......................................................... 72
Figura 3.24 Potencial escalar en el contorno, métodos original y variación de radio. Fuente situada
en el origen y 20 imágenes ............................................................................................................................... 73
Figura 3.25 Variación del coste con el radio, 5 imágenes y fuente en el origen. ............................................ 74
Figura 3.26 Coste en función del radio. 10 cargas imágenes y fuente en (0.2Ȝ ,0.2Ȝ)..................................... 75
Figura 3.27 Situación con 10 imágenes y fuente en (0.2Ȝ ,0.2Ȝ) ...................................................................... 83
Figura 3.28 Representación Potencial escalar en función del ángulo. 10 imágenes y fuente en origen ......... 83
Figura 3.29 Evolución del Coste en función del número de iteración, paso dinámico. Empleo de 10
imágenes y fuente situada en la posición (0.2Ȝ,0.2Ȝ) ....................................................................................... 84
Figura 3.30 Evolución del Coste en función del número de iteración, paso fijo. Empleo de 10 imágenes
y fuente situada en la posición (0.2Ȝ ,0.2Ȝ) ...................................................................................................... 84
Figura 3.31 Representación Potencial escalar en función del ángulo tras optimizar el valor complejo
de las cargas. Empleo de 10 imágenes y fuente situada en el origen............................................................... 85
Figura 3.32 Situación con 10 imágenes y fuente en (0.5Ȝ ,0.5Ȝ) ...................................................................... 95
Figura 3.33 Potencial escalar en función del ángulo. 10 imágenes y fuente en (0.5Ȝ ,0.5Ȝ) ........................... 95
Figura 3.34 Representación de la evolución del coste y de las posiciones de las cargas imagen al
realizar una optimización de las posiciones. 10 imágenes y fuente en (0.5Ȝ ,0.5Ȝ).......................................... 96
Figura 3.35 Potencial escalar en función del ángulo tras optimizar las posiciones de las cargas
imagen. 10 imágenes y fuente en (0.5Ȝ ,0.5Ȝ) ................................................................................................... 97
Figura 3.36 Representación 20 cargas imágenes alrededor del cilindro, fuente en origen coordenadas ..... 107
-XIII-
Figura 3.37 Valor de las condiciones del potencial vector en el contorno del cilindro, 20 imágenes,
fuente situada en el origen. ............................................................................................................................ 108
Figura 3.38 Distribución del error del potencial vector en el contorno del cilindro, 20 imágenes,
fuente situada en el origen. ............................................................................................................................ 108
Figura 3.39 Representación de situación y coste del potencial escalar en un cilindro, empleo de 20
imágenes y fuente situada en (0.6 Ȝ,0.6 Ȝ) ...................................................................................................... 109
Figura 3.40 Detalle del Coste del potencial vector para cada condición. 20 imágenes y fuente situada
en (0.6 Ȝ,0.6 Ȝ) ................................................................................................................................................ 109
Figura 3.41 Diferencias en la ubicación de dipolos imágenes. a) Método general. b) Variación
dinámica del radio. Empleo de 20 dipolos imagen y fuente situada en (0.6Ȝ ,0.6Ȝ)....................................... 112
Figura 3.42 Variación del coste del potencial vector en función del radio. Empleo de 20 dipolos
imagen y fuente situada en (0.6Ȝ ,0.6Ȝ) .......................................................................................................... 113
Figura 3.43 Variación del coste del potencial vector en función del radio en torno al mínimo. Empleo
de 20 dipolos imagen y fuente situada en (0.6Ȝ ,0.6Ȝ).................................................................................... 114
Figura 3.44 Comparación de la distribución de la condición del potencial vector en el contorno de un
cilindro. a) Optimizando el radio. b) Método general. Empleo de 20 dipolos imagen y fuente situada
en (0.6Ȝ ,0.6Ȝ) ................................................................................................................................................. 114
Figura 3.45 Ejemplo de convergencia del coste del potencial vector en el infinito. 15 imágenes, fuente
en el origen de coordenadas........................................................................................................................... 115
Figura 3.46 Variación de la condición del potencial vector en el contorno del cilindro en función del
radio, existiendo 2 mínimos. 15 cargas imágenes y fuente en (0.4Ȝ ,0.0Ȝ)..................................................... 116
Figura 4.1 Dipolo unitario situado en el interior de una cavidad cuadrada ................................................. 122
Figura 4.2 Ubicación de una condición de contorno Gv en la pared de una guía cuadrada ........................ 123
Figura 4.3 Tramo oblicuo y sus ejes de coordenadas .................................................................................... 124
Figura 4.4 Ubicación de una condición de contorno del potencial vector .................................................... 125
Figura 4.5 Estructura cilíndrica, de radio a O , analizada con 20 imágenes y fuente situada en el
centro. Perspectivas de formulación cartesiana y formulación cilíndrica. .................................................... 130
Figura 4.6 Comparación del potencial escalar en un cilindro bajo las formulaciones cartesiana y
cilíndrica. 20 imágenes, fuente en el origen................................................................................................... 131
Figura 4.7 Comparación del potencial vector en un cilindro bajo las formulaciones cartesiana y
cilíndrica. 20 imágenes, fuente en el origen................................................................................................... 132
Figura 4.8 Ejemplos varios de estructuras arbitrarias creadas..................................................................... 135
Figura 4.9 Fichero descriptor de una guía ridge con sus lados redondeados ............................................... 136
Figura 4.10 Ejemplo de guía ridge con lados redondeados........................................................................... 136
Figura 4.11 Distribución uniforme de 15 puntos tangentes en una guía ridge “redondeada” ..................... 139
Figura 4.12 Distribución uniforme de 15 puntos tangentes en una guía triangular..................................... 139
Figura 4.13 Ubicación de 3,4 y 7 puntos tangentes en una guía triangular, empleando una distribución
de los mismos por vértices y mitades de segmentos. ...................................................................................... 140
Figura 4.14 Ubicación de 32 puntos tangentes en una guía ridge redondeada, empleando una
distribución de los mismos por vértices y mitades de segmentos. .................................................................. 140
Figura 4.15 Ubicación de puntos tangentes en guía ridge redondeada en función de la distancia media
de cada tramo a la fuente ............................................................................................................................... 141
Figura 4.16 Ubicación de 20 puntos de forma uniforme y saltando los vértices, sobre una guía
rectangular. .................................................................................................................................................... 142
Figura 4.17 Ubicación de imágenes de forma circular alrededor de una guía ridge redondeada y una
guía rombo. .................................................................................................................................................... 142
Figura 4.18 Ubicación de imágenes mediante posicionamiento especular. Ejemplos con una guía
cuadrada de bordes redondeados y una guía triangular. .............................................................................. 143
Figura 4.19 Ubicación de imágenes en tramos curvos .................................................................................. 144
Figura 4.20 Combinación de técnicas de ubicación de imágenes: especular y recta tangente auxiliar
para tramos curvos. Ejemplo sobre guía cuadrada con 2 lados redondeados............................................... 144
Figura 4.21 Ubicación de imágenes siguiendo el contorno de las guías “ridge” y “estrella” ..................... 144
Figura 4.22 Reconocimiento de tramos rectos situados entre un punto cualquiera y el exterior de la
estructura. Ejemplo con una guía cuadrada. ................................................................................................. 145
Figura 4.23 Reconocimiento de tramos oblicuos situados entre un punto cualquiera y el exterior de la
estructura. Ejemplo con una guía triangular. ................................................................................................ 146
Figura 4.24 Ejemplo de reconocimiento del interior de estructuras arbitrarias. Visualización en 3D. ........ 147
Figura 4.25 Potenciales para una guía cuadrada de 2x2m. 300Mhz. 40 imágenes....................................... 148
Figura 4.26 Potenciales para una guía triangular. 300Mhz. 30 imágenes.................................................... 150
Figura 4.27 Potenciales para una guía rombo. 300Mhz. 30 imágenes.......................................................... 151
-XIV-
Figura 4.28 Potenciales para una guía ridge. 300Mhz. 60 imágenes............................................................ 152
Figura 4.29 Potenciales para una guía estrella. 64 imágenes. ...................................................................... 153
Figura 4.30 Representación del potencial escalar en el contorno de una guía cuadrada. 40 imágenes,
fuente en el centro de la guía. ........................................................................................................................ 155
Figura 4.31 Representación del potencial escalar en el contorno de una guía ridge. 60 imágenes,
fuente en el centro de la guía. ........................................................................................................................ 156
Figura 4.32 Representación del potencial escalar en el interior de una guía ridge. 60 imágenes
situadas a 0.25Ȝ siguiendo el contorno de la guía. ........................................................................................ 156
Figura 4.33 Representación del Coste del potencial vector en una guía rombo, con 40 imágenes. 300
Mhz, fuente situada en el centro de la guía. ................................................................................................... 159
Figura 4.34 Representación del Coste del potencial vector en una guía cuadrada con dos lados
redondeados, con 64 imágenes. 300 Mhz, fuente situada en el centro de la guía.......................................... 160
Figura 4.35 Coste de los potenciales escalar y vector en una guía triangular, 30 imágenes ubicadas de
forma circular sin optimizar su radio. Fuente en el origen de coordenadas. ................................................ 161
Figura 4.36 Evolución del coste en función del radio del círculo que ubica a las 30 imágenes. Guía
triangular. ...................................................................................................................................................... 161
Figura 4.37 Ejemplos de ubicación de cargas mediante el método de la variación de la distancia
especular. Casos de guía ridge (40 imágenes) y guía rombo (30 imágenes) ................................................. 162
Figura 4.38 Ubicación de 64 imágenes en una guía “estrella” mediante un posicionamiento especular
sin optimizar. Representación del coste escalar en su contorno. ................................................................... 163
Figura 4.39 Situación final y representación de Costes en el contorno para potencia escalar y vector.
Guía estrella, 64 imágenes. ............................................................................................................................ 163
Figura 4.40 Situación de partida al analizar una guía rectangular, con 20 imágenes situadas a 0.5Ȝ de
la estructura. Representación del coste de los potenciales escalar en el contorno........................................ 165
Figura 4.41 Evolución del coste del potencial escalar al evaluar una guía rectangular, 20 imágenes......... 166
Figura 4.42 Situación de partida y fin de las 20 imágenes empleadas para analizar una guía ridge. .......... 168
Figura 4.43 Evolución del coste en función del nº de iteración. Optimización de las posiciones de 25
imágenes alrededor de una guía ridge. .......................................................................................................... 168
Figura 4.44 Análisis de una guía cuadrada, 40 imágenes. Fuente en el centro, 15 Ghz. .............................. 171
Figura 4.45 Potenciales normalizados en un corte a 16.5m a lo largo del eje ‘x’ desde el borde
inferior en una guía cuadrada de 66x66mm. Comparativa de método imágenes y método analítico
descrito en [10]. ............................................................................................................................................. 171
Figura 5.1 Numeración de capas y fuentes .................................................................................................... 178
Figura 5.2 Cilindro a analizar. Radio 1 metro, altura 0.5 metros. 40 imágenes situadas de forma
uniforme alrededor de la estructura. Fuente en el centro, a 0.25 metros de altura....................................... 182
Figura 5.3 Representación de potenciales en el interior de una guía cilíndrica, 1 metro de radio, 0.5
de alto, 40 imágenes alrededor. Fuente en el centro, a 0.25metros. 300Mhz. ............................................... 183
Figura 5.4 Representación de potenciales en el interior de una guía cilíndrica, 1 metro de radio, 0.5
de alto, 40 imágenes alrededor. Fuente en el centro, a 0.25metros. 600Mhz. ............................................... 184
Figura 5.5 Corte de los potenciales en el eje X en función de distintas distancias de ubicación de las
50 imágenes empleadas al analizar una guía circular de radio 1metro, 0.5 metros de alto. 300 Mhz.......... 184
Figura 5.6 Relación de ángulos en la formulación cartesiana con transformada Sommerfeld ..................... 186
Figura 4.7 Situación al analizar una guía cuadrada 66x66x30 mm. 50 imágenes, situadas siguiendo su
contorno situadas a 15 mm............................................................................................................................. 187
Figura 5.8 Potenciales en el interior de una guía cuadrada 66x66x30 mm. 50 imágenes, situadas
siguiendo su contorno situadas a 15 mm........................................................................................................ 188
Figura 5.9 Ubicación del corte en el eje X en la guía cuadrada para la visualización de los potenciales.... 188
Figura 5.10 Comparación de los potenciales en una guía cuadrada de 66x66x30mm en un corte en el
eje X a 16.5mm del eje Y, 50 imágenes, distintas posiciones de las mismas. Fuente en el centro. Se
muestra el resultado analítico esperado según [10]. 10 Ghz......................................................................... 189
Figura 5.11 Situación a analizar: guía cuadrada 66x66x30mm, 40 imágenes ubicadas a 33mm del
contorno. Fuente descentrada, situada en (-10.7,0)mm del plano X-Y. ......................................................... 190
Figura 5.12 Comparación de los potenciales de una guía cuadrada 66x66x30mm en un corte en el eje
X a 16.5mm del eje Y, 40 imágenes, ubicadas a 33mm del contorno. Fuente en el centro.Se muestra el
resultado analítico esperado según [10]. 4 Ghz. ........................................................................................... 190
Figura 5.13 Ubicación de imágenes 40 imágenes alrededor de una guía cuadrada 66x66x30mm en
función aumentada de su distancia a la fuente (descentrada)........................................................................ 191
Figura 5.14 Representación de los potenciales de una guía cuadrada 66x66x30 en un corte, al
emplear 40 imágenes situadas en función aumentada de su distancia a la fuente ponderada. 4 Ghz. .......... 191
-XV-
Figura 5.15 Comparación de los potenciales en un corte de una guía cuadrada de 66x66x30mm en
función del número de imágenes empleadas, situadas a 33mm de la guía. 10 Ghz. ...................................... 192
Figura 5.16 Comparación de los potenciales de una guía cuadrada 66x66x30 en un corte, al emplear
un anillo de 40 imágenes variando su distancia a la guía, con el resultado analítico de [10]. 15 Ghz......... 193
Figura 5.17 Situación planteada al analizar con 3 anillos una guía cuadrada ............................................. 194
Figura 5.18 Nomenclatura usada para la ubicación de imágenes y puntos tangentes en análisis
multianillo ...................................................................................................................................................... 195
Figura 5.19 Comparación de los potenciales de una guía cuadrada 66x66x30 en un corte, al emplear
tres anillos de 40 imágenes cada uno, variando su distancia a la guía , con el resultado analítico de
[10]. 15 Ghz. .................................................................................................................................................. 197
Figura 5.20 Comparación de los potenciales de una guía cuadrada 66x66x30 en un corte distinto a la
altura de la fuente, al emplear tres anillos de 40 imágenes cada uno situados a 33mm de la cavidad ,
con el resultado analítico de [10]. 15 Ghz. .................................................................................................... 198
Figura 5.21 Potenciales a media altura de una guía cuadrada de 66x66x30mm, 3 anillos de 50
imágenes cada uno, ubicados a 33mm de la cavidad. Frecuencia de 15Ghz................................................. 199
Figura 5.22 Comparación de los potenciales en un corte del eje X, en una guía cuadrada de 66x66x30
con un corte cuadrado en su esquina superior derecha de 18x18mm. 80 imágenes. 10Ghz. ........................ 201
Figura 5.23 Ubicación de 3 y 5 anillos de 80 imágenes alrededor de una guía cuadrada de 66x66x30
mm que presenta un corte de 18x18mm en su esquina superior derecha....................................................... 202
Figura 5.24 Potencial escalar en un corte para 1,3 y 5 anillos de 80 imágenes, en función de la
distancia de los mismos a la cavidad cuadrada de 66x66x30mm con un corte de 18x18mm. ....................... 202
Figura 5.25 Potencial vector en un corte para 1,3 y 5 anillos de 80 imágenes, en función de la
distancia de los mismos a la cavidad dual cuadrada de 66x66x30mm con un corte de 18x18mm. ............... 203
Figura 5.26 Potencial escalar en un corte para 1,3 y 5 anillos de 80 imágenes, en función de la
distancia cercana de los mismos a la cavidad cuadrada de 66x66x30mm con un corte de 18x18mm........... 204
Figura 5.27 Potencial vector en un corte para 1,3 y 5 anillos de 80 imágenes, en función de la
distancia cercana de los mismos a la cavidad cuadrada de 66x66x30mm con un corte de 18x18mm........... 205
Figura 5.28 Potencial escalar en un corte para 1,3 y 5 anillos en función del número de imágenes por
anillo en una cavidad cuadrada de 66x66x30mm con un corte de 18x18mm. ............................................... 206
Figura 5.29 Potencial vector en un corte para 1,3 y 5 anillos en función del número de imágenes por
anillo en una cavidad cuadrada de 66x66x30mm con un corte de 18x18mm. ............................................... 207
Figura 5.30 Visualización de “picos” de error en el potencial escalar de una guía cuadrada de
66x66x30mm, con un corte en una esquina de 18x18mm. 6.122Ghz. 60 imágenes. ...................................... 207
Figura 5.31 Visualización de las matrices del potencial escalar para las guías cuadrada y dual, al
emplear 80 imágenes...................................................................................................................................... 208
Figura 5.32 Ejemplos de guías rombos, con distintos ángulos internos. ....................................................... 208
Figura 5.33 Potenciales en el corte X central para una guía rombo con ángulos internos de 90º................ 209
Figura 5.34 Potenciales en el corte X central para una guía rombo con ángulos internos de 60º................ 210
Figura 5.35 Potenciales en el corte X central para una guía rombo con ángulos internos de 45º................ 210
Figura 5.36 Potenciales en el corte X central para una guía rombo con ángulos internos de 30º................ 211
Figura 5.37 Potenciales en el eje X de una guía rombo con ángulos internos de 30º. Variación del
número de imágenes de un único anillo. Frecuencia de 10Ghz. .................................................................... 212
Figura 5.38 “Picos” de error al analizar una guía rombo con 2 ángulos internos con 30º a 12 Ghz.
Potencial vector GAyy.................................................................................................................................... 212
Figura 5.39 Situación de la fuente y un punto de observación en el interior de una guía cuadrada de
66x66x30mm, empleo de 60 imágenes situadas a 33mm................................................................................ 214
Figura 5.40 Evaluación de los potenciales en un punto fijo en el interior de una guía cuadrada de
66x66x30mm, al variar la frecuencia de 5.2 a 7.5Ghz................................................................................... 215
Figura 5.41 Potenciales en la primera resonancia de la cavidad cuadrada 66x66x30mm. Comparación
con el campo eléctrico proporcionado por HFSS©. ...................................................................................... 217
Figura 5.42 Potenciales en la segunda resonancia de la cavidad cuadrada 66x66x30mm.
Comparación con el campo eléctrico proporcionado por HFSS©. ............................................................... 218
Figura 6.1 Situación de la fuente y un punto de observación en el interior de una guía rectangular de
132x66x30mm, empleo de 60 imágenes situadas a 12mm.............................................................................. 224
Figura 6.2 Evaluación de los potenciales en un punto fijo en el interior de una guía rectangular de
132x66x30mm, al variar la frecuencia de 5.2 a 7.0 Ghz................................................................................ 225
Figura 6.3 Potenciales en la primera resonancia de la cavidad rectangular 132x66x30mm.
Comparación con el campo eléctrico proporcionado por HFSS©. ............................................................... 227
Figura 6.4 Potenciales en la segunda resonancia de la cavidad rectangular 132x66x30mm.
Comparación con el campo eléctrico proporcionado por HFSS©. ............................................................... 228
-XVI-
Figura 6.5 Potenciales en la tercera resonancia de la cavidad rectangular 132x66x30mm.
Comparación con el campo eléctrico proporcionado por HFSS©. ............................................................... 229
Figura 6.6 Situación de la fuente y un punto de observación en el interior de una guía triangular de
altura 30mm, empleo de 30 imágenes situadas a 12mm. ............................................................................... 230
Figura 6.7 Evaluación de los potenciales en un punto fijo en el interior de una guía triangular de
altura 30mm, al variar la frecuencia de 5 a 9.0 Ghz...................................................................................... 231
Figura 6.8 Evaluación de los potenciales cruzados en un punto fijo en el interior de una guía
triangular de altura 30mm, al variar la frecuencia de 5 a 9.0 Ghz. .............................................................. 231
Figura 6.9 Potenciales en la primera resonancia de una cavidad triangular con 30mm de altura.
Comparación con el campo eléctrico proporcionado por HFSS©. ............................................................... 234
Figura 6.10 Potenciales en la segunda resonancia de una cavidad triangular con 30mm de altura.
Comparación con el campo eléctrico proporcionado por HFSS©. ............................................................... 235
Figura 6.11 Situación de la fuente y un punto de observación en el interior de una guía rombo de
dimensiones 66x66x30mm y sus 4 ángulos internos a 90º. Empleo de 60 imágenes situadas a 12mm.......... 236
Figura 6.12 Evaluación de los potenciales en un punto fijo en el interior de una guía rombo de
dimensiones 66x66x30mm al variar la frecuencia de 5 a 15.0 Ghz. .............................................................. 237
Figura 6.13 Potenciales en la primera resonancia de una cavidad rombo de dimensiones
66x66x30mm. Comparación con el campo eléctrico proporcionado por HFSS©. ........................................ 238
Figura 6.14 Potenciales cruzados en la primera resonancia de una guía rombo de dimensiones
66x66x30mm................................................................................................................................................... 239
Figura 6.15 Situación de 3 anillos con 60 imágenes cada uno al analizar una guía rombo de
dimensiones 66x66x30mm con sus 4 ángulos internos a 90º.......................................................................... 239
Figura 6.16 Potenciales a la frecuencia de resonancia de 14.18Ghz de una cavidad rombo de
dimensiones 66x66x30mm. Comparación con el campo eléctrico proporcionado por HFSS©..................... 240
Figura 6.17 Situación de 1 anillo con 60 imágenes situados sobre una guía en forma de trapecio, de
altura 30mm. .................................................................................................................................................. 241
Figura 6.18 Evaluación de los potenciales en un punto fijo en el interior de una guía trapezoidal de de
altura 30mm al variar la frecuencia de 4.5 a 9.0 Ghz.................................................................................... 242
Figura 6.19 Potenciales a la frecuencia de resonancia de 7.36Ghz de una cavidad trapezoidal de
altura 30mm. Comparación con el campo eléctrico proporcionado por HFSS©. ......................................... 244
Figura 6.20 Potenciales a la frecuencia de resonancia de 7.826Ghz de una cavidad trapezoidal de
altura 30mm. Comparación con el campo eléctrico proporcionado por HFSS©. ......................................... 245
Figura 6.21 Potencial escalar a la frecuencia de resonancia de 8.94Ghz de una cavidad trapezoidal de
altura 30mm. Comparación con el campo eléctrico proporcionado por HFSS©. ......................................... 246
Figura 6.22 Situación de 1 anillo con 90 imágenes situados sobre una guía cuadrada de 66x66x30mm
con un corte de 18x18mm en su esquina superior derecha. ........................................................................... 247
Figura 6.23 Evaluación de los potenciales en un punto fijo en el interior de una guía cuadrada de
66x66x30mm con un corte de 18x18 en una esquina, al variar la frecuencia de 4.5 a 9.0 Ghz..................... 247
Figura 6.24 Potencial escalar en la primera resonancia en función del número de imágenes, guía
cuadrada 66x66x30mm con un corte de 18x18 en su esquina superior derecha. .......................................... 250
Figura 6.25 Campo eléctrico Ez obtenido en una guía dual en la primera resonancia, a 6.122Ghz. ............ 250
Figura 6.26 Potencial escalar a la frecuencia de resonancia de 7.125Ghz de una cavidad dual de
altura 30mm. Comparación con el campo eléctrico proporcionado por HFSS©. ......................................... 251
Figura 6.27 Potencial escalar a la frecuencia de resonancia de 7.629Ghz de una cavidad dual de
altura 30mm. Comparación con el campo eléctrico proporcionado por HFSS©. ......................................... 252
Figura 6.28 Potencial escalar a la frecuencia de resonancia de 8.952Ghz de una cavidad dual de
altura 30mm. Comparación con el campo eléctrico proporcionado por HFSS©. ......................................... 252
Figura 6.29 Situación de 1 anillo con 90 imágenes situados sobre una guía flange de altura 30mm. .......... 253
Figura 6.30 Representación de los potenciales en un corte del eje X en una guía flange a 7 Ghz. 1
anillo con 50 imágenes................................................................................................................................... 254
Figura 6.31 Representación de los potenciales en un corte del eje X en una guía flange a 7 Ghz. 3
anillos con 50 imágenes cada uno.................................................................................................................. 255
Figura 6.32 Representación de los potenciales en un corte del eje X en una guía flange a 7 Ghz. 5
anillos con 50 imágenes cada uno.................................................................................................................. 255
Figura 6.33 Representación de los potenciales en un corte del eje X en una guía flange a 7 Ghz. 5
anillos, variación del número de imágenes por anillo. .................................................................................. 256
Figura 6.34 Evaluación de los potenciales en un punto fijo en el interior de una guía flange de altura
30mm al variar la frecuencia de 5 a 9.0 Ghz. ................................................................................................ 257
Figura 6.35 Potenciales a la frecuencia de resonancia de 7.155Ghz de una cavidad flange de altura
30mm. Comparación con el campo eléctrico proporcionado por HFSS©..................................................... 259
-XVII-
Figura 6.36 Potenciales a la frecuencia de resonancia de 8.419Ghz de una cavidad flange de altura
30mm. Comparación con el campo eléctrico proporcionado por HFSS©..................................................... 260
Figura 6.37 Situación de 1 anillo con 150 imágenes situados sobre una guía estrella de altura 30mm. ...... 261
Figura 6.38 “Picos” de error en los vértices de una guía estrella al ser analizada con 70 imágenes.......... 261
Figura 6.39 Potenciales en ejes ‘x’ e ‘y’ al variar el número de imágenes en el análisis de una guía
estrella............................................................................................................................................................ 262
Figura 6.40 Evaluación de los potenciales en un punto fijo en el interior de una guía estrella de altura
30mm al variar la frecuencia de 5 a 9.0 Ghz. ................................................................................................ 263
Figura 6.41 Primer modo de resonancia en una guía estrella de altura 30mm. Frecuencia 5.723 Ghz ....... 264
Figura 6.42 Potencial escalar a la frecuencia de resonancia de 6.9358Ghz de una cavidad tipo
estrella de altura 30mm. Comparación con el campo eléctrico proporcionado por HFSS©. ....................... 265
Figura 6.43 Situación de 1 anillo con 90 imágenes situados sobre una guía ridge de altura 30mm............. 266
Figura 6.44 “Picos” de error en los vértices de una guía ridge al ser analizada con 70 imágenes ............. 266
Figura 6.45 Potenciales en el eje ‘x’ al variar el número de imágenes en el análisis de una guía ridge. .... 267
Figura 6.46 Evaluación de los potenciales en un punto fijo en el interior de una guía ridge de altura
30mm al variar la frecuencia de 7 a 10.5 Ghz. .............................................................................................. 268
Figura 6.47 Potencial escalar a la frecuencia de resonancia de 7.893Ghz de una cavidad tipo ridge de
altura 30mm. Comparación con el campo eléctrico proporcionado por HFSS©. ......................................... 269
Figura 6.48 Potencial escalar a la frecuencia de resonancia de 9.1127Ghz de una cavidad tipo ridge
de altura 30mm. Comparación con el campo eléctrico proporcionado por HFSS©. .................................... 270
Figura 6.49 Guía cuadrada multicapa y ubicación de 50 imágenes alrededor de la misma......................... 271
Figura 6.50 Potenciales en el eje ‘x’ al variar la distancia de 50 imágenes a una guía cuadrada con
dos dieléctricos en su interior. ....................................................................................................................... 271
Figura 6.51 Potenciales vector al realizar un barrido en un punto fijo entre 2 y 7 Ghz en una guía
cuadrada con 2 dieléctricos. Comparación con el método analítico descrito en [10]................................... 272
Figura 6.52 Comparación de los potenciales en una guía cuadrada de 66x66x30mm con un dieléctrico
de 15mm con H r 2 . Validación con los resultados obtenidos en [10] ...................................................... 273
Figura 6.53 Potencial escalar a 5.5Ghz en una guía cuadrada de 120x120mm, con dos capas de
12mm, una de ellas aire y la otra con permitividad de valor 2. ..................................................................... 273
-XVIII-
-XIX-
Índice de Tablas
-XX-
Tabla 4.3 Posibles valores del parámetro “Ubicación de imágenes” ........................................................... 173
Tabla 5.1 Frecuencias de resonancia de los primeros modos de una cavidad cuadrada de
66x66x30mm................................................................................................................................................... 214
Tabla 5.2 Precisión al obtener las frecuencias de corte de una guía cuadrada de 66x66x30mm.................. 215
Tabla 6.1 Frecuencias de resonancia de los primeros modos de una cavidad rectangular de
132x66x30mm................................................................................................................................................. 224
Tabla 6.2 Precisión al obtener las frecuencias de corte de una guía rectangular de 132x66x30mm. ........... 226
Tabla 6.3 Precisión al obtener las frecuencias de corte de una guía rectangular de 132x66x30mm. ........... 232
Tabla 6.4 Precisión al obtener las frecuencias de corte de una guía rombo de 66x66x30mm, con sus
cuadro ángulos internos a 90º........................................................................................................................ 237
Tabla 6.5 Precisión al obtener las frecuencias de corte de una guía trapecio de 30mm de altura................ 242
Tabla 6.6 Precisión al obtener las resonancias el potencial escalar en una guía dual ................................. 248
Tabla 6.7 Precisión al obtener las resonancias el potencial vector en una guía dual ................................... 248
Tabla 6.8 Precisión al obtener las frecuencias de corte de una guía flange de 30mm de altura................... 257
Tabla 6.9 Precisión al obtener las frecuencias de corte de una guía flange de 30mm de altura................... 263
Tabla 6.10 Precisión al obtener las frecuencias de corte de una guía flange de 30mm de altura................. 268
Tabla AI.1 Coste de distintos métodos de integración por pulso, en función del número de imágenes,
posición fuente en (0.0Ȝ ,0.0Ȝ)........................................................................................................................ 279
Tabla AI.2 Coste de distintos métodos de integración por pulso, en función del número de imágenes,
posición fuente en (0.2Ȝ ,0.2Ȝ)........................................................................................................................ 280
Tabla AI.3 Coste de distintos métodos de integración por pulso, en función del número de imágenes,
posición fuente en (0.4Ȝ ,0.0Ȝ)........................................................................................................................ 281
Tabla AI.4 Coste de distintos métodos de integración por pulso, en función del número de imágenes,
posición fuente en (0.4Ȝ ,0.4Ȝ)........................................................................................................................ 281
Tabla AI.5 Coste de distintos métodos de integración por pulso, en función del número de imágenes,
posición fuente en (0.5Ȝ ,0.0Ȝ)........................................................................................................................ 282
Tabla AI.6 Coste de distintos métodos de integración por pulso, en función del número de imágenes,
posición fuente en (0.7Ȝ ,0.7Ȝ)........................................................................................................................ 282
Tabla AII.1 Coste de la ubicación de cargas en picos de potencial, con distintas situaciones iniciales,
en función del número de imágenes y con posición fuente en (0.0Ȝ ,0.0Ȝ) ..................................................... 283
Tabla AII.2 Coste de la ubicación de cargas en picos de potencial, con distintas situaciones iniciales,
en función del número de imágenes y con posición fuente en (0.2Ȝ ,0.2Ȝ) ..................................................... 284
Tabla AII.4 Coste de la ubicación de cargas en picos de potencial, con distintas situaciones iniciales,
en función del número de imágenes y con posición fuente en (0.4Ȝ ,0.0Ȝ) ..................................................... 284
Tabla AII.5 Coste de la ubicación de cargas en picos de potencial, con distintas situaciones iniciales,
en función del número de imágenes y con posición fuente en (0.4Ȝ ,0.4Ȝ) ..................................................... 285
Tabla AII.6 Coste de la ubicación de cargas en picos de potencial, con distintas situaciones iniciales,
en función del número de imágenes y con posición fuente en (0.5Ȝ ,0.0Ȝ) ..................................................... 285
Tabla AII.7 Coste de la ubicación de cargas en picos de potencial, con distintas situaciones iniciales,
en función del número de imágenes y con posición fuente en (0.7Ȝ ,0.7Ȝ) ..................................................... 286
-XXI-
-XXII-
-XXIII-
Capítulo 1: Introducción
1. Introducción
1.1 Descripción del proyecto
1.1.1 Planteamiento inicial del proyecto
-1-
Capítulo 1: Introducción
1.1.2 Objetivos
En segundo lugar, se extenderá la técnica para poder analizar estructuras con una
geometría mucho más compleja. Así se propone el análisis de cavidades cuadradas,
rectangulares, triangulares, ridge, etc. Otra característica que presenta el método es la
posibilidad de cálculo de las frecuencias de resonancia de estas cavidades. En cuanto a
la validación de los resultados, se podrá realizar de forma numérica para las cavidades
cuadrada y rectangular (pues tienen solución analítica [8]) y se empleará el software
comercial HFSS© para la validación del resto de resultados.
-2-
Capítulo 1: Introducción
-3-
Capítulo 1: Introducción
-4-
Capítulo 1: Introducción
& wU
J (1.1)
wt
&
Además de la corriente de desplazamiento wD wt , ya puesta de manifiesto en la
ley de Ampère, pueden existir en cada punto corrientes de conducción (la denominada
como ley de Ohm):
-5-
Capítulo 1: Introducción
& &
J V E (1.2)
donde H r y P r son los valores relativos mientras que H 0 y P0 son los valores en el vacío,
siendo éstos:
10 9 F
H0
36S m (1.4)
7 H
P0 4S 10 m
En función del tipo de medio en el que nos encontremos, los valores relativos de
la permitividad y la permeabilidad serán reales, complejos, etc. En la tabla siguiente se
muestra el tipo de magnitud que adquieren en función del tipo de medio en el que se
realice el análisis:
-6-
Capítulo 1: Introducción
& &
nˆ ( D1 D2 ) U
& & &
nˆ u ( E1 E 2 ) M
& & (1.5)
nˆ ( B1 B 2 ) 0
& & &
nˆ u ( H 1 H 2 ) J
donde n̂ es el vector normal a la superficie.
Sin embargo, a pesar de los problemas que plantean, son la única manera de
resolver problemas electromagnéticos de forma general y se consigue una enorme
precisión en las soluciones encontradas si el planteamiento y método escogido es el
adecuado.
-7-
Capítulo 1: Introducción
-8-
Capítulo 1: Introducción
& &
Destacar la presencia en todas las ecuaciones anteriores de G (r r ' ) que es la
denominada función de Green del medio homogéneo, definida como respuesta al
impulso de la ecuación de Helmholzt:
& &
jk r r '
& & 1 e (1.10)
G (r r ' ) & &
4S r r '
-9-
Capítulo 1: Introducción
& &
§ jwP0 & & e jk r r ' ·
¨ ¸
& & ¨ 4S ³³S
J (r ' ) & &
4S r r '
dS
¸
nˆ u Eind (r ' ) nˆ u ¨ (1.12)
& & ¸
¨ j '
& & e
jk r r '
¸
¨ 4SwH ³³S S ³³S J (r ' ) 4S r& r& ' dS ¸
© 0 ¹
-10-
Capítulo 2: Formulación espacial en cavidades circulares
-11-
Capítulo 2: Formulación espacial en cavidades circulares
-12-
Capítulo 2: Formulación espacial en cavidades circulares
& &
Siendo Gv r , r ´ el potencial de una carga en condiciones de espacio libre:
-13-
Capítulo 2: Formulación espacial en cavidades circulares
& &
& & P 0 e jko r r ´
Gv r , r ´ & & (2.7)
4S r r ´
Hasta ahora hemos impuesto la condición de potencial nulo en un único punto
del contorno del cilindro, no obstante es deseable la presencia de más puntos con el fin
de poder satisfacer la condición en toda la pared del cilindro. Para ello, empleando la
misma estrategia que ahora, la extenderemos hasta un total de N puntos tangentes en el
contorno del cilindro. De esta forma, evaluaremos de forma numérica el valor complejo
de las cargas imágenes de forma que satisfagan de forma simultánea que el potencial en
todos y cada uno de los puntos tangentes es nulo.
Por ejemplo, la figura 2.4 muestra la situación de dos cargas imagen:
-14-
Capítulo 2: Formulación espacial en cavidades circulares
N
& & & & &
Gvcy (r ) GV (r , ro ' ) ¦ qk Gv (r , rk ' ) (2.10)
k 1
Figura 2.5 Dipolo fuente en el eje ‘x’ y sus dos dipolos imagen asociados a un punto tangente
Al igual que para el caso del potencial escalar, impondremos en primer lugar las
condiciones de contorno a cumplir en un punto discreto ‘i’ de la pared del cilindro.
En este caso, la componente tangencial del campo eléctrico (definido en la
ecuación 2.1) debe ser nulo:
&
E ptoi eˆt 0 (2.11)
donde êt es el vector tangente al punto ‘i’ de la estructura; esto conlleva a dos diferentes
&
condiciones, obtenidas tras analizar el campo eléctrico E definido por la ecuación 2.1
&
eˆ p u A 0 (2.12)
&
A 0 (2.13)
-15-
Capítulo 2: Formulación espacial en cavidades circulares
&
1) eˆ p u A 0
&
Que nos indica que la componente tangencial del vector potencial A
debe ser nula, lo que nos anularía el primer térnimo del campo eléctrico en ese
punto de la ecuación 2.1
Así, esta condición implica:
& & & & & &
eˆ p u ( I1xG Axx ri , r1´ eˆx I1y G Ayy ri , rk ´ eˆ y ) eˆ p u (G Axx ri , r0´ eˆx ) 0 (2.14)
&
Donde hemos tenido en cuenta la presencia de la carga fuente, en r0 .
-16-
Capítulo 2: Formulación espacial en cavidades circulares
Figura 2.6 Descomposición de un dipolo I1 en los dos ejes, con unos pesos determinados.
i 1,2,3...N
Estando definidas las constantes Cix,k y Ciy,k por la ecuación 2.20. Ambas
constantes pueden ser evaluadas para un medio general multicapa mediante el empleo
de la transformada de Sommerfeld [9] (como veremos posteriormente), aunque para el
espacio libre, ambas constantes coinciden a un valor dado.
Para obtener este valor, evaluaremos las partes en las que se compone la
ecuación 2.20 a):
& & & &
G Axx ri , rk ´ P 0 e jk r r ´ 1
0
en esta última ecuación, al estar definido el gradiente únicamente en la dirección êU , los
elementos de las direcciones êT y êM se anulan, por lo que desarrollando obtenemos:
-17-
Capítulo 2: Formulación espacial en cavidades circulares
eˆ U G ri , rk ´ eˆ U « ( 0 & & ) eˆ r »
xx
A eˆ U «( 0
& & & & 2 ) eˆ r
¬« wr 4S r r ´ ¼» 4S «¬ r r´ r r´
& & & &
P 0 e jko r r ´ ª jk 1 º P 0 e jko r r ´ § jk 1 ·
eˆ U «( & &0 & & 2 ) eˆ r » ¨ eˆ U eˆ r ( & &0 & & 2 )
4S «¬ r r ´ r r ´ »¼ 4S ¨ r r´ r r´
© ¹ (2.24)
& &
P0 e
jko r r ´ § ·
¨ >eˆ r& r&´@ ( jk 0 1
) ¸
& & & &
4S r r ´ ¨ U r r ´ r& r&´ 2 ¸
© ¹
De esta forma, ya tenemos todos los elementos para poder resolver la ecuación
2.21, de la que obtendremos el valor de los dipolos para poder recuperar el vector
potencial en el interior de la cavidad cilíndrica, mediante:
N
& & & & &
xx
G Acy (r ) G Axx r , r0 ´ ¦ I kx G Axx ri , rk ´
k 1
N
(2.26)
& & &
G yx
Acy (r ) ¦ I G ri , rk ´
k 1
k
y yx
A
-18-
Capítulo 2: Formulación espacial en cavidades circulares
2.1.2 Resultados
1. Posición (0Ȝ,0 Ȝ)
2. Posición (0.25 Ȝ,0.25 Ȝ)
3. Posición (0.625 Ȝ,0 Ȝ)
4. Posición (0.7 Ȝ,0.7 Ȝ)
-19-
Capítulo 2: Formulación espacial en cavidades circulares
Figura 2.8 Potencial escalar en el interior del cilindro. Fuente en (0 Ȝ,0 Ȝ) y 20 cargas imágenes
-20-
Capítulo 2: Formulación espacial en cavidades circulares
Figura 2.9 Potencial escalar en el contorno del cilindro. Fuente en (0 Ȝ,0 Ȝ) y 20 cargas imágenes
xx
Figura 2.10 Potencial vector G A en el interior del cilindro. Fuente en (0 Ȝ,0 Ȝ) y 20 dipolos imagen.
Comprobamos la orientación del potencial vector a lo largo del eje X, así como
un mayor valor en la posición en la que se encuentra la carga.
-21-
Capítulo 2: Formulación espacial en cavidades circulares
xy yx
Figura 2.11 Representación potencial vector G A y G A . Fuente en (0 Ȝ,0 Ȝ) y 20 dipolos imágenes.
yy
Figura 2.12 Representación potencial vector G A . Fuente en (0 Ȝ,0 Ȝ) y 20 dipolos imágenes.
-22-
Capítulo 2: Formulación espacial en cavidades circulares
Figura 2.13 Representación del potencial escalar en el contorno. Fuente en (0,0) y 5 cargas
imágenes
En este caso observamos como hemos impuestos cinco valores nulos del
potencial escalar en el contorno, pero en el resto de puntos de la pared el valor del
potencial aumenta de forma considerable (valores muy elevados, de hasta 3.3), por lo
que el potencial escalar obtenido con esta representación tendría un error muy
importante. No obstante el coste computacional en este caso sería muy reducido,
pudiéndose emplear para una primera aproximación.
Figura 2.14 Representación del potencial escalar en el contorno. Fuente en (0 Ȝ,0 Ȝ) y 30 cargas
imágenes
-23-
Capítulo 2: Formulación espacial en cavidades circulares
xy
Figura 2.15 Comparación potencial escalar y vectorial G A . Fuente en (0 Ȝ,0 Ȝ) y con 5 y 30
cargas/dipolos imagen
Para terminar con este caso, deberemos indicar que la situación de la fuente en la
posición central (0,0) del cilindro, ha provocado la simetría de todas las gráficas
mostradas hasta ahora. Con una variación de la posición de la fuente eléctrica, esta
simetría desaparecerá, obteniéndose resultados muy diferentes en función de la posición
de la misma.
-24-
Capítulo 2: Formulación espacial en cavidades circulares
En este caso, situaremos la fuente eléctrica en la posición (0.25 Ȝ, 0.25 Ȝ), lo que
provocará en primer lugar, una variación de la situación de las 20 cargas/dipolos imagen
alrededor de la cavidad circular, como podemos observar en la figura 2.16:
Figura 2.17 Potencial escalar en el interior del cilindro. Fuente en (0.25 Ȝ,0.25 Ȝ) y 20 cargas
imágenes
-25-
Capítulo 2: Formulación espacial en cavidades circulares
Figura 2.18 Representación del potencial escalar en el contorno. Fuente en (0.25 Ȝ ,0.25 Ȝ) y 30
cargas imágenes
Así vemos como en este caso, para el ángulo sobre el que se encuentra la fuente
obtenemos un potencial menor, aumentando éste en la dirección opuesta a la fuente.
Este comportamiento es particular, ya que de forma general suele producirse un
aumento del valor del Potencial escalar conforme la carga se acerca a la pared del
cilindro, como observaremos de forma detallada en un caso posterior.
A continuación, mostraremos el potencial vector G Axx para este caso en el
interior del cilindro:
xx
Figura 2.19 Potencial vector G A en el interior del cilindro. Carga en (0.25 Ȝ, 0.25 Ȝ) y 20 cargas imagen
Comprobamos como en este caso no se aprecia de una forma tan clara el efecto
de la orientación de la fuente a lo largo del eje ‘x’ como en el caso anterior,
obteniéndose además una distribución no simétrica.
-26-
Capítulo 2: Formulación espacial en cavidades circulares
Mostraremos en las figuras 2.20 y 2.21 los potenciales vector G Axy y G Ayx
respectivamente:
xy
Figura 2.20 Representación potencial vector G A . Fuente en (0.25 Ȝ, 0.25 Ȝ) y 20 dipolos imágenes.
yx
Figura 2.21 Representación potencial vector G A . Fuente en (0.25 Ȝ, 0.25 Ȝ) y 20 dipolos imágenes.
-27-
Capítulo 2: Formulación espacial en cavidades circulares
yy
Figura 2.22 Representación potencial vector G A . Carga en (0.25 Ȝ, 0.25 Ȝ) y 20 dipolos imágenes.
Figura 2.23 Situación de 20 cargas/dipolos imagen alrededor de la cavidad con la fuente en (0.625,0)
-28-
Capítulo 2: Formulación espacial en cavidades circulares
Figura 2.24 Potencial escalar en el interior del cilindro. Fuente en (0.625 Ȝ,0 Ȝ) y 20 cargas imágenes
xx
Figura 2.26 Potencial vector G A en el interior del cilindro. Carga en (0.625 Ȝ, 0 Ȝ) y 20 cargas
imagen
-29-
Capítulo 2: Formulación espacial en cavidades circulares
xy
Figura 2.27 Representación potencial vector G A . Carga en (0.625 Ȝ, 0 Ȝ) y 20 dipolos imágenes.
yx
Figura 2.28 Representación potencial vector G A . Carga en (0.625 Ȝ, 0 Ȝ) y 20 dipolos imágenes.
Finalmente, representaremos el potencial vectorial G Ayy , que, al igual que el resto
de potenciales, presenta simetría respecto al eje ‘x’:
yy
Figura 2.29 Representación potencial vector GA . Carga en (0.625 Ȝ, 0 Ȝ) y 20 dipolos imágenes.
-30-
Capítulo 2: Formulación espacial en cavidades circulares
Figura 2.30 Situación de 20 cargas/dipolos imagen alrededor de la cavidad con la fuente en (0.7, 0.7)
-31-
Capítulo 2: Formulación espacial en cavidades circulares
Figura 2.32 Potencial escalar en el interior del cilindro. Fuente en (0.7 Ȝ,0.7 Ȝ) y 20 cargas imágenes
Figura 2.33 Representación del potencial escalar en el contorno. Fuente en (0.7 ,0.7) y 20
cargas/dipolos imágenes
La presencia de este pico de potencial, que evita que se cumplan las condiciones
de contorno, es uno de los problemas que tiene esta aproximación: el error introducido
cuando la fuente se encuentra cerca de las paredes. En el capítulo siguiente
introduciremos varias técnicas con las que conseguiremos reducir de forma significativa
este error.
-32-
Capítulo 2: Formulación espacial en cavidades circulares
xx
Figura 2.34 Potencial vector G A en el interior del cilindro. Fuente en (0.7 Ȝ, 0.7 Ȝ) y 20 cargas
imagen
Comprobamos como éste aparece únicamente en el borde del cilindro, en la
situación donde se encuentra la fuente. Para el resto de posiciones del interior de la
cavidad, el valor del potencial es nulo.
xy
Figura 2.35 Representación potencial vector G A . Fuente en (0.7 Ȝ, 0.7 Ȝ) y 20 dipolos imágenes.
-33-
Capítulo 2: Formulación espacial en cavidades circulares
yx
Figura 2.36 Representación potencial vector G A . Fuente en (0.7 Ȝ, 0.7 Ȝ) y 20 dipolos imágenes.
yy
Figura 2.37 Representación potencial vector G A . Fuente en (0.7 Ȝ, 0.7 Ȝ) y 20 dipolos imágenes.
Por último comentar que para este tipo de casos en los que la fuente se encuentra
muy cercana a las paredes del cilindro, un aumento del número de cargas/dipolos
imagen no tiene porqué provocar un mejor cumplimiento de las condiciones de
contorno. Es más, un mayor número de imágenes puede provocar que, a pesar de
cumplirse los nulos en los puntos impuestos, el resto de puntos intermedios se vea
sometido a unos valores de potencial muy elevados dejándose de cumplir por tanto las
condiciones de contorno en esos puntos e introduciendo un mayor error en la
aproximación de los potenciales obtenida.
-34-
Capítulo 2: Formulación espacial en cavidades circulares
5) Conclusiones
-35-
Capítulo 2: Formulación espacial en cavidades circulares
En este caso, la formulación planteada parte de los mismos principios que para
el caso de fuentes eléctricas, pero aparece una diferencia importante: las condiciones de
contorno que se deben satisfacer en las paredes del cilindro son diferentes, pues las
fuentes en este caso son magnéticas.
La distribución de las cargas/dipolos alrededor de la cavidad circular sigue
exactamente la misma forma que en el caso anterior, definiéndose un número de puntos
tangentes alrededor del contorno del cilindro con sus planos tangentes, y asignando a
cada uno de estos puntos una carga/dipolo imagen simétrica respecto a la fuente y los
planos tangentes.
En este caso, el campo magnético creado por la fuente
& será:
& &
H jZF Im jZF j
& F (2.27)
ZP0H 0
La condición de contorno que debe cumplir el campo magnético en el contorno
del cilindro será el valor nulo de su componente normal, es decir:
&
H eˆU 0 en la pared del cilindro (2.28)
Siendo ê p el vector unitario radial, en coordenadas cilíndricas, como podemos
observar en la figura 2.38:
Para el caso del potencial escalar eléctrico, vamos a imponer que la componente
normal de su gradiente sea nula en el contorno del cilindro, como muestra la ecuación
2.29:
(I m ) eˆ U 0 en la pared del cilindro (2.29)
-36-
Capítulo 2: Formulación espacial en cavidades circulares
& &
Siendo Gw r , r ´ el potencial escalar magnético de una carga en condiciones de
espacio libre:
& &
jko r r ´
& & 1 e
G w r , r ´ & & (2.32)
4SP 0 r r ´
Al igual que ocurría para el caso de fuentes eléctricas, será necesario imponer las
condiciones de contorno para el potencial escalar eléctrico en más de un punto de la
pared del cilindro. Para ello, emplearemos la misma estrategia que hasta ahora, pero
extendiéndolo a un total de N puntos tangentes en la pared del cilindro, de tal forma que
exigiremos que se cumplan simultáneamente los requisitos en todos los puntos.
Así planteamos el sistema de ecuaciones:
N
& & & &
¦ q k G w ri , rk ´ eˆ pi G w ri , ro ´ eˆ pi
k 1 (2.33)
i 1,2,3...N
Comprobamos que en este caso tenemos una dificultad añadida, el cálculo del
& &
gradiente del potencial escalar eléctrico, Gw ri , rk ´ , No obstante, para espacio libre
podemos desarrollar este gradiente de una forma sencilla.
Si trabajamos con coordenadas cilíndricas, el gradiente será:
wt 1 wt 1 wt
t eˆr eˆT eˆM (2.34)
wr r wT r sin(T ) wM
-37-
Capítulo 2: Formulación espacial en cavidades circulares
Figura 2.39 Dipolo fuente en el eje ‘x’ y sus dos dipolos imagen asociados a un punto tangente
Al igual que para el caso del potencial escalar, impondremos en primer lugar las
condiciones de contorno a cumplir en un punto discreto ‘i’ de la pared del cilindro.
En este caso, partimos de la ecuación 2.28 que nos definía que la componente
normal del campo magnético en la pared del cilindro debe de ser nula, lo que tras el
&
análisis del campo magnético H (definido en la ecuación 2.27) nos lleva a dos
condiciones diferentes a imponer:
&
eˆ p F 0 (2.37)
&
eˆ p ( F ) 0 (2.38)
Así, impondremos de forma conjunta sobre el punto tangente ‘i’ las dos
condiciones definidas por las ecuaciones 2.37 y 2.38.
Para satisfacer estas condiciones, emplearemos al igual que en el caso de fuentes
eléctricas dos dipolos ortogonales, uno orientado en el eje ‘x’ y el otro en el eje ‘y’, tal y
como podemos observar en la figura 2.39
Cada uno de estos dipolos orientados tendrá un determinado valor, que
deberemos obtener resolviendo el sistema impuesto en el punto tangente ‘i’.
&
Que nos indica que la componente normal del vector potencial F debe
ser nula, lo que nos anularía el primer término del campo eléctrico en ese punto
de la ecuación 2.27.
Así, esta condición implica:
-38-
Capítulo 2: Formulación espacial en cavidades circulares
&
Donde hemos tenido en cuenta la presencia de la carga fuente, en r0 .
-39-
Capítulo 2: Formulación espacial en cavidades circulares
U
(2.47)
& & G yy
r&i , r&k ´
C y
i ,k eˆ p G ri , rk ´
yy
F
F
i 1,2,3...N
Estando definidas las constantes Cix,k y Ciy,k por la ecuación 2.47. Ambas
constantes pueden ser evaluadas para un medio general multicapa mediante el empleo
de la transformada de Sommerfeld [9], aunque para el espacio libre, ambas constantes
coinciden a un valor dado.
Para obtener este valor, evaluaremos las partes en las que se compone la
ecuación 2.47a):
& & & &
G Fxx ri , rk ´ H 0 e jk r r ´ 1 0
Y b):
& &
& & H0 e
jko r r ´ § ·
eˆU G ri , rk ´
xx
& & ¨ >eˆ r& r&´@ ( & jk0& 1 ) ¸ (2.50)
F
¨ U & &
4S r r ´ © r r ´ r r ´ 2 ¸¹
-40-
Capítulo 2: Formulación espacial en cavidades circulares
De esta forma, ya tenemos todos los elementos para poder resolver la ecuación
2.48, de la que obtendremos el valor de los dipolos para poder recuperar el vector
potencial eléctrico en el interior de la cavidad cilíndrica, mediante:
N
xx & & & & &
G Fcy (r ) G Fxx r , r0 ´ ¦ M kx G Axx ri , rk ´
k 1
N
(2.52)
& & &
G yx
Fcy (r ) ¦ M G ri , rk ´
k 1
k
y yx
A
2.2.2 Resultados
-41-
Capítulo 2: Formulación espacial en cavidades circulares
xx
Figura 2.41 Potencial vector eléctrico G F normalizado en el eje x=0
xx
Figura 2.42 Potencial vector eléctrico G F normalizado en el eje y=0
yx
Figura 2.43 Potencial vector eléctrico G F normalizado evaluado con 20 imágenes
-42-
Capítulo 2: Formulación espacial en cavidades circulares
En primer lugar, deberemos introducir los datos iniciales al programa, que serán:
-43-
Capítulo 2: Formulación espacial en cavidades circulares
o ludcmp.f90
Esta rutina se emplea para realizar la factorización LU de una matriz de números
complejos, tal y como se detalla en [12].
o lubksb.f90
Esta rutina se emplea para realizar la sustitución hacia atrás de una matriz de
coeficientes complejos, que ha sido factorizada en su forma LU previamente,
mediante la rutina ludcmp.f90, como se explica en [12]. Lo emplearemos para
poder resolver sistemas de ecuaciones lineales complejas.
o nrtype.f90
Es un módulo que se encarga de definir algunas constantes y parámetros usados
de forma habitual en el cálculo numérico.
o nrutil.f90
Es un módulo que define una gran cantidad de subrutinas numéricas empleadas
de forma frecuente, y que han sido definidas en [12].
o const1r.90
Este módulo se encarga de la definición de algunas constantes físicas que
emplearemos a lo largo del programa principal, como por ejemplo el número pi,
la permitividad del vacío, etc.
-44-
Capítulo 2: Formulación espacial en cavidades circulares
o pgenegeometry.f90
Este módulo define el nombre de variables, para que sean usadas en varios
programas con el mismo nombre, facilitando así la portabilidad entre ellos. Por
ejemplo, para la frecuencia (‘freq’), el número de onda (‘k0’), etc.
o pescalarf_contorno.F90
Esta subrutina se encarga de la situación de N puntos tangentes alrededor del
cilindro, ubicar las cargas imágenes alrededor de los mismos, y calcular el valor
complejo de las cargas para que se satisfagan las condiciones de contorno en
esos N puntos. Para ello, implementa la teoría del potencial escalar bajo fuentes
eléctricas que ha sido desarrollada en el punto 2.1
Recibirá toda la información de entrada necesaria (como el número de puntos
tangentes, el radio del cilindro, ubicación de la carga fuente, etc) y devolverá el
potencial escalar calculado, tanto en todo el interior del cilindro como en el
contorno del mismo (pudiéndose así verificar si se cumplen o no las condiciones
de contorno).
Por último, también grabará para su posterior representación la situación de las
cargas imágenes, puntos tangentes en el cilindro y la fuente, en el subdirectorio
“resultados”.
o pvectorf.F90
Esta subrutina se encarga de la situación de N puntos tangentes alrededor del
cilindro, ubicar los dipolos imágenes alrededor de los mismos, y calcular el
valor complejo de los dipolos para que se satisfagan las condiciones de contorno
en esos N puntos. Para ello, implementa la teoría del potencial vectorial
magnético bajo fuentes eléctricas que ha sido desarrollada en el punto 2.1 Este
cálculo se realizará bajo dos suposiciones: que el dipolo fuente se encuentra en
su ubicación orientado respecto al eje ‘x’ y respecto al eje ‘y’, pudiendo así
obtener todos los potenciales vectoriales.
o Potril.F90
Se trata del programa principal, que se encargará de extraer los parámetros de
análisis del fichero “Datos.in” e ir llamando a las distintas funciones para
calcular todos los potenciales (escalar y vectorial). Finalmente, grabará en el
subdirectorio “resultados” los valores obtenidos para todos los cálculos.
-45-
Capítulo 2: Formulación espacial en cavidades circulares
Hemos empleado dos ficheros script de Matlab, para poder mostrar gráficamente
los resultados:
o representa_situacion.m
Es una función que se encarga de coger toda la información de la ubicación de
las imágenes, los puntos tangentes del cilindro, la posición de la fuente, y la
forma del cilindro, para poder visualizar la situación que estamos analizando.
o inicio.m
Es el programa principal, que además de llamar a la función para representar la
situación en la que nos encontramos, se encarga de extraer de los ficheros toda la
información de los potenciales y visualizarlos de forma gráfica, en 3D y
mostrando el corte con los ejes de todos los potenciales. Para el caso del
potencial escalar, se incluye además el potencial obtenido a lo largo del contorno
del cilindro.
2.3.3.2 Ejecución
Los dos scripts de Matlab se encuentran por defecto dentro del subdirectorio
“resultados” que es donde se guardan todos ficheros de información para el cálculo de
los potenciales. En ese directorio, habrá que ejecutar el programa “inicio.m”,
escribiendo “inicio” en la línea de comandos de Matlab ó abriendo el fichero y
ejecutándolo.
-46-
Capítulo 3: Optimización de la función espacial de Green
o Potencial escalar:
-47-
Capítulo 3: Optimización de la función espacial de Green
o Potencial vectorial:
Tras haber optimizado el Coste, tanto escalar como vectorial, por alguna de las
técnicas descritas anteriormente, vamos a emplearlo para calcular el error obtenido en la
aproximación de los potenciales en función de la situación de las carga/dipolo fuente,
comprobando como el Coste se reduce de forma más que considerable.
Así, podemos limitar el número de cargas/dipolos imagen para conseguir una
determinada exactitud en la aproximación; ello es muy aconsejable debido a los
problemas de sobredimensionamiento que se producen si se imponen demasiadas
condiciones de contorno, apareciendo problemas a la hora de resolver los sistemas de
ecuaciones pues no están bien condicionados.
Representaremos gráficamente como se distribuye cada uno de los Costes
alrededor del cilindro, para observar cuales son las zonas en la que se produce un
aumento del error y tratar de minimizarlas.
En cuanto al coste computacional, éste varía en función de la técnica y del
número de imágenes empleadas, siendo superior al coste que obteníamos en la
implementación inicial (presentada en el capítulo 2). No obstante, conseguimos
contener el coste total, depurando cada una de las técnicas con algoritmos de
parada/detención que se detendrán de forma automática cuando hallamos alcanzado la
precisión deseada (en caso de que sea posible).
En este caso particular, una geometría cilíndrica, en el que los errores son muy
reducidos con un número no muy elevado de imágenes (15-20) la optimización puede
ser suplida con un aumento de imágenes, que consiguen minimizar el error de la misma
forma. No obstante, todas las técnicas aquí descritas no toman como partida la
geometría cilíndrica, sino que pueden ser empleadas para cualquier otro tipo de
estructura arbitraria (como veremos en el Capítulo 4) donde el error obtenido por las
imágenes sin optimizar es demasiado elevado, siendo incorrecta la aproximación de los
potenciales obtenida.
-48-
Capítulo 3: Optimización de la función espacial de Green
-49-
Capítulo 3: Optimización de la función espacial de Green
Figura 3.1 Representación 15 cargas imágenes alrededor del cilindro, fuente en origen coordenadas
Figura 3.2 Representación del potencial escalar en contorno del cilindro en función del ángulo.
Empleo de 15 imágenes, fuente en origen de coordenadas
-50-
Capítulo 3: Optimización de la función espacial de Green
Por ejemplo, basándonos en las suposiciones del caso anterior, pero situando la
carga en la posición (0.7Ȝ, 0.7 Ȝ) obtenemos una situación y una distribución del coste
como muestra la figura 3.3:
Figura 3.3 Representación de situación y coste del potencial escalar en un cilindro, empleo de 15
imágenes y fuente situada en (0.7 Ȝ,0.7 Ȝ)
Así:
-51-
Capítulo 3: Optimización de la función espacial de Green
Nº cargas 5 10 15 20 25 30
Nº cargas 5 10 15 20 25 30
x c) Posición (0.4Ȝ,0Ȝ):
Nº cargas 5 10 15 20 25 30
2.5848e- 4.6492e-
Coste 16.1649 19.9669 0.3775 0.0026
005 007
Tabla 3.3 Coste método original, posición fuente en (0.4Ȝ ,0Ȝ)
x d) Posición (0.4Ȝ,0.4Ȝ):
Nº cargas 5 10 15 20 25 30
-52-
Capítulo 3: Optimización de la función espacial de Green
x e) Posición (0.5Ȝ,0Ȝ):
Nº cargas 5 10 15 20 25 30
9.6659e- 3.7767e-
Coste 17.9904 19.9656 0.9056 0.0070
005 006
Tabla 3.5 Coste método original, posición fuente en (0.5Ȝ ,0Ȝ)
x f) Posición (0.7Ȝ,0.7Ȝ):
Nº cargas 5 10 15 20 25 30
Tras analizar los resultados para el caso inicial, podemos extraer dos
conclusiones interesantes:
-53-
Capítulo 3: Optimización de la función espacial de Green
i 1,2,3...N
Finalmente, debido a la presencia del pulso, la integral únicamente estará
definida dentro de las posiciones del mismo, es decir, la integral de la función de green
en el contorno del segmento que está delimitado por la función pulso.
N
& & & &
¦
k 1
qk ³ Gv (ri , rk ' ) dri ³ Gv (ri , ro ' ) dri
(3.6)
segmento _'i ' segmento _'i '
i 1,2,3...N
Siendo este el sistema que evaluaremos y que nos proporcionará el valor
complejo de las cargas imagen.
-54-
Capítulo 3: Optimización de la función espacial de Green
Por su parte la función de Green que va desde los puntos tangentes (r1) hasta la
fuente (r0’), seran:
De esta forma, podemos partir del sistema de ecuaciones inicial para el caso
Point-Matching que como vimos en el capítulo 2 resultaba:
ª Gv (r1 , r1 ' ) Gv (r1 , r2 ' ) .... Gv(r1 , rN ' ) º ª q1 º ª Gv(r1 , ro ' ) º
« Gv(r , r ' ) » « » « »
« 2 1 .... .... .... » * « q 2 » « Gv(r2 , ro ' ) » (3.10)
« .... .... .... .... » « .... » « ... »
« » « » « »
¬Gv(rN , r1 ' ) .... .... Gv(rN , rN ' )¼ ¬q N ¼ ¬Gv(rN , ro ' )¼
y renombrarlo de forma:
-55-
Capítulo 3: Optimización de la función espacial de Green
ª
³ Gv1 (r1 , r2 ' ) * dr1
« segmento ³ Gv(r , r ' ) * dr
segmento1
1 2 1 .... ³ Gv(r , r
segmento1
1 N' ) * dr1 º
» ªq º
ª
³ Gv1 (r1 , ro ' ) * dr1 º»
« segmento
« » « 1» « »
« ³ Gv(r2 , r1 ' ) * dr2 .... .... .... » « q2 » « ³ Gv(r2 , ro ' ) * dr2 »
« segmento 2 » * « .... » « segmento 2 » (3.12)
« .... .... .... .... » « » « ... »
« » ¬q N ¼ « »
« ³ ³ « ³
Gv(rN , r1 ' ) * drN .... .... Gv ( r , r ' ) * dr Gv ( r , r ' ) * dr
N N N
» N o N
»
¬ segmentoN segmentoN ¼ ¬ segmentoN ¼
-56-
Capítulo 3: Optimización de la función espacial de Green
ª F 2 F 2' .... FN ( n ) º ª q1 º ªG 2º
« » « » « G 3»
« F 3 .... .... .... » « q 2 »
* « » (3.16)
« .... .... .... .... » « .... » « ... »
« » « » « »
¬« F1 .... .... F1( N ) ¼» ¬q N ¼ ¬ G1¼
Una condición necesaria y suficiente para que el sistema tenga solución es que el
determinante de la matriz ‘F’ exista y sea distinto de 0. Es decir, que:
ª F1 F 2 F1' F 2' .... F1( n ) F 2 ( n ) º
« »
«F 2 F 3 .... .... .... » z0 (3.17)
« .... .... .... .... »
« »
¬« FN F1 .... .... FN F1( N ) ¼»
(N)
siendo ‘N’ el rango de la matriz, que en nuestro caso, coincide con el número de puntos
tangentes impuestos sobre el contorno del cilindro.
Por tanto, podemos concluir sobre la aproximación mediante Trapecios que:
-57-
Capítulo 3: Optimización de la función espacial de Green
a) Simpson’s rule
-58-
Capítulo 3: Optimización de la función espacial de Green
siendo h la distancia entre F1 y F2 en el cilindro (es decir, entre dos puntos contiguos de
un mismo segmento).
c) Bode rule
-59-
Capítulo 3: Optimización de la función espacial de Green
siendo h la distancia entre F1 y F2 en el cilindro (es decir, entre dos puntos contiguos de
un mismo segmento).
Este caso es una generalización del caso anterior, empleando un mayor número
de puntos internos, y lo utilizamos para ver la influencia de los mismos en el
cumplimiento de las condiciones de contorno del cilindro.
-60-
Capítulo 3: Optimización de la función espacial de Green
3.2.2.4 Aproximación del Trapecio tomando una gran cantidad de puntos intermedios
Figura 3.9 Detalle de los puntos intermedios tomados en un arco del cilindro al realizar una
aproximación por Trapecio general
Observamos como en este caso tomamos suficientes puntos para conseguir unos
resultados más satisfactorios, quedando cada uno de los arcos completamente
caracterizados.
Para obtener unos resultados comparativos entre las distintas técnicas empleadas
y la técnica original de Point-Matching, fuimos variando la situación de la carga fuente
en el interior del cilindro y el número de cargas imágenes (y por tanto puntos tangentes
sobre los que imponer las condiciones de contorno) con el que abordaremos el
problema, observando el error cometido en las condiciones de contorno en cada caso.
Tras implementar los métodos anteriores, se han obtenido unos resultados no
muy satisfactorios, siendo la técnica inicial de Point-Matching (es decir establecimiento
de las condiciones de contorno en puntos discretos finitos) muy superior en la mayoría
de las ocasiones siendo además su coste computacional más reducido.
-61-
Capítulo 3: Optimización de la función espacial de Green
x Las aproximaciones integrales cerradas (es decir, las de Simpson, Simpson 3/8
y Bode) presentan el mismo problema que la aproximación por Trapecio, siendo
únicamente correctas cuando empleamos un número de cargas imágenes impar,
pues en caso contrario, obtenemos una sistema de ecuaciones incompatible (que
es resuelto de forma errónea por el software, dando por tanto resultados con un
error mucho mayor que en el resto de los casos). Para el caso en el que son
correctas, consiguen un resultado muy similar al obtenido con Point-Matching,
mejorando levemente el error cometido cuando la carga se encuentra situada en
el centro del cilindro, y empeorando los resultados cuando la carga se va
acercando a las paredes. Su empleo no es práctico, pues depende del número de
cargas imagen escogidas (par o impar), aumenta su coste computacional y
proporciona unos resultados ‘oscilantes’ (a veces mejora a veces empeora) frente
a los de Point-Matching.
-62-
Capítulo 3: Optimización de la función espacial de Green
Este método es uno de los más intuitivos que presentaremos, siendo capaz de
proporcionar unos resultados moderados, ligeramente superiores a los de Point-
Matching para situaciones particulares.
A diferencia de los casos anteriores, variaremos de forma dinámica el número de
cargas imágenes empleadas. De esta forma, comenzaremos analizando los potenciales y
la distribución del mismo alrededor del contorno de la estructura (para ver el grado de
cumplimiento de las condiciones impuesta), buscaremos la posición en la que se
produce el máximo de este potencial y lo contrarrestaremos introduciendo una nueva
condición de contorno en ese punto.
Para la aplicación de este método es necesario una situación inicial de partida
(una configuración de cargas imágenes, posicionamiento de las mismas y puntos
tangentes sobre los que imponer las condiciones de contorno), sobre la cual se irá
evolucionando introduciendo más cargas imágenes que controlaran el potencial en
puntos discretos correspondientes a máximos de potencial.
Veamos el funcionamiento del método mediante un ejemplo; supongamos una
configuración con 14 cargas imágenes, situadas alrededor de un cilindro infinitesimal en
la que se ha colocado una fuente puntual en la posición (0.7Ȝ,0.7Ȝ), tal y como indica la
figura siguiente:
Figura 3.10 Representación de situación y coste del potencial escalar en un cilindro, empleo de 14
imágenes y fuente situada en (0.7 Ȝ,0.7 Ȝ)
-63-
Capítulo 3: Optimización de la función espacial de Green
Figura 3.11 Representación de situación y coste del potencial escalar en un cilindro tras ubicar un
punto tangente en un pico de potencial, empleo de 15 imágenes y fuente situada en (0.7 Ȝ,0.7 Ȝ)
Realizando un zoom podemos observar con más detalle la posición del nuevo
punto tangente donde forzar la condición de contorno y de la carga imagen asociada:
Figura 3.12 Detalle de la ubicación de un nuevo punto de condición de potencial nulo donde se
encontraba anteriormente un pico de potencial.
-64-
Capítulo 3: Optimización de la función espacial de Green
De esta forma, se impone alrededor del pico de potencial una gran cantidad de
condiciones de nulidad, forzando a éste a tener valores más reducidos. Gráficamente se
puede observar en la figura 3.13 la distribución del potencial en el caso que presenta un
menor coste, es decir, 25 cargas imágenes habiendo evolucionado las mismas desde una
situación inicial de 19.
-65-
Capítulo 3: Optimización de la función espacial de Green
Figura 3.13 Detalle de la distribución del potencial, 25 cargas imágenes habiendo evolucionado su
colocación en picos de potencial desde 19 cargas. Fuente en (0.7 Ȝ, 0.7Ȝ)
Figura 3.14 Distribución del potencial a lo largo del contorno del cilindro con 15 imágenes.
a) Situando las cargas en picos de potenciales y evolucionando desde una situación inicial de 10
imágenes. b) Empleando Point-Matching.
-66-
Capítulo 3: Optimización de la función espacial de Green
Figura 3.15 Situación 20 cargas imagen alrededor cilindro, fuente en (0.5Ȝ ,0.5Ȝ)
-67-
Capítulo 3: Optimización de la función espacial de Green
-68-
Capítulo 3: Optimización de la función espacial de Green
Comprobamos como en este caso la ubicación de las cargas imagen coincide con
la ubicación de las imágenes con el método inicial. No obstante es en la variación del
radio donde radica la potencia de este método de optimización.
A priori no podemos conocer que valor del radio será el que proporcione un
coste mínimo. Por esta razón, representamos el coste en función del radio, variando el
mismo entre 1 y 7600 aproximadamente (es decir un rango muy amplio, con un paso de
valores entre radio ‘grueso’, para obtener así una primera aproximación)
De esta forma, obtenemos una variación del error en función del radio de:
Figura 3.18 Coste en función del radio (0-8000). 20 imágenes, fuente en (0Ȝ,0Ȝ).
-69-
Capítulo 3: Optimización de la función espacial de Green
Figura 3.19 Coste en función del radio (0-11). 20 imágenes, fuente en (0Ȝ,0Ȝ).
Observamos en la figura 3.8 como el error aumenta cuanto más cerca está la
carga del cilindro, descendiendo hasta valores menores de error conforme situamos el
radio a una mayor distancia. Comprobamos como para la ubicación original (en la que
el radio era igual a 2O ) el coste ya es muy reducido, no obstante con un aumento del
radio obtendremos costes todavía menores, como podemos comprobar en la figura 3.9:
Figura 3.20 Coste en función del radio (2-7). 20 imágenes, fuente en (0Ȝ,0Ȝ).
Así, hemos obtenido una reducción del coste más que notoria, pasando de un
error de 7.1487 10 -5 a uno de 2.7080 10 -7 , sin más que emplear un radio de las cargas
imágenes de valor 7O .
También indicar que el error converge, aunque en este caso al estar situada la
carga en el origen la convergencia se daría teóricamente en el infinito.
-70-
Capítulo 3: Optimización de la función espacial de Green
Figura 3.21 Coste en función del radio, paso grueso y fino. 20 imágenes y fuente en (0Ȝ,0Ȝ)
-71-
Capítulo 3: Optimización de la función espacial de Green
Supongamos por ejemplo la misma situación que hemos estado analizando hasta
ahora, pero empleando un total de 10 cargas imágenes. La variación del coste en
función del radio será:
Figura 3.22 Variación del coste con el radio variable (0-80). 10 cargas imágenes, fuente en (0Ȝ,0Ȝ)
-72-
Capítulo 3: Optimización de la función espacial de Green
Figura 3.24 Potencial escalar en el contorno, métodos original y variación de radio. Fuente situada
en el origen y 20 imágenes
x a) Posición (0Ȝ,0Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Técnica
7.1487e- 1.3369e- 3.07174e-
Point-matching 10.1294 1.3525 0.0061
005 006 008
2.6108e- 1.3858e- 2.7904e-
Error 9.6952 1.1228 0.0013
007 011 014
Variación
1.400Ȝ 3.100 Ȝ
del radio 11Ȝ 11 Ȝ 11 Ȝ 4Ȝ
Radio mínimo mínimo
converge converge converge converge
local local
Tabla 3.8 Coste variación radio para distintos nº de cargas, posición fuente en (0Ȝ ,0 Ȝ)
-73-
Capítulo 3: Optimización de la función espacial de Green
Figura 3.25 Variación del coste con el radio, 5 imágenes y fuente en el origen.
x b) Posición (0.2Ȝ,0.2Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Técnica
5.6664e- 5.9574e- 4.0841e-
Point-matching 12.1475 16.5825 0.0972
004 006 008
3.8533e- 4.9440e- 1.9817e-
Error 10.9688 3.2548 0.0352
005 007 008
Variación
1.200 Ȝ 5.200 Ȝ
del radio 2.400 Ȝ 4 .400 Ȝ 2.500 Ȝ 16.000 Ȝ
Radio mínimo mínimo
converge converge converge converge
local local
Tabla 3.9 Coste variación radio para distintos nº de cargas, posición fuente en (0.2Ȝ ,0.2Ȝ)
-74-
Capítulo 3: Optimización de la función espacial de Green
Figura 3.26 Coste en función del radio. 10 cargas imágenes y fuente en (0.2Ȝ ,0.2Ȝ)
x c) Posición (0.4Ȝ,0.0Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Técnica
2.5848e- 4.6492e-
Point-matching 16.1649 19.9669 0.3775 0.0026
005 007
5.8839e- 1.8230e- 3.2894e-
Error 10.3919 6.8724 0.0862
004 005 007
Variación
1.200 Ȝ 4.400 Ȝ
del radio 2.400 Ȝ 2.300 Ȝ 2.000 Ȝ 2.000 Ȝ
Radio mínimo mínimo
converge converge converge converge
local local
Tabla 3.10 Coste variación radio para distintos nº de cargas, posición fuente en (0.4Ȝ ,0Ȝ)
x d) Posición (0.4Ȝ,0.4Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Técnica
2.8172e- 2.5517e-
Point-matching 26.6829 23.6945 1.4449 0.0132
004 005
6.2641e-
Error 11.1007 9.1644 0.2953 0.0256 0.0019
004
Variación
1.200 Ȝ 1.200 Ȝ
del radio 2.900 Ȝ 16 Ȝ 1.700 Ȝ 1.800 Ȝ
Radio mínimo mínimo
converge converge converge converge
local local
Tabla 3.11 Coste variación radio para distintos nº de cargas, posición fuente en (0.4Ȝ, 0.4Ȝ)
Comprobamos como en esta posición hay tres casos en los que una variación del
radio no proporciona un coste menor que el original. Esto se produce para un número de
imágenes relativamente elevado (20-30) cuando el error ya es reducido. Destacamos la
utilidad del método para reducir de forma muy importante el coste para un número de
imágenes más reducido, en algunos casos de forma espectacular (5-10 imágenes).
-75-
Capítulo 3: Optimización de la función espacial de Green
x e) Posición (0.5Ȝ,0.0Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Técnica
9.6659e- 3.7540e-
Point-matching 17.9904 19.9656 0.9056 0.0070
005 006
2.4757e- 9.1774e-
Error 10.1965 8.5877 0.1721 0.0047
004 006
Variación
1.200 Ȝ 1.200 Ȝ
del radio 2.900 Ȝ 1.900 Ȝ 1.800 Ȝ 1.800 Ȝ
Radio mínimo mínimo
converge converge converge converge
local local
Tabla 3.12 Coste variación radio para distintos nº de cargas, posición fuente en (0.5Ȝ, 0.0Ȝ)
x f) Posición (0.7Ȝ,0.7Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Técnica
Point-matching 22.8267 8.0995 1.6061 4.2177 0.4510 1.1371
Error 24.2562 20.4576 11.5060 11.9342 8.1999 8.5931
Variación 1.1000 1.0900 1.0600 1.0400 1.1300
del radio Radio 1.4300
mínimo mínimo mínimo mínimo minimo
converge
local local local local local
Tabla 3.13 Coste variación radio para distintos nº de cargas, posición fuente en (0.7Ȝ, 0.7Ȝ)
Podemos observar como en casi todos los casos obtenemos un coste mucho
menor que con la técnica original, exceptuando claramente la posición (0.7Ȝ, 0.7Ȝ). En
este caso, la carga fuente se encuentra muy próxima a la pared del cilindro, por lo que
produce un pico de potencial en ella que las cargas imágenes de nuestro método, al estar
situadas de forma simétrica alrededor del cilindro, no pueden contrarrestar. No obstante,
si la carga se encuentra cerca de la pared del cilindro (§0.9 Ȝ) otras técnicas se muestran
mejores que con la anterior ubicación, especialmente al emplear un número reducido de
cargas imágenes.
El error se disminuye de una forma más drástica cuando el número de cargas
imágenes es reducido (ó el coste inicial de partida muy elevado) y de una forma más
moderada cuando el número de cargas imágenes es más elevado.
En cuanto al efecto del radio en el coste, hemos comprobado que siempre se va a
producir un coste mínimo, que puede o no coincidir con el coste al que converge cuando
el radio aumenta mucho su valor (es decir, si el coste mínimo es el valor de
convergencia o es un mínimo absoluto que se presenta anteriormente).
Como principal desventaja de este método, señalaremos que el coste
computacional para conseguir los resultados es relativamente alto, pues debemos
calcular los potenciales para un gran número de situaciones y obtener el Coste para cada
uno de ellos (operación computacionalmente costosa).
-76-
Capítulo 3: Optimización de la función espacial de Green
Este método se basa en modificar el valor complejo de las cargas imágenes con
el fin de optimizarlo, para poder conseguir así un coste mínimo.
Tras resolver el sistema de ecuaciones definido en la ecuación 3.1 ó 3.4 (que son
la misma), obtenemos un conjunto de N valores complejos, uno para cada carga ‘qx’
imagen. Con estos valores complejos de las cargas calcularemos el potencial escalar en
el interior del cilindro (o la estructura arbitraria que hayamos definido), como
estudiamos de forma detallada en el Capitulo 2.
Por tanto deberemos de optimizar el valor de las cargas (que es una magnitud
compleja) empleando para ello un método de descenso de gradiente. Es necesario
definirnos de forma adecuada un Coste a optimizar (que no tiene porqué coincidir con el
coste de la aproximación del potencial escalar definido en la ecuación 3.3) que nos dará
el valor de las cargas buscado en cada iteración.
En general, una optimización por descenso de gradiente sigue la estructura:
w (k 1) w (k ) P w C (k ) (3.27)
Para nuestro caso, vamos a optimizar el valor complejo de las cargas, por lo cada
término de la ecuación se refiere a:
x P Æ Constante de paso.
Por tanto la ecuación que estamos definiendo en 3.27 puede escribirse como:
-77-
Capítulo 3: Optimización de la función espacial de Green
§ w ·
¨ C (k ) ¸
§ q1 (k 1) · § q1 (k ) · ¨ wq1 ¸
¨ ¸ ¨ ¸ ¨ w ¸
C (k ) ¸
¨ q 2 (k 1) ¸ ¨ 2q ( k ) ¸ ¨ wq 2
¨ ¸ ¨ ... ¸ P ¨ ¸
... ¨ ... ¸ (3.28)
¨ ¸ ¨ ¸
¨ q N 1 (k 1) ¸ ¨ q N 1 (k ) ¸ ¨ w C (k ) ¸
¨ q (k 1) ¸ ¨ q (k ) ¸ ¨ wq N 1 ¸
© N ¹ © N ¹ ¨ w ¸
¨¨ C (k ) ¸¸
© wq N ¹
Será por tanto necesario calcular una función de Coste que sea derivable en
función de las cargas complejas.
Para ello realizaremos dos versiones:
y como:
K
d ª º
wp Gvcy
dq p «
¬
Gv ( r , ro ) ¦
k 1
qk Gv(r , rk )» Gv(r , rp )
¼
(3.31)
-78-
Capítulo 3: Optimización de la función espacial de Green
Es por todas estas razones por la que este resultado no es adecuado para el
cálculo del gradiente del coste.
En este punto definiremos el coste a optimizar como la integral del cuadrado del
valor absoluto de la función de Green en el contorno de la estructura, tal y como
observamos en la ecuación 3.33:
2
C oste ³ Gv
Contorno
cy dc (3.33)
*
Donde la notación Gvcy implica el número complejo conjugado de Gvcy .
De esta forma podemos desarrollar la ecuación 3.33 teniendo en cuenta la
definición del potencial escalar Gvcy :
2 *
C oste ³ Gv
Contorno
cy dc ³ Gv
Contorno
cy Gvcy dc
K K * (3.35)
ª º ª º
³ «Gv ( r , ro ) ¦ q k Gv ( r , r )
k » « Gv ( r , ro ) ¦ q k Gv(r , rk )» dc
Contorno ¬ k 1 ¼ ¬ k 1 ¼
-79-
Capítulo 3: Optimización de la función espacial de Green
obtener la derivada del coste definido en la ecuación 3.37 respecto a cada una de las
cargas imagen que hayamos definido en nuestro problema (en general, será un número
N).
Al ser el valor de las cargas complejos, el gradiente que deberemos obtener será
de la forma:
§ wC (q) wC (q ) ·
¨ j ¸
¨ w Re( q1 ) w Im(q1 ) ¸
wC (q ) wC ( q ) ¸
1 ¨¨ j
q C ( q) w Re(q 2 ) w Im(q 2 ) ¸ (3.38)
2¨ ¸
¨ ...... ¸
¨ wC ( q ) wC (q) ¸
¨ w Re(q ) j w Im(q ) ¸
© N N ¹
Calcularemos en este punto la derivada parcial del coste respecto a una carga
compleja arbitraría ‘p’:
K
§ ·
¨ Gv(r , ro ) Gv(r , ro )* Gv(r , ro ) ¦ qk*Gv(r , rk )* ¸
w ¨ k 1 ¸
q p Coste ³
wq p Contorno ¨ K
ª K
º ª K
º ¸dc (3.39)
¨ Gv(r , ro )* ¦ qk Gv(r , rk ) «¦ qk Gv(r , rk )» «¦ qk*Gv(r , rk )* » ¸
¨ ¸
© k 1 ¬k 1 ¼ ¬k 1 ¼¹
w * w
³
Contorno
Gv(r , ro )
wq p
q p Gv( r , r p ) * ³ Gv(r , r ) wq q
Contorno
o
p
pr jq i Gv( r , r p ) *
(3.41)
1 ª w º
³ Gv( r , ro ) « q pr jqi Gv(r , rp ) * j w q pr jqi Gv(r , rp ) * » dc
Contorno
2 «¬ wq pr wq pi »¼
³ Gv(r , r ) Gv(r , r
o p )*
Contorno
-80-
Capítulo 3: Optimización de la función espacial de Green
K K
w w
³ Gv * (r , ro ) ¦ q k Gv(r , rk )dc *
³ Gv (r , ro ) ¦ q Gv(r , r )
k k
wq p Contorno k 1 Contorno
wq p k 1
w w
³
Contorno
Gv * (r , ro )
wq p
q p Gv(r , r p ) ³ Gv(r , r ) wq q
Contorno
o
p
pr jq i Gv(r , rp )
(3.42)
1 ª w º
³
*
Gv (r , ro ) « q pr jqi Gv(r , rp ) j w q pr jqi Gv(r , rp )» dc
Contorno
2 ¬« wq pr wq pi ¼»
1
³ Gv(r , r ) 2 >Gv(r , r
Contorno
o p ) Gv(r , rp ) @ 0
wq p ³ «¦ kr ( q q ki )Gv ( r , rk » «¦ ( q kr q ki )Gv ( r , rk ) » dc
)
Contorno ¬ k 1 ¼ ¬k 1 ¼
§ª K * º w
³ ¨¨ «¦ q k Gv(r , rk ) * » >
(q pr j q pi )Gv(r , rp ) @
Contorno © ¬ k 1 ¼ wq p
ªK º w
>
«¦ q k Gv(r , rk )» wq (q pr j q pi )Gv(r , rp ) dc
*
@
¬k 1 ¼ p
ª w º (3.43)
§ª K * *º 1 w
³ ¨
¨ « ¦ q k Gv ( r , rk ) » « q pr j q i Gv ( r , r p ) j q pr jq i Gv ( r , r p ) »
Contorno © ¬ k 1 ¼ 2 «
¬ wq pr wq pi »¼
ªK º 1 ª w º·
«¦ q k Gv(r , rk )» « q pr j qi Gv * (r , rp ) j w q pr jqi Gv * (r , rp )» ¸¸dc
¬k 1 ¼ 2 ¬« wq pr wq pi ¼» ¹
§ª K * *º 1 ªK º 1 ·
³ ¨
¨« ¦ q k Gv ( r , rk )
» 2 2 Gv ( r , r p ) >
«¦ q k Gv(r , rk )» 2 Gv (r , r p ) Gv (r , rp ) ¸¸dc
* *
@
Contorno © ¬ k 1 ¼ ¬k 1 ¼ ¹
ªK * *º
³ «¦ q k Gv(r , rk ) » Gv(r , rp )dc
Contorno ¬ k 1 ¼
Una vez analizados los cuatro términos en los que habíamos descompuesto la
ecuación 3.39 vamos a recomponerlos con el fin de obtener una expresión para el
gradiente del Coste:
*
§ K ·
¨ Gv(r , r ) Gv(r , r )* ª¦ q Gv(r , r )º Gv(r , r ) ¸dc
q p Coste ³ ¨
Contorno ©
o p «
¬k 1
k k »
¼
p
¸
(3.44)
¹
-81-
Capítulo 3: Optimización de la función espacial de Green
-82-
Capítulo 3: Optimización de la función espacial de Green
Figura 3.28 Representación Potencial escalar en función del ángulo. 10 imágenes y fuente en origen
-83-
Capítulo 3: Optimización de la función espacial de Green
Figura 3.29 Evolución del Coste en función del número de iteración, paso dinámico. Empleo de 10
imágenes y fuente situada en la posición (0.2Ȝ,0.2Ȝ)
Figura 3.30 Evolución del Coste en función del número de iteración, paso fijo. Empleo de 10
imágenes y fuente situada en la posición (0.2Ȝ ,0.2Ȝ)
-84-
Capítulo 3: Optimización de la función espacial de Green
Figura 3.31 Representación Potencial escalar en función del ángulo tras optimizar el valor complejo
de las cargas. Empleo de 10 imágenes y fuente situada en el origen
Podemos comprobar como los valores sobre los que nos movemos son mucho
más reducidos, manteniéndose básicamente los 10 nulos provocados por las 10
condiciones de las cargas (aunque esto no tiene porqué ser así). En general, el nuevo
potencial escalar alrededor del cilindro tras optimizar el valor complejo de las cargas no
tiene porqué guardar ninguna relación con el original, salvo descensos acusados (que no
nulos) en las puntos tangentes del cilindro donde impusimos las condiciones de
contorno.
Mostraremos a continuación cual es el valor final del coste obtenido en este
método en función del número de cargas imagen empleadas (5, 10, 15, 20, 25, 30) y de
la posición de la fuente en el interior del cilindro [ (0Ȝ,0Ȝ) , (0.2Ȝ,0.2Ȝ) , (0.4Ȝ,0Ȝ) ,
(0.4Ȝ,0.4Ȝ) , (0.5Ȝ,0Ȝ) , (0.7Ȝ,0.7Ȝ) ]. Indicaremos en rojo aquellos resultados que
mejoren respecto al valor inicial de cargas imágenes.
x a) Posición (0Ȝ,0Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Técnica
Point- 7.1487e- 1.3369e- 3.07174e-
10.1294 1.3525 0.0061
matching 005 006 008
Optimización
3.85533e- 4.549699E- 8.63487e- 1.9519e-
del valor de 10.12939 0.898036
003 005 007 008
las cargas
Tabla 3.14 Coste optimizando cargas para distintos nº de imágenes, posición fuente en (0.0Ȝ , 0.0Ȝ)
-85-
Capítulo 3: Optimización de la función espacial de Green
x b) Posición (0.2Ȝ,0.2Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Técnica
5.6664e- 5.9574e- 1.204944e-
Point-matching 12.1475 16.5825 0.0972
004 006 007
Optimización del
2.44732 4.40864e- 4.54379e- 8.99507e-
valor de las 10.7727 0.08024
004 006 008
cargas
Tabla 3.15 Coste optimizando cargas para distintos nº de imágenes, posición fuente en (0.2Ȝ, 0.2Ȝ)
x c) Posición (0.4Ȝ,0Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Técnica
Point- 2.5848e- 5.137437e-
16.1649 19.9669 0.3775 0.0026
matching 005 007
Optimización
2.0126e- 2.03818e- 4.24563612e-
del valor de 11.2488 3.734492 0.300117
003 005 007
las cargas
Tabla 3.16 Coste optimizando cargas para distintos nº de imágenes, posición fuente en (0.4Ȝ, 0Ȝ)
x c) Posición (0.4Ȝ,0.4Ȝ):
Nº cargas 5 10 15 20 25 30
Tipo
Técnica
Point- 26.6829 23.6945 1.4449 1.320842e- 2.8172e- 2.55236e-005
matching 002 004
Optimización 10.3198 7.17198 0.97371 9.925287e- 2.6741e- 2.55235e-005
del valor de 003 004
las cargas
Tabla 3.17 Coste optimizando cargas para distintos nº de imágenes, posición fuente en (0.4Ȝ, 0.4Ȝ)
-86-
Capítulo 3: Optimización de la función espacial de Green
x d) Posición (0.5Ȝ,0Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Técnica
Point- 9.6659e- 3.77679e-
17.9904 19.9656 0.9056 0.00700
matching 005 006
Optimización
7.816235e- 3.766583e-
del valor de 12.0109 5.811915 0.66110031 0.005482
005 006
las cargas
Tabla 3.18 Coste optimizando cargas para distintos nº de imágenes, posición fuente en (0.5Ȝ, 0Ȝ)
x e) Posición (0.7Ȝ,0.7Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Técnica
Point-
22.8267 8.0995 1.6061 4.2177 0.4510 1.1371
matching
Optimización
del valor de 22.6211220 8.03442 1.57554 4.214055 0.45094 1.137081
las cargas
Tabla 3.19 Coste optimizando cargas para distintos nº de imágenes, posición fuente en (0.7Ȝ, 0.7Ȝ)
-87-
Capítulo 3: Optimización de la función espacial de Green
Este método se basa en modificar las posiciones de las cargas imagen, con el fin
de optimizarlas para obtener de esta forma un coste mínimo.
Al resolver el sistema de ecuaciones definido en la ecuación 3.1 ó 3.4 (que son
la misma), obtenemos un conjunto de N soluciones complejas con las que podemos
recomponer el potencial en el interior de la estructura a analizar. Para el planteamiento
de este sistema, habíamos supuesto unas posiciones conocidas de las cargas imagen
(como por ejemplo, situarlas de forma especular a un plano tangente a cada punto del
cilindro).
En este método se propone optimizar la posición de las cargas imagen (en sus
ejes ‘x’ e ‘y’) de tal forma que las nuevas posiciones de las cargas calculadas
proporcionen un coste mínimo. Para el cálculo de estas nuevas posiciones emplearemos
un método de descenso de gradiente, de forma similar al utilizado para optimizar los
valores complejos de las cargas. No obstante, en este caso los valores a optimizar son
reales (y no complejos como el valor de las cargas) y además se deben optimizar dos
parámetros de forma conjunta (las posiciones en el eje ‘x’ y en el eje ‘y’).
La actualización de las posiciones de cada imagen se realizará mediante la
conocida fórmula del gradiente descendente:
x (k 1) x (k ) P w C (k ) (3.45)
y (k 1) y (k ) P w C (k )
Donde:
x P Æ Constante de paso.
-88-
Capítulo 3: Optimización de la función espacial de Green
imagen que minimicen el coste. Para completar el algoritmo, será necesario definirnos
un Coste derivable respecto a las posiciones; presentaremos dos versiones similares a
las mostradas en la optimización de las cargas.
C oste
Contorno
³ Gv cy dc (3.46)
De donde:
K
d ª º d
xiGvcy
dxi «
¬
Gv ( r , ro ) ¦
k 1
qk Gv(r , rk )»
¼
qi
dxi
Gv(r , ri ) (3.48)
-89-
Capítulo 3: Optimización de la función espacial de Green
El gradiente del coste empleado para optimizar las posiciones ‘y’ de las cargas
imagen se obtiene de una forma totalmente dual, resultando:
jk r r '
e ( y y ' ) (1 j k r r ' )
yi Coste qi ³ 3
dc (3.54)
Contorno r r'
En este punto definiremos el coste a optimizar como la integral del cuadrado del
valor absoluto de la función de Green en el contorno de la estructura, tal y como
observamos en la ecuación 3.55:
-90-
Capítulo 3: Optimización de la función espacial de Green
2
C oste ³ Gvcy dc
Contorno
(3.55)
K
§ ·
¨ Gv(r , ro ) Gv(r , ro )* Gv(r , ro ) ¦ qk*Gv(r , rk )* ¸
¨ k 1 ¸
Coste ³ ¨ K
ª K
º ª K
*º
¸dc (3.56)
¨ Gv(r , ro ) ¦ qk Gv(r , rk ) «¦ qk Gv(r , rk )» «¦ qk Gv(r , rk ) » ¸
Contorno ¨ * *
¸
© k 1 ¬k 1 ¼ ¬k 1 ¼¹
Calcularemos en este punto la derivada parcial del coste respecto a una posición
x de una carga imagen arbitraria ‘p’:
K
§ ·
¨ Gv(r , ro ) Gv(r , ro )* Gv(r , ro ) ¦ qk*Gv(r , rk )* ¸
w ¨ k 1 ¸
x p Coste ³
wx p Contorno ¨ K
ª K
º ª K
*º
¸dc (3.57)
¨ Gv(r , ro ) ¦ qk Gv(r , rk ) «¦ qk Gv(r , rk )» «¦ qk Gv(r , rk ) » ¸
* *
¨ ¸
© k 1 ¬k 1 ¼ ¬k 1 ¼¹
w *
wx p ³ Gv(r , r ) Gv(r , r ) dc
Contorno
o o 0 (3.58)
-91-
Capítulo 3: Optimización de la función espacial de Green
K K
w w
³ Gv(r , ro ) ¦ qk*Gv(r , rk )*dc ³ Gv(r , ro ) ¦ q Gv(r , r )
*
k k
*
wx p Contorno k 1 Contorno
wx p k 1
w
³ Gv(r , ro ) q*p Gv(r , rp )* (3.59)
Contorno
wx p
jk r r '
* e ( x x' ) (1 j k r r ' )
³ Gv(r , r ) q
Contorno
o p
r r'
3
dc
w
³ Gv* (r , ro ) q p Gv( r , rp ) (3.60)
Contorno
wx p
jk r r '
* e ( x x' ) (1 j k r r ' )
³ Gv (r , ro ) q p
Contorno r r'
3
dc
( x x' ) (1 j k r r ' ) º
K jk r r '
e
q*p (¦ qk Gv(r , rk )) 3 » dc
k 1 r r' »¼
Así, podemos obtener el gradiente del coste respecto la posición ‘x’ de una carga
imagen cualquiera, ‘p’, como:
-92-
Capítulo 3: Optimización de la función espacial de Green
De igual forma que hemos definido el coste respecto a las posiciones en ‘x’ de
las cargas imágenes, podemos definirlo respecto a las posiciones ‘y’; para ello usaremos
la ecuación 3.63 sin más que sustituir la función auxiliar ‘F’ por otra que tenga en
cuenta la posición en ‘y’ de las cargas. Esta función para las posiciones ‘y’ queda
definida en la ecuación 3.64:
& & jk r r '
& & dGv(r , r ' ) e ( y y ' ) (1 j k r r ' )
F (r , r ' ) 3
dy ' r r'
* & & * jk r r '
(3.64)
& & dGv ( r , r ' ) e ( y y ' ) (1 j k r r ' )
F * (r , r ' ) 3
dy ' r r'
-93-
Capítulo 3: Optimización de la función espacial de Green
-94-
Capítulo 3: Optimización de la función espacial de Green
Figura 3.33 Potencial escalar en función del ángulo. 10 imágenes y fuente en (0.5Ȝ ,0.5Ȝ)
-95-
Capítulo 3: Optimización de la función espacial de Green
Criterio de parada primando velocidad (a) Criterio de parada primando coste (b)
Figura 3.34 Representación de la evolución del coste y de las posiciones de las cargas imagen al
realizar una optimización de las posiciones. 10 imágenes y fuente en (0.5Ȝ ,0.5Ȝ)
-96-
Capítulo 3: Optimización de la función espacial de Green
de limitar los valores de pico en la estructura, mientras que el resto se sitúan alrededor
de la cavidad, en una situación más alejada.
Con esta última configuración, obtenemos una distribución del potencial en el
contorno de la estructura circular de:
Figura 3.35 Potencial escalar en función del ángulo tras optimizar las posiciones de las cargas
imagen. 10 imágenes y fuente en (0.5Ȝ ,0.5Ȝ)
x Posición (0Ȝ,0Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Técnica
7.1487e- 1.3369e- 3.07174e-
Point-matching 10.1294 1.3525 0.0061
005 006 008
6.0212e- 7.0595e- 1.3369e-
Optimización de las 1.0771 3.027e-
10.1061 003 005 006
posiciones 008
Tabla 3.20 Coste optimizando posiciones para distintos nº de imágenes, fuente en (0.0Ȝ , 0.0Ȝ)
-97-
Capítulo 3: Optimización de la función espacial de Green
x Posición (0.2Ȝ,0.2Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Técnica
5.6664e- 5.9574e- 1.204944e-
Point-matching 12.1475 16.5825 0.0972
004 006 007
7.8030E- 4.4101e- 4.5577e- 1.012308e-
Optimización de 11.8225 9.21682
002 004 006 007
las posiciones
Tabla 3.21 Coste optimizando posiciones para distintos nº de imágenes, fuente en (0.2Ȝ, 0.2Ȝ)
x Posición (0.4Ȝ,0Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Técnica
Point- 2.5848e- 5.137437e-
16.1649 19.9669 0.3775 0.0026
matching 005 007
Optimización 2.19611e- 2.071509E- 4.587410e-
14.5984 10.9457 0.30708
de las 003 005 007
posiciones
Tabla 3.22 Coste optimizando posiciones para distintos nº de imágenes, fuente en (0.4Ȝ, 0Ȝ)
x Posición (0.4Ȝ,0.4Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Técnica
Point- 1.320842e- 2.8172e-
26.6829 23.6945 1.4449 2.55236e-005
matching 002 004
Optimización 1.23631e- 2.80277e-
19.0104 9.920484 1.07586 2.55236e-005
de las 002 004
posiciones
Tabla 3.23 Coste optimizando posiciones para distintos nº de imágenes, fuente en (0.4Ȝ, 0.4Ȝ)
x Posición (0.5Ȝ,0Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Técnica
Point- 9.6659e- 3.77679e-
17.9904 19.9656 0.9056 0.00700
matching 005 006
Optimización 5.588779e- 8.241987e- 3.776069e-
17.0055 9.167151 0.711393
de las 003 005 006
posiciones
Tabla 3.24 Coste optimizando posiciones para distintos nº de imágenes, fuente en (0.5Ȝ, 0Ȝ)
-98-
Capítulo 3: Optimización de la función espacial de Green
x Posición (0.7Ȝ,0.7Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Técnica
Point-
22.8267 8.0995 1.6061 4.2177 0.4510 1.1371
matching
Optimización 20.8588 3.35552 0.63364 1.473625 0.14394 0.32348
de posiciones
Tabla 3.25 Coste optimizando posiciones para distintos nº de imágenes, fuente en (0.7Ȝ, 0.7Ȝ)
-99-
Capítulo 3: Optimización de la función espacial de Green
Mostraremos a continuación cual es el valor final del coste obtenido por los
distintos métodos (tanto los comentados en apartados anteriores como la combinación
de los mismos) en función del número de cargas imagen empleadas (5, 10, 15, 20, 25,
30) y de la posición de la fuente en el interior del cilindro [ (0Ȝ,0Ȝ) , (0.2Ȝ,0.2Ȝ) ,
(0.4Ȝ,0Ȝ) , (0.4Ȝ,0.4Ȝ) , (0.5Ȝ,0Ȝ) , (0.7Ȝ,0.7Ȝ) ]. Indicaremos en rojo aquellos resultados
que mejoren respecto al valor inicial de cargas imágenes.
x a) Posición (0Ȝ,0Ȝ):
Nº
cargas
5 10 15 20 25
Tipo
Técnica
Point-
10.1294 1.3525 0.0061 7.1487e-005 1.3369e-006
matching
Solo varía 9.6945 1.12238 1.33286-003 2.56416E-007 1.34045E-011
radio
Solo
optim. 10.1293 0.89803 3.8553E-003 4.5499E-005 8.63487E-007
cargas
Variación
radio +
9.6945 0.7715838 8.48775E-004 1.63088E-007 8.5839E-012
optimiza
cargas
Optimiza
10.1061 1.07716 6.0212E-003 7.05951e-005 1.3392e-006
posiciones
Optimiza
4.549704e-
posiciones 10.1061 0.893449 3.8555E-003 8.634874e-007
005
y cargas
Tabla 3.26 Comparativa distintos métodos de optimización en función del número de imágenes,
fuente en (0.0Ȝ , 0.0Ȝ)
-100-
Capítulo 3: Optimización de la función espacial de Green
x b) Posición (0.2Ȝ,0.2Ȝ):
Nº
cargas
5 10 15 20 25 30
Tipo
Técnica
Point- 5.6664e- 5.9574e- 4.0841e-
12.1475 16.5825 0.0972
matching 004 006 008
Solo varía 3.2547 3.5181e- 3.8519e- 4.8903e- 1.9586e-
10.9611
radio 002 005 007 008
Solo
2.4473 8.0247e- 4.4086e- 4.5437e- 8.9950e-
optim. 10.7727
002 004 006 008
cargas
Variación
2.7858e- 3.5059e- 4.8902e-
radio + 10.4712 2.4647 1.9586e-
002 005 007
optimiza 008
cargas
7.8030e- 4.4101e- 4.5577e- 1.0123e-
Optimiza 11.8225 9.2168
002 004 006 007
posiciones
Optimiza
7.6710e- 4.4044e- 4.5415e- 8.9952e-
posiciones 10.7593 3.9511
002 004 006 008
y cargas
Tabla 3.27 Comparativa distintos métodos de optimización en función del número de imágenes,
fuente en (0.2Ȝ , 0.2Ȝ)
x c) Posición (0.4Ȝ,0Ȝ):
Nº
cargas
5 10 15 20 25 30
Tipo
Técnica
Point- 2.5848e- 4.6492e-
16.1649 19.9669 0.3775 0.0026
matching 005 007
6.8721
Solo varía 10.3900 8.6191e- 5.8619e- 1.8183e- 3.2837e-
radio 002 004 005 007
Solo
2.0126e- 2.0381e- 4.2456e-
optim. 11.2488 3.7344 0.3001
003 005 007
cargas
Variación
6.8744e- 4.8472e- 1.7717e- 3.0937e-
radio + 10.0980 4.8108
002 004 005 007
opti
cargas
-101-
Capítulo 3: Optimización de la función espacial de Green
x d) Posición (0.4Ȝ,0.4Ȝ):
Nº
cargas
5 10 15 20 25 30
Tipo
Técnica
Point- 1.3208e- 2.8172e-
26.6829 23.6945 1.4449 2.5523e-005
matching 002 004
Solo varía 11.0951 1.3208e- 2.8172e-
9.0305 0.2952 2.5523e-005
radio 002 004
Solo optim. 9.9252e- 2.6741e-
10.3198 7.1719 0.9737 2.5523e-005
cargas 003 004
Variación
9.9252e-
radio + 10.4546 7.8761 0.2371 2.6741e-
003 2.5523e-005
optimización 004
cargas
1.2363e- 2.8027e-
Optimiza 19.0104 9.9204 1.0758 2.5523e-005
002 004
posiciones
Optimiza
0.8390 1.0916e- 2.6748e-
posiciones y 10.8615 6.5178 2.5523e-005
002 004
cargas
Tabla 3.29 Comparativa distintos métodos de optimización en función del número de imágenes,
fuente en (0.4Ȝ , 0.4Ȝ)
x e) Posición (0.5Ȝ,0Ȝ):
Nº
cargas
5 10 15 20 25 30
Tipo
Técnica
Point- 9.6659e-
17.9904 19.9656 0.9056 0.0070 3.7767e-006
matching 005
Solo varía 10.1459 8.5620 4.6655e- 9.6659e-
0.1721 3.7767e-006
radio 003 005
Solo 7.8162e-
5.4822e-
optim. 12.0109 5.8119 0.6611 005 3.7665e-006
003
cargas
-102-
Capítulo 3: Optimización de la función espacial de Green
Variación
3.7698e- 7.8162e-
radio + 9.9387 6.6690 0.1358
003 005 3.7665e-006
opti
cargas
5.5887e- 8.2419e-
Optimiza 17.0058 9.1671 0.7113 3.7760e-006
003 005
posiciones
Optimiza
5.1803 0.5735 5.4893e- 7.9875e-
posiciones 11.4009 3.7695e-006
003 005
y cargas
Tabla 3.30 Comparativa distintos métodos de optimización en función del número de imágenes,
fuente en (0.5Ȝ , 0.0Ȝ)
x f) Posición (0.7Ȝ,0.7Ȝ):
Nº
cargas
5 10 15 20 25 30
Tipo
Técnica
Point-
22.8267 8.0995 1.6061 4.2177 0.4510 1.1371
matching
Solo varía 22.8266 8.0995 1.6061 4.2176 0.4509 1.1370
radio
Solo
22.6211 8.0344 1.5755 4.2140
optim. 0.4509 1.1370
cargas
Variación
radio + 22.6211 8.0344 1.5755 4.2140
0.4509 1.1370
opti
cargas
OptimCiza 20.8588 3.3555 0.6336 1.4736 0.1439 0.3234
posiciones
Optimiza
posiciones 4.0479 2.3237 0.6335 1.098 0.1401 0.3234
y cargas
Tabla 3.31 Comparativa distintos métodos de optimización en función del número de imágenes,
fuente en (0.7Ȝ , 0.7Ȝ)
-103-
Capítulo 3: Optimización de la función espacial de Green
3.2.8 Conclusiones
-104-
Capítulo 3: Optimización de la función espacial de Green
el del resto, a pesar de que ha sido controlado y moderado con técnicas como una
gestión dinámica del paso de actualización o una detección automática del algoritmo.
El último de los métodos presentados consistió en una combinación de técnicas,
se basaba en variar primero la posición de las cargas imágenes (mediante una técnica de
variación del radio u optimización de las posiciones) y luego optimizar el valor
complejo de las cargas. Esta técnica compuesta ha proporcionado los mejores resultados
en todas las situaciones, disminuyendo el error de forma sustancial respecto al
proporcionado por Point-Matching. Además presenta la ventaja de no depender del tipo
de estructura que se está analizando, siendo fácilmente extensible a geometrías
arbitrarias.
-105-
Capítulo 3: Optimización de la función espacial de Green
i 1,2,3...N
-106-
Capítulo 3: Optimización de la función espacial de Green
N N
& & & & & &
sin(Mi ) ¦ GAxx ri , rk ´ I kx sin(Mi ) G Axx ri , r0´ cos(Mi ) ¦ G Ayy ri , rk ´ I ky 0
k 1 k 1
N N
(3.67)
cos(Mi ) ¦ C I cos(Mi ) C sin Mi ¦ C I
x x x y y
i ,k k 1, 0 i ,k k 0
k 1 k 1
N N
cos(M c ) ¦ Ccx,k I k cos(Mi ) C1x,0 sin Mc ¦ Ccy,k I k dc
x y
F2 ³
Contorno k 1 k 1
0 (3.69)
De tal forma que cuanto más próximo sea el resultado de las ecuaciones 3.68 y
3.69 a cero más preciso será la representación del potencial vector en el interior de la
estructura analizada. Finalmente, sumaremos ambas condiciones con el fin de obtener
un coste total representativo del potencial vector:
Coste F1 F2 (3.70)
Una vez definido el coste podemos aplicarlo a distintas situaciones con el fin de
evaluar la precisión del potencial vector en cada caso.
Concretamente, el valor del coste dependerá del número de imágenes empleadas
(pues cuantas más imágenes usemos, un mayor número de puntos discretos cumplirán la
condición de contorno) y de la posición de la carga fuente en el interior del cilindro.
Supongamos un cilindro de radio a O , (siendo la frecuencia de 300Mhz) en el
que situamos la carga fuente en el interior del cilindro y empleamos un total de 20
imágenes para calcular el potencial escalar. La situación es representada en la
figura 3.36:
Figura 3.36 Representación 20 cargas imágenes alrededor del cilindro, fuente en origen
coordenadas
-107-
Capítulo 3: Optimización de la función espacial de Green
Figura 3.37 Valor de las condiciones del potencial vector en el contorno del cilindro, 20 imágenes,
fuente situada en el origen.
Figura 3.38 Distribución del error del potencial vector en el contorno del cilindro, 20 imágenes,
fuente situada en el origen.
-108-
Capítulo 3: Optimización de la función espacial de Green
Por ejemplo, basándonos en las suposiciones del caso anterior, pero situando la
carga en la posición (0.6Ȝ, 0.6 Ȝ) obtenemos una situación y una distribución del coste
como muestra la figura 3.39:
Figura 3.39 Representación de situación y coste del potencial escalar en un cilindro, empleo de 20
imágenes y fuente situada en (0.6 Ȝ,0.6 Ȝ)
Ello implica que la representación del potencial vector en este caso es mucho
menos exacta que en el caso anterior, no obstante todavía no es un valor muy elevado y
se puede dar por válido.
Podemos observar la descomposición del Coste anterior, comprobando cual de
las dos condiciones impuestas sobre el contorno del cilindro presenta un mayor error:
Figura 3.40 Detalle del Coste del potencial vector para cada condición. 20 imágenes y fuente situada
en (0.6 Ȝ,0.6 Ȝ)
-109-
Capítulo 3: Optimización de la función espacial de Green
Nº cargas 5 10 15 20 25 30
Tabla 3.32 Coste vectorial método original, posición fuente en (0Ȝ ,0Ȝ)
Nº cargas 5 10 15 20 25 30
Tabla 3.33 Coste vectorial método original, posición fuente en (0.2Ȝ ,0.2Ȝ)
x c) Posición (0.4Ȝ,0Ȝ):
Nº cargas 5 10 15 20 25 30
Tabla 3.34 Coste vectorial método original, posición fuente en (0.4Ȝ ,0Ȝ)
-110-
Capítulo 3: Optimización de la función espacial de Green
x d) Posición (0.4Ȝ,0.4Ȝ):
Nº cargas 5 10 15 20 25 30
Tabla 3.35 Coste vectorial método original, posición fuente en (0.4Ȝ ,0.4Ȝ)
x e) Posición (0.5Ȝ,0Ȝ):
Nº cargas 5 10 15 20 25 30
Tabla 3.36 Coste vectorial método original, posición fuente en (0.5Ȝ ,0Ȝ)
x f) Posición (0.7Ȝ,0.7Ȝ):
Nº cargas 5 10 15 20 25 30
Tabla 3.37 Coste vectorial método original, posición fuente en (0.7Ȝ ,0.7Ȝ)
-111-
Capítulo 3: Optimización de la función espacial de Green
El método de la variación del radio del circulo en el que situaremos los dipolos
imagen ha sido implementado debido a la excelentes resultados que proporcionó la
versión desarrollada para el potencial escalar.
Al igual que en éste, el método se basa en una ubicación de los dipolos imagen
diferente a la empleada hasta ahora en el método general (expuesto en el capítulo 2),
manteniendo los dipolos imágenes en un circulo alrededor de la estructura a analizar
(independientemente de la forma de la estructura y de la posición de la fuente en el
interior de la misma).
Supongamos una situación de 20 dipolos imágenes, donde hemos situado el
dipolo fuente en una ubicación distinta del origen, por ejemplo en (0.6Ȝ ,0.6Ȝ); podemos
observar las diferencias entre el método general de ubicación de dipolos y el actual:
Observamos como al situar los dipolos imagen mediante la variación del radio,
éstos se encuentran distribuidos uniformemente a lo largo de un círculo que envuelve a
la estructura a analizar, sin depender de la posición del dipolo fuente (situado en el
interior de la estructura) como pasaba en el método original.
-112-
Capítulo 3: Optimización de la función espacial de Green
cálculo de las distancias de los dipolos imágenes al centro del cilindro, a la posición de
la fuente y a los puntos tangentes en el contorno del cilindro.
Figura 3.42 Variación del coste del potencial vector en función del radio. Empleo de 20 dipolos
imagen y fuente situada en (0.6Ȝ ,0.6Ȝ)
-113-
Capítulo 3: Optimización de la función espacial de Green
Figura 3.43 Variación del coste del potencial vector en función del radio en torno al mínimo.
Empleo de 20 dipolos imagen y fuente situada en (0.6Ȝ ,0.6Ȝ)
a) b)
Figura 3.44 Comparación de la distribución de la condición del potencial vector en el contorno de
un cilindro. a) Optimizando el radio. b) Método general. Empleo de 20 dipolos imagen y fuente
situada en (0.6Ȝ ,0.6Ȝ)
-114-
Capítulo 3: Optimización de la función espacial de Green
x a) Posición (0Ȝ,0Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Técnica
1.3035e- 2.8246e- 7.3812e-
Point-matching 34.9888 14.9860 0.1174
3 5 7
4.7608e- 2.4342e- 2.6383e- 1.4201e-
Error 23.9857 6.4980
2 5 9 13
Variación
1.210Ȝ 7.500 Ȝ
del radio 30.00 Ȝ 30.00 Ȝ 30.00 Ȝ 30.00 Ȝ
Radio Mínimo Mínimo
Converge Converge Converge Converge
local local
Tabla 3.38 Coste vectorial variación radio para distintos nº de cargas, posición fuente en (0Ȝ ,0 Ȝ)
Figura 3.45 Ejemplo de convergencia del coste del potencial vector en el infinito. 15 imágenes,
fuente en el origen de coordenadas.
-115-
Capítulo 3: Optimización de la función espacial de Green
x b) Posición (0.2Ȝ,0.2Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Técnica
1.5561e- 1.5819e- 3.7200e-
Point-matching 31.1104 31.9765 1.3449
2 4 6
2.2103e- 2.7185e- 1.0271e-
Error 22.7601 16.6434 0.7416
3 5 6
Variación
1.280 Ȝ 2.0700 Ȝ 3.1900 Ȝ 2.500 Ȝ 2.2800 Ȝ
del radio 30.00 Ȝ
Radio Mínimo Mínimo Mínimo Mínimo Mínimo
Converge
local local local local local
Tabla 3.39 Coste vectorial variación radio para distintos nº de cargas, posición fuente en (0.2Ȝ ,0.2Ȝ)
x c) Posición (0.4Ȝ,0.0Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Técnica
6.3463e-
Point-matching 64.0995 62.1398 3.6035 0.1014 1.1809e-5
4
1.3308e- 7.3222e-
Error 24.5626 18.3595 1.5567 6.2745e-5
2 4
Variación
1.230 Ȝ 1.230 Ȝ 1.7600 Ȝ 3.3800 Ȝ 2.090 Ȝ 2.3700 Ȝ
del radio
Radio Mínimo Mínimo Mínimo Mínimo Mínimo Mínimo
local local local local local local
Tabla 3.40 Coste vectorial variación radio para distintos nº de cargas, posición fuente en (0.4Ȝ ,0Ȝ)
Figura 3.46 Variación de la condición del potencial vector en el contorno del cilindro en función del
radio, existiendo 2 mínimos. 15 cargas imágenes y fuente en (0.4Ȝ ,0.0Ȝ)
-116-
Capítulo 3: Optimización de la función espacial de Green
x d) Posición (0.4Ȝ,0.4Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Técnica
1.1148e-
Point-matching 57.1065 90.6108 10.4649 0.7906 4.2507e-4
2
2.3337e-
Error 25.3086 18.5186 3.6720 0.2095 9.277e-3
2
Variación
1.190 Ȝ 1.270 Ȝ 1.630 Ȝ 1.7800Ȝ 1.7400Ȝ 1.650Ȝ
del radio
Radio Mínimo Mínimo Mínimo Mínimo Mínimo Mínimo
local local local local local local
Tabla 3.41 Coste vectorial variación radio para distintos nº de cargas, posición fuente en (0.4Ȝ, 0.4Ȝ)
Comprobamos como en esta posición hay dos casos en los que una variación del
radio no proporciona un coste menor que el original. Esto se produce para un número de
imágenes relativamente elevado (25-30) cuando el error ya es reducido. Comprobamos
como conforme se va acercando el dipolo fuente a la pared, la variación del radio va
proporcionando unos resultados más modestos al usar un número de dipolos imagen
elevado.
x e) Posición (0.5Ȝ,0.0Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Técnica
3.1735e-
Point-matching 70.9312 83.0511 7.2793 0.3847 7.6270e-5
3
6.1674e- 5.9330e-
Error 24.0265 18.2668 2.6496 8.4696e-4
2 3
Variación
1.240 Ȝ 1.240 Ȝ 1.630 Ȝ 2.410 Ȝ 1.940 Ȝ 2.160 Ȝ
del radio
Radio Mínimo Mínimo Mínimo Mínimo Mínimo Mínimo
local local local local local local
Tabla 3.42 Coste vectorial variación radio para distintos nº de cargas, posición fuente en (0.5Ȝ, 0.0Ȝ)
x f) Posición (0.7Ȝ,0.7Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Técnica
Point-matching 167.0453 159.1433 85.0425 245.5419 535.3140 378.9134
Error 158.192 156.153 151.083 144.478 144.8139 144.3359
Variación 1.140 Ȝ 1.100 Ȝ 1.050 Ȝ 1.410 Ȝ 1.040 Ȝ 1.160 Ȝ
del radio Radio Mínimo Mínimo Mínimo Mínimo Mínimo Mínimo
local local local local local local
Tabla 3.43 Coste vectorial variación radio para distintos nº de cargas, posición fuente en (0.7Ȝ, 0.7Ȝ)
-117-
Capítulo 3: Optimización de la función espacial de Green
Podemos observar como en casi todos los casos obtenemos un coste mucho
menor que con la técnica original, incluyendo la posición (0.7Ȝ, 0.7Ȝ). No obstante a
pesar de mejorar los resultados originales, los costes obtenidos son muy elevados,
dando lugar a aproximaciones del potencial vector relativamente imprecisas. Ello es
debido a que el dipolo fuente se encuentra muy próximo a la pared del cilindro, por lo
que produce un pico de potencial en ella que los dipolos imágenes de nuestro método, al
estar situadas de forma simétrica alrededor del cilindro, no pueden contrarrestar.
Al igual que ocurría con este método al optimizar el potencial escalar, el error se
disminuye de una forma más drástica cuando el número de cargas imágenes es reducido
(ó el coste inicial de partida muy elevado) y de una forma más moderada cuando el
número de cargas imágenes es más elevado.
Siempre existirá un coste mínimo, que se producirá en un mínimo local o cuando
el algoritmo converja a un valor determinado.
Como principal desventaja de este método, señalaremos que el coste
computacional para conseguir los resultados es relativamente alto, pues debemos de
calcular los potenciales para un gran número de situaciones y obtener el Coste para cada
uno de ellos (operación computacionalmente costosa)
Para concluir, podemos considerar la técnica como muy satisfactoria pues es
muy efectiva en la gran mayoría de los casos, simple y fácilmente extensible a una
geometría arbitraria.
3.3.3 Conclusiones
-118-
Capítulo 3: Optimización de la función espacial de Green
-119-
Capítulo 3: Optimización de la función espacial de Green
-120-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
4.1 Introducción
-121-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
optimización por gradiente del valor complejo de las cargas y de las posiciones. De esta
forma, tratamos de obtener la representación de los potenciales con el menor error
posible.
Finalmente, realizaremos una comparación del potencial en el interior de una
gruía cuadrada con una altura en el eje ‘z’ teóricamente infinita con el potencial
calculado con el presente método (que considera un corte en ‘z’ infinitesimal). A pesar
de que físicamente no estamos representando la misma magnitud, obtenemos unas
similitudes muy importantes, dando a pie a la introducción en el capítulo 5 de guías en
3D (y no en un una estructura bidimensional como hemos analizado hasta la fecha).
-122-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
a) Potencial escalar
Figura 4.2 Ubicación de una condición de contorno Gv en la pared de una guía cuadrada
& &
Siendo Gv r , r ´ el potencial de una carga en condiciones de espacio libre:
& &
& & P 0 e jko r r ´
Gv r , r ´ & & (4.4)
4S r r ´
-123-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
b) Potencial vector
-124-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
&
E ptoi eˆt 0 (4.7)
-125-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
& &
2) A 0
wA x wA y
dx dx
wA y (4.18)
wA x
dy dy
-126-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
dAx
Este desarrollo es totalmente dual para , obteniendo un resultado:
dy
j*k * r r '
wA x P0 e * ( y y ' ) * (1 j * k * r r ' )
3 (4.23)
dy 4 S r r'
En este punto tenemos totalmente definida la ecuación 4.14, que nos permitía
cumplir la condición 4.9.
Una vez analizadas las dos condiciones para un punto cualquiera ‘i’ tangente a
la estructura, vamos a extender la formulación a un total de N puntos, con el fin de que
las condiciones impuestas se verifiquen en el conjunto de la estructura.
-127-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
Para ello analizaremos nuevamente las condiciones 4.8 y 4.9, extendiendo sus
resultados para que satisfagan las condiciones en todos los puntos de la estructura. En
este análisis, consideraremos que existe un dipolo en el interior, orientado a lo largo del
eje ‘x’. La formulación para el caso de que esta orientación sea en el eje ‘y’ es
totalmente dual.
&
1) eˆn u A 0
-128-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
i 1,2,3...N
-129-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
Figura 4.5 Estructura cilíndrica, de radio a O , analizada con 20 imágenes y fuente situada en el
centro. Perspectivas de formulación cartesiana y formulación cilíndrica.
-130-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
Figura 4.6 Comparación del potencial escalar en un cilindro bajo las formulaciones cartesiana y
cilíndrica. 20 imágenes, fuente en el origen.
-131-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
Figura 4.7 Comparación del potencial vector en un cilindro bajo las formulaciones cartesiana y
cilíndrica. 20 imágenes, fuente en el origen.
-132-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
-133-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
Guía cuadrada con dos lados redondeados Guía cuadrada con bordes redondeados
-134-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
-135-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
0
3 4
0
4 5
1
5 6 2 0.5 1
0
6 7
1
7 8 3 0.5 1
0
8 9
0
9 10
0
10 11
1
11 12 4 0.5 1
0
12 1
! Finalmente, indicaremos la posición de la carga fuente.
0 0
Figura 4.9 Fichero descriptor de una guía ridge con sus lados redondeados
Donde observamos que hemos resaltado los vértices para una mejor
visualización de las coordenadas.
Mostraremos a continuación como es el formato del fichero “Geom.in”,
construyendo paso a paso la geometría anterior.
Para empezar a describir la guía, deberemos de indicar el número de vértices
totales de que consta:
! Fichero de configuración de la geometria de la cavidad
------------------------------------------------
! Número de vértices
12
-136-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
Después, indicaremos la posición de los centros desde los cuales vamos a trazar
los tramos curvos, teniendo en cuenta que volveremos a referenciarlos según su
posición:
! Definiremos los centros de los tramos curvos:
1 1 Æ Define el centro nº 1
-1 1 Æ Define el centro nº 2
-1 -1
1 -1
Finalmente, indicaremos como vamos a unir los distintos vértices, distinguiendo
en función de que sean tramos rectos ó curvos:
a) Tramos recto:
-137-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
Por último indicar que todas las medidas introducidas se encuentran metros (m),
y que esta es la única unidad de medida que permite el software desarrollado.
-138-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
Figura 4.11 Distribución uniforme de 15 puntos tangentes en una guía ridge “redondeada”
Observamos como en este caso tras empezar a ubicar los puntos tangentes desde
el vértice superior derecho, ningún otro vértice recibe un punto tangente.
-139-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
a) b)
c)
Figura 4.13 Ubicación de 3,4 y 7 puntos tangentes en una guía triangular, empleando una
distribución de los mismos por vértices y mitades de segmentos.
Aplicando este concepto a una guía ridge con los lados redondeados, con un
total de 32 puntos tangentes, tendremos una ubicación de los mismos de:
Figura 4.14 Ubicación de 32 puntos tangentes en una guía ridge redondeada, empleando una
distribución de los mismos por vértices y mitades de segmentos.
-140-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
En este caso la táctica a seguir será asignar a cada tramo un número de puntos
tangentes y luego distribuirlos en él. La forma de realizar esta asignación tendrá dos
componentes:
o Un número de puntos tangentes fijo (por ejemplo del 70%, aunque es
modificable) que se va a distribuir de forma uniforme por la
estructura, de forma que a todos los tramos le corresponda el mismo
número de puntos tangentes.
o El resto, se asignará en función de la distancia media de cada tramo a
la posición de la carga fuente, de forma inversamente proporcional.
Así cuanto más cerca esté un tramo de la carga fuente, más puntos
tangentes se le asignarán.
Figura 4.15 Ubicación de puntos tangentes en guía ridge redondeada en función de la distancia
media de cada tramo a la fuente
-141-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
Figura 4.16 Ubicación de 20 puntos de forma uniforme y saltando los vértices, sobre una guía
rectangular.
Una vez que tenemos totalmente definida la estructura y ubicados los puntos
tangentes donde forzaremos que se cumplan las condiciones de contorno, el siguiente
paso será situar de forma adecuada las cargas/dipolos imágenes alrededor.
Impondremos el mismo número de cargas imágenes que puntos tangentes se hayan
definido, tal y como lo requiere la formulación.
En este apartado mostraremos aquellos tipos de ubicación de imágenes estáticos,
es decir, que no van a cambiar su posición en función del “coste” de los potenciales.
De entre las distintas posibilidades de situar las imágenes, se han
implementando las siguientes:
a) Situación circular
En este primer caso, vamos a situar las cargas imágenes de forma circular
alrededor de la estructura arbitraria, de forma independiente a como se hayan
situado los puntos tangentes en el contorno.
El radio del círculo sobre el que se sitúan las cargas imágenes se calcula a
partir de la distancia máxima desde el centro hasta la estructura, multiplicándolo por
un factor corrector.
Mostramos a continuación dos ejemplos con una guía ridge “redondeada” y
una guía rombo:
-142-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
b) Posicionamiento especular
-143-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
Figura 4.20 Combinación de técnicas de ubicación de imágenes: especular y recta tangente auxiliar
para tramos curvos. Ejemplo sobre guía cuadrada con 2 lados redondeados.
-144-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
Uno de los problemas que se nos plantea a la hora de trabajar con geometrías
arbitrarias es saber si un punto cualquiera se encuentra en el interior o en el exterior de
la estructura, para poder así determinar si procede calcular el potencial en ese punto o
no.
Es por tanto necesario definirse una estrategia que nos permita conocer con
exactitud la ubicación de cada punto del espacio respecto a nuestra estructura.
Para ello, hemos implementado el siguiente algoritmo:
3) Consideraremos todos los tramos rectos en los que al menos uno de sus
vértices se encuentre en la dirección y sentido que une nuestro punto arbitrario
con el rectángulo, estando los de la dirección contraria uno por arriba y otro por
abajo:
Figura 4.22 Reconocimiento de tramos rectos situados entre un punto cualquiera y el exterior de la
estructura. Ejemplo con una guía cuadrada.
-145-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
Figura 4.23 Reconocimiento de tramos oblicuos situados entre un punto cualquiera y el exterior de
la estructura. Ejemplo con una guía triangular.
5) Otro asunto algo más complejo es cuando nos encontramos ante tramos
curvos. En este caso, hallaremos en primer lugar la distancia desde nuestro
punto arbitrario hasta el centro que define el tramo curvo. Si esta distancia es
mayor que el radio, comprobaremos si en nuestra dirección y sentido de
movimiento vamos a tener delante al menos uno de los vértices del tramo curvo.
En caso afirmativo de alguna de las dos condiciones anteriores, pasaremos a
estudiar más en detalle el caso:
o En primer lugar, comprobaremos si tenemos por delante los dos
vértices del tramo curvo, en cuyo caso es seguro que lo deberemos de
atravesar.
o En caso contrario, discerniremos si estamos en el interior del tramo
curvo o no en función de los ángulos. Para ello, hallaremos el ángulo
desde el centro del tramo curvo a cada uno de los vértices, y desde
nuestro punto arbitrario al centro del cilindro, comprobando
posteriormente si el ángulo respecto al punto arbitrario se encuentra
comprendido entre los ángulos de los vértices.
De esta forma, podemos analizar cualquier punto del espacio pudiendo distinguir
si pertenece o no al interior de cualquier estructura arbitraria, por compleja que sea. Esto
es muy importante a la hora de realizar el cálculo del potencial escalar en el interior de
nuestra estructura, pues ya tenemos definido completamente cual será ese interior.
Mostraremos a continuación algunos ejemplos de cómo el software programado
es capaz de determinar si un punto del espacio está en el interior de la estructura o no:
-146-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
Figura 4.24 Ejemplo de reconocimiento del interior de estructuras arbitrarias. Visualización en 3D.
1) Cuadrada
-147-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
-148-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
2) Triangular
-149-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
3) Rombo
-150-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
4) Ridge
Para analizar la guía tipo ridge, al ser algo más compleja, emplearemos un total
de 60 imágenes distribuidas de forma que sigan el contorno de la guía, relativamente
cerca de la misma. Al no existir tramos oblicuos ni curvos, las componentes del
potencial vector cruzadas son nulas.
Los potenciales obtenidos fueron los siguientes:
-151-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
5) Estrella
-152-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
-153-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
Para concluir este apartado debemos de indicar que los potenciales en el interior
de cada guía dependen de una forma muy significativa de cómo se ubiquen los puntos
tangentes en el contorno de la misma así como de la situación de las imágenes.
Pequeñas variaciones de estos elementos sobre algunos tipos de guías puede dar lugar a
resultados distintos. Es por esto que se necesita una medida para evaluar la corrección
de estos potenciales: el “coste” o error que cometemos al calcularlos. Gracias a este
“coste” podremos saber cuando los potenciales obtenidos son correctos y con qué
margen de error lo son. Además, nos permitirá analizar los puntos en los que se cometen
los errores en las implementaciones con el fin de poder atacarlos directamente para
reducirlos.
-154-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
Figura 4.30 Representación del potencial escalar en el contorno de una guía cuadrada. 40 imágenes,
fuente en el centro de la guía.
-155-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
Otro ejemplo sobre una guía más compleja puede ser la guía tipo “ridge”. Si
empleamos un total de 60 imágenes distribuidas siguiendo el contorno de la guía a una
distancia de 0.25O (una distribución diferente a la empleada en la visualización de los
potenciales en el apartado anterior), el potencial escalar en el contorno de la estructura
será:
Figura 4.31 Representación del potencial escalar en el contorno de una guía ridge. 60 imágenes,
fuente en el centro de la guía.
Figura 4.32 Representación del potencial escalar en el interior de una guía ridge. 60 imágenes
situadas a 0.25Ȝ siguiendo el contorno de la guía.
-156-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
El coste total de la implementación del potencial escalar para esta guía ridge es
de: Coste 17.2302 , mientras que el coste de la implementación mostrada en el
apartado anterior fue de 7.58 (un coste mucho más reducido provocado por una
separación diferente de las imágenes a la estructura).
Comprobamos como para el caso actual el coste es elevado estando propiciado
casi exclusivamente por lo picos de las cuatro esquinas superiores interiores. No
obstante el error es muy localizado, siendo el potencial en el interior de la guía correcto.
i 1,2,3...N
-157-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
N N
& & & & & &
sin(Mi ) ¦ G Axx ri , rk ´ I kx sin(Mi ) G Axx ri , r0´ cos(Mi ) ¦ G Ayy ri , rk ´ I ky 0
k 1 k 1
N N
(4.36)
cos(Mi ) ¦ C I cos(Mi ) C sin Mi ¦ C I
x x x y y
i ,k k 1, 0 i ,k k 0
k 1 k 1
N N
& & & & & &
F1 ³ sin(Mc ) ¦ G Axx rc , rk ´ I kx sin(M c ) G Axx rc , r0 ´ cos(M c ) ¦ G Ayy rc , rk ´ I ky dc
Contorno k 1 k 1 (4.37)
N N
cos(M c ) ¦ Ccx,k I k cos(Mi ) C1x,0 sin Mc ¦ Ccy,k I k dc
x y
F2 ³
Contorno k 1 k 1
0 (4.38)
Así, cuanto más reducido sea el valor de las ecuaciones 4.37 y 4.38 más fiel será
la representación del potencial vector en el interior de la estructura. Finalmente,
sumaremos ambas condiciones con el fin de obtener un coste total representativo del
potencial vector:
Coste F1 F2 (4.39)
A pesar de que las ecuaciones son muy similares a las presentadas en el capítulo
3 para las guías circulares, existe una diferencia fundamental: la evaluación de las
constantes definidas en la ecuación 4.35.
Tras definir el coste del potencial vector, ya podemos aplicarlo a las distintas
estructuras obteniendo así una medida del grado de corrección de los potenciales. Al
igual que en el caso del potencial escalar, presentamos la distribución de este coste para
cada uno de los tramos que constituyan la geometría analizada, distinguiendo además
entre las dos condiciones (4.37 y 4.38) y el coste final (4.39).
-158-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
Figura 4.33 Representación del Coste del potencial vector en una guía rombo, con 40 imágenes. 300
Mhz, fuente situada en el centro de la guía.
En conjunto, obtenemos una distribución del error simétrica, con unos errores
reducidos y un aumento importante de éstos en los picos.
El coste total obtenido ha sido de: 0.3221 , un valor muy reducido para los
niveles de error del potencial vector. Al igual que ocurre con el potencial escalar, una
variación de la distribución de imágenes o puntos tangentes modificará este coste.
Otro ejemplo sobre una guía más compleja puede ser la guía cuadrada pero con
dos lados “redondeados”. Si empleamos un total de 64 imágenes distribuidas siguiendo
el método especular y rectan tangentes para los tramos curvos, el potencial vector en el
contorno de la estructura será:
-159-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
Figura 4.34 Representación del Coste del potencial vector en una guía cuadrada con dos lados
redondeados, con 64 imágenes. 300 Mhz, fuente situada en el centro de la guía.
Una vez definidos el ‘Coste’ del potencial escalar y vector, vamos a utilizarlo
con el fin de optimizar las posiciones de las cargas imágenes al analizar cualquier tipo
de estructura, obteniendo así una representación de los potenciales mucho más precisa.
Los métodos que vamos a presentar en este apartado son una transposición
directa de los métodos estudiados en el capítulo 3, pero generalizados para estructuras
arbitrarias. Conforme vayamos presentado cada método, trataremos de estudiar distintas
estructuras, obteniendo así sus costes y tratándolos de optimizar.
Mostraremos de forma conjunta la optimización de los potenciales escalar y
vector cuando sea posible, mientras que otros tipos de optimizaciones (como las de
descenso de gradiente) únicamente son aplicables al potencial escalar.
-160-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
Figura 4.35 Coste de los potenciales escalar y vector en una guía triangular, 30 imágenes ubicadas
de forma circular sin optimizar su radio. Fuente en el origen de coordenadas.
Figura 4.36 Evolución del coste en función del radio del círculo que ubica a las 30 imágenes. Guía
triangular.
-161-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
Comprobamos como el radio inicial que envolvía la estructura tenía valor muy
cercano al óptimo, por lo que la mejoría de los potenciales no es muy notable: de un
coste del potencial escalar de 0.2176 pasamos a uno de 0.2020, mientras que para el
caso del potencial vector pasamos de un coste de 1.7925 a uno de 1.601. Esto sumado al
coste computacional requerido para una mejoría tan poco significativa, hace
desaconsejar este método en este caso.
No obstante su eficacia variará en función de la estructura que estemos
analizando y de la ubicación de la fuente en el interior de la misma; no obteniendo unos
resultados satisfactorios en la mayoría de las ocasiones.
Los resultados obtenidos con este método son muy satisfactorios (ya que las
imágenes se adaptan a la forma de la geometría a estudiar), especialmente para aquellas
estructuras que no presentan “entrantes”; en el caso de que los tuvieran, el método
comienza a partir de aquel punto en el que todas las imágenes caen fuera de la
estructura.
Como contrapartida, el coste computacional de este método es relativamente
alto, debido a la gran cantidad de operaciones que hay que hacer para evaluar el coste en
el contorno de una estructura arbitraría.
Vamos a mostrar un ejemplo de este método aplicándolo sobre una guía tipo
estrella, empleando para ello 64 imágenes.
Para partir de un resultado general evaluaremos esta guía con esta misma táctica
de ubicación de imágenes sin optimizar; de esta forma podremos apreciar las mejoras
existentes en el caso de que se produzcan. El método general proporciona una ubicación
de imágenes y un coste a la hora de implementar los potenciales de:
-162-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
-163-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
§ ªK º *·
«¦ q k Gv(r , rk )» Gv(r , rp ) ¸¸dc
*
q p C oste ³ ¨
¨ Gv ( r , ro ) Gv ( r , r p ) (4.42)
Contorno © ¬k 1 ¼ ¹
-164-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
Figura 4.40 Situación de partida al analizar una guía rectangular, con 20 imágenes situadas a 0.5Ȝ
de la estructura. Representación del coste de los potenciales escalar en el contorno.
-165-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
Figura 4.41 Evolución del coste del potencial escalar al evaluar una guía rectangular, 20 imágenes.
Donde x (k 1), y (k 1) son dos vectores que representan las posiciones ‘x’ e
‘y’ de cada carga imagen que estamos calculando en la iteración actual,
x (k ), y (k ) representan las posiciones de las imágenes en la iteración anterior, ‘ȝ’ es la
constante de paso y w C (k ) es el gradiente del coste en la iteración k-ésima.
De esta forma, partiendo de unas posiciones iniciales de las imágenes (por
ejemplo cualquiera de las expuestas en este capítulo) alrededor de la estructura arbitraria
-166-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
-167-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
La situación de partida y final de las imágenes tras ejecutar el algoritmo son las
siguientes:
Figura 4.42 Situación de partida y fin de las 20 imágenes empleadas para analizar una guía ridge.
Figura 4.43 Evolución del coste en función del nº de iteración. Optimización de las posiciones de 25
imágenes alrededor de una guía ridge.
-168-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
En primer lugar observamos como para las guías cuadrada y rombo se obtiene
resultados muy parecidos en los dos tipos de optimizaciones empleadas, siendo
ligeramente superior la optimización del radio donde se ubican las imágenes. No
obstante, para el resto de estructuras obtenemos un coste sensiblemente superior al
optimizar la ubicación de las imágenes de forma especular, ya que se adapta a la forma
de la estructura. Es interesante ver el orden de magnitud de los costes en función del
método, y como con una correcta optimización se puede pasar de un coste de más
2492.7 a tan sólo 0.14598,
-169-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
De forma general, podemos analizar cualquier estructura por todos los métodos
presentados, advirtiendo que, en función de las particularidades de cada geometría, unas
técnicas funcionaran de un modo mucho más efectivo que otro.
Para realizar el estudio, nos adecuaremos a los datos que tenemos de los
potenciales, es decir, una guía cuadrada de 66 x66mm . Emplearemos un total de 40
imágenes, distribuidas alrededor del cuadrado siguiendo su estructura. La frecuencia
empleada en este caso la modificamos para adaptarla a los datos de que disponemos,
siendo de 15 Ghz.
-170-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
Figura 4.44 Análisis de una guía cuadrada, 40 imágenes. Fuente en el centro, 15 Ghz.
Figura 4.45 Potenciales normalizados en un corte a 16.5m a lo largo del eje ‘x’ desde el borde
inferior en una guía cuadrada de 66x66mm. Comparativa de método imágenes y método analítico
descrito en [10].
-171-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
Comprobamos como los resultados obtenidos son muy buenos, obteniendo unas
diferencias mínimas entre los dos métodos, a pesar de haber empleado distintas
posiciones de las imágenes para envolver a la estructura.
Indicar que también que se han realizado otras comprobaciones a otras
frecuencias y otros cortes para esta misma guía cuadrada, obteniendo siempre una
convergencia muy buena con los resultados analíticos mostrados en [10].
De esta forma, y a pesar de que como comentamos anteriormente no estamos
representando exactamente el mismo resultado, podemos dar por válidos los resultados
y el cálculo de potenciales realizado hasta ahora así como la formulación presentada
como correcta.
Así, el método presentado es capaz de calcular los potenciales en el interior de
cualquier estructura, con una precisión que también podemos acotar.
En este punto, se abre la puerta a analizar estructuras reales: geometrías en tres
dimensiones donde deberemos de considerar la variación en el eje “z”. De esta forma,
podremos validar los resultados obtenidos en aquellas guías que presentan soluciones
analíticas (por ejemplo la cuadrada) y también en otras más complejas, empleando para
ello como apoyo un software comercial como HFSS©.
En este caso, y a diferencia del capítulo anterior en el que existían una gran
cantidad de programas Fortran 90 realizando cada uno un tipo de optimización, se ha
diseñado un único programa capaz de analizar cualquier estructura definida y de
proporcionar todos los resultados.
Para indicarle al software que estructura queremos analizar así como la posible
optimización de coste que vamos a llevar a cabo sobre ella, emplearemos dos ficheros
de configuración principales:
a) Geom.in: Que ya fue expuesto en el apartado 4.3.2, y en el que se definen todas las
características de la estructura que va a ser analizada: vértices, tramos rectos, tramos
curvos, etc. Además, se añade información acerca de la ubicación de la fuente puntual
en el interior de la guía.
-172-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
5
! Distancia empleada (en función del método)
1
! Optimización gradiente imágenes
0
! Optimización gradiente posiciones
0
Comprobamos como los primeros parámetros son comunes a los que necesitaban
los programas empleados en capítulos anteriores: frecuencia, número de imágenes y
número de puntos de la rejilla.
Posteriormente, aparecen nuevos parámetros como son:
b) Ubicación de imágenes:
Este parámetro hace referencia a la manera en la que ubicaremos las distintas
imágenes alrededor de la estructura. Los posibles valores son:
c) Distancia empleada:
Este parámetro controla la separación de las imágenes a la estructura, en el
caso de que las imágenes sigan a la misma ó la distancia de separación de
especular entre las imágenes y la fuente.
-173-
Capítulo 4: Extensión de la formulación espacial de Green a estructuras arbitrarias
-174-
Capítulo 5: Análisis de cavidades tridimensionales
Una de las ventajas que tiene la técnica desarrollada en los capítulos anteriores
es que puede extenderse fácilmente a cavidades en tres dimensiones, es decir,
considerando el eje ‘z’ que hasta ahora no habíamos tenido en cuenta. Para ello vamos a
situar “tapas” en la parte inferior y superior de la estructura, dejándola totalmente
cerrada.
Además, el interior de la guía a analizar no tiene porqué ser aire (como hemos
considerado hasta ahora al emplear la función de Green en espacio libre), sino que
puede estar compuesto por un medio multicapa, teniendo cada capa un dieléctrico
diferente que caracterizaremos por su permitividad.
De esta forma parámetros como la frecuencia ó dimensión eléctrica van a
adquirir una importancia mucho mayor, pues nos van a definir totalmente la cavidad y
la forma que los potenciales van a adquirir en su interior.
Realizaremos todos los cálculos suponiendo fuentes eléctricas en el interior de
las cavidades; para el caso de fuentes magnéticas el desarrollo es similar (aunque no
dual, pues cambia el tipo de fuente pero no el tipo de pared).
-175-
Capítulo 5: Análisis de cavidades tridimensionales
-176-
Capítulo 5: Análisis de cavidades tridimensionales
Donde:
~
G (k U ) o Función de Green en el dominio espectral
G ( U ) o Función de Green en el dominio espacial.
J n (k U U ) o Función de Bessel de primera especie y orden ‘n’
Uo Distancia entre punto fuente y observación.
Para un medio estratificado, los potenciales base son derivados del dominio
espectral, empleando para ello las corrientes y voltajes evaluadas en un línea de
transmisión transversal equivalente, que describirá al medio estratificado. Un estudio
detallado de estos conceptos se puede encontrar en [14] y [2].
~ ~ VJTE ( z )
G Axx G Ayy (5.3)
jZ
-177-
Capítulo 5: Análisis de cavidades tridimensionales
-178-
Capítulo 5: Análisis de cavidades tridimensionales
También indicar que el tipo de fuente puede ser eléctrica (1) ó magnética (2).
En cuanto al resto de parámetros, se refieren a la precisión de las integrales de la
transformada Sommerfeld.
Comentar finalmente las grandes posibilidades que nos abre el empleo de esta
técnica: podemos evaluar cavidades compuestas por un gran número de “capas” de
distintos dieléctricos, obteniendo así sus potenciales. Además, podemos considerar las
cavidades con o sin “tapas”. Ello nos permitirá analizar estructuras reales, como filtros,
antenas, etc. obteniendo resultados prácticos válidos, como puede comprobarse en [6].
Para el caso del cálculo del potencial escalar, encontrábamos el valor complejo
de las cargas imágenes solventando el siguiente sistema de ecuaciones:
-179-
Capítulo 5: Análisis de cavidades tridimensionales
N
& & & &
¦ q G r , r ´
k 1
k v i k Gv ri , ro ´
(5.4)
i 1,2,3...N
i 1,2,3...N
-180-
Capítulo 5: Análisis de cavidades tridimensionales
& &
wG A ri , rk ´ ~ ~ & &
B oB jk x G A ri , rk ´
wx
& & (5.10)
wG A ri , rk ´ ~ ~ & &
C oC jk y G A ri , rk ´
wy
~
Donde ‘B’ y ‘C’ son las derivadas en el dominio espacial, mientras que ‘ B ’ y
~
‘ C ’ se encuentran en el dominio espectral.
-181-
Capítulo 5: Análisis de cavidades tridimensionales
Figura 5.2 Cilindro a analizar. Radio 1 metro, altura 0.5 metros. 40 imágenes situadas de forma
uniforme alrededor de la estructura. Fuente en el centro, a 0.25 metros de altura.
-182-
Capítulo 5: Análisis de cavidades tridimensionales
Figura 5.3 Representación de potenciales en el interior de una guía cilíndrica, 1 metro de radio, 0.5
de alto, 40 imágenes alrededor. Fuente en el centro, a 0.25metros. 300Mhz.
-183-
Capítulo 5: Análisis de cavidades tridimensionales
Figura 5.4 Representación de potenciales en el interior de una guía cilíndrica, 1 metro de radio, 0.5
de alto, 40 imágenes alrededor. Fuente en el centro, a 0.25metros. 600Mhz.
Observamos como los resultados son coherentes en todo momento con los
esperados, obteniéndose los potenciales de forma correcta.
Indicar que se han realizado múltiples pruebas, obteniendo una concordancia
total al variar el número de imágenes, con convergencia muy rápida.
Figura 5.5 Corte de los potenciales en el eje X en función de distintas distancias de ubicación de las
50 imágenes empleadas al analizar una guía circular de radio 1metro, 0.5 metros de alto. 300 Mhz.
-184-
Capítulo 5: Análisis de cavidades tridimensionales
N
~ & & ~ & &
¦q k > @
S 0 Gv ri , rk ´ > @
S 0 Gv ri , ro ´
(5.14)
k 1
i 1,2,3...N
-185-
Capítulo 5: Análisis de cavidades tridimensionales
En este punto es muy importante tener claro qué ángulo tratamos en cada caso,
que podemos clarificar en la figura 5.6.
El cálculo de la otra ecuación definida en 5.15 es exactamente igual, sin más que
~ ~
sustituir G Axx por GAyy . De esta forma, ya tenemos todos los elementos para resolver el
-186-
Capítulo 5: Análisis de cavidades tridimensionales
sistema de ecuaciones que nos proporcionará los valores de los dipolos imagen,
pudiendo reconstruir el potencial en el interior de cualquier estructura.
Finalmente, esta reconstrucción de los potenciales será llevada a cabo mediante
la ecuación 5.13.
Figura 5.7 Situación al analizar una guía cuadrada 66x66x30 mm. 50 imágenes, situadas siguiendo
su contorno situadas a 15 mm.
-187-
Capítulo 5: Análisis de cavidades tridimensionales
Figura 5.8 Potenciales en el interior de una guía cuadrada 66x66x30 mm. 50 imágenes, situadas
siguiendo su contorno situadas a 15 mm.
Figura 5.9 Ubicación del corte en el eje X en la guía cuadrada para la visualización de los
potenciales
-188-
Capítulo 5: Análisis de cavidades tridimensionales
Figura 5.10 Comparación de los potenciales en una guía cuadrada de 66x66x30mm en un corte en
el eje X a 16.5mm del eje Y, 50 imágenes, distintas posiciones de las mismas. Fuente en el centro. Se
muestra el resultado analítico esperado según [10]. 10 Ghz.
En los resultados mostrados hasta ahora hemos supuesto siempre que la fuente
estaba situada en el centro de la guía cuadrada, y a una altura media. Vamos a mostrar
un ejemplo en el que variaremos la ubicación de la fuente, situándola en una posición de
(-10.7,0) mm en el plano X-Y resultante de observar desde arriba la guía cuadrada, tal y
como observamos en la figura siguiente:
-189-
Capítulo 5: Análisis de cavidades tridimensionales
Figura 5.11 Situación a analizar: guía cuadrada 66x66x30mm, 40 imágenes ubicadas a 33mm del
contorno. Fuente descentrada, situada en (-10.7,0)mm del plano X-Y.
Figura 5.12 Comparación de los potenciales de una guía cuadrada 66x66x30mm en un corte en el
eje X a 16.5mm del eje Y, 40 imágenes, ubicadas a 33mm del contorno. Fuente en el centro.Se
muestra el resultado analítico esperado según [10]. 4 Ghz.
-190-
Capítulo 5: Análisis de cavidades tridimensionales
Figura 5.13 Ubicación de imágenes 40 imágenes alrededor de una guía cuadrada 66x66x30mm en
función aumentada de su distancia a la fuente (descentrada)
Figura 5.14 Representación de los potenciales de una guía cuadrada 66x66x30 en un corte, al
emplear 40 imágenes situadas en función aumentada de su distancia a la fuente ponderada. 4 Ghz.
-191-
Capítulo 5: Análisis de cavidades tridimensionales
Figura 5.15 Comparación de los potenciales en un corte de una guía cuadrada de 66x66x30mm en
función del número de imágenes empleadas, situadas a 33mm de la guía. 10 Ghz.
Para terminar este apartado indicaremos que las estructuras que hemos analizado
hasta ahora son eléctricamente pequeñas. Hemos empleado frecuencias de cómo mucho
10 Ghz, mientras que la guía cuadrada analizada presentaba una altura total de 30mm.
Así:
c 3 108 Altura 0.03
O 9
0.03m o Re lacion 1 (5.18)
f 10 10 O 0.03
-192-
Capítulo 5: Análisis de cavidades tridimensionales
Ghz, este único “anillo” de imágenes no será suficiente para imponer las condiciones de
contorno a lo largo de toda la estructura. Un ejemplo de este caso lo tenemos al analizar
la guía cuadrada tratada hasta ahora con 40 imágenes y una fuente situada en el centro
de la guía. Emplearemos distintas ubicaciones de las imágenes, en función de su
distancia en milímetros al lado de la guía:
Figura 5.16 Comparación de los potenciales de una guía cuadrada 66x66x30 en un corte, al emplear
un anillo de 40 imágenes variando su distancia a la guía, con el resultado analítico de [10]. 15 Ghz.
La única solución para este caso sería situando más de un anillo a lo largo de la
altura de la estructura, de tal forma que la envolviera completamente y se satisficieran
las condiciones de contorno en toda la cavidad. Así introducimos la formulación
“multianillo” que presentamos en el apartado siguiente.
-193-
Capítulo 5: Análisis de cavidades tridimensionales
Figura 5.17 Situación planteada al analizar con 3 anillos una guía cuadrada
Comprobamos como los anillos se han distribuido de tal forma que engloban a
toda la estructura de la guía cuadrada. La forma de ubicar las imágenes en el plano X-Y
no varía, pudiéndose emplear cualquiera de las técnicas presentadas hasta el momento.
En este punto, con la distribución de las imágenes por anillos, vamos a “rehacer”
la formulación empleada, con el fin de tener en cuenta las interacciones entre todas las
imágenes y puntos de contorno. Debido a las múltiples semejanzas que presentan la
formulación en coordenadas cilíndricas y en cartesianas, realizaremos el estudio de
forma genérica, señalando las diferencias en cada caso cuando sea necesario.
-194-
Capítulo 5: Análisis de cavidades tridimensionales
Así, la forma a la que nos referiremos a cada imagen y a cada punto tangente
será como la que se indica en la figura 5.18:
Figura 5.18 Nomenclatura usada para la ubicación de imágenes y puntos tangentes en análisis
multianillo
-195-
Capítulo 5: Análisis de cavidades tridimensionales
i 1,2,3...N ; j 1,2...M
La formulación de las constantes en ‘y’ es la misma, sin más que sustituir las
funciones de Green orientadas según el eje ‘x’ anteriores por otras orientadas según el
eje ‘y’. Al igual que en el caso anterior, hemos considerado que el dipolo fuente se
encuentra orientado a lo largo del eje ‘x’; si estuviera en el eje ‘y’ la formulación es
totalmente dual.
De esta forma, obtenemos un sistema de ecuaciones para las dos formulaciones,
de dimensión 2 M N u 2 M N , cuya solución nos proporcionará los valores
complejos de los dipolos imagen situados en cada anillo ( 2 N en cada anillo, ‘N’
dipolos orientados a lo largo del eje ‘x’ y ‘N’ dipolos orientados a lo largo del eje ‘y’).
Una vez resuelto el sistema de ecuaciones, con un coste computacional cada vez
mayor debido al orden del mismo, podremos recomponer los potenciales en el interior
de la estructura que estemos analizando mediante:
M N
xx & ~ & & ~ & &
G Acy > @
(r ) S0 G Axx r , r0´ ¦¦ I px ,k S 0 G Axx r , rp ,k ´ > @
p 1 k 1
M N
(5.25)
& ~ & &
G yx
Acy (r ) k
y
>
¦¦ I S0 GAyx r , rp,k ´ @
p 1 k 1
-196-
Capítulo 5: Análisis de cavidades tridimensionales
En la figura 5.19 observamos como obtenemos una convergencia total con los
resultados analíticos proporcionados en [10]; igualmente apreciamos como todos los
potenciales convergen entre si, incluyendo las situaciones en las que las imágenes están
muy alejadas de la cavidad (0.35 metros frente a una cavidad de 66x66 milímetros).
Cabe destacar que el coste computacional en este caso es superior al del resto;
ello es debido a la necesidad de obtener las transformadas de Sommerfeld de una
cavidad en la que se han definido varias capas en su interior y al aumento del orden del
sistema de ecuaciones creado.
Figura 5.19 Comparación de los potenciales de una guía cuadrada 66x66x30 en un corte, al emplear
tres anillos de 40 imágenes cada uno, variando su distancia a la guía , con el resultado analítico de
[10]. 15 Ghz.
-197-
Capítulo 5: Análisis de cavidades tridimensionales
Figura 5.20 Comparación de los potenciales de una guía cuadrada 66x66x30 en un corte distinto a
la altura de la fuente, al emplear tres anillos de 40 imágenes cada uno situados a 33mm de la
cavidad , con el resultado analítico de [10]. 15 Ghz.
Vemos como los resultados siguen coincidiendo de forma total con los analíticos
descritos en [10], lo que nos proporciona una validación del método programado.
-198-
Capítulo 5: Análisis de cavidades tridimensionales
Figura 5.21 Potenciales a media altura de una guía cuadrada de 66x66x30mm, 3 anillos de 50
imágenes cada uno, ubicados a 33mm de la cavidad. Frecuencia de 15Ghz.
-199-
Capítulo 5: Análisis de cavidades tridimensionales
5.4.1 Guía dual formada por guía cuadrada con corte en una esquina
-200-
Capítulo 5: Análisis de cavidades tridimensionales
Figura 5.22 Comparación de los potenciales en un corte del eje X, en una guía cuadrada de
66x66x30 con un corte cuadrado en su esquina superior derecha de 18x18mm. 80 imágenes. 10Ghz.
Se realizaron una gran cantidad de pruebas con el fin de corregir este problema:
distribución de las imágenes de forma circular, distribución pegada a la cavidad o muy
separada, e incluso un aumento puntual de las imágenes en la zona con entrante con el
fin de evitar cualquier problema producido en esa zona. Ninguno de los métodos
anteriormente propuestos ha dado el resultado esperado.
Una de las posibles causas de esta falta de convergencia puede ser debida a que
la altura eléctrica que estamos dominando con un único anillo es muy grande; para tratar
de obtener una convergencia, disminuiremos la frecuencia de 10 Ghz a 1Ghz y
realizaremos un estudio similar al anterior (variando la distancia de ubicación de
imágenes) empleando 1,3 y 5 anillos de 80 imágenes cada uno. De esta forma tratamos
de averiguar a que distancia de la cavidad deberemos de ubicar las imágenes para
obtener un resultado óptimo.
-201-
Capítulo 5: Análisis de cavidades tridimensionales
Figura 5.23 Ubicación de 3 y 5 anillos de 80 imágenes alrededor de una guía cuadrada de 66x66x30
mm que presenta un corte de 18x18mm en su esquina superior derecha.
En la siguiente figura se muestra los resultados de los cortes del potencial escalar
en función de la ubicación de las imágenes, para 1,3 y 5 anillos:
1 anillo 3 anillos
5 anillos
Figura 5.24 Potencial escalar en un corte para 1,3 y 5 anillos de 80 imágenes, en función de la
distancia de los mismos a la cavidad cuadrada de 66x66x30mm con un corte de 18x18mm.
-202-
Capítulo 5: Análisis de cavidades tridimensionales
-203-
Capítulo 5: Análisis de cavidades tridimensionales
cavidad. Así esperamos encontrar una zona convergente de ubicación de las imágenes,
con el fin de poder analizar los potenciales en el interior de la guía.
Para el potencial escalar, empleando 1,3 y 5 anillos, obtenemos unos resultados
de:
1 anillo 3 anillos
5 anillos
Figura 5.26 Potencial escalar en un corte para 1,3 y 5 anillos de 80 imágenes, en función de la
distancia cercana de los mismos a la cavidad cuadrada de 66x66x30mm con un corte de 18x18mm.
Realizamos un estudio similar para el caso del potencial vector, con el fin de
obtener la distancia óptima de ubicación de imágenes:
-204-
Capítulo 5: Análisis de cavidades tridimensionales
-205-
Capítulo 5: Análisis de cavidades tridimensionales
1 anillo 3 anillos
5 anillos
Figura 5.28 Potencial escalar en un corte para 1,3 y 5 anillos en función del número de imágenes
por anillo en una cavidad cuadrada de 66x66x30mm con un corte de 18x18mm.
En cuanto al potencial vector, obtenemos unos resultados de:
-206-
Capítulo 5: Análisis de cavidades tridimensionales
Figura 5.30 Visualización de “picos” de error en el potencial escalar de una guía cuadrada de
66x66x30mm, con un corte en una esquina de 18x18mm. 6.122Ghz. 60 imágenes.
-207-
Capítulo 5: Análisis de cavidades tridimensionales
Figura 5.31 Visualización de las matrices del potencial escalar para las guías cuadrada y dual, al
emplear 80 imágenes.
-208-
Capítulo 5: Análisis de cavidades tridimensionales
Una de las formas de mostrar que los resultados son correctos es mostrando su
convergencia al variar la distancia de separación de las imágenes dentro de un
determinado margen, así como mostrando que a partir de un determinado número de
imágenes son resultados son totalmente convergentes.
Así obtenemos:
Figura 5.33 Potenciales en el corte X central para una guía rombo con ángulos internos de 90º
-209-
Capítulo 5: Análisis de cavidades tridimensionales
Figura 5.34 Potenciales en el corte X central para una guía rombo con ángulos internos de 60º
Figura 5.35 Potenciales en el corte X central para una guía rombo con ángulos internos de 45º
-210-
Capítulo 5: Análisis de cavidades tridimensionales
Figura 5.36 Potenciales en el corte X central para una guía rombo con ángulos internos de 30º
Observamos que a pesar de que los resultados del potencial escalar son válidos,
no lo son para el potencial vector, apareciendo diferencias más que significativas entre
las distintas ubicaciones de las imágenes.
De la misma forma, si elegimos una posición de las imágenes (por ejemplo a 10
milímetros de la cavidad) observamos como al variar el número de las mismas no
obtenemos resultados convergentes:
-211-
Capítulo 5: Análisis de cavidades tridimensionales
Figura 5.37 Potenciales en el eje X de una guía rombo con ángulos internos de 30º. Variación del
número de imágenes de un único anillo. Frecuencia de 10Ghz.
Figura 5.38 “Picos” de error al analizar una guía rombo con 2 ángulos internos con 30º a 12 Ghz.
Potencial vector GAyy.
-212-
Capítulo 5: Análisis de cavidades tridimensionales
5.4.3 Conclusiones
-213-
Capítulo 5: Análisis de cavidades tridimensionales
Un ejemplo de esta técnica aplicada en una guía circular con dos dieléctricos
puede ser encontrado en [6].
Aplicaremos el cálculo de las frecuencias de resonancia a una guía cuadrada con
dieléctrico el aire H r 1 , ya que podremos obtener los valores teóricos de las mismas
y nos servirá como una nueva validación de la formulación de las imágenes espaciales
en coordenadas cartesianas presentada.
Se realizará el análisis sobre una guía cuadrada de dimensiones 66 x66 x30mm ,
empleando para ello un total de 60 imágenes, situadas alrededor de la cavidad a una
distancia de 33mm de la misma. La fuente será fijada en fuente en la posición
(15,0,15)mm y el punto fijo de observación en (-15,0,15)mm. La situación que estamos
planteando es la siguiente:
Figura 5.39 Situación de la fuente y un punto de observación en el interior de una guía cuadrada de
66x66x30mm, empleo de 60 imágenes situadas a 33mm.
Modo de Frecuencia de
propagación corte (Ghz)
101 5.4922
111 5.9439
201 6.757
211 7.129
Tabla 5.1 Frecuencias de resonancia de los primeros modos de una cavidad cuadrada de
66x66x30mm
-214-
Capítulo 5: Análisis de cavidades tridimensionales
-215-
Capítulo 5: Análisis de cavidades tridimensionales
-216-
Capítulo 5: Análisis de cavidades tridimensionales
-217-
Capítulo 5: Análisis de cavidades tridimensionales
-218-
Capítulo 5: Análisis de cavidades tridimensionales
1) Geom.in:
2) Datos.in
-219-
Capítulo 5: Análisis de cavidades tridimensionales
3) Separa_anillo.in
4) Anillos_manual.in
-220-
Capítulo 5: Análisis de cavidades tridimensionales
Finalmente destacar que la opción automática distribuirá los anillos y creará las
capas, pero impondrá las mismas características a todas ellas: permitividad y
permeabilidad relativa de valor 1. En la opción manual, se da posibilidad de cambiar los
dieléctricos de la cavidad así como de distribuir los anillos en las posiciones que se
crean más adecuadas.
-221-
Capítulo 5: Análisis de cavidades tridimensionales
-222-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
Indicar también que hemos obviado la inclusión de la guía cuadrada simple, pues
ya fue ampliamente estudiada en el capítulo anterior (hallando sus frecuencias de
resonancia, visualizando sus potenciales, comprobando su convergencia con el número
de imágenes y ubicación de las mismas, etc).
Para finalizar el capítulo, se muestran el estudio de una guía cuadrada con dos
dieléctricos en su interior, con el fin de probar la facilidad que presenta la formulación
presentada para trabajar con medios multicapa.
-223-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
Figura 6.1 Situación de la fuente y un punto de observación en el interior de una guía rectangular
de 132x66x30mm, empleo de 60 imágenes situadas a 12mm.
Al igual que ocurría con la guía cuadrada, la guía rectangular permite calcular de
forma analítica las frecuencias de resonancia (o de corte de los modos) mediante la
ecuación:
2 2 2
c §m· § n · § p·
fc ¨ ¸ ¨ ¸ ¨ ¸ (6.1)
2 © A¹ © B¹ ©C ¹
Donde ‘c’ es la velocidad de la luz en el interior de la guía; ‘m’,’n’ y ‘p’ son el
número de modo excitado (según los ejes ‘x’,’y’,’z’ respectivamente) y ‘A’, ‘B’ y ‘C’
son las dimensiones de la guía a analizar.
Aplicando la ecuación 6.1 podemos encontrar las frecuencias de corte de los
distintos modos que se propagan en la guía rectangular estudiada:
Modo de Frecuencia de
propagación corte (Ghz)
011 5.4922
201 5.4922
111 5.6086
211 5.944
301 6.0516
311 6.4643
021 6.7573
121 6.8522
221 7.1293
Tabla 6.1 Frecuencias de resonancia de los primeros modos de una cavidad rectangular de
132x66x30mm
-224-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
-225-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
Tras obtener las frecuencias de corte y validarlas con las obtenidas de forma
teórica, vamos a mostrar la distribución de los potenciales en el interior de la guía
rectangular, visualizando de forma clara que modo se excita en cada frecuencia de corte.
Tal y como comentamos anteriormente, validaremos la distribución del potencial
escalar mediante la “forma” adquirida por el campo eléctrico E z (en lo modos TM), ya
que la definición de ambas magnitudes está muy relacionada (como se puede observar
en 5.14 para el potencial escalar y en [8] para el campo eléctrico). La representación de
este campo eléctrico la realizaremos mediante HFSS©, al igual que en ejemplos
anteriores.
-226-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
-227-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
Por último y para terminar con el análisis de esta guía, vamos a mostrar los
resultados en una tercera resonancia. Para ello excitaremos el modo 311 a una
frecuencia de 6.458 Ghz:
-228-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
-229-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
De esta forma podemos dar por válidos los resultados, pues el método ha sido
capaz de obtener con una gran precisión las frecuencias de resonancia de la cavidad y ha
obtenido unos potenciales en el interior de la guía cuya distribución hemos validado con
HFSS©.
En este apartado vamos a analizar una guía triangular de altura 30mm y con una
geometría en 2D mostrada en la figura 6.6, empleando para ello un número muy
reducido de imágenes, tan sólo 30.
Figura 6.6 Situación de la fuente y un punto de observación en el interior de una guía triangular de
altura 30mm, empleo de 30 imágenes situadas a 12mm.
-230-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
Tras obtener las resonancia, vamos a realizar una comparación entre las
frecuencias obtenidas y las proporcionadas por HFSS© con el fin de ver la precisión
que conseguimos con el método de las imágenes espaciales.
-231-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
-232-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
-233-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
Figura 6.9 Potenciales en la primera resonancia de una cavidad triangular con 30mm de altura.
Comparación con el campo eléctrico proporcionado por HFSS©.
-234-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
Figura 6.10 Potenciales en la segunda resonancia de una cavidad triangular con 30mm de altura.
Comparación con el campo eléctrico proporcionado por HFSS©.
De esta forma podemos dar por válidos los resultados, pues el método ha sido
capaz de obtener con una gran precisión las frecuencias de resonancia de la cavidad y ha
obtenido unos potenciales en el interior de la guía cuya distribución hemos validado con
HFSS©.
-235-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
segundo lugar para que fueran equivalente deberíamos rotar también la ubicación de la
carga/dipolo fuente situado en su interior (y no situarlos en el eje X e Y).
De esta manera aparecerán diferencias respecto al análisis de la guía cuadrada
tanto en las frecuencias de corte encontradas como en la visualización de los potenciales
obtenidos.
Figura 6.11 Situación de la fuente y un punto de observación en el interior de una guía rombo de
dimensiones 66x66x30mm y sus 4 ángulos internos a 90º. Empleo de 60 imágenes situadas a 12mm.
-236-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
Tras obtener las resonancia, vamos a realizar una comparación entre algunas de
las frecuencias obtenidas y las proporcionadas por HFSS© con el fin de ver la precisión
que conseguimos con el método de las imágenes espaciales.
-237-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
-238-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
En cuanto a los potenciales cruzados, obtenemos unos valores idénticos para las
dos componentes (GAyx y GAxy):
Figura 6.14 Potenciales cruzados en la primera resonancia de una guía rombo de dimensiones
66x66x30mm.
Así la situación a analizar ubicando los 3 anillos con 60 imágenes cada uno será
la siguiente:
Figura 6.15 Situación de 3 anillos con 60 imágenes cada uno al analizar una guía rombo de
dimensiones 66x66x30mm con sus 4 ángulos internos a 90º.
-239-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
-240-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
De esta forma podemos dar por válidos los resultados, pues el método ha sido
capaz de obtener con una gran precisión las frecuencias de resonancia de la cavidad y ha
obtenido unos potenciales en el interior de la guía cuya distribución hemos validado con
HFSS©.
En este apartado vamos a analizar una guía con forma trapecio, con una altura de
30mm. El lado mayor tendrá una longitud total de 66mm, 40mm tendrá el lado menor y
una separación entre ambos de 66mm. Analizaremos la guía con un total de 60
imágenes distribuidas alrededor de la misma a una distancia de 12mm.
La situación que estamos analizando es la siguiente:
Figura 6.17 Situación de 1 anillo con 60 imágenes situados sobre una guía en forma de trapecio, de
altura 30mm.
Esta guía no tiene solución analítica, por lo que emplearemos otra vez HFSS©
para validar las frecuencias de corte obtenidas así como la distribución de los
potenciales calculados.
-241-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
Tras obtener las resonancia, vamos a realizar una comparación entre algunas de
las frecuencias obtenidas y las proporcionadas por HFSS© con el fin de ver la precisión
que conseguimos con el método de las imágenes espaciales.
-242-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
-243-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
-244-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
-245-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
Así, podremos concluir que los potenciales obtenidos por el método de las
imágenes espaciales para esta guía son válidos, pues hemos realizado una validación
tanto en la ubicación de las frecuencias de resonancia (con una gran precisión) como de
la distribución del potencial escalar en el interior de la guía.
-246-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
Figura 6.22 Situación de 1 anillo con 90 imágenes situados sobre una guía cuadrada de
66x66x30mm con un corte de 18x18mm en su esquina superior derecha.
-247-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
Comprobamos como aparecen diferencias mayores que las del resto de guías, no
obstante siguen siendo pequeñas y podemos dar por buenos los resultados.
-248-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
60 imágenes 60 imágenes
80 imágenes 80 imágenes
90 imágenes 90 imágenes
-249-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
Figura 6.25 Campo eléctrico Ez obtenido en una guía dual en la primera resonancia, a 6.122Ghz.
-250-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
-251-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
-252-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
errores numéricos en las esquinas adyacentes al error, que a pesar de encontrarse muy
localizados, si que tienen su cierta influencia en el resto de la guía.
Para concluir los resultados de esta guía, indicar en primer lugar la falta de
convergencia al variar el número de imágenes y las posiciones de las mismas. Ello nos
lleva a que no sepamos que situación es la más propicia para analizar la guía.
Tras múltiples pruebas, observamos como una ubicación de 90 imágenes a 5mm
de la guía proporcionaba un menor error; la visualización de este error es sencilla pues
se manifiesta como “picos” de potencial en las esquinas adyacentes al entrante. Con esta
configuración, conseguimos obtener unas frecuencias de corte muy similares a las
proporcionadas por HFSS© y una distribución de potencial que, salvo el error muy
localizado en las esquinas adyacentes al entrante, se corresponde con el obtenido con
HFSS©.
Indicar que puede existir otra configuración de ubicación de imágenes y número
de las mismas que proporcione un error en la guía mucho menor; una posible forma de
encontrar esta configuración sería mediante el concepto de “Coste” introducido en el
capítulo 3, extendiendo el mismo a estructuras tridimensionales.
Otro ejemplo de guía de ondas cuyo análisis de potenciales puede ser importante
es la denominada guía “flange”. La guía que analizaremos parte de una estructura
cuadrada de dimensiones 66x66x30mm a la que se realiza un corte en diagonal en su
vértice superior derecho. Realizaremos el análisis empleando un total de 50 imágenes
situadas alrededor de la cavidad a una distancia de 5 mm. La situación que estamos
analizando es la siguiente:
Figura 6.29 Situación de 1 anillo con 90 imágenes situados sobre una guía flange de altura 30mm.
Esta guía, al igual que ocurría con la guía dual analizada anteriormente, presenta
un entrante que puede ocasionar “picos” de error, por lo que realizaremos un estudio de
su convergencia en función tanto de las imágenes empleadas como de la posición de las
mismas.
En primer lugar, mostraremos cortes para los potenciales en el eje ‘X’ para un
total de 50 imágenes en función de la distancia de las mismas a la guía. Emplearemos
para ello una frecuencia de 7Ghz. Los resultados que obtenemos son los siguientes:
-253-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
-254-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
Podemos observar como al aumentar el número de anillos los resultados son más
convergentes: en el potencial escalar las diferencias son inexistentes y en el potencial
vector, aunque todavía existen, se han reducido de manera considerable.
Mostraremos finalmente los resultados de convergencia al emplear 5 anillos
alrededor de la cavidad, comprobando como cada vez los resultados son más
convergentes e independientes de la distancia de las imágenes a la guía:
-255-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
De esta forma, vamos a situar las imágenes a 5mm de la cavidad. Sobre esta
distancia, vamos a variar el número de imágenes empleadas, con el fin de estudiar su
convergencia y encontrar un número mínimo de las mismas que nos proporcione un
resultado válido. Emplearemos un total de 5 anillos para esta última comparación, pues
era en la situación en la que la convergencia era máxima. Así, obtenemos unos
resultados de:
-256-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
Tras obtener las resonancia, vamos a realizar una pequeña comparación entre
algunas de las frecuencias obtenidas y las proporcionadas por HFSS© con el fin de ver
la precisión que conseguimos con el método de las imágenes espaciales.
De la tabla anterior podemos destacar como las dos frecuencias introducidas por
el potencial escalar (6.023 y 7.155) obtiene un error mucho más reducido que el resto.
Ello es debido a que en el barrido frecuencia se ha empleado un único anillo, debido al
gran coste computacional que requiere emplear un sistema formado por 3 ó 5 anillos.
No obstante, los resultados siguen siendo bastante satisfactorios, con un error relativo
inferior al 0.82% en todos los casos.
-257-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
-258-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
Figura 6.35 Potenciales a la frecuencia de resonancia de 7.155Ghz de una cavidad flange de altura
30mm. Comparación con el campo eléctrico proporcionado por HFSS©.
-259-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
Figura 6.36 Potenciales a la frecuencia de resonancia de 8.419Ghz de una cavidad flange de altura
30mm. Comparación con el campo eléctrico proporcionado por HFSS©.
Hasta este punto hemos obtenido unos resultados válidos en todas las guías que
hemos analizado, a pesar de tener más o menos problemas para encontrar una ubicación
y número de imágenes capaces de analizar las estructuras propuestas.
-260-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
En este punto vamos a analizar una guía “estrella”, que parte de una guía
cuadrada de 132 x132 x30mm , sobre la que se han realizado cuatro cortes cuadrados en
las esquinas, siendo cada corte de una dimensiones de 44 x44mm .
Observamos como las dimensiones de la guía propuesta en este caso son
diferentes a las analizadas hasta ahora (casi el doble de longitud), es por ello que el
número de imágenes que tenemos que emplear para su análisis también debe de
aumentar. En primer lugar, vamos a mostrar la guía a analizar empleando para ello un
total de 150 imágenes, situadas a una distancia de 10mm de la guía:
Figura 6.37 Situación de 1 anillo con 150 imágenes situados sobre una guía estrella de altura 30mm.
Figura 6.38 “Picos” de error en los vértices de una guía estrella al ser analizada con 70 imágenes
-261-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
Vemos como aparece el mismo efecto que en el caso de la guía dual, debido a la
presencia de no 1 sino 4 entrantes. Esto hace que el análisis de esta guía sea muy
complejo mediante la técnica de las imágenes espaciales.
-262-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
Así por ejemplo, el primer modo se propaga a la frecuencia de 5.7234 Ghz sin
que ésta frecuencia sea detectada. Ello podría ser debido en una primera instancia a un
-263-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
problema del método que calcula las frecuencias de resonancia, sin embargo, tras
realizar múltiples simulaciones no se podido encontrar una ubicación de imágenes que
proporcione la distribución adecuada del potencial escalar para este modo. Podemos
observar en la figura 6.41 como con HFSS© aparece esta primera resonancia de forma
clara, mientras que con el método de las imágenes no.
Hasta este punto todos los resultados son incoherentes, sin embargo aparecen
muestra de un aparecente buen comportamiento del método a otra frecuencia, como es
la 6.9358, donde se consigue una distribución adecuada de los potenciales:
-264-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
Para terminar este apartado, vamos a realizar un análisis de la guía rige. La guía
esta basada en una cuadrada de dimensiones 66 x66 x30mm sobre la que se realizan dos
cortes cuadrados de dimensiones 22 x22mm en los lados superior e inferior.
-265-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
Figura 6.43 Situación de 1 anillo con 90 imágenes situados sobre una guía ridge de altura 30mm.
Debido a los dos entrantes que posee esta guía, nos encontramos ante los
mismos problemas que en la guía estrella: una falta de convergencia entre las distintas
ubicaciones de las imágenes y también en cuanto a su número. Ello nos impide a
primera instancia el cálculo adecuado de los potenciales en el interior de la estructura,
ya que aparecen “picos” de error en las esquinas adyacentes de los entrantes,
invalidando las posibles convergencias de resultados.
Figura 6.44 “Picos” de error en los vértices de una guía ridge al ser analizada con 70 imágenes
En primer lugar, vamos a ubicar las imágenes a 5mm de la cavidad (pues parece
que es una distancia óptima para que las imágenes no se solapen dentro de los entrantes)
y variaremos el número de imágenes empleadas, con el fin de estudiar la convergencia.
Así, a 6Ghz de frecuencia, obtenemos los siguientes cortes de potenciales:
-266-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
-267-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
Al igual que ocurría con la guía estrella aparecen muy pocas resonancias; eso es
debido a que la fuente es incapaz de excitar al modo en estas geometrías complejas.
Compararemos los resultados con los obtenidos con HFSS para evaluar su precisión:
Las frecuencias se han calculado con una precisión bastante pobre, pero lo más
importante es que faltan muchas resonancias por obtener pues los modos no han sido
excitados. A pesar de que puede ocurrir que se sitúe la fuente en un nulo particular de
un modo, no es lógico que únicamente un modo en un barrido tan amplio. No obstante,
también puede ser debido a que la fuente ha sido incapaz de excitar al resto de modos;
por ello se ha realizado el análisis con distintas posiciones de la fuente, obteniendo
siempre el mismo resultado.
-268-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
-269-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
Tras estos resultados podemos concluir que el método no funciona del todo bien
para esta guía, ya que propicia unos “picos” de error en las esquinas adyacentes a los
entrantes que invalidan el resto de resultados. A pesar de ello, se siguen consiguiendo
ciertos resultados puntuales válidos. Los potenciales calculados serán totalmente válidos
cuando se consigan eliminar estos errores, existiendo entonces total convergencia entre
las distintas representaciones.
Hasta ahora, todos los análisis de guías han sido realizados con aire ( H r 1) en
su interior. No obstante, la teoría de las imágenes espaciales puede aplicarse de forma
directa a medios multicapa gracias a la versatilidad de la transformada de Sommerfeld
[9]. Sin más que definir las características de cada capa (como permitividad eléctrica y
permeabilidad magnética, tal y como se expuso en el capítulo anterior), podremos
realizar el cálculo de los potenciales en el interior de guías estratificadas. Ello
proporciona una herramienta muy poderosa de análisis para estructuras reales mediante
la ecuación integral.
-270-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
En primer lugar vamos a realizar un estudio con 50 imágenes (pues para el caso
de espacio libre eran más que suficientes) en función de la distancia de las mismas a la
guía cuadrada multicapa. Emplearemos para ello una frecuencia de 5Ghz. Los
resultados son:
-271-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
Figura 6.51 Potenciales vector al realizar un barrido en un punto fijo entre 2 y 7 Ghz en una guía
cuadrada con 2 dieléctricos. Comparación con el método analítico descrito en [10].
-272-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
Para una estructura como la presentada, con una altura eléctrica considerable
(30mm, para una frecuencia de 6Ghz hace un total de 0.6Ȝ) serían necesarios una gran
cantidad de anillos, siendo el coste computacional muy elevado.
La convergencia aumenta, no obstante, cuando esta altura eléctrica es menor. Si
analizamos una cavidad cuadrada con unas dimensiones diferentes, de 120 x120mm y
una altura total de 24mm (12 mm cada capa) con 8 anillos que rodeen la estructura,
obtendremos una convergencia adecuada. Por ejemplo, mostramos las medidas a
5.5Ghz en la figura 6.53.
Figura 6.53 Potencial escalar a 5.5Ghz en una guía cuadrada de 120x120mm, con dos capas de
12mm, una de ellas aire y la otra con permitividad de valor 2.
-273-
Capítulo 6: Ejemplos de Análisis de guías-onda tipo
Una posible línea de continuación del proyecto sería un estudio entre la cantidad
de anillos que envuelven la estructura (y el número de imágenes que emplea cada uno
de ellos) en función de la altura eléctrica de la cavidad a analizar.
Así, con un único programa, se han logrado estudiar todas las guías presentadas
en este capítulo, así como es posible realizar análisis de otras muchas, sin más que
definirlas de manera adecuada. Es destacable la gran versatilidad del software, siendo
capaz de adaptarse a los distintos casos y proporcionar sus resultados de forma
automática, tanto en el cálculo de los mismo (en Fortran 90) como en su visualización
(Matlab ©).
-274-
Capítulo 7: Conclusiones y Líneas futuras
-275-
Capítulo 7: Conclusiones y Líneas futuras
Finalmente, se mostraron los resultados del cálculo de los potenciales para una
gran cantidad de guías con distintas geometrías: rectangular, rombo, triangular, flange,
etc. Se consiguió validar el método empleando el software comercial HFSS©,
obteniéndose una convergencia casi total de las frecuencias de corte obtenidas y una
visualización del potencial escalar, cuya distribución debe de coincidir con el campo
eléctrico Ez.
En este punto cabe destacar la robustez del software: definiendo una estructura
arbitraría en un fichero de texto, es capaz de realizar de forma automática todo el
estudio de los potenciales en su interior, teniendo en cuenta su altura (estratificada si es
necesario), vértices etc.
Debido a la amplitud del proyecto, han surgido durante la elaboración del mismo
una gran cantidad de líneas de trabajo con las que se puede seguir el mismo.
-276-
Capítulo 7: Conclusiones y Líneas futuras
Por último, comentar que todos los desarrollos presentados a lo largo del
proyecto se han basado en la suposición de fuentes eléctricas; una continuación directa
del proyecto será realizar la formulación basándose en fuentes magnéticas (que no es
dual, pues aunque cambia el tipo de fuente no cambia el tipo de contorno, por lo que la
condición a cumplir varía). No obstante, su obtención debe de ser inmediata, dada su
gran similitud con los desarrollos mostrados.
-277-
Capítulo 7: Conclusiones y Líneas futuras
-278-
ANEXO I: Resultados númericos de la integración por pulsos
ANEXO I
x a) Posición (0Ȝ,0Ȝ):
Nº
Tipo Cargas 5 10 15 20 25
Técnica
7.1487e- 1.3369e-
Point-matching 10.1294 1.3525 0.0061
005 006
1.3369e-
Trapecio 10.1294 74.5772 0.0061 0.1318
006
9.1169e-
Simpson 30.3169 47.4165 0.0041 3.5133
007
8.4151e-
Simpson 3/8 35.1437 31.8513 0.0039 1.1527
007
8.7649e-
Bode 25.4819 28.0553 0.0039 0.5095
007
Semi-abierta 4 4.5510e- 8.7649e-
25.9997 0.8972 0.0039
puntos internos 005 007
Semi-abierta 10 4.5522e- 8.7649e-
26.0353 0.8972 0.0039
puntos internos 005 007
8.6352e-
Trapecio 3.8555e- 4.54981e-
26.6402 0.8973 007
generalizado 003 005
Tabla AI.1 Coste de distintos métodos de integración por pulso, en función del número de
imágenes, posición fuente en (0.0Ȝ ,0.0Ȝ)
-279-
ANEXO I: Resultados númericos de la integración por pulsos
x b) Posición (0.2Ȝ,0.2Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Técnica
Point- 5.6664e- 5.9574e- 4.0841e-
12.1475 16.5825 0.0972
matching 004 006 008
5.9574e-
Trapecio 12.1475 18.3701 0.0972 1.5026 0.0112
006
1.9408e+ 2.3079e+ 4.7875e-
Simpson 17.8651 0.0938 416.8706
005 009 006
2.7496e+ 1.8730e+ 4.6727e- 1.7996e+
Simpson 3/8 26.9861 0.0920
005 009 006 003
3.6733e+ 1.1266e+ 4.6372e-
Bode 40.2739 0.0907 671.2963
005 009 006
Semi-abierta
5.3603e- 4.6498e- 7.5398e-
4 puntos 40.6732 5.5586 0.0912
004 006 009
internos
Semi-abierta
4.6314e- 1.2779e-
10 puntos 33.8777 42.1664 0.0908 0.0016
006 007
internos
2.6125e-
Trapecio 9.0850e- 2.1552e- 4.6282e-
33.1510 65.0231 007
generalizado 002 003 006
Tabla AI.2 Coste de distintos métodos de integración por pulso, en función del número de
imágenes, posición fuente en (0.2Ȝ ,0.2Ȝ)
x c) Posición (0.4Ȝ,0Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Técnica
Point- 2.5848e- 4.6492e-
16.1649 19.9669 0.3775 0.0026
matching 005 007
2.5848e-
Trapecio 16.1649 21.7083 0.3775 0.8872 0.0088
005
2.3114e-
Simpson 21.3260 18.0569 0.4148 0.7403 0.0094
005
2.2760e-
Simpson 3/8 26.2434 13.2931 0.4059 0.6847 0.0140
005
2.2615e-
Bode 37.0349 12.3691 0.3993 0.7613 0.0049
005
Semi-abierta
2.2746e- 5.2649e-
4 puntos 37.9058 6.6327 0.4029 0.0024
005 007
internos
Semi-abierta
0.3996 2.2608e- 5.1975e-
10 puntos 33.4921 6.6674 0.0024
005 007
internos
-280-
ANEXO I: Resultados númericos de la integración por pulsos
x c) Posición (0.4Ȝ,0.4Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Técnica
Point- 2.8172e- 2.5517e-
26.6829 23.6945 1.4449 0.0132
matching 004 005
2.8172e-
Trapecio 26.6829 68.8166 1.4449 0.3749 0.0172
004
3.3688e+ 1.7612e+ 3.6452e-
Simpson 17.1538 1.9235 81.7399
003 011 004
4.6252e+ 2.3064e+ 3.5823e-
Simpson 3/8 19.7604 1.8716 69.2776
003 012 004
5.5227e+ 4.4826e+ 3.5422e-
Bode 21.8704 1.8319 63.6620
003 011 004
Semi-abierta
21.1368 1.8271 3.4845e- 7.5339e-
4 puntos 81.5308 0.2039
004 004
internos
Semi-abierta
20.4663 1.8306 1.6533 3.5342e-
10 puntos 685.1304 0.0061
004
internos
3.5392e- 7.335e-
Trapecio 20.3296 1285.018 1.8339 2.2861
004 003
generalizado
Tabla AI.4 Coste de distintos métodos de integración por pulso, en función del número de
imágenes, posición fuente en (0.4Ȝ ,0.4Ȝ)
x d) Posición (0.5Ȝ,0Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Técnica
Point- 9.6659e- 3.7540e-
17.9904 19.9656 0.9056 0.0070
matching 005 006
9.6659e-
Trapecio 17.9904 32.0972 0.9056 0.7542 0.0121
005
1.2030e-
Simpson 23.9292 13.5542 1.1268 0.5670 0.0200
004
Simpson 1.1790e-
25.1236 14.2283 1.0987 0.4035 0.0038
3/8 004
1.1630e-
Bode 27.9842 14.4025 1.0774 0.8057 0.0129
004
-281-
ANEXO I: Resultados númericos de la integración por pulsos
Semi-
12.7546
abierta 4 28.9088 1.0889 1.1726e- 8.4674e-
0.0101
puntos 004 006
internos
Semi-
abierta 10 27.5673 12.7675 1.0782 1.1636e- 8.1904e-
0.0098
puntos 004 006
internos
Trapecio
8.9901e- 1.16444e- 7.61462e-
generali- 27.2404 12.5845 1.07945
003 004 006
zado
Tabla AI.5 Coste de distintos métodos de integración por pulso, en función del número de
imágenes, posición fuente en (0.5Ȝ ,0.0Ȝ)
x e) Posición (0.7Ȝ,0.7Ȝ):
Nº cargas
Tipo 5 10 15 20 25 30
Integración
Point-
22.8267 8.0995 1.6061 4.2177 0.4510 1.1371
matching
Trapecio 22.8267 51.6313 1.6061 5.9493 0.4510 5.1767
1.7441e+
Simpson 21.7667 6.0057 3.0976 0.8455 2.1369
013
3.6043e+
Simpson 3/8 36.9026 6.7298 2.9544 0.8493 1.3641
012
7.7497e+
Bode 26.3056 36.9438 2.4971 0.7794 0.8441
012
Semi-abierta 42.9424
2.1737
4 puntos 6.0783 7.1962 0.7504 0.7286
internos
Semi-abierta 1.8135
356.1569
10 puntos 41.5098 10.0066 0.7050 0.7609
internos
Trapecio 34.4937 10.7793 1.5838 855.8030 0.6978 0.7804
generalizado
Tabla AI.6 Coste de distintos métodos de integración por pulso, en función del número de
imágenes, posición fuente en (0.7Ȝ ,0.7Ȝ)
-282-
BIBLIOGRAFÍA
ANEXO II
x a) Posición (0Ȝ,0Ȝ):
x b) Posición (0.2Ȝ,0.2Ȝ):
-283-
BIBLIOGRAFÍA
14 6.4884 1.2442
15 4.9207 1.1727 0.0972
16 3.7564 0.8740 0.0359
17 2.6878 0.7288 0.0862
18 2.6345 0.6295 0.0917
19 1.0413 0.2856 0.0770 0.0015
20 0.1475 0.2405 0.0831 0.0038 5.6664e-004
21 0.1066 0.3402 0.0733 0.0080
22 0.0402 0.1910 0.0462 0.0081
23 0.0237 0.1302 0.0428 0.0073
24 0.0142 0.0967 0.0300 0.0053
25 0.0119 0.0591 0.0259 0.0056 5.9574e-006
Tabla AII.2 Coste de la ubicación de cargas en picos de potencial, con distintas situaciones iniciales,
en función del número de imágenes y con posición fuente en (0.2Ȝ ,0.2Ȝ)
x c) Posición (0.4Ȝ,0Ȝ):
x c) Posición (0.4Ȝ,0.4Ȝ):
-284-
BIBLIOGRAFÍA
14 33.7669 11.8571
15 18.2913 9.5892 1.4449
16 17.0453 6.4370 0.5339
17 10.0158 5.3626 1.4715
18 6.2894 6.6233 1.2682
19 2.1613 5.9441 1.0167 0.0348
20 0.1773 5.5390 1.1615 0.1500 0.0132
21 0.1042 5.5099 0.8259 0.1389
22 0.1352 6.2922 0.9196 0.1125
23 0.0798 3.7750 0.7934 0.1062
24 0.1049 4.3934 0.7925 0.1053
25 0.0692 3.3206 0.6214 0.0900 2.8172e-004
Tabla AII.4 Coste de la ubicación de cargas en picos de potencial, con distintas situaciones iniciales,
en función del número de imágenes y con posición fuente en (0.4Ȝ ,0.4Ȝ)
x d) Posición (0.5Ȝ,0Ȝ):
x e) Posición (0.7Ȝ,0.7Ȝ):
-285-
BIBLIOGRAFÍA
13 1.1217 11.2623
14 1.1220 2.2962
15 1.1530 2.0504 1.6061
16 1.0949 2.1947 0.8639
17 0.9710 2.0100 0.6412
18 1.0705 1.9359 0.6069
19 1.8399 1.9205 0.6187 3.7204
20 0.9958 1.8298 0.5594 0.1986 4.2177
21 1.7867 1.8400 0.5878 0.2029
22 1.5244 1.7896 0.4508 0.2289
23 1.5881 1.7708 0.4382 0.1698
24 1.9221 1.7239 0.3882 0.1303
25 0.6174 1.3796 0.3731 0.1441 0.4510
Tabla AII.6 Coste de la ubicación de cargas en picos de potencial, con distintas situaciones iniciales,
en función del número de imágenes y con posición fuente en (0.7Ȝ ,0.7Ȝ)
-286-
BIBLIOGRAFÍA
Bibliografía
[1] O. I. Sukharevsky and S. V. Kukobko and A. Z. Sazonov, “Volume Integral
Equation Analysis of a two-dimensional radome with a sharp nose”, IEEE Transactions
on Antennas and Propagation, Vol. 53, No. 4, pp. 1500-1506, April 2005.
[2] K.A. Michalski and J.R. Mosig, “Multilayered media Green’s functions in integral
equation formulations”, IEEE Trans. Antennas Propag., vol 45, no. 3, pp. 508-519, Mar.
1997.
[3] Yi-Chi Shih, “Design of Waveguide E-Plane Filters with all Metal Inserts”, IEEE
Transactions on Microwave Theory and Techniques”, Vol. 32, No. 7, pp. 695-704, July
1984.
[4] Alejandro Alvarez Melcon and Juan R. Mosig, “Two Techniques for the Efficient
Numerical Calculation of the Green's Functions for Planar Shielded Circuits and
Antennas”, IEEE Transactions on Microwave Theory and Techniques, Vol. 48, No. 9,
pp. 1492-1504, September 2000.
[5] Jind-Yeh Lee and Tzyy-Sheng Horng and Nicolaos G. Alexopoulos, “Analysis of
Cavity-Backed Aperture Antennas With a Dielectric Overlay”, IEEE Transactions on
Antennas and Propagation, Vol. 42, No. 11, pp. 1556-1562, November 1994.
[8] Balanis, C.A. (1989). “Advanced Engineering Electromagnetics”. John Wiley and
Sons.
[9] J.R. Mosig. Integral Equation Technique. New York: Wiley, 1989
[10] Alejandro Álvarez Melcon, “Applications of the integral equation technique to the
analysis and synthesis of multilayered printed shielded microwave circuits and cavity
backed antennas” Thèse EPFL, no 1901 (1998), Département d'électricité. Dir.: Juan
Mosig.
[12] William H. Press, Saul A. Teukolsky, William T. Vetterling and Brian P. Flannery,
“Numerical Recipes in Fortran 90”, Second Edition, 1999, Cambridge University Press.
[13] J.G. de Jalón, J.I. Rodríguez, A. Brazales “Aprenda Matlab 6.1 como si estuviera
en primero”
-287-
BIBLIOGRAFÍA
[14] S.-G. Pan and I. Wolff, “Scalarization of dyadic spectral Green’s function and
network formalism for the three-dimensional full-wave analysis of planar lines and
antennas”, IEEE Trans. Microw. Theory Tech, vol. 42, no. 11, pp 2128-2127, Nov.
1994.
[15] A. A. Melcón, J.R. Mosig, and M. Guglielmi, “Efficient CAD of boxed microwave
circuits based on arbitrary rectangular elements”, IEEE Trans. Microw. Theory Tech.,
vol 47, no. 7, pp. 1045-1058. Jul. 1999.
[18] A.I.Nosich, V.B. Yurchenko and A.Altintas “Numerically exact analysis of a two-
dimensional variable-resistivity reflector fed by a complex point source”, IEEE Trans.
Antennas Propag., vol 45, no.11, pp. 1592-1601, November. 1997.
[20] Yi Huang “Conducting triangular chambers for EMC measurements”, Meas. Sci.
Technol. 10 (1999). L21-L24.
[24] Ángel Cardama, Lluís Jofre Roca “Antenas” Edicions UPC. ISBN: 84-8301-625-7
-288-
BIBLIOGRAFÍA
[28] http://es.wikipedia.org/wiki/James_Clerk_Maxwell
[29] www.upv.es/antenas/Tema_2/Ecuaciones_Maxwell.htm
-289-