Capı́tulo 4

Ecuaciones en derivadas
parciales de primer orden

Objetivos
Resolver problemas de valores iniciales para ecuaciones diferenciales en
derivadas parciales de primer orden cuasilineales.

4.1. Introducción
Hasta el momento nos hemos ocupado de las ecuaciones diferenciales ordi-
narias, que son aquellas en las que las magnitudes que se pretende modelizar
por medio de ecuaciones diferenciales dependen tan sólo de una variable inde-
pendiente, normalmente un tiempo o en algunos casos una coordenada espacial.
Sobre todo en el último caso, planteando problemas en el plano o en el espacio,
es fácil imaginar problemas fı́sicos o ingenieriles en los que intervenga más de
una variable independiente.
Al igual que para las ecuaciones diferenciales ordinarias, el orden de una
ecuación en derivadas parciales es el orden más alto de las derivadas que apare-
cen en la ecuación.
Por ejemplo, la forma más general de una ecuación de primer orden para
una función u(x1 , . . . , xm ) de variables independientes x1 ,. . . , xm es

F (x1 , . . . , xm , u, ux1 , . . . , uxm ) = 0,

donde F es una función con ciertas condiciones de diferenciabilidad.

Notación: para aligerar la notación, muchas veces nos referiremos a las deriva-
das parciales por medio de subı́ndices,
∂u ∂2u ∂2u
ux := , uxx := , uxy := ,
∂x ∂x2 ∂x∂y
siempre que no cause confusión.
También, por aligerar la notación, como hemos hecho en el párrafo prece-
dente, eliminaremos las dependencias de las funciones. Por ejemplo, escribiendo
uxy en lugar de uxy (x, y), siempre que sean conocidas.

89
90 CAPÍTULO 4. EDP DE PRIMER ORDEN

En el caso de dos variables independientes, x, y, la expresión general se

reduce a F (x, y, u, ux , uy ) = 0. Las variables independientes pueden ser dos
coordenadas espaciales, x, y, o una coordenada temporal y una espacial, x, t,
dependiendo del problema.
Un ejemplo de ecuación de primer orden es la ecuación de Burgers,
∂u(t, x) ∂u(t, x)
+ u(t, x) = 0, (ut + uux = 0) ,
∂t ∂x
que se emplea en el estudio de ondas de choque en gases.
Esta ecuación es, como veı́amos en la parte de ecuaciones ordinarias, cua-
silineal, ya que es lineal en las derivadas de orden más alto, en este caso las
derivadas primeras respecto a t y respecto a x.
Una ecuación es lineal cuando es lineal en la variable dependiente y en todas
sus derivadas.
Por ejemplo, la ecuación de Burgers no es lineal, por culpa del producto del
segundo término.
En cambio, la ecuación del transporte,
∂u(t, x) ∂u(t, x)
+ = 0, (ut + ux = 0) ,
∂t ∂x
sı́ que es lineal.
Un ejemplo sencillo de ecuación de orden superior es la ecuación de Poisson,
∂ 2 V (x, y, z) ∂ 2 V (x, y, z) ∂ 2 V (x, y, z)
∆V (x, y, z) = + + = κρ(x, y, z),
∂x2 ∂y 2 ∂z 2
que rige el comportamiento de un potencial escalar V (x, y, z) gravitatorio (o
electrostático) en presencia de densidad de masa (o de carga) ρ(x, y, z). La
constante es κ = 4πG, donde G es la constante de gravitación universal, en el
caso gravitatorio. Y en el caso electrostático κ = −1/ε, siendo ε es la constante
dieléctrica del medio.
Este es un ejemplo de ecuación lineal en derivadas parciales de segundo orden
con tres variables independientes, las coordenadas x, y, z. La única variable
dependiente es el potencial V (x, y, z). En notación simplificada la ecuación de
Poisson se escribirı́a
Vxx + Vyy + Vzz = κρ.
La ecuación de la cuerda vibrante,
∂ 2 u(t, x) 2
2 ∂ u(t, x)


