This article presents the current state-of-the-art understanding of underwater dilative slope failure (breaching). Experimental investigations are reviewed, providing critical insights into the underlying physics of breaching and pointing out knowledge gaps, which underscore the need for further research. Besides, field observations at several locations across the globe are outlined, highlighting the hazard of breaching and the need for effective coastal management strategies to mitigate the associated risks. Furthermore, existing methods for analyzing and predicting the slope failure evolution are discussed and reflected upon, including analytical approaches and numerical models, ranging from simplified 1D models to advanced 3D coupled flow-soil approaches. Lastly, open questions are posed and key future directions are identified to enhance our understanding of the breaching failure. Overall, this review paper provides a valuable resource for researchers and decision makers involved in slope stability and flow slide risk assessment.@en

The chapter gives an overview of the sediment dispersion generated by the mining process. Within the field of dredging engineering, ample experience is available regarding equipment, turbidity generated by equipment, and sediment transport processes. High up the environmental impact mitigation hierarchy are avoidance and minimization. That is where engineering can provide (part of) the solution. It is our aim to predict and consider how we can improve the mining process and equipment. Within this context, our focus is on those processes that are likely to take place close to the seabed. On the one hand, our work focuses on the prediction and reduction of the amount of sediment that might get suspended. On the other hand, considering the conditions under which the suspended sediment might be released in the most optimal way to reduce dispersion, we have performed and analysed small-scale and full-scale laboratory experiments of a hydraulic collector design and various dynamic sedimentation experiments.@en

The presence of mud layers on the bottom of ports and waterways can have negative effects on the hydrodynamic behaviour of marine vessels. This numerical study investigates the effect of muddy seabeds on the full-scale resistance of an oil tanker sailing straight ahead. The objective is to determine the influence of factors such as the densimetric Froude number, UKC and mud rheology at speeds between 3 and 9 knots. The numerical study is conducted using a finite-volume Reynolds-Averaged Navier–Stokes (RANS) flow solver combined with the Volume-Of-Fluid (VOF) method to capture the mud–water interface. At certain critical speeds, the presence of mud increased the ship’s total resistance by up to 15 times compared to the case with solid bottoms. The non-Newtonian rheology of mud was found to influence the ship’s resistance mainly at low speeds and when sailing through the mud layer. This article also shows that, when sailing through mud, the computed resistance at high speeds may be underestimated because of two effects, namely ‘water lubrication’ and ‘numerical ventilation’.@en

Accurate characterization of mechanical perturbations on the seabed is essential for developing models assessing the environmental impacts from physical disturbances. Furthermore, understanding the relationship between (1) seabed resistance and (2) penetration depth, can also facilitate the development of more efficient and less impactful fishing gears. This study examines these two aspects of tickler chain rigged beam trawling via large-scale physical experiments. Three scaled down models (“light,” “medium,” and “heavy” designs) were developed to represent the impacts from typical beam trawl configurations used in the North Sea and were towed at various speeds on a saturated sand bed. Results reveal that increasing the towing speed reduces the mean penetration depth and the steady-state towing resistance of the gears. Smaller scale physical model tests incorporating tickler chains in sand, demonstrate that the towing resistance is significantly influenced by the soil compaction and particle sizes. Moreover, our study offers a simple and efficient method to estimate the penetration depth and towing resistance of prototype beam trawl gears in sand. This approach, along with the associated research, may be valuable for marine scientists assessing trawling impacts and demersal fishing gear designers seeking to optimize efficiency while minimizing seabed disturbance.

@en

