El famoso problema del quitanieve

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:

t0h0=0t1=12:00h1x1=0t2=13:00h2x2=2t3=14:00h3x3=3\begin{array}{l|c|r} t_{0} & h_{0}=0 & \\ \hline t_{1}=12:00& h_{1} & x_{1}=0\\ \hline t_{2}=13:00& h_{2} & x_{2}=2\\ \hline t_{3}=14:00& h_{3} & x_{3}=3 \\ \hline \end{array}

Donde tt: tiempo en horas; hh: altura de la nieve; xx: distancia recorrida por el quitanieve en millas.

Podemos afirmar razonablemente que el quitanieves saca la nieve a caudal constante q˙\dot{q} y que la altura de la nieve se incrementa con el tiempo a una rapidez k˙\dot{k}, haciendo ww como el ancho del quitanieves, podemos establecer las ecuaciones diferenciales.

La ecuación para la altura de la nieve (A)(A) es válida en todo el rango t0t_{0} a t3t_{3}

dhdt=k˙      (A) \dfrac{dh}{dt}=\dot{k}\space\space\space\space\space\space(A)

Y la ecuación del quitanieve (B)(B) es válida en el rango t1t_{1} a t3t_{3}

d(whx)dt=q˙      (B)\dfrac{d\left(w\cdot h\cdot x\right)}{dt}=\dot{q}\space\space\space\space\space\space(B)

Separando variables e integrando la ecuación (A)(A) desde t0t_{0} hasta tt y sabiendo que h0=0h_{0}=0.

h0hdh=k˙t0tdt\int_{h_{0}}^{h}\mathrm{d}h =\dot{k} \int_{t_{0}}^{t}\mathrm{d}t

h=k˙(tt0)h=\dot{k} (t-t_{0})

Evaluando por separado para t2t_{2} y t3t_{3} respectivamente.

h2=k˙(t2t0)h_{2}=\dot{k} (t_{2}-t_{0})

h3=k˙(t3t0)h_{3}=\dot{k} (t_{3}-t_{0})

Dividiendo ambas encuaciones y sabiendo que t2=13t_{2}=13 y t3=14t_{3}=14 obtenemos:

h2h3=13t014t0      (C)\frac{h_{2}}{h_{3}}=\frac{13-t_{0}}{14-t_{0}}\space\space\space\space\space\space(C)

Con la ecuación (B)(B), separando variables e integrando en el intervalo t1t_{1} y t2t_{2} y luego en el intervalo t2t_{2} y t3t_{3} tenemos

Recuerde el teorema fundamental del cálculo abd(F(x))=F(b)F(a)\int _{a}^{b}\mathrm{d}(F(x))= F(b)-F(a)

wh1,x1h2,x2d(hx)=qt1t2dtw\int_{h_{1},x_{1}}^{h_{2},x_{2}}\mathrm{d}(h\cdot x)=q\int_{t_{1}}^{t_{2}}\mathrm{d}t

w(h2x2h1x1)=q(t2t1)w(h_{2}x_{2}-h_{1}x_{1})=q(t_{2}-t_{1})

w(h3x3h2x2)=q(t3t2)w(h_{3}x_{3}-h_{2}x_{2})=q(t_{3}-t_{2})

Diviendo las dos últimas ecuaciones y reemplazando x1=0x_{1}=0, x2=2x_{2}=2, x3=3x_{3}=3; t1=12t_{1}=12, t2=13t_{2}=13, t3=14t_{3}=14:

h2x2h1x1h3x3h2x2=t2t1t3t2\frac{h_{2}x_{2}-h_{1}x_{1}}{h_{3}x_{3}-h_{2}x_{2}} =\frac{t_{2}-t_{1}}{t_{3}-t_{2}}

2h20h13h32h2=13121413\frac{2 h_{2}-0 h_{1}}{3 h_{3}-2 h_{2}} =\frac{13-12}{14-13}

Reordenando convenientemente:

h2h3=34      (D)\frac{h_{2}}{h_{3}}=\frac{3}{4}\space\space\space\space\space\space(D)

Igualando las ecuaciones (C)(C) y (D)(D) y despejando t0t_{0}

34=13t014t0\frac{3}{4}=\frac{13-t_{0}}{14-t_{0}}

t0=10t_{0}=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 dxdt1h\frac{dx}{dt} \varpropto \frac{1}{h} lo cuál no tiene ningún sentido, no considera el aporte de la velocidad de incremento de altura de la nieve dhdt\frac{dh}{dt} a la velocidad del quitanieve. Una ecuación apropiada para describir el fenómeno aproximado se deriva de la ecuación (B)(B):

w(xdhdt+hdxdt)=q˙w\left( x\frac{dh}{dt}+h\frac{dx}{dt} \right) = \dot{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 hdxdt=constanteh\frac{dx}{dt} = constante lo cuál no es cierto porque sabemos que hh varía