SWEpy
Solver GPU de alto rendimiento para las ecuaciones de aguas someras 2D
SWEpy es un solver en Python para las ecuaciones de aguas someras bidimensionales (SWE). Usa el método de volúmenes finitos sobre mallas triangulares no estructuradas y corre en la GPU mediante CuPy, lo que permite simular de forma eficiente fenómenos hidrodinámicos de gran escala, como roturas de presas, tsunamis e inundaciones.
El solver implementa un esquema de volúmenes finitos central-upwind de alto orden, bien balanceado y que preserva la positividad, apto para frentes de mojado y secado y topografía irregular. Para reducir la difusión numérica propia de los volúmenes finitos, usa una reconstrucción WENO de segundo orden y un integrador temporal Runge-Kutta de tercer orden que preserva la estabilidad fuerte (SSP-RK3). Su diseño modular permite combinar distintas condiciones iniciales, condiciones de borde y términos fuente.
SWEpy se validó con casos de referencia, entre ellos la isla cónica de Synolakis y el flujo de Bryson sobre un obstáculo gaussiano, y se aplicó a simulaciones de gran escala de la rotura de la presa de Malpasset (1959) y del tsunami de Maule de 2010 (Mw 8,8). Entrega resultados de alta resolución en hardware de consumo y sirve tanto para investigación como para pronóstico operacional.
SWEpy requiere la arquitectura CUDA de NVIDIA para el cálculo en paralelo. Para asegurar estabilidad numérica y reproducibilidad, recomendamos administrar el entorno de Python con Miniconda.
Para consultas científicas, propuestas de colaboración o soporte técnico, escribe a los autores.