A common way to transport solids in large quantities is by using a carrier fluid to transport the solids as a concentrated solid/liquid mixture or slurry through a pipeline. Typical examples are found in dredging, mining and drilling applications. Dependent on the slurry properties and flow conditions, horizontal slurry pipe flow is either in the fixed-bed, sliding-bed or fully-suspended regime. In terms of non-dimensional numbers, the flow is fully characterized by the bulk liquid Reynolds number (Re), the Galileo number (Ga, a measure for the tendency of particles to settle under gravity), the solid bulk concentration (ϕ_b), the particle/fluid density ratio (ρ_p/ρ_f), the particle/pipe diameter ratio (D_p/D_pipe), and parameters related to direct particle interactions such as the Coulomb coefficient of sliding friction (μ_c). To further our fundamental understanding of the flow dynamics, we performed experiments and interface-resolved Direct Numerical Simulations (DNS) of slurry flow in a horizontal pipe. The experiments were performed in a transparent flow loop with D_pipe=4 cm. We measured the pressure drop along the pipeline, the spatial solid concentration distribution in the cross-flow plane through Electrical Resistance Tomography (ERT), and used a high-speed camera for flow visualization. The slurry consisted of polystyrene beads in water with D_p=2mm, ρ_p/ρ_f=1.02, Ga between 40–45 and ϕ_b between 0.26–0.33. The different flow regimes were studied by varying the flow rate, with Re varying from 3272 till 13830. The simulations were performed for the same flow parameters as in the experiments. Taking the experimental uncertainty into account, the results from the DNS and the experiments are in reasonably good agreement. The results for the pressure drop agree also fairly well with popular empirical models from literature. In addition, we performed a parametric DNS study in which we solely varied Re and Ga. In all flow regimes, a secondary flow of Prandtl's second kind is present, ascribed to the presence of internal flow corners and a ridge of densely packed particles at the pipe bottom during transition towards the fully-suspended regime. In the bulk of the turbulent flow above the bed, secondary flow transport of streamwise momentum dominates over turbulent diffusion in regions where the secondary flow is strong and vice versa where it is weak. The transition between flow regimes appears to be governed by the competition between the net gravity force on the particles and shear-induced particle migration from particle–particle interactions. This competition can be expressed by the Shields number, θ. For θ≲0.75, gravity is dominant and the flow is in the fixed-bed regime. For θ≳0.75, shear-induced migration becomes progressively more important for increasing θ. Low-concentration zones flanking the sliding bed start to form at the top corners of the bed, and gradually expand downwards along the pipe wall till the pipe bottom is reached. For θ≳1.5, shear-induced migration is responsible for lifting the particle bed away from the wall, associated with the onset of the suspended regime. For θ≫1, gravity is of minor importance and the mean flow eventually reaches axi-symmetry with a high-concentration particle core at the pipe center and negligible secondary flow.

@en

As a result of the dilation of soil matrix, dense submarine sand slopes can temporarily be steeper than the natural angle of repose. These slopes gradually fail by the detachment of individual grains and intermittent collapses of small coherent sand wedges. The key question is whether steep disturbances in a submarine slope grow in size (destabilizing breaching) or gradually diminish (stabilizing breaching) and thereby limit the overall slope failure and resulting damage. The ability to predict whether the breaching failure is stabilizing or destabilizing is also crucial for the assessment of safety of submarine infrastructure and hydraulic structures located along rivers, lakes, and coasts. Through a set of large-scale laboratory experiments, we investigate the validity of an existing criterion to determine the failure mode of breaching (i.e., stabilizing or destabilizing). Both modes were observed in these experiments, providing a unique set of data for analysis. It is concluded that the existing method has limited forecasting power. This was quantified using the mean absolute percentage error, which was found to be 92%. The reasons behind this large discrepancy are discussed. Given the complexity of the underlying geotechnical and hydraulic processes, more advanced methodologies are required.@en

Reliable prediction of the erosion rate of sediment beds is important for many applications in coastal and river engineering. Theoretical understanding of empirically derived scaling relations is still lacking. This applies in particular for the scaling anomaly between low and high Shields number conditions. In this work, the erosion process is studied from the perspective of the phase-averaged turbulent kinetic energy (TKE) equations. The multi-phase TKE equations are written in a form that allows for a direct comparison with the TKE equation that appears for a stratified single-phase flow under the Boussinesq approximation. This reveals that next to buoyancy destruction, several other TKE modulation mechanisms become important at high Shields numbers and concentrations. Two scaling laws are derived for both moderate and high Shields numbers, and are tested against a wide range of experimental data.@en