− c = F (t, x), utt − c2 uxx = F (t, x) ,
∂t2 ∂x2
describe la evolución de la separación u de una cuerda de su posición de equilibrio
en presencia de una fuerza f , siendo c la velocidad de propagación de la onda a
lo largo de la cuerda. Obviamente, la separación de la cuerda depende del tiempo
t, pero también depende de la posición x sobre la cuerda, ya que en unos puntos
se separará del equilibrio más que en otros, como sucede en la cuerdas de un
violı́n o una guitarra.
Esta es una ecuación lineal en derivadas parciales de segundo orden con dos
variables independientes, el tiempo t y la coordenada x a lo largo de la cuerda.
Las matemáticas también proporcionan ecuaciones en derivadas parciales de
forma natural:
4.1. INTRODUCCIÓN 91

Pensemos en superficies de revolución en torno al eje Z. Sabemos que las

superficies se pueden parametrizar por medio de gráficas de funciones de dos
variables, (x, y, z(x, y)). En el caso de las superficies de revolución,
p la altura z
se puede expresar tan sólo en función del radio cilı́ndrico ρ = x2 + y 2 . Por
tanto, se pueden describir por una función z(x2 + y 2 ). Usando la regla de la
cadena, denotando u = ρ2 = x2 + y 2 ,
∂z dz ∂u dz ∂z dz ∂u dz
= = 2x , = = 2y ,
∂x du ∂x du ∂y du ∂y du
con lo cual, eliminando dz/du de ambas expresiones obtenemos la ecuación en
derivadas parciales,
∂z(x, y) ∂z(x, y)
y −x = 0, (yzx − xzy = 0) ,
∂x ∂x
que satisfacen las superficies de revolución parametrizadas por gráficas de fun-
ciones (x, y, z(x, y)).
También podemos pensar en sistemas de ecuaciones en derivadas par-
ciales, sistemas de ecuaciones en los que hay más de una variable dependiente.
Un ejemplo de la fı́sica son las ecuaciones para un potencial vectorial A(x, y, z)
de una fuerza F(x, y, z), que, como sabemos, es solución de la ecuación vectorial
F = rot A, que en coordenadas,

∂Az (x, y, z) ∂Ay (x, y, z) x
− = F (x, y, z)  

∂y ∂z 


x z
∂A (x, y, z) ∂A (x, y, z) y
− = F (x, y, z) ,
y
∂z x
∂x 

∂A (x, y, z) ∂A (x, y, z) 

− = F z (x, y, z) 

∂x ∂y
proporciona un sistema de tres ecuaciones diferenciales lineales en derivadas
parciales de primer orden para las tres coordenadas del potencial, Ax , Ay , Az .
Denominaremos solución general de una ecuación en derivadas parciales
de orden n a una familia de soluciones de la ecuación que dependa de n funciones
independientes de las variables independientes.
Ejemplo 4.1.1 Hallar la solución general de la ecuación uxy (x, y) = F (x, y).
Este caso es sencillo, ya que se reduce a integrar en ambas variables,
Z Z
∂ux (x, y)
ux (x, y) = dy + g(x) = F (x, y) dy + g(x),
∂y
donde g(x) es la constante de integración al integrar en la variable y. A su vez,
Z Z Z 
∂u(x, y)
u(x, y) = dx + H(y) = F (x, y) dy + g(x) dx + H(y),
∂x
con lo cual la solución general de la ecuación es
Z Z 
u(x, y) = F (x, y) dy dx + G(x) + H(y),

denotando por G(x) la primitiva de g(x). Como se puede observar, depende

de dos funciones arbitrarias, G(x), H(y), como corresponde a una ecuación de
segundo orden.
92 CAPÍTULO 4. EDP DE PRIMER ORDEN

Ejemplo 4.1.2 Comprobar que la solución general de la ecuación homogénea

de la cuerda vibrante, utt − c2 uxx = 0 es u(x, t) = G(x − ct) + H(x + ct), siendo
F , G funciones arbitrarias.

En el tema siguiente se obtendrá dicha solución general, ası́ que aquı́ nos
limitaremos a comprobar su validez.
Denotemos v(x, t) = x − ct, w(x, t) = x + ct, con lo cual u(v, w) = G(v) +
H(w). Por la regla de la cadena, calculamos las derivadas primeras de la solución
general,

∂u ∂u ∂v ∂u ∂w dG dH
= + = −c +c ,
∂t ∂v ∂t ∂w ∂t dv dw
∂u ∂u ∂v ∂u ∂w dG dH
= + = + ,
∂x ∂v ∂x ∂w ∂x dv dw
que nos sirven a su vez para calcular las derivadas segundas,

