CAPITULO 1:

ECUACIONES DIFERENCIALES PARCIALES DE PRIMER ORDEN

Usando el mtodo de las caractersticas, resolveremos problemas de valor inicial
(PVI) para EDP de primer orden lineales, cuasilineales y no lineales en general,
discutiendo el comportamiento de la solucin en algunos ejemplos concretos. Este
mtodo est basado en una interpretacin geomtrica de la EDP que conduce a un
sistema de ecuaciones diferenciales ordinarias (EDO) no lineal en general, y no siempre
fcil de resolver.

ECUACIONES LINEALES

Sean a, b, c, y d funciones suaves (o regulares) y consideremos la EDP lineal (EDPL)

general de primer orden
(1.1)

a( x, t )

u
u
+ b( x, t )
= c( x, t )u + d( x, t ), x IR, t > 0.
x
t

Idea del mtodo: La solucin u de (1.1) es una superficie en el espacio (x,t,u) que
podemos parametrizar utilizando dos parmetros, digamos s y . Para hallar la
expresin de u=u(x,t) suponemos dada una curva inicial C de parmetro , que ser el
dato inicial o condicin inicial (CI) que acompaa a (1.1).
Si por cada punto de C podemos trazar otra curva de parmetro s, entonces la
superficie as determinada ser la solucin buscada.
El lado izquierdo de (1.1) tiene la forma de una derivada direccional de u(x,t) en la
direccin del vector [a, b] en los puntos (x,t) donde [a, b] est definido y no sea nulo.
Luego las ecuaciones
(1.2)

dx
= a( x, t );
ds

dt
= b( x, t ) ,
ds

s = parmetro,

determinan una familia uniparamtrica de curvas en el plano (x,t):

(1.3)

x=x(s) ; t=t(s),

cuyos vectores tangentes [x '(s), t '(s)] tienen la misma direccin de [a, b] (ver Fig.1).
Por lo tanto la derivada total de u sobre las curvas (1.3) ser
(1.4)

du u dx u dt
=
+
ds x ds
t ds

u
u
a+
b
x
t
= cu + d
Luego tenemos el sistema de EDO
=

Fig. 1 Derivada Direccional.

dt
du
dx
=a;
=b;
= cu + d,
ds
ds
ds

(1.5)

llamado ecuaciones caractersticas de (1.1).

Obviamente las dos primeras ecuaciones se pueden resolver independientemente de la
tercera.
Las soluciones de (1.5) son familias de curvas: x=x(s), t=t(s), u=u(s), llamadas curvas
caractersticas (CC). Estas curvas en el plano (x,t) se llaman curvas caractersticas
bases (CCB). Bajo hiptesis de suavidad sobre los datos a, b, c, y d, la teora de EDO
garantiza una nica solucin (x(s), t(s), u(s)), que ser una CC del espacio (x,t,u) que
pasa por un punto dado (x0, t0, u0).
No estamos interesados en la solucin general de (1.1) sino que en hallar una solucin
particular u = u(x,t) que ser una nica superficie en el espacio y que contenga una
curva dada C. Esta curva ser la CI para (1.1). La EDP y la CI constituyen un PVI o,
ms especficamente, un Problema de Cauchy.
Mtodo. Sea C la curva inicial dada paramtricamente por
(1.6)

C:

x=x() ;

t=t() ;

u=u(), = parmetro.

Para cada valor de tenemos un punto de C por el cual pasa una CC, La familia de
curvas determinadas por los puntos de C est dada por:
(1.7)

x = x(s,) ,

t = t(s,)

, u = u(s,)
2

tal que cuando s = 0 tenemos C, es decir, x(0,) = x() , t(0,) = t() y u(0.) = u().
Las ecuaciones (1.7) son una representacin biparamtrica de una superficie en e
espacio (x,t,u) que contiene a C.
Resolviendo los PVI para las EDO (1.5):
dx
= a( x, t )
ds
x s=0 = x( )

dt
= b( x, t )
ds
t s=0 = t( )

du
= c( x, t )u + d( x, t )
ds
u s=0 = u( )

obtenemos: x=x(s), t=t(s,) y u=u(s,).

Si de las dos primeras ecuaciones podemos despejar s y en funcin de x y t, y
sustituyendo estos valores en u=u(s,), entonces obtenemos u(x,t). Esta funcin u(x, t)
as conseguida satisface (1.1) en al menos una vecindad de C, y la CI (1.6). Adems es
la nica solucin del PVI.

Fig. 2 Curvas caractersticas.

Si este mtodo no conduce a una solucin entonces no existe o existen infinitas
soluciones. Este ltimo caso resulta si, por ejemplo, C coincide con una curva
caracterstica; este tipo de problema se conoce como "Problema de Valor Inicial
Caracterstico".

