W. Luis Mochán Backal
Instituto de Ciencias Físicas, UNAM. Miembro de la ACMor
No es la primera vez; no será la última vez. En época de huracanes no es extraño que caigan lluvias torrenciales en las cuencas de nuestros grandes ríos. Previendo estos meteoros, las presas reguladoras y las hidroeléctricas deben disminuir con anticipación su nivel de agua para poder amortiguar los consecuentes incrementos en el flujo de los caudales que las alimentan. Si no, corren el riesgo de desbordarse o de sufrir daños estructurales. Ha sucedido que, ante las crecidas de los ríos, se abran de emergencia y en demasía sus compuertas, inundando las tierras bajas. Buena parte del litoral del Golfo de México está formada por planicies de poca altura sobre el nivel del mar. Por ejemplo, en la Figura 1 se muestra un mapa topográfico de la zona de Tabasco. La figura muestra grandes franjas de tierra de entre 100 km y 150 km de ancho cuya altura no sobrepasa un par de decenas de metros sobre el nivel del mar. Al tener una planicie prácticamente horizontal, no hay una pendiente que impulse al agua en exceso hacia el mar. Por lo tanto, no es infrecuente que dicha zona se inunde y que el agua permanezca durante muchos días cubriendo tierras y poblados (ver Figuras 2 y 3). Incluso para los fenómenos más evidentes, hay ciencia que nos ayuda a explicarlos, a entenderlos y posiblemente a tomar medidas pertinentes que permitan evitar potenciales desastres como lo son las inundaciones.
Tabasco, sus inundaciones y las ecuaciones diferenciales
La primera vez que viajé en Tabasco, en noviembre de 1999, gran parte de su territorio llevaba más de un mes inundado. Se decía que, por un cambio en la política de generación eléctrica en México, se había privilegiado la compra de energía a empresas privadas, bajando la producción en las presas hidroeléctricas de la CFE, las cuales acumularon tanta agua que no pudieron responder adecuadamente ante las fuertes lluvias de octubre, lo cual condujo a la tragedia. Quizás por eso y por muchos otros motivos, desde entonces ha sido recurrente la noticia de que Tabasco fue inundado una y otra vez.
La inundación de Tabasco de 2017 coincidió con el final de un curso sobre ecuaciones diferenciales (ver referencia (1)) que impartía yo en la Licenciatura en Ciencias de la UAEMor. Se me ocurrió que el estudio de las inundaciones podría ser un tema interesante para presentar a mis alumnos, pues su solución ilustra muchos de los conceptos que cubre el curso. Como no lo hice, y recientemente Tabasco sufrió de nuevo una inundación mayor, aprovecho este espacio para mostrarle a Ud., querido lector, cómo pueden estudiarse estos fenómenos matemáticamente. Antes de que huyendo abandone la lectura de este artículo, permítame decirle que no mencionaré las ecuaciones diferenciales parciales no lineales de segundo orden acopladas que hubiera deseado presentar a mis alumnos, sino que recurriré a una versión simplificada que espero pueda seguir.
La vaca esférica
Estudiar con todo detalle el problema de la dinámica hidráulica en una región geográfica es un problema complicado. Los resultados cuantitativos precisos dependen de detalles como la orografía, la textura del piso, su absorción, la evaporación, la precipitación, etc. Aquí, para obtener resultados más cualitativos y semicuantitativos, recurriré a aproximaciones, como la de la conocida vaca esférica, haciendo referencia a las historias que algunos cuentan para reírse de los físicos (ver referencia (2)).
Cuentan que un granjero, cuyas vacas dejaron de dar leche, pidió ayuda a un equipo de académicos encabezados por un físico. La respuesta decía: Hemos hallado la solución exacta, pero sólo para el caso de una vaca esférica colocada en el espacio vacío…
Pocos granjeros respetarían las predicciones de una teoría sobre la producción lechera de las vacas basada en las premisas anteriores. Los físicos solemos hacer modelos de la realidad que asemejan más una caricatura cruda en que se trazan sólo las líneas esenciales, más que un cuadro que describa todos los detalles de la escena.
El modelo que emplearé está ilustrado en la Figura 4. Una región plana horizontal de longitud L y ancho a se halla cubierta de agua hasta una altura h(t) que depende del tiempo t. En un extremo se elevan montañas que contienen al agua. En el otro, el agua se vierte al mar, desapareciendo del sistema. La longitud del sistema L es miles de veces mayor que la altura del agua h, por lo cual la figura no está dibujada a escala. Mostraremos que su ancho a es irrelevante. Si nos esperamos un pequeñísimo intervalo Δt de tiempo t que vaya de t a t+Δt, el nivel del agua subirá una pequeñísima distancia Δh de h(t) hasta h(t+Δt)=h(t)+Δh. Claro que el nivel bajará, no subirá, lo cual simplemente significa que Δh<0 es una cantidad negativa. El volumen del agua ΔVb que baja es el del prisma ancho y plano en la parte superior de la figura, de dimensiones L, a y |Δh|, ΔVb=La|Δh|=-LaΔh. Durante el mismo tiempo, un volumen Vv de agua se vierte hacia el mar gracias a la corriente que fluye hacia la derecha de la figura con velocidad v. En un intervalo de tiempo Δt el líquido avanza una distancia vΔt, por lo que el volumen del agua vertido ΔVv es el del prisma mostrado del lado derecho de la figura, ΔVv=ahvΔt. Estrictamente, en lugar del símbolo de igual (=) debí haber empleado un símbolo de aproximadamente igual (≅) y escribir ΔVv≅ahvΔt, haciendo notar que la aproximación es mejor mientras más pequeños sean Δt y Δh. La conservación de la materia nos indicaría que los dos volúmenes son iguales entre sí; el nivel del agua baja pues el volumen correspondiente es vertido hacia el mar, ΔVb=ΔVv, de donde se obtiene -LaΔh≅ahvΔt y Δh/Δt≅-hv/L. Esta aproximación es mejor mientras más pequeño sea el intervalo Δt. Quienes hayan tomado un curso de cálculo reconocerán aquí la ecuación diferencial dh/dt=-hv/L.
Energía
La ecuación diferencial previa es insuficiente para resolver el problema, pues no sabemos aún cual es la velocidad v con que fluye el agua. Podemos estimar esta velocidad de hacer algunas consideraciones energéticas. Al llevar una piedra de masa M a una altura z le proporcionamos una energía potencial gravitacional Mgz, donde g es la aceleración de la gravedad, cuyo valor es aproximadamente g≅10m/s2, como podemos verificar soltándola y observando que cae adquiriendo cada vez mayor energía cinética. Por eso, no debemos colocar nuestra cabeza debajo de piedras que podrían caer. Análogamente, volviendo a la Figura 4, vemos que la energía potencial de la capa superficial de agua de ancho Δh es ΔUg=ΔMgh=ρgLah|Δh|, donde ΔM=ρLa|Δh| es su masa y ρ es la densidad del agua, cuyo valor aproximado es ρ≅1000kg/m3. Al bajar el nivel del agua, ¿qué pasa con dicha energía potencial?
Viscosidad
Para contestar, tome un huevo crudo entre sus dedos, colóquelo sobre una mesa y póngalo a girar como si fuera un trompo. Notará que ¡es imposible! El huevo detendrá inmediatamente su rotación. Sin embargo, si repitiera la experiencia con un huevo duro, notaría que el huevo mantendría su giro un tiempo relativamente largo. ¿Por qué se detiene el huevo cuando está crudo? Porque en su interior contiene un líquido viscoso (ver referencia (3)). Al pretender poner al huevo en rotación, en realidad ponemos a girar sólo al cascarón. El líquido en su interior se queda inicialmente en reposo. Lentamente, el ímpetu del ligero cascarón empieza a difundirse hacia el interior, poniendo en movimiento a la relativamente pesada clara y yema. De acuerdo a las leyes de conservación del ímpetu y del ímpetu angular, conforme el interior del huevo empieza a moverse, el cascarón se frena. La fricción entre distintas capas internas que giran con distinta velocidad acaba disipando la energía y convirtiéndola en calor (la clara del huevo es un poco más complicada, pues más que viscoso, es un líquido viscoelástico, ver referencia (4)). En cambio, un huevo duro actúa como un sólido.
Si una capa de un líquido se mueve más rápidamente que otra capa contigua, le transfiere gradualmente su cantidad de movimiento a través de su frontera en proporción al área, la diferencia de velocidades y en proporción inversa al ancho de las capas (ver Figura 5). Similarmente, al fluir sobre el fondo, una corriente de agua siente una fuerza que la frena, disipando su energía. La fuerza está dada por F≅ηLaΔv/Δz, donde Δv es el cambio en la velocidad cuando ascendemos una distancia Δz y η es el coeficiente de viscosidad, que para el agua toma el valor η≅10-3kg/m s. Para quienes saben cálculo, la expresión correcta es F≅ηLa∂v/∂z. Aquí debemos notar que la velocidad es en realidad un campo de velocidades que dependen de la posición y del tiempo.
Una fuerza F que actúa cierta distancia Δx produce un trabajo ΔW=FΔx. De acuerdo a la Figura 4, el agua avanza una distancia Δx=vΔt durante un tiempo Δt. Luego, el trabajo que se realiza sobre el agua, la energía que disipa conforme fluye, es ΔW=FΔx≅ηvLaΔt(Δv/Δz). Aproximando (Δv/Δz)≅v/h, la energía disipada resulta ser ΔW≅ηv2 LaΔt/h. Esta energía se convierte en calor.
Balance de energía
¿De dónde viene esta energía que se convierte en calor? De la energía potencial gravitacional que pierde el líquido conforme baja su nivel. Igualando la energía disipada con el cambio de energía potencial gravitacional del agua llegamos a una ecuación para el cambio del nivel del agua conforme transcurre el tiempo, Δh/Δt=-ρg h4/ηL2. Esta ecuación, válida en el límite de intervalos de tiempo muy pequeños, es en realidad una ecuación diferencial dh/dt=-ρg h4/ηL2 fácilmente soluble. La solución es t=(η/3ρg)L2(1/h3-1/h03), donde h0 es la altura inicial. El lector interesado puede verificar este resultado, aún sin saber cálculo, usando la ley del binomio 1/(h+Δh)3=(h+Δh)-3≅h-3-3h-4Δh+… donde los términos que faltan son despreciables por ser pequeños cuando Δh→0.
¿Por qué Tabasco permanece inundado tanto tiempo?
Para ilustrar el resultado anterior, consideremos un ejemplo que podría aplicarse a regiones de Tabasco. Consideremos una planicie costera con una extensión de 1km inundada hasta una altura h0=1m. Suponiendo que ya ha dejado de llover y que no hay corrientes tierra arriba que traigan agua adicional, ¿cuánto tiempo t tardaría el agua en bajar hasta una altura de, digamos, h=5cm? Primero notamos que 1/h3=8000 m-3, mientras que 1/h03=1 m-3, por lo cual 1/h3-1/h03=7999 m-3≅8000 m-3. Esto ilustra un resultado curioso: Si la altura inicial es varias veces mayor que la altura final, la duración no dependerá prácticamente de la altura inicial, sólo de la altura final. Esto se debe a que en la fase inicial el nivel del agua baja rápidamente, para estancarse mucho tiempo cerca de su nivel final. Entonces podemos aproximar t≅(η/3ρg)(L2/h3). Sustituyendo los valores apropiados para el agua podemos evaluar la constante
η/3ρg≅10-3/(3⋅103⋅10)s/m=3⋅10-8s/m. Como L2=1,000,000m2, nuestra fórmula lleva a
t=3⋅10-8⋅106⋅8⋅103s=240s, 4 minutos. Sin embargo, si la planicie midiera 10km, el tiempo aumentaría a t=24,000s, poco menos de siete horas. Y si fuese más grande aún, una planicie de 100km, el tiempo se incrementaría a t=2,400,000s, ¡casi un mes!
Un poco de autocrítica
El modelo empleado arriba tiene varios errores conceptualmente graves. Por un lado, no toma en cuenta que la velocidad del líquido en el extremo de la planicie, junto a las montañas, es necesariamente cero. Por lo tanto, la velocidad debe depender de x. El flujo es más veloz cerca de la orilla que cerca de la montaña. Por otro lado, hay una fuerza que empuja al agua hacia el mar, la cual proviene de una caída de la presión (la presión debe ser mayor tierra adentro que cerca de la orilla). Eso implica que el nivel h es mayor tierra adentro que cerca del borde, es decir, h también es un campo que depende de la posición además de depender del tiempo. Entonces, la ecuación diferencial que obtuvimos no puede ser correcta y debe ser reemplazada por una ecuación en derivadas parciales (ver referencia 5) que mezcle de forma adecuada la dependencia espacial con la dependencia temporal. Resulta sorprendente que, aun incluyendo estas dificultades, el problema pudo resolverse de manera exacta. Más sorprendente aún es que los resultados mencionados arriba prácticamente no cambian.
Otro tipo de complicaciones tienen que ver con los detalles de la orografía, con los ríos que surcan la planicie, con las corrientes de agua provenientes de río arriba y con la presencia de lluvias continuas durante la inundación. El flujo del agua puede además ser turbulento cuando su velocidad es alta, con vórtices caóticos que disipan la energía mucho más eficientemente que el flujo laminar empleado en nuestro cálculo, y que por tanto reducen la velocidad y aumentan la duración de las fases iniciales de la inundación. Finalmente, el modelo tampoco toma en cuenta que el agua se evapora y se filtra. Incluir todas estas consideraciones complicarían enormemente el modelo y requeriría datos geográficos detallados y modelos matemáticos complejos que sólo se pueden resolver mediante uso masivo de computadoras digitales.
Conclusiones
Consideraciones relativamente sencillas como son la ley de conservación de la materia y la ley de la conservación y disipación de la energía permiten elaborar modelos relativamente simples de problemas complejos, como el de las inundaciones, cuya comprensión puede tener una enorme relevancia social. Estos modelos pueden resolverse para obtener predicciones específicas mediante el uso de algunas técnicas matemáticas simples. En particular, se obtuvo que el tiempo de permanencia del agua en una planicie es proporcional al cuadrado de su tamaño y aproximadamente proporcional al inverso del cubo de la altura final, y prácticamente no depende de la altura inicial. Estos cálculos y conceptos pueden hacer la diferencia entre generar o evitar una catástrofe.
Esta columna se prepara y edita semana con semana, en conjunto con investigadores morelenses convencidos del valor del conocimiento científico para el desarrollo social y económico de Morelos. Desde la Academia de Ciencias de Morelos externamos nuestra preocupación por el vacío que genera la extinción de la Secretaría de Innovación, Ciencia y Tecnología dentro del ecosistema de innovación estatal que se debilita sin la participación del Gobierno del Estado.
Agradecimientos
Este trabajo se hizo con apoyo parcial de la UNAM a través del proyecto IN111119 de DGAPA-PAPIIT.
Referencias
- Ecuación diferencial – Wikipedia, la enciclopedia libre, https://bit.ly/33sDwzr
- Vaca esférica – Wikipedia, la enciclopedia libre https://bit.ly/3fOtFZs
- Viscosidad – Wikipedia, la enciclopedia libre https://bit.ly/2HSHud1
- Viscoelasticidad – Wikipedia, la enciclopedia libre https://bit.ly/39ujurY
- Derivada parcial – Wikipedia, la enciclopedia libre, https://bit.ly/33xvIfD
- Mapa topográfico de Tabasco https://bit.ly/37jbJ5n
- Cambio climático y corrupción, causas de las inundaciones en Tabasco – Firmas.Mx https://bit.ly/3fIQq13
- México: esto sucede cuando se tumba selva en el territorio maya, Animal Político, 28 de noviembre, 2020. https://bit.ly/2JuhdCg
Fuente: Academia de Ciencias de Morelos