Respuesta al impulso de un operador diferencial lineal no homogéneo
Este artículo trata sobre el enfoque clásico de las funciones de Green. Para una discusión moderna, vea la
solución fundamental .
En matemáticas , la función de Green es la respuesta al impulso de un operador diferencial lineal no homogéneo definido en un dominio con condiciones iniciales o condiciones de contorno específicas.
Esto significa que si L es el operador diferencial lineal, entonces
- la función de Green G es la solución de la ecuación LG = δ , donde δ es la función delta de Dirac ;
- la solución del problema de valor inicial Ly = f es la convolución ( G * f ), donde G es la función de Green.
A través del principio de superposición , dada una ecuación diferencial ordinaria lineal (EDO), L (solución) = fuente, primero se puede resolver L (verde) = δ s , para cada s , y teniendo en cuenta que, dado que la fuente es una suma de delta funciones , la solución es una suma de funciones de Green, así, por linealidad de L .
Las funciones de Green llevan el nombre del matemático británico George Green , quien desarrolló el concepto por primera vez en la década de 1820. En el estudio moderno de ecuaciones diferenciales parciales lineales , las funciones de Green se estudian en gran medida desde el punto de vista de las soluciones fundamentales .
En la teoría de muchos cuerpos , el término también se usa en física , específicamente en teoría cuántica de campos , aerodinámica , aeroacústica , electrodinámica , sismología y teoría de campos estadísticos , para referirse a varios tipos de funciones de correlación , incluso aquellas que no se ajustan a la definición matemática. . En la teoría cuántica de campos, las funciones de Green asumen el papel de propagadores .
Definición y usos
Una función de Green, G ( x , s ) , de un operador diferencial lineal que actúa sobre distribuciones sobre un subconjunto del espacio euclidiano , en un punto s , es cualquier solución de
-
|
|
( 1 )
|
donde δ es la función delta de Dirac . Esta propiedad de la función de Green se puede aprovechar para resolver ecuaciones diferenciales de la forma
-
|
|
( 2 )
|
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 una función de Green única. Las funciones de Green pueden clasificarse, 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.
Las funciones de Green también son herramientas útiles para resolver ecuaciones de onda y ecuaciones de difusión . En mecánica cuántica , la función de Green del hamiltoniano es un concepto clave con vínculos importantes con el concepto de densidad de estados .
La función de Green, tal como se usa en física, generalmente se define 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.
Si el operador es invariante a la traslación , es decir, cuando tiene coeficientes constantes con respecto a x , entonces la función de Green puede tomarse como un núcleo de convolución , es decir,
En este caso, la función de Green es la misma que la respuesta al impulso de la teoría del sistema lineal invariante en el tiempo .
Motivación
En términos generales, si se puede encontrar dicha función G para el operador , 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, lo que produce
Esto significa que
-
|
|
( 3 )
|
es una solución a la ecuación
Por lo tanto, se puede obtener la función u ( x ) mediante el conocimiento de la función de Green en la ecuación ( 1 ) y el término fuente en el lado derecho en la ecuación ( 2 ). Este proceso se basa en la linealidad del operador .
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 G de Green que satisfaga la ecuación ( 1 ). Por esta razón, la función de Green también se denomina a veces la solución fundamental asociada al operador .
No todos los operadores admiten la función de Green. La función de Green también se puede considerar como un inverso a la derecha de . Aparte de las dificultades de encontrar la 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 da un resultado teóricamente exacto.
Esto se puede pensar como una expansión de f de acuerdo con una función delta base de Dirac (proyectando f sobre ; y una superposición de la solución en cada proyección . Esta ecuación integral se conoce como una ecuación integral de Fredholm , cuyo estudio constituye Fredholm teoría .
Funciones de Green para resolver problemas de valores de frontera no homogéneos
El uso principal de las funciones de Green en matemáticas es resolver problemas de valores de frontera no homogéneos . En la física teórica moderna , las funciones de Green también se utilizan normalmente como propagadores en los diagramas de Feynman ; el término función de Green se usa a menudo para cualquier función de correlación .
Estructura
Sea el operador de Sturm-Liouville , un operador diferencial lineal de la forma
y sea el operador de
condiciones de contorno con valores vectoriales
Sea una función continua en Supongamos además que el problema
es "regular", es decir, la única solución para todo x es .
Teorema
Hay una y solo una solución que satisface
y está dado por
donde es una función de Green que satisface las siguientes condiciones:
-
es continuo en y .
- Para , .
- Para , .
-
Derivado de "salto": .
- Simetría: .
Funciones de Green avanzadas y retardadas
A veces, la función de Green se puede dividir en una suma de dos funciones. Uno con la variable positiva (+) y el otro con la variable negativa (-). Estas son las funciones de Green avanzadas y retardadas, y cuando la ecuación en estudio depende del tiempo, una de las partes es causal y la otra anticausal. En estos problemas, por lo general, la parte causal es la importante. Estas son con frecuencia las soluciones a la ecuación de ondas electromagnéticas no homogénea .
Encontrar las funciones de Green
Unidades
Si bien no fija de manera única la forma que tomará la función de Green, realizar un análisis dimensional para encontrar las unidades que debe tener la función de Green es una importante verificación de cordura en 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 sólo de las unidades de sino también del número y las unidades del espacio del cual los vectores de posición y son elementos. Esto conduce a la relación:
donde se define como "las unidades físicas de ", 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 de 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 está completo, entonces es posible construir un Función de Green a partir de estos autovectores y autovalores .
"Completo" significa que el conjunto de funciones {Ψ n } satisface la siguiente relación de completitud ,
Entonces lo siguiente es válido,
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 asumió.
El estudio general de la función de Green escrito en la forma anterior, y su relación con los espacios funcionales formados por los vectores propios, se conoce como teoría de Fredholm .
Hay varios otros métodos para encontrar las funciones de Green, incluido el método de imágenes , la separación de variables y las transformadas de Laplace .
Combinando las funciones de Green
Si el operador diferencial se puede factorizar como, entonces la función de Green 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 , de forma análoga a cómo el operador lineal invertible , definido por , está representado por sus elementos de matriz .
Una identidad más 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 pueda factorizar, poniéndolo en la forma:
donde están los ceros de . Tomando la transformada de Fourier de con respecto a ambos y da:
La fracción entonces se puede dividir en una suma usando una descomposición en fracciones parciales antes de volver de transformación de Fourier a y espacio. Este proceso produce identidades que relacionan integrales de las funciones de Green y sumas de las mismas. Por ejemplo, si entonces una forma de la 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).
Tabla de funciones de Green
La siguiente tabla ofrece un resumen de las funciones de Green de los operadores diferenciales que aparecen frecuentemente, en donde , , es la función escalón de Heaviside , es una función de Bessel , es una función de Bessel modificada de primera clase , y es una función de Bessel modificada del segundo tipo . Donde aparece el tiempo ( t ) en la primera columna, se enumera la función de Green avanzada (causal).
Operador diferencial L
|
Función de Green G
|
Ejemplo de aplicacion
|
|
|
|
|
|
|
|
|
|
dónde
|
con
|
Oscilador armónico subamortiguado 1D
|
dónde
|
con
|
Oscilador armónico sobreamortiguado 1D
|
dónde
|
|
Oscilador armónico 1D críticamente amortiguado
|
Operador de Laplace 2D
|
con
|
Ecuación de Poisson 2D
|
Operador de Laplace 3D
|
con
|
Ecuación de poisson
|
Operador de Helmholtz
|
|
Ecuación 3D estacionaria de Schrödinger para partículas libres
|
en dimensiones |
|
Potencial de Yukawa , propagador de Feynman
|
|
|
Ecuación de onda 1D
|
|
|
Ecuación de onda 2D
|
Operador D'Alembert
|
|
Ecuación de onda 3D
|
|
|
Difusión 1D
|
|
|
Difusión 2D
|
|
|
Difusión 3D
|
|
con
|
Ecuación 1D de Klein-Gordon
|
|
con
|
Ecuación 2D de Klein-Gordon
|
|
con
|
Ecuación 3D de Klein-Gordon
|
|
con
|
ecuación del telegrafista
|
|
con
|
Conducción de calor relativista 2D
|
|
con
|
Conducción de calor relativista 3D
|
Funciones de Green para el laplaciano
Las funciones de Green para operadores diferenciales lineales que involucran al laplaciano pueden usarse fácilmente usando la segunda de las identidades de Green .
Para derivar el teorema de Green, comience con el teorema de divergencia (también conocido como teorema de Gauss ),
Dejemos y sustituya en la ley de Gauss.
Calcule y aplique la regla del producto para el operador ∇,
Conectando esto al teorema de divergencia produce el teorema de Green ,
Suponga que el operador diferencial lineal L es el laplaciano , ∇², y que existe una función de Green G para el laplaciano. La propiedad definitoria de la función de Green todavía se mantiene,
Deje entrar la segunda identidad de Green , vea las identidades de Green . Luego,
Usando esta expresión, es posible resolver la ecuación de Laplace ∇ 2 φ ( x ) = 0 o la ecuación de Poisson ∇ 2 φ ( x ) = - ρ ( x ), sujeto a las condiciones de frontera de Neumann o Dirichlet . En otras palabras, podemos resolver para φ ( x ) en todas partes dentro de un volumen donde (1) el valor de φ ( x ) se especifica en la superficie límite del volumen (condiciones de frontera de Dirichlet), o (2) la derivada normal de φ ( x ) se especifica en la superficie delimitante (condiciones de contorno de Neumann).
Suponga que el problema consiste en resolver φ ( x ) dentro de la región. Entonces la integral
se reduce a simplemente φ ( x ) debido a la propiedad definitoria de la función delta de Dirac y tenemos
Esta forma expresa la propiedad bien conocida de las funciones armónicas , que si el valor o la derivada normal se conoce en una superficie límite, entonces el valor de la función dentro del volumen se conoce en todas partes .
En electrostática , φ ( x ) se interpreta como potencial eléctrico , ρ ( x ) como densidad de carga eléctrica y la derivada normal como componente normal del campo eléctrico.
Si el problema es resolver un problema de valor límite de Dirichlet, la función de Green debe elegirse de manera que G ( x , x ′) desaparezca cuando x o x ′ está en la superficie límite. Por tanto, sólo queda uno de los dos términos de la integral de superficie . Si el problema es resolver un problema de valor en la frontera de Neumann, la función de Green se elige de manera que su derivada normal se desvanezca en la superficie limítrofe, ya que parecería ser la opción más lógica. (Ver Electrodinámica clásica de Jackson JD, página 39). Sin embargo, la aplicación del teorema de Gauss a la ecuación diferencial que define la función de Green produce
lo que significa que la derivada normal de G ( x , x ′) no puede desaparecer en la superficie, porque debe integrarse a 1 en la superficie. (Nuevamente, vea Electrodinámica clásica de Jackson JD, página 39 para este y el siguiente argumento).
La forma más simple que puede tomar la derivada normal es la de una constante, a saber, 1 / S , donde S es el área de la superficie. El término superficial en la solución se convierte en
donde es el valor medio 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í.
Sin condiciones de contorno, la función de Green para el Laplaciano ( función de Green para la ecuación de Laplace de tres variables ) es
Suponiendo que la superficie limítrofe llega al 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
Ejemplo
Encuentre la función de Green para el siguiente problema, cuyo número de función de Green es X11:
Primer paso: la función de Green para el operador lineal en cuestión se define como la solución a
-
|
|
(Ec. * )
|
Si , entonces la función delta da cero y la solución general es
Porque , la condición de frontera en implica
si y .
Porque , la condición de frontera 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 .
Asegurar la continuidad en la función del Green en implica
Uno puede garantizar discontinuidad apropiada en la primera derivada mediante la integración de la ecuación diferencial que define (es decir, la Ec. * ) A partir de a y tomando el límite cuando tiende a cero. Tenga en cuenta 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 y obtener
Entonces, la función de Green para este problema es:
Más ejemplos
- Deje que n = 1 y dejar que el subconjunto sea todos R . Deje que L sea . Entonces, la función escalón de Heaviside H ( x - x 0 ) es una función de Green de L en x 0 .
- Sea n = 2 y sea el subconjunto el cuarto de plano {( x , y ): x , y ≥ 0} y L el laplaciano . Además, suponga que se impone una condición de límite de Dirichlet en x = 0 y que se impone una condición de límite de Neumann en y = 0 . Entonces la función del X10Y20 Green es
- Sea , y los tres son elementos de los números reales. Entonces, para cualquier función de reales a reales , 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 .
Ver también
Referencias
-
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: Publicaciones de Dover. ISBN 0-486-63947-9.
El capítulo 5 contiene una descripción muy legible del uso de las funciones de Green para resolver problemas de valores límite en electrostática.
-
Polianina, 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.
-
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 Fourier Analysis y sus aplicaciones . Serie de Matemáticas. Wadsworth y Brooks / Cole.
-
Cole, KD; Beck, JV; Haji-Sheikh, A .; Litkouhi, B. (2011). "Métodos para obtener las funciones de Green". Conducción de calor usando las funciones de Green . Taylor y Francis. págs. 101-148. ISBN 978-1-4398-1354-6.
-
Verde, 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 verdes diádicas de espacio infinito en electromagnetismo . Londres, Reino Unido / San Rafael, CA: IoP Science (Reino Unido) / Morgan y Claypool (Estados Unidos). Bibcode : 2018idgf.book ..... F .
enlaces externos