Ejemplo 1. Movimiento de ondas unidireccional. La EDP de primer orden

u
u
(1.8)
+c
= 0 , c = cte., x IR, t > 0
t
x
representa un movimiento de ondas unidireccional en IR (conocida tambin como
ecuacin de conveccin). Resolveremos un PVI para (1.8) con la condicin inicial
(1.9)

u(x,0) = F(x), x IR

siendo F una funcin dada.

Sol:
(1.10)

Observando que la CI est dada a t=0, parametrizamos C como sigue:

C: x= , t=0 , u=F().

Las ecuaciones caractersticas son:

(1.11)

dx
=c,
ds

dt
du
= 1,
=0
ds
ds

que, junto a las CI conforman los tres PVI para estas EDO, cuyas soluciones son
(1.12)

x(s,) = cs + , t(s,) = s ,

u(s,) = F().

Despejando y s en funcin de x y t en las dos primeras expresiones

(1.13)

s = t , = x - ct

y sustituyendo en u(s, ), obtenemos

(1.14)

u(x, t) = F(x - ct).

Si F C1(IR) entonces (1.14) es la solucin del PVI (1.8), (1.9). Notar que la solucin
est determinada slo en trminos del dato inicial.
Interpretacin fsica de la solucin.
Si F representa la forma inicial de una onda, entonces u(x, t) = F(x - ct) muestra que la
onda no cambia de forma cuando t crece (e.d. a medida que el tiempo transcurre).
Adems, si c >0, entonces la onda "viaja" a la derecha con velocidad dx/dt = c; y si c<0
la onda "viaja" a la izquierda. Las CC son x-ct = . De este modo, para cada valor de ,
u es constante sobre las curvas caractersticas. Esto significa que el dato inicial se
transmite a lo largo de las CCB, en el sentido que cualquiera sea el valor de u a t=0 en
algn x, ese valor es el mismo en todos los puntos x pertenecientes a las CCB.

Fig. 3 Propagacin unidireccional de ondas.

Idea grfica.
Si el dato F(x) es la delta de Dirac (x), entonces la solucin "formal" ser u(x,t)=(x-ct),
es decir, la solucin es idnticamente nula salvo sobre las caractersticas sobre las
cuales se propaga el dato inicial.
Con todo lo anterior, ya podemos afirmar que las perturbaciones se propagan a lo
largo de las curvas caractersticas bases. (Fig. 4)

Fig. 4 Propagacin de una perturbacin (la funcin ).

5

Si el dato F(x) es una funcin discontinua, por ejemplo u(x,0) = H(x+1) - H(x-1) siendo
H(x) la funcin de Heaviside:

(1.15)

1, x > 0
H( x ) =
,
0, x 0

entonces la solucin de (1.8) con c=1 tiene la forma u(x, t) = H(x+1-t) - H(x-1-t).

Fig. 5 Propagacin de discontinuidades.

Del grfico de esta ltima funcin, deducimos que no slo las perturbaciones se
propagan sobre las CCB, sino que tambin las discontinuidades (ver Fig. 5). Lo mismo
ocurre en los puntos donde F(x) no es derivable. En efecto, basta considerar por
ejemplo el dato inicial
(1.16)

0, x > 1
u(x,0)= R(x)=
1 x , x 1

entonces u(x, t) = R(x-ct), grficamente indica que la propagacin de los saltos de las
derivadas ocurre sobre las CCB (ver Fig.6).
De las ecuaciones caractersticas bases x-ct = = constante, deducimos que el
dominio de dependencia de u(x,t) sobre el dato inicial est representado por u punto
x= sobre el eje de las x. La influencia del dato inicial en un punto particular sobre la
solucin u(x,t), justamente son todos los puntos de la CC x-ct = .

Fig. 6 Propagacin de los saltos de discontinuidad.

NOTA: Este es un buen ejemplo para iniciar el estudio de la solucin numrica por el
mtodo de las diferencias finitas. (Ver F John, " Partial Differential Equations", Springer
Int. 1978).
Ejemplo 2. Solucin D'Alembert de la ecuacin de ondas.
El mtodo de las caractersticas permite resolver un Problema en EDP para una
ecuacin lineal de segundo orden que rige el movimiento de ondas sobre IR con
velocidad >0:
(1.17)

2
2u
2 u

= 0 , x IR , t > 0
t 2
x 2

Para tener un PVI agregamos las condiciones iniciales

(1.18)
(1.19)

u (x,0) = f(x)
ut(x,0) = g(x)

siendo f y g funciones dadas.

