Feasibility of the porous zone approach to modelling vegetation in CFD

Vegetation within stormwater ponds varies seasonly and its presence affects the flow field, which in turn affects the pond’s Residence Time Distribution and its effectiveness at pollutant removal. Vegetated flows are complex and, as a result, few suitable tools exist for evaluating realistic stormwater pond designs. Recent research has suggested using a porous zone to represent vegetation within a CFD model, and this paper investigates the feasibility of this approach using ANSYS Fluent. One of the main benefits of using a porous zone is the ability to derive the relevant parameters from the known physical characteristics of stem diameter and porosity using the Ergun equation. A sensitivity analysis on the viscous resistance factor 1/α1/α and the inertial resistance factor C2C2 has been undertaken by comparing model results to data collected from an experimental vegetated channel. Best fit values of C2C2 were obtained for a range of flow conditions including emergent and submerged vegetation. Results show the CFD model to be insensitive to 1/α1/α but very sensitive to values of C2C2. For submerged vegetation, values of C2C2 derived from the Ergun equation are under-predictions of best-fit C2C2 values as only the turbulence due to the shear layer is represented. The porous zone approach does not take into account turbulence generated from stem wakes such that no meaningful predictions for emergent vegetation were obtained. C2C2 values calculated using a force balance show better agreement with best-fit C2C2 values than those derived from the Ergun equation. Manually fixing values of kk and εε within the porous zone of the model shows initial promise as a means of taking stem wakes into account.