La solución del problema de valor inicial es la convolución ( ).
A través del principio de superposición , dada una ecuación diferencial ordinaria (EDO) lineal , primero se puede resolver , para cada s , y darse cuenta de que, dado que la fuente es una suma de funciones delta , la solución es también una suma de funciones de Green, por linealidad de L .
donde δ es la función delta de Dirac . Esta propiedad de una función de Green se puede aprovechar para resolver ecuaciones diferenciales de la forma
Si el núcleo de L no es trivial, entonces la función de Green no es única. Sin embargo, en la práctica, alguna combinación de simetría , condiciones de contorno y/u otros criterios impuestos externamente darán como resultado una función de Green única. Las funciones de Green pueden categorizarse, por el tipo de condiciones de contorno satisfechas, por un número de función de Green . Además, las funciones de Green en general son distribuciones , no necesariamente funciones de una variable real.
La función de Green, tal como se utiliza en física, suele definirse con el signo opuesto. Es decir,
esta definición no cambia significativamente ninguna de las propiedades de la función de Green debido a la uniformidad de la función delta de Dirac.
En términos generales, si se puede encontrar una función G para el operador L , entonces, si multiplicamos la ecuación 1 para la función de Green por f ( s ) y luego integramos con respecto a s , obtenemos,
Debido a que el operador es lineal y actúa solo sobre la variable x (y no sobre la variable de integración s ), se puede sacar al operador de la integración, obteniendo
Esto significa que
es una solución a la ecuación
De esta manera, se puede obtener la función u ( x ) a través del conocimiento de la función de Green en la ecuación 1 y el término fuente en el lado derecho de la ecuación 2. Este proceso se basa en la linealidad del operador L .
En otras palabras, la solución de la ecuación 2 , u ( x ) , se puede determinar mediante la integración dada en la ecuación 3. Aunque se conoce f ( x ) , esta integración no se puede realizar a menos que también se conozca G. El problema ahora radica en encontrar la función de Green G que satisface la ecuación 1. Por esta razón, la función de Green también se denomina a veces solución fundamental asociada al operador L.
No todos los operadores admiten una función de Green. Una función de Green también puede considerarse como una inversa derecha de L . Además de las dificultades para encontrar una función de Green para un operador en particular, la integral en la ecuación 3 puede ser bastante difícil de evaluar. Sin embargo, el método proporciona un resultado teóricamente exacto.
La función de Green no es necesariamente única, ya que la adición de cualquier solución de la ecuación homogénea a una función de Green da como resultado otra función de Green. Por lo tanto, si la ecuación homogénea tiene soluciones no triviales, existen múltiples funciones de Green. En algunos casos, es posible encontrar una función de Green que no se anule solo para , que se denomina función de Green retardada, y otra función de Green que no se anule solo para , que se denomina función de Green avanzada. En tales casos, cualquier combinación lineal de las dos funciones de Green también es una función de Green válida. La terminología avanzada y retardada es especialmente útil cuando la variable x corresponde al tiempo. En tales casos, la solución proporcionada por el uso de la función de Green retardada depende solo de las fuentes pasadas y es causal , mientras que la solución proporcionada por el uso de la función de Green avanzada depende solo de las fuentes futuras y es acausal. En estos problemas, a menudo ocurre que la solución causal es la físicamente importante. El uso de la función de Green avanzada y retardada es especialmente común para el análisis de soluciones de la ecuación de onda electromagnética no homogénea .
Encontrar las funciones de Green
Unidades
Si bien no fija de manera única la forma que adoptará la función de Green, realizar un análisis dimensional para encontrar las unidades que debe tener una función de Green es una comprobación de cordura importante para cualquier función de Green encontrada por otros medios. Un examen rápido de la ecuación definitoria
muestra que las unidades de dependen no solo de las unidades de sino también del número y las unidades del espacio del que son elementos los vectores de posición y . Esto conduce a la relación:
donde se define como "las unidades físicas de " [ se necesita más explicación ] , y es el elemento de volumen del espacio (o espacio-tiempo ).
Por ejemplo, si y el tiempo es la única variable entonces:
Si , el operador d'Alembert , y el espacio tiene 3 dimensiones entonces:
Expansiones de valores propios
Si un operador diferencial L admite un conjunto de vectores propios Ψ n ( x ) (es decir, un conjunto de funciones Ψ n y escalares λ n tales que L Ψ n = λ n Ψ n ) que es completo, entonces es posible construir una función de Green a partir de estos vectores propios y valores propios .
"Completo" significa que el conjunto de funciones {Ψ n } satisface la siguiente relación de completitud ,
Entonces se cumple lo siguiente:
donde representa la conjugación compleja.
La aplicación del operador L a cada lado de esta ecuación da como resultado la relación de completitud que se supuso.
El estudio general de la función de Green escrita en la forma anterior, y su relación con los espacios de funciones formados por los vectores propios, se conoce como teoría de Fredholm .
Si el operador diferencial se puede factorizar como entonces la función de Green de se puede construir a partir de las funciones de Green para y :
La identidad anterior se sigue inmediatamente de tomar como la representación del operador derecho inverso de , análogo a cómo para el operador lineal invertible , definido por , se representa por sus elementos de matriz .
Otra identidad sigue para los operadores diferenciales que son polinomios escalares de la derivada, . El teorema fundamental del álgebra , combinado con el hecho de que conmuta consigo mismo , garantiza que el polinomio se puede factorizar, poniendo en la forma:
donde son los ceros de . Tomando la transformada de Fourier de con respecto a ambos y da:
La fracción se puede dividir en una suma utilizando una descomposición en fracciones parciales antes de volver a la transformada de Fourier en el espacio y . Este proceso produce identidades que relacionan integrales de funciones de Green y sumas de las mismas. Por ejemplo, si entonces una forma para su función de Green es:
Si bien el ejemplo presentado es manejable analíticamente, ilustra un proceso que funciona cuando la integral no es trivial (por ejemplo, cuando es el operador en el polinomio).
Las funciones de Green para operadores diferenciales lineales que involucran al laplaciano se pueden utilizar fácilmente utilizando la segunda de las identidades de Green .
Supongamos que el operador diferencial lineal L es el laplaciano , ∇ 2 , y que existe una función de Green G para el laplaciano. La propiedad definitoria de la función de Green sigue siendo válida,
Usando esta expresión, es posible resolver la ecuación de Laplace ∇ 2 φ ( x ) = 0 o la ecuación de Poisson ∇ 2 φ ( x ) = − ρ ( x ) , sujetas a condiciones de contorno de Neumann o de Dirichlet . En otras palabras, podemos resolver φ ( x ) en todas partes dentro de un volumen donde (1) el valor de φ ( x ) está especificado en la superficie límite del volumen (condiciones de contorno de Dirichlet), o (2) la derivada normal de φ ( x ) está especificada en la superficie límite (condiciones de contorno de Neumann).
Supongamos que el problema consiste en resolver φ ( x ) dentro de la región. Entonces la integral
se reduce simplemente a φ ( x ) debido a la propiedad definitoria de la función delta de Dirac y tenemos
Esta forma expresa la conocida propiedad de las funciones armónicas , que si se conoce el valor de la derivada normal en una superficie límite, entonces el valor de la función dentro del volumen se conoce en todas partes .
Si el problema consiste en resolver un problema de valor límite de Dirichlet, la función de Green debe elegirse de modo que G ( x , x ′) se anule cuando x o x′ se encuentren en la superficie límite. Por lo tanto, solo permanece uno de los dos términos en la integral de superficie . Si el problema consiste en resolver un problema de valor límite de Neumann, podría parecer lógico elegir la función de Green de modo que su derivada normal se anule en la superficie límite. Sin embargo, la aplicación del teorema de Gauss a la ecuación diferencial que define la función de Green da
como resultado que la derivada normal de G ( x , x ′) no puede anularse en la superficie, porque debe integrarse a 1 en la superficie. [3]
La forma más simple que puede adoptar la derivada normal es la de una constante, es decir, 1/ S , donde S es el área de la superficie. El término de superficie en la solución se convierte
en donde es el valor promedio del potencial en la superficie. Este número no se conoce en general, pero a menudo no es importante, ya que el objetivo suele ser obtener el campo eléctrico dado por el gradiente del potencial, en lugar del potencial en sí.
Suponiendo que la superficie límite se extiende hasta el infinito y conectando esta expresión para la función de Green finalmente se obtiene la expresión estándar para el potencial eléctrico en términos de densidad de carga eléctrica como
Primer paso: La función de Green para el operador lineal en cuestión se define como la solución de
Si , entonces la función delta da cero y la solución general es
Para , la condición de contorno en implica
si y .
Para , la condición de contorno en implica
La ecuación de se omite por razones similares.
Para resumir los resultados hasta ahora:
Segundo paso: La siguiente tarea es determinar y .
Garantizar la continuidad de la función de los Verdes implica
Se puede asegurar la discontinuidad adecuada en la primera derivada integrando la ecuación diferencial definitoria (es decir, Ec. * ) de a y tomando el límite cuando tiende a cero. Nótese que solo integramos la segunda derivada ya que el término restante será continuo por construcción.
Las dos ecuaciones de (dis)continuidad se pueden resolver para y para obtener
Entonces la función de Green para este problema es:
Más ejemplos
Sea n = 1 y el subconjunto sea todo R . Sea L . Entonces, la función escalón de Heaviside Θ( x − x 0 ) es una función de Green de L en x 0 .
Sean , y los tres elementos de los números reales. Entonces, para cualquier función con una derivada -ésima que sea integrable en el intervalo : La función de Green en la ecuación anterior, , no es única. ¿Cómo se modifica la ecuación si se suma a , donde satisface para todos (por ejemplo, con )? Además, compare la ecuación anterior con la forma de una serie de Taylor centrada en .
^ En la jerga técnica, "regular" significa que sólo existe la solución trivial ( ) para el problema homogéneo ( ).
Referencias
^ Cole, KD; Beck, JV; Haji-Sheikh, A.; Litkouhi, B. (2011). "Métodos para obtener funciones de Green". Conducción de calor mediante funciones de Green . Taylor y Francis. págs. 101–148. ISBN 978-1-4398-1354-6.
^ algunos ejemplos tomados de Schulz, Hermann (2001). Physik mit Bleistift: das analytische Handwerkszeug des Naturwissenschaftlers (4. Aufl ed.). Fráncfort del Meno: alemán. ISBN 978-3-8171-1661-4.
^ Jackson, John David (14 de agosto de 1998). Electrodinámica clásica . John Wiley & Sons. pág. 39.
Bayin, SS (2006). Métodos matemáticos en ciencia e ingeniería . Wiley. Capítulos 18 y 19.
Eyges, Leonard (1972). El campo electromagnético clásico . Nueva York, NY: Dover Publications. ISBN 0-486-63947-9. El capítulo 5 contiene una explicación muy legible del uso de las funciones de Green para resolver problemas de valores límite en electrostática.
Polyanin, AD; Zaitsev, VF (2003). Manual de soluciones exactas para ecuaciones diferenciales ordinarias (2.ª ed.). Boca Raton, FL: Chapman & Hall/CRC Press. ISBN 1-58488-297-2.
Barton, Gabriel (1989). Elementos de las funciones de Green y propagación: potenciales, difusión y ondas . Publicaciones científicas de Oxford. Oxford: Nueva York: Clarendon Press; Oxford University Press. ISBN 978-0-19-851988-1. Libro de texto sobre la función de Green con pasos resueltos.
Polyanin, AD (2002). Manual de ecuaciones diferenciales parciales lineales para ingenieros y científicos . Boca Raton, FL: Chapman & Hall/CRC Press. ISBN 1-58488-299-9.
Mathews, Jon; Walker, Robert L. (1970). Métodos matemáticos de la física (2.ª ed.). Nueva York: WA Benjamin. ISBN 0-8053-7002-1.
Folland, GB Análisis de Fourier y sus aplicaciones . Serie de matemáticas. Wadsworth y Brooks/Cole.
Green, G (1828). Ensayo sobre la aplicación del análisis matemático a las teorías de la electricidad y el magnetismo . Nottingham, Inglaterra: T. Wheelhouse. páginas 10-12.
Faryad y, M.; Lakhtakia, A. (2018). Funciones de Green diádicas en el espacio infinito en electromagnetismo. Londres, Reino Unido / San Rafael, CA: IoP Science (Reino Unido) / Morgan y Claypool (EE. UU.). Código Bibliográfico : 2018idgf.book.....F. ISBN 978-1-68174-557-2.
Şeremet, VD (2003). Manual de funciones y matrices de Green . Southampton: WIT Press. ISBN 978-1-85312-933-9.
Introducción a la técnica de la función de Green en el desequilibrio de Keldysh por AP Jauho
Biblioteca de funciones de Green
Tutorial sobre las funciones de Green
Método de elementos de contorno (para obtener una idea de cómo se pueden utilizar las funciones de Green con el método de elementos de contorno para resolver problemas potenciales numéricamente) Archivado el 7 de febrero de 2012 en Wayback Machine.
En Citizendium
Videoconferencia del MIT sobre la función de Green