∂ut ∂ut ∂v ∂ut ∂w d2 G d2 H

= + = c2 2 + c2 ,
∂t ∂v ∂t ∂w ∂t dv dw2
∂ux ∂ux ∂v ∂ux ∂w d2 G d2 H
= + = + ,
∂x ∂v ∂x ∂w ∂x dv 2 dw2
y, al sustituirlas en la ecuación de la cuerda,
 
2 d2 G d2 H d2 G d2 H
utt − c uxx = c2 2 + c2 −c 2
+ ≡ 0,
dv dw2 dv 2 dw2

comprobamos que efectivamente es una familia de soluciones de la ecuación.

Como, además, la familia depende de dos funciones arbitrarias, G y H y la
ecuación es de orden dos, efectivamente es una solución general de la ecuación.
La solución general de una ecuación en derivadas parciales no es sencilla
de conseguir en la mayorı́a de los casos, por lo que en general abordaremos
problemas de valores iniciales, problemas de contorno y problemas mixtos, tal
como hicimos con las ecuaciones ordinarias.

4.2. Ecuaciones cuasilineales de primer orden

Como ya hemos mencionado, se puede abordar el problema de Cauchy o
de valores iniciales para cualquier ecuación en derivadas parciales de primer
orden. En este curso, no obstante, nos limitaremos al caso de las ecuaciones
cuasilineales, comenzando con dos variables independientes, x, y, y una variable
dependiente u(x, y). La forma más general de ecuación cuasilineal para u es

a(x, y, u)ux + b(x, y, u)uy = c(x, y, u),

donde a, b, c son funciones de tres variables, con ciertas propiedades de diferen-

ciabilidad.
Podemos representar las funciones u(x, y) como superficies en R3 que sean
gráficas de funciones de dos variables z = u(x, y). La superficie estará formada
por puntos {(x, y, z) ∈ R3 : z = u(x, y)}. Es decir, a cada punto (x, y) del plano
le asignamos la altura dada por el valor u(x, y).
4.2. ECUACIONES CUASILINEALES DE PRIMER ORDEN 93

z=u(x,y)

z
u(x,y)

x y
(x,y)

Figura 4.1: Gráfica de la función u(x, y)

La superficie tiene por ecuación implı́cita F (x, y, z) = u(x, y) − z = 0 y, por

tanto, grad F (x, y, z) = (ux , uy , −1) proporciona un vector perpendicular a la
superficie en cada uno de sus puntos.

Consideremos curvas parametrizadas por δ(τ ) = x(τ ), y(τ ), z(τ ) en algún
intervalo y que estén contenidas en la superficie dada por la solución. Es decir,


F x(τ ), y(τ ), z(τ ) ≡ 0.

Derivando esta expresión con respecto al parámetro τ de la curva,


dF x(τ ), y(τ ), z(τ ) ∂F (x, y, z)  dx(τ )
0 = =
dτ ∂x 
(x(τ ),y(τ ),z(τ )) dτ
 
∂F (x, y, z)  dy(τ ) ∂F (x, y, z)  dz(τ )
+ +
∂y 
(x(τ ),y(τ ),z(τ )) dτ ∂z 
(x(τ ),y(τ ),z(τ )) dτ



= ẋ(τ )ux x(τ ), y(τ ), z(τ ) + ẏ(τ )uy x(τ ), y(τ ), z(τ ) − ż(τ ),

denotando con un punto la derivada respecto al parámetro τ .

Comparando esta expresión con la ecuación diferencial, aux + buy − c = 0,
observamos que las curvas que verifican el sistema autónomo

dx(τ )

= a x(τ ), y(τ ), z(τ ) ,
dτ
dy(τ )

= b x(τ ), y(τ ), z(τ ) ,
dτ
dz(τ )

= c x(τ ), y(τ ), z(τ ) , (4.1)
dτ
o en notación simplificada, donde el punto denota derivación respecto a τ ,

ẋ = a, ẏ = b, ż = c,

están contenidas en la superficie dada por alguna solución de la ecuación dife-

