Docstoc

Derivadas_parciales

Document Sample
Derivadas_parciales Powered By Docstoc
					                            Ecuaciones en Derivadas Parciales.

1. Introducción.
Una ecuación diferencial en derivadas parciales (PDE), por su semejanza con las ODE, es una
ecuación donde una cierta función incógnita u viene definida por una relación entre sus derivadas
parciales con respecto a las variables independientes.

Si u = u(x,y,z), una ecuación diferencial en derivadas parciales sería

                               ∂u ∂u ∂u ∂ 2 u ∂ 2 u ∂ 2 u ∂ 2 u 
                 F x, y, z, u , , , , 2 , 2 , 2 ,              ,... = 0
                               ∂x ∂y ∂z ∂x ∂y ∂z ∂x∂y 

Se denomina orden de la PDE al más alto grado de derivación parcial que aparece en la expresión.

Así

                 ∂2 u ∂2 u
                     +     =0                                                            (1)
                 ∂x 2 ∂y 2

es una PDE de 2 orden, mientras que

                     ∂u ∂u
                u·     +   =0                                                            (2)
                     ∂x ∂y

es una PDE de primer orden.

La ecuación (1) es lineal ya que u y sus derivadas aparecen sin multiplicarse y no aparecen
elevadas a potencias. La ecuación (2) es, en cambio, no lineal.

Se podrían hacer algunas consideraciones acerca de las PDE. Mientras que el número de
constantes que hay que eliminar en una familia de curvas definen el orden de la ecuación
diferencial ordinaria de la que es solución, aquí la génesis se puede considerar vista de otro modo.

Sea una función de (x,y) que verifica que es de la forma

                u(x,y) = f(x + y)·g(x - y)

donde f y g son funciones arbitrarias. Entonces

                 ∂u
                    = f ′ · g + g′ · f
                 ∂x
                 ∂u
                    = f ′ · g − g′ · f
                 ∂y

a su vez

                 ∂2 u
                      = f ′′ · g + 2 · f ′ · g′ + f · g′′
                 ∂x 2
                ∂2 u
                     = f ′′ · g − 2 · f ′ · g′ + f · g ′′
                ∂y 2

                 ∂ 2u
                      = f ′′ · g − f ′ · g ′ − f · g′′ + f ′ · g ′ = f ′′ · g − f · g′′
                ∂x∂y

de donde se deduce que

                ∂2 u ∂2 u
                    −     = 4 · f ′ · g′
                ∂x 2 ∂y 2

Pero

                ∂u ∂u
                  +   = 2· f ′·g
                ∂x ∂y

                ∂u ∂u
                  −   = 2 · f · g′
                ∂x ∂y
                                     2
                 ∂u 
                       2
                          ∂u 
                  −   = 4 • f ′ • g′ • u
                 ∂x     ∂y 

                                                                2
                 ∂ 2 u ∂2 u  ∂u     ∂u 
                                     2

                 2 − 2 •u=   − 
                 ∂x    ∂y    ∂x     ∂y 

Se ha llegado a una ecuación diferencial de 2° orden en derivadas parciales cuya solución tiene la
forma

                u(x,y) = f(x + y)·g(x - y)

Considérese otro ejemplo:

                u = x·f(y)

                ∂u
                   = f (x )
                ∂x
                         ∂u
                u = x·
                         ∂x
es la ecuación diferencial de primer orden cuya solución tiene la forma

                u = x·f(y)

donde f es una función arbitraria.

Otro ejemplo es el siguiente:

                u = f (x2 + y2 )




                                                                                                2
                ∂u
                   = 2· x·f ′
                ∂x
                ∂u
                   = 2· y·f ′
                ∂y

                     ∂u      ∂u
                y·      − x·    =0
                     ∂x      ∂y

Otro ejemplo:

                u = f(x + y) + g(x - y)

                ∂u
                   = f ′ + g′
                ∂x
                ∂u
                   = f ′ − g′
                ∂y

                ∂2 u
                     = f ′′ + g ′′
                ∂x 2

                ∂2 u
                     = f ′′ + g ′′
                ∂y 2

                ∂2 u ∂ 2 u
                    =
                ∂x 2 ∂y2

Repárese que según los resutados obtenidos existen infinitas soluciones posibles de la PDE. Pero
ahora la arbitrariedad de la solución general viene dada en términos de funciones, apareciendo
tantas como el orden de la ecuación.

