Earth
Mareas en un estuario
Mira la marea de la tarde subir por un estuario en embudo junto al faro: el cauce se estrecha hacia el río, así que la onda crece al avanzar y la carrera en la cabecera es mucho mayor que en la boca, la llenante dura menos que la vaciante, y el agua salada del mar penetra decenas de kilómetros tierra adentro contra el río. Sube el caudal del río y la sal retrocede hacia el mar; afila el embudo o enderézalo; cierra la cabecera y el cauce resuena como un cuarto de onda, con la marea de la cabecera a / cos(kL) de la del mar. El agua quieta y cerrada sigue exactamente quieta, un estuario salado sin río sigue exactamente salado, y agua y sal cuadran con lo que cruzó la boca y la cabecera.

Qué ocurre
Un cauce de 60 km y 8 m de profundidad cuya anchura se estrecha desde 3 km como exp(−x/25 km), en 120 celdas sobre una malla escalonada: la superficie en los centros, la velocidad en las caras. Pasos adelante-atrás de 20 s actualizan la superficie con los caudales de las caras y luego la velocidad con la nueva pendiente, con advección descentrada aguas arriba y fricción de Manning implícita. La marea M2 del mar de 1,5 m entra por la boca a lo largo de dos periodos; el río aporta 100 m³/s en la cabecera. La sal viaja con los mismos caudales, descentrada aguas arriba, con una dispersión mareal de 150 m²/s, así que un estuario salado sin río sigue salado y el agua quieta y cerrada sigue quieta. Un cauce prismático cerrado comprueba el cuarto de onda a/cos(kL).
Qué probar
Mira crecer la carrera embudo arriba y la llenante durar menos que la vaciante; sube el río y empuja la sal fuera; endereza o afila el embudo; cierra la cabecera y compara con el cuarto de onda.
El modelo y sus límites
Un cauce unidimensional bien mezclado de profundidad uniforme, sin secado de las llanuras mareales, salinidad sin efectos de densidad, una sola componente de marea. Las carreras son las de este estuario idealizado.
Ecuaciones del módulo
Ecuaciones b·∂η/∂t + ∂(b·H·u)/∂x = 0, b = b0·e^(−x/L_b) ∂u/∂t + u·∂u/∂x = −g·∂η/∂x − g·n²·|u|·u/H^(4/3) ∂(b·H·S)/∂t + ∂(b·H·u·S − b·H·K·∂S/∂x)/∂x = 0 closed prismatic channel: a_head = a/cos(kL), k = ω/√(g·h)