Para resolver (1.17) lo reescribimos como un sistema equivalente de EDP de primer
orden. Para ello introducimos los operadores diferenciales x, t definidos por
2
2u
2 u

= ( 2t 2 2x )u = ( t + x )( t x )u = 0
2
2
t
x
Sea (t - x)u = v, luego (t + x)v = 0. De este modo tenemos el sistema acoplado

(1.20)

u
u

= v( x, t )
t
x
(1.21)
v
v
+
=0
t
x
Debemos hallar la condicin inicial para v(x, t). De la primera ecuacin de (1.21)
tenemos, a t=0, v(x, 0)= g(x) -f ' (x). Luego el PVI para v(x, t) ser
v
v
+
=0
t
x
v(x ,0) = g(x) -f ' (x),

cuya solucin, de acuerdo al ejemplo 1, es v(x, t) = g(x -t) - f ' (x -t).

Sustituyendo en la primera ecuacin de (1.21), tenemos el PVI para u(x, t):
u
u

= g( x t ) f ' ( x t )
t
x
u( x, 0) = f ( x )

(1.22)

Ahora resolvemos (1.22) por el mtodo de las caractersticas:

dx
= ,
EDO: ds

C:
Luego,

x=

dt
=1 ,
ds

t=0 ,

x = -s+ , t = s ,

du
= g f '
ds

u = f().
u=