In submerged sandy slopes, soil is frequently eroded as a combination of two main mechanisms: breaching, which refers to the retrogressive failure of a steep slope forming a turbidity current, and instantaneous sliding wedges, known as shear failure, that also contribute to shape the morphology of the soil deposit. Although there are several modes of failures, in this paper we investigate breaching and shear failures of granular columns using the two-fluid approach. The numerical model is first applied to simulate small-scale granular column collapses (Rondon et al., Phys. Fluids, vol. 23, 2011, 073301) with different initial volume fractions to study the role of the initial conditions in the main flow dynamics. For loosely packed granular columns, the porous medium initially contracts and the resulting positive pore pressure leads to a rapid collapse. Whereas in initially dense-packing columns, the porous medium dilates and negative pore pressure is generated stabilizing the granular column, which results in a slow collapse. The proposed numerical approach shows good agreement with the experimental data in terms of morphology and excess of pore pressure. Numerical results are extended to a large-scale application (Weij, doctoral dissertation, 2020, Delft University of Technology; Alhaddad et al., J. Mar. Sci. Eng., vol. 11, 2023, 560) known as the breaching process. This phenomenon may occur naturally at coasts or on dykes and levees in rivers but it can also be triggered by humans during dredging operations. The results indicate that the two-phase flow model correctly predicts the dilative behaviour and the subsequent turbidity currents associated with the breaching process.

@en

When dredging rock using a Cutter Suction Dredger the high amount of spillage is problematic, since it prevents an energy efficient removal process. This papers presents a coupled DEM-FVM method to simulate spillage, that can be used for optimizing the design and working method of the Cutter Suction Dredger. In these simulations, the challenge was to model relatively large particles in a complex and rotating geometry. To ensure stability and reduce computational time we used smoothing kernels to map the forces and the concentration between the discrete elements and the fluid mesh. The method is validated for the fluid flow in the rotating cutter head. This model incorporates all physical processes to predict flow induced spillage in cutter heads within feasible calculation times.@en

This article presents a new turbulence closure based on the k-ω SST model for predicting turbulent flows of Herschel–Bulkley fluids, including Bingham and power-law fluids. The model has been calibrated with direct numerical simulations (DNS) data for fully-developed pipe flow of shear-thinning and viscoplastic fluids. The new model shows good agreement in the mean velocity, average viscosity, mean shear stress budget and friction factor. The latter compares well also against correlations from the literature for a wide range of Reynolds numbers. With the new model, improvements are also observed in the iterative convergence, which is often difficult for calculations with yield-stress fluids. Additionally, three eddy-viscosity models for Newtonian fluids, namely the k-ω SST, k-kL and Spalart–Allmaras model, have been tested on turbulent Herschel–Bulkley flows. Results show that (i) the new model produces the best prediction; (ii) the standard SST model may be considered for simulations of weakly shear-thinning/viscoplastic fluids at high Reynolds numbers; (iii) the k-kL and the Spalart–Allmaras models appear to be unsuitable for turbulent Herschel–Bulkley flows. The new model is simple and appealing for engineering applications concerned with turbulent wall-bounded flows and is presented in a formulation that can be easily adapted to other generalised Newtonian fluids.

@en

When investigating the effect of muddy seabeds on marine vessels using Computational Fluid Dynamics (CFD) software, one challenge is to adequately describe the complex non-Newtonian fluid behaviour of mud. Although a number of rheological models have been proposed in the past, mud sediments are often simply regarded either as highly viscous Newtonian fluids or as Bingham fluids in many engineering applications. In this study, we investigate the accuracy of the Bingham model for numerical predictions of the viscous forces on a plate moving through fluid mud in laminar regime. In this context, a plate could be regarded as the flat bottom of a ship hull. The aim is to provide CFD practitioners with information about the accuracy of the Bingham model for the prediction of the frictional resistance of a ship sailing through fluid mud. This work presents a comparison of experimental and numerical data on the resistance of a plate moving through fluid mud from the Europoort area (Netherlands). Results suggest that the regularised Bingham model can be a reasonable compromise between simplicity and accuracy for CFD simulations to investigate the effect of muddy seabeds on marine vessels. A comparison between CFD data and analytical formulas is also presented.@en

