Search for: All records

Abstract. In this paper, an Eulerian two-phase flow model, sedFoam, is extended to include an air phase together with its water and sediment phases. The numerical model called sedInterFoam is implemented using the open-source library OpenFOAM. sedInterFoam includes the previous features of sedFoam for sediment transport modeling and also solves the air–water interface using the volume-of-fluid method coupled with the waves2Foam toolbox for free-surface wave generation and absorption. Using sedInterFoam, four test cases are successfully reproduced to validate the free-surface evolution algorithm's implementation, mass conservation of sediment and fluid phases, and predictive capabilities and to demonstrate its potential in modeling a broader range of coastal applications with sediment transport dominated by surface waves. 

The present study evaluates the performance of two numerical approaches in estimating non-equilibrium scour patterns around a non-slender square structure subjected to a transient wave, by comparing numerical findings with experimental data. This study also investigates the impact of the structure’s positioning on bed evolution, analyzing configurations where the structure is either attached to the sidewall or positioned at the centerline of the wave flume. The first numerical method treats sediment particles as a distinct continuum phase, directly solving the continuity and momentum equations for both sediment and fluid phases. The second method estimates sediment transport using the quadratic law of bottom shear stress, yielding robust predictions of bed evolution through meticulous calibration and validation. The findings reveal that both methods underestimate vortex-induced near-bed vertical velocities. Deposits formed along vortex trajectories are overestimated by the first method, while the second method satisfactorily predicts the bed evolution beneath these paths. Scour holes caused by wave impingement tend to backfill as the flow intensity diminishes. The second method cannot sufficiently capture this backfilling, whereas the first method adequately reflects the phenomenon. Overall, this study highlights significant variations in the predictive capabilities of both methods in regard to the evolution of non-equilibrium scour at low Keulegan–Carpenter numbers. 

Abstract. In this paper, an Eulerian two-phase-flow model sedFoam is extended to include an air phase together with the water and sediment phases. The numerical model called sedInterFoam is implemented using the open source library OpenFOAM. SedInterFoam includes the previous features of sedFoam for sediment transport modeling and also solves the air/water interface using the volume of fluid method coupled with the waves2Foam toolbox for free surface wave generation and absorption. Using sedInterFoam, four test cases are successfully reproduced to validate the free-surface evolution algorithm implementation, mass conservation of sediment and fluid phases, predictive capabilities and demonstrate its potential in modelling a broader range of coastal applications with sediment transport dominated by surface waves. 

A Reynolds-averaged two-phase Eulerian model for sediment transport, SedFoam, is utilized in a twodimensional domain for a given sediment grain size, flow period, and mobility number to study the asymmetric and skewed flow effects on the sediment transport over coarse-sand migrating ripples. First, the model is validated with a full-scale water tunnel experiment of orbital ripple driven by acceleration skewed (asymmetric) oscillatory flow with good agreement in the flow velocity, net sediment transport, and ripple migration rate. The model results showed that the asymmetric flow causes a net onshore sediment transport of both suspended and near-bed load (the conventional bed load and part of the near-bed suspended load, responsible for ripple migration). The suspended load transport is driven by the “positive phase-lag” effect, while the near-bed transport is due to the large erosion of the boundary layer on the stoss flank, sediment avalanching on the lee flank, and the returning flux induced by the stoss vortex. Together, these processes result in a net onshore transport rate. In contrast, for an energetic velocity skewed (skewed) flow, the net transport rate is offshore directed. This is due to a larger offshore-directed suspended load transport rate, resulting from the “negative phase-lag” effect, compared to the onshore-directed near-bed load transport rate. Compared to the asymmetric flow, the onshore near-bed load transport (and migration) rate is limited by the larger offshore directed flux associated with returning flow on the lee side, due to a stronger lee vortex generation during the onshore flow half-cycle. In the combined asymmetric-skewed case, the near-bed load and migration rate are higher than in the asymmetric flow case. Moreover, the offshore-directed suspended load is much smaller compared to the skewed flow case due to a competition between the negative (due to velocity skewness) and positive (due to acceleration skewness) phase-lag effects. As a result, the net transport rate is onshore directed but slightly smaller than the asymmetric flow case. 

