BIBLIOTECA CIENCIAS

EXACTAS Y NATURALES
Q*324
.V34

11111111111 IH
15/T130 JRSIDAD DE SONORA
DIVISIÓN DE CIENCIAS EXACTAS Y NATURALES
DEPARTAMENTO DE MATEMÁTICAS

"Diste • Green

73/1 NI

7b. 19,

presenta

Horacio Valenzuela Martínez

Hermosillo, Sonora, México Octubre del 2001.

UNIVERSIDAD DE SONORA
DIVISIÓN DE CIENCIAS EXACTAS Y NATURALES
DEPARTAMENTO DE MATEMATICAS

"Distr Green

as
presenta

Horacio Valenzuela Martínez

Hermosillo, Sonora, México Octubre del 2001.
 DEDICATORIA

A DIOS Nuestro señor por hacerme de mi un hombre de bien.

A mis Padres ALEJANDRO y ALICIA que con su empeño, dedicación y

comprensión supieron guiarme a lo largo de esta etapa de la vida.

A mis Hermanos ALEJANDRO y YUNIVA con quiénes escribí una de las

más bellas páginas de mi vida.

A mis Abuelos JUAN y REFUGIO por su apoyo moral y espiritual

A todos mis Tíos que de alguna o de otra forma me apoyaron y me ayudaron

para salir adelante.

A todos mis Maestros por su ayuda y dedicación que me brindaron en el

transcurso de mis estudios.
Contenido

PREFACIO iii

1 Elementos de la Teoría de Distribuciones 1

1.1 Introducción 2
1.1.1 Motivación 2
1.1.2 La Delta de Dirac 4
1.2 Distribuciones Definidas sobre la Recta Real 6
1.2.1 Definición y Ejemplos 7
1.2.2 Derivadas, Multiplicación por una Función Suave e Integrales 17
1.2.3 Argumentos Compuestos y Convolución 23
1.3 Distribuciones Multidimensionales 30
1.3.1 La Delta de Dirac Multidimensional 34
1.3.2 Operaciones Generales 38
1.3.3 Interpretación de 5(a) 40
1.3.4 Producto de Deltas 42
1.3.5 Parte Finita de una Integral Divergente 43
*Comentario 1 50

2 Funciones de Green 51
2.1 Introducción 52
2.1.1 Motivación 52
2.1.2 Solución Fundamental 57
2.1.3 Construcción y Propiedades 59
2.2 Funciones de Green en Dimensiones Mayores 65
2.2.1 Funciones de Green en dos Dimensiones 66
2.2.2 Método de la Función de Green 69
2.2.3 Funciones de Green en tres Dimensiones. 70
*Comentario 2 72
 CONTENIDO

•
3 Aplicaciones 73
3.1 Soluciones Generalizadas de Ecuaciones Diferenciales 74
3.1.1 Ecuaciones Diferenciales Ordinarias 74
3.1.2 Ecuaciones Diferenciales Parciales 78
3.2 Método de las Imágenes 81
3.2.1 Función de Green sobre un Semi-Plano 82
3.2.2 El Problema de Robin 84
3.2.3 Funcion de Green sobre un Círculo 85
3.2.4 Función de Green sobre una Esfera 88
3.3 Algunas Aplicaciones Interesantes 89
3.3.1 La Fórmula de Euler—Leibnitz 89
3.3.2 Teorema de la Divergencia Generalizado 92
3.3.3 Entre Espejos Te Veas 94
3.3.4 Diferentes Representaciones para la Solución de la Ecuación de Onda 95
3.3.5 La Fórmula de Kirchhoff 97
*Conclusión 101

BIBLIOGRAFIA 102
 IOTECA
5.;,:u..1¿.,)
SA,IEk Is1
14 fit A MI GRAN11FIA

PREFACIO

En el presente trabajo se cubre lo básico que se debe de saber acerca la Teoría

de Distribuciones y Funciones de Green, desde un punto de vista aplicado. Esto es,
en el material no se abordan los conceptos y resultados con el rigor requerido por
las matemáticas abstractas; las construcciones teórico-formales abstractas se reducen
prácticamente al mínimo Quizá esto pueda significar la desaprobación de algunos
matemáticos puros quienes tengan una formación más o menos buena en la teoría, pero
no es relevante puesto que el enfoque sigue más o menos la filosofía de las matemáticas
aplicadas y desde mi punto de vista puedo asegurar que hacerlo de esta forma no resta
méritos ni mucho menos ya que el objetivo es presentar de la manera más clara y sencilla
el material elegido, esperando sea de utilidad tanto para consulta como para referencia
futura. Tal vez este comentario parezca un poco adelantado o fuera de lugar, sin
embargo lo hago en base a que aún cuando hasta la fecha se han publicado gran cantidad
de artículos y textos (incluyendo trabajos de tesis) sobre el tema, muchos de los cuales
teniendo presentaciones densas y extremadamente abstractas, a final de cuentas además
de ser de difícil lectura para los principiantes prácticamente no enseñan gran cosa. En
particular, me refiero a que en muchas de estas publicaciones las funciones generalizadas
multidimensionales, presentes en diversas aplicaciones reales, son tratadas como objetos
raros y curiosos, la mayoría de las veces difíciles de comprender tanto por la notación
utilizada como por el enfoque riguroso con el que se tratan, omitiendo en la mayoría
de los casos aspectos verdaderamente importantes. Por supuesto que hay excepciones,
pero también es cierto que muchos de estos materiales son difíciles de conseguir como
sucede a menudo al no estar disponibles en bibliotecas cercanas o accesibles.

iii
iv

e
Además, una buena razón para seguir el enfoque elegido se debe a que la teoría
de distribuciones es interesante, pero en toda su extensión, es una muy compleja,
sútil y difícil rama del análisis matemático, que para su formalización requiere bases
realmente sofisticadas, comúnmente difíciles de tener a nivel licenciatura. Sin embargo,
las funciones generalizadas pueden tratarse de forma relativamente fácil de comprender
y utilizar si inicialmente la teoría no se trata con mucha abstracción; con el enfoque
aplicado, para estudiarlas con cierta libertad solo es necesario tener conocimiento
elemental de geometría diferencial, análisis (real y complejo, vectorial, funcional y
tensorial) y quizá un poco de topología, la cual es necesaria, pero no un obstáculo. Por
consiguiente, el presente trabajo se dedica a presentar rudimentariamente la teoría de
distribuciones, enfatizando la aplicación, intentando con esto despertar el interés hacia
esta rama "olvidada" de las matemáticas, revelando de cierto modo la potencia de
la teoría mediante ejemplos concretos e ilustrativos, ayudando así a entender aunque
sea un poco los conceptos abstractos inherentes en las funciones generalizadas. El
tema clave o de enlace es el manejo de la distribucion delta de Dirac y distribuciones
multidimensionales relacionadas; de hecho estas son las distribuciones que aparecen
con más frecuencia en la mayoría de los problemas generalizados encontrados en la
práctica. Además, el camino a seguir hasta llegar a resolver ciertos problemas clásicos
de cálculo y de ecuaciones diferenciales es el estudio de las Funciones de Green, las
cuales aparecen al introducir la función delta de Dirac como término no homogéneo en
ciertas ecuaciones diferenciales, al momento de determinar la solución fundamental de
una ecuación diferencial.
En la exposición, trato de ser claro y sencillo, probando algunos resultados, esta-
bleciendo otros y formulando suposiciones explícitamente de tal forma que el trabajo
resulte lo más simple, pero ilustrativo posible.
Agradezco a la Universidad de Sonora, en específico al Departamento de Matemá-
ticas, la oportunidad de presentar este trabajo para obtener el grado de licenciatura. En
especial, agradezco al M.C. Rodrigo González González el apoyo tanto en las materias
que impartió y de las que formé parte, las cuales contribuyeron bastante a la realización
del presente trabajo, así como el apoyo brindado al dirigir el mismo.

Hermosillo, Sonora, México Horacio Valenzuela Martínez.

Octubre del 2001.
 Capítulo 1
Elementos de la Teoría de
Distribuciones

Desde los 1950's, cuando Laurent Schwartz publicó su Theories des Distributions,
las funciones generalizadas han tenido mucha aplicación en varios campos de las ciencias
e ingeniería y tal vez uno de los aspectos más útiles de esta teoría sea el hecho que
las funciones discontinuas pueden manipularse con la misma facilidad y libertad que
las funciones continuas o diferenciables, lo cual proporciona una poderosa herramienta
tanto en la formulación como en la solución de muchos problemas de aplicación.
En este capítulo, abusando de la terminología definiremos a las funciones generali-
zadas como funcionales lineales y continuos sobre cierto espacio especial de funciones,
que llamaremos funciones de prueba, después introduciremos el concepto de diferen-
ciación generalizada, de fundamental importancia para muchos aspectos aplicados y
finalmente presentaremos con el nuevo enfoque algunos resultados del análisis clásico,
los cuales tienen un gran impacto sobre muchas áreas de las matemáticas, en particular
sobre las ecuaciones diferenciales parciales.
Para comprender la teoría de distribuciones es necesario un cambio de pensamiento
en la forma usual de ver a las funciones ordinarias; funciones localmente (Lebesgue)
integrables —funciones que tienen integral finita sobre cualquier conjunto finito. Esto
implica la necesidad de ampliar el marco teórico de las funciones de tal forma que
ahora las funciones ordinarias sean un subconjunto de las nuevas funciones al que
denominaremos funciones generalizadas regulares. Cualquier otro tipo de "función"
que se introduzca pertenecerá a la categoria de las funciones generalizadas singulares.

1
2 CAPÍTULO 1. ELEMENTOS DE LA TEORÍA DE DISTRIBUCIONES

1.1 Introducción

Esta sección simplemente se utiliza como plataforma de visualización de la teoría

de distribuciones, también llamadas funciones generalizadas. En la primer parte, como
un primer intento de introducción, se relaciona esta teoría con las distribuciones de
carga físicas y fuentes concentradas o instantáneas presentes en varias áreas; fuerzas
y momentos concentrados de la mecánica de sólidos, fuentes de calor localizadas en la
conducción de calor, masas puntuales en la teoría del potencial, fuerzas impulsivas en
la acústica y mecánica de impactos, fuentes y vórtices en la mecánica de fluidos, cargas
puntuales, por capas o volumétricas en la electrostática, entre otras. En la segunda
parte, se introduce la más importante y popular de las funciones generalizadas, la
"función" delta de Dirac. De hecho, esta entidad matemática recibe un tratamiento
especial a lo largo de todo el trabajo ya que es la que tiene el más amplio rango de
aplicación. Además, muchas de las propiedades importantes de las "nuevas funciones"
se pueden ejemplificar e ilustrar fácilmente con la función simbólica delta de Dirac.

1.1.1 Motivación

Con el fin de preparar el terreno para empezar de lleno con la teoría de distribu-
ciones, es buena idea introducir primeramente el concepto de distribución de carga, el
cual es muy utilizado por físicos e ingenieros. En la aplicación, frecuentemente se habla
de carga puntual, línea de carga, superficie de carga, volumen de carga y dipolos. Pero,
exactamente ¿cómo se puede dar un significado matemático preciso a tales términos?,
considerando la física de la situación. Es decir, ¿cómo se puede modelar y/o mani-
pular matemáticamente situaciones o procesos físicos que involucren términos de esta
naturaleza? Por ejemplo, cuando físicamente se trata de medir una carga de inmediato
se piensa en hacerlo mediante un experimento en el que se utilice algún instrumento
cuya lectura proporcione la medida de la carga contenida en una región cuyo tamaño
dependa precisamente de la clase de instrumento que se utiliza. Si la carga es puntual,
podría utilizarse un instrumento o un conjunto de instrumentos con la capacidad de
medir la carga en regiones cada vez más pequeñas conteniendo al punto de interés.
Pero, resulta que prácticamente llega un momento en el que ya no es posible tener un
instrumento físico que logre tal objetivo. Lo mismo sucede con cualquier otro tipo de
1.1. Introducción 3

carga. Por consiguiente, es conveniente buscar una forma indirecta de hacerlo y una
posible alternativa podría ser: proponer un conjunto determinado de instrumentos, de
tal forma que correspondiendo a cada uno de estos instrumentos, 0, exista un número,
Q(0) (la lectura del instrumento) que proporcione precisamente la medida de la carga.
Por ejemplo, para medir la carga en un sólido, con densidad volumétrica de carga p(x),
el número Q(0) podría obtenerse mediante la integral triple
Q(q5) = f p(x)0(x) dx, (1.1.1)

proporcionando una especie de promedio, para funciones 0 apropiadas, de la carga

volumétrica en el sólido. Las propiedades inmediatas que se observan en este proceso
son
(1.1.2)
Q(a0) = f P( x ) [a0(x)] dx = aQ(0),
donde a es un escalar y
Q(01 + 02 ) = f p(x)(01 (x) + 02 (x)] dx = Q(02 ) + Q(02 ), (1.1.3)
para cualquier par de instrumentos 02 y 02 . De hecho, en general siempre será necesario
suponer que los instrumentos a utilizar mínimamente tengan las propiedades (1.1.2) y
(1.1.3), lo cual matemáticamente significa que Q(0) es un funcional lineal sobre el con-
junto de instrumentos 0, los cuales pertenecerán a un conjunto especialmente elegido
de forma tal que si se utiliza un funcional integral por lo menos se asegure su convergen-
cia. Por supuesto, además deberán de cumplirse otras condiciones y restricciones tanto
físicas como matemáticas. De esta forma, la elección del conjunto de instrumentos y sus
operaciones lleva al establecimiento de un nuevo conjunto de "objetos" matemáticos
que originan cierta estructura, realmente digna de estudio; la teoría de distribuciones,
la cual dará respuesta a la pregunta anteriormente planteada, pero mucho más aún,
representará una poderosa herramienta en la solución de varios problemas importantes.
En efecto, esta construcción de nuevos objetos representará una ampliación del marco
teórico de las funciones ordinarias utilizadas en el análisis clásico y como tal deberá
incluir todas las operaciones y propiedades previamente establecidas Además, una vez
establecidos los principales elementos y resultados de la nueva estructura, el siguiente
paso será desarrollar un nuevo cálculo con sus elementos, inmediatamente tratar de
formalizarlo, para finalmente buscar su posible aplicación, y este precisamente es el
esquema que se sigue en el presente trabajo.


4 CAPÍTULO 1. ELEMENTOS DE LA TEORÍA DE DISTRIBUCIONES

•
1.1.2 La Delta de Dirac

La teoría de distribuciones es una rama de las matemáticas relativamente reciente,

muy importante en áreas de la aplicación física, la cual tomó su forma final aproxi-
madamente hacia la mitad del Siglo XX. Aún cuando Paul Dirac y Sergey Sobolev
(así como varios antecesores) contribuyeron bastante, la teoría de distribuciones fue
desarrollada de forma elegante y definitiva por Laurent Schwartz, con el fin de dar una
base matemática rigurosa para las "funciones" delta, las cuales durante años habían
estado siendo utilizadas de forma incorrecta y con impunidad por físicos e ingenieros,
quienes hablaban de y usaban el término función delta de Dirac en un sentido intuitivo.
Pero, de todos es bien sabido que utilizar algo intuitivamente lleva en la mayoría de
los casos a conclusiones falsas o paradojas extrañas
De hecho uno de los errores en la aplicación fue introducir la función delta como
un límite puntual de ciertas funciones regulares, lo cual por supuesto no es correcto.
En específico, se consideraba (y tal vez aún se considere) a la función delta de Dirac
como el límite, cuando e -4 0, de la familia de funciones rectangulares

f (x)

para 1x1 <e

para 1x1 >e €1 = €0=1 x

Fig.1 Representación Simple (INGENUA)

de la "función" delta de Dirac.
La observación inmediata es que cuando e O, los rectángulos cada vez son más
angostos, pero más altos, conservando áreas constantes, ya que para cada e,

fe (x) dx = 1. (1.1.4)
fut
Esto llevó a definir la función delta de Dirac como el límite puntual

I
8(x) ci-g Ó f e (x), (1.1.5)

el cual, como se puede ver fácilmente, es cero dondequiera excepto en el punto x = O,

donde es infinito.


1.1. Introducción 5

Así, incluso en nuestros días, la definición ordinaria de la función delta de Dirac que
frecuentemente se maneja es
para x O,
(1.1.6)
para > O.

Además, por supuesto con l a condición que f 5(x) dx = 1, que tenga la bien probada

propiedad de convolución f 5(x — a)0(x) dx = 0(a), que tenga unidades de frecuencia

cuando la variable sea el tiempo, [[ 5(t) ] = 1/T, y que sea la derivada de la función
escalón o salto unitario, llamada función de Heaviside, 5(x) = c±h(x),
dx
h
si x > 0,
h(x) { 01 :
si x < 0.
x
Fig.2 Función de Heaviside.

Todas estas propiedades lucen muy atractivas, sin embargo antes de poder utilizarlas
es necesario tener una definición precisa de la "función" delta de Dirac ya que a esta
altura surgen las preguntas de rigor: ¿cómo puede una función tomar como valor a
oo y anularse en todos los demás puntos? y más aún, ¿cómo puede ser posible que
sea la derivada de una función con salto (discontinua) y además ser diferenciable?
Pero, el verdadero dilema no son estas preguntas, incluso podrían hasta omitirse, sino
que el asunto crítico surge desde el mismo momento en el que se intenta establecer el
límite puntual (1.1.5) como definición, ya que de ser así la integral (1.1.4) no está bien
definida si se considera en el sentido de Riemann y es igual a cero si se extiende como
una integral de Lebesgue. Sin embargo, observemos que para cada E > O, la integral

71[0] = f fe (x)0(x) dx (1.1.7)

existe para cualquier función continua 0, y cuando E O+, converge al valor de

esta función en cero, lim Te [0] = 71 [0] = 0(0), donde se denomina función de
prueba (el instrumento de medición que se mencionó anteriormente) la cual pertenence
a cierto conjunto por establecer. Así, una posible definición matemática correcta de
la "función" delta de Dirac puede ser mediante un proceso de este tipo. Es decir,
6 CAPÍTULO 1. ELEMENTOS DE LA TEORÍA DE DISTRIBUCIONES

como limite de una integral del tipo (1.1.7), en vez de un límite puntual de funciones
ordinarias como en (1.1.5). Esto es, utilizar para definirla un funcional lineal sobre el
conjunto de las funciones de prueba, generado por una función que determine todas
las propiedades del funcional a la que se le denominará el kernel del funcional.
Paul A. M. Dirac (1902-1984), uno de los principales creadores de la mecánica
cuántica y "padre" de la función delta, reconoció antes que muchos de sus contem-
poráneos la necesidad de emplear un enfoque funcional para este tipo de "funciones"
a las que llamó distribuciones. Por esta razón, de aquí en adelante la versión rigurosa
para la intuitiva función delta de Dirac será la de distribución delta de Dirac o sim-
plemente delta de Dirac. Esto para enfatizar que la "delta" no es una función. Por
su parte, el matemático ruso Sergey Sobolev empezó la creación formal de la teoría de
distribuciones en los 1930's (alrededor de 1936) al buscar soluciones generalizadas para
ecuaciones diferenciales parciales y contribuyó bastante. Poco tiempo después (1945),
el matemático Francés Laurent Schwartz completó el trabajo de estos fundamentos
al construir una estructura precisa para la teoría de distribuciones basándose en los
conceptos de espacios vectoriales topológicos localmente convexos. De hecho, Schwartz
ganó en 1950 la Field Medal (versión matemática del Premio Nobel) por este trabajo.
De esta forma, por primera vez desde Newton, las ideas acerca de las funciones y su
diferenciabilidad vinieron a tener una mayor revisión.

1.2 Distribuciones Definidas sobre la Recta Real

Las distribuciones, también conocidas con el nombre de funciones generalizadas
o funciones simbólicas, se introducen en esta parte de una manera un tanto infor-
mal, utilizando una dimensión e ilustrando con algunos ejemplos la razón del nombre
"generalizadas o simbólicas". Como hay mucho que decir y es necesario entender los
principales conceptos y resultados de esta teoría, esta sección resulta extensa. Además,
resulta interesante hacer una "copia" del cálculo elemental utilizando estos nuevos ob-
jetos matemáticos y una vez logrado esto, inmediatamente establecer la generalización
de una forma mucho más rápida y natural. Habrá muchas cosas que se escapen como
sucede en todo escrito y habrá otras que no queden del todo claras, sin embargo para
complementar ambas se puede encontrar una buena referencia en la bibliografía
1.2. Distribuciones Definidas sobre la Recta Real 7

1.2.1 Definición y Ejemplos

En Cálculo se relaciona a una función f (x) con una tabla de pares ordenados
(x, f (x)) y, aún cuando comúnmente esta tabla tiene un número infinito de elementos,
se establece que esta tabla representa a una curva plana a la que se le llama gráfica
de la función. En la teoría de las funciones generalizadas, también se describe a f (x)
con una tabla de números, los cuales son producidos por un funcional lineal T [O] que
se define mediante el proceso integral

T[P] Tf [q5] f f (x)q5(x) dx, (1.2.1)

generado por el kernel f(x), actuando sobre las funciones de prueba 0, las cuales
pertenecen a cierto conjunto D, por seleccionar. Una vez que se elije el conjunto de
funciones de prueba, V, el conjunto de funcionales lineales sobre este, llamado el espacio
dual de D, y denotado por D', se determina automáticamente. Precisamente, este tipo
de funcionales cuando son continuos son los que se identifican como distribuciones.
Por lo tanto, una función f(x) se puede describir ahora por una tabla de sus valores
funcionales sobre un conjunto dado de funciones de prueba. Así que primeramente se
debe de especificar el espacio de las funciones de prueba.

Afirmación: Entre más reducido sea el conjunto de funciones de prueba, más grande
será el conjunto de funcionales lineales definido sobre este y viceversa.
Como resultado es posible identificar todos y cada uno de los valores de f (x) a
partir de la tabla de valores {T[q5], E D}. Además, esta tabla debe ser tal que se
pueda distinguir dife- rentes funciones ordinarias en el sentido de Lebesgue; dos fun-
ciones ordinarias f y g son iguales casi dondequiera si difieren solo sobre un conjunto
de medida cero. Por lo tanto, como regla, para obtener un amplio conjunto de dis-
tribuciones tenemos que imponer restricciones estrictas y específicas sobre el conjunto
de funciones de prueba. Sin embargo, al mismo tiempo, el conjunto de funciones de
prueba no debe ser demasia- do reducido, ya que esto restringiría el rango de proble-
mas donde la herramienta de la teoría de distribuciones podría utilizarse, entre otras
complicaciones que implicaría.
8 CAPÍTULO 1. ELEMENTOS DE LA TEORÍA DE DISTRIBUCIONES

Ejemplo 1: El funcional que asigna a cada función de prueba su valor en cero se

identifica con la distribución delta de Dirac.

Nótese que en el sentido usual de funciones no existe ninguna función f (x) tal que
ff (x)0(x) dx 0(0), para toda función de prueba q5. No obstante, la representación
es muy conveniente y útil en el sentido que podemos introducir el término función
simbólica 8(x) de tal forma que

f5(x)0(x) dx = 0(0), (1.2.2)

quedando claro con esto que 8(x) no es una función de x (por ejemplo, no tiene ningún
sentido hablar de 8(5)). Lo mismo aplica para otras distribuciones más generales, que
más tarde identificaremos con el nombre de distribuciones singulares.

Ahora, respecto al conjunto de funciones de prueba D, existen algunas demandas

naturales. En particular, es conveniente identificar kerneles continuos f vía la definición
integral (1.2.1). En otras palabras, una vez que los valores del funcional T[0] son
conocidos, para toda tb E D, el kernel f debe determinarse de forma única. Es decir,
se requiere que el conjunto de distribuciones D' sea lo suficientemente rico como para
incluir por lo menos a todas las funciones continuas.

Espacio de las funciones generalizadas D'

funciones generalizadas
singulares
QE D
números reales
0(0) T[0] o complejos

Fig.3 Las Funciones Generalizadas D' son Funcionales Lineales y Continuos

sobre el Espacio Ti de las Funciones de Prueba.
1.2. Distribuciones Definidas sobre la Recta Real 9

Un buen candidato (tal como S. Sobolev lo eligió) para el conjunto de funciones

de prueba /3 es

D = E Cr}
= familia de funciones infinitamente diferenciables con soporte compacto.
(1.2.3)

Definición 1.1 Una función se dice ser de soporte compacto si es igual a cero fuera
de cierto conjunto acotado, denotado por Sopt 0, el cual se denomina soporte de cb y
se define como la cerradura del conjunto donde la función es diferente de cero. Esto es,

Sopt = {x ok(x) 0}. (1.2.4)

Nota: D tiene la estructura de espacio vectorial (lineal) y no es vacío trivialmente,

ya que por lo menos contiene a la función

C el/(x2-1) para ixj < 1
w(x) =

para x > 1.

Seleccionando C 2.25 (condición de

normalización), f w(x) dx = 1.
Fig.4 Función en Forma de Campana.
Si se cambia de escala para w, para cada e > O, se produce una nueva función
{1}, la cual está en D. Aplicando la fórmula de cambio de variable se
we (x) = 1w (É),
¡
tiene que we (x) dx = 1. A partir del Teorema del Valor Medio Generalizado para
integrales se tiene que si f (x) es una función continua, entonces cuando E —> O,

TE [fl = f f (x)o.),(x) dx 3 f (0).

Así, el valor de f en O puede recobrarse evaluando funcionales TE . El valor de f

en cualquier otro punto y puede recobrarse evaluando la integral sobre funciones de
prueba u) s trasladadas por y. Esto muestra la afirmación de que el valor de una función
continua en un punto cualquiera es determinado por los valores de un funcional sobre
las funciones de prueba elegidas y viceversa.