rencial.
El teorema de existencia y unicidad para sistemas de ecuaciones ordinarias
afirma que si las funciones a, b, c son de clase C 1 en un punto (x0 , y0 , z0 ), por
dicho punto pasará una sola curva solución del sistema.
94 CAPÍTULO 4. EDP DE PRIMER ORDEN

Obsérvese que un cambio de parámetro τ = τ (σ),


dδ τ (σ) dδ(τ )  dτ (σ)
= ,
dσ dτ τ (σ) dσ


simplemente conduce a multiplicar la velocidad δ̇(τ ) de la curva por la función

dτ (σ)/dσ y, por tanto, a multiplicar por esta misma función el lado derecho del
sistema 4.1.
Por ello, como el parámetro empleado es irrelevante, para no comprome-
ternos de antemano con una parametrización de las curvas, muchas veces se
reescribe el sistema como
dx dy dz
= = ,
a b c
de modo que desaparezca la referencia explı́cita al parámetro.
Las curvas que verifican el sistema 4.1 se denominan curvas caracterı́sticas
de la ecuación diferencial aux + buy − c = 0.
En el caso de que las funciones a, b no dependan de u,
a(x, y)ux + b(x, y)uy = c(x, y, u),
como sucede, por ejemplo, en el caso de las ecuaciones lineales, las dos primeras
ecuaciones del sistema caracterı́stico se desacoplan de la tercera,
dx(τ )

= a x(τ ), y(τ ) ,
dτ
dy(τ )

= b x(τ ), y(τ ) , (4.2)
dτ
y permiten definir una familia de curvas en el plano XY llamadas proyecciones
caracterı́sticas, o más confusamente caracterı́sticas, que, al no depender de
u, son comunes para todas las soluciones y son, por tanto, una propiedad de la
ecuación.
También podemos hablar de proyecciones caracterı́sticas en el caso cuasili-
neal como la proyección de las curvas caracterı́sticas sobre el plano XY , aunque
en este caso, en general, dependerán de la solución u considerada.
Ejemplo 4.2.1 Obtener las curvas caracterı́sticas de la ecuación ut + kux = 0.
En este ejemplo, a = k, b = 1, c = 0, con lo cual el sistema es
ẋ = k, ṫ = 1, ż = 0,
o equivalentemente, dividiendo todas las ecuaciones por la segunda,
dx(t) dz(t)
= k, = 0,
dt dt
y su solución general es
x(t) = kt + x0 , z(t) = z0 ,
donde x0 , z0 son constantes.
Las curvas caracterı́sticas en este ejemplo son rectas horizontales parame-
trizadas por δ(t) = (x0 + kt, t, z0 ), todas ellas paralelas, ya que su velocidad es
δ̇(t) = (k, 1, 0).
4.3. PROBLEMA DE VALORES INICIALES 95

Ejemplo 4.2.2 Obtener las curvas caracterı́sticas de la ecuación ut + uux = 0.

Este ejemplo no es lineal, a = z, b = 1, c = 0, y su sistema caracterı́stico es

ẋ = z, ṫ = 1, ż = 0,

o dividiendo todas las ecuaciones por la segunda,

dx(t) dz(t)
= z, = 0,
dt dt
y su solución general es

x(t) = x0 + z0 t, z(t) = z0 ,

donde x0 , z0 son constantes.

Las curvas caracterı́sticas en este ejemplo vuelven a ser rectas horizontales
parametrizadas por δ(t) = (x0 + z0 t, t, z0 ), pero no son paralelas, ya que su
velocidad es δ̇(t) = (z0 , 1, 0).

4.3. Problema de valores iniciales

Las curvas caracterı́sticas de una ecuación nos permiten construir la solución
de los problemas de valores iniciales. Supongamos que conocemos los valores
(dato inicial) u|Γ de la variable dependiente u a lo largo de una curva Γ (ver
Figura 4.2). Queremos conocer la solución u de la ecuación

a(x, y, u)ux + b(x, y, u)uy = c(x, y, u),

que toma los valores de u|Γ sobre la curva Γ. Es decir

Trazando las caracterı́sticas que pasan por los puntos de Γ, podemos cons-
truir la superficie de la solución de la ecuación que tiene por dato inicial u|Γ ,
siempre que Γ corte transversalmente a las caracterı́sticas.

z=u(x,y)
Γ

