Articulo de referencia

Método de matriz T

El método de la matriz de transición ( método de la matriz T o TMM ) es una técnica computacional de dispersión de luz por partículas no esféricas formulada originalmente por Pe...

El método de la matriz de transición ( método de la matriz T o TMM ) es una técnica computacional de dispersión de luz por partículas no esféricas formulada originalmente por Peter C. Waterman en 1965. [ 1 ] [ 2 ] La técnica también se conoce como método de campo nulo y método de condición de contorno extendida (EBCM). [ 3 ] En el método, los elementos de la matriz se obtienen al igualar las condiciones de contorno para las soluciones de las ecuaciones de Maxwell . Se ha extendido ampliamente para incorporar diversos tipos de medios lineales que ocupan la región que encierra al dispersor. [ 4 ] El método de la matriz T demuestra ser altamente eficiente y se ha utilizado ampliamente en el cálculo de la dispersión electromagnética de partículas simples y compuestas. [ 5 ]

Definición de la matriz T

El campo eléctrico incidente y dispersado se expande en funciones de onda vectoriales esféricas (SVWF), que también se encuentran en la dispersión de Mie . Son las soluciones fundamentales de la ecuación vectorial de Helmholtz y pueden generarse a partir de las soluciones fundamentales escalares en coordenadas esféricas , las funciones de Bessel esféricas de primer tipo y las funciones de Hankel esféricas. En consecuencia, hay dos conjuntos de soluciones linealmente independientes denotados comoMETRO1,norte1{\displaystyle \mathbf {M} ^{1},\mathbf {N} ^{1}}yMETRO3,norte3{\displaystyle \mathbf {M} ^{3},\mathbf {N} ^{3}}, respectivamente. También se les llama SVWF regulares y salientes, respectivamente. Con esto, podemos escribir el campo de incidentes como

miinortedo=norte=1metro=nortenorte(ametronorteMETROmetronorte1+bmetronortenortemetronorte1).{\displaystyle \mathbf {E} _{inc}=\sum _{n=1}^{\infty }\sum _{m=-n}^{n}\left(a_{mn}\mathbf {M} _{mn}^{1}+b_{mn}\mathbf {N} _{mn}^{1}\right).}

El campo disperso se expande en SVWF radiantes:

misdoat=norte=1metro=nortenorte(FmetronorteMETROmetronorte3+gramometronortenortemetronorte3).{\displaystyle \mathbf {E} _{scat}=\sum _{n=1}^{\infty }\sum _{m=-n}^{n}\left(f_{mn}\mathbf {M} _{mn}^{3}+g_{mn}\mathbf {N} _{mn}^{3}\right).}

La matriz T relaciona los coeficientes de expansión del campo incidente con los del campo dispersado.

(Fmetronortegramometronorte)=T(ametronortebmetronorte){\displaystyle {\begin{pmatrix}f_{mn}\\g_{mn}\end{pmatrix}}=T{\begin{pmatrix}a_{mn}\\b_{mn}\end{pmatrix}}}

La matriz T está determinada por la forma y el material del dispersor y, para un campo incidente dado, permite calcular el campo dispersado.

Cálculo de la matriz T

El método estándar para calcular la matriz T es el método del campo nulo , que se basa en las ecuaciones de Stratton-Chu. [ 6 ] Estas ecuaciones establecen que los campos electromagnéticos fuera de un volumen dado pueden expresarse como integrales sobre la superficie que encierra dicho volumen, considerando únicamente las componentes tangenciales de los campos en la superficie. Si el punto de observación se encuentra dentro de este volumen, las integrales se anulan.

Al utilizar las condiciones de contorno para las componentes del campo tangencial en la superficie del dispersor,

norte×(misdoat+miinortedo)=norte×miinortet{\displaystyle \mathbf {n} \times (\mathbf {E} _{scat}+\mathbf {E} _{inc})=\mathbf {n} \times \mathbf {E} _{int}}

y

norte×(Hsdoat+Hinortedo)=norte×Hinortet{\displaystyle \mathbf {n} \times (\mathbf {H} _{scat}+\mathbf {H} _{inc})=\mathbf {n} \times \mathbf {H} _{int}},

dóndenorte{\displaystyle \mathbf {n} }Siendo el vector normal a la superficie del dispersor, se puede obtener una representación integral del campo disperso en términos de las componentes tangenciales de los campos internos en la superficie del dispersor. Se puede obtener una representación similar para el campo incidente.