(g f ' )d + f ()
0

Despejando y s en funcin de x y t de las dos primeras ecuaciones y considerando

que f y g son funciones de x-t , resulta : x+s = x-s = -2s,
s

u(s , ) = [g( 2) f ' ( 2)]d + f ( ) .

0

Sea = - 2, d = -2d. Luego,

u(s, ) =

1
2

2 s +

g( )d +

1 2 s +
f '( )d + f ( )
2

1 2 s+
1
g( )d + [ f ( 2s + ) f ( )] + f ( )

2
2

1 2 s+
1
g( )d + f ( 2s + ) + f ( )

2
2

es decir,
8

x t
1
[f ( x t) + f ( x + t)] 1 x+ t g()d
2
2
y de aqu resulta la solucin D'Alembert de la ecuacin de ondas

u(s, t) =

1
[f ( x t) + f ( x + t)] + 1
2
2
Interpretacin fsica de la solucin.
(1.23)

Sea G(z) =

u(x, t) =

x + t

x t

g( )d .

g( )d la antiderivada de g, luego (1.23) puede escribirse como

1
1

1
1
G( x + t ) + f ( x t ) G( x t ) .
u(x , t) = f ( x + t ) +
2
2
2
2

El primer sumando de (1.24) representa ondas que viajan a la izquierda con velocidad ,
y el segundo sumando, ondas que viajan a la derecha con la misma velocidad.
La expresin (1.23) (o bien (1.24)) nos dice que u(x, t) depende de los extremos del
intervalo x t a travs del dato f(x), y de todos los puntos de ese intervalo a travs del
dato g(x). Luego, si por ejemplo g 0, la onda que se propaga slo depende de los
puntos x t. Esta propiedad se expresa diciendo que no hay difusin de la onda, o
que la onda se propaga sin difusin obedeciendo el Principio de Huygens.
Interpretamos grficamente esta doble propagacin con los datos: g 0, =100, y f(x)
es una onda suave con soporte en [-2, 2], tal como muestra la Fig. 7.

(1.24)

Fig. 7 Propagacin bidireccional de ondas.

Observamos que la onda no cambia de forma, tal como en el ejemplo 1. Sin embargo,
hay una superposicin al comienzo. Por cada dcima de segundo la onda se traslada a
9

derecha e izquierda en una unidad de longitud y se suman las ondas obtenidas.

Volveremos sobre esta ecuacin ms adelante.
Ejemplo 3. Una ecuacin con singularidades.
Consideremos el PVI
u
u
x
+t
= cu , c = cte. , x IR, t > 0
x
t
u( x, 1) = f ( x )
Observamos que el vector [x, t] es nulo en (0,0); es natural preguntarse entonces cul
ser el efecto de esta singularidad en la solucin ? (notar que la condicin inicial no est
dada a t=0).
(1.25)

Sol. Tenemos los problemas en EDO:

dx
=x
,
ds
x( 0 ) =
cuyas soluciones son

(1.26)

(1.27)

dt
du
=t
= cu
,
ds
ds
t(0) = 1 u(0) = f ( )

x(s, )= es , t(s, )=es , u(s, )=f()ecs

(notar que no puede darse la condicin inicial a t=0, pues en tal caso no podemos
resolver el PVI para dt/ds=t).
x
y s = ln t, luego
De aqu resulta =
t
x
u(x, t) = f t c
t
es la solucin de (1.25), pero no es vlida para todo x y todo t. En efecto, (1.28) tiene
una singularidad en t=0; basta considerar que si f es una constante k y c>0, entonces la
solucin u = ktc no ofrece problemas, pero si c<0 entonces la solucin explota a t=0.

(1.28)

Ejemplo 4. Un problema con infinitas soluciones.

Consideremos la EDP (1.25) con c=1,
(1.29)

u
u
+t
=u
x
t

dx
dt
dx x
=x y
= t , resulta
= , t 0 , de
ds
ds
dt t
donde x = t, =cte., son las curvas caractersticas bases (CCB), esto es, un haz de
rectas que pasan por el origen (excepto el eje x). La ecuacin caracterstica para u ser

Combinado las ecuaciones caractersticas

10

du
ds = du = u
dt dt
t
ds
de donde u = t, =cte.
Debemos exigir u x =t = t con , fijados. Parametrizamos C como:

, u = .

De (1.26) y (1.27), para esta C, obtenemos las soluciones

C:

(1.30)

x = ,

t=

x(s, ) = es

t(s, ) = e s ,

u(s, ) = e s .

De las dos primeras ecuaciones de (1.30) no podemos despejar s y en funcin de x y

t. Sin embargo, u = (/)x (=u(0, )) y u = t (=u(t, t)) son soluciones de (1.25). Luego
no existe una nica solucin. Pero hay ms, como puede cambiar de una curva
caracterstica a otra, de u = t se sigue que
x
u(x, t) = g t
t
es una solucin general de la EDP xux + tut = u, cualquiera sea gC1(IR); y como
x
g( z) z= = concluimos que g t tambin es solucin del PVI (1.25), es decir,
t
tenemos infinitas soluciones.
Por otro lado, si las condiciones iniciales se dan sobre las CCB x = t tal u no puede
ser expresada como u =t, entonces no existe solucin.

(1.31)

NOTA: Tcnicamente este problema no es un PVI-caracterstico pues la curva C en el

espacio (x,t,u) no es una curva caracterstica.

ECUACIONES CUASILINEALES.

Ahora estudiaremos una forma especial de EDP no lineal cuya forma general es
(1.32)

a(x,t,u)

u
x

b(x,t,u)

u
= c(x,t,u),
t

donde a, b y c son funciones dadas con cierta regularidad, digamos de C1 en x , en t y

en u. En particular, dado que la no linealidad est dada slo en los trminos de mayor
orden, la ecuacin se llama cuasilineal.
Un ejemplo tpico de PVI asociado a una ecuacin cuasilineal es el que modela el
trfico en carreteras:

11

ut + g(u)ux = 0
u(x, 0) = f(x),
sobre el que volveremos ms adelante.
Tal como en el caso lineal, hallaremos la solucin interpretando geomtricamente la
ecuacin diferencial.
Supongamos que existe una solucin u=u(x,t) de (1.32), e.d., que existe una superficie
integral que expresamos implcitamente como
(1.33)

F(x,t,u) = u(x,t) - u = 0.

Tomando diferenciales obtenemos la ecuacin de la tangente a F en todo punto

F
F
F
dx +
dt +
du = 0 ,
x
t
u
e.d., uxdx+utdt-1du=0, de donde, tenemos la diferencial total du = uxdx+utdt. Luego, el
vector normal a la superficie F(x,t,u)=0 en todo punto ser:

F n = [u x , u t , 1]

Por lo tanto (1.32) puede escribirse como

(1.34)

u
u
u u

+b
c = [a, b, c ] , , 1 = 0,
x
t
x t

lo que implica que [a,b,c] pertenece al plano tangente a la superficie u=u(x,t) en todo
punto (x,t,u) donde F 0 .
De este modo, la ecuacin (1.32) puede considerarse como el requerimiento geomtrico
que toda superficie integral o solucin u(x,t) que pasa por el punto (x,t,u) debe ser
tangente el vector [a,b,c]. Si en algn punto de la superficie empezamos a movernos en
la direccin del vector tangente conocido [a,b,c], lo hacemos sobre una curva que est
completamente contenida en la superficie F(x,t,u)=0. Esta curva se llama caracterstica,
y hallando el conjunto de todas estas curvas caractersticas, hallamos la solucin.

Sea r = x i + t j + u k el vector posicin de un punto sobre una curva caracterstica.

Siguiendo la idea de 1.1, expresamos paramtricamente r en trminos de s (ver Fig. 8).
Por lo tanto, el vector tangente a la curva en un punto particular ser:
dr dx dt du
=
i+
j+
k.
ds ds
ds
ds
Como los vectores [xs, ts, us] y [a,b,c] son tangentes a la curva, sus componentes son
proporcionales, e.d.,

(1.35)

dx dt du
ds = ds = ds ,
a
b
c

o bien
12

Fig. 8 Parametrizacin de una curva caracterstica.

(1.36)

dx
= a( x, t, u) ;
ds

dt
= b( x, t, u) ;
ds

du
= c( x, t, u) .
ds

Las ecuaciones diferenciales (1.36) son la ecuaciones caractersticas de (1.32) y sus

soluciones son las curvas caractersticas de (1.32).
NOTA: De la teora de EDO sabemos que si a, b, c son suficientemente suaves y no se
anulan en algn punto, entonces existe una nica curva que pasa por cada punto (x0, t0,
u0).
La condicin inicial para (1.32) resulta de especificar el valor de u(x,t) sobre una curva
dada del plano (x,t), determinndose una curva C en el espacio (x,t,u), llamada curva
inicial.
Resolucin del PVI: Para hallar la solucin, pasamos una curva caracterstica por cada
punto de la curva inicial C. Si estas curvas generan una superficie, sta superficie es la
solucin del PVI.
Para reafirmar la existencia y unicidad de solucin, damos el siguiente resultado:
TEOREMA: " Sean a,b,c C1 en x,t,u. Sea C: x=x(), t=t(), u=u(), tal que C tenga
vector tangente continuo y que

(1.37)

dt

dx

( ) = d a[x( ), t(), u()] d b[ x(), t(), u( )] 0

sobre C. Entonces existe una nica solucin u=u(x,t) definida en alguna vecindad de C
satisfaciendo el PVI:

13

(1.38)

u
u

a( x, t, u)
+ b( x, t, u) = c( x, tu)
.

x
t

u[ x( ), t( )] = u( ) "

Idea de la demostracin: De la teora de EDO, sabemos que el sistema caracterstico

dx
= a( x, t, u)
ds
u s=0 = x( )

tiene una nica solucin

(1.39)
x = x(s,);

dt
= b( x, t, u)
ds
t s=0 = t( )

du
= c( x, t, u)
ds
u s=0 = u( )

t = t(s,);

u = u(s,)

con derivada continua en s y tal que

(1.40)

x(0,) = x()

t(0,) = t()

u(0,) = u().

x = x(s, )
en s=0 es
El jacobiano de la transformacin
t = t(s, )

x x
( x, t )
s = dt a dx b
(1.41)
s =0 =
s =0 = ( )
t t d
(s, )
d
s
que es distinto de cero por hiptesis, y por continuidad no se anula en una vecindad de
C. Luego, por el teorema de la funcin implcita, podemos despejar s y como
funciones de x y t cerca de C, e.d.,
(1.42)

u(x,t) = u[s(x,t), (x,t)]

es una solucin del PVI (1.38).

La unicidad resulta del hecho que dos superficies integrales cualesquiera, que
contienen la misma curva inicial, coinciden sobre todas las CC que pasan por la curva
inicial.
NOTA: La hiptesis (1.37) significa que la curva inicial C no es caracterstica ni es la
envolvente de CC.
Si C es una curva caracterstica, entonces existe solucin pero no es nica. Si C es una
envolvente entonces u puede o no ser diferenciable sobre C. Si C no es caracterstica
salvo en un conjunto discreto de puntos, entonces surgen problemas especiales en la
vecindad de esos puntos.
Ejemplo 5. Movimiento unidireccional no lineal de ondas.
Consideremos el siguiente PVI

14

(1.43)

u
u

EDP :
+u
=0
,

t
x
CI : u( x,0) = f ( x )

siendo f(x) una funcin dada.

SOL.: Parametrizamos C como sigue:
(1.44)
x=;
t=0 ;
u = f()

Aplicamos (1.37): () = -1.

Las ecuaciones caractersticas son:
dx
dt
du
= u;
= 1;
= 0;
ds
ds
ds
du
= 0 resulta que u=cte. sobre C, luego u(s,) = u(0,) = f().
Del hecho que
ds
Por lo tanto,
dx
dt

= u = f ( )
=1

x
(
s
,

)
=

+
sf
(

)
y
ds
ds

t(s, ) = s .

x s =0 =
t s =0 = 0

(1.45)

Si u verifica la hiptesis del teorema, el sistema anterior puede invertirse cerca de C

para dar: s = t ; = (x,t). Luego la solucin ser
(1.46)

u = f((x,t))

y como s = t y = x- sf()=x - tu, resulta

(1.47)

u = f(x - tu).

Discusin fsica de la solucin: Si u=u(x,t) representa una onda, entonces la forma de la

onda, en todo instante, est dada por la curva u = u(x,t) en el plano (x, u) con t=
parmetro y u = f(x) la forma inicial de la onda.
Claramente el valor numrico de u(x,t) da la altura de la onda en x y a tiempo t.
Cmo de mueve la onda?. Primero debemos hallar la velocidad dx/dt en cada punto.
Las dos primeras ecuaciones de (1.45) dan
(1.48)

dx
=u,
dt

e.d., la velocidad de propagacin no es constante. Por lo tanto, a mayor amplitud de la

onda, u(x,t), mayor velocidad (distinto al caso lineal !). Luego, si u>0 x va a la
derecha; si u=0 x permanece fijo y si u<0 x va a la izquierda. Por lo tanto si u cambia de
signo hay partes de la onda que van hacia la izquierda y otras a la derecha. En tal caso
el viaje no es unidireccional.

15

Supongamos que u>0; luego los puntos x tal que u(x,t) es grande y la onda el "alta", se
movern ms rpido que los puntos x donde u(x,t) tiene un valor menor y la onda es
ms "baja". Por lo tanto si a t=0 hay partes de altas de la onda a la izquierda de otros
puntos, stos pueden alcanzar los puntos bajos y pasarlos; ah se produce el quiebre
de la onda .

Fig. 9 Quiebre de ondas.

En el tiempo de quiebre (tq), u(x, t) es multvoca y deja de ser una solucin vlida para
el PVI (pensar que si u representa una densidad, la multivocidad no tiene sentido).
Por qu se produce tan desagradable suceso?. En general, porque en el modelo (1.43)
se han despreciado trminos con derivadas de mayor orden (existe un modelo mejor: la
+ 1
ecuacin de Burger: u t + 1 +
uu x u xx = 0 .
2
Observacin importante: El PVI
v
v D 2 v
= c
+
t
x 2 x 2
v(x,0)= (x),
donde c y D son constantes y es la "delta de Dirac" , tiene como solucin
( x ct )2

v(x,t) = kexp
2Dt

que es muy regular para todo t >0. Pero si D=0 (e.d., consideramos la EDP de primer
orden) la solucin es v(x,t) = (x-ct) que es singular para todo t>0. Es decir, al incluir
trminos con derivadas de mayor orden en la EDP, la ecuacin tiende a suavizar la
solucin.

16

NOTA: Existe otro mtodo, llamado "ondas de choque discontinuo", que extiende la
validez de la solucin ms all del tiempo de quiebre (la ecuacin de la onda de choque
depende fuertemente del fenmeno en estudio).
Volvamos al PVI (1.43).
Volvamos al PVI (1.43): la forma de la onda depende del dato inicial, en particular si la
onda se quiebra, tambin depende del dato inicial. Cabe preguntarse si es posible hallar
el tiempo de quiebre tq. La respuesta es afirmativa y veremos dos mtodos:
1er Mtodo. Determinemos la pendiente de la forma de la onda u= u=(x,t)=f(x-tu) a
tiempo t. Derivando implcitamente

f ' ( x tu)
u
=
.
x 1 + tf ' ( x tu)
De x(s,) = +sf() ; t(s,) = s resulta que las curvas caractersticas bases (CCB) para
(1.43) es la familia de rectas: x tf()= , siendo el parmetro.
u
y la onda
Cuando 1+tf '(x-tu)=1+tf '() se anula por primera vez, entonces
x
empieza a quebrarse; luego
(1.49)

tq = -

(1.50)

1
f ' ( )

es el tiempo de quiebre.
2do Mtodo. Recordemos que si F(x, t,)=0 es una familia uniparamtrica de curvas con
parmetro , la envolvente se encuentra eliminando del sistema:
F( x, t, ) = 0

F ( x, t, ) = 0

Considerando las curvas caractersticas x=+sf(); t=s. Resulta que u=f() sobre las
CCB
(1.51)
x tf()= , =parmetro.
Si para dos o ms valores de las rectas (1.51) se cortan., entonces u es multivaluada
en el punto de interseccin (pues u(x, t) debe ser igual a f() en ese punto).
Para hallar los puntos de interseccin debemos determinar la envolvente de (1.51).
F = x-tf() - = 0
Tenemos as:

F
= -tf ' () 1 =0.

17

x tf ( ) = 0
1
, e.d., las
x-tf() = , lo que da tf '(x-tf())=1, de donde t = tf ' ( ) 1 = 0
f ' ( )
1
.
veces que las CC tocan la envolvente est dado por t = f ' ( )
1
siempre es negativo y la onda no se quiebra.
De (1.50), si f ()>0 , entonces f ' ( )
Geomtricamente esto significa que u = f(x) es una funcin montona creciente en x
(ver Fig. 10).

Fig. 10 No existe quiebre de la onda.

Si f '()<0 , entonces u=f(x) es montona decreciente y los puntos ms altos de la

onda se mueven ms rpido que los puntos bajos y pueden sobrepasarlos,
producindose as el quiebre d la onda. (si f ' cambia de signo, generalmente hay
quiebre de onda).

Fig. 11 Existe quiebre de la onda.

Hay otra interpretacin de los resultados anteriores sobre el quiebre de la onda: De

18

1
y x - tf() = siendo =parmetro, resulta que el tiempo de quiebre est
f ' ( )
determinado por el punto sobre la curva

t=-

f ( )
1
;t=; u = f(),
f ' ( )
f ' ( )
Donde t alcanza el mnimo valor, es decir, cuando el hessiano

Cs=0 : x = -

x
s
t
s

x
t
t

= 0.
s=0

Pero de (1.37) esto implica que ()=0 en ese punto; luego si C se elige como curva
inicial C para la EDP ut+uux=0, el teorema no garantiza solucin nica. Adems,
f ' ( x tu) C f ' ( )
u
=
=
u no es diferenciable sobre C.
x 1 + tf ' ( x tu) 1 1
Luego, aun cuando el dato inicial garantice solucin nica, sta solucin se quiebra si la
condicin ()0 fuese violada. Esta curva quebrada no es una CC para la EDP pues u
no es constante sobre ella.
Finalmente, si la curva inicial es una CC entonces el PVI tiene infinitas soluciones. En
efecto, recordemos que las ecuaciones caractersticas para nuestro problema son:
dx
dt
du
=u ;
=1 ;
=0.
ds
ds
ds
Luego la recta x=ct es una caracterstica para la EDP pues con t=s, x=cs satisfacen las
ECP. Sea u=c=cte. Sobre x=cxt. Entonces u=f(x-tu) es una solucin del PVI, pero
tambin los es f(0)=c.

Tres ejemplos de f(x) para ut + uux =0.

Ejemplo 6. Consideremos el dato inicial

(1.52)

u(x,0) = f(x) = -x.

Como f '(x) = -1, la onda debe quebrarse, e.d. u ser multvoca en algn punto. En
efecto, como f() = -, tq=1 es el tiempo de quiebre.
x
Adems, u=f(x-tu) = -(x-tu) u = -x +tu, es decir, u =
es la solucin del PVI
t 1
u t + uu x = 0
, solucin que explota en t=1.

u( x,0) = x
En la Fig. 12 observamos el comportamiento de la solucin de t=0 a t=t1, concluyendo:

19

i)
existe una rotacin alrededor del origen en el plano (x,u)
ii)
u crece linealmente respecto de x
iii)
t=1 u = eje u (multivaluada all)
La multivocidad de u en t=1 tambin puede visualizarse considerando las CCB:x= (1-t).
iv)