Figura 4.2: Gráfica de la solución u(x, y) como unión de caracterı́sticas trans-

versales a una curva Γ

Por contra, si la curva Γ fuese una caracterı́stica, la construcción anterior

serı́a imposible, ya que no podrı́amos continuar el dato inicial más allá de Γ. No
podemos en general dar el dato inicial sobre una curva caracterı́stica.
96 CAPÍTULO 4. EDP DE PRIMER ORDEN

De hecho, otra manera de construir las curvas caracterı́sticas es plantear el

problema de valores iniciales para la ecuación aux + buy = c en el que como
dato se dan los valores u|Γ de u a lo largo
de una curva Γ. Es decir,

si Γ está
parametrizada por γ(s) = f (s), g(s), h(s) , conocemos u f (s), g(s) = h(s).
Haciendo uso de la regla de la cadena, podemos conocer también los valores
de ux , uy a lo largo de Γ,

 
dh(s) du f (s), g(s) ∂u(x, y)  df (s) ∂u(x, y)  dg(s)
= = + ,
ds ds ∂x  f (s),g(s)
ds ∂y  f (s),g(s)
ds

con lo cual, podemos conocer ux |Γ , uy |Γ , resolviendo el sistema lineal,

h′ = u x f ′ + u y g ′ , aux + buy = c,

en el que las funciones a, b, c, ux , uy están restringidas a tomar valores sobre la

curva Γ. Con este procedimiento podemos obtener, no sólo las derivadas prime-
ras, sino las derivadas de orden superior, iterando la aplicación de la regla de la
cadena. Y conocidas todas las derivadas de u, podrı́amos expresar formalmente
la solución como una serie de Taylor.
Este sistema lineal está determinado y tiene solución única si y sólo si


 a f (s), g(s), h(s) f ′ (s) 

 6= 0.
 b f (s), g(s), h(s) g ′ (s) 

Es decir si los vectores (a, b), (f ′ , g ′ ) son paralelos, lo que equivale a que
los vectores (a, b, c), γ ′ = (f ′ , g ′ , h′ ) sean paralelos (recordemos que ambos tie-
nen que ser perpendiculares a grad F (x, y, z) = (ux , uy , −1)), el sistema no está
determinado. Pero precisamente (a, b, c) nos define la velocidad de las curvas ca-
racterı́sticas, con lo cual lo que estamos enunciando es que, para que el problema
tenga solución, debemos dar el dato sobre una curva Γ que no sea tangente (es
decir, que sea transversal) en ningún punto a ninguna caracterı́stica.

Teorema 4.3.1 Sea la ecuación a(x, y, u)ux + b(x, y, u)uy = c(x, y, u). El pro-
blema de valores iniciales a lo largo de una curva Γ tiene solución unica en un
entorno de Γ si a, b, c son funciones de clase C 1 y


 a f (s), g(s), h(s) f ′ (s) 

 6= 0, (4.3)
 b f (s), g(s), h(s) g ′ (s) 


donde γ(s) = f (s), g(s), h(s) es una parametrización de la curva Γ en un
intervalo.

Otra manera de interpretar la condición de transversalidad 4.3 es pensar

en un cambio de variables x = x(τ, s), y = y(τ, s) para resolver el problema
de valores iniciales. La condición para que el cambio de variables sea invertible
(difeomorfismo local) en un entorno de la curva Γ, que situamos en τ = 0, es
precisamente que el jacobiano del cambio no se anule,
 
 ∂x ∂x 
  

∂(x, y)  ∂τ ∂s   a f (s), g(s), h(s) f ′ (s) 
 
0 6= (0, s) =  = 
.
g ′ (s) 

∂(τ, s)  ∂y ∂y  b f (s), g(s), h(s)
 
 
∂τ ∂s (0,s)
4.3. PROBLEMA DE VALORES INICIALES 97

A efectos prácticos, si el
dato inicial está dado por una curva Γ parametrizada
por γ(s) = f (s), g(s), h(s) , tendremos que resolver el sistema caracterı́stico,

ẋ = a, ẏ = b, ż = c.

Las condiciones iniciales en τ = 0 son
x(s, 0) = f (s), y(s, 0) = g(s), z(s, 0) =