Desde el punto de vista de la Matemática puede parecer más preciso obtener en cualquier caso la
solución general, sin embargo, se van a buscar soluciones dentro del campo de la Física por lo que
sólo interesará una solución particular concreta. Estas soluciones particulares van a satisfacer unas
determinadas condiciones de contorno y de valor inicial.

Es decir, se va a tratar de obtener la solución de una cierta PDE que verifique unas condiciones en
el contorno del dominio en que está definida (condiciones de contorno), y si además una variable es
el tiempo "t" las condiciones en t = 0 se darán como dato (condiciones iniciales).

Por último, y por lo que respecta a la clasificación, cuando cada término de la ecuación diferencial
contiene la función o sus derivadas esta ecuación se dice homogénea.

Algunos ejemplos típicos de ecuaciones en derivadas parciales son:

Ecuación de difusión:

                ∂u       ∂2 u
                   = c2 · 2
                ∂t       ∂x




                                                                                                   3
Es la clásica ecuación unidimensional de difusión del calor, de segundo orden, lineal, homogénea y
de coeficientes constantes.

Ecuación de onda:

                ∂2u        ∂2u
                     = c2 · 2
                ∂t 2       ∂x

Es la clásica ecuación de onda unidimensional, que describe fenómenos de tipo oscilatorios y es
también de segundo orden, lineal, homogénea y de coeficientes constantes.

Ecuación de Laplace:

                ∂2 u ∂2 u
                    +     =0
                ∂x 2 ∂y 2

Esta es una ecuación bidimensional, de segundo orden, lineal, homogénea y de coeficientes
constantes, describiendo potenciales eléctricos o gravitatorios o procesos de difusión en los que se
ha alcanzado un equilibrio térmico.

Ecuación de Poisson:

                ∂2 u ∂2 u
                    +     = f ( x, y )
                ∂x 2 ∂y 2

Es también una ecuación bidimensional, de segundo orden, lineal, de coeficientes constantes, pero
no homogénea.

Este curso se va a centrar exclusivamente en el estudio de las ecuaciones diferenciales de 2
orden lineales con coeficientes constantes, que son las más habituales en distintos campos de la
física.




2. ECUACIONES DE SEGUNDO                      ORDEN        LINEALES       HOMOGÉNEAS            DE
COEFICIENTES CONSTANTES
Una ecuación de 2 orden en derivadas parciales lineal, homogénea, con coeficientes constantes
tiene la forma (supuesto dos variables independientes):

                  ∂2u        ∂2u     ∂2u        ∂u       ∂u
                a· 2 + 2·h·      + b· 2 + 2·f ·    +2·g·    +c· u = 0
                  ∂x        ∂x∂y     ∂y         ∂x       ∂y

donde

a,h,b,f,g,c son constantes. Por comparación con una cónica

                       a        h f   x
                                     
                (x y 1) h       b g  y = 0
                       f
                                g c   1
                                      




                                                                                                  4
                  a · x2 + b · y2 + 2 · h · x · y + 2 · f · x + 2 · g · y + c = 0

se puede decir que estas ecuaciones se clasifican en elípticas, parabólicas e hiperbólicas de igual
modo que las cónicas. Esto es, si

                  a · b − h 2 > 0 la ecuación es elíptica

                  a · b − h 2 = 0 la ecuación es parabólica

                  a · b − h 2 < 0 la ecuación es hiperbólica

Según esto, las clá sicas ecuaciones de difusión, de ondas y de Laplace pertenecen a los tipos

Ecuación de difusión: parabólica

Ecuación de onda: hiperbólica

Ecuación de Laplace: elíptica
Nota: Esta clasificación sigue siendo válida incluso cuando los coeficientes de la ecuación
a, b, h, f, g, c so funciones variables de x e y. En estos casos la ecuación puede cambiar de
tipo al pasar de un cuadrante a otro. Por ejemplo la ecuación

                     ∂2 u         ∂2u     ∂2 u
                  y · 2 + 2· x ·      + y· 2 = 0
                     ∂x          ∂x∂y     ∂y