x
.
t 1
Para t = 1, todas las rectas se cortan en (x,t)=(0,1). A cada una de estas rectas le
corresponde un valor de y como u(x,t) es constante sobre cada una de estas rectas, u
ser multvoca en t = 1 (ver Fig. 13)

Fig. 12 Movimiento de la onda u =

Fig. 13 Multivocidad de u en t=1

20

NOTA: Obviamente, si f(x)=x entonces u(x,t) =

x
est bien definida para todo t>0.
t +1

Ejemplo 7. Consideremos ahora la CI: u(x,0)= f(x) = 1-x2.

Usando la forma implcita u=f(x-tu) la solucin es u=1-(x-tu)2 = 1-x2 +2xtu u2t2, e.d.,
u2t2+u(1-2xt)+(x2-1)=0, de donde
x 1 + 1 + 4 t( t x )

t
2t 2
Luego la solucin tiene dos ramas.
x 1 1 + 4 t( t x )
Si calculamos lim
,debemos considerar el signo ms en (1.53) para
t 0 t
2t 2
que as u(x,t) verifique la CI u(x,0)=1-x2, luego

(1.53)

(1.54)
de donde

u(x,t) =

x 1 + 1 + 4 t( t x )

,
t
2t 2
u 1
1
= +
si 1+4t(t-x)0, e.d.,
t t t 1 + 4t( t x )
u=