Background: Both influenza and SARS-CoV-2 viruses show a strong seasonal spreading in temperate regions. Several studies indicated that changes in indoor humidity could be one of the key factors explaining this. Objective: The purpose of this study is to quantify the association between relevant epidemiological metrics and humidity in both influenza and SARS-CoV-2 epidemic periods. Methods: The atmospheric dew point temperature serves as a proxy for indoor relative humidity. This study considered the weekly mortality rate in the Netherlands between 1995 and 2019 to determine the correlation between the dew point and the spread of influenza. During influenza epidemic periods in the Netherlands, governmental restrictions were absent; therefore, there is no need to control this confounder. During the SARS-CoV-2 pandemic, governmental restrictions strongly varied over time. To control this effect, periods with a relatively constant governmental intervention level were selected to analyze the reproduction rate. We also examine SARS-CoV-2 deaths in the nursing home setting, where health policy and social factors were less variable. Viral transmissibility was measured by computing the ratio between the estimated daily number of infectious persons in the Netherlands and the lagged mortality figures in the nursing homes. Results: For both influenza and SARS-CoV-2, a significant correlation was found between the dew point temperature and the aforementioned epidemiological metrics. The findings are consistent with the anticipated mechanisms related to droplet evaporation, stability of virus in the indoor environment, and impairment of the natural defenses of the respiratory tract in dry air. Significance: This information is helpful to understand the seasonal pattern of respiratory viruses and motivate further study to what extent it is possible to alter the seasonal pattern by actively intervening in the adverse role of low humidity during fall and winter in temperate regions. Impact: A solid understanding and quantification of the role of humidity on the transmission of respiratory viruses is imperative for epidemiological modeling and the installation of non-pharmaceutical interventions. The results of this study indicate that improving the indoor humidity by humidifiers could be a promising technology for reducing the spread of both influenza and SARS-CoV-2 during winter and fall in the temperate zone. The identification of this potential should be seen as a strong motivation to invest in further prospective testing of this non-pharmaceutical intervention.

@en

The ship’s resistance and manoeuvrability in shallow waters can be adversely influenced by the presence of fluid mud layers on the seabed of ports and waterways. Fluid mud exhibits a complex non-Newtonian rheology that is often described using the Herschel–Bulkley model. The latter has been recently implemented in a maritime finite-volume CFD code to study the manoeuvrability of ships in the presence of muddy seabeds. In this paper, we explore the accuracy and robustness of the CFD code in simulating the flow of Herschel–Bulkley fluids, including power-law, Bingham and Newtonian fluids as particular cases. As a stepping stone towards the final maritime applications, the study is carried out on a classic benchmark problem in non-Newtonian fluid mechanics: the laminar flow around a sphere. The aim is to test the performance of the non-Newtonian solver before applying it to the more complex scenarios. Present results could also be used as reference data for future testing. Flow simulations are carried out at low Reynolds numbers in order to compare our results with an extensive collection of data from the literature. Results agree both qualitatively and quantitatively with literature. Difficulties in the convergence of the iterative solver emerged when simulating Bingham and Herschel–Bulkley flows. A simple change in the interpolation of the apparent viscosity has mitigated such difficulties. The results of this work, combined with our previous code verification exercises, suggest that the non-Newtonian solver works as intended and it can be thus employed on more complex applications.@en

Dredging is the relocation of soil. Before the soil can be transported, it has to be loosened. This can be done hydraulically (jetting) or mechanically (cutting). Often, water jets are used to erode the soil layer. Over time, pickup functions have been derived to predict the amount of erosion corresponding to the flow conditions. However, existing pickup functions are inaccurate at high flow velocities. During the current study, erosion experiments have been done at high flow velocities (up to 4.7 m/s) corresponding to a bed shear stress of up to 60 Pa and a Shields parameter (θ) of up to 30. The results of these experiments were compared with a number of well-known data sets and pickup functions.

