RESUMO A discretização da rede fluvial é uma etapa crítica para o cálculo da propagação de vazões em modelos hidrológicos. No entanto, quando se trata de modelagem hidrológico-hidrodinâmica mais complexa, adaptações na representação espacial das unidades de cálculo do modelo são necessárias para permitir simulações eficientes, especialmente para aplicações em grande escala. O objetivo deste artigo foi avaliar os impactos da discretização dos rios nas vazões simuladas, níveis de água e estabilidade numérica de um modelo hidrológico-hidrodinâmico baseado em divisões por sub-bacias, usando um método de segmentação de comprimento do rio (Δx) fixo. O estudo de caso foi a bacia do rio Purus, um afluente do rio Amazonas, que abrange uma área que conta desde regiões de respostas rápidas a montante até rios de várzea a jusante. Os resultados indicam que as descargas máximas e mínimas são menos afetadas pelo Δx adotado, enquanto os níveis de água são mais influenciados por essa seleção. Mostra-se que, para a propagação unidimensional inercial local usando um modelo explícito, o Δx e o alfa da condição CFL (Courant-Friedrichs-Lewy) devem ser cuidadosamente escolhidos para evitar erros de balanço de massa. Além disso, um limitador de fluxo baseado em número Froude simples é proposto e testado.
ABSTRACT The discretization of river networks is a critical step for computing flow routing in hydrological models. However, when it comes to more complex hydrologic-hydrodynamic models, adaptations in the spatial representation of model calculation units are further required to allow cost-effective simulations, especially for large scale applications. The objective of this paper is to assess the impacts of river discretization on simulated discharge, water levels and numerical stability of a catchment-based hydrologic-hydrodynamic model, using a fixed river length (Δx) segmentation method. The case study was the Purus river basin, a sub-basin of the Amazon, which covers an area that accounts for rapid response upstream reaches to downstream floodplain rivers. Results indicate that the maximum and minimum discharges are less affected by the adopted Δx (reach-length), whereas water levels are more influenced by this selection. It is showed that for the explicit local inertial one-dimensional routing, Δx and the α parameter of CFL (Courant-Friedrichs-Lewy) condition must be carefully chosen to avoid mass balance errors. Additionally, a simple Froude number-based flow limiter to avoid numerical issues is proposed and tested.