El problema del quitanieve
Un día comenzó a nevar en forma intensa y constante. Un quitanieve comenzó a medio día, y avanzó 2 millas la primera hora y una milla la segunda ¿A qué hora comenzó a nevar?
Fuente: Ralph Palmer Agnew(1942), Differential Equations
Resolución del problema
Del enunciado anotamos nuestros datos y variables:
t0t1=12:00t2=13:00t3=14:00h0=0h1h2h3x1=0x2=2x3=3
Donde t: tiempo en horas; h: altura de la nieve; x: distancia recorrida por el quitanieve en millas.
Podemos afirmar razonablemente que el quitanieves saca la nieve a caudal constante q˙ y que la altura de la nieve se incrementa con el tiempo a una rapidez k˙, haciendo w como el ancho del quitanieves, podemos establecer las ecuaciones diferenciales.
La ecuación para la altura de la nieve (A) es válida en todo el rango t0 a t3
dtdh=k˙ (A)
Y la ecuación del quitanieve (B) es válida en el rango t1 a t3
dtd(w⋅h⋅x)=q˙ (B)
Separando variables e integrando la ecuación (A) desde t0 hasta t y sabiendo que h0=0.
∫h0hdh=k˙∫t0tdt
h=k˙(t−t0)
Evaluando por separado para t2 y t3 respectivamente.
h2=k˙(t2−t0)
h3=k˙(t3−t0)
Dividiendo ambas encuaciones y sabiendo que t2=13 y t3=14 obtenemos:
h3h2=14−t013−t0 (C)
Con la ecuación (B), separando variables e integrando en el intervalo t1 y t2 y luego en el intervalo t2 y t3 tenemos
Recuerde el teorema fundamental del cálculo ∫abd(F(x))=F(b)−F(a)
w∫h1,x1h2,x2d(h⋅x)=q∫t1t2dt
w(h2x2−h1x1)=q(t2−t1)
w(h3x3−h2x2)=q(t3−t2)
Diviendo las dos últimas ecuaciones y reemplazando x1=0, x2=2, x3=3; t1=12, t2=13, t3=14:
h3x3−h2x2h2x2−h1x1=t3−t2t2−t1
3h3−2h22h2−0h1=14−1313−12
Reordenando convenientemente:
h3h2=43 (D)
Igualando las ecuaciones (C) y (D) y despejando t0
43=14−t013−t0
t0=10
Empezó a nevar a las 10:00 en la mañana.
So that’s was all folks!
Aclaraciones finales
Nota: En muchas de las resoluciones para este problema, se sacan del bolsillo, de que la velocidad del quitanieve es inversamente proporcional a la altura de la nieve es decir dtdx∝h1 lo cuál no tiene ningún sentido, no considera el aporte de la velocidad de incremento de altura de la nieve dtdh a la velocidad del quitanieve. Una ecuación apropiada para describir el fenómeno aproximado se deriva de la ecuación (B):
w(xdtdh+hdtdx)=q˙
En la ecuación anterior se observa el aporte diferencial en dos dimensiones, notese que si reemplaza la proporcionalidad sugerida, esta ecuación pierde su sentido haciendo que hdtdx=constante lo cuál no es cierto porque sabemos que h varía