GK

G.H. Keetels

info

Please Note

36 records found

Journal article (2025) - S.M.S. Alhaddad, A.L. Snyder, S.V. Bult, G.H. Keetels
Very limited research has been carried out to investigate sediment erosion caused by subaqueous inclined water jets, despite the fact that such water jets are used in subsea engineering (e.g., dredging, trenching, and deep sea mining). Therefore, we conducted a set of novel small-scale experiments to primarily study the effect of jetting inclination on cohesive sediment erosion. The experimental results reveal that vertical jetting results in the largest cavity depth (or ’erosion depth’), but not in the largest cavity size (sediment production). The erosion depth increases with the jetting angle reaching its maximum at 90° and then begins to decrease with further increase in the jetting angle. The results also indicate that the cavity width (or ’erosion width’) is not necessarily correlated with the impingement region but is instead associated with the erosion-effective jet width—the width of the jet where flow velocities are high enough to penetrate the bed. Analysis of the cavity size showed that the largest sediment production was achieved at a 65° jetting angle among the tested jetting angles (25°, 45°, 65°, 90°, 115°, 135°, and 155°). The erosion depth was found to be highly proportional to the impingement force exerted by the flow on the clay. ...
This study investigates the influence of multiple jet parameters on the flow field of translating impinging inclined water jets. We conducted full-scale stereoscopic particle image velocimetry and pressure measurements and three-dimensional computational fluid dynamics simulations for Reynolds numbers in the range of. Considering the complex mechanism of a translating impinging jet, a good concordance is observed between the experimental and numerical results. The translation-to-jet velocity ratio is identified as a critical parameter in determining whether the jet flow predominantly exhibits impinging characteristics or behaves as a jet in cross-flow. It is found that, for, jet impingement is minimal. The stand-off distance to nozzle diameter ratio determines the relative influence of the cross-flow on the jet flow. The effect of is similar to a stationary impinging jet, with the potential core extending up to, but entrainment is enhanced by the relative cross-flow. For an inclined jet, i.e. jet angle, the direction of the jet, either backward or forward, governs the deflection of the flow. Higher pressures are recorded for a backward directed jet compared with a forward directed jet for supplementary angles. ...
Journal article (2025) - Mosaab Mahgoub, G.H. Keetels, Said Alhaddad
The cutter suction dredger (CSD) is one of the main vessels utilized in the dredging industry. The dynamic actions related to its rotating cutter head are the main trigger for sediment release and turbidity generation by this vessel. The ability to predict the evolution of this turbidity and suspended sediment concentrations is imperative for effective environmental management. This predictive capability allows ecologists to estimate potential damage, enabling environmental managers to propose appropriate mitigation measures. In this study, we conducted a qualitative numerical assessment of the characteristics of turbidity currents generated as a result of cutter suction dredging of densely-packed sand. To this end, we developed a one-dimensional physics-based model, providing the order of magnitude of sediment fluxes and concentration levels. In addition, a quantitative sensitivity analysis is performed to unravel the relative influence of key operational parameters on the generated turbidity currents. The results of this research reveal that breaching (dilative underwater slope failure) is a major source of sediment release by CSDs and should be incorporated in the source-term estimation. It is also found that the cut ratio is the most influential operational parameter on the generated turbidity. ...

Insights from experiments and interface-resolved Direct Numerical Simulations

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 (ρpf), the particle/pipe diameter ratio (Dp/Dpipe), 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 Dpipe=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 Dp=2mm, ρpf=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. ...
Journal article (2024) - S. Lovato, S.L. Toxopeus, J.W. Settels, G.H. Keetels, A. Kirichek
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’. ...
Journal article (2024) - Bithin Ghorai, Justin Tiano, Pieke Molenaar, Karline Soetaert, Geert Keetels
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. ...
Review (2024) - S.M.S. Alhaddad, G.H. Keetels, Dick Mastbergen, C. van Rhee, Cheng-Hsien Lee, Eduard Puig Montellà, Julien Chauchat
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. ...
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. ...
Journal article (2023) - S.M.S. Alhaddad, Dave Weij, C. van Rhee, G.H. Keetels
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. ...
Journal article (2023) - E. P. Montellà, J. Chauchat, C. Bonamy, D. Weij, G. H. Keetels, T. J. Hsu
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. ...
Journal article (2023) - B.J. Nieuwboer, C. van Rhee, G.H. Keetels
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. ...
Journal article (2023) - G.H. Keetels, Julien Chauchat, W.P. Breugem
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. ...
Journal article (2022) - G. H. Keetels, L. Godderis, B. J.H. van de Wiel
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. ...

CFD modelling and towing tank experiments

Journal article (2022) - S. Lovato, Alex Kirichek, S.L. Toxopeus, J.W. Settels, G.H. Keetels
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. ...
Journal article (2022) - S. Lovato, S.L. Toxopeus, J.W. Settels, G.H. Keetels
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. ...
Journal article (2022) - Oliver A. Heijmeijer, Arno J. Nobel, Geert Keetels, Cees Van Rhee
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. ...
Journal article (2022) - S. Lovato, G. H. Keetels, S. L. Toxopeus, J. W. Settels
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. ...
Journal article (2021) - S.L. Lovato, Serge Leon Toxopeus, Just W. Settels, G.H. Keetels, Guilherme Vaz
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. ...