Aplicaciones Ecuaciones Diferenciales
Aplicaciones Ecuaciones Diferenciales
Aplicaciones Ecuaciones Diferenciales
orden.
Prof. Leonardo Sáenz Baez
LEY DE ENFRIAMIENTO DE NEWTON.
El enunciado de la ley de Newton para el enfriamiento de un cuerpo es:
La razón o velocidad con que cambia la temperatura de un cuerpo, con
respecto al tiempo, es proporcional a la diferencia de temperaturas entre el
cuerpo y el medio ambiente
Si denotamos con T a la temperatura del cuerpo en un tiempo t cualquiera y con T m a
la temperaura del medio ambiente, entonces la ley de enfriamiento se puede formular
mediante la siguiente ecuación diferencial
dT −kT − T m
dt
donde k es una constante de proporcionalidad positiva y T T m , el signo menos en la
ecuación es necesario ya que la temperatura es una función decreciente y la derivada
debe ser negativa. La T , es la variable dependiente y t la variable independiente, T m y
k son constantes.
Esta ecuación diferencial se puede escribir
dT kT kT m
dt
la cual es una ecuación diferencial lineal de primer orden, cuya solución, como ya
hemos visto se puede encontrar escribiendo:
pt k
gt kT m
siendo el factor de integración
t e
ptdt e kdt e kt
y por lo tanto la solución es
T e −kt e kt kT m dt Ce −kt T m e −kt e kt Ce −kt
T T m Ce −kt
Ahora bien, cuando t 0 la temperatura tiene un cierto valor inicial T i , esta condición
nos permite encontrar el valor de C
Ti Tm C C Ti − Tm
de manera que la solución se puede escribir
T T m T i − T m e −kt
debe hacerse notar que cuando t → un tiempo muy grande , e −kt → 0 y , entonces,
T → T m , lo cual significa que para un tiempo largo, la temperatura del cuerpo en
enfriamiento se aproxima a la temperatura del medio ambiente.
1
La ecuación diferencial
dT −kT − T m
dt
tambien es separable y se puede escribir como
dT −kdt
T − T m
integrando ambos lados se tiene
lnT − T m −kt C
y tomando antilogaritmos se escribe
T − T m e −ktC e C e −kt C 1 e −kt
T T m C 1 e −kt
y si para t 0, T T i , entonces
Ti Tm C1 C1 Ti − Tm
y la solución será
T T m T i − T m e −kt
donde k , la constante de proporcionalidad se puede determinar conociendo la
temperatura T en algún otro valor del tiempo t ≠ 0, digamos T T 1, para cuando
t t 1 ≠ 0. , en cuyo caso tendriamos
T 1 T m T i − T m e −kt 1
e −kt 1 T 1 − T m
Ti − Tm
− kt 1 ln T 1 − T m
Ti − Tm
k − 1 ln T 1 − T m 1 ln T i − T m
t1 Ti − Tm t1 T1 − Tm
Ejemplo concreto.
Una barra metálica a una temperatura de 100 0 F se pone en un cuarto a una
temperatura constante de 0 0 F. Si después de 20 minutos la temperatura de la barra
es de 50 0 F, hallar:
a) El tiempo que que debe transcurrir para que la barra adquiera una temperatura de
25 0 F y
b) La temperatura de la barra después de 10 minutos.
Solución.
Los datos del problema son:
T i 100 0 F, Tm 00F
T 1 50 0 F para t 1 20min
por lo tanto
2
k 1 ln 100 − 0 1 ln2 ln2 20
1
20 50 − 0 20
1 t
kt t ln2 20 ln2 20
y la solución será
t
T 0 100 − 0e −ln2 20
T 100 1 t 100 1 t 1002 − 20
t
e ln2 20 2 20
T 1002 − 20
t
a) Si T 25, entonces
25 1002 − 20
t
25 2 − 20t
100
1 2 − 20t
4
tomando logaritmo natural
ln 1 − t ln2
4 20
− ln4 − t ln2
20
2 ln2 t ln2
20
2 t
20
t 40 min
Después de transcurridos 40 min. la barra tendrá una temperatura de 25 0 F.
b) Si t 10, entonces
T 1002 − 20 1002 − 2 100 50 2 ≈ 70. 71 0 F
10 1
2
T ≈ 70. 71 0 F
La temperatura después de transcurridos 10 min sera aproximadamente de 70. 71 0 F.
HORA DE UN DECESO.
En la investigación de un homicidio o de una muerte accidental, a menudo es
importante estimar el momento de la muerte.
Supongamos que en t 0 se descubre un cadáver y que su temperatura es T i , se
supone que en el instante t d , del fallecimiento o del deceso, la temperatura del
cuerpo T d tenía la temperatura normal de 37 0 C 98. 6 0 F.
Considerando que se cumple la ley de enfriamiento de Newton
dT −kT − T m
dt
cuya solución es
3
T T m T i − T m e −kt
donde T m es la temperatura del medio ambiente y T i es la temperatura inicial cuando
t 0.
donde k se determina con
k 1 ln T i − T m
t1 T1 − Tm
siendo t 1 otro tiempo distinto de cero y T 1 su temperatura correspondiente
Entonces empleando
T T m T i − T m e −kt
que al despejar t se tiene
T − T m e −kt
Ti − Tm
t − 1 ln T − T m
k Ti − Tm
t d − 1 ln T d − T m
k Ti − Tm
con esta expresión calculamos t d cuando T T d 37, y k se calcula con la expresión
ya indicada anteriormente.
Por ejemplo, supóngase que la temperatura del cadáver es de 29. 4 0 C, cuando es
descubierto y dos horas mas tarde su temperatura es de 23. 3 0 C, siendo la
temperatura ambiente de 20 0 C. ¿Cuál fué la hora del deceso?.
Temperatura inicial T i 29. 4 0 C
Temperatura ambiente T m 20 0 C
Solución. Los datos que se tienen son Temperatura medida T 1 23. 3 0 C
tiempo de medición t 1 2 horas
Temperatura estándar T d 37 0 C
con estos valores se tiene
k 1 ln T i − T m 1 ln 29. 4 − 20
t1 T1 − Tm 2 23. 3 − 20
k ≈ 0. 523 393 611
t d − 1 ln T d − T m − 1 ln 37 − 20
k Ti − Tm . 523 393 611 29. 4 − 20
t d ≈ −1. 132 042 2
esto es, la muerte se estima que ocurrió 1 h con 8 minutos antes de encontrar el
cadáver.
MODELOS DE CRECIMIENTO.
Si N Nt es una población que crece (o decrece) al transcurrir el tiempo y si su
velocidad de crecimiento es proporcional a la población presente en un tiempo t ,
4
entonces se puede escribir
dN kN
dt
donde k es una constante positiva y el signo de la ecuación debe ser positivo para
fenómenos de crecimiento y negativo para los casos de decrecimiento.
Estamos asumiendo que Nt es una función derivable y por lo tanto continua, sin
embargo para problemas de población, donde Nt es un entero y discreto, esta
suposición es incorrecta, sin embargo la ecuación diferencial dá una buena
aproximación a las leyes de crecimiento que rigen tal sistema (modelo exponencial).
Tratemos uno de los casos digamos el de crecimiento cuya ecuación es
dN kN
dt
esta ecuación es separable y se puede escribir
dN kdt
N
integrando
lnN kt C
tomando antilogarítmos
N e ktC e C e kt C 1 e kt
N C 1 e kt
si t 0 se tiene un cierto valor de N N 0 llamado población inicial, entonces
N0 C1
y la ecuación quedara
N N 0 e kt
en el caso de fuera un fenómeno de decrecimiento, la ecuación quedaría
N N 0 e −kt
Ejemplo. Se sabe que la población de cierto país aumenta a una velocidad
proporcional al número de habitantes en un tiempo t. Si después de 2 años la
población se ha duplicado y después de tres años la población es de 20 000
habitantes, hallar el número de habitantes que había inicialmente en el país, esto es
cuando t 0.
Sol. El modelo es de crecimiento por lo que la ecuación que lo modela es
N N 0 e kt
Para t 2 sabemos que N 2N 0 de manera que
2N 0 N 0 e 2k
2 e 2k
1
2 2 ek
con lo cual la ecuación quedaría
5
t 1 t
N N 0 e k N 0 2 2
t
N N02 2
tanbién sabemos que para t 3, N 20000 de manera que
3
210 4 N 0 2 2
210 4 4 10 4 2
N0 3 10 5 2 10 3
22 2 2
N 0 ≈ 7071
POBLACIONES LIMITADAS.
En el modelo exponencial o de Malthus
dP kP
dt
cuya solución es
P P 0 e kt
ocurre cuando en una población se considera que los índices de natalidad y
mortalidad son constantes, ya que en tales casos la tasa de crecimiento por
individuo sería
1 dP − k
P dt
dP kP
dt
Si se considera que el índice de nacimientos es una función lineal decreciente
b − mP y el índice de mortalidad es constante se tiene
1 dP − b − mP − b − − mP a − mP
P dt
dP Pa − mP
dt
Ecuación que es un modelo de crecimiento, llamado Modelo logístico de Poblaciones
propuesto por el Matemático y Biólogo Belga P.F. Verhulst.
La solución de esta ecuación separable será
dP dt
Pa − mP
Pa dP
− mP
dt C
6
1
a
1
P
m
a − mP
tC
1
a dP − −mdP
P a − mP
tC
Pa − mP 0
ln at
P 0 a − mP
Pa − mP 0
e at
P 0 a − mP
Pa − mP 0 P 0 a − mPe at
Pa − mP 0 aP 0 e at − mP 0 Pe at
Pa − mP 0 mP 0 e at aP 0 e at
P aP 0 e at
a − mP 0 mP 0 e at
P aP 0
a − mP 0 e −at mP 0
Esta solución P Pt se conoce como la función logística, la cual está acotada ya
que
a
lim Pt m
t→
7
En este tipo de curvas cuando t → , P → ma , y si t → −, P → 0, estas curvas tienen
un punto de inflexión, donde la concavidad cambia de signo, en el cual la velocidad de
crecimiento adquiere su valor máximo, dicho punto se encuentra pues, en el máximo
de la función
dP Pa − mP
dt
y como esta representa una parábola que se abre hacia abajo, y corta al eje
a , su máximo se encuentra en el vértice o sea cuando P a ,
horizontal en 0 y en m
2m
que es la mitad de la población máxima.
200
y
150
100
50
-5 -4 -3 -2 -1 0 1 2 3 4 5
x
8
dP dt
Pa − m ln P
dP
P
dt, integrando quedaría
a − m ln P
1 lna − m ln P t C
− m
y como para t 0, P P 0 se tiene
C −m 1 lna − m ln P 0
de manera que
1
m lna − m ln P 0 − lna − m ln P t
a − m ln P 0
ln mt
a − m ln P
tomando antilogarítmos
a − m ln P 0
e mt
a − m ln P
a − m ln P 0 e −mt a − m ln P
m ln P a − a − m ln P 0 e −mt
ln P ma − a − ln P 0 e −mt
m
a − a −ln P e −mt
P eme m 0
a a
y cuando t → ; P → e m , lo cual nos indica que la población está limitada por e m .
¿Para que valor de P se tiene el punto de inflexión?
En este caso como
dP Pa − m ln P
dt
se tiene
′
d 2 P P−m P a − m ln PP ′
dt 2 P
d 2 P P ′ a − m ln P − m Pa − m ln Pa − m ln P − m
dt 2
al igualar la segunda derivada con cero tenemos que los posibles valores de P en los
que puede haber puntos de inflexión son:
P 0, que se descarta por ser una asíntota horizontal
a
a − m ln P 0 P e m , que se descarta por ser también asíntota horizontal
a−m
a − m ln P − m P e m
a−m
el punto de inflexión se encuentra entonces para cuando P e m .
9
APLICACIONES EN LA MECÁNICA.
Descripción del movimiento de los cuerpos. En un cuerpo en movimiento se
encuentran involucrados ciertos parámetros fundamentales como son: la velocidad
v del cuerpo en un tiempo t cualquiera, el desplazamiento x en cualquier tiempo o
bién la aceleración a del cuerpo. En la descripción del movimiento se emplearán
principios básicos de la mecánica, como lo es la segunda ley de Newton que
establece que la fuerza F que actúa en un cuerpo de masa m es directamente
proporcional a la variación instantánea de la cantidad de movimiento la cual se define
como el producto de la masa por la velocidad, esto es:
F k d mv
dx
donde k es una constante de proporcionalidad que empleando un sistema adecuado
de unidades tendremos k 1, quedando entonces la segunda ley como
F d mv
dx
o si la masa es constante
F m dv ma
dx
donde
a dv
dt
es la aceleración del cuerpo o variación instantánea de la velocidad con respecto al
tiempo, o bien, como la velocidad es la variación instantánea del desplazamiento se
tiene
2
a dv d dx d 2x
dt dt dt dt
y como
a dv dv dx v dv
dt dx tt dx
tambien se puede escribir la segunda ley de Newton como
F ma mv dv
dx
Ejemplo 1. (Descripción del movimiento uniformemente acelerado).
Describiremos el movimiento de un cuerpo de masa m el cual se mueve con una
aceleración constante a, siendo v 0 su velocidad inicial (al tiempo t 0). Supondremos
que el cuerpo parte del reposo (el desplazamiento inicial x 0 0) y que la dirección
positiva es la usual para el eje-x.
Solución. A un tiempo t cualquiera, la única fuerza que actúa en el cuerpo es ma , de
modo que
10
m dv ma
dt
dv a
dt
dv adt
v at C
y como para t 0, v v 0 tendremos
C v0
∴ v v 0 at
o si se emplea la integral definida se puede escribir
m dv ma
dt
dv a
dt
dv adt
v t
v dv a
0
0
v − v 0 at
v v 0 at
y como v dx
dt
dx v 0 at
dt
x t
0 dx 0 v 0 atdt
x v 0 t 1 at 2
2
si consideramos la velocidad al cuadrado se tiene
v 2 v 0 at 2 v 20 2av 0 t a 2 t 2 v 20 2av 0 t 1 at 2 v 20 2ax
2
v v 0 2ax
2 2
11
m dv −mg
dt
dv −g
dt
v dv −g
dx
vdv −gdx
Como en x máx tenemos que v 0, entonces
0 x máx
v vdv −g
0
dx
0
− 1 v 20 −gx máx
2
v2
x máx 0
2g
Por otro lado como
dv −gdt
se tiene que
v t
v dv −g dt
0
0
v − v 0 −gt
v v 0 − gt
Cuando v 0 obtenemos el tiempo de ascenso
t A vg0
y en este momento el cuerpo deja de moverse y principia a descender.
Para calcular el tiempo de descenso t D hagamos lo siguiente. Sabemos que en
suretorno el cuerpo recorre una distancia
v2
x máx 0
2g
habiendo comenzado su movimiento con una velocidad nula, entonces como
dv −gdt
v t
0 dv −g 0 dt
v −gt
dx −gt
dt
de donde
12
dx −gtdt
0 tD
x dx −g tdt
0
máx
g
− x máx − t 2D
2
2
v0 g 2
t
2g 2 D
v 20
t 2D
g2
∴ t D vg0 t A
esto es, la velocidad con que se llega al punto de partida es la misma velocidad inicial
v 0 , el signo menos es por el sentido positivo que seleccionamos hacia arriba.
Ejemplo 3. Un cuerpo de masa m es lanzado verticalmente hacia arriba con una
velocidad v 0 . La resistencia del aire es proporcional a la velocidad del cuerpo.
Encuentre la altura máxima alcanzada y el tiempo que tarda en alcanzar tal altura.
El problema es como el anterior salvo que ahora se toma en cuenta la resistencia que
ofrece el aire al movimiento del cuerpo.
Durante el ascenso, las fuerzas que actúan sobre el cuerpo son W mg y R kv,
ambas en sentido negativo, de manera que la ecuación de movimiento será
m dv −mg − kv
dt
dv k v −g
dt m
k y ft −g
la cual es una ecuación lineal, para la cual pt m
siendo su factor de integración
t e
ptdt e mk t
siendo la solución
v e − m t e m t −gdt Ce − m t
k k k
13
mg − mk t mk t
e e Ce − m t
k
v−
k
mg
Ce − m t
k
v−
k
y como para t 0, v v 0
mg
v0 C
k
de modo que
mg mg − mk t
v− v0 e
k k
el tiempo de ascenso, para el cual se tiene que v 0, será
mg mg − mk t
0− v0 e
k k
mg mk t mg
e v0
k k
k t
e m mg k v 0 mg kv 0 1
k mg
k kv 0
m t ln mg 1
t A m ln kv mg 1
0
k
Para calcular la altura máxima pongamos
dv k v −g
dt m
v dv m k v −g
dx
vdv m k v g dx 0
vdv dx 0
mk v g
m m v g − g dv dx 0
k
k mk v g
m 1− g
dv dx 0
k vg
k
m
m v − m mg ln kv mg x C
k k k m
m v − mg ln kv mg x C
k k m
y como para x 0 se tiene v v 0
mg kv mg
C m v0 − ln 0 m
k k
de manera que
14
mg kv mg mg kv mg
x m v0 − ln 0 m − m v− ln m
k k k k
mg kv mg mg kv mg
x m v0 − ln 0 m − v− ln m
k k k
mg kv mg kv mg
x m v 0 − v ln m − ln 0 m
k k
mg kv mg
x m v 0 − v ln
k k kv 0 mg
cuando v 0, x x máx
mg mg
x máx m v 0 ln
k k kv 0 mg
Tarea. Calcular el tiempo de descenso t D .
Veamos como calcular el tiempo de descenso.
Para evitar problemas con los signos, cambiaremos el sistema de referencia,
tomando ahora t 0 y x 0 en el punto máximo de la trayectoria y consideraremos
la direccion positiva del movimiento hacia abajo.
Durante el descenso las fuerzas que actuan son la resistencia del aire R kv, opuesta
al movimiento y la del peso del cuerpo W mg , en favor del movimiento De manera
que la ecuación de movimiento será:
m dv mg − kv
dt
separando variables
mdv dt
mg − kv
integrando
v t
0 mdv
mg − kv
0 dt
− m lnmg − kv| v0 t
k
− m lnmg − kv − lnmg t
k
mg − kv
− m ln mg t
k
mg − kv k t
ln mg −m
mg − kv
e− m t
k
mg
despejando a la v se tiene
mg − kv mge − m t
k
mg1 − e − m t kv
k
mg
1 − e − m t
k
v
k
15
y como v dx se tiene
dt
mg
1 − e − m t dt
k
dx
k
e integrando
x mg t
0 dx
k 0
1 − e − m t dt
k
mg − k t
x t m e m | t0
k k
mg m − k t
x t e m − m
k k k
en esta ecuación cuando
x x max , se tiene que t t D
o sea
mg − k t
x max tD m e m D − m
k k k
esta ecuación es necesario resolverla numéricamente, pero de cualquier modo es
evidente que no se cumple para t D t A ,
por lo que concluimos que t D ≠ t A .
En efecto, si t D t A se tendría
kv
mg m
x max ln kvmg
0
1 m e − mk mk ln mg0 1 − m
k k k k
mg m −1
x max ln kvmg
0
1 m kv 0 1 − m
k k k mg k
m2g kv 0 mg mg
x max ln mg −1
0 mg
2 kv
k
m2g mg mg
x max − ln −1
k 2 kv 0 mg kv 0 mg
m2g mg kv 0
x max − ln −
k 2 kv 0 mg kv 0 mg
m2g mg m2g kv 0
x max − 2 ln − 2
k kv 0 mg k kv 0 mg
m2g mg m 2 gv 0
x max − 2 ln −
k kv 0 mg kkv 0 mg
pero, en realidad
m2g mg
x máx m v 0 2 ln
k k kv 0 mg
de manera que
tD ≠ tA
Si tomamos un caso concreto con m 1 kg, v 0 10 m/seg, g 9. 8 m/seg 2 ,
entonces manejando algunos valores de k, que se aproximen a cero, tendremos
16
k x max tA tD vf
1 3. 107664 0. 703299 0. 917653 5. 8853
0. 5 3, 84 0. 82449 0. 95566 7. 4455
0. 1 4. 77952 0. 97164 1. 00442 9. 36317
0. 01 5. 0676 1. 01523 1. 01189 9. 9345
donde
mg
1 − e − m t D
k
vf
k
En la tabla se puede notar que cuando k → 0, t A → t D y v f → v 0
De manera que si no hay resistencia del aire (k 0), entonces t A t D y v f v 0 .
PROBLEMA DE MEZCLAS
Problema. El lago Erie tiene un volumen de 458 km 3 y el flujo de entrada y salida se
realizan ambos a razón de 175 km 3 por año. Suponga que para el tiempo t 0 (años)
su concentración de contaminación es de 0.05% y que un tiempo después la
concentración de contaminantes que ingresa en el agua es de 0.01%. Suponiendo
que el agua se mezcla perfectamente dentro del lago, calcule el tiempo que debe
transcurrir para que la concentración de contaminantes se reduzca al 0.02%.
17
la cual es una ecuación diferencial lineal para la cual pt 0. 3821, ptdt 0. 3821t
siendo el factor de integración
t e 0.3821t
y por lo tanto la solución es
xt e −0.3821t e 0.3821t 0. 0175dt Ce −0.3821t
Solución.
dx 500. 0005 − x 500 . 25 −. 0 625x
dt 8000
dx . 0 625x . 25
dt
la ecuación diferencial es lineal con
18
pt 0. 0625
t e
0.0625dt e 0.0625t
x 0. 25 Ce −0.0625t 4 Ce −0.0625t
0. 0625
x 4 Ce −0.0625t
si t 0, x 0. 00258000 20. 0
de manera que
20 − 4 C
C 16
x 4 16e −0.0625t
ahora que tiempo debe transcurrir para que
x 0. 0018000 8
8 4 16e −0.0625t
resolviendo para t se tiene
1 e −0.0625t
4
ln 1 −0. 0625t
4
t ln 4 22. 180 709 8
0. 0625
t ≈ 22. 18 dias
19