h(s). O, lo que es lo mismo, u f (s), g(s) = h(s).
Resolver el problema de valores iniciales para funciones arbitrarias f, g, h
equivale a obtener la solución general de la ecuación, en forma paramétrica,


u x(s, τ ), y(s, τ ) = z(s, τ ),

ya que la solución dependerá de los parámetros s, τ , en lugar de las coordenadas

x, y.

Ejemplo 4.3.1 Resolver el problema de valores iniciales para la ecuación lineal

ut + kux = 0, donde k es una constante, para u(x, 0) = h(x).

Nos están dando el dato inicial a lo largo del eje X, que será la curva Γ, por
lo que parece razonable tomar como parámetro s la propia coordenada x,

x(s, 0) = f (s) = s, t(s, 0) = g(s) = 0, z(s, 0) = h(s).

Los coeficientes de la ecuación son a = k, b = 1, c = 0.

Comprobamos la condición de transversalidad,
   
 a f ′ (s)   k 1 
 
 b g ′ (s)  =  1 0  = −1 6= 0.
 
(f (s),g(s),h(s))

Planteamos el sistema caracterı́stico,

ẋ = k, ṫ = 1, ż = 0,

cuya solución es

x(s, τ ) = kτ + α(s), t(s, τ ) = τ + β(s), z(s, τ ) = γ(s),

a la que imponemos la condición inicial,

s = x(s, 0) = α(s), 0 = t(s, 0) = β(s), h(s) = z(s, 0) = γ(s),

y concluimos que la solución del problema de valores iniciales, en forma pa-

ramétrica es

x(s, τ ) = kτ + s, t(s, τ ) = τ, z(s, τ ) = h(s).

No siempre podremos, pero en este caso es factible eliminar los dos paráme-
tros s, τ ,
τ = t, s = x − kτ = x − kt,
y obtener la solución del problema en función exclusivamente de las coordenadas
x, t, u(x, y) = z,
u(x, t) = h(x − kt),
98 CAPÍTULO 4. EDP DE PRIMER ORDEN

que corresponde a una onda viajera de velocidad k que se desplaza hacia la

derecha a lo largo del eje X, considerando que t es un tiempo.
Al dejar libre el dato inicial, ya que no hemos facilitado la expresión de h(s),
esta es la solución general de la ecuación, que depende de una función arbitraria,
la propia h.
Con esta expresión podemos resolver cualquier problema de valores iniciales.
Por ejemplo, si el dato inicial fuera u(x, 0) = sin x, es decir, h(s) = sin s, la
solución del problema serı́a simplemente u(x, t) = sin(x − kt).

Ejemplo 4.3.2 Resolver el problema de valores iniciales para la ecuación cua-

silineal ut + uux = 0, para u(x, 0) = h(x).

Esta es la ecuación de Burgers, en la que u representa la velocidad de un

fluido que recorre el eje X, con velocidad dependiente tanto del tiempo, como
de la posición.
La derivada ut representa, pues, la aceleración en un punto x en un instante
t. Si queremos conocer la aceleración de un partı́cula de fluido, tendremos que
realizar la derivada total,
 
du(x(t), t) ∂u(x, t)  ∂u(x, t)  dx(t)
= + = ut + uux |(x(t),t) .
dt ∂t 
(x(t),t) ∂x 
(x(t),t) dt

Por tanto, la ecuación de Burgers expresa la condición de que la aceleración

total de las partı́culas individuales es nula.
Nos están dando el dato inicial a lo largo del eje X, por lo que tomamos
como parámetro s la propia coordenada x,

x(s, τ ) = f (s) = s, t(s, τ ) = g(s) = 0, z(s, τ ) = h(s).

Los coeficientes de la ecuación son a = z, b = 1, c = 0.

Comprobamos la condición de transversalidad,
   
 a f ′ (s)   h(s) 1 
 
 b g ′ (s)  = 
 1
 = −1 6= 0.
(f (s),g(s),h(s))
0 

Planteamos el sistema caracterı́stico,

ẋ = z, ṫ = 1, ż = 0,

del cual podemos resolver inmediatamente las dos últimas ecuaciones,

t(s, τ ) = τ + β(s), z(s, τ ) = γ(s),

e imponerles las condiciones iniciales,

0 = t(s, 0) = β(s), h(s) = z(s, 0) = γ(s) ⇒ t(s, τ ) = τ, z(s, τ ) = h(s),

