2 METODO DE FOURIER. Problema de Sturm Lioville
2 METODO DE FOURIER. Problema de Sturm Lioville
2 METODO DE FOURIER. Problema de Sturm Lioville
. (2.1.1)
(2.1.6)
. (2.1.7)
(2.1.8)
(2.1.9)
o sea
(2.1.15)
Los valores de son llamados valores característicos y siempre se obtienen como las raíces
positivas de la ecuación característica. Con lo anterior las soluciones (2.1.8) y (2.1.12) para
y X(x) son:
, (2.1.16)
, (2.1.17)
, (2.1.18)
El principio de superposición establece que la solución de una ecuación diferencial es la
combinación lineal de soluciones linealmente independientes, por lo tanto, la solución es
(2.1.19)
(2.1.20)
La ecuación anterior indica que puede ser expresada como una serie de Fourier de senos y
que los valores de las , n = 1,2,3,…, deben satisfacer dicha igualdad. Para determinar dichos
valores se aplica la propiedad de ortogonalidad de funciones (similar a la ortogonalidad de
vectores).
La ortogonalización se efectúa multiplicando toda la ecuación, por la misma función característica,
con diferente valor característico y se integra en el dominio (0 ≤ x ≤ l ).
(2.1.21)
Hay que notar, que en la anterior ecuación se define un vector de ecuaciones, una ecuación para
cada valor de m. Se puede mostrar por integración directa
De donde
En general
(2.1.22)
Al sustituir los valores de las constante An de la serie de Fourier de senos en la ecuación (2.1.19) se
tiene la solución completamente determinada para el problema planteado
(2.1.23)
Esta serie es convergente por lo que se puede truncar la sumatoria hasta un valor determinado de n,
dependiendo de la precisión deseada.
Ejemplo
Como caso particular se considerará una barra de hierro de 10 cm de longitud y difusividad
térmica a2 = 0.15 cgs sujeta a una temperatura inicial dada por la ecuación
T(x,0) = f(x) = 100, 0 < x < l.
Para simplificar la sumatoria se puede proceder a efectuar la integración
De donde
Al sustituir los valores de An, en la solución (2.1.23) se obtiene que la temperatura está dada por
Notar que solo se calcula para la mitad de la barra, ya que este físicamente es un problema
simétrico y se predice una distribución parabólica. Además ña convergencia de la serie debida al
rápido decremento del exponencial.
2.2 Problema de Sturm-Liouville.
(2.2.1)
, (2.2.2)
. (2.2.3)
Dependiendo del valor particular de las constantes A1, A2, B1 y B2, se obtienen condiciones de
primero, segundo y tercer tipo, homogéneas.
Sean yn(x) y ym(x) dos funciones características que se corresponden a dos valores característicos n
y m y que satisfacen al problema de Sturm-Liouville. Sustituyendo estas dos funciones y valores
característicos en la ecuación (2.2.1) se tiene las siguientes ecuaciones.
Multiplicando la primera ecuación por ym(x), la segunda por yn(x) y restando las ecuaciones
resultantes se obtiene
Integrando la anterior en el dominio, y cambiando el orden se tiene
(2.2.4)
Pero, las funciones características también deben satisfacer las condiciones frontera (2.2.5) y
(2.2.6), por lo cual se tiene
La, igualdad a cero de los determinantes anteriores se explica a causa de que un sistema de
ecuaciones homogéneo tiene solución si y solo si su determinante es igual a cero. Por lo tanto, la
ecuación (2.2.7) se reduce a
(2.2.8)
n m (2.2.9)
La ecuación anterior se le conoce como integral de ortogonalidad con función de peso, p(x).
Una generalización del concepto anterior se obtiene al considerar la expansión de una función
seccionalmente continua, f(x), en una serie generalizada de Fourier
, (2.2.10)
(2.2.11)
de acuerdo a la ecuación (2.2.9), todas las integrales del lado derecho son cero, a excepción de
una, cuando m = n, por lo tanto, la constante de combinación lineal es
(2.2.12)
(2.1.13)
hay que notar la semejanza entre la norma definida en la ecuación (2.2.13) y la norma para un
conjunto de vectores n-dimensional con componentes vk, dada por
. (2.2.14)
2.1.13
3. Problema que genera una Serie de Fourier de Cosenos
El siguiente ejemplo genera la serie de Fourier de Cosenos (Del libro de Smirnov)
Una sustancia disuelta con concentración inicial constante, C0, encerrada entre los planos x = 0 y x
= h, se empieza a difundir a t = 0, en un solvente localizado entre los planos x = h y x = l. Se
desea determinar el proceso de compensación de concentraciones, suponiendo que los planos x = 0
y x = l son impermeables a la sustancia.
C0 C=0
|-------------|----------------------------|
, (2.3.1)
, (2.3.2)
y CF (fronteras aisladas)
, (2.3.3)
. (2.3.4)
Matemáticamente la ecuación (2.3.1) es una EDP de segundo orden, homogénea, y las CF son de
segundo tipo homogéneas, con lo cual se tiene un problema homogéneo que puede ser resuelto
mediante el MÉTODO DE FOURIER.
Se propone solución
, (2.3.5)
con la cual, y de forma similar al procedimiento seguido en el problema anterior, se llega a las
soluciones
(2.3.6)
. (2.3.7)
Mediante la sustitución de (2.3.5) en (2.3.3) y (2.3.4) se obtiene
, (2.3.8)
. (2.3.9)
de donde, y/o .
Primero, si , la ecuación (2.3.7) es
, (2.3.10)
Con la condición (2.3.9) en la anterior se obtiene
(2.3.11)
(2.3.12)
y de (2.3.6) es
(2.3.13)
, (2.3.14)
.(2.3.16)
. (2.3.17)
(2.3.18)
y cuando
(2.3.17)
Hay que notar que la solución (2.3.17), es la suma de una solución en estado estacionario y una
solución transitoria, esta última tiende a cero cuando t tiende a infinito.
Solución de la ecuación de onda
. (2.3a.1)
(2.3a,3)
y condiciones frontera
, (2.3a.4)
. (2.3a.5)
(2.3a,7)
(2.3a,8)
(2.3ª,9)
(2.3ª,10)
(2.3ª.11)
La solución de /2.3ª,7) es
3a12
3ª13
Si se aplica la condición inicial 3ª,9 nos da
De donde
Si queda 3ª.14
Aplicando la condición 3ª.10 en 3ª,13
y
3ª,15
Aplicando la condición 3ª,11
De donde
En consecuencia
3ª.16
Ortogonalizando
Si m = n
Y la solución es
Cuando la ecuación de difusión (1.1.22) esta dada en coordenadas esféricas y el transporte solo
ocurre en dirección del radio se tiene
(2.4.1)
Lo problemas que se modelan con este tipo de ecuación pueden ser simplificados mediante la
transformación
(2.4.2)
para sustituir la anterior en la ecuación (2.4.1) es necesario obtener las siguientes derivadas
parciales
(2.4.3)
(2.4.4)
(2.4.5)
(2.4.6)
hay que notar que esta ecuación es idéntica a la ecuación de calor (o difusión) unidimensional en
coordenadas cartesianas.
Esferas de uranio de reciente formación, en las que ocurren reacciones nucleares, llevan a
cabo autogeneración de calor, Q, de densidad constante. Son enfriadas por medio de un fluido con
coeficiente de transferencia de calor h, de manera que el balance de calor para una sola partícula es
, (2.4.7)
, (2.4.9)
, (2.4.10)
que es equivalente a
. (2.4.11)
. (2.4.13)
La anterior ecuación puede ser separada en dos ecuaciones diferenciales, una para cada variable
dependiente.
. (2.4.14)
(2.4.15)
Físicamente, el término de generación de calor existe aun en el estado estacionario, lo que implica
su presencia en la ecuación para Ts(r). Por otra parte, matemáticamente, es justificable dejar el
término no homogéneo en la EDO, haciendo homogénea la EDP.
De igual manera, se elige que las condiciones frontera para el problema parcial sean homogéneas,
por lo luego de sustituir la definición (2.4.12) en las CF (2.4.10) y (2.4.9) se obtiene dos
condiciones para el estado transitorio
, (2.4.16)
, (2.4.17)
, (2.4.18)
, (2.4.19)
donde H = h/k.
Solución para la temperatura en estado estacionario. Tomando la EDO (2.4.14) e integrando
(2.4.20)
Al aplicar la condición (2.4.19) en la solución anterior, se llega a que C1 = 0, así que luego de otra
integración la temperatura en estado estacionario es
, (2.4.21)
(2.4.22)
Solución para la temperatura en estado transitorio. El problema descrito por las ecuaciones
(2.4.15), (2.4.16) y (2.4.17) es homogéneo, y se simplifica utilizando la definición (2.4.2) ,
, para dar
(2.4.23)
, (2.4.24)
, (2.4.25)
. (2.4.28)
con condiciones obtenidas a partir de (2.24) y (2.4.28),
, (2.4.29)
, (2.4.30)
, (2.4.32)
C3 = 0 nos llevaría a la solución trivial, por lo cual se llega a la ecuación característica
(2.4.33)
4R
3R
1R 2 R
R
, n = 1,2,3,… (2.4.34)
. n = 1,2,3,… (2.4.35)
. (2.4.36)
. (2.4.34)
Para encontrar los valores de las constantes de combinación lineal, An, es necesario aplicar la
condición inicial,
(2.4.36)
Ahora es necesario ortogonalizar, por lo que se multiplica la ecuación (2.4.36) por la ecuación
característica y la función de peso, y se integra en el dominio,
Para el caso en que m n de acuerdo al problema de Sturm-Louville, la integral del lado derecho
se vuelve cero, mientras que cuando n = m, esta integral toma la forma
obtiene que
(2.4.37)
La solución final es
(2.4.39)
(2.4.40)
2.5. Método de Separación de Variables en Varias Dimensiones
El método de separación de variables puede ser aplicado a ecuaciones diferenciales que
contiene más de dos variables independientes. En este tipo de problemas se generan dos o
más problemas de Sturm - Louville, así que para aplicar el método, es necesario que todas
las condiciones menos una sean homogéneas. A continuación, se presenta un ejemplo que
ilustra la manera de resolver este tipo de problemas.
(2.5.1)
Frontera x = 0: (2.5.2)
Frontera x = l: (2.5.3)
Frontera y = 0: (2.5.4)
Frontera y = h: (2.5.5)
Como se puede apreciar, este problema es homogéneo, así que su solución se puede obtener
mediante el método de separación de variables, en el cual se supone que la temperatura,
T(x,y,t), es el producto de tres funciones independientes entre si
(2.5.7)
Al sustituir esta definición en la ecuación (2.5.1) se tiene
Dividiendo la anterior por , y despejando el término que contiene x se logra
la primera separación
En este caso, como el término del lado derecho solo depende de x y el del lado izquierdo no,
ambos lados son independientes, por lo que se pueden igualar a una constante de separación
(2), por lo que a partir de esta ecuación se obtiene una ecuación diferencial ordinaria para
X(x), y otra ecuación que depende de y y de t.
(2.5.8)
(2.5.9)
Como se observa, el lado izquierdo de la ecuación (2.5.9) solo depende de t, mientras que el
lado derecho solo depende de y, así que ambos lados son independientes y se pueden igualar a
una constante de separación con lo cual se obtienen otras dos ecuaciones diferenciales
ordinarias:
(2.5.10)
(2.5.11)
(2.5.13)
(2.5.14)
(2.5.15)
(2.5.16)
las anteriores son definidas a causa de que las funciones X(x), Y(y) y deben ser diferentes
de cero para que la solución propuesta también lo sea.
Las condiciones (2.5.13) y (2.5.14) junto con la ecuación diferencial (2.5.8), al igual que las
condiciones (2.5.15) y (2.5.16) junto con la ecuación diferencial (2.5.11) forman dos problemas
de Sturm Louville cuyas funciones y valores característicos se obtuvieron en el primer
problema de este capítulo. Funciones y valores característicos que también pueden ser
obtenidos de la tabla 2.1, con base a estos:
(2.5.18)
(2.5.19)
(2.5.20)
(2.5.21)
Para tener una solución única, se utiliza el principio de superposición de soluciones, y debido a
que tanto n como m inician de uno hasta infinito, es necesario realizar dos sumatorias, una para
cada contador, es decir,
(2.5.22)
(2.5.23)
Debido a que se cuenta con una doble sumatoria, se debe ortogonalizar dos veces, una para x y
otra para y, por lo tanto es necesario multiplicar por las dos funciones características con otro
valor característico, e integrar dos veces, una en el dominio de x y otra en el dominio de y, es
decir
La primera integral del lado derecho es igual a cero cuando n p, mientras que sí n = p, es igual
a l/2. De igual manera, la segunda integral del lado derecho es igual a cero cuando m q,
mientras que sí m = q, es igual a h/2. Así que se eliminan todos los términos de la sumatoria
menos uno, para despejar así el valor de Anm
(2.5.24)
(2.5.25)
en donde w y v son variables mudas de integración.
, (3.1.1)
Esta Ecuación diferencial lineal homogénea es conocida como ecuación de Bessel y sus
soluciones son llamadas funciones de Bessel o funciones cilíndricas.
La mayoría de los problemas con valores en la frontera generan una ecuación de Bessel
con = n, n = 0,1,2,….., por lo cual es de interés encontrar la solución de la ecuación de
Bessel
x2 y’’(x) + x y´(x) + (x2 – n2) y(x) = 0, n = 0,1,2,…. (3.1.2)
Para encontrar el valor de la variable dependiente y(x), se utiliza el método de Frobenius
proponiendo una serie generalizada como la solución buscada
(3.1.3)
O efectuando la multiplicación
= -n y con estos valores forzosamente A1 = 0. Además, el coeficiente de x k+s para k 2 debe ser
igual a cero, o sea
(s – n + k)(s – n + k)Ak + Ak-2 = 0, k = 2,3,….
, k = 2,3,…. (3.1.4)
………
La sustitución sucesiva de las anteriores da
, k = 1,2,3,….
, k = 0,1,2,… (3.1.6)
(3.1.7)
(3.1.9)
(3.1.10)
(3.1.11)
La cual es una particularización de la ecuación
(3.1.12)
Además, J0(x) y J1(x) recuerdan las series de las funciones coseno y seno, las cuales son
, - < x <
, - < x <
J0(x) = 0, en x = J1(x) = 0, en x =
2.4048 3.8317
5.5201 7.0156
8.6537 10.1735
11.7915 13.3237
14.9309 16.4706
18.0711 19.6159
, (0) = 1 (3.1.16)
(3.2.1)
(3.2.2)
(3.2.3)
se observa que tiene la misma forma y se determina que la función de peso p() = . Así, la
integral de ortogonalidad para la función de Bessel esta dada por:
(3.2.4)
Dicho de otra manera, para ortogonalizar se requiere de multiplicar por , que es la
función de peso para las funciones de Bessel. Sin embargo también es necesario conocer el
valor de la integral (3.2.4) cuando j = k, es decir, que es necesario encontrar la norma de las
funciones de Bessel, ||y(jx)||, dada por la integral
(3.2.5)
Para encontrar el valor de esta integral se toma la ecuación (3.2.2) y se multiplica por
, y queda
(3.2.6)
se llega a
al sustituir los límites y despejar se obtiene la norma o magnitud de las funciones de Bessel
(3.2.7)
, 0<ρ<R (3.3.2)
(3.3.6)
(3.3.7)
(3.3.8)
(3.3.10)
Para que se mantenga finita la temperatura cuando → 0 (condición frontera (3.3.3) es
necesario que C3 = 0
(3.3.11)
(3.5) en (3.4) da T(R,t) = P(R)θ(t) = 0
P(R) = 0 (3.3.12)
Aplicando la anterior en (3.11) se tiene
(3.3.13)
de la anterior se obtiene la ecuación característica
J0(λR) = 0 (3.3.14)
Las raíces positivas de la anterior dan los valores característicos
λ1 = 2.4048/R
λ2 = 5.5201/R
, i = 1,2,3,... (3.3.16)
, i = 1,2,3,... (3.3.17)
Luego, de acuerdo a la solución (3.5) y aplicando el principio de superposición se obtiene
(3.3.18)
(3.3.19)
Serie de Fourier-Bessel
Ortogonalizando la anterior se tiene
, j = 1,2,3,.... (3.3.20)
i=j
Para integrar la integral del lado izquierdo de la anterior, de tablas (Del libro de Spiegel Manual
de Formulas y tablas matemáticas o del libro de Operational Mathematics de Churchill) se tiene
la identidad
Para la segunda integral del lado izquierdo también se recurre a tablas de donde
por lo tanto
y si i = j de (3.2.7)
Ya que
de donde
, i = 1,2,3,....
(3.3.21)
, 0<ρ<R (3.2)
C(ρ,t) < M, t>0 (3.3)
C(R,t) = C1, t>0 (3.4)
Se propone que la solución sea la suma de dos soluciones una en estado estacionario y otra
transitoria
(3.5)
, (3.6)
, t>0 (3.7)
, t>0 (3.8)
Por su similitud con el problema anterior se puede tomar la solución (3.3.18) como la solución
para nuestra variable
(3.9)
(3.10)
Aplicando la condición inicial (3.2) se tiene
, j = 1,2,3,....
Aplicando la ortogonalización se obtiene
Al sustituir en (3.10)
Funciones de Legendre.
Las soluciones de la ecuación diferencial
(1)
o su equivalente
(2)
Son conocidas como funciones de Legendre de orden p, donde p se supone real y positiva.
Estas aparecen en problemas con fronteras esféricas, cuando existe simetría rotacional.
Se propone solución por el método de Frobenius
(3)
Sustituyendo en la ecuación diferencial se tiene
Y de los coeficientes de
Si s = 0
k=1
k=2
Por inducción se determina que
Si s =1
k =1
Por inducción
Lueg
o se obtienen solución
(4)
(6)
Y si n es non
(7)
Además, en x = 1
(8)
(9)
(10)
(13)
.
Al comparar la ecuación diferencial (2) se determina que
, , (función de peso) y
se determina que la solución de (1) que es finita en el intervalo existe solo si p es
entero positivo o cero, que es , y se corresponde a . De acuerdo a
la teoría del problema de Sturm Liouville la integral de ortogonalidad consiste en
(14)
y alguna función de x puede ser representada como una combinación lineal de polinomios de
Legendre
(15)
Ortogonalizando se obtiene
Para obtener el valor de la norma o magnitud de los polinomios de Legendre se puede utilizar la
fórmula de recurrencia (12)
(16)
Al multiplicar toda la anterior por y ortogonalizando se obtiene
(17)
Donde el paréntesis se utiliza para indicar la integral de las funciones que abarca
(19)
de (18) y (19) se obtiene
Y el conjunto
es ortonormal en el intervalo .
y es independiente de .
satisface la ecuación de Laplace en coordenadas esféricas
, (1)
(2)
puede representar la temperatura al interior de una esfera cuya superficie tiene una
temperatura que se puede expresar como una función de , o puede representar el potencial
electrostático dentro de una esfera libre de cargas, cuando
Si se define, , (3)
Para sustituir en el segundo termino de (1)
Sustituyendo
Con lo cual la ecuación diferencial (1) es
(4)
(6)
Sea
(7)
Sustituyendo en (4) se tiene
(8)
(9)
(11)
Para la e.d. (8) se tiene
(12)
(16)
Ortogonalizando
Si m = n
Recordando la definición de
-------------------------------------------------------------------------------------------------------------------
La base de un hemisferio sólido se encuentra aislada. El flujo
o con
(1)
De donde
(2)
(3)
(4)
(5)
(6)
ortogonalizado
Si m = n
INTEGRALES DE FOURIER
Transferencia de calor en una región semi infinita
(1)
(2)
(3)
(4)
Sea
(5)
Aplicando (5) en (1), (3) y (4) se obtiene
(6)
(7)
,
donde , es continua a lo largo del dominio x > 0
(8)
(10)
En la ultima integral se invirtió el orden de integración
Por ortogonalidad si
y si
En consecuencia de (10) se obtiene
(11)
(12)