Mate
Mate
Trabajo Grupal
Ecuaciones diferenciales parciales (EDP)
11 de julio de 2018
1
MATEMÁTICA SUPERIOR
1.1.1. Ejemplo
A continuación se presentan una serie de ecuaciones importantes:
El orden de una ecuación diferencial parcial, es el orden de la mayor derivada que aparece
en la ecuación.
1.1.2. Ejemplo
La ecuación ut + uxxx = 0, es de tercer orden; la ecuación ut + uux = 0, es de primer orden;
la ecuación ut + uxx = 0, es de segundo orden.
Se recuerda del curso de ecuaciones diferenciales ordinarias, que dada una ecuación diferen-
cial se le puede asociar un operador L, el cual está definido sobre un espacio de funciones, por
ejemplo, con la ecuación:
dl2 u(t) du(t)
2
+3 − 5u(t) = t2 + 1
dt dt
se asocia el operador
2
MATEMÁTICA SUPERIOR
1.1.3. Ejemplo
∂u ∂ 2 u ∂u ∂ 2 u
Con la ecuación, − 2 = cos x, se asocia el operador L[u] = − 2
∂t ∂x ∂t ∂x
∂ 2u ∂ 2u ∂ 2u ∂ 2u ∂ 2u
Con la ecuación de Laplace, 2 + 2 + 2 = 0, se asocia el operadorL[u] = + +
∂x ∂y ∂z ∂x2 ∂y 2
∂ 2u
.
∂z 2
La ecuación de Burgers ut + uux = 0 tiene asociado el operador no lineal L[u] = ut + uux .
Es decir, en forma análoga a las ecuaciones diferenciales ordinarias, las ecuaciones en
derivadas parciales, se dicen lineales, si el operador asociado es lineal.
I = B 2 − 4AC.
Una vez que se tiene la invariabilidad del signo del discriminante I bajo el cambio de coor-
denadas, tiene sentido clasificar de acuerdo con tal signo a las ecuaciones en derivadas parciales
lineales de segundo orden con coeficientes constantes. Se tiene la siguiente definición:
Se dice que la ecuación (1.2) es:
Parabólica, si I = 0
Elı́ptica, si I ¡0
Hiperbólica, si I ¿0
3
MATEMÁTICA SUPERIOR
0
x = xy − z, x(0) = cosh s,
y 0 = y 2 − 1, y(0) = 0,
0
z = yz − x, z(0) = senh s.
Con lo cual, la solución viene dada en forma implı́cita por la expresión del hiperboloide
parabólico de una hoja
x2 + y 2 − z 2 = 1,
4
MATEMÁTICA SUPERIOR
Dado que dtd u(x(t), y(t)) = 0 se sabe del cálculo diferencial que u(x(t), y(t)) = α, donde
α es una constante. Cabe aclarar que posiblemente se requiera una constante diferente para
diferentes valores x0 , y0 . Resumiendo, la función u es constante sobre cada recta paralela a la
ecuación cartesiana de la recta parametrizada es bx − ay = c, donde c = bx0 − ay0 por lo que
u es constante cuando bx-ay es constante, es decir,
1.4.1. Ejemplo
Encuentre la solución de la ecuación 2ux + uy = 0, la cual satisface la condición inicial
u(0, y) = cosy
Solución. Usamos directamente la fórmula anterior con lo cual tenemos u(x, y) = f (x − 2y)
donde f es una función por determinar a partir de la condición inicial u(0, y) = cosy. Tenemos
u(0, y) = f (−2y) = cosy. Si ponemos w = −2y, entonces y = −w/2 y ası́ f (w) = cos(−w/2) =
5
MATEMÁTICA SUPERIOR
cos(w/2), dado que la función coseno es función par. De donde se obtiene la solución.
x − 2y
u(x, y) = f (x − 2y) = cos( )
2
El lector debe verificar derivando directamente que la u obtenida satisface la ecuación 2ux +uy =
0 y también que u satisface la condición inicial u(0, y) = cosy
1.5.1. Ejemplo
Resuelva el problema
ux − yuy = 0,
u(0, y) = y 3 ,
Solución. Las caracterı́sticas de la ecuación están dadas por el sistema
0
x (t) = 1,
0
y (t) = −y,
el cual puede ser integrado fácilmente para obtener
x(t) = t + c1 ,
y(t) = c2 e−t ,
6
MATEMÁTICA SUPERIOR
De esta forma puede definirse G(x, y) = yex = c como las curvas caracterı́sticas de la ecua-
ción.Por lo tanto,
u(x, y) = f (yex ),
donde f es una función por determinar a partir de la condición u(0, y) = y 3 , es decir f (y) = y 3 .
Finalmente se llega a la solución del problema: u(x, y) = (yex )3 = y 3 e3x .
1.5.2. Ejemplo
Encuentre las solución del problema
3
yx ux + uy = 0, x 6= 0
u(x, 0) = cosx,
Solución. En este caso la ecuación cartesiana de las caracterı́sticas puede encontrarse directa-
mente integrando la ecuación
dx
= yx3
dy
Tenemos entonces separando variables e integrando 0 − x−2 = y 2 + c o bien
1
−c = + y2
x2
Por lo tanto, las soluciones de la edp en el problema son de la forma
u(x, y) = f (x−2 + y 2 )
donde f se debe determinar a partir de la condición u(x, 0) = cosx. Claramente, para x 6= 0, la
ecuación y la condición dada llevan a
cosx = f (x−2 )
Y al definir w = x−2 se concluye que f (w) = cos(1/w1/2 ) y por lo tanto la solución del problema
es
7
MATEMÁTICA SUPERIOR
curva parametrizada
∂u ∂u
(a, b, c).( , , 1) = 0
∂x ∂y
El lector debe recordar de sus cursos de cálculo que para una superficie dada por la gráfica
∂u ∂u
dez = u(x, y),el vector( , , 1)representa una normal a la superficie. De esta manera el vec-
∂x ∂y
tor (a,b, c) se encuentra en el plano tangente a una superficie solución en cada punto. Por lo
tanto se puede construir una curva contenida en la superficie solución al resolver el sistema
1.6.1. Ejemplo
Resuelva el problema
ux + uy = x + y, (2,22)
u(x, 0) = x2
Solución:En este caso, de la edp tenemos que dx/dt = 1, dy/dt = 1,du/dt = x + y y la condición
inicial x(s, 0) = s,y(s, 0) = 0,u(s, 0) = s2 . Al integrar se obtiene x = t + c1(s),y = t + c2 (s).De
la condición inicial t = 0 se llega a x = s+t, y = t,u = t2 + st + s2 . En este caso es posible
escribir u como función de x, y para obtener
u(x, y) = x2 ?xy + y 2 .
El lector puede comprobar que u satisface la edp y la condición inicial.
8
MATEMÁTICA SUPERIOR
ut + g(u)ux = 0
puede ser resuelta por el método de las caracterı́sticas estudiado en este capı́tulo. Para esta
ecuación las caracterı́sticas están dadas por
dt(s, τ )/dτ = 1,
Ası́ u(s, τ ) = c1, τ (s, τ ) = τ + c2 y x(s, τ ) = g(c1)τ + c3. Si se dan las condiciones iniciales
para τ = 0: x(s,0) = x0(s), t(s,0) = t0(s), y u(s,0) = u0(s). Se obtiene la solución.
1.7.1. Ejemplo
Resuelva el problema
ut + uux = 0
0, x < 0,
u(x, 0) =
1, x > 0.
Solución. Para este problema la condición inicial (2.29) puede ser representada paramétrica-
mente como
t(s, 0) = 0, x(s, 0) = s, u(s, 0) = u(x(s, 0), 0).
Para las caracterı́sticas tenemos el sistema:
dt(s, τ )/dτ = 1,
9
MATEMÁTICA SUPERIOR
Al integrar se obtiene
t(s, τ ) = τ,
ut + uux = 0
0, s < 0,
u(s, 0) =
1, s > 0.
0, s < 0,
x(s, τ ) = s + τ
1, s > 0.
Las caracterı́sticas pueden construirse ahora partiendo de las condiciones iniciales por ejemplo
la caracterı́stica que pasa por el punto (−2, 0) es x = −2 y la caracterı́stica que pasa por (2,0)
es x = 2+t. Entonces tenemos u(−2, t) = 0 y u(2+t,t) = 1, para t > 0. De esta manera se ha
obtenido la solución del problema en forma paramétrica.
De la figura 2.3 el lector puede darse cuenta que no hay solución dada por este método en la
cuña 0 < x < t. Esta situación inquietante puede resolverse adentrándose en el estudio de las
ondas de choque (shocks) y puede consultarse por ejemplo en el artı́culo clásico de Lax [21]. El
siguiente ejemplo tomado del mismo artı́culo de Lax da la pauta para el estudio de soluciones
generalizadas, las cuales pueden ser discontinuas, pero cuyo estudio va más allá de los objetivos
de este libro.
Definicion
La ecuación del calor es una importante ecuación diferencial en derivadas parciales del tipo
parabólica que describe la distribución del calor (o variaciones de la temperatura) en una re-
gión a lo largo del transcurso del tiempo. Para el caso de la variable temporal t, y de una
función de tres variables
(x, y, z) ∈ R3
el enunciado clásico es:
∂u ∂ 2u ∂ 2u ∂ 2u
− α( 2 + 2 + 2 ) = 0
∂t ∂x ∂y ∂z
10
MATEMÁTICA SUPERIOR
La ecuación de difusión, es una versión más general de la ecuación del calor, y se relaciona
principalmente con el estudio de procesos de difusión quı́mica.
En el caso especial de la difusión de calor en un medio isotrópico y homogéneo en un espacio
de tres dimensiones, la ecuación es:
∂u ∂ 2u ∂ 2u ∂ 2u
= α( 2 + 2 + 2 )
∂t ∂x ∂y ∂z
∂u
= α(uxx + uyy + uzz )
∂t
En donde:
u = u (x, y, z, t), es la temperatura como una función del espacio y del tiempo;
∂u
, es la razón de cambio de temperatura en un punto respecto del tiempo;
∂t
uxx , uyy , uzz , son derivadas segundas parciales (conducciones térmica) de la temperatura res-
pecto de x, y, z, respectivamente;
11
MATEMÁTICA SUPERIOR
k
α = , es la difusividad térmica, una cantidad especı́fica del material que depende de la
cp ρ
conductividad térmica k, la densidad de masa ρ, y la capacidad de calor especı́fica cp
La ecuación del calor es una consecuencia de la ley de Fourier de conducción (ver conduc-
ción del calor). Si el medio no es todo el espacio, entonces para resolver la ecuación de calor
se necesita solo que se especifiquen las condiciones de contorno para u. Para determinar la
unicidad de las soluciones en todo el espacio se necesita suponer una cota exponencial sobre
el crecimiento de las soluciones.3? Las soluciones de la ecuación del calor se caracterizan por
una suavidad gradual de la distribución de temperatura inicial por el flujo de calor desde las
áreas más cálidas hacia las más frı́as de un objeto. Generalmente, muchos estados diferentes
y condiciones iniciales tienden al mismo equilibrio estable. La ecuación del calor es un ejem-
plo prototipo de una ecuación diferencial en derivadas parciales del tipo parabólico Usando el
operador de Laplace, la ecuación se puede simplificar, y generalizar a ecuaciones similares en
el espacio de un número arbitrario de dimensiones:
2
ut = α∇ u = ∆u
donde el operador de Laplace puede ser ∆o∇2 , la divergencia del gradiente, se toma en las
variables espaciales. La ecuación del calor gobierna la difusión de calor, es decir, es otro proceso
de difusión, tal como la difusión de partı́culas o la propagación del potencial de acción en células
neuronales. Sin embargo no son de naturaleza difusiva, algunos problemas de mecánica cuántica
también están gobernadas por una ecuación análoga a la ecuación del calor (ver más adelante).
También puede utilizarse para modelar algunos problemas en finanzas, como por ejemplo en
los procesos de Black-Scholes o Ornstein-Uhlenbeck.
12
MATEMÁTICA SUPERIOR
8. Resuelva la ecuación
ut + 3ux = x2 t + 1
con la condición inicial u(x, 0) = x + 2
Solución. x−?3t + 2 + t + (x2 t2 )/2? − xt3 + 3/4t4
13
MATEMÁTICA SUPERIOR
10. Cuáles de las siguientes ecuaciones son lineales y cuáles no lineales. Determı́nelo encon-
trando el operador asociado y en caso que sea lineal demuéstrelo.
a)ux u2xy + 2xuuyy −?3xyuy −?u = 0
b)a(x, y, ux , uxy )uxyy + b(x, y, uyy )uyyy − f (x, y) = 0
c)cos(x + y)ux + cosuy = 0
Solucion:2a, 2b, 2c, no lineales
14