# Revisión técnica ronda A1 — física y matemática (nivel publicación)

Eres revisor de un artículo de hidráulica de flujos de detritos. Revisa SOLO física y matemática del modelo (no estilo de código).
Caso: aluvión del 8-oct-2026 en San Carlos de Apoquindo (Las Condes, Santiago). Lee `sim.py` (modelo) y `prep.py` (terreno).

Modelo implementado:
- ANUGA 4.0.1 (aguas someras 2D, volúmenes finitos, algoritmo DE1), mezcla homogénea "bulked": se conserva el volumen de mezcla.
- Fricción de mezcla por paso fraccionado semi-implícito (operador `Reologia`): S_f = τy/(γm h) + K η |u|/(8 γm h²) + n_td² |u|²/h^{4/3}
  con τy = α2 e^{β2 Cv}, η = α1 e^{β1 Cv} (O'Brien & Julien 1988; tabla FLO-2D: glenwood1/2, aspen1/2), ρm = ρw + Cv(ρs−ρw).
  Actualización: |u|' = max(0, |u| − Δt g τy/(γm h)) / (1 + Δt g [K η/(8 γm h²) + n² |u|/h^{4/3}]).
- Fricción Manning de ANUGA = 0. Entrada: hidrograma triangular de mezcla (Qp, volumen V, tp) en un círculo r = 14 m en el cauce abierto
  aguas arriba de Av. La Plaza (supuesto: el ducto subterráneo quedó obstruido y todo el caudal escurrió por superficie).
- Terreno: FABDEM (DTM 30 m) remuestreado a 4 m; edificios de OSM elevados como bloques (altura por DSM−DTM, building:levels, o 6 m);
  calles rebajadas 0,15 m. Malla: 60–120 m² en el corredor, 600 m² fuera.
- Calibración prevista: barrido de Qp, V, Cv y suelo; métrica = coincidencia de la zona con h > 5 cm con el polígono municipal
  (El Alba / La Plaza / Quebrada Honda / Padre Hurtado) y con las calles cerradas.

Preguntas obligatorias (responde cada una con: correcto / incorrecto / dudoso, justificación y corrección propuesta con ecuación):
1. ¿Es correcta la forma de S_f y la conversión de unidades (poise→Pa·s, dyn/cm²→Pa)? ¿Son correctos los coeficientes α, β citados para
   glenwood1/2 y aspen1/2? Si no los recuerdas con certeza, dilo y propone valores con fuente.
2. ¿Es consistente aplicar S_f como fricción por paso fraccionado sobre la cantidad de movimiento de ANUGA (que resuelve agua)?
   ¿Falta algún término por la densidad de la mezcla (gradiente de presión, término de fuente de pendiente)? ¿Cómo debe tratarse el umbral τy
   (condición de parada y reinicio) para no generar flujo espurio ni detener artificialmente el frente?
3. ¿El esquema semi-implícito es estable y consistente (orden, positividad)? ¿Hay sesgo con Δt variable de ANUGA?
4. Valores razonables de K (resistencia laminar) y n_td para calles urbanas pavimentadas y para cauce natural.
5. Representar edificios como bloques elevados vs. agujeros en la malla vs. porosidad: ¿qué es más defendible con DTM de 30 m?
6. ¿Qué incertidumbre introduce un DTM de 30 m remuestreado a 4 m en la dirección del flujo urbano? ¿Cómo cuantificarla
   (p. ej. Monte Carlo con ruido espacialmente correlacionado)?
7. Diseño de la calibración inversa (Qp, V, Cv): métrica recomendada (índice de Jaccard / F-score / distancia de alcance),
   identificabilidad de parámetros y cómo reportar el intervalo de confianza del recorrido.
8. Cualquier error físico/matemático adicional que veas.
Entrega un informe en markdown con secciones numeradas y una lista final "Cambios obligatorios" y "Cambios recomendados".