la pendiente

cuando
(1.55) x = t +

1
.
4t

El radical en (1.54) es real si

x<t+

1
, luego u est bien definida ah. Ahora, t>0;
4t

1
x < 1, es decir, u est bien definida para x<1. En x=1, el radical se anula
4t
para t= y cuando x +el tiempo t donde anula 0 (de (1.55)). Por lo tanto si la CI se
da sobre todo IR, la onda empieza a quebrarse inmediatamente a x= +.
Las CCB para este problema son:
x<t+

(1.56)

x-t(1-2) = , = parmetro

y la envolvente de estas curvas es

1
=0
(idntica a (1.56))
4t
Puesto que dos curvas cercanas de (1.56) se cortan sobre la envolvente para valores
diferentes de , entonces u es doblemente valuada en (1.57). En efecto, la raz en (1.54)
se anula sobre (1.57) en el plano (x,t), es decir, u se divide en dos ramas dadas por
(1.54) sobre (1.57). Una idea grfica est dada en la Fig. 14.
Un ejercicio interesante es estudiar el caso cuando el dato inicial se restringe al
intervalo - < x a, a>1.

(1.57)

xt-

21

Movimiento de la onda 1-x2

Fig. 14

Ejemplo 8. Consideremos ahora el dato inicial u(x,0)= f(x) = senx. Obviamente la