es elíptica en la región y 2 − x 2 > 0, parabólica a lo largo de las rectas y 2 − x 2 = 0, e
hiperbólica en la región y 2 − x 2 < 0.
En el caso de más variables independientes la forma general de una ecuación diferencial lineal de
segundo orden en derivadas parciales es:
 n   n                           ∂ 2u     n                               ∂u
∑ ∑ ark ( x1, x 2,......x n ) ⋅ ∂x      + ∑ b k ( x1, x 2 ,......x n ) ⋅
                                    ∂x k k =1                            ∂x k
                                                                              + c( x 1, x 2 ,......x n ) ⋅ u = f ( x1, x 2 ,......xn )
r =1 k =1                         r


Se denomina parte principal de la ecuación diferencial al primer sumando simbólico. Considérese la
matriz cuadrada nxn cosntituida por los coefiecientes de la parte principal de la ecuación
diferencial.
Por ser una matriz simétrica en el campo real, tiene autovalores reales.
    a) Si todos los autovalores fueran del mismo signo (ninguno nulo), la ecuación se denomina
        elíptica.
    b) Si un autovalor fuera de signo opuesto a los otros, no siendo nulo ninguno, la ecuación es
        hiperbólica.
    c) Si algún autovalor es nulo, se denomina parabólica.
    d) En el resto de casos se denomina ultrahiperbólica.

3. ECUACIONES DE EULER
Se llama ecuación de Euler a una ecuación de la forma

                       ∂2 u           ∂2u     ∂2u
                  a·        + 2· h ·      + b· 2 = 0
                       ∂x 2          ∂x∂y     ∂y

La solución general se puede obtener del siguiente modo, haciendo el cambio




                                                                                                                               5
                                    ∂α         ∂α
                α = p· x+ q ·y         =p         =q
                                    ∂x         ∂y

                                   ∂β        ∂β
                β = r · x + s· y      =r        =s
                                   ∂x        ∂y

donde p,q,r y s son constantes

                ∂u ∂u ∂α ∂u ∂β      ∂u      ∂u
                  =  ·  +  ·   = p·    + r·
                ∂x ∂α ∂x ∂β ∂x      ∂α      ∂β

                ∂2 u    ∂  ∂u  ∂α ∂  ∂u  ∂β
                   2 =     •     +    •      =
                ∂x     ∂α  ∂x  ∂x ∂β  ∂x  ∂x

                       ∂2 u        ∂2 u            ∂ 2u    ∂2 u 
                = p • p • 2 + r •        + r • p •      +r• 2
                       ∂α         ∂α∂β          ∂α∂β       ∂β 

                         ∂2u              ∂2 u       ∂2 u
                = p2 ·        + 2· r ·p·       + r2 · 2
                         ∂α 2            ∂α∂β        ∂β

de igual modo

                ∂2 u       ∂2 u             ∂2 u       ∂2 u
                     = q2 · 2 + 2 · q · s·       + s2 · 2
                ∂y 2       ∂α              ∂α∂β        ∂β

por último

                 ∂2u   ∂  ∂u   ∂  ∂u  ∂α ∂  ∂u  ∂β
                     =     =      •     +  •
                ∂x ∂y ∂y  ∂x  ∂α  ∂x  ∂y ∂β  ∂x  ∂y

                     ∂     ∂u    ∂u    ∂      ∂u      ∂u 
                =      p •    +r•  •q+    p •    + r •  •s =
                    ∂α  ∂α       ∂β    ∂β  ∂α         ∂β 

                         ∂2 u                     ∂2u         ∂2 u
                = p · q · 2 + ( q · r + p · s) ·      + r · s· 2
                         ∂α                      ∂α∂β         ∂β

Sustituyendo en la ecuación diferencial

                     2 ∂ 2u               ∂2u     ∂2u 
                a • p •      + 2• r •p •      +r • 2 +
                                                 2

                        ∂α 2             ∂α∂β     ∂β 

                                  ∂2u                        ∂ 2u        ∂2 u 
                + 2 • h • p • q •      + ( q • r + p • s) •       + r •s• 2  +
                                  ∂α 2                      ∂α∂β         ∂β 

                           ∂ 2u              ∂2u       ∂2u 
                + b • q 2 • 2 + 2 • q • s •      + s2 • 2  = 0
                           ∂α               ∂α∂β       ∂β 




                                                                                   6
                 ∂2 u                                         ∂2u
                      • (a • p 2 + 2 • h • p • q + b • q 2 ) + 2 • (a • r 2 + 2 • h • r • s + b • s 2 ) +
                 ∂α 2                                         ∂β

                                                                           ∂2 u
                 + 2 • ( a • p • r + h • q • r + h • p • s + b • q • s)         =0              (3)
                                                                          ∂α∂β

