Publications

@article{sMartinezA_2024_1,

title = {Mixing-phase model for shear-induced contractive/dilative effects in unsteady water-sediment mixture flows},

author = {S. Martínez-Aranda and J. Fernández-Pato and P. García-Navarro},

url = {https://www.sciencedirect.com/science/article/pii/S0309170824000976},

doi = {https://doi.org/10.1016/j.advwatres.2024.104710},

issn = {0309-1708},

year  = {2024},

date = {2024-01-01},

journal = {Advances in Water Resources},

volume = {188},

pages = {104710},

abstract = {Among the geophysical surface processes, mud and debris flows show one of the most complex and challenging behaviour for scientists and modellers. These flows consist of highly-unsteady gravity-driven movements of water-sediment mixtures with non-Newtonian rheology where the solid concentration could be about 40%–80% of the flow volume and which occur along steep and irregular terrains. Furthermore, the appearance of dynamic pressures in the fluid filling the intergranular pores increases the complexity and dominates the behaviour of the fluidized water-sediment material, leading to the appearance of significant density gradients during the movement. The dynamic pressure in the pore-fluid changes the effective normal stress within the mobilized material, affecting the frictional shear stress between grains and leading to the solid phase dilation/contraction. This must be properly accounted for when developing realistic models for water-sediment surface flows. In this work, a novel physically-based approach for modelling multi-grain dense-packed water-sediment flows is presented. A novel closure formulation for the pressure distribution within the pore-fluid during the movement of dense-packed water-sediment materials has been derived. This closure allows to relate the appearance of shear-induced dynamic pore pressures to the contractive/dilative behaviour of the solid aggregate. The resultant system of depth-averaged conservation laws includes continuity of the density-variable water-sediment material and the different solid classes transported in the flow, as well as the linear momentum equation for the fluidized bulk material, and it is solved using a well-balanced fully-coupled Finite Volume (FV) method. The resultant simulation tool is faced to synthetic, laboratory and real-scale benchmark cases to test its robustness and accuracy. The presence of dynamic pore pressures within the pore-fluid leads to the appearance of a deviatoric contribution to the solid flux, which causes the shear-induced separation of the solid and liquid phases and sustains the flow mobility for long distances, as it has been observed in real mud and debris events.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{diagnostics14111067,

title = {A New Quantitative Tool for the Ultrasonographic Assessment of Tendons: A Reliability and Validity Study on the Patellar Tendon},

author = {Isabel Albarova-Corral and José Segovia-Burillo and Miguel Malo-Urriés and Izarbe Ríos-Asín and Jesús Asín and Jorge Castillo-Mateo and Zeus Gracia-Tabuenca and Mario Morales-Hernández},

url = {https://www.mdpi.com/2075-4418/14/11/1067},

year  = {2024},

date = {2024-01-01},

journal = {Diagnostics},

volume = {14},

pages = {DOI: 10.3390/diagnostics14111067},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{valles2024_01,

title = {A 2D shallow water flow model with 1D internal boundary condition for subgrid-scale topography},

author = {P. Vallés and J. Fernández-Pato and M. Morales-Hernández and I. Echeverribar and P. García-Navarro},

url = {https://www.sciencedirect.com/science/article/pii/S0309170824001039},

doi = {https://doi.org/10.1016/j.advwatres.2024.104716},

year  = {2024},

date = {2024-01-01},

journal = {Advances in Water Resources},

volume = {189},

pages = {104716},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{ORTEGAMOYA2024113378,

title = {A vertically non-uniform temperature approach for the friction term computation in depth-averaged viscoplastic lava flows},

author = {J. Ortega-Moya and S. Martínez-Aranda and J. Fernández-Pato and P. García-Navarro},

url = {https://www.sciencedirect.com/science/article/pii/S0021999124006260},

doi = {https://doi.org/10.1016/j.jcp.2024.113378},

issn = {0021-9991},

year  = {2024},

date = {2024-01-01},

journal = {Journal of Computational Physics},

volume = {519},

pages = {113378},

abstract = {Recently, depth-averaged shallow flow models have been adapted to modelling liquefied lava flows, generally characterized by a marked temperature-dependent non-Newtonian rheology. Modelling these complex flows requires to include the effects of the depth-averaged temperature gradients in the governing equations. The most significant term to correctly predict the lava mobility is the flow resistance term, which is widely estimated using the linear viscoplastic Bingham model. This non-Newtonian model allows to relate the bed shear stress to the depth-averaged lava flow features by means of a cubic equation with analytical solution when assuming a uniform temperature distribution along the vertical profile. Nevertheless, the lava temperature is non-uniform along the vertical due to the heat transfer at the bottom and the free surface, and hence the classical cubic Bingham model is not valid anymore. In this work, a depth-averaged shallow flow model is adapted for realistic lava flows considering influence of the non-uniform vertical temperature profile in the non-Newtonian resistance. This requires to modify the rheological viscoplastic models for ensuring the coupling between flow dynamics and temperature evolution. Three non-uniform temperature vertical distributions are considered: linear, piece-wise and diffusion profiles. Synthetic tests are used to show the influence of the temperature vertical profile on the numerical results. Furthermore, laboratory experimental data are used to validate this novel viscoplastic resistance formulation and to show that the calibration of its parameters is possible.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{NAVASMONTILLA2024106349,

title = {Exploring the potential of TENO and WENO schemes for simulating under-resolved turbulent flows in the atmosphere using Euler equations},

author = {A. Navas-Montilla and J. Guallart and P. Solán-Fustero and P. García-Navarro},

url = {https://www.sciencedirect.com/science/article/pii/S0045793024001816},

doi = {https://doi.org/10.1016/j.compfluid.2024.106349},

issn = {0045-7930},

year  = {2024},

date = {2024-01-01},

journal = {Computers & Fluids},

volume = {280},

pages = {106349},

abstract = {This paper focuses on the design and analysis of very high-order finite volume methods for the computation of simplified meso- and micro-scale atmospheric flows. In a dry atmosphere, these flows can be represented by the Euler equations with a gravitational source term. Two different approaches are considered here. While one of the approaches is fully conservative for the total energy, the other is formulated in a non-conservative form. The main focus of the paper is to analyze the performance of such models in combination with the traditional WENO reconstruction and the novel TENO reconstruction by examining the spectral properties of these reconstruction methods. The overarching goal is to determine whether the combination of these models and numerical schemes can be used to build an implicit Large Eddy Simulation framework, shedding light on their potential advantages or limitations in representing under-resolved atmospheric processes in the meso- and micro-scales.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{mairal2023_01,

title = {The entropy fix in augmented Riemann solvers in presence of source terms: Application to the Shallow Water Equations},

author = {J. Mairal and J. Murillo and P. García-Navarro},

url = {https://www.sciencedirect.com/science/article/pii/S0045782523005352},

year  = {2023},

date = {2023-01-01},

journal = {Computer Methods in Applied Mechanics and Engineering},

volume = {417A},

pages = {DOI: 10.1016/j.cma.2023.116411},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{murillo2023_02,

title = {Numerical coupling of 0D and 1D models in networks of vessels including transonic flow conditions. 

Application to short-term transient and stationary hemodynamic simulation of postural changes},

author = {J. Murillo and P. García-Navarro},

url = {https://onlinelibrary.wiley.com/doi/full/10.1002/cnm.3751},

year  = {2023},

date = {2023-01-01},

journal = {Int J Numer Meth Biomedical Engineering},

pages = {doi:10.1002/cnm.3751},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{jfpato2023_01,

title = {A 2D Hydraulic Simulation Model Including Dynamic Piping and Overtopping Dambreach},

author = {J. Fernández-Pato and S. Martínez-Aranda and P. García-Navarro},

url = {https://www.mdpi.com/2073-4441/15/18/3268},

year  = {2023},

date = {2023-01-01},

journal = {Water},

volume = {15(18)},

pages = {https://doi.org/10.3390/w15183268},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{cav2023_01,

title = {SERGHEI (SERGHEI-SWE) v1.0: a performance-portable high-performance parallel-computing shallow-water solver for hydrology 

and environmental hydraulics},

author = {D. Caviedes-Voullieme and M. Morales-Hernández and M. R. Norman and I. Ozgen-Xian},

url = {https://gmd.copernicus.org/articles/16/977/2023/},

year  = {2023},

date = {2023-01-01},

journal = {Geosci. Model Dev.},

volume = {16},

pages = {977-1008},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{mar2023_01,

title = {Impacts of climate change on future hurricane induced rainfall and flooding in a coastal watershed: a case study on Hurricane Harvey.},

author = {Xudong Li and Dan Fu and John Nielsen-Gammon and Sudershan Gangrade and Shih-Chieh Kao and Ping Chang and Zhe Zhang Nathalie Voisin Mario Morales Hernández and Huilin Gao},

url = {https://www.sciencedirect.com/science/article/pii/S0022169422013440},

year  = {2023},

date = {2023-01-01},

journal = {Journal of Hydrology},

volume = {616},

pages = {128774},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{solan2023_01,

title = {Development of POD-based Reduced Order Models applied to shallow water equations using augmented Riemann solvers},

author = {P. Solán-Fustero and J. L. Gracia and A. Navas-Montilla and P. García-Navarro},

url = {https://www.sciencedirect.com/science/article/pii/S0045782523001627},

year  = {2023},

date = {2023-01-01},

journal = {Computer Methods in Applied Mechanics and Engineering},

volume = {410},

pages = {116038},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{martinez2023_01,

title = {Non-Equilibrium Bedload Transport Model Applied to Erosive Overtopping Dambreach},

author = {S. Martínez-Aranda and J. Fernández-Pato and P. García-Navarro},

url = {https://www.mdpi.com/2073-4441/15/17/3094},

year  = {2023},

date = {2023-01-01},

journal = {Water},

volume = {15(17)},

pages = {3094},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{valles2023_01,

title = {2D numerical simulation of floods in Ebro River and analysis of strategies to model the Mequinenza Reservoir},

author = {Pablo Vallés and Isabel Echeverribar and Juan Mairal and Sergio Martínez-Aranda and Javier Fernández-Pato and Pilar García-Navarro},

url = {https://www.mdpi.com/2624-795X/4/2/9},

year  = {2023},

date = {2023-01-01},

journal = {Geohazards},

volume = {1},

pages = {1-23},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{isa01,

title = {Extension of a Roe-type Riemann Solver Scheme to model Non-Hydrostatic Pressure Shallow Flows},

author = {I. Echeverribar and P. Brufau and P. García-Navarro},

year  = {2023},

date = {2023-01-01},

journal = {Applied Mathematics and Computation},

volume = {440},

pages = {127642},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{isa02,

title = {A GPU-based 2D viscous flow model with variable density and heat exchange},

author = {I. Echeverribar and S. Martínez-Aranda and J. Fernández-Pato and P. García-Navarro},

year  = {2023},

date = {2023-01-01},

journal = {Advances in Engineering Software},

volume = {175},

pages = {103340},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{GANGRADE2023130157,

title = {Unraveling the 2021 Central Tennessee flood event using a hierarchical multi-model inundation modeling framework},

author = {Sudershan Gangrade and Ganesh R. Ghimire and Shih-Chieh Kao and Mario Morales-Hernández and Ahmad A. Tavakoly and Joseph L. Gutenson and Kent H. Sparrow and George K. Darkwah and Alfred J. Kalyanapu and Michael L. Follum},

url = {https://www.sciencedirect.com/science/article/pii/S0022169423010995},

doi = {https://doi.org/10.1016/j.jhydrol.2023.130157},

issn = {0022-1694},

year  = {2023},

date = {2023-01-01},

journal = {Journal of Hydrology},

volume = {625},

pages = {130157},

abstract = {Flood prediction systems need hierarchical atmospheric, hydrologic, and hydraulic models to predict rainfall, runoff, streamflow, and floodplain inundation. The accuracy of such systems depends on the error propagation through the modeling chain, sensitivity to input data, and choice of models. In this study, we used multiple precipitation forcings (hindcast and forecast) to drive hydrologic and hydrodynamic models to analyze the impacts of various drivers on the estimates of flood inundation depth and extent. We implement this framework to unravel the August 2021 extreme flooding event that occurred in Central Tennessee, USA. We used two radar-based quantitative precipitation estimates (STAGE4 and MRMS) as well as quantitative precipitation forecasts (QPF) from the National Weather Service Weather Prediction Center (WPC) to drive a series of models in the hierarchical framework, including the Variable Infiltration Capacity (VIC) land surface model, the Routing Application for Parallel Computation of Discharge (RAPID) river routing model, and the AutoRoute and TRITON inundation models. An evaluation with observed high-water marks demonstrates that the framework can reasonably simulate flood inundation. Despite the complex error propagation mechanism of the modeling chain, we show that inundation estimates are most sensitive to rainfall estimates. Most notably, QPF significantly underestimates flood magnitudes and inundations leading to unanticipated severe flooding for all stakeholders involved in the event. Finally, we discuss the implications of the hydrodynamic modeling framework for real-time flood forecasting.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{adrian01,

title = {Numerical characterization of seiche waves energy potential in river bank lateral embayments},

author = {C. Juez and A. Navas-Montilla},

year  = {2022},

date = {2022-01-01},

journal = {Renewable Energy},

volume = {186},

pages = {143-156},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{adrian03,

title = {A POD-based ROM strategy for the prediction in time of advection-dominated problems},

author = {P. Solán-Fustero and J. L. Gracia and A. Navas-Montilla and P. Garc'ia-Navarro},

year  = {2022},

date = {2022-01-01},

journal = {Journal of Computational Physics},

volume = {471},

pages = {111672},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{sergioLaserProfile,

title = {Laser Profilometry to measure the bed evolution in a dam-break flow},

author = {R. Meurice and S. Martínez-Aranda and M. Ebrahimi and P. García-Navarro and S. Soares-Frazão},

doi = {https://doi.org/10.1080/00221686.2022.2059584},

year  = {2022},

date = {2022-01-01},

journal = {Journal of Hydraulic Research},

volume = {(60)},

number = {5},

pages = {725-737},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{sergio2Dresistance,

title = {Novel discretization strategies for the 2D non-Newtonian resistance term in geophysical shallow flows},

author = {S. Martínez-Aranda and J. Murillo and M. Morales-Hernández and P García-Navarro},

doi = {https://doi.org/10.1016/j.enggeo.2022.106625},

year  = {2022},

date = {2022-01-01},

journal = {Engineering Geology},

volume = {302},

pages = {106625},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{sergiomudflowGPU,

title = {A GPU-accelerated Efficient Simulation Tool (EST) for 2D variable-density mud/debris flows over non-uniform erodible beds},

author = {S. Martínez-Aranda and J. Murillo and P García-Navarro},

doi = {https://doi.org/10.1016/j.enggeo.2021.106462},

year  = {2022},

date = {2022-01-01},

journal = {Engineering Geology},

volume = {296},

pages = {106462},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{psolanfus01,

title = {A POD-based ROM strategy for the prediction in time of advection-dominated problems},

author = {P. Solán-Fustero and J. L. Gracia and A. Navas-Montilla and P García-Navarro},

doi = {https://doi.org/10.1016/j.jcp.2022.111672},

year  = {2022},

date = {2022-01-01},

journal = {Journal of Computational Physics},

volume = {471},

pages = {111672},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{sym13091658,

title = {A Solution of the Junction Riemann Problem for 1D Hyperbolic Balance Laws in Networks including Supersonic Flow Conditions on Elastic Collapsible Tubes},

author = {Javier Murillo and Pilar García-Navarro},

url = {https://www.mdpi.com/2073-8994/13/9/1658},

doi = {10.3390/sym13091658},

issn = {2073-8994},

year  = {2021},

date = {2021-01-01},

journal = {Symmetry},

volume = {13},

number = {9},

abstract = {The numerical modeling of one-dimensional (1D) domains joined by symmetric or asymmetric bifurcations or arbitrary junctions is still a challenge in the context of hyperbolic balance laws with application to flow in pipes, open channels or blood vessels, among others. The formulation of the Junction Riemann Problem (JRP) under subsonic conditions in 1D flow is clearly defined and solved by current methods, but they fail when sonic or supersonic conditions appear. Formulations coupling the 1D model for the vessels or pipes with other container-like formulations for junctions have been presented, requiring extra information such as assumed bulk mechanical properties and geometrical properties or the extension to more dimensions. To the best of our knowledge, in this work, the JRP is solved for the first time allowing solutions for all types of transitions and for any number of vessels, without requiring the definition of any extra information. The resulting JRP solver is theoretically well-founded, robust and simple, and returns the evolving state for the conserved variables in all vessels, allowing the use of any numerical method in the resolution of the inner cells used for the space-discretization of the vessels. The methodology of the proposed solver is presented in detail. The JRP solver is directly applicable if energy losses at the junctions are defined. Straightforward extension to other 1D hyperbolic flows can be performed.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{isa03,

title = {Efficient Reservoir Modelling for Flood Regulation in the Ebro River (Spain)},

author = {P. Echeverribar and P. Vallés and J. Mairal and P. García-Navarro},

year  = {2021},

date = {2021-01-01},

journal = {Water},

volume = {12},

pages = {3160},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{mario01,

title = {TRITON: A Multi-GPU open source 2D hydrodynamic flood model},

author = {M. Morales-Hernández and Md B. Sharif and A. Kalyanapu and S. K. Ghafoor and T.T.Dullo and S. Gangrade and S. -C. Kao and M. R. Norman and K. J. Evans},

doi = {https://doi.org/10.1016/j.envsoft.2021.105034},

year  = {2021},

date = {2021-01-01},

journal = {Environmental Modelling and Software},

volume = {141},

pages = {105034},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{mario02,

title = {Evaluating precipitation, streamflow, and inundation forecasting skills during extreme weather events: A case study for an urban watershed},

author = {Xudong Li and Cheryl Rankin and Sudershan Gangrade and Gang Zhao and Kris Lander and Nathalie Voisin and Manqing Shao and Mario Morales-Hernández and Shih-Chieh Kao and Huilin Gao},

doi = {https://doi.org/10.1016/j.jhydrol.2021.127126},

year  = {2021},

date = {2021-01-01},

journal = {Journal of Hydrology},

volume = {603(Part D)},

pages = {127126},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{mario03,

title = {Assessing climate-change-induced flood risk in the Conasauga River watershed: an application of ensemble hydrodynamic inundation modeling},

author = {T. T. Dullo and G. K. Darkwah and S. Gangrade and M. Morales-Hernández and M. B. Sharif and A. J. Kalyanapu and S. -C. Kao and S. Ghafoor and M. Ashfaq},

doi = {https://doi.org/10.5194/nhess-21-1739-2021},

year  = {2021},

date = {2021-01-01},

journal = {Nat. Hazards Earth Syst.},

volume = {21},

pages = {1739-1757},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{mario04,

title = {Transport Control at Channel Junctions Using Adjoint Sensitivity},

author = {G. Gordillo and M. Morales-Hernández and P. García-Navarro},

doi = {https://doi.org/10.3390/math10010093},

year  = {2021},

date = {2021-01-01},

journal = {Mathematics},

volume = {10},

pages = {93},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{FernandezPato2021,

title = {An efficient GPU implementation of a coupled overland-sewer hydraulic model with pollutant transport},

author = {J. Fernández-Pato and P. García-Navarro},

doi = {doi.org/10.3390/hydrology8040146},

year  = {2021},

date = {2021-01-01},

journal = {Hydrology},

volume = {8(4)},

pages = {146},

abstract = {Numerical simulation of flows that consider interaction between overland and drainage networks has become a practical tool to prevent and mitigate flood situations in urban environments, especially when dealing 

with intense storm events, where the limited capacity of the sewer systems can be a trigger for flooding. Additionally, in order to prevent any kind of pollutant dispersion through the drainage network, it is very interesting 

to have a certain monitorization or control over the quality of the water that flows in both domains. In this sense, the addition of a pollutant transport component to both surface and sewer hydraulic models would benefit the 

global analysis of the combined water flow. On the other hand, when considering a realistic large domain with complex topography or streets structure, a fine spatial discretization is mandatory. Hence the number of grid cells 

is usually very large and, therefore, it is necessary to use parallelization techniques for the calculation, the use of Graphic Processing Units (GPU) being one of the most efficient due to the leveraging of thousands of processors 

within a single device. In this work, an efficient GPU-based 2D shallow water flow solver (RiverFlow2D-GPU) is fully coupled with EPA�s Storm Water Management Model (SWMM). Both models are able to develop a transient water quality 

analysis taking into account several pollutants. The coupled model, referred to as RiverFlow2D-GPU UD (Urban Drainge) is applied to three real-world cases, covering the most common hydraulic situations in urban hydrology/hydraulics. 

A UK Environmental Agency test case is used as model validation, showing a good agreement between RiverFlow2D-GPU UD and the rest of the numerical models considered. The efficiency of the model is proven in two more complex domains, 

leading to a >100x faster simulations compared with the traditional CPU computation.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{morales2021c,

title = {TRITON: A Multi-GPU open source 2D hydrodynamic flood model},

author = {M. Morales-Hernández and Md B. Sharif and A. Kalyanapu and S. K. Ghafoor and T. T. Dullo and S. Gangrade and S. -C. Kao and M. R. Norman and K. J. Evans},

doi = {10.1016/j.envsoft.2021.105034},

issn = {1364-8152},

year  = {2021},

date = {2021-01-01},

journal = {Environmental Modelling & Software},

volume = {141},

pages = {105034},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{spectral_analysis,

title = {Application of approximate dispersion-diffusion analyses to under-resolved Burgers turbulence using high resolution WENO and UWC schemes},

author = {P. Solán-Fustero and A. Navas-Montilla and E. Ferrer and J. Manzanero and P. García-Navarro},

doi = {10.1016/j.jcp.2021.110246},

year  = {2021},

date = {2021-01-01},

journal = {Journal of Computational Physics},

volume = {435},

pages = {110246},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{morale2021b,

title = {Assessing climate-change-induced flood risk in the Conasauga River watershed: an application of ensemble hydrodynamic inundation modeling},

author = {Tigstu T. Dullo and George K. Darkwah and Sudershan Gangrade and Mario Morales-Hernández and M. Bulbul Sharif and Alfred J. Kalyanapu and Shih-Chieh Kao and Sheikh Ghafoor and Moetasim Ashfaq},

doi = {10.5194/nhess-21-1739-2021},

year  = {2021},

date = {2021-01-01},

volume = {21},

number = {6},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{bedload_nonCapacity,

title = {Comparative analysis of HLLC- and Roe-based models for the simulation of a dambreak flow in an erodible channel with a 90 bend},

author = {S. Martínez-Aranda and R. Meurice and S. Soares-Frazão and P. García-Navarro},

doi = {10.3390/w13131840},

year  = {2021},

date = {2021-01-01},

journal = {Water},

volume = {13},

pages = {1840},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{navasmontilla2021a,

title = {Hyperbolic Reformulation Approach to Enable Efficient Simulation of Groundwater Flow and Reactive Transport},

author = {Ilhan Özgen-Xian and Adrian Navas-Montilla and Dipankar Dwivedi and Sergi Molins},

doi = {10.1089/ees.2020.0363},

year  = {2021},

date = {2021-01-01},

journal = {Environmental Engineering Science},

volume = {38},

number = {3},

pages = {181-191},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{bedload_fully_coupled,

title = {Comparison of new efficient 2D models for the simulation of bedload transport using the augmented Roe approach},

author = {S. Martínez-Aranda and J. Murillo and P. García-Navarro},

doi = {10.1016/j.advwatres.2021.103931},

year  = {2021},

date = {2021-01-01},

journal = {Advances in Water Resources},

volume = {153},

pages = {103931},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{morales2021a,

title = {Simulation of Hurricane Harvey flood event through coupled hydrologic-hydraulic models: Challenges and next steps},

author = {Tigstu T. Dullo and Sudershan Gangrade and Mario Morales-Hernández and Md Bulbul Sharif and Shih-Chieh Kao and Alfred J. Kalyanapu and Sheikh Ghafoor and Katherine J. Evans},

doi = {10.1111/jfr3.12716},

year  = {2021},

date = {2021-01-01},

journal = {Journal of Flood Risk Management},

volume = {n/a},

pages = {e12716},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{kinect_cavity,

title = {2D experiments and numerical simulation of the oscillatory shallow flow in an open channel lateral cavity},

author = {A. Navas-Montilla and S. Martínez-Aranda and A. Lozano and I. García-Palacín and P. García-Navarro},

doi = {10.1016/j.advwatres.2020.103836},

year  = {2021},

date = {2021-01-01},

journal = {Advances in Water Resources},

volume = {148},

pages = {103836},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{10.2166/hydro.2020.163,

title = {High-performance computing in water resources hydrodynamics},

author = {M. Morales-Hernández and M. B. Sharif and S. Gangrade and T. T. Dullo and S. -C. Kao and A. Kalyanapu and S. K. Ghafoor and K. J. Evans and E. Madadi-Kandjani and B. R. Hodges},

doi = {10.2166/hydro.2020.163},

year  = {2020},

date = {2020-01-01},

journal = {Journal of Hydroinformatics},

volume = {22},

number = {5},

pages = {1217-1235},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{waterquality,

title = {A GPU-based 2D Shallow Water quality model},

author = {G. Gordillo and M. Morales-Hernández and I. Echeverribar and J. Fernández-Pato and P. García-Navarro},

year  = {2020},

date = {2020-01-01},

journal = {Journal of Hydroinformatics},

volume = {22(5)},

pages = {1182-1197},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{FERNANDEZPATO2020b,

title = {Analysis of the performance of different culvert boundary conditions in 2D shallow flow models},

author = {J. Fernández-Pato and S. Martínez-Aranda and M. Morales-Hernández and P. García-Navarro},

doi = {10.2166/hydro.2020.025},

year  = {2020},

date = {2020-01-01},

journal = {Journal of Hydroinformatics},

volume = {22(5)},

pages = {1093-1121},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{geovanny01,

title = {A gradient-descent adjoint method for the reconstruction of boundary conditions in a river flow nitrification model},

author = {G. Gordillo and M. Morales-Hernández and P. García-Navarro},

year  = {2020},

date = {2020-01-01},

journal = {Environmental Science: Processes & Impacts},

volume = {22},

pages = {381-397},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{navas2020a,

title = {Discontinuous Galerkin well-balanced schemes using augmented Riemann solvers with application to the shallow water equations},

author = {A. Navas-Montilla and P. Sol'n-Fustero and J. Murillo and P. García-Navarro},

year  = {2020},

date = {2020-01-01},

journal = {Journal of Hydroinformatics},

volume = {22(5)},

pages = {1038-1058},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{murillo2020c,

title = {Adaptation of flux-based solvers to 2D two-layer shallow flows with variable density including numerical treatment of the loss of hyperbolicity and drying/wetting fronts},

author = {J. Murillo and S. Martínez-Aranda and A. Navas-Montilla and P. García-Navarro},

doi = {10.2166/hydro.2020.207},

year  = {2020},

date = {2020-01-01},

journal = {Journal of Hydroinformatics},

volume = {22(5)},

pages = {972-1014},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{martinez2020b,

title = {A 1D shallow-flow model for two-layer flows based on FORCE scheme with wet-dry treatment},

author = {S. Martínez-Aranda and A. Ramos-Pérez and P. García-Navarro},

doi = {10.2166/hydro.2020.002},

year  = {2020},

date = {2020-01-01},

journal = {Journal of Hydroinformatics},

volume = {22(5)},

pages = {1015-1037},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{martinez2020a,

title = {A robust two-dimensional model for highly sediment-laden unsteady flows of variable density over movable beds},

author = {S. Martínez-Aranda and J. Murillo and P. García-Navarro},

doi = {10.2166/hydro.2020.027},

year  = {2020},

date = {2020-01-01},

journal = {Journal of Hydroinformatics},

volume = {22(5)},

pages = {1038-1060},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{echeverribar2020a,

title = {Analysis of the performance of a hybrid CPU/GPU 1D2D coupled model for real flood cases},

author = {I. Echeverribar and M. Morales-Hernández and P. Brufau and P. García-Navarro},

year  = {2020},

date = {2020-01-01},

journal = {Journal of Hydroinformatics},

volume = {22(5)},

pages = {1198-1216},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{FERNANDEZPATO2020103617,

title = {A 2D finite volume simulation tool to enable the assessment of combined hydrological and morphodynamical processes in mountain catchments},

author = {J. Fernández-Pato and S. Martínez-Aranda and P. García-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S0309170819308334},

doi = {https://doi.org/10.1016/j.advwatres.2020.103617},

issn = {0309-1708},

year  = {2020},

date = {2020-01-01},

journal = {Advances in Water Resources},

volume = {141},

pages = {103617},

abstract = {Nowadays, the great power of modern computers allows to develop computational models able to deal with simulations of several coupled phenomena over detailed complex topography. An efficient and properly calibrated computational model represents a useful tool to provide insight into the catchment dynamics at hydrological and geomorphological levels. In addition, it allows to develop detailed risk management and conservation plans. In this work, we present a coupled surface-groundwater distributed flow model with hydrological (rainfall and infiltration) and geomorphological (suspended and bed load sediment transport) components. The coupled model is applied to well characterized experimental catchments that are used as realistic test cases. The calibration of the water flow model response to rainfall is performed by means of the fitting to experimental outlet hydrographs of the results supplied by a coupled formulation of 2D Shallow Water Equations and 2D Darcy’s law for saturated porous media connected via suitable infiltration laws. The calibration of a suspended and bed load model is also addressed by means of the fitting to experimental outlet sedigraphs. The numerical results show a good agreement between numerical and observed hydrographs and sedigraphs, significantly improving previous published simulations. Additionally, the need to repeat the simulations in the calibration processes is no longer an unapproachable problem.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{CAVIEDESVOULLIEME2020124663,

title = {Performance assessment of 2D Zero-Inertia and Shallow Water models for simulating rainfall-runoff processes},

author = {Daniel Caviedes-Voullième and Javier Fernández-Pato and Christoph Hinz},

url = {http://www.sciencedirect.com/science/article/pii/S0022169420301232},

doi = {https://doi.org/10.1016/j.jhydrol.2020.124663},

issn = {0022-1694},

year  = {2020},

date = {2020-01-01},

journal = {Journal of Hydrology},

volume = {584},

pages = {124663},

abstract = {Rainfall-runoff simulations are increasingly being performed with physically-based and spatially distributed solvers. The current computational and numerical technology enables the use of full shallow water equations solvers to be applied for these type of flow problems. Nonetheless, Zero-Inertia (diffusive wave) solvers have been historically favoured due to their conceptual and mathematical simplicity in comparison to shallow water solvers, with the working assumption that the simplifications introduced by Zero-Inertia will have some assumable impact on accuracy but will also allow for computational efficiency. Since both types of solvers have been primarily developed, benchmarked and compared to each other for fluvial and floodplain simulations, it is relevant to assess t-he relative performance for rainfall-runoff problems. In this work, both solvers are applied to a set of six well known test cases with reference solutions. The performance of the solvers is assessed in terms of global signatures such as hydrographs and flooded areas, but also in terms of spatial distributions of depth and velocity, as well as computational cost. Furthermore, the comparisons are performed across different spatial resolutions. The results show that for rainfall-runoff problems explicit, finite volumes solvers for both equations provide a similar accuracy, but the shallow water solver requires less computational time. The Zero-Inertia solver was found to be less sensitive to mesh refining than the full shallow water solver.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{navas2020b,

title = {Computational hemodynamics in arteries with the one-dimensional augmented fluid-structure interaction system: viscoelastic parameters estimation and comparison with in-vivo data},

author = {G. Bertaglia and A. Navas-Montilla and A. Valiani and M. I. M. García and J. Murillo and V. Caleffi},

year  = {2020},

date = {2020-01-01},

journal = {Journal of Biomechanics},

volume = {100},

pages = {109595},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{navas2020c,

title = {Numerical ability of hyperbolic flux solvers to compute 2D shear layers in turbulent shallow flows},

author = {A. Navas-Montilla and C. Juez},

year  = {2020},

date = {2020-01-01},

journal = {Advances in Water Resources},

volume = {135},

pages = {103482},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{10.2166/hydro.2019.080,

title = {Finite volume model for the simulation of 1D unsteady river flow and water quality based on the WASP},

author = {Geovanny Gordillo and Mario Morales-Hernández and Pilar García-Navarro},

url = {https://doi.org/10.2166/hydro.2019.080},

doi = {10.2166/hydro.2019.080},

year  = {2020},

date = {2020-01-01},

journal = {Journal of Hydroinformatics},

volume = {22(2)},

pages = {327-345},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{García-Navarro2019,

title = {The shallow water equations and their application to realistic cases},

author = {P. García-Navarro and J. Murillo and J. Fernández-Pato and I. Echeverribar and M. Morales-Hernández},

url = {https://doi.org/10.1007/s10652-018-09657-7},

doi = {10.1007/s10652-018-09657-7},

issn = {1573-1510},

year  = {2019},

date = {2019-01-05},

journal = {Environmental Fluid Mechanics},

abstract = {The numerical modelling of 2D shallow flows in complex geometries involving transient flow and movable boundaries has been a challenge for researchers in recent years. There is a wide range of physical situations of environmental interest, such as flow in open channels and rivers, tsunami and flood modelling, that can be mathematically represented by first-order non-linear systems of partial differential equations, whose derivation involves an assumption of the shallow water type. Shallow water models may include more sophisticated terms when applied to cases of not pure water floods, such as mud/debris floods, produced by landslides. Mud/debris floods are unsteady flow phenomena in which the flow changes rapidly, and the properties of the moving fluid mixture include stop and go mechanisms. The present work reports on a numerical model able to solve the 2D shallow water equations even including bed load transport over erodible bed in realistic situations involving transient flow and movable flow boundaries. The novelty is that it offers accurate and stable results in realistic problems since an appropriate discretization of the governing equations is performed. Furthermore, the present work is focused on the importance of the computational cost. Usually, the main drawback is the high computational effort required for obtaining accurate numerical solutions due to the high number of cells involved in realistic cases. However, the proposed model is able to reduce computer times by orders of magnitude making 2D applications competitive and practical for operational flood prediction. Moreover our results show that high performance code development can take advantage of general purpose and inexpensive Graphical Processing Units, allowing to run almost 100 times faster than old generation codes in some cases.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{MARTINEZARANDA201991,

title = {A comparative analysis of capacity and non-capacity formulations for the simulation of unsteady flows over finite-depth erodible beds},

author = {S. Martínez-Aranda and J. Murillo and P. García-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S0309170819302453},

doi = {https://doi.org/10.1016/j.advwatres.2019.06.001},

issn = {0309-1708},

year  = {2019},

date = {2019-01-01},

journal = {Advances in Water Resources},

volume = {130},

pages = {91 - 112},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{FernandezPato2019b,

title = {Simulación de avenidas mediante un modelo hidráulico/hidrológico distribuido en un tramo urbano del río Ginel (Fuentes de Ebro)},

author = {J. Fernández-Pato and A. Sánchez and P. García-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S0045793019301185},

doi = {10.1080/23863781.2019.1622473},

year  = {2019},

date = {2019-01-01},

journal = {Ribagua},

volume = {6},

number = {1},

pages = {49-62},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{10.4995/ia.2019.12169,

title = {Modelos de simulación de alto orden para la resolución de fenómenos de propagación de ondas en flujos de lámina libre con turbulencia},

author = {A. Navas-Montilla and J. Murillo and P. García-Navarro},

url = {https://doi.org/10.4995/ia.2019.12169},

doi = {10.4995/ia.2019.12169},

issn = {1134-2196},

year  = {2019},

date = {2019-01-01},

journal = {Ingeniería del agua},

volume = {23},

number = {4},

pages = {275-287},

abstract = {En este trabajo se presenta una herramienta de simulación basada en la resolución transitoria de grandes remolinos (URANS o DA-LES) para flujos turbulentos de aguas poco profundas en los que la turbulencia es predominantemente horizontal. El aspecto fundamental del modelo es la combinación de una discretización de alto orden en espacio y tiempo con una modelización de los efectos en el flujo promedio de las escalas turbulentas no resueltas. El modelo propuesto garantiza con precisión de máquina el equilibrio hidrostático (propiedad well-balanced) gracias a la utilización de una formulación del flujo numérico que incluye los términos fuente en la resolución del problema de Riemann derivativo en las paredes de las celdas. Se presenta una validación del modelo utilizando datos de literatura para un experimento de laboratorio que involucra un flujo de aguas poco profundas sobre una isla cónica, que da lugar a la generación de una calle de vórtices aguas abajo de la isla. Los resultados numéricos muestran que el modelo propuesto es capaz de reproducir fenómenos turbulentos bidimensionales, proporcionando un mayor nivel de detalle que la aproximación RANS tradicional.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{BERTAGLIA2019109595,

title = {Computational hemodynamics in arteries with the one-dimensional augmented fluid-structure interaction system: viscoelastic parameters estimation and comparison with in-vivo data},

author = {Giulia Bertaglia and Adrián Navas-Montilla and Alessandro Valiani and Manuel Ignacio Monge García and Javier Murillo and Valerio Caleffi},

url = {http://www.sciencedirect.com/science/article/pii/S0021929019308589},

doi = {https://doi.org/10.1016/j.jbiomech.2019.109595},

issn = {0021-9290},

year  = {2019},

date = {2019-01-01},

journal = {Journal of Biomechanics},

pages = {109595},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{MURILLO201974,

title = {Formulation of exactly balanced solvers for blood flow in elastic vessels and their application to collapsed states},

author = {J. Murillo and A. Navas-Montilla and P. García-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S0045793019301185},

doi = {https://doi.org/10.1016/j.compfluid.2019.04.008},

issn = {0045-7930},

year  = {2019},

date = {2019-01-01},

journal = {Computers & Fluids},

volume = {186},

pages = {74 - 98},

abstract = {In this work, numerical solvers based on extensions of the Roe and HLL schemes are adapted to deal with test cases involving extreme collapsing conditions in elastic vessels. To achieve this goal, the system is transformed to provide a conservation–law form, allowing to define Rankine–Hugoniot conditions. The approximate solvers allow to describe the inner states of the solution. Therefore, source term fixes can be used to prevent unphysical values of vessel area and, at the same time, the eigenvalues of the system control stability. Numerical solvers of different order are tested using a wide variety of Riemann problems, including extreme vessel collapse and blockage. In all cases, the robustness of the approximate solvers presented here is checked using first and third order methods in time and space, using the WENO reconstruction scheme in combination with the TVDRK3 method.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{ECHEVERRIBAR2019103444,

title = {2D numerical simulation of unsteady flows for large scale floods prediction in real time},

author = {I. Echeverribar and M. Morales-Hernández and P. Brufau and P. García-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S0309170819304786},

doi = {https://doi.org/10.1016/j.advwatres.2019.103444},

issn = {0309-1708},

year  = {2019},

date = {2019-01-01},

journal = {Advances in Water Resources},

volume = {134},

pages = {103444},

abstract = {The challenge of finding a compromise between computational time and level of accuracy and robustness has traditionally expanded the use simplified models rather than full two-dimensional (2D) models for flood simulation. This work presents a GPU accelerated 2D shallow water model for the simulation of flood events in real time. In particular, an explicit first-order finite volume scheme is detailed to control the numerical instabilities that are likely to appear when used in complex topography. The model is first validated with the benchmark test case of the Toce River (Italy) and numerical fixes are demonstrated to be necessary. The model is next applied to reproduce real events in a reach of the Ebro River (Spain) in order to compare simulation results with field data. The second case deals with a large domain (744 km2) and long flood duration (up to 20 days) allowing an analysis of the performance and speed-up achieved by different GPU devices. The high values of fit between observed and simulated results as well as the computational times achieved are encouraging to propose the use of the model as forecasting system.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{MARTINEZARANDA2019,

title = {A 1D numerical model for the simulation of unsteady and highly erosive flows in rivers},

author = {S. Martínez-Aranda and J. Murillo and P. García-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S004579301930009X},

doi = {https://doi.org/10.1016/j.compfluid.2019.01.011},

issn = {0045-7930},

year  = {2019},

date = {2019-01-01},

journal = {Computers & Fluids},

abstract = {This work is focused on a numerical finite volume scheme for the coupled shallow water-Exner system in 1D applications with arbitrary geometry. The mathematical expressions modeling the hydrodynamic and morphodynamic components of the physical phenomenon are treated to deal with cross-section shape variations and empirical solid discharge estimations. The resulting coupled equations can be rewritten as a non-conservative hyperbolic system with three moving waves and one stationary wave to account for the source terms discretization. Moreover, the wave celerities for the coupled morpho-hydrodyamical system depend on the erosion-deposition mechanism selected to update the channel cross-section profile. This influence is incorporated into the system solution by means of a new parameter related to the channel bottom variation celerity. Special interest is put to show that, even for the simplest solid transport models as the Grass law, to find a linearized Jacobian matrix of the system can be a challenge in presence of arbitrary shape channels. In this paper a numerical finite volume scheme is proposed, based on an augmented Roe solver, first order accurate in time and space, dealing with solid transport flux variations caused by the channel geometry changes. Channel cross-section variations lead to the appearance of a new solid flux source term which should be discretized properly. The stability region is controlled by wave celerities together with a proper reconstruction of the approximate local Riemann problem solution, enforcing positive values for the intermediate states of the conserved variables. Comparison of the numerical results for several analytical and experimental cases demonstrates the effectiveness, exact well-balancedness and accuracy of the scheme.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{NAVASMONTILLA2019445,

title = {Improved Riemann solvers for an accurate resolution of 1D and 2D shock profiles with application to hydraulic jumps},

author = {A. Navas-Montilla and J. Murillo},

url = {http://www.sciencedirect.com/science/article/pii/S0021999118307496},

doi = {https://doi.org/10.1016/j.jcp.2018.11.023},

issn = {0021-9991},

year  = {2019},

date = {2019-01-01},

journal = {Journal of Computational Physics},

volume = {378},

pages = {445 - 476},

abstract = {From the early stages of CFD, the computation of shocks using Finite Volume methods has been a very challenging task as they often prompt the generation of numerical anomalies. Such anomalies lead to an incorrect and unstable representation of the discrete shock profile that may eventually ruin the whole solution. The two most widespread anomalies are the slowly-moving shock anomaly and the carbuncle, which are deeply addressed in the literature in the framework of homogeneous problems, such as Euler equations. In this work, the presence of the aforementioned anomalies is studied in the framework of the 1D and 2D SWE and novel solvers that effectively reduce both anomalies, even in cases where source terms dominate the solution, are presented. Such solvers are based on the augmented Roe (ARoe) family of Riemann solvers, which account for the source term as an extra wave in the eigenstructure of the system. The novel method proposed here is based on the ARoe solver in combination with: (a) an improved flux extrapolation method based on a previous work, which circumvents the slowly-moving shock anomaly and (b) a contact wave smearing technique that avoids the carbuncle. The resulting method is able to eliminate the slowly-moving shock anomaly for 1D steady cases with source term. When dealing with 2D cases, the novel method proves to handle complex shock structures composed of hydraulic jumps over irregular bathymetries, avoiding the presence of the aforementioned anomalies.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{García-Navarro2019b,

title = {Use of internal boundary conditions for levees representation: application to river flood management},

author = {I. Echeverribar and M. Morales-Hernández and P. Brufau and P. García-Navarro},

doi = {10.1007/s10652-018-09658-6},

year  = {2019},

date = {2019-01-01},

journal = {Environmental Fluid Mechanics},

abstract = {River floods can be simulated with the 2D shallow water system of equations using finite volume methods, where the terrain is discretized in cells that form the computational mesh. Usually a proper treatment of wet/dry fronts is required. River levees can be modelled as part of the topography by means of sufficiently small cells of higher elevation than the rest of the bed level in locally refined meshes. This procedure is associated with a large computational time since the time step depends directly on the cell size. The alternative proposed in this work includes the levees as internal boundary conditions in the 2D numerical scheme. In particular, levees have been defined by a weir law that, depending on the relative values of water surface levels on both sides, can formulate the discharge for different situations (i.e. free flow and submerged flow). In addition, having identified numerical difficulties in cases of low discharge under free flow conditions, a novel procedure to avoid oscillations has been developed and called volume transport method. The validation and comparison between methods has been carried out with benchmark test cases and, in addition, with a real flood event in the Ebro River (Spain)},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{García-Navarro2019,

title = {Use of internal boundary conditions for levees representation: application to river flood management},

author = {I. Echeverribar and M. Morales-Hernández and P. Brufau and P. García-Navarro},

url = {https://doi.org/10.1007/s10652-018-09657-7},

doi = {10.1007/s10652-018-09658-6},

year  = {2019},

date = {2019-01-01},

journal = {Environmental Fluid Mechanics},

abstract = {River floods can be simulated with the 2D shallow water system of equations using finite volume methods, where the terrain is discretized in cells that form the computational mesh. Usually a proper treatment of wet/dry fronts is required. River levees can be modelled as part of the topography by means of sufficiently small cells of higher elevation than the rest of the bed level in locally refined meshes. This procedure is associated with a large computational time since the time step depends directly on the cell size. The alternative proposed in this work includes the levees as internal boundary conditions in the 2D numerical scheme. In particular, levees have been defined by a weir law that, depending on the relative values of water surface levels on both sides, can formulate the discharge for different situations (i.e. free flow and submerged flow). In addition, having identified numerical difficulties in cases of low discharge under free flow conditions, a novel procedure to avoid oscillations has been developed and called volume transport method. The validation and comparison between methods has been carried out with benchmark test cases and, in addition, with a real flood event in the Ebro River (Spain)},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{Morales-Hernández2018,

title = {Diffusion–dispersion numerical discretization for solute transport in 2D transient shallow flows},

author = {M. Morales-Hernández and J. Murillo and P. García-Navarro},

url = {https://doi.org/10.1007/s10652-018-9644-2},

doi = {10.1007/s10652-018-9644-2},

issn = {1573-1510},

year  = {2018},

date = {2018-11-01},

journal = {Environmental Fluid Mechanics},

abstract = {The 2D solute transport equation can be incorporated into the 2D shallow water equations in order to solve both flow and solute interactions in a coupled system of equations. In order to solve this system, an explicit finite volume scheme based on Roe's linearization is proposed. Moreover, it is feasible to decouple the solute transport equation from the hydrodynamic system in a conservative way. In this case, the advection part is solved in essence defining a numerical flux, allowing the use of higher order numerical schemes. However, the discretization of the diffusion–dispersion terms have to be carefully analysed. In particular, time-step restrictions linked to the nature of the solute equation itself as well as the numerical diffusion associated to the numerical scheme used are question of interest in this work. These improvements are tested in an analytical case as well as in a laboratory test case with a passive solute (fluorescein) released from a reservoir. Experimental measurements are compared against the numerical results obtained with the proposed model and a sensitivity analysis is carried out, confirming an agreement with the longitudinal coefficients and an underestimation of the transversal ones, respectively.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{Morales-Hernández2018,

title = {Diffusion–dispersion numerical discretization for solute transport in 2D transient shallow flows},

author = {M. Morales-Hernández and J. Murillo and P. García-Navarro},

url = {https://doi.org/10.1007/s10652-018-9644-2},

doi = {10.1007/s10652-018-9644-2},

issn = {1573-1510},

year  = {2018},

date = {2018-11-01},

journal = {Environmental Fluid Mechanics},

abstract = {The 2D solute transport equation can be incorporated into the 2D shallow water equations in order to solve both flow and solute interactions in a coupled system of equations. In order to solve this system, an explicit finite volume scheme based on Roe's linearization is proposed. Moreover, it is feasible to decouple the solute transport equation from the hydrodynamic system in a conservative way. In this case, the advection part is solved in essence defining a numerical flux, allowing the use of higher order numerical schemes. However, the discretization of the diffusion–dispersion terms have to be carefully analysed. In particular, time-step restrictions linked to the nature of the solute equation itself as well as the numerical diffusion associated to the numerical scheme used are question of interest in this work. These improvements are tested in an analytical case as well as in a laboratory test case with a passive solute (fluorescein) released from a reservoir. Experimental measurements are compared against the numerical results obtained with the proposed model and a sensitivity analysis is carried out, confirming an agreement with the longitudinal coefficients and an underestimation of the transversal ones, respectively.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{10.2166/hydro.2018.085,

title = {Calibration of a dynamic Eulerian-lagrangian model for the computation of wood cylinders transport in shallow water flow},

author = {Elisabetta Persi and Gabriella Petaccia and Stefano Sibilla and Pilar Brufau and Pilar García-Navarro},

url = {https://doi.org/10.2166/hydro.2018.085},

doi = {10.2166/hydro.2018.085},

issn = {1464-7141},

year  = {2018},

date = {2018-01-01},

journal = {Journal of Hydroinformatics},

volume = {21},

number = {1},

pages = {164-179},

abstract = {A computational Eulerian–Lagrangian model (ORSA2D_WT) is used for modelling the movement of floating rigid bodies on the water surface. The two-dimensional transport is computed with a dynamic approach, modifying existing formulations for the transport of bodies within fluid flows for the case of floating bodies, by adopting suitable added mass, drag and side coefficients. An original formulation for planar rotation is proposed, which includes the effect of the hydrodynamic torque and a resistance term, named added inertia, based on the difference between the angular velocity of the flow and that of the body. The value of the added inertia coefficient is calibrated against experiments made on purpose, involving the transport of a cylinder in a flume with two side obstacles. The calibrated code is applied to a slightly larger set of experiments for its preliminary evaluation. The outcome of the simulations shows that the streamwise and transversal displacements are well modelled, while some inaccuracies arise when considering the cylinder orientation. The effects of the initial conditions on the cylinders' trajectory and rotation are discussed, showing their influence on the evolution of the rotation angles.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{MORALESHERNANDEZ2018301,

title = {Assessing zebra mussel colonization of collective pressurized irrigation networks through pressure measurements and simulations},

author = {Mario Morales-Hernández and Enrique Playán and Yolanda Gimeno and Alfredo Serreta and Nery Zapata},

url = {http://www.sciencedirect.com/science/article/pii/S0378377418304219},

doi = {https://doi.org/10.1016/j.agwat.2018.04.025},

issn = {0378-3774},

year  = {2018},

date = {2018-01-01},

journal = {Agricultural Water Management},

volume = {204},

pages = {301 - 313},

abstract = {Zebra mussel (Dreissena polymorpha) colonies are becoming a real problem in pressurized irrigation networks. The zebra mussel infestation of the 45 Water Users Associations (WUAs) of the Riegos del Alto Aragón (RAA) irrigation project (121 thousand hectares located in northeastern Spain) was assessed during the period of 2013��?017. Maps of WUA infestation stages were produced. A survey of the WUAs made it possible to assess the relevance of certain structural and management practices in the control of zebra mussels. A method to monitor zebra mussel colonization of pressurized collective irrigation networks was presented. The method is based on the combination of pressure measurements at network hydrants and hydraulic simulations. Normalized pressure, estimated as the difference between simulated and observed pressure, should approach zero in all hydrants in a properly characterized, non-infested network. A positive normalized pressure can indicate the presence of zebra mussel colonies. The methodology was validated using two different test cases located in two RAA WUAs: the first case involved a discrete chemical treatment, while the second case was based on the analysis of three years of telemetry pressure data and remote operation of network hydrants. The existence of an infested reservoir upstream of the WUAs was the most likely source of zebra mussel colonization of the WUA pressurized networks in the RAA project. The desiccation and chemical treatment of the small WUA reservoirs was associated with pest control. The hybrid (measurement-simulation) methodology is able to characterize the presence of zebra mussel colonies in specific reaches of pressurized irrigation networks.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{camara3D,

title = {Towards transient experimental water surfaces: A new benchmark dataset for 2D shallow water solvers},

author = {S. Martínez-Aranda and J. Fernández-Pato and D. Caviedes-Voullième and I. García-Palacín and P. García-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S0309170818303658},

doi = {10.1016/j.advwatres.2018.08.013},

issn = {0309-1708},

year  = {2018},

date = {2018-01-01},

journal = {Advances in Water Resources},

volume = {121},

pages = {130 - 149},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{NAVASMONTILLA2018316,

title = {2D well-balanced augmented ADER schemes for the Shallow Water Equations with bed elevation and extension to the rotating frame},

author = {A. Navas-Montilla and J. Murillo},

url = {http://www.sciencedirect.com/science/article/pii/S0021999118304170},

doi = {https://doi.org/10.1016/j.jcp.2018.06.039},

issn = {0021-9991},

year  = {2018},

date = {2018-01-01},

journal = {Journal of Computational Physics},

volume = {372},

pages = {316 - 348},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{geosciences8080288,

title = {Development of a New Simulation Tool Coupling a 2D Finite Volume Overland Flow Model and a Drainage Network Model},

author = {Javier Fernández-Pato and Pilar García-Navarro},

url = {http://www.mdpi.com/2076-3263/8/8/288},

doi = {10.3390/geosciences8080288},

issn = {2076-3263},

year  = {2018},

date = {2018-01-01},

journal = {Geosciences},

volume = {8},

number = {8},

pages = {288},

abstract = {Numerical simulation of mixed flows combining free surface and pressurized flows is a practical tool to prevent possible flood situations in urban environments. When dealing with intense storm events, the limited capacity of the drainage network conduits can cause undesirable flooding situations. Computational simulation of the involved processes can lead to better management of the drainage network of urban areas. In particular, it is interesting to simultaneuously calculate the possible pressurization of the pipe network and the surface water dynamics in case of overflow. In this work, the coupling of two models is presented. The surface flow model is based on two-dimensional shallow water equations with which it is possible to solve the overland water dynamics as well as the transformation of rainfall into runoff through different submodels of infiltration. The underground drainage system assumes mostly free surface flow that can be pressurized in specific situations. The pipe network is modeled by means of one-dimensional sections coupled with the surface model in specific regions of the domain, such as drains or sewers. The numerical techniques considered for the resolution of both mathematical models are based on finite volume schemes with a first-order upwind discretization. The coupling of the models is verified using laboratory experimental data. Furthermore, the potential usefulness of the approach is demonstrated using real flooding data in a urban environment.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{CaviedesFernandezPato2018,

title = {Cellular Automata and Finite Volume solvers converge for 2D shallow flow modelling for hydrological modelling},

author = {Daniel Caviedes-Voullième and Javier Fernández-Pato and Christoph Hinz},

url = {http://www.sciencedirect.com/science/article/pii/S0022169418304438},

doi = {10.1016/j.jhydrol.2018.06.021},

issn = {0022-1694},

year  = {2018},

date = {2018-01-01},

journal = {Journal of Hydrology},

volume = {563},

pages = {411 - 417},

abstract = {Surface flows of hydrological interest, including overland flow, runoff, river and channel flow and flooding have received significant attention from modellers in the past 30 years. A growing effort to address these complex environmental problems is in place in the scientific community. Researchers have studied and favoured a plethora of techniques to approach this issue, ranging from very simple empirically-based mathematical models, to physically-based, deductive and very formal numerical integration of systems of partial-differential equations. In this work, we review two families of methods: cell-based simulators ��?later called Cellular Automata ��?and Finite Volume solvers for the Zero-Inertia equation, which we show to converge into a single methodology given appropriate choices. Furthermore, this convergence, mathematically shown in this work, can also be identified by critically reviewing the existing literature, which leads to the conclusion that two methods originating from different reasoning and fundamental philosophy, fundamentally converge into the same method. Moreover, acknowledging such convergence allows for some generalisation of properties of numerical schemes such as error behaviour and stability, which, importantly, is the same for the converging methodology, a fact with practical implications. Both the review of existing literature and reasoning in this work attempts to aid in the effort of synchronising and cross-fertilizing efforts to improve the understanding and the outlook of Zero-Inertia solvers for surface flows, as well as to help in clarifying the possible confusion and parallel developments that may arise from the use of different terminology originating from historical reasons. Moreover, synchronising and unifying this knowledge-base can help clarify model capabilities, applicability and modelling issues for hydrological modellers, specially for those not deeply familiar with the mathematical and numerical details.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{FernandezPato2018a,

title = {A fractional-order infiltration model to improve the simulation of rainfall/runoff in combination with a 2D Shallow Water model},

author = {J. Fernández-Pato and J. L. Gracia and P. García-Navarro},

doi = {10.2166/hydro.2018.145},

year  = {2018},

date = {2018-01-01},

journal = {Journal of Hydroinformatics},

volume = {20},

number = {4},

pages = {898-916},

abstract = {In this work, a distributed 2D Shallow Water (SW) flow model is combined with a fractional-order version of the Green-Ampt (FOGA) infiltration law to improve rainfall/runoff simulation in real catchments. The surface water model is based on a robust finite volume method on triangular grids that can handle flow over dry bed and multiple wet/dry fronts. When supplied with adequate infiltration laws, this model can provide useful information in surface Hydrology. The classical Green-Ampt law is generalized by using a Caputo fractional derivative of order less than or equal to 1 in Darcy's law. The novelty of this combination is that, on the one hand, the distributed SW simulation provides a detailed surface water distribution and, on the other hand, the FOGA model offers the possibility to model infiltration rates not monotonically decreasing. In order to obtain the best results, a non-uniform order of the fractional derivative depending on the cumulative infiltration and the existence of available surface water is proposed for realistic cases. This allows to improve significantly previous published numerical results in the literature for several storm events on catchments where the infiltration process occurs.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{doi:10.1080/00221686.2017.1300196,

title = {Application of an adjoint-based optimization procedure for the optimal control of internal boundary conditions in the shallow water equations},

author = {Asier Lacasta and Mario Morales-Hernández and Pilar Brufau and Pilar García-Navarro},

url = {https://doi.org/10.1080/00221686.2017.1300196},

doi = {10.1080/00221686.2017.1300196},

year  = {2018},

date = {2018-01-01},

journal = {Journal of Hydraulic Research},

volume = {56},

number = {1},

pages = {111-123},

publisher = {Taylor & Francis},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{Juez2018,

title = {Application of a distributed 2D overland flow model for a rainfall/runoff and erosion simulation in a Mediterranean watershed},

author = {C. Juez and A. Tena and J. Fernández-Pato and R. J. Batalla and P. García-Navarro},

doi = {https://publicaciones.unirioja.es/ojs/index.php/cig/issue/view/196},

year  = {2018},

date = {2018-01-01},

journal = {Cuadernos de Investigación Geográfica},

volume = {44},

number = {2},

pages = {615 - 640},

abstract = {Soil erosion has reemerged as an environmental problem associated with climate change that 

requires the help of simulation tools for forecasting future consequences. This topic becomes even more 

relevant in Mediterranean catchments due to the highly variable and irregular rainfall regime. Hence, an 

approach that includes the rainfall/runoff and erosion phenomena is required for quantifying the amount 

of soil the catchments are transferring to the rivers. As the calibration process of the infiltration and erosion 

parameters can become cumbersome in terms of iterations to the optimal values to fit experimental data, a 

Simplified Catchment Model (SCM) is introduced as a first approach. The set of tuning constants that 

provides the best fit are used as input for re-calibrating the parameters by means of the simulation of the 

real catchment. The modeling effort here presented opens its application to the analysis of the hydro-sedimentary processes at larger temporal and spatial scales.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{FernandezPato2017a,

title = {Implicit finite volume simulation of 2D shallow water flows in flexible meshes},

author = {J. Fernández-Pato and M. Morales-Hernández and P. García-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S0045782517302724},

doi = {10.1016/j.cma.2017.08.050},

year  = {2018},

date = {2018-01-01},

journal = {Computer Methods in Applied Mechanics and Engineering},

volume = {328},

pages = {1 - 25},

abstract = {In this work, an implicit method for solving 2D hyperbolic systems of equations is presented, focusing on the application to the 2D shallow water equations. It is based on the first order Roe’s scheme, in the framework of finite volume methods. A conservative linearization is done for the flux terms, leading to a non-structured matrix for unstructured meshes thus requiring iterative methods for solving the system. The validation is done by comparing numerical and exact solutions in both unsteady and steady cases. In order to test the applicability of the implicit scheme to real world situations, a laboratory scale tsunami simulation is carried out and compared to the experimental data. The implicit schemes have the advantage of the unconditional stability, but a quality loss in the transient solution can appear for high CFL numbers. The properties of the scheme are well suited for the simulation of unsteady shallow water flows over irregular topography using all kind of meshes.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{Juez2017,

title = {Experimental and numerical simulation of bed load transport over steep slopes},

author = {Carmelo Juez and Sandra Soares-Frazao and Javier Murillo and Pilar García-Navarro},

url = {https://doi.org/10.1080/00221686.2017.1288417},

doi = {10.1080/00221686.2017.1288417},

year  = {2017},

date = {2017-01-01},

journal = {Journal of Hydraulic Research},

volume = {55},

number = {4},

pages = {455-469},

publisher = {Taylor & Francis},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{NavasMontilla2017575,

title = {Overcoming numerical shockwave anomalies using energy balanced numerical schemes. Application to the Shallow Water Equations with discontinuous topography},

author = {A. Navas-Montilla and J. Murillo},

url = {http://www.sciencedirect.com/science/article/pii/S0021999117302589},

doi = {https://doi.org/10.1016/j.jcp.2017.03.057},

issn = {0021-9991},

year  = {2017},

date = {2017-01-01},

journal = {Journal of Computational Physics},

volume = {340},

pages = {575 - 616},

abstract = {When designing a numerical scheme for the resolution of conservation laws, the selection of a particular source term discretization (STD) may seem irrelevant whenever it ensures convergence with mesh refinement, but it has a decisive impact on the solution. In the framework of the Shallow Water Equations (SWE), well-balanced STD based on quiescent equilibrium are unable to converge to physically based solutions, which can be constructed considering energy arguments. Energy based discretizations can be designed assuming dissipation or conservation, but in any case, the STD procedure required should not be merely based on ad hoc approximations. The STD proposed in this work is derived from the Generalized Hugoniot Locus obtained from the Generalized Rankine Hugoniot conditions and the Integral Curve across the contact wave associated to the bed step. In any case, the STD must allow energy-dissipative solutions: steady and unsteady hydraulic jumps, for which some numerical anomalies have been documented in the literature. These anomalies are the incorrect positioning of steady jumps and the presence of a spurious spike of discharge inside the cell containing the jump. The former issue can be addressed by proposing a modification of the energy-conservative STD that ensures a correct dissipation rate across the hydraulic jump, whereas the latter is of greater complexity and cannot be fixed by simply choosing a suitable STD, as there are more variables involved. The problem concerning the spike of discharge is a well-known problem in the scientific community, also known as slowly-moving shock anomaly, it is produced by a nonlinearity of the Hugoniot locus connecting the states at both sides of the jump. However, it seems that this issue is more a feature than a problem when considering steady solutions of the SWE containing hydraulic jumps. The presence of the spurious spike in the discharge has been taken for granted and has become a feature of the solution. Even though it does not disturb the rest of the solution in steady cases, when considering transient cases it produces a very undesirable shedding of spurious oscillations downstream that should be circumvented. Based on spike-reducing techniques (originally designed for homogeneous Euler equations) that propose the construction of interpolated fluxes in the untrustworthy regions, we design a novel Roe-type scheme for the SWE with discontinuous topography that reduces the presence of the aforementioned spurious spike. The resulting spike-reducing method in combination with the proposed STD ensures an accurate positioning of steady jumps, provides convergence with mesh refinement, which was not possible for previous methods that cannot avoid the spike.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{Lacasta2017BoundaryAdjoint,

title = {Application of an adjoint-based optimization procedure for the optimal control of internal boundary conditions in the shallow water equations},

author = {Asier Lacasta and Mario Morales-Hernández and Pilar Brufau and Pilar García-Navarro},

url = {http://dx.doi.org/10.1080/00221686.2017.1300196},

doi = {10.1080/00221686.2017.1300196},

year  = {2017},

date = {2017-01-01},

journal = {Journal of Hydraulic Research},

volume = {0},

number = {0},

pages = {1-13},

abstract = {The shallow water equations have been extensively studied and used to model unsteady open channel flows in different applications. They belong to the category of hyperbolic partial differential equations and their treatment has recently been benefited from many numerical contributions leading to robust, accurate and well-balanced solutions. The correct formulation of external and internal boundary conditions is required to achieve a useful model in practical application. Control of internal boundary conditions may be useful in water distribution facilities. Therefore, development of a control strategy based on the fully dynamical mathematical model becomes attractive and justified. This work is devoted to the implementation of an adjoint based sensitivity analysis for the control of sluice gates in open-channel flow, formulated as internal boundary conditions. This is one of the most complex tasks in multiple regulated water delivery situations. Based on gradient method optimizers, the control of the whole channel to satisfy different requirements at several points of the channel is discussed. One of the achievements in this work is that the whole optimization process is performed two orders of magnitude faster than in real time. Moreover, the numerical results show this promising technique is a feasible way to obtain a robust and efficient control method.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{MoralesHernandez2017294,

title = {A Large Time Step explicit scheme (CFL>1) on unstructured grids for 2D conservation laws: Application to the homogeneous shallow water equations},

author = {M. Morales-Hernández and A. Lacasta and J. Murillo and P. García-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S0307904X17301397},

doi = {http://dx.doi.org/10.1016/j.apm.2017.02.043},

issn = {0307-904X},

year  = {2017},

date = {2017-01-01},

journal = {Applied Mathematical Modelling},

volume = {47},

pages = {294 - 317},

abstract = {In this work, a Large Time Step (LTS) explicit finite volume scheme designed to allow CFL > 1 is applied to the numerical resolution of 2D scalar and systems of conservation laws on triangular grids. Based on the flux difference splitting formulation, a special concern is put on finding the way of packing the information to compute the numerical solution when working on unstructured grids. Not only the cell areas but also the length of the interfaces and their orientation are questions of interest to send the information from each edge or interface. The information to update the cell variables is computed according to the local average discrete velocity and the orientation of the edges of the cells involved. The performance of these ideas is tested and compared with the conventional explicit first order and second order schemes in academic configurations for the 2D linear scalar equation and for 2D systems of conservation laws (in particular the shallow water equations) without source terms. The LTS scheme is demonstrated to preserve or even gain accuracy and save computational time with respect to the first order scheme.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{Echeverribar2016,

title = {Simulación numérica con RiverFlow2D de posibles soluciones de mitigación de avenidas en el tramo medio del río Ebro},

author = {Isabel Echeverribar and Mario Morales-Hernández and Asier Lacasta and Pilar Brufau and Pilar García-Navarro},

url = {http://polipapers.upv.es/index.php/IA/article/download/6550/pdf},

doi = {10.4995/ia.2017.6550},

year  = {2017},

date = {2017-01-01},

journal = {Ingeniería del agua},

volume = {21},

number = {1},

pages = {53-70},

abstract = {En este trabajo se presenta un estudio de medidas orientadas a la mitigación de avenidas en el tramo medio del río Ebro: limpieza de vegetación del cauce, uso de zonas de inundación controlada y construcción o re-adaptación de motas. Para ello se utiliza el software RiverFlow2D que resuelve las ecuaciones conservativas del flujo de superficie libre con un método de volúmenes finitos realizando los cálculos sobre GPU. Se comparan los resultados con medidas en estaciones de aforo e información extraída de ortofotos. La medida más efectiva, de las analizadas, ha resultado ser la eliminación de la vegetación en el cauce. Se demuestra que no sólo el área máxima inundada es menor en todo el tramo sino que también reduce la altura de agua hasta en 1 m. El resto de medidas tienen consecuencias locales y de poca entidad cuando los caudales pico son altos, aunque podrían resultar de utilidad para avenidas con caudales más bajos.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{Lacastajh2017021,

title = {Calibration of the 1D shallow water equations: a comparison of Monte Carlo and gradient-based optimization methods},

author = {Asier Lacasta and Mario Morales-Hernández and Javier Burguete and Pilar Brufau and Pilar García-Navarro},

url = {http://jh.iwaponline.com/content/early/2017/01/06/hydro.2017.021},

doi = {10.2166/hydro.2017.021},

issn = {1464-7141},

year  = {2017},

date = {2017-01-01},

journal = {Journal of Hydroinformatics},

publisher = {IWA Publishing},

abstract = {The calibration of parameters in complex systems usually requires a large computational effort. Moreover, it becomes harder to perform the calibration when non-linear systems underlie the physical process, and the direction to follow in order to optimize an objective function changes depending on the situation. In the context of shallow water equations (SWE), the calibration of parameters, such as the roughness coefficient or the gauge curve for the outlet boundary condition, is often required. In this work, the SWE are used to simulate an open channel flow with lateral gates. Due to the uncertainty in the mathematical modeling that these lateral discharges may introduce into the simulation, the work is focused on the calibration of discharge coefficients. Thus, the calibration is performed by two different approaches. On the one hand, a classical Monte Carlo method is used. On the other hand, the development and application of an adjoint formulation to evaluate the gradient is presented. This is then used in a gradient-based optimizer and is compared with the stochastic approach. The advantages and disadvantages are illustrated and discussed through different test cases.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{Murillo201670,

title = {A comprehensive explanation and exercise of the source terms in hyperbolic systems using Roe type solutions. Application to the 1D-2D shallow water equations},

author = {J. Murillo and A. Navas-Montilla},

url = {http://www.sciencedirect.com/science/article/pii/S0309170816305917},

doi = {http://dx.doi.org/10.1016/j.advwatres.2016.10.019},

issn = {0309-1708},

year  = {2016},

date = {2016-01-01},

journal = {Advances in Water Resources},

volume = {98},

pages = {70 - 96},

abstract = {Abstract Powerful numerical methods have to consider the presence of source terms of different nature, that intensely compete among them and may lead to strong spatiotemporal variations in the flow. When applied to shallow flows, numerical preservation of quiescent equilibrium, also known as the well-balanced property, is still nowadays the keystone for the formulation of novel numerical schemes. But this condition turns completely insufficient when applied to problems of practical interest. Energy balanced methods (E-schemes) can overcome all type of situations in shallow flows, not only under arbitrary geometries, but also with independence of the rheological shear stress model selected. They must be able to handle correctly transient problems including modeling of starting and stopping flow conditions in debris flow and other flows with a non-Newtonian rheological behavior. The numerical solver presented here satisfies these properties and is based on an approximate solution defined in a previous work. Given the relevant capabilities of this weak solution, it is fully theoretically derived here for a general set of equations. This useful step allows providing for the first time an E-scheme, where the set of source terms is fully exercised under any flow condition involving high slopes and arbitrary shear stress. With the proposed solver, a Roe type first order scheme in time and space, positivity conditions are explored under a general framework and numerical simulations can be accurately performed recovering an appropriate selection of the time step, allowed by a detailed analysis of the approximate solver. The use of case-dependent threshold values is unnecessary and exact mass conservation is preserved.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{NavasMontilla2016108,

title = {Asymptotically and exactly energy balanced augmented flux-ADER schemes with application to hyperbolic conservation laws with geometric source terms},

author = {A. Navas-Montilla and J. Murillo},

url = {http://www.sciencedirect.com/science/article/pii/S0021999116301024},

doi = {http://dx.doi.org/10.1016/j.jcp.2016.04.047},

issn = {0021-9991},

year  = {2016},

date = {2016-01-01},

journal = {Journal of Computational Physics},

volume = {317},

pages = {108 - 147},

abstract = {Abstract In this work, an arbitrary order HLL-type numerical scheme is constructed using the flux-ADER methodology. The proposed scheme is based on an augmented Derivative Riemann solver that was used for the first time in Navas-Montilla and Murillo (2015) [1]. Such solver, hereafter referred to as Flux-Source (FS) solver, was conceived as a high order extension of the augmented Roe solver and led to the generation of a novel numerical scheme called AR-ADER scheme. Here, we provide a general definition of the FS solver independently of the Riemann solver used in it. Moreover, a simplified version of the solver, referred to as Linearized-Flux-Source (LFS) solver, is presented. This novel version of the FS solver allows to compute the solution without requiring reconstruction of derivatives of the fluxes, nevertheless some drawbacks are evidenced. In contrast to other previously defined Derivative Riemann solvers, the proposed FS and LFS solvers take into account the presence of the source term in the resolution of the Derivative Riemann Problem (DRP), which is of particular interest when dealing with geometric source terms. When applied to the shallow water equations, the proposed HLLS-ADER and AR-ADER schemes can be constructed to fulfill the exactly well-balanced property, showing that an arbitrary quadrature of the integral of the source inside the cell does not ensure energy balanced solutions. As a result of this work, energy balanced flux-ADER schemes that provide the exact solution for steady cases and that converge to the exact solution with arbitrary order for transient cases are constructed.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{Lacasta2016371,

title = {A GPU accelerated adjoint-based optimizer for inverse modeling of the two-dimensional shallow water equations},

author = {A. Lacasta and P. García-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S0045793016302092},

doi = {10.1016/j.compfluid.2016.06.024},

issn = {0045-7930},

year  = {2016},

date = {2016-01-01},

journal = {Computers & Fluids},

volume = {136},

pages = {371 - 383},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{doi:10.1080/00221686.2016.1143042,

title = {An efficient GPU implementation for a faster simulation of unsteady bed-load transport},

author = {Carmelo Juez and Asier Lacasta and Javier Murillo and Pilar García-Navarro},

url = {http://dx.doi.org/10.1080/00221686.2016.1143042},

doi = {10.1080/00221686.2016.1143042},

year  = {2016},

date = {2016-01-01},

journal = {Journal of Hydraulic Research},

volume = {54},

number = {3},

pages = {275-288},

abstract = {ABSTRACTComputational tools may help engineers in the assessment of sediment transport during the decision-making processes. The main requirements are that the numerical results have to be accurate and simulation models must be fast. The present work is based on the 2D shallow water equations in combination with the 2D Exner equation. The resulting numerical model accuracy was already discussed in previous work. Regarding the speed of the computation, the Exner equation slows down the already costly 2D shallow water model as the number of variables to solve is increased and the numerical stability is more restrictive. In order to reduce the computational effort required for simulating realistic scenarios, the authors have exploited the use of Graphics Processing Units in combination with non-trivial optimization procedures. The gain in computing cost obtained with the graphic hardware is compared against single-core (sequential) and multi-core (parallel) CPU implementations in two unsteady cases.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{FernándezPato2016b,

title = {2D Zero-Inertia Model for Solution of Overland Flow Problems in Flexible Meshes},

author = {J. Fernández-Pato and P. García-Navarro},

url = {http://dx.doi.org/10.1061/(ASCE)HE.1943-5584.0001428},

doi = {10.1061/(ASCE)HE.1943-5584.0001428},

year  = {2016},

date = {2016-01-01},

journal = {Journal of Hydrologic Engineering},

volume = {0},

number = {0},

pages = {04016038},

abstract = {A study of the efficiency of a Zero-Inertia model for 2D overland flow simulation is presented in this work. An upwind numerical scheme is used for the spatial discretization in the frame of finite volume methods and an implicit formulation is chosen to avoid numerical instability. The scheme is applied in both structured and unstructured meshes, focusing in the last ones, due to their good adaptability. The Zero-Inertia equation has a non-linear character, hence a linearization is required in the implicit procedure. This is carried out by means of Picard iterations method as a previous step to the system matrix resolution, characteristic of implicit techniques. The BiConjugate Gradient Stabilized (BiCGStab) method combined with sparse storage strategies is selected for the system resolution. A dual threshold Incomplete LU factorization (ILUT) is chosen as matrix preconditioner. Computational efficiency of the implicit temporal discretization for Zero-Inertia model is explored under both steady and unsteady flow conditions by comparing the CPU times against the explicit version of the same model.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{FernándezPato2016c,

title = {Simulación numérica de inundación de valles fluviales mediante un modelo difusivo implícito},

author = {Javier Fernández-Pato and Pilar García-Navarro},

url = {http://dx.doi.org/10.4995/ia.2016.4548},

doi = {10.4995/ia.2016.4548},

year  = {2016},

date = {2016-01-01},

journal = {Ingeniería del agua},

volume = {20},

number = {3},

pages = {115-126},

abstract = {En este trabajo se presenta un modelo difusivo de flujo superficial para la resolución eficiente de problemas de inundación de valles fluviales. La discretización espacial se realiza mediante un esquema upwind de volúmenes finitos, aplicado en una malla triangular no estructurada. Para la discretización temporal se ha empleado un esquema implícito, lo que conlleva la generación de un sistema de tantas ecuaciones como celdas de cálculo tenga la malla computacional. Para su resolución, se ha empleado el método del Gradiente Biconjugado Estabilizado (BiCGStab). La eficiencia computacional se mide realizando una comparación de tiempos de CPU de las versiones explícita e implícita del mismo esquema numérico aplicado al modelo de onda difusiva. El estudio demuestra que, en general, el modelo difusivo se beneficia de una discretización temporal implícita, aumentando su eficiencia en gran medida frente al uso de un esquema explícito. Debido al carácter no lineal de la ecuación difusiva, un paso de tiempo mayor no siempre garantiza un menor coste computacional y resulta necesario encontrar el paso de tiempo óptimo para cada problema. Se ha aplicado el modelo difusivo a un caso test de inundación de valle fluvial propuesto por la UK Environmetal Agency para evaluar su capacidad predictiva en relación a modelos comerciales. Por último, se han comparado los resultados de calado y nivel de agua proporcionados por el modelo propuesto con los generados por una formulación matemática basada en las ecuaciones completas de aguas poco profundas no encontrándose diferencias relevantes en el ejemplo analizado.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{SeguraBeltran2016,

title = {Using post-flood surveys and geomorphologic mapping to evaluate hydrological and hydraulic models: The flash flood of the Girona River (Spain) in 2007},

author = {F. Segura-Beltrán and C. Sanchis-Ibor and M. Morales-Hernández and M. González-Sanchis and G. Bussi and E. Ortiz},

url = {http://www.sciencedirect.com/science/article/pii/S002216941630230X},

doi = {10.1016/j.jhydrol.2016.04.039},

issn = {0022-1694},

year  = {2016},

date = {2016-01-01},

journal = {Journal of Hydrology},

pages = {-},

abstract = {Summary This paper analyzes the Girona River (Spain) flash flood, occurred on the 12th of October 2007, combining hydrological and hydraulic modeling with geomorphologic mapping and post-flood survey information. This research aims to reproduce the flood event in order to understand and decipher the flood processes and dynamics on a system of prograding alluvial fans. The hydrological model TETIS was used to characterize the shape and dimension of the October 2007 Girona River hydrograph. Subsequently, the flood event was reproduced using the free surface flow module of the model RiverFlow2D. The combination of hydrological and hydraulic models was evaluated using post-flood surveys defining maximum flooded area and flood depths. Then, simulations with different peak discharges were carried out to estimate the hydro-geomorphologic response of the Girona River floodplain, through the identification of the activation thresholds in different geomorphic elements. Results showed that the unit peak discharge of the October 2007 flood event (5 m3 s��? km��?) was among the largest ever recorded in the area, according to the existing literature. Likewise, the hydraulic model showed a good performance in reproducing the flood event (FitA = 76%, RMSE = 0.65 m and NSE = 0.6), despite the complexity of the case, an ephemeral and ungauged river. The model simulation revealed the existence of an activation pattern of paleochannels and alluvial fans, which was altered by the presence of some anthropogenic disturbances. This multidisciplinary approach proved to be a useful strategy for understanding flash flood processes in ungauged catchments. It allowed understanding the mechanisms governing floods in alluvial fans systems and it represented a solid contribution for early warning plans and risk mitigation policies.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{FernándezPato2016496,

title = {Rainfall/runoff simulation with 2D full shallow water equations: Sensitivity analysis and calibration of infiltration parameters},

author = {Javier Fernández-Pato and Daniel Caviedes-Voullième and Pilar García-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S0022169416301263},

doi = {doi:10.1016/j.jhydrol.2016.03.021},

issn = {0022-1694},

year  = {2016},

date = {2016-01-01},

journal = {Journal of Hydrology},

volume = {536},

pages = {496 - 513},

abstract = {Summary One of the most difficult issues in the development of hydrologic models is to find a rigorous source of data and specific parameters to a given problem, on a given location that enable reliable calibration. In this paper, a distributed and physically based model (2D Shallow Water Equations) is used for surface flow and runoff calculations in combination with two infiltration laws (Horton and Green–Ampt) for estimating infiltration in a watershed. This technique offers the capability of assigning a local and time-dependent infiltration rate to each computational cell depending on the available surface water, soil type or vegetation. We investigate how the calibration of parameters is affected by transient distributed Shallow Water model and the complexity of the problem. In the first part of this work, we calibrate the infiltration parameters for both Horton and Green–Ampt models under flat ponded soil conditions. Then, by means of synthetic test cases, we perform a space-distributed sensitivity analysis in order to show that this calibration can be significantly affected by the introduction of topography or rainfall. In the second part, parameter calibration for a real catchment is addressed by comparing the numerical simulations with two different sets of experimental data, corresponding to very different events in terms of the rainfall volume. We show that the initial conditions of the catchment and the rainfall pattern have a special relevance in the quality of the adjustment. Hence, it is shown that the topography of the catchment and the storm characteristics affect the calibration of infiltration parameters.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{NavasMontilla2016108b,

title = {Asymptotically and exactly energy balanced augmented flux-ADER schemes with application to hyperbolic conservation laws with geometric source terms},

author = {A. Navas-Montilla and J. Murillo},

url = {http://www.sciencedirect.com/science/article/pii/S0021999116301024},

doi = {http://dx.doi.org/10.1016/j.jcp.2016.04.047},

issn = {0021-9991},

year  = {2016},

date = {2016-01-01},

journal = {Journal of Computational Physics},

volume = {317},

pages = {108 - 147},

abstract = {Abstract In this work, an arbitrary order HLL-type numerical scheme is constructed using the flux-ADER methodology. The proposed scheme is based on an augmented Derivative Riemann solver that was used for the first time in Navas-Montilla and Murillo (2015) [1]. Such solver, hereafter referred to as Flux-Source (FS) solver, was conceived as a high order extension of the augmented Roe solver and led to the generation of a novel numerical scheme called AR-ADER scheme. Here, we provide a general definition of the FS solver independently of the Riemann solver used in it. Moreover, a simplified version of the solver, referred to as Linearized-Flux-Source (LFS) solver, is presented. This novel version of the FS solver allows to compute the solution without requiring reconstruction of derivatives of the fluxes, nevertheless some drawbacks are evidenced. In contrast to other previously defined Derivative Riemann solvers, the proposed FS and LFS solvers take into account the presence of the source term in the resolution of the Derivative Riemann Problem (DRP), which is of particular interest when dealing with geometric source terms. When applied to the shallow water equations, the proposed HLLS-ADER and AR-ADER schemes can be constructed to fulfill the exactly well-balanced property, showing that an arbitrary quadrature of the integral of the source inside the cell does not ensure energy balanced solutions. As a result of this work, energy balanced flux-ADER schemes that provide the exact solution for steady cases and that converge to the exact solution with arbitrary order for transient cases are constructed.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{Juez20161,

title = {A model based on Hirano-Exner equations for two-dimensional transient flows over heterogeneous erodible beds},

author = {C. Juez and C. Ferrer-Boix and J. Murillo and M. A. Hassan and P. García-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S0309170815002584},

doi = {10.1016/j.advwatres.2015.10.013},

issn = {0309-1708},

year  = {2016},

date = {2016-01-01},

journal = {Advances in Water Resources},

volume = {87},

pages = {1 - 18},

abstract = {Abstract In order to study the morphological evolution of river beds composed of heterogeneous material, the interaction among the different grain sizes must be taken into account. In this paper, these equations are combined with the two-dimensional shallow water equations to describe the flow field. The resulting system of equations can be solved in two ways: (i) in a coupled way, solving flow and sediment equations simultaneously at a given time-step or (ii) in an uncoupled manner by first solving the flow field and using the magnitudes obtained at each time-step to update the channel morphology (bed and surface composition). The coupled strategy is preferable when dealing with strong and quick interactions between the flow field, the bed evolution and the different particle sizes present on the bed surface. A number of numerical difficulties arise from solving the fully coupled system of equations. These problems are reduced by means of a weakly-coupled strategy to numerically estimate the wave celerities containing the information of the bed and the grain sizes present on the bed. Hence, a two-dimensional numerical scheme able to simulate in a self-stable way the unsteady morphological evolution of channels formed by cohesionless grain size mixtures is presented. The coupling technique is simplified without decreasing the number of waves involved in the numerical scheme but by simplifying their definitions. The numerical results are satisfactorily tested with synthetic cases and against experimental data.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{MoralesHernández20162087,

title = {Conservative 1D–2D coupled numerical strategies applied to river flooding: The Tiber (Rome)},

author = {M. Morales-Hernández and G. Petaccia and P. Brufau and P. García-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S0307904X15005375},

doi = {10.1016/j.apm.2015.08.016},

issn = {0307-904X},

year  = {2016},

date = {2016-01-01},

journal = {Applied Mathematical Modelling},

volume = {40},

number = {3},

pages = {2087 - 2105},

abstract = {Abstract Coupled 1D�D numerical strategies are presented in this work for their application to fast computation of large rivers flooding. Both 1D and 2D models are built using explicit upwind finite volume schemes, able to deal with wetting-drying fronts. The topography representation is described via cross sections for the 1D model and with quadrilateral/triangular structured/unstructured meshes for the 2D model. The coupling strategies, free of hydraulic structures and tuning parameters, are firstly validated in a laboratory test dealing with a levee break and its flooding into a lateral plane. The numerical results are compared with a fully 2D model as well as with measurements in some gauge points giving satisfactory results. The simulation of a real flooding scenario in the Tiber river near the urban area of Rome (Italy) is then performed. A lateral coupling configuration is provided, in which the flood wave propagation in the main channel is simulated by means of a 1D model and the inundation of the riverside is simulated by means of a 2D model. On the other hand, a frontal coupling, in which the flood wave is simulated in a 1D model first and then it is propagated into a 2D model, is also performed. The flooding extension is almost well captured by all the schemes presented, being the 1D�D lateral configuration the most confident with speed-ups of around 15x.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{FernándezPato2016b,

title = {2D Zero-Inertia Model for Solution of Overland Flow Problems in Flexible Meshes},

author = {J. Fernández-Pato and P. García-Navarro},

url = {http://dx.doi.org/10.1061/(ASCE)HE.1943-5584.0001428},

doi = {10.1061/(ASCE)HE.1943-5584.0001428},

year  = {2016},

date = {2016-01-01},

journal = {Journal of Hydrologic Engineering},

volume = {0},

number = {0},

pages = {04016038},

abstract = {A study of the efficiency of a Zero-Inertia model for 2D overland flow simulation is presented in this work. An upwind numerical scheme is used for the spatial discretization in the frame of finite volume methods and an implicit formulation is chosen to avoid numerical instability. The scheme is applied in both structured and unstructured meshes, focusing in the last ones, due to their good adaptability. The Zero-Inertia equation has a non-linear character, hence a linearization is required in the implicit procedure. This is carried out by means of Picard iterations method as a previous step to the system matrix resolution, characteristic of implicit techniques. The BiConjugate Gradient Stabilized (BiCGStab) method combined with sparse storage strategies is selected for the system resolution. A dual threshold Incomplete LU factorization (ILUT) is chosen as matrix preconditioner. Computational efficiency of the implicit temporal discretization for Zero-Inertia model is explored under both steady and unsteady flow conditions by comparing the CPU times against the explicit version of the same model.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{FernándezPato2016c,

title = {Simulación numérica de inundación de valles fluviales mediante un modelo difusivo implícito},

author = {Javier Fernández-Pato and Pilar García-Navarro},

url = {http://dx.doi.org/10.4995/ia.2016.4548},

doi = {10.4995/ia.2016.4548},

year  = {2016},

date = {2016-01-01},

journal = {Ingeniería del agua},

volume = {20},

number = {3},

pages = {115-126},

abstract = {En este trabajo se presenta un modelo difusivo de flujo superficial para la resolución eficiente de problemas de inundación de valles fluviales. La discretización espacial se realiza mediante un esquema upwind de volúmenes finitos, aplicado en una malla triangular no estructurada. Para la discretización temporal se ha empleado un esquema implícito, lo que conlleva la generación de un sistema de tantas ecuaciones como celdas de cálculo tenga la malla computacional. Para su resolución, se ha empleado el método del Gradiente Biconjugado Estabilizado (BiCGStab). La eficiencia computacional se mide realizando una comparación de tiempos de CPU de las versiones explícita e implícita del mismo esquema numérico aplicado al modelo de onda difusiva. El estudio demuestra que, en general, el modelo difusivo se beneficia de una discretización temporal implícita, aumentando su eficiencia en gran medida frente al uso de un esquema explícito. Debido al carácter no lineal de la ecuación difusiva, un paso de tiempo mayor no siempre garantiza un menor coste computacional y resulta necesario encontrar el paso de tiempo óptimo para cada problema. Se ha aplicado el modelo difusivo a un caso test de inundación de valle fluvial propuesto por la UK Environmetal Agency para evaluar su capacidad predictiva en relación a modelos comerciales. Por último, se han comparado los resultados de calado y nivel de agua proporcionados por el modelo propuesto con los generados por una formulación matemática basada en las ecuaciones completas de aguas poco profundas no encontrándose diferencias relevantes en el ejemplo analizado.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{FernándezPato2016496,

title = {Rainfall/runoff simulation with 2D full shallow water equations: Sensitivity analysis and calibration of infiltration parameters},

author = {Javier Fernández-Pato and Daniel Caviedes-Voullième and Pilar García-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S0022169416301263},

doi = {doi:10.1016/j.jhydrol.2016.03.021},

issn = {0022-1694},

year  = {2016},

date = {2016-01-01},

journal = {Journal of Hydrology},

volume = {536},

pages = {496 - 513},

abstract = {Summary One of the most difficult issues in the development of hydrologic models is to find a rigorous source of data and specific parameters to a given problem, on a given location that enable reliable calibration. In this paper, a distributed and physically based model (2D Shallow Water Equations) is used for surface flow and runoff calculations in combination with two infiltration laws (Horton and Green–Ampt) for estimating infiltration in a watershed. This technique offers the capability of assigning a local and time-dependent infiltration rate to each computational cell depending on the available surface water, soil type or vegetation. We investigate how the calibration of parameters is affected by transient distributed Shallow Water model and the complexity of the problem. In the first part of this work, we calibrate the infiltration parameters for both Horton and Green–Ampt models under flat ponded soil conditions. Then, by means of synthetic test cases, we perform a space-distributed sensitivity analysis in order to show that this calibration can be significantly affected by the introduction of topography or rainfall. In the second part, parameter calibration for a real catchment is addressed by comparing the numerical simulations with two different sets of experimental data, corresponding to very different events in terms of the rainfall volume. We show that the initial conditions of the catchment and the rainfall pattern have a special relevance in the quality of the adjustment. Hence, it is shown that the topography of the catchment and the storm characteristics affect the calibration of infiltration parameters.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{MoralesHernández20162087,

title = {Conservative 1D–2D coupled numerical strategies applied to river flooding: The Tiber (Rome)},

author = {M. Morales-Hernández and G. Petaccia and P. Brufau and P. García-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S0307904X15005375},

doi = {10.1016/j.apm.2015.08.016},

issn = {0307-904X},

year  = {2016},

date = {2016-01-01},

journal = {Applied Mathematical Modelling},

volume = {40},

number = {3},

pages = {2087 - 2105},

abstract = {Abstract Coupled 1D��?D numerical strategies are presented in this work for their application to fast computation of large rivers flooding. Both 1D and 2D models are built using explicit upwind finite volume schemes, able to deal with wetting-drying fronts. The topography representation is described via cross sections for the 1D model and with quadrilateral/triangular structured/unstructured meshes for the 2D model. The coupling strategies, free of hydraulic structures and tuning parameters, are firstly validated in a laboratory test dealing with a levee break and its flooding into a lateral plane. The numerical results are compared with a fully 2D model as well as with measurements in some gauge points giving satisfactory results. The simulation of a real flooding scenario in the Tiber river near the urban area of Rome (Italy) is then performed. A lateral coupling configuration is provided, in which the flood wave propagation in the main channel is simulated by means of a 1D model and the inundation of the riverside is simulated by means of a 2D model. On the other hand, a frontal coupling, in which the flood wave is simulated in a 1D model first and then it is propagated into a 2D model, is also performed. The flooding extension is almost well captured by all the schemes presented, being the 1D��?D lateral configuration the most confident with speed-ups of around 15x.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{lacasta2015calibracion,

title = {Calibración y simulación de un sistema regulado de suministro de agua a través de técnicas de Monte Carlo},

author = {A Lacasta and M Morales-Hernández and M Tejero-Juste and J Burguete and P Brufau and P García-Navarro},

url = {http://polipapers.upv.es/index.php/IA/article/view/3350},

doi = {10.4995/ia.2015.3350},

issn = {1886-4996},

year  = {2015},

date = {2015-01-01},

journal = {Ingeniería del agua},

volume = {19},

number = {3},

pages = {117–133},

abstract = {El flujo en canales se caracteriza por cubrir largas distancias y obedecer a patrones temporales variables. Suele estar regulado por elementos hidráulicos, como compuertas laterales, para asegurar un correcto abastecimiento de agua. La dinámica de este flujo viene gobernada por un sistema de ecuaciones diferenciales en derivadas parciales llamado de aguas poco profundas. Junto a ellas es necesario establecer una formulación simplificada del funcionamiento de las compuertas. El conjunto forma un sistema no lineal que sólo se puede resolver numéricamente. Aquí se propone un esquema descentrado explícito de volúmenes finitos con el fin de resolver todo tipo de regímenes. La modelización de las estructuras hidráulicas (compuertas laterales) introduce parámetros con incertidumbre. Por ello, serán calibradas mediante el algoritmo de Monte Carlo, obteniendo como resultado unos coeficientes asociados a cada una de ellas. Posteriormente serán verificadas utilizando casos reales proporcionados por el equipo de monitorización del canal de Pina de Ebro (Zaragoza).},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{NavasMontilla2015188,

title = {Energy balanced numerical schemes with very high order. The Augmented Roe Flux ADER scheme. Application to the shallow water equations},

author = {A. Navas-Montilla and J. Murillo},

url = {http://www.sciencedirect.com/science/article/pii/S0021999115001217},

doi = {10.1016/j.jcp.2015.03.002},

issn = {0021-9991},

year  = {2015},

date = {2015-01-01},

journal = {Journal of Computational Physics},

volume = {290},

pages = {188 - 218},

abstract = {Abstract In this work, an ADER type finite volume numerical scheme is proposed as an extension of a first order solver based on weak solutions of RPs with source terms. The type of source terms considered here are a special but relevant type of source terms: their spatial integral is discontinuous. The relevant difference with other previously defined ADER schemes is that it considers the presence of the source term in the solutions of the DRP. Unlike the original ADER schemes, the proposed numerical scheme computes the RPs of the high order terms of the DRP departing from time derivatives of the fluxes as initial conditions for these RPs. Weak solutions of the RPs defined for the DRP are computed using an augmented version of the Roe solver that includes an extra wave that accounts for the contribution of the source term. The discretization done over the source term leads to an energy balanced numerical scheme that allows to obtain the exact solution for steady cases with independence of the grid refinement. In unsteady problems, the numerical scheme ensures the convergence to the exact solution. The numerical scheme is constructed with an arbitrary order of accuracy, and has no theoretical barrier. Numerical results for the Burger's equation and the shallow water equations are presented in this work and indicate that the proposed numerical scheme is able to converge with the expected order of accuracy.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{alacastagpuhydro,

title = {GPU implementation of the 2D shallow water equations for the simulation of rainfall/runoff events},

author = {Asier Lacasta and Mario Morales-Hernández and Javier Murillo and Pilar García-Navarro},

url = {http://dx.doi.org/10.1007/s12665-015-4215-z},

doi = {10.1007/s12665-015-4215-z},

issn = {1866-6280},

year  = {2015},

date = {2015-01-01},

journal = {Environmental Earth Sciences},

volume = {74},

number = {11},

pages = {7295-7305},

publisher = {Springer Berlin Heidelberg},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{MoralesHernandezRCE1D2D,

title = {A Riemann coupled edge (RCE) 1D��?D finite volume inundation and solute transport model},

author = {M. Morales-Hernández and A. Lacasta and J. Murillo and P. Brufau and P. García-Navarro},

url = {http://dx.doi.org/10.1007/s12665-015-4754-3},

doi = {10.1007/s12665-015-4754-3},

issn = {1866-6280},

year  = {2015},

date = {2015-01-01},

journal = {Environmental Earth Sciences},

volume = {74},

number = {11},

pages = {7319-7335},

publisher = {Springer Berlin Heidelberg},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{Lacasta201563,

title = {An efficient solution for hazardous geophysical flows simulation using GPUs},

author = {A. Lacasta and C. Juez and J. Murillo and P. Garcia-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S0098300415000321},

doi = {10.1016/j.cageo.2015.02.010},

issn = {0098-3004},

year  = {2015},

date = {2015-01-01},

journal = {Computers & Geosciences},

volume = {78},

number = {0},

pages = {63 - 72},

abstract = {Abstract The movement of poorly sorted material over steep areas constitutes a hazardous environmental problem. Computational tools help in the understanding and predictions of such landslides. The main drawback is the high computational effort required for obtaining accurate numerical solutions due to the high number of cells involved in the calculus. In order to overcome this problem, this work proposes the use of GPUs for decreasing significantly the CPU simulation time. The numerical scheme implemented in GPU is based on a finite volume scheme and it was validated in previous work with exact solutions and experimental data. The computational cost time obtained with the Graphical Hardware technology, GPU, is compared against Single-Core (sequential) and Multi-Core (parallel) CPU implementations. The GPU implementation allows to reduce the computational cost time in two orders of magnitude.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{Murillo2015149,

title = {A Roe type energy balanced solver for 1D arterial blood flow and transport},

author = {J. Murillo and P. Garcia-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S0045793015001504},

doi = {10.1016/j.compfluid.2015.05.003},

issn = {0045-7930},

year  = {2015},

date = {2015-01-01},

journal = {Computers & Fluids},

volume = {117},

number = {0},

pages = {149 - 167},

abstract = {Abstract The approximate solver presented in this work is based on the upwind discretization of the source terms and a genuinely Roe solver for the one-dimensional blood flow equations in arteries. This augmented solver involves the presence of the source terms, ensuring convergence to the exact solution by including an extra wave associated to the change in the material properties and the friction term. The resulting numerical scheme is energy-balanced, that is, ensures equilibrium in rest conditions and is able to ensure numerically a constant level of energy in steady cases with velocity. The resulting numerical solver allows simulating directly mass transport without creating non-physical oscillations. The numerical scheme is assessed using steady and unsteady problems with exact solutions and is compared with models of the systemic arterial tree published in literature.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{HYP:HYP10152,

title = {Modelling sediment deposition and phosphorus retention in a river floodplain},

author = {M. González-Sanchis and J. Murillo and A. Cabezas and J. E. Vermaat and F. A. Comin and P. Garcia-Navarro},

url = {http://dx.doi.org/10.1002/hyp.10152},

doi = {10.1002/hyp.10152},

issn = {1099-1085},

year  = {2015},

date = {2015-01-01},

journal = {Hydrological Processes},

volume = {29},

number = {3},

pages = {384-394},

abstract = {Phosphorus (P) is one of the major limiting nutrient in many freshwater ecosystems. During the last decade, attention has been focused on the fluxes of suspended sediment and particulate P through freshwater drainage systems because of severe eutrophication effects in aquatic ecosystems. Hence, the analysis and prediction of phosphorus and sediment dynamics constitute an important element for ecological conservation and restoration of freshwater ecosystems. In that sense, the development of a suitable prediction model is justified, and the present work is devoted to the validation and application of a predictive soluble reactive phosphorus (SRP) uptake and sedimentation models, to a real riparian system of the middle Ebro river floodplain. Both models are coupled to a fully distributed two-dimensional shallow-water flow numerical model. The SRP uptake model is validated using data from three field experiments. The model predictions show a good accuracy for SRP concentration, where the linear regressions between measured and calculated values of the three experiments were significant (r2 ≥��?.62; p ≤��?.05), and a Nash–Sutcliffe coefficient (E) that ranged from 0.54 to 0.62. The sedimentation model is validated using field data collected during two real flooding events within the same river reach. The comparison between calculated and measured sediment depositions showed a significant linear regression (p ≤��?.05; r2��?��?.97) and an E that ranged from 0.63 to 0.78. Subsequently, the complete model that includes flow dynamics, solute transport, SRP uptake and sedimentation is used to simulate and analyse floodplain sediment deposition, river nutrient contribution and SRP uptake. According to this analysis, the main SRP uptake process appears to be the sediment sorption. The analysis also reveals the presence of a lateral gradient of hydrological connectivity that decreases with distance from the river and controls the river matter contribution to the floodplain. Copyright © 2014 John Wiley & Sons, Ltd.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{CaviedesVoullième2014177,

title = {2D dry granular free-surface flow over complex topography with obstacles. Part I: experimental study using a consumer-grade RGB-D sensor},

author = {Daniel Caviedes-Voullième and Carmelo Juez and Javier Murillo and Pilar García-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S0098300414002210},

doi = {10.1016/j.cageo.2014.09.009},

issn = {0098-3004},

year  = {2014},

date = {2014-01-01},

journal = {Computers & Geosciences},

volume = {73},

number = {0},

pages = {177 - 197},

abstract = {Abstract Avalanches, debris flows and other types of gravity-driven granular flows are a common hazard in mountainous regions. These regions often have human settlements in the lower parts of valleys, with human structures dangerously exposed to the destructive effects of these geophysical flows. Therefore a scientific effort has been made to understand, model and simulate geophysical granular flows. In order for computer models and simulations to be of predictive value they need to be validated under controlled, yet nature-like conditions. This work presents an experimental study of granular flow over a simplified mountain slope and valley topography. The experimental facility has a rough bed with very high slope at the upstream end and adverse slope on the downstream end, following a parabolic profile. Obstacles are present in the lower regions. Transient measurements of the moving granular surfaces were taken with a consumer-grade RGB-D sensor, providing transient 2D elevation fields around the obstacles. Three experimental configurations were tested, with semispheres of different diameters and a square dike obstacle. The experimental results are very consistent and repeatable. The quantitative, transient and two-dimensional data for all three experiments constitute excellent benchmarking tests for computational models, such as the one presented in a companion paper.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{Juez2014142,

title = {2D dry granular free-surface transient flow over complex topography with obstacles. Part II: Numerical predictions of fluid structures and benchmarking},

author = {C. Juez and D. Caviedes-Voullième and J. Murillo and P. García-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S0098300414002222},

doi = {10.1016/j.cageo.2014.09.010},

issn = {0098-3004},

year  = {2014},

date = {2014-01-01},

journal = {Computers & Geosciences},

volume = {73},

number = {0},

pages = {142 - 163},

abstract = {Abstract Dense granular flows are present in geophysics and in several industrial processes, which has lead to an increasing interest for the knowledge and understanding of the physics which govern their propagation. For this reason, a wide range of laboratory experiments on gravity-driven flows have been carried out during the last two decades. The present work is focused on geomorphological processes and, following previous work, a series of laboratory studies which constitute a further step in mimicking natural phenomena are described and simulated. Three situations are considered with some common properties: a two-dimensional configuration, variable slope of the topography and the presence of obstacles. The setup and measurement technique employed during the development of these experiments are deeply explained in the companion work. The first experiment is based on a single obstacle, the second one is performed against multiple obstacles and the third one studies the influence of a dike on which overtopping occurs. Due to the impact of the flow against the obstacles, fast moving shocks appear, and a variety of secondary waves emerge. In order to delve into the physics of these types of phenomena, a shock-capturing numerical scheme is used to simulate the cases. The suitability of the mathematical models employed in this work has been previously validated. Comparisons between computed and experimental data are presented for the three cases. The computed results show that the numerical tool is able to predict faithfully the overall behavior of this type of complex dense granular flow.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{Lacasta20141,

title = {An optimized GPU implementation of a 2D free surface simulation model on unstructured meshes},

author = {A. Lacasta and M. Morales-Hernández and J. Murillo and P. García-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S0965997814001331},

doi = {10.1016/j.advengsoft.2014.08.007},

issn = {0965-9978},

year  = {2014},

date = {2014-01-01},

journal = {Advances in Engineering Software},

volume = {78},

number = {0},

pages = {1 - 15},

abstract = {Abstract This work is related with the implementation of a finite volume method to solve the 2D Shallow Water Equations on Graphic Processing Units (GPU). The strategy is fully oriented to work efficiently with unstructured meshes which are widely used in many fields of Engineering. Due to the design of the GPU cards, structured meshes are better suited to work with than unstructured meshes. In order to overcome this situation, some strategies are proposed and analyzed in terms of computational gain, by means of introducing certain ordering on the unstructured meshes. The necessity of performing the simulations using unstructured instead of structured meshes is also justified by means of some test cases with analytical solution.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{CaviedesVoullième2014206,

title = {Reconstruction of 2D river beds by appropriate interpolation of 1D cross-sectional information for flood simulation},

author = {Daniel Caviedes-Voullième and Mario Morales-Hernández and Ibai López-Marijuan and Pilar García-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S1364815214002308},

doi = {10.1016/j.envsoft.2014.07.016},

issn = {1364-8152},

year  = {2014},

date = {2014-01-01},

journal = {Environmental Modelling & Software},

volume = {61},

number = {0},

pages = {206 - 228},

abstract = {Abstract The 2D numerical simulation of river flow requires a large amount of topographic data to build an accurate Digital Terrain Model which must cover the main river channel and the area likely to be flooded. DTMs for large floodplains are often generated by LiDAR flights. However, it is often impossible to obtain LiDAR data of permanently inundated river beds. These areas are often surveyed and discrete cross-sections of the river channel are obtained. This work presents an algorithm to generate the missing information for the areas between cross-sections. The algorithm allows to generate a river bed which preserves important morphological features such as meanders and thalweg trajectory. Two benchmark cases are studied: a synthetic river-floodplain system and a real case application on a reach of the Ebro river in Spain. The cases are analyzed from a geometry and hydrodynamics perspective by performing 2D simulations with good results.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{Juez201493,

title = {A 2D weakly-coupled and efficient numerical model for transient shallow flow and movable bed},

author = {C. Juez and J. Murillo and P. García-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S0309170814001080},

doi = {10.1016/j.advwatres.2014.05.014},

issn = {0309-1708},

year  = {2014},

date = {2014-01-01},

journal = {Advances in Water Resources},

volume = {71},

number = {0},

pages = {93 - 109},

abstract = {Abstract Recent advances in the simulation of free surface flows over mobile bed have shown that accurate and stable results in realistic problems can be provided if an appropriate coupling between the shallow water equations (SWE) and the Exner equation is performed. This coupling can be done if using a suitable Jacobian matrix. As a result, faithful numerical predictions are available for a wide range of flow conditions and empirical bed load discharge formulations, allowing to investigate the best option in each case study, which is mandatory in these type of environmental problems. When coupling the equations, the SWE are considered but including an extra conservation law for the sediment dynamics. In this way the computational cost may become unrealistic in situations where the application of the SWE over rigid bed can be used involving large time and space scales without giving up to the adequate level of mesh refinement. Therefore, for restoring the numerical efficiency, the coupling technique is simplified, not decreasing the number of waves involved in the Riemann problem but simplifying their definitions. The effects of the approximations made are tested against experimental data which include transient problems over erodible bed. The simplified model is formulated under a general framework able to insert any desirable discharge solid load formula.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{FLD:FLD3892,

title = {Numerical modelling of bridges in 2D shallow water flow simulations},

author = {H. Ratia and J. Murillo and P. García-Navarro},

url = {http://dx.doi.org/10.1002/fld.3892},

doi = {10.1002/fld.3892},

issn = {1097-0363},

year  = {2014},

date = {2014-01-01},

journal = {International Journal for Numerical Methods in Fluids},

volume = {75},

number = {4},

pages = {250–272},

abstract = {Two representations of bridges are presented to fit into simulations based on the 2D SWEs. They include a novel and simple empirical formulation for total head loss caused by bridges which covers all flow regimes: free water surface, partially submerged flow, and fully submerged flow. The performances of the different representations of bridges are assessed using experimental results of steady and unsteady flows, presented for the first time in this work. A case study of a real flood is presented giving insight into the application of the representation techniques to real engineering problems. The proposed solution of a head loss representation together with the novel formulation is shown to be accurate for both steady and unsteady flows, and to perform efficiently for solving the flood case. Copyright © 2014 John Wiley & Sons, Ltd.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{Burguete201491,

title = {SURCOS: A software tool to simulate irrigation and fertigation in isolated furrows and furrow networks},

author = {J. Burguete and A. Lacasta and P. García-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S0168169914000398},

doi = {10.1016/j.compag.2014.02.004},

issn = {0168-1699},

year  = {2014},

date = {2014-01-01},

journal = {Computers and Electronics in Agriculture},

volume = {103},

number = {0},

pages = {91 - 103},

abstract = {Abstract A software tool useful for the numerical computation of surface irrigation and fertigation in furrows and furrow networks was developed. The model solves the complete one-dimensional St-Venant equations together with the transport equation of a passive solute. The flow equations and the solute advection are solved with a high resolution TVD explicit Eulerian scheme. The solute dispersion is solved with a centered implicit Eulerian scheme to avoid further restriction in the allowable time step. The computational speed of the model is high in isolated furrows. In cases of large furrow networks over extended irrigation times the model is slower but affordable computational speed is achieved. The computational model has been designed to be robust, intuitive and able to supply useful visual results. Both the executable and the source code, as well as the examples presented can be downloaded, edited and distributed under a BSD type license.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{Lacasta2014978,

title = {Simulation of PID Control Applied to Irrigation Channels},

author = {A. Lacasta and M. Morales-Hernández and P. Brufau and P. García-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S1877705814001118},

doi = {10.1016/j.proeng.2014.02.109},

issn = {1877-7058},

year  = {2014},

date = {2014-01-01},

journal = {Procedia Engineering},

volume = {70},

number = {0},

pages = {978 - 987},

abstract = {Abstract Open-channel flow usually includes many hydraulic elements to help with the regulation of water supply in terms of automatic control. On the other hand, the one-dimensional Shallow Water Equations (SWE) are widely used to model and predict the flow dynamics in this kind of configurations. In this work, the unsteady SWE are used to model the water motion and they are solved using a finite volume upwind scheme able to cope with all flow regimes. Furthermore, the regulation of hydraulic structures at channels is frequently based on the PID controller. In this work, the implementation and coupling of the channel flow simulation with hydraulic elements and PID regulation is performed.},

note = {12th International Conference on Computing and Control for the Water Industry, CCWI2013},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{FernándezPato2014641,

title = {A Pipe Network Simulation Model with Dynamic Transition between Free Surface and Pressurized Flow},

author = {J. Fernández-Pato and P. García-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S1877705814000721},

doi = {10.1016/j.proeng.2014.02.070},

issn = {1877-7058},

year  = {2014},

date = {2014-01-01},

journal = {Procedia Engineering},

volume = {70},

number = {0},

pages = {641 - 650},

abstract = {Abstract Water flow numerical simulation in urban sewer systems is a topic that combines surface flows and pressurized flows in steady and transient situations. A numerical simulation model is developed in this work, capable of solving pipe networks mainly unpressur- ized, with isolated peaks of pressurization. For this purpose, a reformulation of the mathematical model through the Preissmann slot method is proposed. The numerical model is based on the first order Roe's scheme, in the frame of finite volume methods. The validation has been done by means of several cases with analytic solutions or empirical laboratory data.},

note = {12th International Conference on Computing and Control for the Water Industry, CCWI2013},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{FernándezPatoDWES,

title = {Finite volume simulation of unsteady water pipe flow},

author = {J. Fernández-Pato and P. García-Navarro},

url = {http://www.drink-water-eng-sci.net/7/83/2014/dwes-7-83-2014.html},

doi = {doi:10.5194/dwes-7-83-2014},

year  = {2014},

date = {2014-01-01},

journal = {Drinking Water Engineering and Science},

volume = {7},

number = {0},

pages = {83 - 92},

abstract = {The most commonly used hydraulic network models used in the drinking water community exclusively consider fully filled pipes. However, water flow numerical simulation in urban pipe systems may require to model transitions between surface flow and pressurized flow in steady and transient situations. The governing equations for both flow types are different and this must be taken into account in order to get a complete numerical model for solving dynamically transients. In this work, a numerical simulation tool is developed, capable of simulating pipe networks mainly unpressurized, with isolated points of pressurization. For this purpose, the mathematical model is reformulated by means of the Preissmann slot method. This technique provides a reasonable estimation of the water pressure in cases of pressurization. The numerical model is based on the first order Roe's scheme, in the frame of finite volume methods. The novelty of the method is that it is adapted to abrupt transient situations, with subcritical and supercritical flows. The validation has been done by means of several cases with analytic solutions or empirical laboratory data. It has also been applied to some more complex and realistic cases, like junctions or pipe networks.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{Morales,

title = {A 2D extension of a Large Time Step explicit scheme (CFL > 1) for unsteady problems with wet/dry boundaries},

author = {M. Morales-Hernandez and M. E. Hubbard and P. Garcia-Navarro},

doi = {10.1016/j.jcp.2014.01.019},

year  = {2014},

date = {2014-01-01},

journal = {JOURNAL OF COMPUTATIONAL PHYSICS},

volume = {263},

pages = {303-327},

abstract = {A 2D Large Time Step (LTS) explicit scheme on structured grids is presented in this work. It is first detailed and analysed for the 2D linear advection equation and then applied to the 2D shallow water equations. The dimensional splitting technique allows us to extend the ideas developed in the 1D case related to source terms, boundary conditions and the reduction of the time step in the presence of large discontinuities. The boundary conditions treatment as well as the wet/dry fronts in the case of the 2D shallow water equations require extra effort. The proposed scheme is tested on linear and non-linear equations and systems, with and without source terms. The numerical results are compared with those of the conventional scheme as well as with analytical solutions and experimental data.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{Murillo,

title = {Accurate numerical modeling of 1D flow in channels with arbitrary shape. Application of the energy balanced property},

author = {J. Murillo and P. Garcia-Navarro},

doi = {10.1016/j.jcp.2013.12.040},

year  = {2014},

date = {2014-01-01},

journal = {JOURNAL OF COMPUTATIONAL PHYSICS},

volume = {260},

pages = {222-248},

abstract = {This work focuses on the numerical treatment of 1D flow in channels with arbitrary shape using energy balanced arguments. The system of equations is defined using the mass and momentum conservation equations, allowing the resolution of hydraulic jumps where energy conservation arguments are not valid. When necessary, conservation of mechanical energy takes part actively in the numerical scheme when evaluating the source terms. The numerical scheme is based on an augmented Roe solver that involves the presence of source terms by adding an extra stationary wave. The characteristics of the numerical scheme include the energy balanced property, and being only first order accurate in time and space, leads to exact numerical solutions for steady solutions with independence of the grid refinement in channels with general geometries. Riemann problems considered here involve non-prismatic channels, bed variations and the resonance regime, including the limiting situation when the Riemann data belong to the resonance hypersurface. Numerical results point out that the finite volume numerical scheme with nonconservative terms presented here, converges to the exact solution. The well balanced property is ensured, as it is a particular case of the energy balanced property in cases of quiescent equilibrium.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{CaviedesVoullième2014177,

title = {2D dry granular free-surface flow over complex topography with obstacles. Part I: experimental study using a consumer-grade RGB-D sensor},

author = {Daniel Caviedes-Voullième and Carmelo Juez and Javier Murillo and Pilar García-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S0098300414002210},

doi = {10.1016/j.cageo.2014.09.009},

issn = {0098-3004},

year  = {2014},

date = {2014-01-01},

journal = {Computers & Geosciences},

volume = {73},

number = {0},

pages = {177 - 197},

abstract = {Abstract Avalanches, debris flows and other types of gravity-driven granular flows are a common hazard in mountainous regions. These regions often have human settlements in the lower parts of valleys, with human structures dangerously exposed to the destructive effects of these geophysical flows. Therefore a scientific effort has been made to understand, model and simulate geophysical granular flows. In order for computer models and simulations to be of predictive value they need to be validated under controlled, yet nature-like conditions. This work presents an experimental study of granular flow over a simplified mountain slope and valley topography. The experimental facility has a rough bed with very high slope at the upstream end and adverse slope on the downstream end, following a parabolic profile. Obstacles are present in the lower regions. Transient measurements of the moving granular surfaces were taken with a consumer-grade RGB-D sensor, providing transient 2D elevation fields around the obstacles. Three experimental configurations were tested, with semispheres of different diameters and a square dike obstacle. The experimental results are very consistent and repeatable. The quantitative, transient and two-dimensional data for all three experiments constitute excellent benchmarking tests for computational models, such as the one presented in a companion paper.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{CaviedesVoullième2014206,

title = {Reconstruction of 2D river beds by appropriate interpolation of 1D cross-sectional information for flood simulation},

author = {Daniel Caviedes-Voullième and Mario Morales-Hernández and Ibai López-Marijuan and Pilar García-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S1364815214002308},

doi = {10.1016/j.envsoft.2014.07.016},

issn = {1364-8152},

year  = {2014},

date = {2014-01-01},

journal = {Environmental Modelling & Software},

volume = {61},

number = {0},

pages = {206 - 228},

abstract = {Abstract The 2D numerical simulation of river flow requires a large amount of topographic data to build an accurate Digital Terrain Model which must cover the main river channel and the area likely to be flooded. DTMs for large floodplains are often generated by LiDAR flights. However, it is often impossible to obtain LiDAR data of permanently inundated river beds. These areas are often surveyed and discrete cross-sections of the river channel are obtained. This work presents an algorithm to generate the missing information for the areas between cross-sections. The algorithm allows to generate a river bed which preserves important morphological features such as meanders and thalweg trajectory. Two benchmark cases are studied: a synthetic river-floodplain system and a real case application on a reach of the Ebro river in Spain. The cases are analyzed from a geometry and hydrodynamics perspective by performing 2D simulations with good results.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{FernándezPato2014641,

title = {A Pipe Network Simulation Model with Dynamic Transition between Free Surface and Pressurized Flow},

author = {J. Fernández-Pato and P. García-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S1877705814000721},

doi = {10.1016/j.proeng.2014.02.070},

issn = {1877-7058},

year  = {2014},

date = {2014-01-01},

journal = {Procedia Engineering},

volume = {70},

number = {0},

pages = {641 - 650},

abstract = {Abstract Water flow numerical simulation in urban sewer systems is a topic that combines surface flows and pressurized flows in steady and transient situations. A numerical simulation model is developed in this work, capable of solving pipe networks mainly unpressur- ized, with isolated peaks of pressurization. For this purpose, a reformulation of the mathematical model through the Preissmann slot method is proposed. The numerical model is based on the first order Roe's scheme, in the frame of finite volume methods. The validation has been done by means of several cases with analytic solutions or empirical laboratory data.},

note = {12th International Conference on Computing and Control for the Water Industry, CCWI2013},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{FernándezPatoDWES,

title = {Finite volume simulation of unsteady water pipe flow},

author = {J. Fernández-Pato and P. García-Navarro},

url = {http://www.drink-water-eng-sci.net/7/83/2014/dwes-7-83-2014.html},

doi = {doi:10.5194/dwes-7-83-2014},

year  = {2014},

date = {2014-01-01},

journal = {Drinking Water Engineering and Science},

volume = {7},

number = {0},

pages = {83 - 92},

abstract = {The most commonly used hydraulic network models used in the drinking water community exclusively consider fully filled pipes. However, water flow numerical simulation in urban pipe systems may require to model transitions between surface flow and pressurized flow in steady and transient situations. The governing equations for both flow types are different and this must be taken into account in order to get a complete numerical model for solving dynamically transients. In this work, a numerical simulation tool is developed, capable of simulating pipe networks mainly unpressurized, with isolated points of pressurization. For this purpose, the mathematical model is reformulated by means of the Preissmann slot method. This technique provides a reasonable estimation of the water pressure in cases of pressurization. The numerical model is based on the first order Roe's scheme, in the frame of finite volume methods. The novelty of the method is that it is adapted to abrupt transient situations, with subcritical and supercritical flows. The validation has been done by means of several cases with analytic solutions or empirical laboratory data. It has also been applied to some more complex and realistic cases, like junctions or pipe networks.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{ISI:000314801500010,

title = {Energy balance numerical schemes for shallow water equations with 

 discontinuous topography},

author = {J. Murillo and P. Garcia-Navarro},

doi = {10.1016/j.jcp.2012.11.003},

year  = {2013},

date = {2013-03-01},

journal = {JOURNAL OF COMPUTATIONAL PHYSICS},

volume = {236},

pages = {119-142},

abstract = {The well-balanced property that ensures quiescent equilibrium when 

 solving the shallow-water equations with varying topography is extended 

 in this work to ensure numerically a constant level of energy in steady 

 cases with velocity when necessary. This is done in the context of 

 augmented solvers that consider in their definition the presence of a 

 discontinuous bed. In order to guarantee a constant energy state a 

 proper integral approach of the bed source term is presented. This 

 approach is systematically assessed via a series of steady test cases 

 and Riemann problems including the resonance regime. (C) 2012 Elsevier 

 Inc. All rights reserved.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{ISI:000314801500010b,

title = {Energy balance numerical schemes for shallow water equations with 

 discontinuous topography},

author = {J. Murillo and P. Garcia-Navarro},

doi = {10.1016/j.jcp.2012.11.003},

year  = {2013},

date = {2013-03-01},

journal = {JOURNAL OF COMPUTATIONAL PHYSICS},

volume = {236},

pages = {119-142},

abstract = {The well-balanced property that ensures quiescent equilibrium when 

 solving the shallow-water equations with varying topography is extended 

 in this work to ensure numerically a constant level of energy in steady 

 cases with velocity when necessary. This is done in the context of 

 augmented solvers that consider in their definition the presence of a 

 discontinuous bed. In order to guarantee a constant energy state a 

 proper integral approach of the bed source term is presented. This 

 approach is systematically assessed via a series of steady test cases 

 and Riemann problems including the resonance regime. (C) 2012 Elsevier 

 Inc. All rights reserved.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{ISI:000314801500010c,

title = {Energy balance numerical schemes for shallow water equations with 

 discontinuous topography},

author = {J. Murillo and P. Garcia-Navarro},

doi = {10.1016/j.jcp.2012.11.003},

year  = {2013},

date = {2013-03-01},

journal = {JOURNAL OF COMPUTATIONAL PHYSICS},

volume = {236},

pages = {119-142},

abstract = {The well-balanced property that ensures quiescent equilibrium when 

 solving the shallow-water equations with varying topography is extended 

 in this work to ensure numerically a constant level of energy in steady 

 cases with velocity when necessary. This is done in the context of 

 augmented solvers that consider in their definition the presence of a 

 discontinuous bed. In order to guarantee a constant energy state a 

 proper integral approach of the bed source term is presented. This 

 approach is systematically assessed via a series of steady test cases 

 and Riemann problems including the resonance regime. (C) 2012 Elsevier 

 Inc. All rights reserved.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{ISI:000315008300007,

title = {Verification, conservation, stability and efficiency of a finite volume 

 method for the 1D Richards equation},

author = {D. Caviedes-Voullieme and P. Garcia-Navarro and J. Murillo},

doi = {10.1016/j.jhydrol.2012.12.008},

year  = {2013},

date = {2013-02-01},

journal = {JOURNAL OF HYDROLOGY},

volume = {480},

pages = {69-84},

abstract = {Richards equation is a non-linear partial differential equation that 

 describes flow in porous media. To solve it, numerical methods that need 

 to be conservative, stable, accurate and efficient, are required. Finite 

 volume methods have not been reported as widely as finite difference or 

 finite element methods. This work is focused on providing numerical 

 results and analysis of several finite volume schemes in one dimension. 

 The formulation of time discretization schemes is revisited, a stability 

 analysis is performed, and analytical and experimental benchmakrs are 

 provided. (C) 2012 Elsevier B.V. All rights reserved.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{ISI:000315008300007b,

title = {Verification, conservation, stability and efficiency of a finite volume 

 method for the 1D Richards equation},

author = {D. Caviedes-Voullieme and P. Garcia-Navarro and J. Murillo},

doi = {10.1016/j.jhydrol.2012.12.008},

year  = {2013},

date = {2013-02-01},

journal = {JOURNAL OF HYDROLOGY},

volume = {480},

pages = {69-84},

abstract = {Richards equation is a non-linear partial differential equation that 

 describes flow in porous media. To solve it, numerical methods that need 

 to be conservative, stable, accurate and efficient, are required. Finite 

 volume methods have not been reported as widely as finite difference or 

 finite element methods. This work is focused on providing numerical 

 results and analysis of several finite volume schemes in one dimension. 

 The formulation of time discretization schemes is revisited, a stability 

 analysis is performed, and analytical and experimental benchmakrs are 

 provided. (C) 2012 Elsevier B.V. All rights reserved.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{ISI:000315008300007c,

title = {Verification, conservation, stability and efficiency of a finite volume 

 method for the 1D Richards equation},

author = {D. Caviedes-Voullieme and P. Garcia-Navarro and J. Murillo},

doi = {10.1016/j.jhydrol.2012.12.008},

year  = {2013},

date = {2013-02-01},

journal = {JOURNAL OF HYDROLOGY},

volume = {480},

pages = {69-84},

abstract = {Richards equation is a non-linear partial differential equation that 

 describes flow in porous media. To solve it, numerical methods that need 

 to be conservative, stable, accurate and efficient, are required. Finite 

 volume methods have not been reported as widely as finite difference or 

 finite element methods. This work is focused on providing numerical 

 results and analysis of several finite volume schemes in one dimension. 

 The formulation of time discretization schemes is revisited, a stability 

 analysis is performed, and analytical and experimental benchmakrs are 

 provided. (C) 2012 Elsevier B.V. All rights reserved.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{Juez2013166,

title = {2D simulation of granular flow over irregular steep slopes using global and local coordinates},

author = {C. Juez and J. Murillo and P. García-Navarro},

url = {http://www.sciencedirect.com/science/article/pii/S0021999113005275},

doi = {10.1016/j.jcp.2013.08.002},

issn = {0021-9991},

year  = {2013},

date = {2013-01-01},

journal = {Journal of Computational Physics},

volume = {255},

number = {0},

pages = {166 - 204},

abstract = {Abstract In this work approximate augmented Riemann solvers are formulated providing appropriate numerical schemes for mathematical models of granular flow on irregular steep slopes. Fluxes and source terms are discretized to ensure steady state configurations including correct modeling of start/stop flow conditions, both in a global and a local system of coordinates. The weak solutions presented involve the effect of bed slope in pressure distribution and frictional effects by means of the adequate gravity acceleration components. The numerical solvers proposed are first tested against 1D cases with exact solution and then their results are compared with experimental data in order to check the suitability of the mathematical models described in this work. Comparisons between results provided when using global and local system of coordinates are presented. The obtained results point out that both the global and the local system of coordinates can be used to predict faithfully the overall behavior of the phenomena considered in this work.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{Lacasta2013225,

title = {Preprocess static subdomain decomposition in practical cases of 2D unsteady hydraulic simulation},

author = {A. Lacasta and P. García-Navarro and J. Burguete and J. Murillo},

url = {http://www.sciencedirect.com/science/article/pii/S0045793012000977},

doi = {10.1016/j.compfluid.2012.03.010},

issn = {0045-7930},

year  = {2013},

date = {2013-01-01},

journal = {Computers & Fluids},

volume = {80},

number = {0},

pages = {225 - 232},

abstract = {Explicit finite volume methods are frequently used and widely accepted in hydraulic models based on the shallow water approximation. The main drawback of the approach is the time step size limit imposed by the Courant–Friedrichs–Lewy numerical stability constraint. This leads to excessively long computational times in large scale cases of practical interest. At the same time, the accuracy of the numerical results is associated to the use of fine computational meshes able to achieve enough spatial resolution. Taking into account that hydraulic modelers do not have access, in general, to large computational facilities, suitable and useful parallelization techniques are required. Furthermore, if high performance computing facilities are used, it is usually necessary to provide an estimation of the requirements of computational load to cover the length of the simulation. In this work the suitability of a preprocess static subdomain decomposition is explored and presented as a promising strategy to improve the efficiency of 2D unsteady shallow water computational models over dry bed in medium scale computational facilities and, at the same time, is useful to provide a preprocess computational time estimation if large scale computational facilities are going to be used.},

note = {Selected contributions of the 23rd International Conference on Parallel Fluid Dynamics ParCFD2011},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{WRCR:WRCR20062,

title = {The formulation of internal boundary conditions in unsteady 2-D shallow water flows: Application to flood regulation},

author = {M. Morales-Hernández and J. Murillo and P. García-Navarro},

url = {10.1002/wrcr.20062},

doi = {10.1002/wrcr.20062},

issn = {1944-7973},

year  = {2013},

date = {2013-01-01},

journal = {Water Resources Research},

volume = {49},

number = {1},

pages = {471-487},

abstract = {This work presents a two-dimensional hydraulic model that includes gates as internal structures. The flow is modeled using the two-dimensional shallow water equations and the gates are formulated as internal boundary conditions to provide a simulation tool for water flood management. When open channel flow in a river passes through a gate, the shallow water equations are no longer valid and energy conservation laws are required. The change in the set of equations is avoided by modeling gates as a spatial discontinuity or internal boundary condition, providing an alternative algorithm to the one used in the rest of the flooded computational domain. In the first part of this work, the requirements of an adequate discretization for gate modeling are provided in the context of a finite volume numerical scheme able to handle all kind of flow regimes over complex bed topography. In the second part of this work, the formulation of the internal boundary conditions is verified by means of a test case with exact solution. A benchmark test case is then proposed as a synthetic river reach with lateral storage areas controlled by gates. Dimensional analysis is used to establish the regulation parameters influencing the attenuation of the outlet peak discharge. It is also shown that the peak outflow discharge can be reduced by coupling the present simulator with a proportional-integral-derivative regulation algorithm. Finally, a river reach of the Ebro River is simulated with a real flooding scenario.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{Juez1,

title = {Numerical assessment of bed-load discharge formulations for transient flow in 1D and 2D situations},

author = {C. Juez and J. Murillo and P. Garcia-Navarro},

doi = {10.2166/hydro.2013.153},

year  = {2013},

date = {2013-01-01},

journal = {JOURNAL OF HYDROINFORMATICS},

volume = {15},

number = {4},

pages = {1234-1257},

abstract = {Two-dimensional (2D) transient flow over an erodible bed can be modelled using shallow-water equations and the Exner equation to describe the morphological evolution of the bed. Considering the fact that well-proven capacity formulae are based on one-dimensional (1D) experimental steady flows, the assessment of these empirical relations under unsteady 1D and 2D situations is important. In order to ensure the reliability of the numerical experimentation, the formulation has to be general enough to allow the use of different empirical laws. Moreover, the numerical scheme must handle correctly the coupling between the 2D shallow-water equations and the Exner equation under any condition. In this work, a finite-volume numerical scheme that includes these two main features will be exploited here in 1D and 2D laboratory test cases. The relative performances of Meyer-Peter and Ashida and Michiue, Engelund and Fredsoe, Fernandez Luque and Van Beek, Parker, Smart, Nielsen, Wong and Camenen and Larson formulations are analysed in terms of the root mean square error. A new discretization of the Smart formula is provided, leading to promising predictions of the erosion/deposition rates. The results arising from this work are useful to justify the use of an empirical sediment bed-load discharge formula among the ones studied, regardless of the hydrodynamic situation.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{MoralesHernández201326,

title = {A conservative strategy to couple 1D and 2D models for shallow water flow simulation},

author = {M. Morales-Hernandez and P. Garcia-Navarro and J. Burguete and P. Brufau},

doi = {10.1016/j.compfluid.2013.04.001},

year  = {2013},

date = {2013-01-01},

journal = {Computers & Fluids},

volume = {81},

number = {0},

pages = {26 - 44},

abstract = {Abstract A 1D-2D coupled numerical model is presented in this work. 1D and 2D models are formulated using a conservative upwind cell-centred finite volume scheme. The discretization is based on cross-sections for the 1D model and with triangular unstructured grid for the 2D model. The resulting element of discretization for the coupled model is analysed and two different coupling techniques based on mass conservation and mass and momentum conservation respectively are explored, considering both frontal and lateral configurations. The interaction with the boundaries in each model is highlighted and the necessity of using the appropriate strategy according to the flow regime is also justified. The coupled model is tested through academic test cases where the numerical results are compared with a fully 2D model as well as with experimental measurements in steady and unsteady scenarios. It is also applied to a real world configuration, where the flood wave propagation in the river bed is simulated by means of a 1D model and the inundation of the riverside is dealt with a 2D model. The computational gain is also analysed.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{MoralesHernández201326,

title = {A conservative strategy to couple 1D and 2D models for shallow water flow simulation},

author = {M. Morales-Hernandez and P. Garcia-Navarro and J. Burguete and P. Brufau},

doi = {10.1016/j.compfluid.2013.04.001},

year  = {2013},

date = {2013-01-01},

journal = {Computers & Fluids},

volume = {81},

number = {0},

pages = {26 - 44},

abstract = {Abstract A 1D-2D coupled numerical model is presented in this work. 1D and 2D models are formulated using a conservative upwind cell-centred finite volume scheme. The discretization is based on cross-sections for the 1D model and with triangular unstructured grid for the 2D model. The resulting element of discretization for the coupled model is analysed and two different coupling techniques based on mass conservation and mass and momentum conservation respectively are explored, considering both frontal and lateral configurations. The interaction with the boundaries in each model is highlighted and the necessity of using the appropriate strategy according to the flow regime is also justified. The coupled model is tested through academic test cases where the numerical results are compared with a fully 2D model as well as with experimental measurements in steady and unsteady scenarios. It is also applied to a real world configuration, where the flood wave propagation in the river bed is simulated by means of a 1D model and the inundation of the riverside is dealt with a 2D model. The computational gain is also analysed.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{ISI:000312711400004,

title = {Diffusive-Wave Based Hydrologic-Hydraulic Model with Sediment Transport. 

 I: Model Development},

author = {D. Lopez-Barrera and P. Garcia-Navarro and P. Brufau and J. Burguete},

doi = {10.1061/(ASCE)HE.1943-5584.0000552},

year  = {2012},

date = {2012-10-01},

journal = {JOURNAL OF HYDROLOGIC ENGINEERING},

volume = {17},

number = {10},

pages = {1093-1104},

abstract = {In this paper, a distributed numerical model is proposed on the basis of 

 (1) a hydrologic model for the water exchange laws; (2) a surface runoff 

 model based on a hybrid two-dimensional (2D) diffusive/kinematic wave 

 approximation able to calculate flow over all kinds of wet/dry ground 

 surface slopes; (3) a groundwater flow model based on the 2D Darcy law 

 for both saturated and partly saturated zones; (4) a 2D hillslope 

 erosion model for the sediment transport; and (5) a explicit 

 finite-volume discretization with specific schemes according to the 

 characteristics of the flow equations, upwind for the hyperbolic 

 equations and centered for the parabolic equations. The resulting model 

 offers a variable time step, ensuring numerical stability with the time 

 step size sensitive to the grid cell size in the diffusive wave case and 

 an entropy correction of the upwind fluxes to ensure conservative 

 solutions near local maxima in the slopes controlling the water 

 movement. The validation and practical application of the model is 

 presented in Paper II, in which the potential usefulness of the proposed 

 model is demonstrated. DOI: 10.1061/(ASCE)HE.1943-5584.0000552. (C) 2012 

 American Society of Civil Engineers.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{ISI:000312711400005,

title = {Diffusive-Wave Based Hydrologic-Hydraulic Model with Sediment Transport. 

 II: Validation and Practical Application},

author = {D. Lopez-Barrera and P. Garcia-Navarro and P. Brufau and J. Burguete},

doi = {10.1061/(ASCE)HE.1943-5584.0000551},

year  = {2012},

date = {2012-10-01},

journal = {JOURNAL OF HYDROLOGIC ENGINEERING},

volume = {17},

number = {10},

pages = {1105-1122},

abstract = {The development of a distributed two-dimensional (2D) 

 hydrologic-hydraulic simulation model was presented in Paper I. The 

 simulation model combined overland flow (kinematic/diffusive wave 

 models), hillslope sediment transport, and groundwater flow apart from 

 the water exchange mechanisms between zones. Particular attention was 

 paid to the upwind discretization of the surface flow equations. In this 

 paper, the proposed model is validated by using four test cases with 

 exact solutions, one academic test case, and two laboratory test cases. 

 The model adequately reproduced front advance over dry beds of any slope 

 and water table evolution in simple cases. As practical application of 

 the model, the simulation of real events in two experimental basins is 

 also presented. The work is focused on the influence of the choice of 

 the empirical parameters on the model results concerning solid and 

 liquid discharges. Also, because of the lack of information referring 

 the boundary and initial conditions of the groundwater flow in real 

 basins, it is difficult to evaluate the accuracy of the complete model. 

 DOI: 10.1061/(ASCE)HE.1943-5584.0000551. (C) 2012 American Society of 

 Civil Engineers.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{ISI:000312711400004b,

title = {Diffusive-Wave Based Hydrologic-Hydraulic Model with Sediment Transport. 

 I: Model Development},

author = {D. Lopez-Barrera and P. Garcia-Navarro and P. Brufau and J. Burguete},

doi = {10.1061/(ASCE)HE.1943-5584.0000552},

year  = {2012},

date = {2012-10-01},

journal = {JOURNAL OF HYDROLOGIC ENGINEERING},

volume = {17},

number = {10},

pages = {1093-1104},

abstract = {In this paper, a distributed numerical model is proposed on the basis of 

 (1) a hydrologic model for the water exchange laws; (2) a surface runoff 

 model based on a hybrid two-dimensional (2D) diffusive/kinematic wave 

 approximation able to calculate flow over all kinds of wet/dry ground 

 surface slopes; (3) a groundwater flow model based on the 2D Darcy law 

 for both saturated and partly saturated zones; (4) a 2D hillslope 

 erosion model for the sediment transport; and (5) a explicit 

 finite-volume discretization with specific schemes according to the 

 characteristics of the flow equations, upwind for the hyperbolic 

 equations and centered for the parabolic equations. The resulting model 

 offers a variable time step, ensuring numerical stability with the time 

 step size sensitive to the grid cell size in the diffusive wave case and 

 an entropy correction of the upwind fluxes to ensure conservative 

 solutions near local maxima in the slopes controlling the water 

 movement. The validation and practical application of the model is 

 presented in Paper II, in which the potential usefulness of the proposed 

 model is demonstrated. DOI: 10.1061/(ASCE)HE.1943-5584.0000552. (C) 2012 

 American Society of Civil Engineers.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{ISI:000312711400005b,

title = {Diffusive-Wave Based Hydrologic-Hydraulic Model with Sediment Transport. 

 II: Validation and Practical Application},

author = {D. Lopez-Barrera and P. Garcia-Navarro and P. Brufau and J. Burguete},

doi = {10.1061/(ASCE)HE.1943-5584.0000551},

year  = {2012},

date = {2012-10-01},

journal = {JOURNAL OF HYDROLOGIC ENGINEERING},

volume = {17},

number = {10},

pages = {1105-1122},

abstract = {The development of a distributed two-dimensional (2D) 

 hydrologic-hydraulic simulation model was presented in Paper I. The 

 simulation model combined overland flow (kinematic/diffusive wave 

 models), hillslope sediment transport, and groundwater flow apart from 

 the water exchange mechanisms between zones. Particular attention was 

 paid to the upwind discretization of the surface flow equations. In this 

 paper, the proposed model is validated by using four test cases with 

 exact solutions, one academic test case, and two laboratory test cases. 

 The model adequately reproduced front advance over dry beds of any slope 

 and water table evolution in simple cases. As practical application of 

 the model, the simulation of real events in two experimental basins is 

 also presented. The work is focused on the influence of the choice of 

 the empirical parameters on the model results concerning solid and 

 liquid discharges. Also, because of the lack of information referring 

 the boundary and initial conditions of the groundwater flow in real 

 basins, it is difficult to evaluate the accuracy of the complete model. 

 DOI: 10.1061/(ASCE)HE.1943-5584.0000551. (C) 2012 American Society of 

 Civil Engineers.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{ISI:000312711400004c,

title = {Diffusive-Wave Based Hydrologic-Hydraulic Model with Sediment Transport. 

 I: Model Development},

author = {D. Lopez-Barrera and P. Garcia-Navarro and P. Brufau and J. Burguete},

doi = {10.1061/(ASCE)HE.1943-5584.0000552},

year  = {2012},

date = {2012-10-01},

journal = {JOURNAL OF HYDROLOGIC ENGINEERING},

volume = {17},

number = {10},

pages = {1093-1104},

abstract = {In this paper, a distributed numerical model is proposed on the basis of 

 (1) a hydrologic model for the water exchange laws; (2) a surface runoff 

 model based on a hybrid two-dimensional (2D) diffusive/kinematic wave 

 approximation able to calculate flow over all kinds of wet/dry ground 

 surface slopes; (3) a groundwater flow model based on the 2D Darcy law 

 for both saturated and partly saturated zones; (4) a 2D hillslope 

 erosion model for the sediment transport; and (5) a explicit 

 finite-volume discretization with specific schemes according to the 

 characteristics of the flow equations, upwind for the hyperbolic 

 equations and centered for the parabolic equations. The resulting model 

 offers a variable time step, ensuring numerical stability with the time 

 step size sensitive to the grid cell size in the diffusive wave case and 

 an entropy correction of the upwind fluxes to ensure conservative 

 solutions near local maxima in the slopes controlling the water 

 movement. The validation and practical application of the model is 

 presented in Paper II, in which the potential usefulness of the proposed 

 model is demonstrated. DOI: 10.1061/(ASCE)HE.1943-5584.0000552. (C) 2012 

 American Society of Civil Engineers.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{ISI:000312711400005c,

title = {Diffusive-Wave Based Hydrologic-Hydraulic Model with Sediment Transport. 

 II: Validation and Practical Application},

author = {D. Lopez-Barrera and P. Garcia-Navarro and P. Brufau and J. Burguete},

doi = {10.1061/(ASCE)HE.1943-5584.0000551},

year  = {2012},

date = {2012-10-01},

journal = {JOURNAL OF HYDROLOGIC ENGINEERING},

volume = {17},

number = {10},

pages = {1105-1122},

abstract = {The development of a distributed two-dimensional (2D) 

 hydrologic-hydraulic simulation model was presented in Paper I. The 

 simulation model combined overland flow (kinematic/diffusive wave 

 models), hillslope sediment transport, and groundwater flow apart from 

 the water exchange mechanisms between zones. Particular attention was 

 paid to the upwind discretization of the surface flow equations. In this 

 paper, the proposed model is validated by using four test cases with 

 exact solutions, one academic test case, and two laboratory test cases. 

 The model adequately reproduced front advance over dry beds of any slope 

 and water table evolution in simple cases. As practical application of 

 the model, the simulation of real events in two experimental basins is 

 also presented. The work is focused on the influence of the choice of 

 the empirical parameters on the model results concerning solid and 

 liquid discharges. Also, because of the lack of information referring 

 the boundary and initial conditions of the groundwater flow in real 

 basins, it is difficult to evaluate the accuracy of the complete model. 

 DOI: 10.1061/(ASCE)HE.1943-5584.0000551. (C) 2012 American Society of 

 Civil Engineers.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{ISI:000307299400015,

title = {Augmented versions of the HLL and HLLC Riemann solvers including source 

 terms in one and two dimensions for shallow flow applications},

author = {J. Murillo and P. Garcia-Navarro},

doi = {10.1016/j.jcp.2012.06.031},

year  = {2012},

date = {2012-08-01},

journal = {JOURNAL OF COMPUTATIONAL PHYSICS},

volume = {231},

number = {20},

pages = {6861-6906},

abstract = {Shallow water flows are found in a variety of engineering problems 

 always dominated by the presence of bed friction and irregular 

 bathymetry. These source terms determine completely the possible 

 evolution of the flooded area in time. It is well known that appropriate 

 numerical schemes for this type of flows must be well-balanced. 

 Well-balanced numerical schemes are based on the preservation of cases 

 of quiescent equilibrium over variable bed elevation. Commonly they are 

 formulated as an adaptation of numerical solvers defined for cases 

 without source terms. This procedure is insufficient when applied to 

 real situations. Then, it is possible to argue that appropriate 

 numerical schemes cannot arise directly from those derived from the 

 simplest homogeneous case without source terms. New solutions are 

 presented in this work by defining weak solutions that include the 

 presence of source terms. To do that, the solvers presented in this work 

 extend the number of waves in the well known HLL and HLLC solvers 

 involving a stationary jump in the solution. This is done without 

 modifying the original solution vector of conserved quantities. The 

 resulting approximate Riemann solvers include variable bed level surface 

 and friction. Solvers are systematically assessed via a series of test 

 problems with exact solutions for one and two dimensions, including 

 steady and unsteady flow configurations, variation of the flooded area 

 in time and comparisons with experimental data. The obtained results 

 point out that the new method is able to predict faithfully the overall 

 behavior of the solution and of any type of waves. (C) 2012 Elsevier 

 Inc. All rights reserved.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{ISI:000306606900017,

title = {A large time step 1D upwind explicit scheme (CFL > 1): Application to 

 shallow water equations},

author = {M. Morales-Hernandez and P. Garcia-Navarro and J. Murillo},

doi = {10.1016/j.jcp.2012.06.017},

year  = {2012},

date = {2012-08-01},

journal = {JOURNAL OF COMPUTATIONAL PHYSICS},

volume = {231},

number = {19},

pages = {6532-6557},

abstract = {It is possible to relax the Courant-Friedrichs-Lewy condition over the 

 time step when using explicit schemes. This method, proposed by Leveque, 

 provides accurate and correct solutions of non-sonic shocks. 

 Rarefactions need some adjustments which are explored in the present 

 work with scalar equation and systems of equations. The non-conservative 

 terms that appear in systems of conservation laws introduce an extra 

 difficulty in practical application. The way to deal with source terms 

 is incorporated into the proposed procedure. The boundary treatment is 

 analysed and a reflection wave technique is considered. In presence of 

 strong discontinuities or important source terms, a strategy is proposed 

 to control the stability of the method allowing the largest time step 

 possible. The performance of the above scheme is evaluated to solve the 

 homogeneous shallow water equations and the shallow water equations with 

 source terms. (C) 2012 Elsevier Inc. All rights reserved.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{ISI:000307299400015b,

title = {Augmented versions of the HLL and HLLC Riemann solvers including source 

 terms in one and two dimensions for shallow flow applications},

author = {J. Murillo and P. Garcia-Navarro},

doi = {10.1016/j.jcp.2012.06.031},

year  = {2012},

date = {2012-08-01},

journal = {JOURNAL OF COMPUTATIONAL PHYSICS},

volume = {231},

number = {20},

pages = {6861-6906},

abstract = {Shallow water flows are found in a variety of engineering problems 

 always dominated by the presence of bed friction and irregular 

 bathymetry. These source terms determine completely the possible 

 evolution of the flooded area in time. It is well known that appropriate 

 numerical schemes for this type of flows must be well-balanced. 

 Well-balanced numerical schemes are based on the preservation of cases 

 of quiescent equilibrium over variable bed elevation. Commonly they are 

 formulated as an adaptation of numerical solvers defined for cases 

 without source terms. This procedure is insufficient when applied to 

 real situations. Then, it is possible to argue that appropriate 

 numerical schemes cannot arise directly from those derived from the 

 simplest homogeneous case without source terms. New solutions are 

 presented in this work by defining weak solutions that include the 

 presence of source terms. To do that, the solvers presented in this work 

 extend the number of waves in the well known HLL and HLLC solvers 

 involving a stationary jump in the solution. This is done without 

 modifying the original solution vector of conserved quantities. The 

 resulting approximate Riemann solvers include variable bed level surface 

 and friction. Solvers are systematically assessed via a series of test 

 problems with exact solutions for one and two dimensions, including 

 steady and unsteady flow configurations, variation of the flooded area 

 in time and comparisons with experimental data. The obtained results 

 point out that the new method is able to predict faithfully the overall 

 behavior of the solution and of any type of waves. (C) 2012 Elsevier 

 Inc. All rights reserved.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{ISI:000306606900017b,

title = {A large time step 1D upwind explicit scheme (CFL > 1): Application to 

 shallow water equations},

author = {M. Morales-Hernandez and P. Garcia-Navarro and J. Murillo},

doi = {10.1016/j.jcp.2012.06.017},

year  = {2012},

date = {2012-08-01},

journal = {JOURNAL OF COMPUTATIONAL PHYSICS},

volume = {231},

number = {19},

pages = {6532-6557},

abstract = {It is possible to relax the Courant-Friedrichs-Lewy condition over the 

 time step when using explicit schemes. This method, proposed by Leveque, 

 provides accurate and correct solutions of non-sonic shocks. 

 Rarefactions need some adjustments which are explored in the present 

 work with scalar equation and systems of equations. The non-conservative 

 terms that appear in systems of conservation laws introduce an extra 

 difficulty in practical application. The way to deal with source terms 

 is incorporated into the proposed procedure. The boundary treatment is 

 analysed and a reflection wave technique is considered. In presence of 

 strong discontinuities or important source terms, a strategy is proposed 

 to control the stability of the method allowing the largest time step 

 possible. The performance of the above scheme is evaluated to solve the 

 homogeneous shallow water equations and the shallow water equations with 

 source terms. (C) 2012 Elsevier Inc. All rights reserved.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{ISI:000307299400015c,

title = {Augmented versions of the HLL and HLLC Riemann solvers including source 

 terms in one and two dimensions for shallow flow applications},

author = {J. Murillo and P. Garcia-Navarro},

doi = {10.1016/j.jcp.2012.06.031},

year  = {2012},

date = {2012-08-01},

journal = {JOURNAL OF COMPUTATIONAL PHYSICS},

volume = {231},

number = {20},

pages = {6861-6906},

abstract = {Shallow water flows are found in a variety of engineering problems 

 always dominated by the presence of bed friction and irregular 

 bathymetry. These source terms determine completely the possible 

 evolution of the flooded area in time. It is well known that appropriate 

 numerical schemes for this type of flows must be well-balanced. 

 Well-balanced numerical schemes are based on the preservation of cases 

 of quiescent equilibrium over variable bed elevation. Commonly they are 

 formulated as an adaptation of numerical solvers defined for cases 

 without source terms. This procedure is insufficient when applied to 

 real situations. Then, it is possible to argue that appropriate 

 numerical schemes cannot arise directly from those derived from the 

 simplest homogeneous case without source terms. New solutions are 

 presented in this work by defining weak solutions that include the 

 presence of source terms. To do that, the solvers presented in this work 

 extend the number of waves in the well known HLL and HLLC solvers 

 involving a stationary jump in the solution. This is done without 

 modifying the original solution vector of conserved quantities. The 

 resulting approximate Riemann solvers include variable bed level surface 

 and friction. Solvers are systematically assessed via a series of test 

 problems with exact solutions for one and two dimensions, including 

 steady and unsteady flow configurations, variation of the flooded area 

 in time and comparisons with experimental data. The obtained results 

 point out that the new method is able to predict faithfully the overall 

 behavior of the solution and of any type of waves. (C) 2012 Elsevier 

 Inc. All rights reserved.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{ISI:000306606900017c,

title = {A large time step 1D upwind explicit scheme (CFL > 1): Application to 

 shallow water equations},

author = {M. Morales-Hernandez and P. Garcia-Navarro and J. Murillo},

doi = {10.1016/j.jcp.2012.06.017},

year  = {2012},

date = {2012-08-01},

journal = {JOURNAL OF COMPUTATIONAL PHYSICS},

volume = {231},

number = {19},

pages = {6532-6557},

abstract = {It is possible to relax the Courant-Friedrichs-Lewy condition over the 

 time step when using explicit schemes. This method, proposed by Leveque, 

 provides accurate and correct solutions of non-sonic shocks. 

 Rarefactions need some adjustments which are explored in the present 

 work with scalar equation and systems of equations. The non-conservative 

 terms that appear in systems of conservation laws introduce an extra 

 difficulty in practical application. The way to deal with source terms 

 is incorporated into the proposed procedure. The boundary treatment is 

 analysed and a reflection wave technique is considered. In presence of 

 strong discontinuities or important source terms, a strategy is proposed 

 to control the stability of the method allowing the largest time step 

 possible. The performance of the above scheme is evaluated to solve the 

 homogeneous shallow water equations and the shallow water equations with 

 source terms. (C) 2012 Elsevier Inc. All rights reserved.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}