Subaqueous vortex ripples in equilibrium are characterized by their unique geometry and dimensions. Motivated by the recent direct numerical simulation study of oscillatory turbulent flow over a wavy bottom by Önder & Yuan ( J. Fluid Mech. , vol. 858, 2019, pp. 264–314), the objective of this study is to further investigate the fluid dynamical controls that determine the distinctive equilibrium dimensions of vortex ripples. We use direct numerical simulations to investigate the differences in flow kinetics between sinusoidal oscillatory flow over equilibrium and out-of-equilibrium vortex ripples. In comparison with the equilibrium case, the spanwise coherent vortices, the averaged bottom shear stress on overlying flow and the shear stress distribution on the ripple surface are identified as the key fluid dynamical controls on equilibrium dimensions. Based on these controls, we propose mechanisms in the selection of vortex ripple dimensions. We observe that the flow adjusts in such a way that the interaction between overlying flow and vortex ripples tends to generate the strongest coherent vortices while the ripple surface (or overlying flow) experiences the smallest shear stress averaged over ripple wavelength during the selection process. Through a triple decomposition of the flow, the component of the ripple-induced fluctuation is found to dictate these fluid dynamical controls, which implies that this component plays an important role in the evolution of vortex ripples. 

Abstract A 3D large eddy simulation coupled with a free surface tracking scheme was used to simulate cross‐shore hydrodynamics as observed in a large wave flume experiment. The primary objective was to enhance the understanding of wave‐backwash interactions and the implications for observed morphodynamics. Two simulation cases were carried out to elucidate key processes of wave‐backwash interactions across two distinct stages: berm erosion and sandbar formation, during the early portion of a modeled storm. The major difference between the two cases was the bathymetry: one featuring a berm without a sandbar (Case I), and the other, featuring a sandbar without a berm (Case II) at similar water depth. Good agreement (overall Willmott's index of agreement greater than 0.8) between simulations and measured data in free surface elevation, wave spectrum, and flow velocities validated the model skill. The findings indicated that the bottom shear stress, represented by the Shields parameter, was significant in both cases, potentially contributing substantial sediment transport. Notably, the occurrence of intense wave‐backwash interactions were more frequent in the absence of a sandbar. These intense wave‐backwash interactions resulted in a pronounced horizontal pressure gradient, quantified by high Sleath parameters, exceeding the criteria for momentary bed failure. Additionally, a more vigorous turbulence‐bed interaction, characterized by near‐bed turbulent kinetic energy, was observed in the case lacking a sandbar, potentially augmenting sediment suspension. These insights are pivotal in understanding the mechanisms underlying berm erosion and how sandbar formation serves to protect further beach erosion. 

The flocculation behavior of clay minerals in aquatic environments is an important process in estuarine and riverine dynamics, where strong gradients in salinity can locally occur. Various contradicting observations have been reported in the literature on the impact of salt concentration on the settling process of cohesive sediments. To address this issue in a systematic manner, we investigate the settling behavior of clay minerals as a function of the salt concentration of the ambient water. Specifically, we focus on montmorillonite as a prototype clay mineral with a high cation exchange capacity (CEC). To this end, we study suspensions of Wyoming bentonite (Volclay SPV) as a very important constituent for many constructional and industrial purposes. We perform an experimental campaign to study the settling behavior of moderately turbid montmorillonite concentrations in monovalent salt solutions with different salinities (sodium chloride) to represent different environments ranging from deionized to ocean water, respectively. The subsequent settling process was monitored by taking pictures by a camera in regular time intervals over a total observation time up to 48 h. In addition, a modified hydrometer analysis is conducted to determine the grain size distribution (in terms of an equivalent diameter) of the flocculated clay suspension in salt water. Despite the rather high cation exchange capacity of the investigated clay (CEC=88.1), our results show that the settling speed drastically increases within a range of 0.6–1.0 PSU and stays approximately constant for higher salinities. This critical salt concentration is defined here as the critical coagulation concentration (CCC) and lies well below the salinity of natural open water bodies. The hydrometer analysis revealed that 60% of the agglomerates exceed the equivalent grain size of 20 μm. Finally, the findings of this study are supplemented with experiments studying the effect of Extracellular Polymeric Substances (EPS) on the flocculation behavior of bentonite in salt water. Our results demonstrate that salinity is the original trigger for flocculation, whereas EPS allows for even larger floc size but it does not play a significant role for the settling processes of bentonite in estuarine environments.