Ahora se elige p = r = 1 de modo que q y s sean raíces de la ecuación

                 a + 2 · h · x + b · x2 = 0

                                                    ∂2 u ∂ u
                                                               2
es decir, de modo que los coeficientes de la            y     sean cero.
                                                    ∂α 2 ∂β
                                                            2



Por tanto, llamando a las raíces x 1, x 2 quedaría la ecuación:

                                                         ∂2 u
                  a + h · ( x1 + x 2 ) + b · x1 · x 2         =0
                                                        ∂α∂β

Ahora bien

                              −2 · h                                  a
                 x1 + x2 =                     y         x1 · x 2 =
                               b                                      b
por lo que la ecuación puede expresarse

                                     ∂2u
                   • (a • b − h ) •
                 2             2
                                         =0
                 b                  ∂α∂β

Si a · b − h2 ≠ 0 , es decir la ecuación es elíptica o hiperbólica ⇒

                  ∂2u
                      =0
                 ∂α∂β

cuya solución general se reduce a

                u = F(α) + G(β)

donde F y G son funciones arbitrarias, pero

                 α = x + x1 · y

                β = x + x2 · y

luego la solución general de las ecuaciones elípticas e hiperbólicas es de la forma:

                 u = F( x + x1 • y ) + G( x + x 2 • y )

x 1 y x 2 son reales si la ecuación es hiperbólica, pero si es elíptica, son complejas.




                                                                                                            7
Si la ecuación es parabólica:

                  a · b − h2 = 0

volviendo a la ecuación (3) y haciendo sólo p = 1

                  ∂2 u                               ∂2 u
                       • (a + 2 • h • q + b • q 2 ) + 2 • (a • r 2 + 2 • h • r • s + b • s2 ) +
                  ∂α 2                               ∂β

                                                                    ∂2u
                  + 2 • ( a • r + h • q • r + h • s + b • q • s)        =0            (4)
                                                                   ∂α∂β

Se busca q tal que

                  a + 2 · h · q + b · q2 = 0

                       − h ± h2 − a · c    h
                  q=                    = − (raíz doble)
                             b             b

Llevando este valor a (4)

                                           ∂2 u               h2 • r             ∂2 u
        (a • r + 2 • h • r • s + b • s ) • ∂β 2 + 2 •  a • r + b + h • s − h • s ∂α∂β = 0
              2                         2
                                                                                

pero como h 2 = ab, sustituyendo el valor de "a" despejado de esta igualdad, la ecuación queda

                   h2 • r 2                       ∂ 2u
                            + 2 •h • r •s+ b •s  • 2 = 0
                                                2
                   b                              ∂β

                  1                       ∂ 2u
                    • ( h • r + b • s) 2 • 2 = 0
                  b                       ∂β

con tal que r y s no sean cero simultáneamente, la ecuación resultante es

                  ∂2 u
                       =0
                  ∂β 2

cuya solución general es de la forma

                  u = β · F(α ) + G ( α )

con F y G funciones arbitrarias, pero

                           h
                  α = x−     ·y
                           b
                  β = r · x +s·y

con r y s arbitrarios, pero no simultáneamente ceros.




                                                                                                  8
Luego, la solución general de una ecuación parabólica es

                                             h            h 
                u = ( r • x + s • y ) • F x − • y + G  x − • y
                                             b            b 

Aunque se ha resuelto, desde el punto de vista matemático, las ecuaciones de Euler, estas
soluciones tienen muy poco valor cuando se imponen unas condiciones de contorno dadas y una
condición de valor inicial. Suele resultar muy difícil obtener la expresión de las funciones F y G.
Por ello, este procedimiento más académico que útil va a dar paso a otro más eficaz que,
además, nos va a ayudar a ver el sentido físico de lo que se trata de resolver. Este método se
conoce con el nombre de método de separación de variables.




                                                                                                  9

				
DOCUMENT INFO
Shared By:
Categories:
Tags:
Stats:
views:1258
posted:5/31/2010
language:Spanish
pages:9