10 CAPÍTULO 1. ELEMENTOS DE LA TEORÍA DE DISTRIBUCIONES

Definición 1.2 Un funcional T[0] lineal y continuo sobre el conjunto D se llama

distribución y el conjunto de todas las distribuciones sobre D, denotado por D', el
espacio dual de 29, se denomina espacio de Sobolev-Schwartz.

(*) Un funcional T se dice que es lineal sobre 79 si T[aq5+0], ceT[0]+T[0], para

todo par 0, 2P E D y todo escalar a. Si el funcional está dado por la integral (1.2.1),
entonces su linealidad es inmediata.
(**) Un funcional T sobre D se dice que es continuo si para cualquier sucesión
de funciones {0 k (x)} c D que converge a q5(x) E D, los números T[h] convergen al
número T[q5]. La convergencia de la sucesión de funciones de prueba 0k en D a la
función de prueba O es en el siguiente sentido:
nü

i) los soportes de las O k 's están todos contenidos en un mismo conjunto acotado
del eje x,
ü) cuando k ce, las funciones q5k , y todas sus derivadas 0(: ) (x), n = 1, 2, ... ,
convergen uniformente a las derivadas correspondientes de la función límite 0(x).
Esto es, para cada n = 0,1, 2, ...
-+ 0(n)
0(kn)

Lo anterior genera un tipo especial de topología, que sería interesante estudiar,

así como la topología respectiva (más especial aún) del espacio dual generada por
, sucesiones débiles, a la cual se le llama topología débil. Sin embargo, ni el tiempo y ni
el espacio, disponibles para este trabajo, permiten entrar en este tipo de detalles.
Ejemplo 2: Algunas Distribuciones Básicas.

i) Consideremos el conjunto Lol = { f(x) I f es localmente integrable}, donde

localmente integrable significa que es integrable (en el sentido de Lebesgue) sobre
cualquier intervalo finito del eje x. Las funciones continuas sobre IR son la mejor
ilustración de estas. Así,

TÍO] (LO) = f f(x )O( x ) dx

está bien definido para cualquier q5 E D y es claro que define un funcional lineal
sobre D. También la continuidad de este funcional es fácil de verificar.
1.2. Distribuciones Definidas sobre la Recta Real 11

Una distribución definida de esta forma, con la ayuda de una función estándar
f E 4, se llama distribución regular y se denota por Tf , como en (1.2.1). Es
común también llamar al kernel f (x) función densidad y utilizar la notación
( f , ya que se puede considerar como el producto interior ordinario entre f y 0.
En general, podemos decir que las distribuciones son generalización de funciones
y productos interiores. Así, toda función localmente integrable f genera una
distribución a través del producto interior usual, (f, O) f (x)q5(x) dx. En este
contexto, decimos que las funciones localmente integrables pueden identificarse
con ciertas distribuciones en el espacio DI . Sin embargo, algunos funcionales
lineales y continuos (distribuciones) sobre D no pueden identificarse con kerneles
localmente integrables y entonces son denominados distribuciones singulares.

ü) El ejemplo más simple de una distribución singular es el funcional que asigna

a cada función de prueba 0 E D su valor en x = O (Ejemplo 1). Esto es, la
distribución delta de Dirac, a la cual nos referimos por 5[0] = 0(0). Cabe
señalar, para este caso, que algunas veces (1.2.1) se escribe heurísticamente como
f8(x)0(x) dx = 0(0), aunque formalmente ésto no representa la integración
en algún sentido (Riemann o Lebesgue). Sin embargo, sirve como una regla
nemotécnica intuitiva que, si se usa con cuidado, facilita muchos cálculos en
los que este involucrada la delta de Dirac, la cual no es una función, sino una
distribución singular y hablar del argumento x es conveniente ya que se puede
hablar de la distribución i5(x — e) definida por la fórmula

f(5(x — e)(b(x) dx = Çb(e),

algo así como trasladar e unidades la delta de Dirac (ver Subsección 1.2.3 para el
caso más general), la cual denotaremos por 5e y diremos la delta centrada en e.
Otros ejemplos importantes de distribuciones singulares son los dipolos o multi-
polos, los cuales se introducen posteriormente en la Subsección 1.2.3.

iii) En muchos contextos aparece el concepto de valor principal de Cauchy (Augustin

Louis Cauchy (1789-1857), fue uno de los matemáticos más prolíficos y creativos
de los "tiempos modernos") el cual se puede establecer de diferentes formas.
 12 CAPÍTULO 1. ELE
MENTOS DE LA TEORÍA DE DISTRIBUCIONES

Para empezar, se define la di stribución valor

p rincipal de Cauchy como el límite
simétrico de la integral singular f 0(7)
— dr. Por ejemplo, para
T—X x O,
PVf q—
5(7) dr lim [f e 4_ i ce] 0÷
40 r) dr, (1.2.5)
cuando el límite existe, con una definición
análoga para otros valores de x.

P
odríamos seguir dando ejemplos de distribuciones
imp ortantes, pero no es el
momento adecuado para hacerlo ya que se requiere
teoría. En el punto iii) de la siguiente pr ofundizar un poco más en la
o bservación se ilustra esta necesidad.

Obs ervación: i)
En (1.2.5) podemos escribir para fE
CO

Si 95 E D, entonces existe 1?
o. fe
tal que 0(x) O para lxj > I?
y

PV f d T V p f
Ahora, escribiendo q5(x) [95(x) — 0(0)]
+ 0(0), se obtiene un proceso

PV f &)_ dr f R 0 ) — 00 ) dr
-R
que remueve la si ngularidad en O. Es decir,
hecho que 0 E D satisface la c
regulariza al integrando, lo cual sigue del
ondición de L ipschitz,
aún cuando q5 no sea de soporte lq5(a)— q5(b)i < kja — bi. También,
compacto, tenemos que

l_ r) dr = f R ° 7) Ø(o) dr
(
o PV1 T E -*ce
-R

PV f
P IE) dr =
Jo
s er) dr
lo cual sigue de la fórmula de cambio de variable y del teorema del valor medio del
cálculo.
1.2. Distribuciones Definidas sobre la Recta Real 13

iii) Otro enfoque conveniente para el valor principal es utilizar una integral de
Cauchy,
J °(r) T —
dT = lim
y-.0+
f Çb(r)
T—Z
dT, (1.2.6)

donde z = x + iy es un parámetro complejo. Separando las partes real e imaginaria,

1 T —x iy
z (7- x)2 + y2 x)2 + y2
(r
y sustituyendo esta última relación en (1.2.6), obtenemos

hm
y—>0+ f 4)(T) dT = 19V
T— 1 q5(T)
T —
dT + i hm

Como el término que multiplica a en la integral del lado derecho es una curva de

(T — X) 2 + y2
0(r) dr.

Lorentz multiplicada por 7r, que converge débilmente a 8(x — 7) (ver Ejemplo 3b)),
entonces
°(T) = rv f °(7)
f— Z
dT
T—Z
dT f iluk(x),

donde el signo (+) significa que y -+ 0+ (a través del semi-plano superior) y el signo ( —)
cuando y 0_ (a través del semi-plano inferior). Esto determina las dos distribuciones

1 1
=PV x + iird(x — T),
— x — iO T—
(1.2.7)
1 1
T — X + i0 =PV
—
x a-8(x — 7).

Por ejemplo, una vez que se introduce el concepto de derivada generalizada (Sub-
sección 1.2.2), es posible mostrar que PV ( i= 1 ln Ix I.
( x dx

iv) También, es posible definir el concepto de soporte de una distribución; una

distribución T se dice que se anula en una región abierta 8 del eje x si T[01= 0, VOED
con Sopt c B. El complemento de la región más grande para la que la distribución
T se anula es a lo que se llama el soporte de la distribución. T y se denota por Sopt T.
Por ejemplo,
Sopt ó = {0}.
 14 CAPÍTULO 1. ELEMENTOS DE LA TEORÍA DE DISTRIBUCIONES

v) Además, como ya se observó en iii), para definir ciertas dsitribuciones es

necesario introducir el concepto de convergencia débil. Aún cuando la distribución
delta de Dirac formalmente no puede representarse como un funcional integral, esta
puede obtenerse como un límite de una sucesión de funcionales integrales,

Tk [0] f fk (x)0(x) dx, (1.2.8)

respecto a kerneles funciones regulares, como por ejemplo las funciones rectangulares
vistas anteriormente. En este sentido, muchas de las propiedades de la delta 5 pueden
11
derivarse de las propiedades de la sucesión {fk }. La aproximación de la 5 por fk en el
sentido que para cada E V,
Tkkb) --> b[1b] (1.2.9)
cuando k cc, es llamada aproximación débil y la convergencia correspondiente,
convergencia débil.
La elección de una sucesión de distribuciones {Tk }, representadas por funciones
regulares {fk }, no es única y se puede seleccionar de diferentes formas, incluso se
puede seleccionar de forma tal que las funciones fk (x) sean infinitamente diferenciables
(suaves, pero no necesariamente con soporte compacto).

Ejemplo 3: a) Considerar la familia de funciones Gaussianas,

1 e-x2 pe
Ye(x)
f2 ire
(paramentrizada por un parámentro e > O),
con la propiedad f eye (x) dx = 1.

Fig.5 Funciones Gaussianas, las cuales

convergen débilmente a la 8.

La sucesión fk (x) = 71/k (x), k = 1, 2, ... , cuando k co, aproxima en forma

débil a la delta de Dirac. Obsérvese que las campanas tienen picos cada vez más altos,
concentrándose alrededor de x O, mientras conservan el área bajo sus gráficas.

1.2. Distribuciones Definidas sobre la Recta Real 15

b) . Las curvas de Lorentz (físicos) o densidades de Cauchy (matemáticos),

As(x)

1 E
A e (x) = 7! x2 + E2

Fig.6 Curvas de Lorentz o Densidades de Cauchy.

son elementos de otra sucesión de funciones regulares que converge débilmente
a la delta de Dirac, al definir de nuevo los kerneles fk = Aiik . Nótese que
aunque a primera vista las funciones de Lorentz se parecen mucho a las funciones
Gaussianas, tienen ligeras diferencias; una de estas es que una función Gaussiana
decrece exponencialmente a O, cuando x —› loo, mientras que una función de
Lorentz decrece más lentamente y las áreas bajo las Gaussianas son mucho más
concentradas alrededor de x O.

c) Las funciones regulares que convergen débilmente a 5 no necesariamente concen-

tran sus valores diferentes de cero alrededor del origen x O y no necesariamente
tienen valores positivos, como en a) y b). Incluso existen funciones con argumen-
tos complejos cuyas partes real e imaginaria tienden débilmente a la delta de
Dirac, f (x)
1

