GeoscienceAustralia
Simulating modified porous shallow-water equations
I want to change the ANUGA source codes by adding a porosity parameter into the classical shallow-water equation to simulate urban flooding. This governing equation was described in Velickovic et al.(2010): “The porosity approach leads to modified shallow-water equations including a porosity parameter to account for the area available to the flow in a given control volume”. I don't know if ANUGA can simulate this modified shallow-water equation with porosity. Do I need to configure this new parameter in config.py?
@wpdengzheng, could you give more detail about your reference. I couldn't find the paper with the info you have given.
My guess is that the modification would involve updating the flux/evolve code which is not as simple as adding an operator or forcing term. If it is just a modification to the time derivative term then the modification could be simple.
@stoiver Thank you for your reply. I have a detailed description of the porous equation here, and upload these two references. Referance [1] Isotropic porosity.pdf Referance [2] Anisotropic porosity.pdf
The porosity parameter includes isotropy and anisotropy. When this parameter is added into the governing equation, there are some changes. Isotropic porous shallow-water equation is described in Eq. (1)(2) (Reference [1]), where the conserved vector variable U, the flux F and G have one more parameter than the classical shallow-water equation. Likewise, The source term S also has two additional terms, which are described in Eq. (3)(4) (Reference [1]). I would like to use anisotropic porous shallow-water equation, which is similar to the above isotropic one, but the single porosity becomes two different porous values: edge-based porosity and cell-based porosity, which are shown in Eq.(6)(7)(8) (Reference [2]). I want to change the ANUGA source codes based on the anisotropic porous equation. In the process of running, the porosity values, as a known file, will be called in every time step for each grid cell. But I don't know if ANUGA can simulate this modified governing equation with the new parameter.
Reference [1]: Soares-Frazão S, Lhomme J, Guinot V, et al. Two-dimensional shallow-water model with porosity for urban flood modelling[J]. Journal of Hydraulic Research, 2008, 46(1): 45-64. Reference [2]: Schubert J E, Sanders B F. Building treatments for urban flood inundation models and implications for predictive skill and modeling efficiency[J]. Advances in Water Resources, 2012, 41: 49-64.
So the aim is to allow flow to pass through solid objects such as Buildings ? Can this be done using a Porous River Wall ?
