Numerical modelling of turbulence and sediment concentrations under breaking waves using OpenFOAM®

Abstract

This thesis presents the development of a novel numerical model capable of evaluating suspended sediment dynamics under breaking waves, and is based in the open source Computational Fluid Dynamics software, OpenFOAM®. The hydrodynamics were determined by solving the incompressible, Reynolds-Averaged Navier-Stokes equations for a two-phase fluid using the Finite Volume method, along with a Volume of Fluid scheme that modelled the interface between the air and water phases.

A new library of five turbulence models was developed to include weakly compressible effects through the introduction of density variations in the conservation equations. This library was thoroughly evaluated against existing physical data for surf zone dynamics. A skill score was applied, based on the MSE, to rank the models, with the nonlinear k−ε performing best overall, and the k−ω predicting turbulent kinetic energy most accurately. Furthermore, the numerical model was shown to predict the near-bed hydrodynamics well, through comparison with new in-house physical data obtained in the COAST laboratory.

Suspended sediment concentrations were determined using an advection-diffusion methodology, with near-bed processes modelled using a flux based approach that balances entrainment and deposition. The model was validated against existing experimental data for steady state flow conditions, as well as for regular and breaking waves. The agreement was generally good, with the results indicating that the model is capable of capturing complicated processes such as sediment plumes under plunging breakers.

The validated model was applied to investigate the properties of the sediment diffusivity, which is a vital parameter in suspended sediment dynamics. In physical experiments, sediment diffusivity is commonly estimated implicitly, based on the vertical concentration profile. In this work, this approach was applied to the numerical concentration predictions, and compared with the value directly determined within the model. The estimated value was generally acceptable providing that large horizontal concentration gradients were not present, and diffusion dominated flow advection. However, near the breaking point of waves, large errors were observed at mid-depth of the water column, which strongly correlates with a region of large flow advection relative to diffusion. Therefore, when using this estimation, caution is recommended since this approach can potentially lead to substantial discrepancies.