This paper describes the application of a commercially available, three-dimensional computational fluid dynamic (CFD) model to simulate the flow structure in an upland river that is prone to flooding. Simulations use a rectangular channel geometry, smooth sidewalls and a bed topography obtained from the field site that contains a subdued pool–riffle sequence. The CFD model uses the RNG κ– ε turbulence closure scheme of Yakhot and Orszag (J. Sci. Comput. 1 (1986) 1), as implemented in FLUENT 4.4.4, with a free surface. Results are shown for numerical runs simulating a 1:100 year return interval flood. Output from the numerical model is compared to a physical model experiment that uses a 1:35 scale fibreglass mould of the field study reach and measures velocity using ultrasonic Doppler velocity profiling (UDVP). Results are presented from the numerical and flume models for the water surface and streamwise velocity pattern and for the secondary flows simulated in the numerical model. A good agreement is achieved between the CFD model output and the physical model results for the downstream velocities. Results suggest that the streamwise velocity is the main influence on the flow structure at the discharge and channel configuration studied. Secondary flows are, in general, very weak being below the resolution of measurement in the physical model and less than 10% of the streamwise velocity in the numerical model. Consequently, there is no evidence for a ‘velocity dip’. It is suggested that the subdued topography or inlet morphology may inhibit the development of secondary flows that have been recorded in previous flat-bed, rectangular open channel flows. A significant corollary of these results is that the morphological evolution of the pool–riffle sequence at high discharges may be controlled primarily by the downstream distribution of velocity and sediment transport with little role for lateral sorting and sediment routing by secondary flows. This paper also raises a number of issues that may be of use in future CFD modelling of three-dimensional flow in open channels within the geomorphological community.
Read full abstract