y sustituyendo en la primera obtenemos

ẋ = h(s) ⇒ x(s, τ ) = h(s)τ + α(s),

y la condición inicial nos permite despejar α(s),

s = x(s, 0) = α(s) ⇒ x(s, τ ) = h(s)τ + s,

4.3. PROBLEMA DE VALORES INICIALES 99

con lo cual ya tenemos la solución de todos los problemas de valores iniciales y,

por tanto, la solución general de la ecuación, en forma paramétrica,

x(s, τ ) = h(s)τ + s, t(s, τ ) = τ, z(s, τ ) = h(s).

Podemos eliminar los parámetros s, τ ,

τ = t, s = x − zt,

y obtener la solución en función de las coordenadas, x, t, u(x, y) = z, exclusiva-

mente,
u = h(x − ut),
que en este caso queda en forma implı́cita, al no poder despejarse, en general,
u en función de x, t.
Se trata de una onda viajera, pero de velocidad variable u.
Por ejemplo, en el caso en el que el dato inicial sea u(x, 0) = x, es decir,
h(s) = s, la solución del problema de valores iniciales es
x
u = x − ut ⇒ u(x, t) = .
1+t
La ecuación de Burgers tiene una peculiaridad: dado que sus proyecciones
caracterı́sticas son rectas de la forma x(t) = x0 + z0 t, cabe la posibilidad de que
estas rectas se corten, lo que no sucede con la ecuación del transporte, para la
cual las caracterı́sticas son paralelas.

(X,T)

~
x0 x0 x

Figura 4.3: Proyecciones caracterı́sticas de la ecuación de Burgers

Que haya dos proyecciones caracterı́sticas, x(t) = x0 + z0 t, x(t) = x̃0 + z̃0 t,

que se corten en un punto (X, T ) es grave, ya que la primera asigna a la solución
el valor u(X, T ) = z0 en dicho punto, mientras que la segunda le asigna el valor
u(X, t) = z̃0 , lo que es incompatible si z0 6= z̃0 .
Resolviendo la ecuación,
x̃0 − x0 z0 x̃0 − z̃0 x0
x0 + z0 T = X = x̃0 + z̃0 T ⇒ T = − , X= ,
z̃0 − z0 z0 − z̃0
observamos que esto sucede en un instante T > 0 en el futuro si el signo de
x̃0 − x0 es opuesto al signo de z̃0 − z0 . Teniendo en cuenta que u(x, t) es la
velocidad de las partı́culas del fluido, esto se puede evitar si la velocidad en
t = 0 es una función creciente de x.
100 CAPÍTULO 4. EDP DE PRIMER ORDEN

Esto es lógico, ya que, si las partı́culas más adelantadas (con mayor x inicial)
tienen una velocidad mayor, nunca chocarán con las partı́culas más rezagadas
(con menor x inicial).
Por ello la ecuación de Burgers se utiliza para modelizar ondas de choque,
que aparecen cuando unas partı́culas de fluido alcanzan a las que están situadas
por delante.

4.4. Solución general de la ecuación cuasilineal

Si queremos obtener directamente la solución general de una ecuación cua-
silineal sin pasar por la solución en forma paramétrica, podemos reescribirlo
formalmente como
dx dy du
= = ,
a b c
y reducirlo a un sistema de dos ecuaciones de primer orden, por ejemplo, to-
mando x como variable independiente,
dy b du c
= , = ,
dx a dx a
o cualquier otra.
La ventaja fundamental es que no aparece ningún parámetro. Como hemos
visto, esto es equivalente a un cambio de variable independiente de τ a x.
Como son dos ecuaciones de orden uno, su solución general podrá expresarse
en forma implı́cita por medio de dos constantes,

F (x, y, u) = C1 , G(x, y, u) = C2 ,

que se fijarán con el dato inicial para la curva caracterı́stica.

Si en vez de una sola curva caracterı́stica tuviéramos un problema de valores
iniciales con dato Γ, parametrizado por γ(s), las constantes dependerı́an del
parámetro s, C1 (s), C2 (s) y podrı́amos eliminarlo despejando, a través de o de
manera implı́cita,

C1 = f (C2 ), C2 = g(C1 ), h(C1 , C2 ) = 0,