solucin es u= sen(x-tu) y las CC son x-tsen - = 0, y el tiempo de quiebre (usando
1
(1.50)) resulta ser tq = .
cos
As, el menor valor de t ocurre a t = 1, es decir, tq=1. Una idea grfica est dada en la
Fig. 15.

Fig. 15 Movimiento de la onda senx.

22

NOTA: La solucin de
u
u
+ c(u)
= 0,
t
x
donde c(u) es una funcin arbitraria, no tiene un comportamiento muy diferente de la
solucin de

(1.58)

u
u
+u
=0
t
x
ya que si (1.58) se multiplica por c '(u), resulta una ecuacin para c(u) que tiene la
misma forma de (1.59) con u reemplazada por c(u).

(1.59)

Ejemplo 9. Propagacin de ondas de agua.

De la discusin de los ejemplos anteriores puede deducirse que una ecuacin
cuasilineal puede ser un buen modelo para estudiar la propagacin de ondas de agua
superficiales. En efecto, sea u(x,t) la elevacin de la superficie del agua respecto del
nivel u=0 del agua en reposo. Si suponemos que esta elevacin es pequea como
tambin ux y ut, la EDP que gobierna la propagacin de ondas es
(1.60)

3
1 + u u x + u t = 0 , < x < +, t > 0 .
2

