Earth
Los glaciares modelan valles
Llena el valle en V de un río con un glaciar de seiscientos metros de espesor y mira pasar cien mil años. El hielo fluye más deprisa en el centro de su superficie y nada contra el lecho; el esfuerzo que arrastra sobre el lecho es máximo en las laderas de la V y se anula en su fondo estrecho, así que el deslizamiento y la abrasión recortan las laderas y el valle se ensancha hacia una U - la V antigua queda dibujada para comparar. El hielo limpio corta mucho menos que el cargado de roca; sin deslizamiento no corta nada; haz avanzar y retroceder el glaciar en ciclos. Cada metro cúbico que pierde el lecho va a la morrena. En un canal semicircular la solución es la de Poisson cerrada.

Qué ocurre
Una sección de valle de 2400 m de ancho en una malla de 61 por 40. La velocidad valle abajo del hielo newtoniano (10¹³ Pa·s, pendiente 0,05) resuelve eta grad²u = −rho g sin alpha con u = 0 en el lecho y du/dz = 0 en la superficie; el lecho irregular se trata con diferencias de Shortley-Weller y la solución por sobrerrelajación desde la anterior. El esfuerzo sobre el lecho impulsa un deslizamiento proporcional - una simplificación declarada, el deslizamiento no entra en la solución - y el deslizamiento desgasta el lecho (Hallet), escalado por los derrubios del hielo y la resistencia de la roca. Cada metro cúbico va a la morrena. Pasos de 1 milenio; un segundo son mil años a escala de tiempo 1.
Qué probar
Mira la V ensancharse hacia una U y lee el exponente del perfil; compara hielo limpio y sucio; detén el deslizamiento; haz avanzar y retroceder el glaciar; comprueba el semicírculo contra la fórmula cerrada.
El modelo y sus límites
Una sección sin flujo a lo largo del valle, hielo newtoniano, deslizamiento proporcional al esfuerzo calculado, una sola ley de abrasión; no se modelan el arranque ni la fractura de la roca. Los perfiles son de este valle idealizado.
Ecuaciones del módulo
Ecuaciones η·∇²u = −ρ·g·sin α, u = 0 on the bed, ∂u/∂z = 0 at the surface τ_b = η·|∇u| at the bed, u_b = c_s·τ_b dz_b/dt = −K·D/S·u_b·√(1 + (∂z_b/∂y)²) semicircle: u = ρ·g·sin α·(R² − r²)/(4η), τ_b = ρ·g·sin α·R/2