a través de una función que no conocemos.

Cualquiera de estas tres expresiones nos da una forma de la solución general,

F (x, y, u) = f (G(x, y, u)) , G(x, y, u) = g (F (x, y, u)) ,

h (F (x, y, u), G(x, y, u)) = 0,

y usaremos la que más nos convenga en cada caso.
Lo vemos con un ejemplo:

Ejemplo 4.4.1 Hallar la solución general de la ecuación cuasilineal ut +uux =

0.

El sistema caracterı́stico en este caso quedarı́a como

dx
dt = , du = 0,
u
4.5. SOLUCIÓN GENERAL DE LA ECUACIÓN LINEAL 101

que resolvemos fácilmente, ya que la segunda ecuación se integra inmediatamen-

te, u = C1 , lo que nos permite integrar la primera,

C2 + ut = x ⇒ C2 = x − ut,

con lo cual podemos expresar la solución general de esta ecuación como C1 =

f (C2 ),
u = f (x − ut),
en forma implı́cita, como hemos visto. ✷
Obviamente, también se podrı́an haber empleado las otras expresiones equi-
valentes,
x − ut = g(u), h(u, x − ut) = 0,
pero resultan más incómodas.

4.5. Solución general de la ecuación lineal

Las ecuaciones lineales de primer orden y, en general, las ecuaciones de la
forma
a(x, y)ux + b(x, y)uy = c(x, y, u), (4.4)
se pueden reducir con un cambio de variables a ecuaciones ordinarias de primer
orden.
Para ello, hacemos uso de las proyecciones caracterı́sticas,

ẋ = a(x, y), ẏ = b(x, y),

que en estos casos son propiedades de la ecuación, ya que no dependen de una

solución concreta, puesto que la variable dependiente u no aparece en el sistema
de ecuaciones.
Tomando, por ejemplo, como parámetro τ la propia coordenada x (se podrı́a
tomar la coordenada y indistintamente), las proyecciones caracterı́sticas son
soluciones de la ecuación ordinaria de primer orden,
dy b(x, y)
= ,
dx a(x, y)

que tendrá una solución general de la forma K = F (x, y), en forma implı́cita,
siendo K una constante arbitraria.
Realizamos el cambio de variables U (x, y) = x, V (x, y) = F (x, y), suponien-
do que el cambio tiene jacobiano distinto de cero. En caso contrario, se podrı́a
tomar U (x, y) = y.
Por la regla de la cadena,
∂u ∂u ∂U ∂u ∂V ∂u ∂u ∂F
= + = + ,
∂x ∂U ∂x ∂V ∂x ∂U ∂V ∂x
∂u ∂u ∂U ∂u ∂V ∂u ∂F
= + = ,
∂y ∂U ∂y ∂V ∂y ∂V ∂y
podemos expresar la ecuación en las nuevas variables independientes U, V ,

(uU + uV Fx )a + uV Fy b = c,
102 CAPÍTULO 4. EDP DE PRIMER ORDEN

pero como a = ẋ, b = ẏ, el término

dF
uV Fx a + uV Fy b = uV Fx ẋ + uV Fy ẏ = uV = 0,
dτ
se anula, ya que F (x, y) = K es constante a lo largo de las caracterı́sticas.
Por tanto, la ecuación 4.4 se reduce a una ecuación ordinaria,

auU = c,

ya que no aparece V en las nuevas variables. La solución general de esta ecua-

ción contendrá una constante, que dependerá de V , la variable que etiqueta las
caracterı́sticas de la ecuación. ✷
Si hubiéramos tomado U (x, y) = y, la ecuación ordinaria resultante habrı́a
sido, como es fácil comprobar, buU = c.

Ejemplo 4.5.1 Hallar la solución general de la ecuación del transporte, ut +

kux = 0.

Esta ecuación es lineal y sabemos que sus proyecciones caracterı́sticas veri-

fican
dx
ẋ = k, ṫ = 1 ⇒ = k ⇒ x = kt + K.
dt
Tomamos como variables nuevas U (x, t) = x, V (x, t) = x − kt y la ecuación
se reduce a uU = 0, que tiene por solución u = h(V ). Es decir,

u(x, t) = h(x − kt),

como ya habı́amos obtenido previamente.