@en

In dredging, high pressure water jets are commonly applied to assist the mobilization of soil. This work considers the excavation of cohesive soil. The key objective is to predict the development of the cavity in the soil as a function of the undrained shear strength, translation velocity and hydrodynamic pressure of a single nozzle. A generic computational fluid dynamics (CFD) model has been developed that captures both the jet flow and the soil failure in a single framework. The results are compared with data from a previous experimental study. The CFD model predicts the cavity dimensions with reasonable accuracy. In addition the model provides detailed data to study the cyclic nature of the soil failure process. The CFD model is promising and can be applied for more complex nozzle configurations to assist the design process of dragheads and improve production estimates.

@en

The presence of complex fluids in nature and industrial applications combined with the rapid growth of computer power over the past decades has led to an increasing number of numerical studies of non-Newtonian flows. In most cases, non-Newtonian models can be implemented in existing Newtonian solvers by relatively simple modifications of the viscosity. However, due to the scarcity of analytical solutions for non-Newtonian fluid flows and the widespread use of regularization methods, performing rigorous code verification is a challenging task. The method of manufactured solutions (MMS) is a powerful tool to generate analytical solutions for code verification. In this article, we present and discuss the results of three verification exercises based on MMS: (i) steady single-phase flow; (ii) unsteady two-phase flow with a smooth interface; (iii) unsteady two-phase flow with a free surface. The first and second exercises showed that rigorous verification of non-Newtonian fluid solvers is possible both on single- and two-phase flows. The third exercise revealed that “spurious velocities” typical of free-surface calculations with the Volume-of-Fluid model lead to “spurious viscosities” in the non-Newtonian fluid. The procedure is illustrated herein on a second-order finite volume flow solver, using the regularized Herschel-Bulkley fluid model as an example. The same methodology is however applicable to any flow solver and to all the rheological models falling under the class of generalized Newtonian fluid models.@en

In dredging applications, deep sea mining and land reclamation projects typically large amounts of sediments are transported through pipes in the form of hyper concentrated (40% sediment or more) sediment-water mixtures or slurries. In this paper it is investigated how well a generic Euler-Euler CFD-model is capable to model velocity, concentration profiles and the pressure gradient of sediment above deposition limit velocity in a pipeline. This Euler-Euler solver treats both phases as a continuum with its own momentum and continuity equations. The full kinetic theory for granular flows is accounted for (no algebraic form is used) and is combined with a buoyant k-ε turbulence model for the fluid phase. The influence of the mesh size has been checked and grid convergence is achieved. All numerical schemes used are of second-order accuracy in space. The pressure gradient was calibrated by adjusting the specularity coefficient in one calibration case and kept constant afterwards. Simulations were carried out in a wide range of slurry flow parameters, in situ volume concentration (9–42%), pipe diameter (0.05–0.90 m), particle diameter (90–440 μm) and flow velocity of (3–7 m/s). The model shows satisfactory agreement to experimental data from existing literature.

@en

In dredging applications, deep sea mining and land reclamation projects typically large amounts of sediments are transported through pipes in the form of hyper concentrated (40% sediment or more) sediment-water mixtures or slurry. These slurries can flow at three different regimes. 1: fully suspended (homogeneous and stratified) 2: partially suspended with a sliding bed 3: partially suspended with a fixed bed. At the moment it is hard to predict the transport regime, the volume flux of particles and the pressure drop (friction factor) of these slurries within these regimes. The goal is to establish a model 3D continuum model that is able to predict the aforementioned aspects of slurry flow in a wide range of slurry flow conditions. In this paper it is investigated how well an existing CFD-model is able to model velocity and concentration profiles of suspended sediment in a pipeline. The CFD-model that is used is TwoPhaseEulerFoam which is part of OpenFOAM. This Euler-Euler solver treats both the phases as a continuum with its own momentum and continuity equations. The phases are coupled with coupling terms such as the drag force.@en