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

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

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

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