Numerical modeling of cohesive sediment transport and bed morphology in estuaries

Two major lines of investigation have been pursued in this thesis: (1) More efficient, robust and realistic numerical techniques are designed for the simulation of complex turbulent fluid flows; (2) A new algorithm and its analysis is performed in the context of multiphasic fluid flow, for a cohesive fine-grained sediment (fluid mud) transport in estuaries. Estuaries exist between marine and freshwater system where waters of different physical, chemical and biological composition meet, combine and disperse primarily due to tidal influences. In the present thesis, the behavior of cohesive sediment in estuaries is reviewed based on the existing literature. Basic theories and recent developments are introduced to describe the formation of fluid mud from a very dilute suspension and its motion down a natural river bed with complex bathymetry. The present work contributes to the numerical simulation of complex turbulent multiphasic fluid flows encountered in estuarine channels, with the aim of the better understanding of the underlying physical processes as well as predicting realistically the cohesive sediment transport and bed morphology in such a zone. The model is based on the mass preserving method by using the so-called Raviart-Thomas finite element on the unstructured mesh in the horizontal plane. In the vertical, the computational domain is divided into number of layers at predefined heights and the method uses a conventional conforming P1 finite element scheme, with the advantage that the lowermost and uppermost layers variable height allow a faithful representation of the time-varying bed and free surface, respectively. Concerning the modeling of turbulence, the research effort focuses on the turbulence two-equation k - ε closure for the vertical parameterization of eddy viscosity. More precisely, a robust up-to-date algorithm is used for this issue. The new methodology is developed with the aim to account for more general relevant effects in the closure. The proposed model offers the capability to cope with the stiffness problem introduced by the large difference between the solid phase flow time scale and the hydrodynamic one, by using a sub-cycling strategy, whereas the splitting scheme is adopted with the aim of stability and the positivity of the relevant turbulent variables. The flexibility of the model and its performance are evaluated on several free-surface flow configurations with increasing complexity : homogeneous unsteady non-uniform flows in plane open channel flows, U-shaped (193°) curved open channel flow. Concerning the cohesive sediment transport, most of the existing models in the literature assume the analogous transport characteristics with that of the coarse sediment and adopt the relevant developed sediment transport for the latter to treat the former. Moreover, these existing models do not account for the consolidation of the mud-bed. The present research effort focused on a novel methodology based on the realistic em