f, (x) = e -2.22E
`Lile

e
fs (x) oscila a razones cada vez más grandes mientras que e se hace pequeño,
dando como resultado que Re {Te m} —> 0(0) cuando e -3 0.
16 CAPÍTULO 1. ELEMENTOS DE LA TEORÍA DE DISTRIBUCIONES

d)* Todavía más aún, las funciones regulares no necesariamente son pares o simétricas
Considerar
D(x)
-2 -1/x
X e para x > O,
D(x)
{ O para x < O.

1 x
La sucesión de funciones De (x) = — D E )
E

también aproxima débilmente a la delta.

Fig.7 Curvas de Dawson.

e) Así como estos ejemplos, existen muchas sucesiones de funciones regulares que
convergen débilmente a la delta de Dirac o a otro tipo de distribuciones singulares
(ver Ejemplo 4). Otros ejemplos importantes son algunas de las sucesiones que
Heaviside utilizó en sus cálculos,

{ n2 (! — lx0 , ixl 5 1/n

Tn (x) =
O, lx1 > 1/n

sen nt xa
Sn (t)= rt (n —> cx)), Fa =a —> —1),
F(a + 1)
las cuales también convergen débilmente a la delta de Dirac.

A continuación se enuncia uno de los principales resultados de la teoría de dis-

tribuciones (para su demostración referirse a [6]).

Teorema 1.2.1 El espacio TY es completo.

Este resultado nos dice que toda sucesión convergente de funciones generalizadas en
D' siempre convergerá a una función generalizada también en .
1.2. Distribuciones Definidas sobre la Recta Real 17

•
1.2.2 Derivadas, Multiplicación por una Función %ave e Integrales

La diferenciabilidad infinita de las funciones de prueba elegidas D nos permite

definir derivadas de orden arbitrario para cualquier distribución T E D', sin ningún
temor por la existencia o no existencia de la diferenciabilidad como el que se tiene para
las funciones ordinarias. Esta es una de las muchas ventajas que tienen las distribu-
ciones sobre las funciones elementales del cálculo. También, es posible definir nuevas
distribuciones a partir de otras distribuciones previamente conocidas, como por ejemplo
multiplicar por una función suave y determinar la primitiva de una distribución.

Antes de dar una definición general observemos que la familiar fórmula de inte-
gración por partes aplicada a una función diferenciable f (x) y una función de prueba
0(x) e D se reduce a

I(x)¢(x) dx = — f f(x)01 (x) dx (1.2.10)

ya que f (x)0(x) ree = O debido a que la función de prueba q5 es cero fuera de cierto
conjunto acotado en el eje x. Si consideramos la distribución asociada a la función
regular f, entonces (1.2.10) nos dice que

(711)195] = (1.2.11)

la cual es válida para cualquier función de prueba 0. Nótese que (1.2.11) no depende
de la diferenciabilidad de f . Por lo tanto, esta idea se puede extender a cualquier
distribución, dando como resultado la siguiente definición de diferenciabilidad.

Definición 1.3 Si T es una distribución en D', entonces su derivada T' es definida

como una distribución en D', determinada por sus valores (como un funcional) sobre
las funciones de prueba q5 E D por la relación

T'[0] = —T[0 11, (1.2.12)

la cual está bien definida ya que representa a un funcional lineal y continuo sobre D.
 18 CAPÍTULO 1. ELEMENTOS DE LA TEORÍA DE DISTRIBUCIONES

Las derivadas de orden superior se definen mediante la aplicación recursiva de la

operación de la primer derivada, dando como resultado que la n—ésima derivada de T
es de nuevo una distribución sobre 7) determinada por sus valores sobre las funciones
de prueba E D por la fórmula
71(n) [0] = ( - 1)nT[0(n)]. (1.2.13)
Así, en el sentido distribucional, siempre se tienen derivadas de todos los ordenes.
En mi opinión esta es una de las mejores propiedades que tienen las funciones gene-
ralizadas. De hecho, otra de las razones por la que surge la teoría de distribuciones es
precisamente por el deseo tanto de físicos como de matemáticos de emplear derivadas
más generales de funciones no diferenciables en el sentido del cálculo.

Ejemplo 4: Consideremos la distribución 6(x—e), definida por 5(x— e)[q5] 0(e). De

acuerdo a (1.2.13), su n—ésima derivada es 8(x — 0 (71) (0] = (-1) n Q5(n) (e). En particular,
su primera derivada 5' es un funcional sobre D, el cual con signo negativo describe a
un dipolo, otro ejemplo de distribución singular, el cual también puede obtenerse como
límite débil de sucesiones regulares, por ejemplo, si se consideran las derivadas de las
funciones Gaussianas -ye , se tiene que

— f 7É(x)¢(x) dx = x)0'(x)
f -ye( dx —> 5[q51 = cuando e —› O.

Ve(x)

Fig.8 Funciones aproximadoras al negativo del dipolo

I
Similarmente, se obtiene el mismo resultado si se utiliza cualquier otro tipo de aproxi-
madores de la delta.
1.2. Distribuciones Definidas sobre la Recta Real 19

Otra operación lineal, la cual produce una nueva distribución a partir de una
distribución dada T y una función infinitamente diferenciable a, es la multiplicación
de T por a.

Definición 1.4 El producto de una función a E C°° y una distribución T E D' es una
distribución en D' determinada por

(aT)[0]= T[(a0)], Çb E D. (1.2.14)

En el caso particular de una función constante a(x) = c, se tiene la nueva
distribución en D',
(1.2.15)
(cT)[0] = 71(0)1.

Ejemplo 5: Multiplicación de la delta (5(x — e) por una función suave,

(a (x) 6 (x ))[01 = 6(x = a (e)95 (1) = a (e) 5 ( x e)[01.

Es decir, a(x)5(x — e) , a(e)(5(x — e).

La diferenciación de distribuciones multiplicadas por funciones suaves satisface la

clásica regla de Leibniz para la derivada del producto de dos funciones.

Teorema 1.2.2 Si a E C°° y T E D', entonces (aT)' = a'T + aT'.

Dem. Tenemos que para cualquier función de prueba q5 E D,

(aT) 1 [0] = —(an[q511= = —TRacb)' — a/q51 = —T[(aq5) 1] +T[akb]

= T'[a41 + (a'T)[0]= (aT')[q5] + (a'T)[0] = (a'T + aT9[0].

Por lo tanto, (aT)' = a'T + aT'.

Similarmente, es fácil probar la fórmula de Leibniz generalizada, para distribuciones,

n a(k)T(n-k)
(aT) (n) =
k=0
 20 CAPÍTULO 1. ELEMENTOS DE LA TEORÍA DE DISTRIBUCIONES

Ejemplo 6: a) (a(x)8(x — e))/[0] = —a(e)01(e), ya que

( a ( x ) 8 ( x — e))1[0] = — ( a ( x ) 8 ( x e))[ 95 1 = — ( a (C) 8 ( x — IDA

= —a (C) 8 (x = —a(1)01(e).
En particular, (c8(x—e))1[0] = —0'(e) y en general (a(x)8(x —1)) (n) = a(e)8 (n) (x — e).
Es decir, cualquier función suave que multiplique a la delta de Dirac puede considerarse
como una constante que sale fuera de la operación diferenciación.

b) (a89[0] = 61a0] = —8[(a0)1= —8[ a'q5+ a0' = — a'(0-0(C) — a(e)0'(e). En

,111
particular, para a(x) = x y e = O, tenemos que

x61 (x)= —8(x). (1.2.16)

Es decir, la delta de dirac es una solución generalizada de la ecuación diferencial

ordinaria xy' = —y.
o

Nota: i) La ecuación (1.2.16), así como la ecuación más general

xn8(n) (x) = ( - 1)nn!5(x), (1.2.17)

no es posible establecerlas sin el enfoque de la delta de Dirac en el sentido distribucional.

A partir de aquí ya podemos observar la diferencia y la potencia de las herramientas
matemáticas modernas como lo es la teoría de distribuciones.
—6(x) 1
ii) 6'(x) ya que — no es una función de clase C°° Orlo está permitido
dividir a una distribución por una función que se anula en ciertos puntos!). Esto es, la
multiplicación de distribuciones no está bien definida. Por ejemplo, la multiplicación
de la delta 8(x) por si misma (es decir, 82 (x)) no tiene ningún significado en el sentido
de las distribuciones.
n-1

Ejemplo 7: Es buen ejercicio verificar que T E ck8 (k) (x) es una solución
k=o
generalizada de la ecuación xn 7" 0 (solución del problema de la división del cero).
1.2. Distribuciones Definidas sobre la Recta Real 21

Ahora, al igual como se hace en cálculo integral, consideremos el problema de

determinar una distribución S tal que S' = T. Es decir, encontrar la primitiva o
integral indefinida de una distribución. Por ejemplo, la función de Heaviside genera la
distribución de Heaviside
00
Tk [0] = f h(x)q5(x) dx = fo q5(x) dx.

En el sentido del análisis clásico h(x) no tiene derivada en cero. Sin embargo, su
derivada distribu-cional la cual de aquí en adelante denotaremos con una barrita () está
bien definida y está dada por
h'(x) = 5(x), (1.2.18)
00
ya que 711,[0] = — Tk [q51 ] = — f h(x)0'(x) dx = — f 0 1 (x) dx = q5(0) = 5[0]. Así,
se dice que la función de Heaviside es la primitiva de la delta de Dirac, en el sentido
generalizado o de distribuciones.
Una vez determinada la derivada de la función de Heaviside, inmediatamente se
puede calcular la derivada distribucional de cualquier función f (x) continua por piezas,
con discontinuidades de salto en los puntos xk , k = 1, 2, ... , m. Es evidente que una
función de este tipo siempre se puede expresar en la forma

f(x ) = f s (x) + ETro

k=1
lfk ih(x — x k ), (1.2.19)

donde [fk l = f(x k +0)— f (xk — O) denota el tamaño del salto en xk y f, es una función
continua suave por piezas.
f (x)

(x)

Xy Xy

Fig.9 La parte continua fs (x) puede obtenerse de f(x)

removiendo de esta los saltos.


22 CAPÍTULO 1. ELEMENTOS DE LA TEORÍA DE DISTRIBUCIONES

Por lo tanto como la derivada distribucional es una operación lineal, tenemos que
771

f' = { h}' + E Ifk1 .5 ( x - xk), (1.2.20)
k=1

donde {•} significa "únicamente la parte suave". Si consideramos el orden contrario

obtenemos la antiderivada de cualquier distribución que consta de una función local-
mente integrable más una combinación lineal de deltas centradas en los puntos de salto.
En el mismo sentido, tenemos la fórmula
711

f(n) = {fs}(n) + E E[fP-15(,i-i)(x-xk) (1.2.21)

k=1 j=1
para la n-ésima derivada de una función f continua por piezas.

Ejemplo 8: Como parte del trabajo, pero no se incluye el procedimiento para no

agotar espacio, fue interesante mostrar que

sgn' (x) = 28(x), donde sgn denota a la función signo,

para x > O
sgn (x) = {10 para x = O
—1 para x < O.
si f (x) = x 2 sgn (x), entonces

f' (x ) = 21x1, f" = 2 sgn (x), fm(x) = 48(x), etc.

c) f (x) = satisface a la ecuación diferencial ordinaria —y" + y = 8(x) y

v(x,()=
_2ilkeikix-C} es solución generalizada de y" + k2v = 8(x — ().

Ahora, vamos a finalizar esta sección enunciando el siguiente resultado (su de-
mostración se puede encontrar en [6, 11]), el cual es de gran importancia y de una
enorme potencia.
Teorema 1 2 3 (Teorema Estructural de la Teoría de Distribuciones): Las
funciones generalizadas en TY son derivadas generalizadas de orden finito de funciones
continuas.
1.2. Distribuciones Definidas sobre la Recta Real 23

4.
Observación: El concepto de integral definida para una distribución, es mucho más
complejo de manipular, está presente solo en ciertos contextos; para trabajar con éste
se necesitan algunas suposiciones y restricciones adicionales y por lo mismo, no se
considera en este trabajo.

1.2.3 Argumentos Compuestos y Convolución

En esta parte consideraremos otras dos operaciones con las cuales se generan nuevas
distribuciones a partir de otras ya dadas o construidas: la composición y convolución
de distribuciones.

Si f (x) es una función regular, entonces la fórmula de cambio de variable usual

nos establece que

f f (a (x))0(x) dx = f f (Y)0(,3 (0)101 (01 dy,

donde y = a(x), x = /3(y). Esto es, fi

Una suposición que garantiza la válidez de la fórmula de cambio de variable es

que a(x) sea una función estrictamente monótona y que aplique el eje x completo
en el eje y completo (función 1 - 1 y sobre). En el sentido funcional, se necesita que
.0([3(y))1fii(y) E V, lo cual se asegura si fi E C°°. Así, con esta condición,

71(a(x))[0(x)] 71[0(fi(Y))ifil(Y)11 (1.2.22)

determina la composición de una distribución T con una función a(x).

Ejemplo 9: a) Si a(x) = x — e (una traslación), entonces

5(a(x))[(b] = (5 [(b ( y + C)) = 95(e) = 5 (x e)[95]•

Es decir, 8(a(x)) = 8(x — e), la cual actúa sobre las funciones de prueba de la siguiente
manera

(5(x — e)H = 0(E).



24 CAPÍTULO 1. ELEMENTOS DE LA TEORÍA DE DISTRIBUCIONES

En cambio, 5(a(x) - e) = 5(x - fi())) ya que

la, (0(e)),

8 ( a ( x ) e)[0] = 5 (Y - e)[45 ( 3 (Y))1fii (Y)ii = 0(fi(e))ifii(e)i,

donde fi' (e) = 1/ci t (NI)) (por el teorema de la función inversa). En específico tenemos
por ejemplo que 8(-7x + 5) = Ti5(x -
1
También, tenemos la propiedad de auto-similaridad de la delta 8(cx) = — 6(x)
ici
y con esto vemos que 8(-x) 8(x), algo así como que la delta es una "función" par
o simétrica.

En general, si T[b(-x)] = T[0(x)], se dice que la distribución es par y si T[0(-x)] =

-T[0(x)], entonces impar.
Bajo las mismas condiciones para a(x), la fórmula de cambio de variable (1.2.22)
es equivalente a

f f (a(x))q5(x) dx = f f(y ) [4 faw<y dej dy, (1.2.23)

la cual tiene sentido si 0(e) de E D, siendo este el caso si la función inversa

dy lY(E)<Y
de a(x) = y es C°° y de lento crecimiento en el infinito Por ejemplo, si a(x) = cx-b,
la integral (1.2.23) se reduce a

1,,
f f (cx - b)0(x) dx = f f (y)[ (5 L_bcy 0(0 del dy.
fx
J
Poniendo . icb(e) de = 1. (x), tenemos que

(y+b)/c
+ b
0(e) de = 11( ) 11 ( 00) si c > O,
f 00
0(e) de =
fce—b<y
y b
0(e) de 100) <1( ) si c < O.
r
ti(y+
b)lc
Así, para cualquier valor de c O, tenemos

0(e) de = (Y ± b)
dy ce_b<1,
1.2. Distribuciones Definidas sobre la Recta Real 25

Sustituyendo esta relación en la integral correspondiente a (1.2.23), obtenemos

J f (cx — a)¢(x) dx = icr 1 ff (y)0 ( Y +e a ) dy,

resultado obtenido al aplicar directamente la fórmula de cambio de variable del cálculo.

La ventaja de la fórmula (1.2.23) es que se puede aplicar aún si el argumento es una
función no monótona. Esto ejemplifica que hay ocasiones en las que es conveniente
abandonar el enfoque riguroso de la teoría de distribuciones. Por ejemplo, (1.2.23) se
2 2
puede utilizar para establecer 8(x — a ). Calcular

f8(x 2 — a 2 )q5(x) dx = f 8 ( y ) [dy ft2 a2<ii 0(e) de] dy,

pero

pvt4-Fa2
f ea' <y
owde=
-N/Y+a2
«ocie= 4. ( ‘,/y + a2 'D ( — Vy + a2),

de donde
d 1 r
«I) = 2vy ai vfi(Vy + a2 ) — 0(-4 ± az)]
dy je2_ 6,2 <y
y por lo tanto,

f8(x2 — a 2 )0(x) dx = I 5(y)

r
[0(N/Y + a2 ) + 0(—A4 + a2 )] dy
\,/y + a2

= 211a1[0(a) +'P(—a)]

Como 1 1
[0(a) + O( —a)] = 2{4 [8(x — a) + 8(x + a)] 0(x) dx,
21a1 f

se concluye que

6(x2 a2 ) = 2 1 ai [5(x — a) + 5(x + a)] .

 26 CAPÍTULO 1. ELEMENTOS DE LA TEORÍA DE DISTRIBUCIONES

• En general, cuando se opera la delta de Dirac con argumentos no monótonos,

8(a(x) — y), donde la función y = a(x) tiene raíces simples o múltiples, x k =
k 1, 2, ... , n(y),

Fig.10 Argumento para la delta de Dirac no monótono.

necesariamente se tiene que abandonar el enfoque riguroso de la teoría de distribuciones,

resultando de gran utilidad la fórmula (1.2.23). Este es el caso cuando nos encontramos
con la fórmula
= 8(x — xk)
8(a(x)

o en vista de las propiedades de la delta de Dirac,

8(a(x) — y) = 1 (5(x — xk), (1.2.24)

atx
l()i k=1

válida siempre que la función a(x) tenga derivadas continuas no cero en vecindades de
los puntos xk •

Ahora, si a(x) es estrictamente creciente, entonces

8 (x — /3 (y)) 5(z — 0(0)

8(a(x) y) = c4/3 (Y»
= al(x)

1111
de donde a' (x)8(a(x) — y) = 8(x — )3(y)). Diferenciando esta igualdad con respecto a
y, obtenemos
a a
ai(x) — 8(a(x) — y) =6(a(y) — x).
ay
1.2. Distribuciones Definidas sobre la Recta Real 27

Por otro lado,

a
8(a(x) — y) = (x)— 6(a(x) — y).
ax ay
Combinando estas dos últimas expresiones, obtenemos la relación
a a
8(a(x) — y) + —5(3(y) — x) = O.
Ox

Nótese que las variables x y y juegan el mismo papel y al aplicarse a las funciones
de prueba, producen

—
dy sb ( N y)) = 01(Y)-
43 0(3)

dx 95(a (x )) = á (x) da 0(a),

la regla de la cadena usual. Estas expresiones junto con la multiplicación por una
función suave implican la regla de la cadena para distribuciones,

d---T(a(x)) = da T(a(x)) á (x). (1.2.25)

Nota: También, podemos utilizar la fórmula (1.2.23) para establecer este resultado.
De la definición (1.2.12) y la fórmula (1.2.23), tenemos que

f {f (a(x))}5(x) dx — f f (a(x))41(x) dx = — f f (y) o' (e) dl;] dy

f(E)<y

= ff t ( y) [f 01 (e) de] dy.

0,(1)<Y

Como a(C) es una función creciente de y, entonces tenemos que a(C0) y. Por con-
siguiente,

Le y k
01(e) de = feo 951 (e) de = 95(e0)


28 CAPÍTULO 1. ELEMENTOS DE LA TEORÍA DE DISTRIBUCIONES

Por otro lado, si ponemos n = a(C), entonces

d
L dy
ukv d
(e)0(e) de = dy- f Yeo 0(1) cin = 0(e.),

así que

O i ( e) de = fc,(1)<, a'(e)0(1) de
fa (f)<v dy
Y

f {f (a(x))} 1 q5(x) dx =f f (y)[ciy- Lekv ai (e)0(e) del dy

f fi (a(x)) á (x)0(x) dx,

lo que demuestra (1.2.25).

Ejemplo 10: Supongamos que f (x) = h(x), la función de Heaviside unitaria y

a(x) = x 2 — a2 , entonces

2 — a2 ) 5(x 2 — a2) 2x = — [5(x — a) + 5(x + a)]

dx h(x a

= 8(x — a) — S(x + a).

Observación: La función considerada tiene un salto de magnitud -1 en x = —a y

otro de magnitud +1 en x = a.

h(x2 — a2)

—a a
x

Por lo que su derivada es —8(x + a) + S(x — a).

1.2. Distribuciones Definidas sobre la Recta Real 29

Una combinación de traslación e integración da como resultado otra operación

lineal importante sobre las distribuciones: la convolución de una distribución T con
una función 4' E D o en general con otra distribución S E DI . A continuación vamos a
dar la definición de convolución de una manera muy general, solo por mencionar este
tipo de operación.

La convolución T *O de una distribución T con una función suave 0, es una función

regular C°° dada por la fórmula

(T * 0)(x) = Tt [0(x — t)] = f f (t)0(x — t) dt (1.2.26)

definida puntualmente para cada x por separado. En particular, (T * 0)(0) = Tr(11,

donde q!)(t) = 0(—t). En especial á * 4' = 0. Es decir, la delta se comporta como el
elemento unidad para la "multiplicación" convolución.

También, es posible convolucionar dos distribuciones, pero con un enfoque más

débil. Si S es una distribución con soporte compacto y T es una distribución arbitraria
en D', entonces su convolución T * S es una distribución en Di que actúa sobre las
funciones de prueba E To de la siguiente manera:

(T * S)[0] = Tx[Sy[0(x + y)]] (1.2.27)

o equivalentemente

(T * S)[q5] = (T * (S * 0)(0). (1.2.28)

De nuevo, 5* S = S. Además, se puede diferenciar la convolución de dos distribu-

ciones, obteniendo

(S *T) (k) = S(k) *T = T(k) *S. (1.2.29)

Si las dos distribuciones involucradas en la convolución son de soporte compacto,

entonces la operación además de ser conmutativa es asociativa, lo cual da al espacio
de las distribuciones D' "cierta" estructura algebraica. En general, la convolución de
distribuciones no es asociativa, como se puede observar con el caso 1 * (5' * 1i).


30 CAPÍTULO 1. ELEMENTOS DE LA TEORÍA DE DISTRIBUCIONES

1.3 Distribuciones Multidimensionales

Por analogía con las distribuciones en IR, las distribuciones sobre lr se definen
como funcionales lineales y continuos sobre el espacio 79 = (lRn ) de las funciones
de prueba 0(x) 0(x l , x2 , . , x Ti ) infinitamente diferenciables (existen las derivadas
parciales de cualquier orden) con soporte compacto en Ift n , donde el soporte de una
función con argumento vectorial se define igual que en (1.2.4), con la diferencia que
ahora la función se anula fuera de un subconjunto n dimensional.

Ejemplo 11: De nuevo, al igual que en la Definición 1.1, el ejemplo inmediato de

función de prueba es
a2
a2 ) 114 < a,
1/2
0(x) =
Il x =
ll X1)
=1
o 114 a.

Es fácil verificar que D tiene la estructura de espacio vectorial y por consiguiente,

(X — X0)
si 01 y 02 E V, también c1 O1 + c202 E D y si 0 E D, también q5 E D, la
E
cual se anula fuera de una bola de radio e centrada en x° . Además, si 0(x) E D
y a(x) E C ce (Ir , entonces a(x)0(x) E D y si q5( xi, x2, , x,,,) E V' (IR') y
x„„,±2 , , xn ) E D (IR') , entonces q51P E D (Iftn).

Antes de continuar, es preciso introducir una notación especial que servirá para
compactar algunas expresiones utilizadas posteriormente. Sean k 1 , k2 , , k, enteros
positivos de tal forma que k = ( k 1 , k 2 , . , k,) representa a un multi-índice n dimen-
sional, Ikl Eig k1 + k 2 + • • + k r, y
aikj
k
D = ,g,k2
(1.3.1)
. axk.
axl "2 n

ki = O significa que la diferenciación con respecto a la variable correspondiente j no

está presente.
1.3. Distribuciones Multidimensionales 31

En el mismo contexto, cualquier operador diferencial lineal L de orden r en n

variables puede expresarse en la forma

L ak(x)Dk, (1.3.2)
Iklcr
donde los coeficientes ak (x) k2 (Xi, x2, xr,) son funciones arbitrarias. Por
ejemplo, si el orden es 2 y el número de variables es también dos, entonces

L= E ak(x)Dk ak(x)Dk + E ak (x)Dk + E ak(x)Dk

iki<2 ikko iki=2

, a ,a
ao,o(x, y) + ai,o(x, y )— + ao ' 1. (x, yj—
ax Oy
8 2 y
+a2,o(x, y) + (x, y) + ao,2(x, y) 9y2
axay

Definición 1.5 La sucesión de funciones de prueba {0„,(x)}'

m=i se dice que es una
sucesión nula en 13 si y solo si

existe una región común acotada fuera de la cual toda 0„,(x) se anula. Es decir,
los soportes de todos los miembros de la sucesión estan contenidos en una misma
bola simple, lo suficientemente grande.

para todo multi-índice k de dimensión n,

lim max D k tfini (x) =0.

,n-->co xElFr

Esto es, {0,,(x)} tiende uniformemente a cero en Ir, lo mismo para la sucesión
{Dk 0,(x)}

Por lo tanto, decir que {e„,(x)} es una sucesión nula en 1) significa que la conver-
gencia a O es muy fuerte (se requiere que O m y todas sus derivadas parciales tiendan a
cero uniformemente en ir).
32 CAPÍTULO 1. ELEMENTOS DE LA TEORÍA DE DISTRIBUCIONES

Definición 1.6 Un funcional lineal sobre D es continuo si siempre que {95,.„(x)} sea
una sucesión nula en D, la sucesión numérica
7;[0] = f f (x).0„,(x) dx
tiende a O cuando m oo. La integral es múltiple con respecto al volumen usual. En
este caso, el funcional lineal y continuo
f
T[0] = fr f (x)0(x) dx, (1.3.3)
recibe el nombre de distribución n-dimensional. El conjunto de distribuciones n-
dimensionales se denota por D i el cual tiene la estructura de espacio vectorial. Se
podría decir mucho más acerca del espacio de las funciones de prueba D y espacio de
las distribuciones D'" , incluso se podría hablar de la topología presente en estos espa-
cio y establecer el análisis de los mismos, pero además de desviar el sentido aplicado
del presente trabajo, esto representaría bastante espacio de discusión. Por esta razón
es que se introducen únicamente las distribuciones n-dimensionales más importantes,
algunas de sus principales propiedades y conceptos de cálculo, sin restar la importancia
del aspecto riguroso que poseé la teoría de distribuciones.

Teorema 1.3.1 Si f : IR" -Y IR es localmente integrable (1.1 f dx existe, en el

sentido de Lebesgue, para todo dominio acotado fi c IR"), entonces f genera una
distribución n-dimensional mediante la definición del producto interior usual

T f [O] (f, fur f (x )0( x) dx. (1.3.4)

Dem. La linealidad de (f, 0) sobre D sigue directamente de la definición (1.3.4). Para

probar la continuidad, sea {0„,(x)} una sucesión nula tal que todos sus elementos se
anulan fuera de la bola finita n. Entonces,
0„,)¡ < (max 10,7,(x)i) f I (x)i dx.
€11

Como {0„,} es una sucesión nula, implica que hm max 14,1 = O. La integrabilidad
971-› CO

local de f garantiza que f dx es finito, así que hm ( f , 0„,) = O. Por lo tanto, el

M-4 CO

funcional (1.3.4) es continuo y por consiguiente una distribución.

1.3. Distribuciones Multidimensionales 33

Observación: i) Por medio de (1.3.4) toda función f localmente integrable puede

considerarse también como una distribución. Como función puntual, f tiene el valor
f(x) en cada punto x E IR"; como distribución, f tiene una acción (promedio) (LO)
sobre una función de prueba elegida O E D. Además, el valor de una función en
cualquier punto puede recuperarse (de forma única) de los valores del funcional sobre
las funciones de prueba que se elijan.
ü) Si A (x) y f2 (x) coinciden excepto en un número finito de puntos, entonces
generan la misma distribución. En general, dos funciones se dice que son iguales
casi dondequiera si fn lh -121 dx O para todo conjunto acotado SZ. Dos funciones
localmente integrables iguales casi dondequiera generan la misma distribución. En caso
contrario, generan distribuciones diferentes. Una vez que se muestra esta afirmación, se
considera a funciones iguales casi dondequiera como la misma distribución. Para hacer
esto más formal y consistente es necesario introducir la noción de clases de equivalencia,
no necesario debido al objetivo de este trabajo.

Definición 1.'T Una distribución se llama regular si es generada por un kernel f (x)
localmente integrable, mediante la expresión (1.3.4).

Ejemplo 12: a) (c, = I c g5(x) dx (distribución regular generada por la función

constante f (x) = e).

b) Sea SI un dominio en IR". El funcional

(hl], 45)
f <15(x) dx (1.3.5)
a
es una distribución regular, llamada distribución de Heaviside general. La función hn
es la función de Heavíside, también llamada la función indicador del conjunto

def 1 para x E fi,

hn = (1.3.6)
Sl O en otra parte.

e) Cualquier función continua o continua por piezas define, por medio de (1.3.4),
una distribución regular


34 CAPÍTULO 1. ELEMENTOS DE LA TEORÍA DE DISTRIBUCIONES

Idas distribuciones que no son regulares se llaman distribuciones singulares. Los

ejemplos más comunes de tal tipo de distribuciones son las deltas de Dirac, sus derivadas
y los valores principales de Cauchy que se desprenden de integrales cuyos integrandos
tienen singularidades simples.

1.3.1 La Delta de Dirac Multidimensional

La más importante de las distribuciones singulares, desde el punto de vista del

presente trabajo, es la delta de Dirac, ya que para determinar una función de Green
(Capítulo 2) es de mucha utilidad conocer sus propiedades, además de otras aplicaciones
interesantes en las que está presente También, sin restar importancia se introducen
otras distribuciones singulares, con el fin de establecer que el marco teórico de las
funciones elementales se extiende de manera considerable.

La delta de Dirac multidimensional es el funcional lineal y continuo del tipo (1.3.3)

que asigna a la función de prueba su valor en un punto específico E IR", el cual es
llamado polo o singularidad de la delta,

5(x - 4)kbl r- 0(4) (1.3.7)

y tiene la propiedad que f 5(x - 4) dx = 1, si E fi

Otro ejemplo de distribución singular es el dipolo multidimensional, el cual vamos

denotar por
T(x - d-41 • V8(x 4), (1.3.8)

donde C es el vector unitario en la dirección del dipolo, y el operador derivada direccional

V actúa sobre 5 de tal manera que T(x - CO) =-- .e • Vc/iix=k Por ejemplo,
en IR3 la configuración de una fuente electrostática que consiste de una carga positiva
1
- localizada en el punto + (1) £, y una carga negativa de la misma magnitud en
e
1.3. Distribuciones Multidimensionales 35

La distribución correspondiente a estas cargas tiene la acciót

95 e2 el. (1 t)

sobre cualquier E D. El dipolo unitario con eje P se obtiene cuando e -+ O (algo

así como que se juntan las cargas a la vez, incrementando la magnitud de ambas
indefinidamente). La función generalizada correspondiente
dcfr
1- [0] = di (I),
derivada distribucional de la delta, define una distribución singular, que se identifica
con el dipolo.

En el caso multidimensional, 8(x) también se puede obtener mediante limites

débiles de distribuciones representadas por kerneles regulares fk : IR" 1R. Por
ejemplo, en IR3 es conveniente tomar fk (x) = g k ( x i ) gk (x2 ) g k (x3 ), con gk función regular
de una variable aproximadora de la delta. En este contexto,

8(x) = 5(x i )5(x 2 )6(x3 )• (1.3.9)

Esta separación como un producto se conoce como propiedad de isotropia de la delta
3—dimensional (o invariancia respecto a el grupo de rotaciones de IR3 ). Una prueba
de esta propiedad es considerar la sucesión de funciones Gaussianas, Ejemplo 3a) y su
producto coordenada a coordenada, lo cual lleva a

fe (x) l E ) 3 er2/2e2,
Grir

que depende únicamente de la magnitud del vector, r = y no de su orientación

en el espacio. Lo mismo sucede con otro tipo de aproximadores débiles. Similarmente,
podemos considerar a la delta de Dirac en el espacio—tiempo como el producto
8(x,t) = 8(x)8(t).

Ahora, si se considera el cambio de coordenadas 1-1, y = a(x), entonces

5(a(x) — e ) = 8 (,3 (1) — z) (1.3.10)
1J (x)I
donde x = /3(y) son las coordenadas inversas a a(x) y J es el Jacobiano de la
transformación de las coordenadas y a las coordenadas x.


36 CAPÍTULO 1. ELEMENTOS DE LA TEORÍA DE DISTRIBUCIONES

Ejemplo 13: En particular, es interesante mostrar que la delta de Dirac en coorde-

í esféricas y cilíndricas tiene la expresión
nadas

8(u) = 6( 9) (P, 45,

47rp2 COORDENADAS ESFÉRICAS

(x , y , z) 6(x) = 8(x)6(y)8(z)
COORDENADAS CARTESIANAS

6(ii) = 8(r)6(z) (r , , z)
2irr COORDENADAS CHADRICAS

i 1 Únicamente vamos a considerar el caso del cambio de coordenadas rectangulares

kir
J1 x, y, z a cilíndricas r, 9, z. Tenemos que para este caso, la ecuación (1.3.10) tiene la
' 11 11
I expresión
I
' ' 8(r — r 0)8(61 — 00)6(z — zo)
1,I 6(x — E) = ro > 0. (1.3.11)

De aquí podemos observar que, debido a que el sistema cilíndrico es singular en el

origen, no se puede considerar directamente la situación cuando la fuente está localizada
en el origen. Sin embargo, vamos a considerar el problema mediante un proceso de
límite,

considerando un anillo uniforme de fuentes de magnitud total unitaria, a una profun-

didad h (r = ro , z = —h). Como la porción de anillo entre Bo y eo + M0 representa
AB,
a un fuente de magnitud —, su representación distribucional está dada por (1.3.11)
27r
multiplicada por A00/27r. Así, la representación de todo el anillo se obtiene integrando
desde O, = O hasta 00 = 27r, dando como resultado

6(r — ro)8(z + h)
2irr
1.3. Distribuciones Multidimensionales 37

41

Tomando el límite cuando ro —> O, tenemos que para una fuente puntual unitaria
localizada a una profundidad h tiene la representación

5(x + h) = 5(r)5(z +
2irr
donde h = (O, O, —h), de tal forma que si la fuente se localiza en el origen, entonces
tiene la representación

5(r)8(z)
5(x) —
2irr
o

Por otro lado, la discusión en torno a la delta de Dirac de hecho puede extenderse
a distribuciones delta cuyo soporte singular sea más general que un simple punto. Por
ejemplo, si a es una superficie en IR3 y p(x) es una función localmente integrable
sobre a (densidad de área, que por simplicidad se toma p(x) -a- 1), entonces la delta de
Dirac superficie (o capa superficial simple) 6, -a- 8(c) se establece como

&Y] k 5(c)p(x)q5(x) dx = I p(x)0(x) dS, E 0113 ) . (1.3.12)

Un caso especial ocurre cuando es la esfera r Y P(x) — 47rr2:

1 0(x) dS.
lukbl ilirr2
De la misma forma, se establece que la delta de Dirac línea 51 -a 5(1) para una curva
en IR? es
54 01 d" 5(2) 0(x) dx = f 0(x) ds. (1.3.13)

En la Subsección 1.3.3 consideraremos con más detalle la deducción e interpretación

de estas expresiones de la delta para este tipo de soportes, las cuales son importantes
porque se utilizan en física en muchos contextos. Por ejemplo, se utilizan como mode-
lo matemático para superficies y alambres eléctricamente cargados en electrostática,
fuerzas concentradas y momentos en mecánica de sólidos, vórtices y fuentes en mecánica
de fluídos, entre otras.
 38 CAPÍTULO 1. ELEMENTOS DE LA TEORÍA DE DISTRIBUCIONES

1.3.2 Operaciones Generales

Ahora, si a(x) y f (x) son localmente integrables, el producto a(x) f (x) no
necesariamente es localmente integrable (aún cuando a f sea localmente integrable, su
acción sobre 95 puede no estar relacionada con la acción individual de a y f). Sin
embargo, si a(x) E C' (IR"), entonces aq5 E D y

aT[d.] (a f , f a(x) f (x)qb(x) dx = ( f , ad.) T[a(b].

tan
Así, a(x) f (x) genera una nueva distribución, la multiplicación de una distribución
por una función infinitamente diferenciable. Un caso específico es

a(x)6(x — = a(4)5(x — 4). (1.3.14)

En general, bajo las condiciones apropiadas, es válida la fórmula de cambio de

variable multidimensional para distribuciones

f f ( y - 1)0( y) dy f f (a (x) — C)0( a (x ))P ( x )i dx. (1.3.15)

Generalizando el concepto de derivada, tenemos que

8 _ [50]
Tkb] — —T (1.3.16)
axi axi

ti
y aplicando repetidas veces (1.3.16),

Dk T[0] = (-1) 1k1 [Dk d.] . (1.3.17)

Si L es un operador diferencial lineal de orden r en n variables (ver (1.3.2)), con

;J,
k
coeficientes a infinitamente diferenciables, entonces

L T[¢] (L f , kk
( E a D f, q5) = (f, E (-1)i k iDk (a05))
iki<r Ik15r
(1.3.18)
= E
T [iki<r
(-4)1kLok (ak
1.3. Distribuciones Multidimensionales 39

Nota: si) El operador de orden r que aparece en el lado derecho de (1.3.18) es conocido
como el operador formal adjunto de L y se denota por L. Así,

L T[0] = T[L*q5]. (1.3.19)

ii) Siempre se tiene que (L*) * = L y si L = L*, entonces L se llama operador

formal auto-adjunto.

Por último, consideremos ahora una función f (x) con discontinuidades multidi-
mensionales. En específico, sobre la superficie a(x) = O, con discontinuidad de salto
[f] dado por la relación

Lfl = f (0- = 0± ) - f = (1.3.20)

Nota: a = 0 + representa al lado derecho de la superficie (hacia donde apunta Va).

f
Al igual que en la ecuación (1.2.20), nos gustaría determinar . Para hacer esto es
oxi
necesario poner sobre la superficie un sistema de coordenadas (u 1 , u 2 ) y extender el
sistema al espacio en una vecindad de la superficie a lo largo de sus normales. Sea
u 3 = a la tercer coordenada, la cual está bien definida en la vecindad de a = 0.
Debemos señalar que f es continua en las variables u 1 y u 2 , pero discontinua en la
tercer variable u3 . Por lo tanto, tenemos que

bf _
i-=- 1, 2
pu, au,
(1.3.21)
f
bu
82,13 = — + 115(u3).

Así, por la regla de la cadena, tenemos

Buss

+ Lf115.8(a)« (1.3.22)
au; ari au;afax,
afaf
ax, aui + Lfl 1-
°ax,'6
3 (u3) =

En forma vectorial tenemos que

= Vf Ifl Va5(a). (1.3.23)



40 CAPÍTULO 1. ELEMENTOS DE LA TEORÍA DE DISTRIBUCIONES

Análogamente, podemos definir la divergencia y el rotacional generalizados como

V•f=V•f+Va• [fi 8(a),

(1.3.24)
V xf=Vxf+Vcrx [fi 8(a).

La derivación rigurosa de ambos resultados requiere cierto manejo de la invariancia

de la divergencia y rotacional en sistemas de coordenadas curvilineas en general, parte
importante del análisis vectorial--tensorial clásico. Los resultados (1.3.23) y (1.3.24) se
pueden compactar, utilizando * para las tres operaciones, en la fórmula
.111
4111 V *f=Gr *f+ V a* [fi 8(a). (1.3.25)

Lo más importante de esta parte es que la ecuación (1.3.17) implica la excelente

y definitiva conclusión acerca de la diferenciabilidad.

Teorema 1.3.2 Toda distribución es diferenciable a orden arbitrario.

1.3.3 Interpretación de 8(a)

El interés principal que nos mueve en esta parte es referente al significado de 8(a)
y 5'(a), donde a = O es una superficie en el espacio 3-dimensional, bajo la suposisión
que 'Tall = 1 en todo punto sobre la superficie. Si a no tiene esta propiedad, entonces
al = a 111Voll si la tiene y se redefine la superficie. Con estas condiciones, consideremos
de nuevo la integral (1.3.12),

5,401= f 8(a)0(x) dx (1.3.26)

introduciendo para su análisis un sistema de coordenadas curvilineas (u 1 , u2 ) sobre la

superficie a = 0, extendiendo localmente estas variables al espacio a lo largo de las
normales a la superficie. La tercer coordenada, u 3 = a, representa la distancia local
desde la superficie ya que =—da = 1. Así, a = u3 = const O es una superficie
du3
paralela a a = O, con u3 pequeña.
1.3. Distribuciones Multidimensionales 41

•
De la Geometría Diferencial (ver por ejemplo [7]), se tiene que si x x(u l , u2,u3)
es uno a uno sobre la superficie, entonces

dx = U2, U3) d'U/ dU2dU3 (1.3.27)

donde g2 es el determinante de la primera forma fundamental correspondiente a la
superficie a = u3 = const. Utilizando (1.3.27) en (1.3.26) e integrando respecto a u3,
obtenemos

(5, [0] f (fi[ x ( u i, u2, u 3)1 8 ( u 3)1/42( u 2, u2, u 3 ) duiclu2du3

= (fi[ x ( 24., u 2, 0 )] 1/42(4 , u2, O) duidu2 (1.3.28)

= f 0(x) dS

Ahora, respecto a 8'(c), tenemos

IVO] = f 0(x)6I (u) dx

= f OIx(u2, u2, u3)]8 1 (u3) 3g2(u2, U2, 2i3) dUidt2dU3 (1.3.29)

a
= f
au3 3g2 u i , U2, 1£31
[499(x ) y
1,3=o
duidu2.

De nuevo, utilizando otro resultado de Geometría Diferencial, tenemos que

fi
n 1/42N., u2, u3) = — 2K„.(u2 , u2, u3) 1/42(24 , u2, u3),
ou3
donde Ka denota a la curvatura media de la superficie u = u3 = const. Así,
ao
(5a , [0] = —
utt3
( ( X (U1, U2, 0 )) 1/42 ( UD U2, O) dUldU2

f21<a(u1, uz, 0 ) 1b [x ( u 1, u2, 0 )] 1/42 (U 1 / U2/ 0) dUldu2 (1.3.30)

= f {-- + 2Ka (x)0(x)] dS,

fin
an
es la derivada normal usual de 0. Esto muestra la importancia de la Geometría
8n
Diferencial en la teoría de distribuciones multidimensionales y que estas no son fáciles
de tratar.
 42 CAPÍTULO 1. ELEMENTOS DE LA TEORÍA DE DISTRIBUCIONES

4. 1.3.4 Producto de Deltas

En general, en el espacio de distribuciones la multiplicación no está definida.

Sin embargo, este hecho debe de esclarecerse puesto que por un lado, obviamente las
funciones ordinarias, las cuales también son funciones generalizadas, siempre pueden
multiplicarse y por consiguiente también deben poder multiplicarse en el sentido dis-
tribucional. Pero, puede ser que la multiplicación de una distribución regular por si
misma lleve a a cierta ambigüedad, como por ejemplo, en el caso de multiplicar la
función de Heaviside por si misma n veces, h n (x) = h(x)h(x) • • • h(x) h(x). Tomar
la derivada en ambos lados lleva a nhn-1 (x)hi (x) = h'(x), de donde supuestamente
se obtiene que hn-1 (x) = 1/n lo cual es una contradicción ya que si se toma n = 2,
entonces h(x) = 1/2. La multiplicación de una distribución regular y una no singular
o dos distribuciones singulares puede que no siempre este definida. Por ejemplo, 82 no
tiene ningún significado en el sentido de distribuciones, así como tampoco f (x)5(x - e)
cuando f(x) tiene discontinuidad de salto en x = e. Sin embargo, es posible que este
problema de la multiplicación se pueda remediar considerando otros espacios de dis-
tribuciones, pero ya esto es material que no vamos a abordar en el presente trabajo.
Únicamente vamos a considerar un caso muy especial; dos deltas de Dirac multidimen-
sionales se pueden multiplicar siempre que sus soportes sean apropiados.

Sean al = O y a2 O dos superficies que se intersectan a lo largo de una curva

1' (ver Fig.11)
752 = n2

Dal =

Cr2 = O

0-1 = O
Fig.11 Soporte del Producto de dos Deltas.
1.3. Distribuciones Multidimensionales 43

Supongamos que Do-1 = n 1 y Val = n2, donde lin]. = = 1. El asunto

ahora es interpretar

óol 6* 2[0] =f 6(01)8(52)q5(x) dx. (1.3.31)

Sobre el plano normal local a la 1"-curva definimos u 1 = al, uz = az y u 3 = s, donde
s es la distancia a lo largo de la r-curva. Extendiendo u 1 , u2 al espacio en una vecindad
del plano a lo largo de una de sus normales, tenemos
duidu2dus
dx =' (1.3.32)
sen
donde sen O = x n2 11. Sustituyendo (1.3.32) en (1.3.31) e integrando con respecto
a u1 y u2 , tenemos que

5,,(5,[0]
= I 1)x)
(
se n
8(ui )5(u2) duidu2du3 = (fr( x)
o1=0se n O
ds. (1.3.33)

a =0

Este resultado podría resultar muy útil en la aplicación, así como otros más generales
derivados de multiplicar distribuciones Delta adecuadas.

1.3.5 Parte Finita de una Integral Divergente

La parte finita de una integral divergente es sin duda un punto clave para muchos
problemas aplicados. Previamente en la Sección 1.2 se dió la definición y una breve
introducción de este tópico con funciones de una variable, denominando al resultado
"valor principal de Cauchy". Por ejemplo, mencionábamos que es interesante mostrar
que

PV ( 1)
x= dx
lnixi

y precisamente aprovecharemos el espacio de esta parte para obtener esta relación.

La función f (x) = lnixi es localmente integrable, mientras que su derivada ordinaria
= 1 no lo es sobre cualquier intervalo que contenga x = O. Sin embargo,
sabemos que como una función generalizada, lnixi tiene derivada de todos los ordenes.
Pero, ¿cuál es entonces la relación entre la derivada generalizada y la derivada ordinaria
de f (x)?
44 CAPÍTULO 1. ELEMENTOS DE LA TEORÍA DE DISTRIBUCIONES

Trabajemos con el funcional Tf [0], que representa a f (x), como sigue

T[q]= f lnix1q5(x) dx.

Utilizando la definición de derivada generalizada (1.2.12), tenemos

711 [0] — f 1M01 (x) dx.

La integración por partes sería una opción para obtener el término 1/x en est
representación integral, pero no se puede efectuar debido a que 1/x no es localmen-
integrable. Otra opción es introducir otro funcional auxiliar, dependiendo de E del cu
71q5] se pueda obtener como límite Sea he (x) la función ayuda

—e < X < as
hs(x) = {°1
otra parte

definida para alguna constante a > O y un parametro e > O (ver Fig.12)

he(x)
1

-E ae

Fig.12 Función de ayuda para obtener la parte finita de la integral divergente f q5(x)/x

Es obvio que 1Mx¡ puede escribirse como el límite de funciones

lim he (r)lnIxj = lnixi.

-30

Si definimos
T:(4)) — f he (x)lnixl 61 (x) dx,

entonces por el Teorema de Completes de Di (Teorema 1.2.1) tenemos que

lim n [0] = r [0].

e--10


1.3. Distribuciones Multidimensionales 45

Ira función h,(x)lnix I tiene dos discontinuidades de salto en x = — e y x = as.

Podemos usar la integración por partes clásica, partiendo la recta real en dos intervalos
o usar la derivada generalizada (mucho más conveniente)

[he (x)lnix 0(x) dx,

f dx
integrando sobre Sopt 0, sin tener que preocuparnos de los límites en la integración
por partes puesto que 0 se anula en estos puntos límite.
Tomando la derivada del término entre parentesis cuadrados en la expresión ante-
rior y utilizando a(x)5(x — e) = a(e)5(x — C) (Ejemplo 5), tenemos que

d he(x)
— [he (x)lnixl ] = ln s5(x + e) + ln (as)5(x — as).
dx

y por consiguiente, tenemos que

TE [0] = —ln eq5(—e) + ln (ae)q5(ae)+ f he (x) 0(x) dx

x

= 0(0)1n a + f 116(x) 0(x) dx + o(e),

donde o(s) representa a los términos en e de orden grande. Tomando el límite, cuando
E —Y O, obtenemos

T'ko] = lime,0 T;[0]

= 0(0)1n a + lim f he(x) 0(x) dx

e—>0 x

= 0(0)1n a + lim [f °(xx) dx + f °( x) dx]

e—>0
—oo x

Para a = 1, ln a = O y el límite que queda precisamente es lo que se define como

valor principal de Cauchy, relación (1.2.5). Esto es,

T'[0] = PV f 95(x) dx.



46 CAPÍTULO 1. ELEMENTOS DE LA TEORÍA DE DISTRIBUCIONES

Tontar a .==_ 1 en el procedimiento anterior es por conveniencia, pero no es necesario

ya que numéricamente el resultado es el mismo para cualquier valor de a. En general,
a este procedimiento de límite se le conoce como "tomar la parte finita de una integral
divergente o regularizar una integral con singularidades" (regularización canónica de
Gelfand, [6]). Así,
1
— lnixj = -
dx x
sobre cualquier intervalo abierto que no incluya a x = O, pero cuando x = O pertenece
al intervalo considerado, entonces
pv ( 1 \ = d
x) II
Ahora, ¿cómo se utiliza este procedimiento en la aplicación? Supongamos, por
ejemplo, que hemos reducido la solución de un problema a la evaluación de la expresión

u(x) = —
dx 0(y)lnlx – yl dy, x El.

Si la integral es continua como función de x, entonces

d
u(x) = f 0(y)— — y dy = f q5(y)PV ( 1 dy.
dx x — y)
la cual se interpreta como la parte finita de la integral divergente por medio del pro-
cedimiento descrito anteriormente.

Nota: Si no se utilizara el valor principal de Cauchy, entonces aplicando integración

por partes implicaría que
u(x) = – f q5 1 (y)lnix — y dy

también sería un resultado legítimo en el caso de existir la integral. Sin embargo, el

problema es que aún cuando esta ecuación integral se desprendiera de algún problema
aplicado la función 0(x) podría tener singularidades en la frontera del intervalo I, de
tal forma que el procedimiento de integración por partes no sería válido y la ecuación
integral sería divergente. Por lo tanto, la única alternativa que nos queda es introducir
el proceso de regularización mediante el valor principal de Cauchy de 1/(x – y), el cual
es un kernel bien conocido en la teoría de las ecuaciones integrales singulares, y ahora
lo podemos identificar como una distribución singular


1.3. Distribuciones Multidimensionales 47

El concepto equivalente para dimensiones mayores es de difícil manipulación, de-

pende del problema y puede tener serias implicaciones. Por ejemplo, en tres dimen-
siones, consideremos la integral
52
I(x) (1.3.34)
axl

donde r = II x — y11, y x E 9 (alguna región en el espacio). Suponiendo que la integral

82 8-2
es una función C i. en x, podemos reemplazar — por — e intercambiar integración
8x21 axl.
y derivada,

(1 a2 1
I(x) 82
0(y)— - ) dy = f 0(y)— dy. (1.3.35)
ax? r n ay? r

Nota: i) Utilizamos diferenciación generalizada en lugar de diferenciación ordinaria

porque el resultado es una integral divergente.

Tenemos que

a (1.} = 7.1
=- X1 - Y1
7') 74'

a2 (1 ) 3r? — r2
ayi r r5

iii) Como rj /r 3 es localmente integrable, escribimos

2 (1) ra (ri
(1.3.36)
a y? =
Con estas observaciones, procedemos a encontrar la parte finita de la integral
divergente (1.3.35). Sea
f (x, y, e) = g ( r i, r 2, r3) — e = O (1.3.37)

una superficie suave por pedazos que contiene al punto y = x, donde ri = Xi - yi,
i = 1, 2, 3 y g una función homogénea de orden 1; g(ar i , are, are) = a g ( r i, r z, rs),
condición que asegura que la superficie g(ari , are, are) — e = O corresponde a


48 CAPÍTULO 1 ELEMENTOS DE LA TEORÍA DE DISTRIBUCIONES

g(ri , r2 , r3 ) — e/a = O para a 0. Así, las superficies g e = O para diferentes

valores de a son similares en forma. Además, de la homogeneidad de g se sigue que
f (x, y, O) = g(ri , r2 , r3 ) = O consiste de un solo punto (y = x). Por ejemplo, para una
esfera con centro en y = x y radio e, tenemos que

f (x, y, E) = Vr? + r? + rg — e = 0.

Ahora, sea Vy f = n, donde rt es el vector normal unitario a la superficie. Además,

suponiendo que f > O fuera y f < 0 dentro de esta superficie, respectivamente, intro-
ducimos una función auxiliar He (y) mediante la relación

HE(y) = 1 (f (1.3.38)
O < O).

La derivada generalizada requerida en (1.3.36) se define por la relación

lim a 114(Y)ril
6—>0 r3
(1.3.39)
[ri ni 5(f) + 371 —
5 r2nÁy)}
6—+0 r3

donde n 1 es la componente de n en la dirección del eje yl . Por consiguiente, 1(x) puede

escribirse como

— r2
1(x) = fi m c/5(y) dS + lim f He(Y)lb (Y) dy, (1.3.40)
s —r0 f1=
0 e—>0 r5

donde hemos usado (1.3.28) para integrar 8(f) en (1.3.40). Utilizando una expansión
en series de Taylor para 0(y) en y = x, encontramos que

lim If=0 rir31 0(y) dS a f0(x),

e—a)

donde afe es una constante que depende de la forma de la superficie f = 0. Por

47r
ejemplo, para la esfera dada por ir? + + = e, tenemos que af = — mientras
3 '
1.3. Distribuciones Multidimensionales 49

que si f O es un cilindro circular con eje paralelo al eje y l , radio de la base e y altura
-y, eh > 1, entonces af 4/r. Así, la ecuación (1.3.40) puede finalmente escribirse
como

1(x) = a/ 0(x) + lim f 3r1 r2 14(y)0(y) dy. (1.3.41)

e—ko r

Numéricamente, I(x) debe ser la misma sin importar la forma de la superficie

f = 0. Como (37-? —r 2 )/r5 cerca de y = x toma valores tanto positivos como negativos,
la forma de f = O cuando e --> 0 afecta al valor de la integral de volumen en el
proceso de suma (efecto similar al del conocido resultado para series condicionalmente
convergentes, las cuales pueden hacerse convergentes a cualquier valor reordenando los
términos). El término a f 0(x) en la ecuación (1.3.41) compensa el cambio en el valor
de la integral cuando la superficie f = O cambia de tal forma que 1(x) numéricamente
es el mismo valor sea cual sea esta superficie. Precisamente en esto consiste el proceso
de regularización.

En la aplicación, la principal implicación es en la forma de obtener el valor de la

integral de volumen, ya que comúnmente en la práctica esta se obtiene numéricamente.
La región de integración en la integral de volumen tiene un hoyo en el que se encuentra
y = x, cuya frontera está dada por f = O. El valor de a f debe, por consiguiente,
corresponder a la malla utilizada en la integración de volumen numérica. Si el hoyo
es rectangular, caso más frecuente, entonces ninguno de los valores mencionados para
af es apropiado para el problema. Si -en lugar de tomar una derivada generalizada

dentro de la integral impropia (1.3.35) se utiliza una derivada ordinaria, entonces nos
lleva a una integral divergente y parecería que el problema físico del cual proviene esta
integral estaría mal o no tendría solución. Por lo tanto, en muchas ocasiones, no es
que un problema sea difícil o éste mal planteado lo que lo hace complicado, sino los
procedimientos no adecuados que se utilizan para atacarlo. Así, el analista debe revisar
siempre la causa de la apariencia errónea de un problema real y utilizar la técnica más
apropiada para resolverlo, y la teoría de las funciones generalizadas en algunos casos
se convierten en un buen candidato.
 50 CAPÍTULO 1 ELEMENTOS DE LA TEORÍA DE DISTRIBUCIONES
151

Comentario 1
En cualquier dimensión, tenemos la relación

licSce TYDS'DE',

V= {funciones infinitamente diferenciables con soporte compacto}

= espacio dual de D, Distribuciones de Schwartz.
= {funciones infinitamente diferenciables de rápido decrecimiento}
S' = espacio dual de S, Distribuciones Moderadas o Naturales (deT emple).
E = {funciones infinitamente diferenciables C"}
= espacio dual de E, Distribuciones con Soporte Compacto.

En nuestro contexto es más que suficiente considerar las distribuciones comunes, las
Distribuciones de Schwartz 73' puesto que este es el más general (el espacio de las dis-
tribuciones fuertes). De hecho la distribución que más se utiliza en el presente trabajo
es la delta de Dirac la cual es una distribución de soporte compacto, lo que implica que
pertenece a los tres espacios más importantes de distribuciones e', S' y D'. Por con-
siguiente, podemos trabajar en cualquiera de estos espacios de funciones generalizadas.
Sin embargo, no está demas señalar que el espacio E' es suficiente para el aspecto
algebraico formal, mientras que el espacio S' es necesario para el análisis riguroso,
en específico para el estudio de la Transformada de .7" ourier Generalizada, la cual es
herramienta de fundamental importancia para tratar con libertad un gran número de
problemas en Ecuaciones Diferenciales Parciales. Estudiar y analizar con detalle todo
el material lleva tiempo, espacio y esfuerzo, y por lo mismo, a manera de sugerencia es
posible desarrollar un nuevo tema de tesis paralelo al presente trabajo.
 Capítulo 2

Funciones de Green

A traves del tiempo el hombre ha tenido que modelar procesos naturales que ha
observado a su alrededor; física y matemáticamente. Una buena parte de las leyes
de la naturaleza se describen mediante ecuaciones diferenciales. Por ejemplo, algunas
de estas provienen de procesos de difusión, movimiento de partículas, flujo de fluidos,
ondas elásticas, acústicas, electromagnéticas, entre otras. Algunas ecuaciones que sir-
ven como modelo son la Ecuación de Laplace, Poisson y de Onda. Las ecuaciones de
Laplace y de Poisson son fundamentales en el estudio de mecánica de fluidos, campos
gravitacionales y electrostáticos, y también muy útiles en el análisis de películas de
jabón y cristales líquidos. Por su parte la ecuación de onda permite analizar todo tipo
de problemas de vibración. Aún cuando las soluciones de estas ecuaciones se consi-
deran en detalle hasta en el Capítulo 3, se mencionan en el presente capítulo porque el
objetivo es establecer una forma de como resolver problemas con condiciones de fron-
tera; el enfoque mediante las funciones de Green es un método de solución importante
para este tipo de problemas.
En la primer parte del capítulo se introduce una función de Green de una manera
motivadora al considerar el problema clásico de una cuerda vibrante. Posteriormente,
se define formalmente el concepto de función de Green, se establecen los pasos para
su construcción, se consideran las propiedades más importantes y se presentan algunos
ejemplos simples en dos y tres dimensiones.

51
52 CAPÍTULO 2. FUNCIONES DE GREEN

2.1 Introducción
Un problema con valores en la frontera (PVF) consiste de
Una ecuación Diferencial Lu f,
(2.1.1)
Condiciones de Frontera Bu = h,
donde el par f, h se conoce como los datos de entrada del problema y u la respuesta,
o en otras palabras la función desconocida o solución por determinar.
Ahora, ¿cómo proceder si se tiene que f representa a una fuente concentrada,
como por ejemplo una carga eléctrica puntual, una fuente de calor intensa o una fuerza
externa que impacta a una masa intempestivamente? ¿Se puede determinar alguna
solución que satisfaga tal condición? Si es así, ¿cómo se determina?
Para los casos descritos, se introduce la delta de Dirac en (2.1.1) como término no
homogéneo y al resolver el nuevo problema aparece el concepto matemático llamado
función de Green, la cual representa la respuesta inmediata o solución fundamental del
problema planteado y depende del espacio, las condiciones de frontera y su geometría.
En términos de las funciones de Green el PVF con datos arbitrarios se puede resolver
de una forma explícita en la que se muestra claramente la dependencia de la solución
de los datos involucrados, lo cual viene dado mediante un proceso de convolucion.

2.1.1 Motivación

Para motivar la introducción de las funciones de Green, primeramente considera-

remos el problema de determinar la curva de deflexión de una cuerda elástica sujeta a
cierta tensión producida por una carga. Esto es, consideraremos una cuerda elástica
flexible que es estirada hasta una longitud 2 bajo una tensión T como se muestra en la
siguiente figura

Fig.13 Cuerda estirada y deformada.

2.1. Introducción 53

Supongamos que la cuerda soporta una carga distribuida por unidad de longitud w(x)
de tal forma que bajo esta carga todas las deflexiones son perpendiculares a la longitud
de la cuerda y para cualesquier valor de x la carga que actúa sobre esa posición en la
cuerda es la misma antes y despues de que la cuerda se deforma. Consideremos un
elemento arbitrario de la cuerda estirada y deformada (Fig.14) por la acción de una
carga distribuida para describir las fuerzas que actúan sobre el sistema
F2

ay

w(x)A(x)
al

co(x) <
Fig.14 Fuerzas que actúan sobre la cuerda.
Cuando la cuerda deformada está en equilibrio tanto la fuerza horizontal como la fuerza
neta vertical sobre cada elemento deben ser cero,
F2 cos Ce2 = F1 cos a1,
(2.1.2)
F2 sen a2 — F1 sen al = —w(x)Ax.

La primera ecuación nos dice que suma de las componentes horizontales de la fuerza en
la cuerda es cero, y además como las deflexiones son tan pequeñas, esta componente
horizontal constante no difiere de manera apreciable de la tension T en la cuerda.
Entonces, dividiendo los términos de la primera ecuación de (2.1.2) entre las cantidades
iguales Fi. cos a1 , F2 cos a2 y T, tenemos
w(x)Ax
tan a2 = tan ai (2.1.3)
T
Como tan a2 es la pendiente de la curva de deflexión en el punto x + ¿x y tan al es la
pendiente de la curva de deflexion en el punto x, la ecuación (2.1.3) se puede escribir
como
y'(x + x) — y'(x) — —w(x)
(2.1.4)
Ax T
Cuando el límite Ax O, se obtiene la ecuación diferencial que describe a la curva
de deflexión de la cuerda,
Ty" = —w(x). (2.1.5)
54 CAPÍTULO 2. FUNCIONES DE GREEN

Ahora, la cuestión es determinar la deflexión de la cuerda bajo la influencia de una

carga concentrada de "magnitud" P en algún lugar de la cuerda. La representación
matemática de una carga concentrada es un concepto ficticio, ya que físicamente no es
posible; cualquier carga diferente de cero concentrada en un solo punto implicaría una
presión infinita que cortaría inmediatamente a la cuerda. Sin embargo, el uso de las
cargas concentradas al analizar los sistemas físicos, tal como en el caso de vigas, cuerdas
y cargas eléctricas es común y natural, y el objetivo es mostrar que estas se pueden
representar mediante la distribución delta de Dirac. Si en el sistema se considera una
carga concentrada en algún punto sobre la cuerda, entonces la ecuación (2.1.5) es cero
en todos los puntos de la cuerda excepto en el punto donde se localiza la carga. Esto
es, la ecuación no homogénea (2.1.5) se reduce a la ecuación homogénea, y" = O, cuya
solución es una función lineal, y por lo tanto, se concluye que la curva de deflexión de la
cuerda bajo la influencia de una sola carga concentrada en algún punto sobre la cuerda
consta de dos secciones lineales (triángulo) como se muestra en la figura siguiente,

Fig.15 Carga concentrada en x = e.

Las condiciones sobre la cuerda estan dadas de la siguiente forma,

Fl cos al = F2 cos a2 T,
Fi sin az + F2 sin a2 = —P.

d d
De aquí, tan al + tan a2 = ——
7 o bien
P -= —
T'
de donde se sigue que
—

d "e— e)e (2.1.6

Ti


2.1. Introducción 55

Conocida la deflexion d en el punto de carga, es sencillo determinar la deflexión

de la cuerda en cualquier punto, utilizando triángulos semejantes,
P(P— e)x
O<x<e
y(x) = p(inx)e (2.1.7)
e < x <
TI
Por lo tanto, podemos decir que la deflexión de una cuerda en un punto x, debida a
una carga concentrada de magnitud P aplicada en un punto x = e, es igual a la deflexión
producida en el punto e por una carga igual aplicada en el punto x. Cuando se tiene
una carga unitaria, P = 1, es común usar la notacion G(x, ¿) para y(x) para indicar
que la deflexión depende tanto del punto e donde está aplicada la carga, como del punto
x donde se observa la deflexión. En este ejemplo, directamente se observa que G(x, e)
es simétrica en las dos variables x y e. Es decir, G(x, e) = G(e, x). Comúnmente
G(x, e) se conoce como función de influencia, ya que describe la influencia que una
carga concentrada en el punto e tiene en cualquier punto de la cuerda x. La simetría
de G(x, e) es un ejemplo del Principio de Reciprocidad de Maxwell-Rayleigh, el cual
está presente en muchos sistemas mecánicos y eléctricos.
Es fácil checar que por medio de la función de influencia, G(x, ¿), puede hallarse
una expresión para la deflexión de una cuerda bajo una carga con distribución más
general sin utilizar algún método convencional de ecuaciones diferenciales. Una forma
sencilla de hacerlo es la siguiente. Se divide el intervalo [0, El en n subintervalos, por
medio de los puntos lo = O, £1 , e9 ' ' en = 1, con Aei — la longitud del
subintervalo [ei_1 — ¿d. Además, sea s i un punto arbitrario en el i-ésimo subintervalo
y supongamos que la porción de carga distribuida que actúa sobre este subintervalo,
co(si )Aei , es concentrada en el punto e = s i . Así, la deflexión producida en un punto
x por esta carga es la magnitud de la carga multiplicada por la deflexión producida en
x por una carga unitaria en el punto e = s i . Esto es,
[w(si)Aed G(x, si ) (2.1.8)

Si sumamos las deflexiones infinitesimales producidas en el punto x por las cargas

infinitesimales concentradas (equivalentes a la carga real distribuida), obtenemos

w(si )G(x, (2.1.9)

i=1
56 CAPÍTULO 2. FUNCIONES DE GREEN

Torpando el límite cuando Al; -+ 0, la suma de Riemann (2.1.9) se convierte en la

integral,

y(x) = f w(e)G(e,x) de (2.1.10)

o
la cual representa a la deflexión de la cuerda en un punto arbitrario x. La función de
influencia, G(x, e), tiene la expresión
(1 - e)x
{ G <(x ,e) O<x< e
= T2
G(x, C) = (2.1.11)
x )e
G>(x, e) -= e <x<
TP
Una de las principales propiedades es que G(x,e) debe satisfacer las condiciones
en la frontera del problema; en este caso de la cuerda, las condiciones son de deflexión
cero en x = 0, .1. Claramente, de (2.1.11), tenemos que G(0, e) = e) = O para
todos los valores de e en [0,4 Además, cabe mencionar que G(x, e) es una función
continua de x sobre todo el intervalo [0,1]; el único punto de conflicto es cuando x
pero tomando los límites por la izquiera y derecha de G(x, C),

fim G(x,e)= lim ( t e) x - e)e

x-›r. Ti TI
( 2. 1.12)
(1 - x)e (1 - e)e
lim G (x , e) = fim
-›e+ x-+I-F TI TI
tenemos evidentemente que dichos límites son iguales entre si e iguales a G(e, e), lo
cual muestra que incluso en el punto problemático x = e la función G es continua.
Por otro lado, la derivada respecto a x de G(x, e) es discontinua en el punto x =
donde se tiene un salto; en este caso de magnitud -1/T (hacia abajo). Para verificar
esto debemos checar que G(x, e) es diferenciable en todos los puntos del intervalo [0,1],
excepto en x e, donde se tiene que
dG(x,e) (1- e) (P - e)
lim _ him
x->t- dx TI Ti
(2.1.13)
dG(x, e) -e = -e
11 1/1 - -1.1111
x -› e+ dx x-›e + Ti TP
2.1. Introducción 57

Estos valores de los límites laterales no son iguales; su diferencia es

e)
—e (e —e) 1
(2.1.14)
TI TI

Por lo tanto, G(x, e) consta de dos expresiones lineales (el trángulo de la Fig.15)
continuas que satisfacen la ecuación diferencial homogénea Ty" = O en todos los puntos
del intervalo [0,1], excepto en x e. En el punto x = la segunda derivada no existe
en el sentido ordinario, pero si en el sentido distribucional, la cual de acuerdo a los
resultados del Capítulo 1, es
d2G(x, e) = s(x 1).
dx2 (2.1.15)
Las propiedades en las que hemos hecho enfásis no son accidentales de la función
de influencia específica obtenida, sino todo lo contrario; las funciones de influencia, que
de ahora en adelante llamaremos funciones de Green, son un tipo especial de funciones
asociadas con las ecuaciones diferenciales que precisamente satisfacen todas y cada una
de las propiedades descritas.

2.1.2 Solución Fundamental

Consideraremos ecuaciones diferenciales sobre todo el espacio completo y un
método poderoso para resolverlas basado en el concepto solución fundamental.

Sea L un operador diferencial lineal arbitrario de orden r, (ecuación (1.3.2)), y sea

Lu = f (x), x E S2 C (2.1.16)
una ecuación diferencial no homogénea en el sentido distribucional; u y f son funciones
generalizadas.

Definición 2.1 Una distribución u es solución de (2.1.16) sobre Q si

Tu [L* 0] = T f [0] para cada O E T; (2.1.17)

donde L* es el operador adjunto de L, fórmulas (1.3.18) y (1.3.19). Cuando f es
continua, u es una solución clásica de (2.1.16) si u E C'(Q). Mientras que si f es
localmente integrable (o una distribución más general), entonces u se dice ser solución
débil o generalizada de (2.1.16) en a
58 CAPÍTULO 2. FUNCIONES DE GREEN

Definición 2.2 Una solución fundamental para un operador diferencial L con polo o
singularidad en 4 es una solución de la ecuación no homogénea
Lu = 8(x — 4), (2.1.18)
donde 4 es considerado como un vector de parámetros.

Observaciones:
La ecuación (2.1.18) tiene que ser forsozamente considerada en el sentido dis-
tribucional. La solución de (2.1.18) se denota por G(x, 4) y se refiere como la
función de Green asociada con la ecuación diferencial producida por el opera-
dor L. Frecuentemente (pero no necesariamente siempre) G corresponde a una
función localmente integrable de x. En cada caso, de acuerdo a (2.1.17), G
satisface (2.1.18) si y solo si

TG [Ls O] = 6[0] para cada E D. (2.1.19)

La ecuación (2.1.18) usualmente tendrá muchas soluciones, difiriendo una de
otra por una solución de la ecuación homogénea y del tipo de condiciones dadas.
Para problemas en los cuales IR" sea considerado como un espacio isotrópico,
frecuentemente se selecciona una solución particular sobre la base de la simetría
y comportamiento en el infinito Por su parte, si una de las variables es el tiempo,
se tiene que utilizar el principio de causalidad como criterio apropiado respecto
a tal coordenada.
3. Si L es de coeficientes constantes, es suficiente determinar la solución fundamen-
tal con polo en el origen, G(x, O) = G(x), y trasladar para obtener la solución
fundamental con singularidad en 4, G(x, 4) = G(x — 4).

Son dos las partes principales para determinar una solución fundamental. Primero
debemos construir, frecuentemente por medios intuitivos (no existe nada global y for-
mal, depende de la dimensión, del problema y sus condiciones) un candidato y segundo,
verificar que se satisface (2.1.19). La primera parte a su vez puede dividirse en dos
pasos: resolver la ecuación homogénea para x 4, tomando en cuenta las consi-
deraciones físicas y entonces analizar en una vecindad alrededor de la singularidad en
x = 4, usando la forma integral de (2.1.18).


2.1. Introducción 59

22.3 Construcción y Propiedades

Como el problema de determinar una función de Green depende del caso conside-
rado, vamos a describir los pasos de la construcción de una función de Green para la
ecuación diferencial de segundo orden

Ly = {—d (p(x)—
d + q(x)jy(x) = —f (x), (2.1.20)
dx dx
sujeta a las condiciones homogéneas
a i y(a) + a2y'(a) = 0,
(2.1.21)
b i y(b) + b2 y'(b) = 0,

ya que a final de cuentas ésta es una de las ecuaciones más importantes en las apli-
caciones y de hecho muchas ecuaciones más generales (orden superior) se reducen a
ella.

Definición 2.3 Una función G(x, e) que cumple con las propiedades:

G(x, 1) satisface la ecuación diferencial homogénea (2.1.20) en a < x <

e < x < b.
G(a, = G(b, e) = 0.
G(x, e) es una función continua de x en todo el intervalo a < x < b.
dG(x , e)
es continua en a < x < e, e < x < b, pero tiene una discontinuidad de
dx
dG(x,e)re+ 1
salto de magnitud — en x = e. Esto es,
P(e) dx i x,e _ p(e)

se llama funcion de Green del problema definido por la ecuación diferencial (2.1.20) y
las condiciones de frontera (2.1.21). Además, la función de Green tiene la propiedad
de simetría, G(x, e) = G(5, x) (Principio de Reciprocidad de Maxwell).

Observación: Por definición, G(x, e) es la solución de la ecuación Ly = — fe(x),

donde fe (x) es cero fuera del intervalo ix — SI < e y en dicho intervalo satisface
1.11-c
e_, ) dx = 1
Jf,(x
 CAPÍTULO 2. PUNCIONES DE GREEN

por lo que fc (x) = 6(x — E). Entonces, para x G(x, e) satisface la ecuación
homogenea Ly = 0 y cuando x = e, se tiene que

d +e €
— p(x) dG(x, C) dx + f q(x)G(x, e) dx = f — f< (x) dx.
dx dx fe
Esto es,
dG(x,e)
p(x)dx + f q(x)G(x, e) dx -= — 1
X=1 — Z —E
Si asumimos que G(x, e) es continua y diferenciable en x = e, entonces en el límite
cuando e -+ 0, tenemos que
x=e+
dG(x, e) —1
=--
dx P(e)
lo cual confirma que la derivada de G(x, e) tiene una discontinuidad de salto en x = e.

La solución de la ecuación no homogénea (2.1.20) está dada en términos de la

función de Green sobre el intervalo [a, b] mediante el proceso de convolución.

Teorema 2.1.1 La solución y(x) del problema con valores en la frontera, ecuación
(2.1.20), y condiciones (2.1.21) está dada por la integral

y(x) = f bG(x, e) f(e) de. (2.1.22)

Nota: Como G(x, e) G (e , x) al momento de integrar se debe intercambiar el papel

de la variable y el parametro en la función de Green G(x, E).

Dem Para probar el teorema, utilizaremos la Fórmula de Euler-Leibnitz de diferen-

ciación respecto a un parámetro, ecuación (3.3.6), escribiendo (2.1.22) en la forma

y(x) = I
a
G(x, e) f (e) + f G(x, f (e) de. (2.1.23)
obteniendo

dG(x , e) Lb
(x dGE)
y1(x) = J x f (e) de + f (e) de + [G (x,.- ) f (x) — G(x, x+ )f (x)]
dx dx
(2.1.24)
 2.1. Introducción 61

Como G(x, e) es continua en x = e, se tiene que G (x , x + ) = G (x , x- ) por lo que

dG(x, e)
(x) = dx f (e) de (2.1.25)
a

Nota: Es claro que G(x, O satisface las condiciones de frontera (2.1.21):

aiy(a) + azy i (a) = f [aiG(0, e) + az dG d( O; y) f (e) de = 0,

a

,
( b y)
dGdbx,
biy(b) + b2 y 1 (b) = f [cli G(b, e) + az
a x
f (o)de0.
Ahora diferenciando y' (x), ecuación (2.1.25), se tiene que

d (2, [dG ( d ,
x xx—) f (x) dG( x , x +)
(2.1.26)
Y" (X) = f b d2 GdX f (e) de + dx f (x) I •

z={+ ¿=x+
dG (x , —1 dG , x) —1
Pero, como es equivalente a = , entonces
dx x=z- 19(e) dx e=a- p(x)
(2.1.26) se reduce a
b d(2,1) f (x)
y" (x) = (2.1.27)
a d2 f (e) de p(x)

Tomando el operador L de (2.1.20) en la forma Ly p(x)y" (x) + p' (x)y i (x) + q(x)y(x)
y sustituyendo (2.1.25)y (2.1.27), tenemos

,e) + q(x)(7(x,e)] f(e) de

, d(7(x
L(y) = — f (x) + fb [p ( x) d2G(X,I)
a
d2 P
„
dx

que podemos reescribir como

rb
L ( y ) = (x) L (G(x , e)) (e) de-

Pero por definicion se tiene que L(G (x , e)) = O para x e, obteniendo entonces
L(y) = — f (x). (2.1.28)
62 CAPÍTULO 2. FUNCIONES DE GREEN

La expresión de la función de Green para dos intervalos separados por x = e se escribe

como
G<(x,e) a < x < e,

G(x, = (2.1.29)
G>(x,e) x < b.

La condición de continuidad es expresa como G<(e- , e) = G> (e+ , e) y la condicion de

dG,(1± dG<(e- 1—
salto
dx dx 13(1)
Ejemplo 14: Consideremos el problema

y" = —x

con condiciones de extremo y(0) = 0, y(1) = O.

Para un valor fijo de e la función de Green debe satisfacer la ecuación homogénea

G"(x, e) = O. Además, debe cumplir con la condición de continuidad y de salto en
x = e. Tal función de Green (caso particular del problema de la cuerda, ecuación
(2.1.5)) está dada por

x(1— e) O < x < e,

G(x,e) =
e(1 — x) e < x < 1.
Así, como f (x) = x en (2.1.20), entonces el Teorema 2.1.1 mediante la ecuación (2.1.23)
implica que la solución está dada por

y(x) = f G(x, f (e) + f G(x, f (e) de

ci
1 x(1 — x2)
= e(1 x)e de + f x(1 — e)e =
fo 6


2.1. Introducción 63

Tina forma usual de construcción de una función de Green para la ecuación de

segundo orden (2.1.20), es mediante la utilización del Wronskiano del sistema (Técnicas
de Álgebra Lineal), con una ligera modificación. Consideremos la ecuación homogénea
asociada al problema (2.1.20) y condiciones de frontera (2.1.21). Supongamos que
111( x ) es solución de tal ecuación homogénea, la cual satisface a i y i (a) + (a) = O
y similarmente y2 (x) satisfaciendo la condición b i y2 (b) + b2 y5(b) = O, ambas válidas
en el intervalo [a, b] y yi , y2 linealmente independientes, de tal forma que la solución
general de la ecuación homogénea es y(x) = ci yi (x)+ c2y2 (x). Así, la función de Green
buscada se puede escribir como

Ci(e)Yi( x ) para x < e,

G(x, e ) = (2.1.30)
C2(e)Y2( x ) para e < x.

De la condición de continuidad en x = e se tiene la siguiente relación

C2 (e)y2 (e) — C2 (e)yi (e) = 0 (2.1.31)

y de la condición de discontinuidad de la derivada en x = e se tiene que

1 (2.1.32)
C2(e)yi(e) — ci(e)1(e) = —p(e)
Las des ecuaciones lineales en las incoginitas C 11,(1) y C2(e) tienen por solución

y2 (e)
c' (e) — rcerw(yl, y2; e)
(2.1.33)
yi(C)
C2(C) = pww(yi, y2; e)
donde W(y i , y2 ; e) es el Wronskiano (determinante del sistema) dado por

W (Yi, Yz, e) = ( x )V2( x ) — (x)Y1(x) ( 2.1.34)

Nota: Las soluciones yi y y2 son linealmente independiente si el Wronskiano es

diferente de cero.
 64 CAPÍTULO 2. FUNCIONES DE GREE

Por lo tanto, la funcion de Green tiene la expresión

(x )Y2 (e)
para x <
P(e) W ( yi, Yz; e)
G(x, e) ={ (2.1.35)
Yi(e)Y2(x)
para x > e.
P(e) W ( yi, yz, e)

Ejemplo 15: Consideremos el siguiente PVF

y" + y = —1
y(0) = O, y(1) = O.
Las soluciones estan dadas por

Yi(x) = sen x para O < x < e,

y2 (x) = cos x
ri
para e < x < 1.

Como W(yi , y2 ; e) = y]. (e)yl(e) — y2 (e)yl (e) = —1, entonces

sen x cos e para x < e,

sen e cos x para x > e.

De la misma manera que antes se puede ver que G(x, C) es simétrica en x = e.

Esto es, G(x, e) = G(e,x). La solución del PVF está dada por

y(x) = f G(x,e)f (e) + f i; G(x, e) f (e) de

= —1 + sen x + cos x.
o
En el próximo capítulo se darán otros ejemplos de la solución fundamental para
ecuaciones ordinarias ya que por ahora es más importante y quizá hasta más interesante
considerar el caso de dimensiones mayores.
2.2. Funciones de Green en Dimensiones Mayores 65

2.2 Funciones de Green en Dimensiones Mayores

Resolver un problema con condiciones en la frontera significa determinar una función
que satisfaga una ecuación diferencial dada y condiciones de frontera específicas. El
problema es independiente del tiempo, es decir solo involucra las coordenadas es-
paciales. Los problemas con condiciones iniciales (hay dependencia del tiempo) se
describen con ecuaciones diferenciales parciales de tipo hiperbólico, mientras que los
problemas con condiciones en la frontera estan asociados con ecuaciones diferenciales
parciales de tipo elípticas. La diferencia o el contraste entre los problemas con condi-
ciones iniciales y los problemas con condiciones en la frontera es que estos últimos son
mucho mas difíciles de resolver. Así, para estudiar funciones de Green en dimensiones
mayores a uno es necesario diferenciar entre unos y otros e introducir conceptos que
nos ayuden en su solución.
Existen varias formas de resolver problemas con condiciones en la frontera, depen-
diendo del tipo de condiciones que se tengan. Cuatro son los problemas fundamentales
los cuales a continuación se enumeran.

* Primer Problema con Valores en la Frontera (Problema de Dirichlet):

Consiste en determinar una función u, que además de satisfacer la ecuación diferencial
considerada, en una región Q, satisface
u = gi (x) sobre OQ. (2.2.1)
donde gi (x) es una función continua sobre la frontera OQ del dominio Q, el interior
de una superficie suave por piezas. Geométricamente, podemos dar una interpretación
simple a la solución u del problema de Dirichlet, como un estado estacionario con
valores fijos en la frontera.

* Segundo Problema con Valores en la Frontera (Problema de Neunaann)-

Consiste en determinar una función u definida en O, tal que
Ou
92(x) sobre DO con f f (x) dS = O, (2.2.2)
an
u
donde denota la derivada direccional a lo largo de la normal de la frontera DO. La
On
condición fan f(x) dS = O es conocida como condición de compatibilidad. Física y
geométricamente, la condición de compatibilidad se interpreta como el flujo neto que
cruza por la frontera (no hay fuentes ni sumideros).
 66 CAPÍTULO 2. FUNCIONES DE GREEN

* Tercer Problema con valores en la Frontera ondiciones Mixtas): Consiste

en determinar una función u en 52 que satisface:

Bu
— + q(x)u = g(x) sobre a (2.2.3)
On
donde g y q son funciones continuas. La interpretación física de esta condición es de
radiación a través de la frontera.

* Cuarto Problema con Valores en la Frontera (Problema de Robin): Consiste

en determinar una función u que satisfaga condiciones en la frontera de diferente tipo
sobre diferentes porciones de la frontera de S2, por ejemplo

u = (x) sobre 091 y u = f2 (x) sobre 022, (2.2.4)

donde 812 as-21 U 0Q2.

El problema que con mayor frecuencia se aborda en la aplicación es el de Dirichlet.

1h'
Sin embargo, no todos los problemas de Dirichlet tienen solución; la existencia de una
solución depende de la geometría de la región considerada.

Definición 2.4 Una función u se dice ser una función armónica en un dominio
si satisface la ecuación de Laplace Au = O, con sus dos primeras derivadas continuas
en ft.

Esta es una de las principales ecuációnes en las aplicaciones, junto con la ecuación
no homogénea relacionada, Au = f (x), llamada ecuación de Poisson.

2.2.1 Funciones de Green en dos Dimensiones.

Es momento ya de construir funciones de Green para problemas modelados por

ecuaciones diferenciales parciales, considerando primeramente el caso 2-dimensional,
enfocandonos en específico al problema de Dirichlet,

Au = f (x, y) en SZ c 1R2,
(2.2.5)
u = g(x, y) sobre an.
2.2. Funciones de Green en Dimensiones Mayores 67

La solución en este caso está dada por

aG
u(x, y) = f G(x, y, e, Ti me, n) dedn + g— ds, (2.2.6)
52 On
donde G es la función de Green asociada y n el vector normal a la frontera an de la
región 12. Como la solución u(x, y) se determina una vez que la función de Green es
conocida, el primer paso se concreta a encontrar dicha función. La funcion de Green
para el problema de Dirichlet satisface

i) AG 5(x e)5(y — ri ) en SI

ü) G = O sobre az.
En el caso cuando se considera el operador de Laplace (2.2.5) la función de Green
además satisface

G(x, y, e, n) = G(1, rl, x, y).

8G
G es continua en sus argumentos, tiene discontinuidad en el punto
pero
( e, 77) , lo cual se especifica por

lim
e-10
19G ds = 1,
—
8n
donde n/ es la normal unitaria al círculo Ce : (x — 1)2 + (y = €2.

La función de Green en este caso se interpreta como la respuesta del sistema en el

conjunto de puntos (x, y) debido a la función 8, fuente en el punto (e, G es continua
en S2 casi dondequiera, y su primera y segunda derivada son continuas en ft excepto
en (e, n) Las últimas dos propiedades para la función de Green se establecen en base
a los dos resultados siguientes, los cuales también son válidos en dimensiones mayores.

Teorema 2.2.1 (Principio de Reciprocidad de Maxwell-Rayleigh) La función

de Green asociada al operador de Laplace (2.2.5) es simetrica,

G(x, 1) = G(1, x). (2.2.7)

 68 CAP ÍT ULO 2. F UNCIONES D E GREE N

D QM. Partiendo de la fórmula de Green, identidad conocida del Cálculo Vectorial,

ao
(Çbao,b — 7,b00) dS = fa ( 0 ))0— ds,
0 8n 8n
an
considerando que q5 = G(x, y, e, n) y IP = G(x, y, t'', n'), donde ( e', V) es arbitrario,
tenemos que

(G( x, y , e n)áG(x, y, e', '7) — G(x, y, e', n'),AG(x, y, e, 27)) dS

aG ( x ,:n, e ' ni) ,)0G(x, y, e, 77)\

) ds.
= (G(x,Y,e,r1) ci( , Y, E , an
an

Pero, como G(x, y, e, n) y G(x, y, e, n) se anula en la frontera alZ, entonces

1,11
fL (G(x, y, e, WAG(x, y, e', 71) — G(x, y, e, n)AG(x, y, e, n)) dS = O

I
como AG(x, y, e, 77) = 8(x—e)c5(y—n) y áG(x, y, e', 77) 5(x —19)(y—ng, entonces

fiz G(x, y, e, n)8(x — es)6(y — 77) dS = f£G(x, y, el , 77)5(x —1)5(y — dS,

de donde tenemos que G(e , , e, G(1,7 7 , e y por lo tanto G(x, 1) G(g, x).

0G
Teorema 2.2.2 — es discontinua en el punto j = (1,n), si
n
aG
lim — ds = 1, (2.2.8)
->0 8n
donde n' es la normal del circulo CE : (x — y e)2 + ( = E2.

Dem. Sea RE la región acotada por C. Integrando áG(x, y, e, 77) = 8(x — e)8(y — n)
sobre el disco R6 , tenemos

fiR, AG(x, y, e, n) dS = ffRe AG(x, y, e, 77) dS = ffRe 8(x — e)6.(y — 77) dS = 1.



2.2. Funciones de Green en Dimensiones Mayores 69

Si tomamos el límite cuando e —> O, tenemos que

lim ff AG(x, y, e, n) dS = 1.
E—>0 Re

Aplicando la fórmula de Green con = 1 y b = G(x, y, e, n), obtenemos

lim AG(x, y, e, 77) dS liM

aG fr, Y, e, n) ds,
E—>0
Re E—>0
ceani

de donde li m aG(x, Y, e, 77) ds = 1. Lo cual implica que aG

— es discontinua por el
—>0 ce ans an
teorema fundamental de integrales de linea.

2.2.2 Método de la Función de Green

Al buscar la función de Green de un problema determinado, es conveniente hacerlo

considerando la suma de una integral particular de la ecuación no homogénea y de la
solución de la ecuación homogénea asociada. Esto es,

G(x, = G i(x ,E) + G2(x, 4), (2.2.9)

donde G1 satisface .G 1 = 6(x — e) en a Mientras que G2 satisface AG 2 = O en f2.

Es claro que G satisface AG = 6(x — e). Para que G = O sobre 89 se requiere que
Gy = —G1 sobre 05-2 y G1 no necesariamente tiene que satisfacer condiciones en la

frontera.
Con estas propiedades satisfechas, podemos determinar la solución del problema
de Dirichlet (2.2.5) mediante el método de la función de Green.

Tomando 0(e, 17 ) = c(e, 77, 3c, y) y tP(e, 77) u(C, y) en la fórmula de Green, tenemos

fi n(G(e,n, x, y ) Ao(e, y) — u(E, n) A e G (1, y , x, y)) dAe

= (G(e, x , v) 8u(e ' 11) u ( e , 77) aG(1an

' 77 ' x '
Oft an Y))
70 CAPÍTULO 2. FUNCIONES DE GREEN

Como Aeu((, n) = f (e, n) y A1 G((,n, x, y) = 8(1 - x)5(n - y), entonces

fin (G (e n x Y).f (e, n) — u(e, n) 6 (e — x)( (n - y)) dechl

(G((,77, x, y) au(I ' 71) u((, n) OG (( , x , y)) ds

- n On On
Además, de las condiciones de frontera, se tiene que G = 0 y u = g sobre 85-2. Esto y
la simetría implican

, x,
OG(1en x y) (2.2.10)
u ( x , y) = G(1, 71; x , y ) f(e, de s*, + f y ( e, n) ds.
[f s-

2.2.3 Funciones de Green en tres Dimensiones.

El método de las funciones de Green inmediatamente se extiende al caso de tres

y hasta más dimensiones. Consideraremos solo el caso 3-dimensional ya que este es el
más frecuente en las aplicaciones reales.

Al igual que en dos dimensiones, consideraremos el problema de Dirichlet. La

función de Green para la ecuación de Poisson es la función que satisface

AG = 5(x - 05(y - y)8(z - () en 2 C Ift3.

G = O sobre an.

G(x, y , z , e7 77, = G ((, n, (7 X ) YI Z)•

9G
G es continua todos sus argumentos, pero — tiene discontinuidad en el punto
an
x = = ((, n, (), lo cual se especifica por

aG ,r,
hm — = (2.2.11)
E-H) an'
donde n' es la normal unitaria a la esfera S, : (x - ()2 + (y n + (z ()2 62.
2 2 Funciones de Green en Dimensiones Mayores 71

(procediendo como en el caso de dos dimensiones, la solución del problema de

Dirichlet
Au = f (x) en
(2.2.12)
u = g(x) sobre 3Q
está dada por
u(x, y, z) = ff f ( e, n, ()G (x, y, z; e, n, cledndC
(2.2.13)
0G(x y x; e n ()
+ g((, C) ' ' ' dS,
00 OIL

donde G(x , y, z ; 1, 77, = (x , Y, z; 1, 71, () G 2 (x, y, z; 77, (), con G1 satisfaciendo

AGi = 6(( — x)5(71 — y)5(C — z) y G2 satisfaciendo AG 2 = O en s-z.
 72 CAPÍTULO 2. FUNCIONES DE GILEÉ:

Comentario 2

iPp
A pesar de la importancia y la enorme potencia con la que nos equipa la uti-
lización de las funciones de Green en la determinanción de la solución de ecuaciones
diferenciales, es lamentable que en ningún curso de la carrera de la Licenciatura en
Matemáticas de la Universidad de Sonora ni siquiera se mencione. Quiero o creo en-
lV
im tender que esto principalmente se debe al hecho que quizá hace algún tiempo, cuando
la tecnología no estaba tan desarrollada como lo está actualmente, estudiar esta técnica
no era muy atractivo que digamos debido a que nos deja "a medio camino" al tener
expresiones para las soluciones en forma integral, muhcas de las cuales difíciles de eva-
luar (por lo menos manualmente) o tal vez se deba a la falta de interés o apatía tanto
de estudiantes como de profesores por estudiar conceptos relacionados con la Física.
Sin embargo, el enfoque de las funciones de Green es interesante por si mismo, por
las matemáticas involucradas y por supuesto por el aspecto físico. Ahora en nuestros
tiempos, el método de las funciones de Green ya no nos deja tan a medias, puesto que
existen paquetes de cálculo simbólico como por ejemplo Maple o Mathematica los
cuales representan una poderosa herramienta numérica, auxiliar a técnicas analíticas,
y tal es el caso del enfoque de las funciones de Green. Este trabajo lo elaboré pen-
sando precisamente en eso y espero sirva de punto de partida para un análisis mucho
más profundo de estas técnicas, incorporándolas en alguno de los cursos de ecuaciones
diferenciales.
 1BLIOTECA
DE CIENCIA& EXACTAS
Y NATURALES
EL SABER TAL ME ilIJOR
RARA MI GRANDEZA

Capítulo 3

Aplicaciones

En este capítulo se presentan algunas aplicaciones de la teoría de las funciones gen-

eralizadas con el fin de mostrar la potencia y la elegancia de sus métodos, los cuales,
en ocasiones son más simples y directos que otros enfoques. De hecho para ciertos
problemas no hay método alternativo para determinar la solución respectiva, como
por ejemplo algunos que involucran ecuaciones diferenciales parciales. Esto sucede con
frecuencia en aerodinámica, mecánica de fluidos, electrodinámica y acústica. Muchas
situaciones físicas (problemas reales) se modelan con ecuaciones o sistemas de ecua-
ciones diferenciales no homogéneas cuyos términos no homogéneos se identifican con
fuentes (sumideros) o singularidades que convienen ser descritos mediante distribu-
ciones. Resolver esta clase de problemas significa determinar soluciones generalizadas,
lo cual solo es posible con la herramienta proporcionada por la teoría de distribuciones.
Además, también es posible encontrar funciones generalizadas en otro contexto; desde
la derivación de leyes de conservación hasta en problemas más sofisticados como los de
turbulencia en dinámica de fluidos.

73
 74 CAPÍTULO 3. APLICACIONES

3.1 Soluciones Generalizadas de Ecuaciones Diferenciales

Uno de los orígenes de la teoría de distribuciones fue la necesidad de determinar
soluciones de ecuaciones diferenciales tanto ordinarias como parciales fuera del marco
teorico del análisis clásico; funciones altamente discontinuas o singulares, las cuales son
llamadas soluciones débiles o generalizadas. En esta parte del trabajo se consideran
algunos de estos problemas. No se presenta por completo o con detalle todo el proce-
dimiento en la obtención de las soluciones. Simple y sencillamente pueden considerarse
ti
como ejercicios ya que nuestro interés principal es dar una introducción de la teoría,
esperando motivar un análisis más profundo de la misma.

3.1.1 Ecuaciones Diferenciales Ordinarias

.119 11 a) Sabemos, de las ecuaciones diferenciales elementales [151, que la solución general
11„;
de la más simple de las ecuaciones de orden n, O') = O, es un polinomio de grado
11„
n, P„(x), mientras que si consideramos la ecuación no homogénea, entonces la
solución general es y(x) = yg(x) + yp(x), donde yH es la solución de la ecuación
homogénea y yp una solución particular Para la ecuación no homogénea

y(n) = (5(x - a)

es fácil verificar que una solución particular es

yp(x)=- (( h(x - a),

x n- -
a1)!

donde h es la la función de Heaviside, ya que para n > 1,

(x - a)n-2
Yip (x) = (n - 2)! x > a,

=0 x < a,

puesto que yp es continua para n > 1.

3.1. Soluciones Generalizadas de Ecuaciones Diferenciales 75

I» la misma forma, para O < k < n,

ynx) = (x —
x > a,
(n — k — 1)!
=0 x < a,

pero para k n, yp(n ) (x) = 15(x — a). Estas ideas se pueden resumir de la
siguiente manera:

"La solución general de la ecuación diferencial no homogénea O.) 8(x — a)

es una función ordinaria de x que satisface la ecuación homogénea y(10 = O en
x < a, x > a. Esta función y sus primeras (n — 2) derivadas son continuas en
x = a, pero la (n — 1)-ava derivada tiene un salto de magnitud unitaria". Si
n = 1, la función misma tiene un salto de magnitud uno en x = a. Esto es, la
función de Heaviside se considera como la antiderivada de la delta de Dirac.

Resultados similares son válidos para la solución de la ecuación diferencial de

orden n,
dn-1
Ly ={a„(x)c7F, + (x)dxn_i cie(x)] y = (5(x — e),

si las funciones an, , an son continuas en x = e. Fácilmente se puede verificar

que para x < e, x > e, la función y(x) es la solución de la ecuación homogénea
Ly = O y que para x = e la función y sus primeras (n-2) derivadas son continuas,
pero su (n — 1)-ava derivada tiene discontinuidad de salto en ese punto. La
solución de la ecuación Ly = 8(x — C) se conoce como solución fundamental o
función de Green (Capítulo 2). La solución de la ecuación diferencial ordinaria
no homogénea de orden n más general, Ly = -f(x), está dada en términos de
la función de Green. Resultaría bastante interesante contrastar el Método de
Variación de Parámetros contra el enfoque de las funciones de Green.

b) Como ilustración de la ecuación general de 2 d° orden del inciso a), consideremos

el problema de encontrar la solución fundamental de la ecuación
(3.1.1)
Y" + A2Y = +8 (x —
76 CAPÍTULO S. APLICACIONES

Resolveremos primeramente con la combinación de signos,

G" — A 2G = —5(x — e).

Para x e, tenemos que G satisface la ecuación homogénea

G" — A 2 G = O.

Así, el problema se divide en dos partes:

G> (x, e) = A(e)eAx + B(e)e-Ax para x > e, G< G>

G<(x, e) = C(e)eAx + D(e)e-Az para x <
z=e

Cuando x --> co, G> (x, e) = 0 lo cual implica que A(C) = O y cuando x —› —oo,
G<(x, e) = O por lo que D(C) = O. Así,

G,(x, e) = B(e)e-Áz para x >

G<(x, e) = C(e)e2' para x <

La continuidad, G >(x ,Oix=e+ = G<(x,e)1.,1_, implica que B(C) = C(e)&x{,

mientras que la condicion de salto

dG> dG<
= —1
dxdxx=e+ x=e-
implica que C(C) = e-412A. Por lo tanto, tenemos que

G(x, > (x , e)
= G
G<(x, e)

la cual finalmente se puede escribir como

G(x, e) = 1 . (3.1.2)
2A
3.1. Soluciones Generalizadas de Ecuaciones Diferenciales 77

En caso de considerar los signos positivos (solución fundamental para el oscilador

armónico),

Gil + A2 G = 8(x — e)

tenemos que la función de Green es

G(x, e) = (3.1.3)
2iA

Si se considera el problema en un dominio finito, entonces se necesita también

especificar las condiciones de frontera. Por ejemplo, para condiciones de frontera
cero en x = 0, 1, para la última ecuación considerada, "construyamos" una
e
solución generalizada de la siguiente forma: para x la solución de la ecuación
homogénea correspondiente que satisface la condición de frontera y(0) = O es
y = sen Ax, mientras que la solución que satisface y(1) = O es y = sen A(1— x).
Esto es, tenemos que

sen Ax
y(x) {A
B sen A(1 — x)

donde A y 13 son constantes que se determinan de las condiciones que se tienen

en x e; se requiere que y(x) sea continua en x e y que su primer derivada
tenga un salto de magnitud uno en x = e, lo cual implica que
sen A(1 — e)
A—
sen
B—
AD
D sen Ae cos A(1 — e) + sen A(1 — e) cos Ae sen A.
Por lo tanto,
sen Ax sen A(1 — e)
A sen A
y (x , e) = sen A(1 — x) sen Ae
A sen A

Observemos que y(x, e) = y(e, x), lo cual ilustra el Teorema de Reciprocidad de

Maxwell (Sección 2.2).


78 CAPÍTULO 3. APLICACIONES

3.1.2 Ecuaciones Diferenciales Parciales

a) Consideremos el problema de resolver la ecuación de Poisson,
áu f (x), (3.1.4)
sobre todo el espacio infinito
Para obtener su solución explícitamente (representación integral), primeramente
determinaremos la función de Green asociada G(ae,1), la cual satisface
áG = —5(x — ¿). (3.1.5)
Como los resultados son diferentes dependiendo de la dimensión del espacio con-
siderado, se presenta el resultado simutáneamente para dos y tres dimensiones y luego
se generaliza.
Sea r (p) la distancia radial desde x en dos (tres) dimensiones. Esto es,
r = Il x —
(P =x
ii — x 2
= Vfri — ei.) 2 + (xz —e2)2
2 x -
= NA i — 11) + (x2 — e2) + ( s — 1 3) 2 )•

Supongamos que G(x, 4) depende solo de r (p):

G (x , 1) = G(r) (G(x G(P) ) •
Lejos de la fuente (r O o p O) tenemos AG(x, 4) = O. En dos dimensiones, bus-
caremos soluciones de la ecuación de Laplace con simetría circular (en tres dimensiones
con simetría esférica). Así, en coordenadas polares (esféricas), tenemos que resolver la
ecuación

1 d dG
(r
— — — ) = O,
r dr dr
rO
2 d 2dG\
c ) O,

cuya solución general se obtiene por integración directa,

G(r)= cija r + c2 (G(p)= + c 4) •

P
En cada caso, una de las constantes absorbe la singularidad de la fuente, la cual
se obtiene integrando (3.1.5) alrededor de una pequeña circunferencia (esfera),

ff áG dA = — 1
(iff akG dV = —1)


3.1. Soluciones Generalizadas de Ecuaciones Diferenciales 79

e
Aplicando el Teorema de la Divergencia, tenemos que

ff V . (VG) dA = f VG•n ds = —1 (fll V • (VG) dV = f VG • n dS = —1) •

En dos dimensiones la derivada de la función de Green en la dirección normal a una

circunferencia, VG • n, es 8G /8r (en tres dimensiones acíap), ya que solo hay depen-
dencia radial. Como sobre la circunferencia (esfera) el radio es constante, entonces

8G ( G
27r-- = —1 47,11— = —1) .
ap
En la mayoría de los problemas que involucran singularidades es necesario conside-
rar el límite de una circunferencia (esfera) infinitamente pequeña. Así, la condición de
singularidad se expresa como

, OG 1 (hm _1
Jim = — —
r->0 8r 27- a p 47)
Por consiguiente, tenemos que

1 ( 1)
— c3 =
271-
Las constantes c2 y c4 son arbitrarias e indican que la función de Green para la ecuación
de Poisson para el espacio no acotado es única, excepto por una constante aditiva. Por
conveniencia e2 = c4 = O, y por lo tanto,

1
G(x, 1) = r (G(x,1) = .

En general, la función de Green para la ecuación de Poisson está dada por

C In r, n = 2,
G(x, = (3.1.6)
K r2-n,
n > 2,

donde C =- 4
-, K = 11 (a)
2
(área de la esfera unitaria en IR") y r IIa —
2(n — 2) n/2

80 CAPÍTULO 3. APLICACIONES

Para obtener la solución de la ecuación de Poisson (3.1.4) sobre todo el espacio

infinitas mediante la función de Green asociada, se utiliza la fórmula de Green inte-
grando sobre una gran circunferencia (esfera) y luego se toma el límite cuando el radio
tiende a infinito, dando como resultado

u(x) = ff f (1)G(x, dAe (u(x) = f (1)G(x, dV e) . (3.1.7)

Es posible también resolver la ecuación de Poisson (3.1.4) para un dominio semi-

infinito utilizando el método de las imágenes o incluso sobre un dominio acotado,
pero resulta que el problema depende de la geometría y del tipo de condiciones de
frontera, por lo que se tiene que considerar en particular cada problema en específico.
Dedicaremos la siguiente sección para ilustrar este tipo de problemas.

b) Consideremos ahora el problema de determinar la solución de la ecuación de

Helmholtz,
k2 u = f (x). (3.1.8)
Con signo negativo la ecuación (3.1.8) describe procesos de difusión estacionaria con
fase específica, mientras que con signo positivo es una reducción de la ecuación de onda
al considerar el problema con una sola frecuencia (frecuencia monocromática). Este es
un ejemplo clásico de problema radial.
Consideraremos el problema sobre todo el espacio 3-dimensional, con fuente esta-
cionaria concentrada en x = y tomaremos por convenciencia primeramente el signo
negativo. El flujo de salida a través de una pequeña esfera alrededor de la fuente es el
mismo que el de entrada. Esto es,
lim —
E-40
fLeike

8u
en
dS 1. (3.1.9)

La función u se anula lejos de la fuente y depende solo de la distancia radial desde la

fuente. Así, en coordenadas esféricas u = u(p) satisface la ecuación
1 d (2du „,2
P dio) U ° p > 0, (3.1.10)
P2 dP
lo cual hace que (3.1.9) se reduzca a
du
—1 = lim 42re2 (— dS. (3.1.11)
dp)
3.2. Método de las Imágenes 81

Haciendo la sustitución u = v I p, tenemos que (3.1.10) se transforma en

- q2 V =O

cuya solución es de la forma y = cle-qP + c2e9P . Debido al comportamiento cuando

p -> oo, tenemos que u = Ae-9P . La condición (3.1.11) implica que A = 1/4/r. Así,
u = e-9P/4rrp. Por lo tanto, la función de Green para la ecuación de Helmholtz (3.1.8)
con signo negativo es
e-qiix-Iii
G(x, e). -
(3.1.12)

1
Nótese que si q -› O, entonces G -› lo cual concuerda con el resultado
II

471-11x -
obtenido en a). Es decir, se obtiene de nuevo (3.1.6) la función de Green para la
ecuación de Poisson.
Similarmente, considerando el signo positivo en (3.1.8), obtenemos
eigllx-ell
G(x, 4) — (3.1.13)
- ¿II"
En dos dimensiones, obtener la función de Green para la ecuación de Helmholtz
es más complicado; el resultado queda en términos de funciones de Bessel. En general,
la función de Green depende de la paridad de la dimensión; para dimensión impar, la
función de Green es una combinación de funciones de Bessel esféricas, mientras que
para dimensión par, la función de Green es una combinación de funciones de Bessel
cilíndricas.

3.2 Método de las Imágenes

El método de las imágenes es una técnica que se utiliza para obtener una función de
Green (construcción de) para un dominio finito o semi-infinito, a partir de un dominio
infinito El procedimiento es una buena alternativa cuando se requiere resolver una
ecuación diferencial sobre un dominio finito sobre el que se especifican condiciones
de frontera, solo que tiene la desventaja de que es aplicable únicamente cuando la
frontera tiene simetría simple, por ejemplo respecto a un punto o con respecto a un
eje. A continuacion consideraremos algunos ejemplos.
82 CAPÍTULO 3. APLICACIONES

3.2.1 'Función de Green sobre un Semi-Plano

Consideremos el ejemplo clásico en el que se construye una funcion de Green en

dos dimensiones para la conocida ecuación de Poisson en el semi-plano superior

Lu = f (x), (3.2.1)

x2 > O, sujeta a condiciones no homogéneas en x2 = O,

u(x i3 O ) = p(x). (3.2.2)

La función de Green correspondiente satisface

LIG(x , 1) = —6(x — e)
(3.2.3)
G(xi, o; 12) = o.

En el semi-plano x2 > O únicamente se tiene una fuente concentrada en x = 4 y la

1
función de Green para el espacio completo, — 2 7r In Ilx — no satisface la condición
de frontera en x2 = O, pero el método de las imágenes precisamente proporciona el
camino que permite que la solución satisfaga las condiciones de frontera especificadas.

a) b)
Fig.16 a) Fuente real (positiva) y su imagen b)

Se considera el problema sobre la totalidad del espacio, proponiendo una fuente imagen
(negativa) de la fuente 8(x -4) en x = 4. Sea —5(x —4*) la fuente imagen en x =
donde (* = Fig.16b). Esto implica que la función de Green buscada satisface
la ecuación

G (x , = — [8(x — — (5(x — •)] . (3.2.4)



3.2. Método de las Imágenes 83

Por erprincipio de superposición, la respuesta de tales fuentes es la suma de las dos

respuestas individuales. Así, la función de Green está dada por
1
27r in lix —
G(x ,= — — 27rln
— lix — (3.2.5)
Por simetría, la respuesta en x 2 = O debido a la fuente en x = ¿* debe ser el. negativo
de la respuesta debido a la fuente original en x = 1, de tal forma que se satisface la
condición homogénea en x 2 = O. Por lo tanto, la función de Green para el espacio
semi-infinito x 2 > O puede escribirse como
1
G(x, C) —
II x — Cil
ln „ „ =
1
ln
— + (x2 + e2)2 6)2 (3.2.6)
27r ilx — ihr (xl _ ei)2 ,(_
. (x2 _ 12)2.
Ahora claramente se puede ver que en x 2 = O se tiene que G(x, ¿) = O.

En general, la solución de la ecuación de Poisson, con condiciones de frontera no

homogéneas, está dada en términos de la función de Green,
8G(x, e)
u(x) = ff f (1)G (x, dAt + f g(e) dS. (3.2.7)
a n an
En particular, para el caso cuando Q es el semi-plano x 2 > O, tenemos que
o. a
u(x) = ff f (Z)G(x, 4) dAt — f _co
g(1) —G(x, e)
12=0

De (3.2.6) tenemos que

aG(x,1) 1 X2

ae2 12=0 (xi — e1)2

Así, finalmente obtenemos la solución de la ecuación de Poisson (3.2.1), sujeta a la

condición de frontera (3.2.2)

u(x) = ff f (x)G(x , 1) dA. 1 — —

1 r g(1)
X2
(3.2.8)
_ce (xi — + x2 dei
Partiendo de este hecho, si considera la ecuación de Laplace para el semi-plano
infinito x 2 > O con condicion de frontera u(x i , O) = g(x i ), la solución está dada por
X2 r- 2 del.
u (x i, x2) = s(ei) (xi — ei)2 + x2 (3.2.9)

84 CAPÍTULO 3. APLICACIONES

3.2.2 El Problema de Robin

ilf1;" :1
Otro ejemplo clásico de problemas con condiciones de frontera en espacios semi-
infinitos es el llamado problema de Robin el cual se resuelve también usando el método
de las imágenes.

Consideremos la ecuación de Poisson (3.2.1) en el espacio 2-dimensional en la

región x i > O y x2 > O, sujeta a las condiciones no homogéneas
u(x1,x2)1„ = g1(x2),
au (3.2.10)
==g2(x 1) ,
On x2=-0

e.

a) b)
Fig.17 a) Fuente real y sus imágenes b)
La fuente se localiza en (xl , x2 ), sus imágenes en (x l , —x2 ), (—x l , —x 2 ), (x l , —x2),
como se ilustra en la Fig.17. Por inspección se puede construir la función de Green
para este problema,
1 [(x i + 6) 2 + (x 2 — e2) 2][(xi — ei) 2 + (x 2 + e2)2]
G(x, E) = 47r ln )2 (x2 [( x i e1 )2 + ( x 2 + e2)2] (3.2.11)
[(xi — ¿2)2]
la cual satisface condiciones homogéneas en la frontera.
Como la solución está expresada en términos de la función de Green,

u(x) = if f(x,E)G(x,c) dA l [G(x,c)—

au OG
an — u—
an j dS (3.2. 1 2)

entonces

u (x) = ff f (x, C)G(x, 1) dAe + / cc g2(W c(x, E) del —

o o
gi(6)
aG (xE)
a dez,
n (3.2.13)
donde G(x,) está dada por (3.2.11).
3.2. Método de las Imágenes 85

3.i3 Funcion de Green sobre un Círculo

La función de Green para la ecuación de Poisson sobre un círculo de radio a y con
condiciones de frontera cero sobre la circunferencia se establece como
AG (x , = 5(x — Z),
(3.2.14)
G (x, = O para ilx = a,
y por supuesto se determina por el método de las imágenes, ya que debido a una razón
geométrica presente, existe un punto imagen x = e tal que la respuesta a lo largo
de la circunferencia es constante. Cuando el punto de observacion de la fuente se
encuentra en el interior de la circunferencia, llamaremos a este caso problema interior
(ver Fig.18), en cambio si el punto de observación se encuentra fuera, entonces le
llamaremos problema exterior.

Fig.18 Fuente real e imagen para el círculo (Problema Interior).

Los problemas interior y exterior se resuelven de manera similar, la única diferencia
es al momento de evaluar las condiciones de frontera ya que el vector normal cambia
en cada caso su orientación. Determinaremos una expresion para el caso interior e
inmediatemente a partir de ésta determinaremos la expresión para el caso exterior.
Consideremos la función de Green sobre todo el espacio correspondiente a una
fuente en x = y una fuente imagen en x =
iáG(x, (5(x — 1) — 5(x — 1*) (3.2.15)
la cual, por el principio de superposicion, es la suma de las fuentes individuales,

G ( x , E) = —
27r — 21r
— ln — = 1 ln Ilx + (3.2.16)
+c 47r —
86 CAPÍTULO 3. APLICACIONES

donde c es una constante de homogeneidad que se determina a partir de la existencia

del pinto imagen e y del hecho que G(x, 4) = O sobre la circunferencia 114 = a. Para
que esto ocurra se necesita que

li x 111 2 = kiix (3.2.17)

1
lo cual implica que c -- In k. El punto imagen e * está a lo largo de la misma linea
zbr
radial que el punto fuente 4,
(3.2.18)

El ángulo e entre x y 4 (el mismo entre x y 4*) satisface la Ley de los Cosenos:

Il x -1I1 2 11 x 11 2 + 11411 2 — 2 114 IIcII cose (3.2.19)

itx -11 2 = IIxlI 2 + Per - 2114 iteii cose
La relación (3.2.17) es válida sobre la circunferencia H4 = a, usando (3.2.18), si

a 2 + re2 — 2are cos e = k(a 2 +72r£ — 2a7re cos e),

donde re = 11411, para todos los ángulos e. Esto lleva a dos ecuaciones para 7 y k
a2 + k(a2 + 72r1),
—2are = k( -2127r£),
a2
las cuales tienen por solución k = 1/7 y 7 = Así, el punto imagen se localiza en
e
a2
1* = 72-4'
e
a2
Nótese que III*11= —. Es decir, el producto de los radios de la fuente y su imagen
re
es el radio al cuadrado de la circunferencia. Por consiguiente, entre más cerca esté la
fuente del centro de la circunferencia, más lejos deberá estar el punto imagen.

Por lo tanto, la función de Green buscada está dada por

1 ln ( Ilx — 11 2 a2
G(x, e) — x 127.1 (3.2.20)
3.2. Método de las Imágenes 87

la cual, usando la Ley de los Cosenos, se puede expresar como

1 r2 +r — 2rre cos
G(x, 4) = — 1n a2 (3.2.21)
47r r2r2 + a4 — 2rrea2 cos e

donde e es el ángulo entre x y 4, r = 11x11 y re = III. De aquí se verifica que

G(x, 4) = O, sobre la circunferencia r = a.
Una vez obtenida la función de Green (3.2.21), se puede obtener la solución para
la ecuación de Poisson directamente, utilizando la fórmula general,

u(x) = ff f (4)G(x, 4) dAe + g(E)VeG(x, 4) • n ds. (3.2.22)

Debido a la geometría del problema, para evaluar las integrales se utilizan coordenadas
polares Para la integral de línea se tiene que ds = adO y así
a
g(I)VeG(x, 4) • n ds = f2n g(01) , .e- G(x, 4) a dOe. (3.2.23)
13 rt=G

Por lo tanto, la expresión de la solución de la ecuación de Poisson queda expresada

COMO

1 2n (i„ a2 — r 2
u(r, O) = ff f (r Oe)G (x 4) dAi.e + 91) dee,
27r ja gytel r2 + a2 — 2ra cos(0 —
(3.2.24)

donde O — = e. Esta expresión es valida para el caso interior; la normal apunta en

la dirección radial. Para el caso exterior se tiene que la normal apunta en direccion
opuesta, por lo que solo cambia el signo de la integral de línea en (3.2.24).

Partiendo de esta expresión, ahora determinaremos la solución de la ecuación de

Laplace con condición de frontera g(8) en el círculo r = a,
a2 r2
1 /27r
u(r, O) = ±-
27r g(01) clO e (3.2.25)
r2 + a2 — 2ra cos(61 — Be)

con signo positivo cuando el punto de observación de la fuente está en el interior y

signo negativo en el caso exterior.


88 CAPÍTULO 3. APLICACIONES

A la expresión (3.2.25) se le conoce como fórmula de Poisson, la cual usualmente

se eiqtresa en la forma más compacta,
1 f 2r
u(r, 9) — P(r, O — 01)11(9e) dOe, (3.2.26)
27r o
donde a
r2 — a2
P(r, O — Be) r2 (3.2.27)
a2 — 2ra cos(0 — Be)
se le llama Kernel de Poisson.
Similarmente, obtenemos que la función Green para el semi-círculo, O < r < a, O <
O < 7r, con G(x, ¿) O en la frontera, es

1 r2 + — 2rre cos(0 — Be)

Ina2r27,/
G(x,C) = —
47r a4
2rrea2 cos(0 — Be)
(3.2.28)
_ 1 in a2 r2 + re — 2rre cos(0 + Be)
47r r2r1. + a4 — 2rrea2 cos(0 + Be)

3.2.4 Función de Green sobre una Esfera

La función de Green para la ecuación de Poisson sobre una esfera se establece de

!J.
la misma forma que para el caso del círculo, ecuación (3.2.14), y tiene la forma
1 1 1 1
G(x , (3.2.29)
47r lix — 47r lix —
La relación entre y 1* se determina utilizando una propiedad geométrica simple.
Sea x un punto arbitrario de observación, el punto donde se localiza la fuente y E*
el punto donde se localiza la fuente imagen (Ver Fig.18). Por semejanza de triángulos
a2
se tiene la relación =- 114 . Como lixil = a y = re, entonces II = — y por
114 1111 1
lo tanto
1 re 1 1
G(x,) = (3.2.30)
47r Vr2 r/ a4
— 2rrea2 cos O 47r V r2 re2 — 2rre cos e

donde
cos O = • — cos cos Be + sen sen Be cos(0 — q5e), (3.2.31)
a re
li
donde O y son ángulos esféricos.
3.3. Algunas Aplicaciones Interesantes 89

3.3 Algunas Aplicaciones Interesantes

Para concluir el trabajo vamos a considerar algunos problemas interesantes, espe-
rando con esto ilustrar aún más el uso que tiene o puede tener la teoría de distribuciones.
El primero de los problemas que se presenta es la justificación de una fórmula de
derivación típica de cálculo elemental; la fórmula de derivación de Euler-Leibnitz y
su generalización. El segundo problema, también un resultado clásico de cálculo; el
Teorema de la Divergencia generalizado. El tercero de los problemas una propiedad que
se puede experimentar físicamente; un efecto de óptica. Además, se presentan algunas
representaciones de la solución de la ecuación de onda y por último, se presenta también
una fórmula; la llamada Fórmula de Kirchhoff, la cual aún cuando es muy necesaria
en las aplicaciones físicas es dificil de obtener y utilizar.

3.3.1 La Fórmula de Euler-Leibnitz

En cálculo se obtiene, mediante la utilización de la Regla de la Cadena multidimen-
sional la conocida Fórmula de Euler-Leibnitz para diferenciar una integral que depende
de un parámetro. Ahora veremos como hacerlo utilizando el concepto de diferenciación
generalizada. También, veremos que sólo con la teoría de las funciones generalizadas
es posible establecer el resultado de manera más general.

El problema es determinar la derivada con respecto al parámetro a de la integral

d 1 B(a)
I(a) f (x, a) dx, (3.3.1)
da A(a)
donde A, 13 y f son funciones continuas y B(a) > A(a) para a E [a, b]. Para hacerlo,
se introduce la función de ayuda
H(x, a) = h[x - A(a)]1i[B (a) - (3.3.2)
con h(•) es la función de Heaviside.

Observemos que H(x, a) = 1 cuando A(a) < x < B(a) y H(x, a) = O en otra parte.

Usando esta función de ayuda H(x, a), podemos escribir 1(a) como

d " aH
d fc: H (x , a) f (x, a) dx = J I Ca +
al f:t) dx,
y (3.3.3)
I(a) = Tc'e J_

90 CAPÍTULO 3. APLICACIONES

donde
•
-61/
- (x, a) = — (a)h [B (a) — x] 8[x — A(a)1 + B' (a)h — A(a)] 6[B (a) — x]
Da
= — A' (a)8[x — A(a)] + B' (a)5[B (a) — x].
(3.3.4)
El primer sumando en (3.3.4) se reduce utilizando la propiedad

h[B(a) — x]8[x — A(a)] = h[B(a) — A(a)]8[x — A(a)1 6[x — A(a)], (3.3.5)

ya que B (a) — A(a) > O y similarmente el segundo. Por lo tanto, usando la ecuación
(3.3.5) en la ecuación (3.3.3) e integrando con respecto a x, obtenemos la fórmula de
Euler-Leibnitz

1(a) f B(a)
A(a) L'a
(x, a) dx + B' (a) f [B (a), a] — A' (a) f [A(a), a]. (3.3.6)

Ahora, consideremos el problema más general de determinar la derivada respecto

a un parametro de una integral de volumen. Esto es,

1„ f
ut n(t)
Q(x ,t) dx, (3.3.7)

donde 51(t) es una región en el espacio la cual depende del tiempo y Q(x, t) una
función C'. Supongamos además que asi(t), la frontera de la región fk(t), es suave por

r
piezas y es descrita por la superficie f = O, con f > O en O. Supongamos también
que V f = ni, donde n' es el vector normal unitario hacia dentro de la superficie.
Por último, supongamos que el integrando es continuo respecto al tiempo. Entonces,
podemos reemplazar d/dt por cildt y llevar la diferenciación dentro de la integral,

ft, = —
_, f
d
Ud n(t)
h(f)Q(x,t) dx = I { a8t f 0Q
8(f)Q(x,t) + h(f)- ] dx
I*
(3.3.8)
= f
f
—L
a Q (x ,t) dS
9t
+ f 82
S1(t) U 5
dx,

donde h(f) es la función indicador y hemos usado (1.3.28) para integrar 8(f) en el
segundo paso en (3.3.8).
3.3. Algunas Aplicaciones Interesantes 91

Sabemos que
af
= = vn, (3.3.9)
at
y„
donde vn, y son las velocidades locales normales en la, dirección hacia dentro y hacia
afuera, respectivamente. Así,

I„ =ami) I vn 1Q(x,t) dS + I OQ
at
dx.
1-1(t)
(3.3.10)

Esta ecuación es una generalización de la fórmula de Euler-Leibnitz (3.3.6).

En el mismo contexto, también es posible calcular la derivada de la integral

(3.3.11)
=dt J Q(x,t) dS.

Para hacerlo, primeramente convertiremos la integral de superficie (3.3.11) en la inte-

gral de volumen,

= dt f 5(f)*(x,t) dx, (3.3.12)

n
donde f = O describe a 012(0, V f = es el vector normal unitario hacia afuera de
la superficie y Ej es la restricción de Q a la superficie f = O. Por lo tanto,

I 8 = f [--6 9
at (f)Q(x,t)+8(f)---91
at
d r. (3.3.13)

Usando (1.3.28) y (1.3.30), obtenemos

I, = fa {-2-
8ut — 2vn IC fQ(x,t)] dS, (3.3.14)
n(t)

80 aQ acj
donde ICT es la curvatura media de 81.2(t) y 79—
t = 791 + v,2797-1.

La deducción de la fórmula (3.3.14) por cualquier otro método no es nada trivial.

En la aplicación, las relaciones (3.3.10) y (3.3.14) son conocidas como Teoremas de
Transporte, y son utilizadas en la derivación de leyes de conservación.

92 CAPÍTULO 3. APLICACIONES

3.3.2 Teorema de la Divergencia Generalizado

El Teorema de la Divergencia es un resultado clásico de Cálculo Vectorial que
relaciona una integral de volumen con una integral de superficie de un campo vectorial
continuo y su derivada. Este resultado es utilizado en la derivación de algunas leyes de
conservación físicas en forma diferencial. Lo relevante es que podemos considerar este
resultado desde un mejor ángulo o punto de vista, (más generalizado); permanece válido
para campos vectoriales discontinuos o singulares lo cual implica que también tales leyes
de conservación permanecen válidas cuando todas las derivadas son consideradas como
derivadas generalizadas.

Sea S2 un volumen finito en el espacio y sea V(x) un campo vectorial C'(9) .

Definamos el campo vectorial discontinuo

V(x)
17(x) = S x 112 (3.3.15)
x

Sea f = 0 la superficie frontera an de la región 12 de tal forma que n = V f apunte

hacia afuera de 09 y 11V f = 1 sobre f = 0. Tenemos que

divV=V•V=V•V+LV1 • ná(f) = V • :5-7. —V•n8(f), (3.3.16)

usando el hecho que [Vi = 5- 7 (f = 0+ ) — V(f = 01 = —V (f = 0). Integrando sobre

todo el espacio, obtenemos

.
I f I dxi dx i dx2 dx3 = f vir -ce dx 2 dx3 = O.

Similarmente, obtenemos integrales iguales a cero para 8IT 2 /5x2 y 5173/3x3 , donde Vi
(3.3.17)

es la i-ésima componente del campo vectorial V. Por lo tanto,

fdiv V dx = 0. (3.3.18)

Por otro lado, la integración del lado derecho de (3.3.16), usando la relación (1.3.28),
definiendo V,2 = V n y usando el hecho que

V • V(x) x E2
V = { (3.3.19)
x ft
3.3. Algunas Aplicaciones Interesantes 93

da como resultado el Teorema de la Divergencia Clásico,

rz V • V dx — fon dS O. (3.3.20)
f
Ahora, notemos que (3.3.20) es válida si V tiene una discontinuidad sobre la
superficie a = O en sk,
Va= n

Fig.19 Superficie de discontinuidad del campo vectorial V.

Entonces, la ecuación (3.3.20) es válida si V • V en la integral de volumen se reemplaza
por V • V, donde el único salto de V en la divergencia generalizada proviene de la
discontinuidad en a = 0. Esto es, tenemos que
Vi •V=V•V+ yi • n'5(a), (3.3.21)

con ni = Va la normal unitaria a la superficie a = 0. Así, la ecuación (3.3.20) puede

ahora escribirse como
ir • V dx = f VdS, (3.3.22)
811

la cual a su vez, al usar (3.3.21), se puede escribir como

in V • V dx = f Vn dS — f [Vn,1 dS, (3.3.23)

s„
donde 1V,ill IV] • n' y So es la parte de la superficie a = O encerrada en la región
ft(Ver figura 19).


94 CAPÍTULO 3. APLICACIONES

3.3.3 Entre
• Espejos Te Veas

Un Problema interesante es encontrar la función de Green correspondiente a una

fuente localizada en medio de dos planos paralelos a la misma distancia (espejos), con
potencial específico en cada uno de ellos, como se muestra en la figura siguiente

Fig.20 Espejos paralelos encontrados.

Como se puede observar, determinar la fuente imagen correspondiente a alguno de los
planos es relativamente fácil puesto que esta satisface las condiciones sobre ese plano.
Sin embargo, no satisface las correspondientes al otro plano y por consiguiente se tiene
que considerar otra fuente imagen que las satisfaga, pero nuevamente encontramos que
ésta no va a satisfacer las condiciones del primero y así sucesivamente hasta tener una
infinidad de fuentes imagen. Por lo tanto, la función de Green es la suma infinita de las
contribuciones individuales de cada fuente imagen. Físicamente, el problema es como
si uno se mirara en dos espejos encontrados, en los cuales se venia nuestra imagen refle-
jada un número infinito de veces. Matemáticamente, el problema establecido significa
determinar la función de Green que satisface la ecuación
áG = —6(x — e)
sobre la banda (columna) infinita O < xl < L y —oo < xi < oo, i = 2, 3, ... y las
condiciones de frontera G = O en x1 = O y x1 = L. Para el caso 2-dimensional, por
construción la funcion de Green esta dada por
c°
G(x, e) ix — an l — lnIx — 13„)
= 2- n=—co
donde an = e — 2nL y fin = —e + 2nL, mientras que para tres dimensiones tenemos
que la funcion de Green asociada es
1 1
G (x, = —4)tr És' — a n lj
n=—co — Onll
donde an = — 2nD Y Ott = 2nD, con D (L, O, , O).
3.3. Algunas Aplicaciones Interesantes 95

3.3.4 Diferentes Representaciones para la Solución de la Ecuación de Onda

La solución de la ecuación de onda admite diferentes formas, dependiendo del
problema y del objetivo. A continuación se presenta su solución fundamental (función
de Green), la cual es un ejemplo para problemas no estacionarios. Después, se introduce
la solución en la forma de Poisson y de Kirchhoff, respectivamente, para concluir con
la generalización de la fórmula de Kirchhoff para superficies móviles.

Consideremos la ecuación de onda

2 w = Q(x,t) (3.3.24)

donde q 2 denota al operador de D'alembert, dado por q 2 tí —1 82 — A, sujeta a la

a 2 at2
8
condición inicial cp(z, O) f (x), — w(x, O) = g(x).
8t
Si el problema se analiza sobre una región finita o semi-infinita, entonces además
es necesario considerar condiciones de frontera no homogeneas. Para determinar la
solución de este problema introduciremos la función de Green G(x, t; r) como la
solución correspondiente a una fuente concentrada en x = 1, actuando instantáneamente
en t = T. Esto es,
2G = --8(x — Z)5(t — r).
La función de Green es la respuesta en x al tiempo t debido a una fuente localizada
en en el tiempo T. La respuesta es cero cuando la fuente actúa en t < T ("principio
de causalidad"); G(x,t; 1, T) = O para t < T.
Introduciendo la Transformada de Fourier y algunas de sus propiedades, tenemos que
a2b-
+ a2w2a = 5(t T),
at2 (270n

donde co2 = ca • ca. Por el principio de causalidad tenemos que

02a
=o para t > -r,
at2
cuya solución es de la forma

.= O Si t <

A cos aco(t — r) .B sen aw(t — si t >

96 CAPÍTULO 3. APLICACIONES

&ALI
De la continuidad tenemos que A = O y de la condición de salto B = Así, la
(27)naw
solución a está dada por

G(w, t; x 0 , to) sen aw(t — r).

(2/r)naco
Por lo tanto, la función de Green asociada a la ecuación de onda es
1 j" sen aw(t — r)
Gxt;
( , C , 7) = (27r)" dm
acv
En el caso de una dimensión sin considerar el efecto de fuente alguna, llegamos a
la conocida solución de D'Alembert,
1
cp(x,t) = — De (x + at) + f (x — at)] + j‘t-Fat g(e) de,
2 z_
siendo en este caso

G(x,t;e,r) — e) — a(t — 7)] + h[(x — e) + a(t — r)]).

la función de Green asociada.
En el caso tridimensional, considerando condiciones iniciales iguales a cero, se
obtiene
(p(x,t) = fe: f f f G(x,t; /-)Q(Z, 7) dedr

1 t 1
= 4car fo f ff 5[R — a(t — r)]Q(1,r) decir.

La funcion de Green para este caso específico tiene la expresión

G(x,t; 1, r) = 5[R— a(t — 7)].

4alt
Las únicas fuentes que contribuyen satisfacen lix — xoll = a(t — to). El efecto en
x al tiempo t se debe a todas las fuentes involucradas las cuales emiten ondas que
se propagan a velocidad a. A un tiempo t — r posterior, el efecto se localiza sobre
una cáscara esférica a una distancia a(t — r). Esto es parte de lo que se conoce como
"principio de Huygens" el cual solo es válido en tres dimensiones. En dos dimensiones
el problema no es trivial.
3.3. Algunas Aplicaciones Interesantes 97

3.3.5 La Fórmula de Kirchhoff

Consideremos primeramente el problema de Cauchy para la ecuación de onda
en tres dimensiones. Esto es, determinar 9 = cp(x, y, z, t) que satisfaga la ecuación
(3.3.24), Q = O y las condiciones iniciales

Wi t=0 = Wo(x , y , z),

(3.3.25)
at t=o cPi (x, y, z).

Una solución para este problema fue dada por Poisson, la cual se obtiene de la siguiente
forma: se considera una esfera S con centro en el punto p(1, y, 7) y radio variable
p = al, de tal forma que un punto sobre S tiene coordenadas

+ a/ sen cos O, y + al sen sen O, z + alcas 9),

donde O y 9 son dos de las coordenadas esféricas. El valor promedio Mp (9) de una
función sobre esta esfera está dado por

1
M°(co) = - 7r ff cp(± + al sen cos O, y + al sen sen O, z + al cos t,b) sen 9 MI
4
(3.3.26)

Entonces, la solución de la ecuación de onda en la forma de Poisson se escribe como

d
=— (3.3.27)
dl P((Po)] + IMp(Wi),
la cual proporciona el valor de la función c,o en el punto p en el momento l en términos
de sus valor inicial y el de su derivada con respecto al tiempo a una distancia p =
desde p.
Otra forma del problema de Cauchy fue obtenida por Kirchhoff, la cual se puede
determinar aún para casos más generales. Consideremos la ecuación de onda no
homogénea (3.3.24), donde Q se interpreta como la distribución de densidad de la
fuente. Para obtener el valor de yo, en cualquier punto p en términos de sus valores en
una región 9, se tiene que suponer que en esta región so y sus primeras derivadas son
continuas y que las segundas derivadas y Q son finitas e integrables. Denotando por R
la distancia pq, donde q es cualquier punto en la región 2, por ase la frontera cerrada
de 9 y por la diferenciación en la dirección de la normal exterior. La notación [1
n


98 CAPÍTULO 3. APLICACIONES

indica que el valor de la función se calcula en el tiempo t — Entonces, la fórmula de

Kirchhtff, para p en el interior de Q, está dada por
8cp aR
47rcpp = lif R
dV—ff
aa a ( R ) R L On J aR On at J
199Jan dS. (3.3.28)

Para p fuera de Q, la integral es cero. Para la ecuación de onda homogénea, la integral

de volumen (la cual representa al así llamado "'potencial retardado") no aparece en la
expresión. Además, observemos para este caso que wp depende solo de los valores de w y
sus drivadas en puntos q sobre la superficie frontera de 12 en un tiempo que antecede al
instante t por la cantidad 2, el cual es el tiempo necesario para que una perturbación
viaje de q a p con velocidad a. Así, la fórmula de Kirchhoff se considera como la
forma matemática general del Principio de Huygens, en el que q se considera como una
fuente secundaria que manda perturbaciones al punto p. Si suponemos que la región
12 es una bola de radio R = at, con centro en p, entonces la fórmula de Kirchhoff
8
se reduce a la fórmula de Poisson, observando que 1,o(t — = cp(0), — = y
8n 9R
Ocp 8cp
= —a- . Para una función periodica en el tiempo se obtiene la fórmula conocida
at 9R
como de Helmholtz. La obtención de la fórmula de Kirchhoff en base al análisis clásico
es laboriosa y compleja. Sin embargo, desde el punto de vista generalizado el trabajo
se reduce significativamente. En los últimos años de los 1980's Farassat y Myers [4]
dieron una moderna deducción de la fórmula de Kirchhoff para superficies móviles,
basándose en la teoría distribucional, la cual realmente es corta y evita la dificultad
de tratar con espacios y funciones de Green en cuatro dimensiones. Este resultado es
muy utilizado en aerodinámica, hidrodinámica, electrodinámica y acústica.
Supongamos que la superficie en movimiento sobre la cual se especifican las condi-
ciones de una función cp(x,t) se denota por a(x, t), la cual es una superficie deformable,
definida de tal forma que = 1 sobre esta. Supongamos que 99 satisface la ecuación
de onda en el exterior SZ de la superficie a = 0, la cual es la región definida por o- > 0.
Extendiendo ;o al espacio completo,
a>0
Cp(x,t) = {w(x,t) (3.3.29)
0 o- < 0
tenemos que (75 satisface la ecuación de onda sobre todo el espacio. Sin embargo, (7- ,
tiene discontinuidades sobre la superficie a = 0 las cuales aparecen como fuente en la
3.3. Algunas Aplicaciones Interesantes 99

ecuación de onda y pueden ser de diferente naturaleza. Por ejemplo,

Q(x,t) = S (x,t)
= S (x, t)8(o-)
(3.3.30)
= ái [S(x, t)15(al
=y [S(x, t)8(o-)] ,
donde S (x , t) puede ser una función ordinaria, una función con discontinuidades o bien
una función generalizada, dependiendo del tipo de singularidad.

Notemos que los saltos de 9- 5 y sus derivadas dependen de los valores correspondientes
aw
de w en la superficie ya que vp-i cp(a = 0+) y at1= — (0- = 0+). Aplicando las
at
reglas de la diferenciación generalizada (el operador de onda) a fi), tenemos

1 18
E2 ;2 = ( So n + —Mnsot) 8(c) — - -
at [mn925(0.)]– [9716(a)] (3.3.31)

donde y„ = -- aalet es la velocidad local normal de la superficie a = 0, Mn = vnice

ay, 899
el número local normal de Mach sobre a = 0, wt = — y cpn = — . Como siempre,
Ot 8n
suponemos que Va = n es la normal unitaria positiva a a = O. Los tres tipos de fuente
en el lado derecho de (3.3.24) son del tipo estandar descrito en (3.3.30). Aplicando
el método de las funciones de Green para el operador de onda sobre todo el espacio,
encontramos la fómula de Kirchhoff para movimiento subsónico, ya que debido a una
de las singularidades en (3.3.31), llamada de Doppler, no es aplicable a movimiento
supersónico. La solución para una superficie deformable dada por Farassat [4], es

CP(x,t) = { f (s) duid2+ fa,s, [7(7_2.v1

mr , L. duid2}, (3.3.32)
donde 17(8) es una región en el espacio u 1 , u2 que depende del tiempo sobre la cual
es transformada la superficie a = 0, g2 es el determinante de los coeficientes de la
primera forma fundamental y r* es el tiempo de emisión del punto (u i , u2 ) sobre la
superficie a = 0. La expresión E1 depende de (P, pn , V2W (gradiente sobre la superficie
de w) y de los parámetros geométricos y cinéticos de la superficie a = O. La expresión
E2 únicamente depende de los parámetros geométricos y cinéticos.

100 CAPÍTULO 3. APLICACIONES

En este preoblema se supone que la superficie cerrada a = O en movimiento es

impenetrable. Sin embargo, es posible que se tenga ruido (contaminación por sonido,
penetración de fluido, etc.) dentro de la superficie y como los principios de conservación
continuan siendo válidos en el sentido generalizado, entonces se sigue satisfaciendo la
ecuación de onda, solo que ahora este ruido aparece en los terminos no homogéneos
de una forma mucho muy complicada, bastante interesante como para estudiar con
detalle.

Nota: Las fórmulas de Kirchhoff resultan atractivas y de fundamental importancia

para la aeroacústica en el modelado de alta resolución aerodinámica y simulación bajo
regimen de turbulencia. Las matemáticas utilizadas en la derivación de estas fórmulas
son producto de la última década y ciertamente dignas de estudio.
3.3. Algunas Aplicaciones Interesantes 101

De esta forma termino mi exposición "rudimentaria" de la teoría de distribuciones

y algunas de sus aplicaciones, las cuales básicamente giran alrededor del concepto de
la diferenciación generalizada, con el enfoque de las funciones de Green. En realidad
el poder de la teoría de las funciones generalizadas proviene de sus propiedades opera-
cionales con las cuales es posible obtener fácilmente soluciones discontinuas (soluciones
generalizadas) de ecuaciones diferenciales lineales utilizando el método de las funciones
de Green, al poner el problema en el nuevo marco teórico. Además, la teoría de dis-
tribuciones de Schwartz unifica muchos métodos empíricos y da respuesta a preguntas
fundamentales de las ecuaciones diferenciales parciales lineales y tal vez con un poco de
cuidado también sea posible atacar problemas no lineales, siempre y cuando la multi-
plicación de distribuciones presente así lo permita. Esto es lo que alcanzo a vislumbrar,
pero sin duda lo más importante es que la teoría de distribuciones que nos heredo Lau-
rent Schwartz equipa a cualquier científico o ingeniero con un arma más para hacer
frente o resolver problemas de aplicación. En toda su extensión, el material de las
funciones generalizadas es muy abstracto y de un nivel matemático muy sofisticado, en
contrándose solamente en "excelentes" libros clásicos, difíciles de encontrar en nuestras
bibliotecas. En nuestro contexto, la abstracción se dejó de lado parcialmente ya que el
interés era el trabajar con la delta de Dirac multidimensional, sus propiedades y sus
derivadas con el fin de utilizarla junto con las funciones de Green. El marco teórico
utilizado fue el espacio de las distribuciones de Schwartz 79' por ser este el espacio
más fuerte. Sin embargo, cuando la teoría de distribuciones se anda sobre el conjunto
específico de funciones de prueba, S, el espacio de las distribuciones correspondiente re-
sulta especial para el estudio de la Transformada de Fourier en su versión generalizada,
con la cual es posible resolver aún muchos más problemas de ecuaciones diferenciales.
Espero sinceramente que se considere y se retome el presente trabajo como base para
un nuevo tema de tesis sobre la Transformada de Fourier Generalizada para de esta
forma completar el material de la teoría de las distribuciones.
102
CAPÍTULO 3. APLICACIONES
Bibliografía
ARFKEN G., "Mathematical Methods for Physicists",

Academic Press, 1985

COURANT H. & HILBERT D., "Methosd of Mathematical Physics",

Interscience, 1953 (Vol.I), 1962 (Vol.II)

FARASSAT F., "Introduction to Generalized Functions with Applications in Aerodyna-

mics and Aeroacoustics",
NASA Technical Paper, 1994

FARASSAT F. & Myers M.K., "Extension of Kirchhoff 's Formula to Radiation from
Moving Surfaces",

Journal of Sound and Vibration, Vol.123, No.3, 1998

FRIEDMAN B., "Lectures on Applications-Oriented Mathematics",

John Wiley Sons, Inc., 1991

GEL'FAND I. M. Shilov G.E., "Generalized Functions",

Vol.1—Properties and Operations,

Academic Press, 1964

GOETZ A., "Introduction to Differential Geornetry",

Addison — Wesley, 1970

HILDEBRAND F. B., "Advanced Calculus for Applications",

Prentice-Hall, Inc., 1976

103
104 BIBLIOGRAFÍA

HABERMAN R., "Elementary Applied Partial Differential Equations",

•
Prentice-Hall, 1983

RAUCH J., "Partid Differential Equations",

Springer—Verlag, 1991

RICHARDS I. & Youn H., "Theory of Distributions: A Non-Theoretical Introduction",

Cambridge University Press, 1990

SAICHEV A. I. & WOYCZYI(ISKI W. A., "Distributions in the Physical and Engineering

Sciencies", Vol.I-Distributional and Fractal Calculus, Integral Transforrns and Wavelets,
Birkhüser, 1997

STAKGOLD I., "Green's Functions and Boundary Value Problems",

John Wiley Sons, Inc., 1979

STEWART J., "Multivariable Calculus",

Brooks/Cole Puiblishing Company, 1991

[15] ZILL D.G., "Ecuaciones Diferenciales con Aplicaciones de Modelado",

Thompson Editores, 1997

3!..!OTECA
,)E1(4)15E`...Ç:ACTAS
YliATU 'MIES
PRI HIJOS
I , GRANDEZA