Al expandir el campo interno en términos de SVWF y explotar su ortogonalidad en superficies esféricas, se llega a una expresión para la matriz T. La matriz T también puede calcularse a partir de datos de campo lejano. [ 7 ] Este enfoque evita los problemas de estabilidad numérica asociados con el método de campo nulo. [ 8 ]

En línea se pueden encontrar varios códigos numéricos para la evaluación de la matriz T.Archivado el 16 de diciembre de 2022 en Wayback Machine .

La matriz T se puede encontrar con métodos distintos al método del campo nulo y al método de condiciones de contorno extendidas (EBCM); por lo tanto, el término "método de la matriz T" es desafortunado.

La mejora de la matriz T tradicional incluye el Método de Matriz T de Incrustaciones Invariantes (IITM) de BR Johnson. [ 9 ] El código numérico de IITM fue desarrollado por Lei Bi, basado en el código EBCM de Michael I. Mishchenko . [ 3 ] [ 10 ] Es más potente que EBCM, ya que es más eficiente y aumenta el límite superior del tamaño de partícula durante el cálculo.

Referencias

  1. Waterman, PC (1965). "Formulación matricial de la dispersión electromagnética". Actas del IEEE . 53 (8). Instituto de Ingenieros Eléctricos y Electrónicos (IEEE): 805– 812. doi : 10.1109/proc.1965.4058 . ISSN 0018-9219 . 
  2. Waterman, Peter C. (1971). "Simetría, unitaridad y geometría en la dispersión electromagnética". Physical Review D. 3 ( 4): 825– 839. Bibcode : 1971PhRvD...3..825W . doi : 10.1103/PhysRevD.3.825 .
  3. 1 2 Mishchenko, Michael I. ; Travis, Larry D.; Mackowski, Daniel W. (1996). "Cálculos de la matriz T de la dispersión de la luz por partículas no esféricas: una revisión". Journal of Quantitative Spectroscopy and Radiative Transfer . 55 (5). Elsevier BV: 535– 575. doi : 10.1016/0022-4073(96)00002-7 . ISSN 0022-4073 . 
  4. Lakhtakia, Akhlesh (2018). El teorema de extinción de Ewald-Oseen y el método de condiciones de contorno extendidas, en: El mundo del electromagnetismo aplicado . Cham, Suiza: Springer. doi : 10.1007/978-3-319-58403-4_19 . ISBN 978-3-319-58403-4.
  5. Mishchenko, Michael I .; Travis, Larry D.; Lacis, Andrew A. (2002). Dispersión, absorción y emisión de luz por partículas pequeñas . Cambridge, Reino Unido: Cambridge University Press. ISBN 9780521782524.
  6. Stratton, JA; Chu, LJ (1939-07-01). "Teoría de la difracción de las ondas electromagnéticas". Physical Review . 56 (1). American Physical Society (APS): 99– 107. Bibcode : 1939PhRv...56...99S . doi : 10.1103/physrev.56.99 . ISSN 0031-899X . 
  7. Ganesh, M.; Hawkins, Stuart C. (2010). "Cálculos de la matriz T de dispersión electromagnética tridimensional" . Journal of Computational and Applied Mathematics . 234 (6): 1702– 1709. doi : 10.1016/j.cam.2009.08.018 .
  8. Ganesh, M.; Hawkins, Stuart C. (2017). "Algoritmo 975: TMATROM - Un software de modelo de orden reducido de matriz T". ACM Transactions on Mathematical Software . 44 : 9:1–9:18. doi : 10.1145/3054945 . S2CID 24838138 . 
  9. Johnson, BR (1988-12-01). "Enfoque de matriz T de incrustación invariante para la dispersión electromagnética" . Applied Optics . 27 (23): 4861– 4873. doi : 10.1364/AO.27.004861 . ISSN 2155-3165 . PMID 20539668 .  
  10. Bi, Lei; Yang, Ping; Kattawar, George W.; Mishchenko, Michael I. (2013-02-01). "Implementación eficiente del método de matriz T de incrustación invariante y el método de separación de variables aplicado a partículas no esféricas grandes e inhomogéneas" . Journal of Quantitative Spectroscopy and Radiative Transfer . 116 : 169–183 . doi : 10.1016/j.jqsrt.2012.11.014 . hdl : 2060/20140010884 . ISSN 0022-4073 . S2CID 11722624 .