A Multi-dimensional Two-phase Flow Modeling Framework for Sediment Transport Applications


Book Description

Studying coastal processes is essential for the sustainability of human habitat and vibrancy of coastal economy. Coastal morphological evolution is caused by a wide range of coupled cross-shore and alongshore sediment transport processes associated with short waves, infra-gravity waves, and wave-induced currents. One of the key challenges was that the major transport occurs within bottom boundary layers and it is dictated by turbulence-sediment interactions and inter-granular interactions. Therefore, this study focuses on numerical investigations of sediment transport in the bottom wave boundary layers on continental shelves and nearshore zones, with emphasis on both fine sediment (mud) and sand transports. On the continental shelves, the sea floor is often covered with fine sediments (with settling velocity no more than a few mm/s). Wave-induced resuspension has been identified as one of the major mechanisms in the offshore delivery for fine sediments. A series of turbulence-resolving simulations were carried out to study the role of sediment resuspension/deposition on the bottom sediment transport. Specifically, we focus on how the critical shear stress of erosion and the settling velocity can determine the transport modes. At a given wave intensity associated with more energetic muddy shelves, three transport modes, namely the well-mixed transport (mode I), two-layer like transport with the formation of lutocline (mode II) and laminarized transport (mode III), are obtained by varying the critical shear stress of erosion or the settling velocity. A 2D parametric map is proposed to characterize the transition between transport modes as a function of the critical shear stress and the settling velocity at a fixed wave intensity. In addition, the uncertainties due to hindered settling and particle inertia effects on the transport modes were further studied. Simulation results confirmed that the effect of particle inertia is negligible for fine sediment in typical wave condition on continental shelves. On the other hand, the hindered settling with low gelling concentration can play a key role in sustaining a large amount of suspended sediments and results in the laminarized transport (mode III). Low gelling concentrations can also trigger the occurrence of gelling ignition, a state in which the erosion rate always exceeds the deposition rate. A sufficient condition for the occurrence of gelling ignition is hypothesized for a range of wave intensities as a function of sediment/floc properties and erodibility parameters. In the more energetic nearshore zones, the sea floor is often covered with sand (with settling velocity exceeds 1 cm/s). Based on the open-source CFD toolbox OpenFOAM, a multi-dimensional Eulerian two-phase modeling framework is developed for sediment transport applications. With closures of particle stresses and fluid-particle interactions, the model is able to resolve full sediment transport profiles without conventional bedload/suspended load assumptions. The turbulence-averaged model is based on a modified k-epsilon closure for the carrier flow turbulence and it was used to study momentary bed failure under sheet flow conditions. Model results revealed that the momentary bed failure and the resulting large transport rate were associated with a large erosion depth, which was triggered by the combination of large bed shear stresses and large horizontal pressure gradients. In order to better resolve turbulence-sediment interactions, the modeling framework was also extended with a 3D turbulence-resolving capability, where most of the turbulence-sediment interactions are directly resolved. The model is validated against a steady sheet flow experiment for coarse light particles. It is found that the drag-induced turbulence damping effect was more significant than the well-known density stratification for the flow condition and grain properties considered. Meanwhile, the turbulence-resolving model is able to reproduce bed intermittency, which was driven by turbulent ejection and sweep motions, similar to the laboratory observation. Finally, simulations for fine sand transport in oscillatory sheet flow demonstrate that the turbulence-resolving model is able to capture the enhanced transport layer thickness for fine sand, which may be related to the burst events near flow reversal. Several future research directions, including further improvements of the present modeling framework and science issues that may be significantly benefited from the present turbulence-resolving sediment transport framework, are recommended.













Modeling of Fluid-Solid Two-Phase Geophysical Flows


Book Description

Fluid-solid two-phase flows are frequently encountered in geophysical flow problems such as sediment transport and submarine landslides. It is still a challenge to the current experiment techniques to provide information such as detailed flow and pressure fields of each phase, which however is easily obtainable through numerical simulations using fluid-solid two-phase flow models. This chapter focuses on the Eulerian-Eulerian approach to two-phase geophysical flows. Brief derivations of the governing equations and some closure models are provided, and the numerical implementation in the finite-volume framework of OpenFOAM® is described. Two applications in sediment transport and submarine landslides are also included at the end of the chapter.




Three-dimensional Numerical Analysis of Flow Structure and Sediment Transport Process in Open Channels


Book Description

This research project focuses on the analysis and prediction of flow structures and sediment transport process in open channels by using three-dimensional numerical models. The numerical study was performed using the open source computational fluid dynamics (CFD) solver based on the finite volume method (FVM) – OpenFOAM. Turbulence is treated by means of the two main methodologies; i.e. Large Eddy Simulation (LES) and Reynolds-Averaged Navier–Stokes (RANS). The free surface is tracked using the Volume of Fluid method (VOF). In addition, a new multi-dimensional model for sediment transport based on the Eulerian two-phase mathematical formulation is applied. The results obtained from the different numerical configurations are verified and validated against experimental data sets published in important research journals. The main characteristics of the flow structures are studied by using three set-up cases in steady and unsteady-state (transient) hydraulic flow conditions. On the other hand, the new multi-dimensional model for sediment transport is applied to predict the local scour caused by submerged wall jet test-case. Non-uniform structured elements are used in the grid configuration of the computational domains. A mesh sensitivity analysis is performed in each test-case study in order to obtain independent grid results. This analysis provides a balance between accuracy and optimal computational time. The results demonstrate that the three-dimensional numerical configurations satisfactorily reproduce the temporal variation of the different variables under study with correct trends and high correlation with the experimental values. Regarding the analysis and prediction of the flow structures, the results show the importance of the turbulence approach in the numerical configuration. On the other hand, the results of the new multi-dimensional two-phase model allow to analyze the full dynamics for sediment transport (concentration profile). Although the numerical results are satisfactory, the application of three-dimensional numerical models in field-scale cases requires a high computational resource.