Una CI apropiada es

(1 + cos x ), x , 0 < << 1

0, x >
Las ecuaciones caractersticas son

(1.61)

dx
3
dt
du
= (1 + u) ;
= 1;
= 0.
ds
2
ds
ds
En lugar de parametrizar la CI como la curva C, seguimos otro camino. De (1.62)
podemos escribir las ecuaciones caractersticas como

(1.62)

dt
dx
y du = 0
=
3
1
1+ u
2
Por lo tanto, sobre una caracterstica

(1.63)

(1.64)

dx
3
= 1+ u
dt
2

u = constante .

dx
tambin es constante ah. Luego las
dt
caractersticas son rectas no paralelas. De (1.64),

Como u es constante sobre una caracterstica,

23

(1.65)

x (1 + u)t = c 1

u = c2

3
sobre una caracterstica; por lo tanto, en general u = f(x-(1+ u)t).
2
De la CI (1.61) s sigue que

3
3
x 1 + u t
1 + cos( x 1 + u t ),

2
2
(1.66)
u(x,t)=
3
0 ,
x 1 + u t > .

2
3
Luego, c2 = (1+cosc1) sobre x 1 + u t = c 1 , por ejemplo. La Fig. 16 muestra las
2
3

CCB (1.65). Los valores de las constantes (c1, c2) =
,0 , (-,0), , ,
2
2
....determinan rectas caractersticas particulares.

Fig. 16 Plano caracterstico.

La Fig. 17 muestra la propagacin de los perfiles de ondas e ilustran el fenmeno no

lineal del quiebre del frente de ondas. Para terminar, calculemos el tiempo de quiebre
(por el primer mtodo)

24


sen
3 3

ux = -sen x 1 + u t 1 u x t u x =
3
2 2

1 t sen
2
3
3
con = x 1 + u t . Claramente ux = cuando 1 t sen = 0 .
2
2
Como es constante a lo largo de cualquier caracterstica, esto ocurre primero sobre la

2
recta cuando = y as sen=1. Por lo tanto t q =
.
2
3
Estudie el caso para <0.

Fig. 17 Propagacin de perfiles

Prof. Ral F Jimnez

Departamento Matemtica
Universidad de Concepcin
Marzo 2000.

25