Next Article in Journal
Design of a Self-Organizing Routing Protocol for Underwater Wireless Sensor Networks Based on Location and Energy Information
Next Article in Special Issue
Exploring the Impact of Palladium Loading on Pd-Based Three-Way Catalyst Performance and Propane Reactivity for Emission Reduction in Liquefied Petroleum Gas Engines
Previous Article in Journal
Hybrid Model Predictive Control of a Two-Body Wave Energy Converter with Mechanically Driven Power Take-Off
Previous Article in Special Issue
Update on Risk Criteria for Crude Oil Tanker Fleet
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sediment Resuspension Distribution Modelling Using a Ship Handling Simulation along with the MIKE 3 Application

Faculty of Maritime Studies and Transport, University of Ljubljana, Pot Pomorščakov 4, 6320 Porotorož, Slovenia
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2023, 11(8), 1619; https://doi.org/10.3390/jmse11081619
Submission received: 10 July 2023 / Revised: 16 August 2023 / Accepted: 17 August 2023 / Published: 19 August 2023
(This article belongs to the Special Issue Sustainable and Safe Maritime Transportation)

Abstract

:
The environmental effects of ship propellers were not even close to fully examined before the current massive ships were introduced to sea trade. Larger ships, result in greater length, beam, draft and propulsion power. Of concern here is the under-keel clearance (UKC) and applied power, the most important parameters causing sea bottom sediment resuspension and, consequently, the transport and deposition of washed sediments. The problems are multifarious: shorelines could be contaminated with heavy metals, petroleum hydrocarbons and other organic chemicals, which are sometimes buried deep in the sediment bed. The effects of resuspension on marine life have been well documented by marine biologists. Further, a ship passing through a flow field may have a significant hydrodynamic effect on the shipping channel: waves generated by moving vessels can accelerate shoreline erosion; erosion around quay piles have a negative impact on sea flora. Waves can also affect other manoeuvring vessels or ships at berth. Available empirical models are applicable for a steady state condition, addressing velocity and, consequently, shears at the sea bottom for defined UKC and states of applied power. The idea here is to calculate material resuspension dynamically in the water column based on realistic manoeuvring conditions, which can be a matter of some complexity. During a manoeuvre, the pilot must bring the ship into or out of the harbour in the safest possible way, operating the telegraph, rudder, thrusters and possibly tugs, and also co-ordinating the work of the linesmen. The jet speed powering the vessel is not only a function of the speed of the propeller, but also of the present speed of the ship, which has an effect on the propeller’s constantly changing torque. Additionally, the bathymetry is constantly changing, and the streamlines hit not only the seabed, but also the bank and other structures of the harbour basin. The resuspended material remains in the column long after the ship has finished manoeuvring, moving slowly through the entire water column and being transported not only by the remaining streamlines of the ship but also by general currents. Realistic manoeuvring parameters can be obtained from real-time simulations with a real crew using state-of-the-art Full Mission Bridge Simulators (FMBS); eddies and the like contribute to the distribution and material resuspension and can be calculated by applying numerical modelling. In this paper, a container ship departure manoeuvre is simulated dynamically using Wartsila FMBS obtained data, which are postprocessed and coupled with the MIKE 3 FM hydrodynamic modelling application to which we add the precise port of Koper bathymetry to gain ship propeller spatial jet velocity distribution in specific time domains. Obtained jet velocity distribution is further coupled with the MIKE 3 MT particle tracking application to visualize total resuspended sediment transport patterns, etc. Container ships were selected to amplify the urgency of this phenomenon; they are the most intrusive in terms of resuspending and scouring the seabed given their powerful engines and larger propellers. Passenger ships could have been used, car carriers, or even tankers; but the fear among scientists is that the issue will not be taken seriously enough by certain stakeholders.

1. Introduction

The profitability of larger ships is blinding the community of maritime commerce to the numerous problems that form a chain trailing from the inexorable drive to unquestionably float ever larger ships. Ports around the world are struggling to accommodate deep draft vessels, and though this paper studies the northern Adriatic near the Port of Koper, which sees more than its share of pollution, we are also keenly aware of the safety issues caused by the premature predominance of overly large vessels (collisions, groundings, near misses, delays, increased effects of weather). The trend of the future is toward ships that call at ports with low under-keel clearance (UKC), with less and less distance between the seabed and the propeller tip and greater main engine and thrusters’ power. Propeller thrust, propeller diameter and the other factors mentioned above create a propeller jet velocity field that results in sediment resuspension (SR), which has known negative effects on sea flora and fauna, berth walls, quay piles (erosion) and the deposits of resuspended sediment [1]. This paper deals with sediment resuspension in the Bay of Trieste. This location is specific because of the Idrijca River, which deposits mercury in this region [2]. The process of sediment resuspension by ship propellers raises mercury particles that enter the marine food chain and, interestingly, deliver some mercury to the very mouth of the river that flows into the gulf.
The study area is shown in Figure 1, where three basins of the port area are visible; the container terminal is located in the first basin, close to the city centre.
The velocity field of the propeller jet reaches the seafloor at some distance from (behind) the vessel, depending on the sea depth and the topography of the seabed. Critical jet velocity is the speed at which the particles overcome their own shear stress and resedimentation occurs; obviously, this varies according to the particle characteristics. Additionally, metocean conditions such as currents, density gradients and prevailing winds affect sediment transport and deposition. After the resuspension process, sediment much of which is too dispersed to be tracked and microparticles may travel some miles before settling is deposited in visible masses at predictable locations within the port basins. The bathymetry of the port is constantly changing (this is darkly ironic, as it may require dredging where at much expense dredging was required that the port might accept these very ships, as will be mentioned) and may prevent vessels from docking and manoeuvring at a particular point due to insufficient UKC. Port authorities must perform dredging to allow deep draft vessels to safely enter and berth in the port. The work is time-consuming, causes congestion in the port and creates additional costs due to demurrage costs. The first step in finding a solution to minimize SR is to understand the velocity distribution of the propeller jet. Tools to evaluate this thrust include the Wartsila NTPro 5000 v 5.35 full mission bridge simulator, and now the introduction of hydrodynamic and mud transport models MIKE 3 FM and MT [3].
Recent empirical modelling studies have included ship propeller jet distribution in water columns at specific distance from propeller orifices, but not effects on sea bottom sediment resuspension [4]. Some existing approaches have been limited to coupling resuspended material transport in specific time domains, sea bed roughness, bathymetry, port basin shape and metocean impacts on sediment transport [5].
Several authors performed laboratory experiments to validate propeller rotation impact on the sea bottom. The research aim of the author [6] was to determine scouring patterns based on twin propeller configuration. The measurements were done with 2D Laser Doppler Velocimetry to determine scouring patterns in x and y directions. The limitation of the research is the two-propeller set-up, both left turning; vessels in reality have opposite turning propellers due to elimination of propeller side walk. The authors [7] investigated maximum scour depth, width, length and deposit height, with experimental (Dp = 8.2 cm) ship four blade propeller. Their experiment included three different propeller axial distances from the bed and three propeller rotational speeds. They scanned the bed surface with the photogrammetry technique combined with a Terrestrial Laser Scanner (TLS). The main contributions of this paper were that propellers generate two holes; one smaller directly under the propeller and one bigger downstream. The propeller scouring is not symmetrical in a lateral direction but is deeper on one side with higher resuspended particles deposition on the other side. These authors [8] presented experimental results of the propeller jet flow field around an open quay using the Particle Image Velocimetry (PIV) technique. The experiment included four longitudinal distances from the propeller face to quay toe slope. The results showed that maximum scour depth appears at the toe of the slope and decreases with longitudinal distance from the propeller face. This study was useful for understanding and predicting erosion processes in narrow channel approaches and along port quay infrastructure. The authors presented ship propeller scour patterns using an experimental set-up and stationary position of the propeller. An upgrade could be repeating the work with a propeller moving through the ambient water. Results in this research could be used only in the first step of a ship departure manoeuvre, when ship velocity is zero or close to it.

2. Tools and Methods for the Determination of Ship Propeller Jet Velocity Field Distribution

There are many ways to obtain the kinematic parameters of a ship’s motion. The best would be from the vantage of the ship’s team onboard taking VDR data; but of course this is logistically complex. Another option is to reproduce the manoeuvre using the state-of-the-art ship handling simulators operated by actual marine pilots, which is what we did. The simplest option which is insufficiently accurate is to use Automatic Identification System data transmitted by every vessel in the area.

2.1. Automatic Identification System and Full Mission Bridge Simulator

The automatic identification system (AIS) is a standard tool in the maritime industry, primarily intended for easy identification of ship traffic. Other applications of AIS include collision prevention, emission awareness, tracking oil spills, monitoring fishing activities, spatial planning and much more [9]. Data transmitted by AIS is used in scientific research because the system provides a large amount of statical and dynamic data: vessel name, MMSI number and call sign, type and dimension of vessels, vessel position, course over ground, speed over ground, heading; from this data distance and time of closest point of approach can be determined with any particular object in the vicinity. Extensive data from AIS were analysed for the port of Koper, Basin I, where container ships berth. The decision for selecting this particular location was taken regarding the high container ship traffic density (which is expected to grow in the near future) and low UKC, as little as 0.5 m. Based on acquired AIS data, a complete sailing departure route leg can be reproduced in the simulator. According to the preliminary analysed data, the outbound manoeuvre of a container ship was selected and simulated; in order to simplify the first modelling, only the propulsion system of the main engine was analysed, without tug assistance and bow thruster. The following table (Table 1) contains the relevant data from simulator container ship 1.
The figure below (Figure 2) presents the procedure for obtaining dynamical data of a container ship departure manoeuvre needed for further numerical modelling with MIKE 3. The FMBS simulation program was used to reproduce the departure manoeuvre in real time and place.
Figure 3 shows a screenshot of the simulator instructor station, showing the harbour basin 1 with high resolution bathymetry and the position of the ship at berth. During the manoeuvre, relevant ship dynamic data were recorded at a frequency of 0.5 Hz: vessel position (lon., lat.), speed over ground (SOG), course over ground (COG), heading direction, propeller rates, thrust force and rudder angle.
The thrust force (Te) of the ship is then calculated in the simulator according to the following equation [10]:
T e = C t J , p D P ρ w n 2 D p 4 1 t
Parameters in the equation are: thrust coefficient (Ct), which is a function of speed of advance (J), pitch to propeller diameter ratio (p/DP), water density (ρw) propeller rotation per second (n) and thrust deduction coefficient (t). Speed of advance is calculated with the equation:
J = V n D p
Parameter (V) is vessel speed. The thrust deduction coefficient takes into account the underwater shape of the vessel. For conventionally-shaped single screw propelled vessels, it can be obtained with the following expression [11]:
t = 0.001979 L w l B B C P 1 + 1.0585 C 10 0.00524 0.1418 D p 2 B d + 0.0015 C s t e r n
The parameters in the equation are: length on the waterline (Lwl); ship’s breadth (B); draft (d); coefficient (C10), which for most of the ship is equal to the ratio between the ship’s beam and waterline length, while coefficient (Cstern) varies from −25 to 10 based on the shape of the ships stern; and coefficient (CP1), which is a function of the prismatic coefficient and longitudinal center of buoyancy. The thrust was compared with the outflow velocity of the propeller jet at a nonzero ship speed (V0;v≠0). The propeller rates recorded by the simulator are further processed using semiempirical methods in order to obtain the applied power (Papp) and to determine the propeller jet exit velocity and the jet velocity distribution behind the propeller.

2.2. Semi-Empirical Equations to Calculate Propeller Jet Efflux Velocity and Velocity Distribution behind the Ship Propeller

The propeller revolutions per minute (RPMapp) are recorded during ship manoeuvres with a time step of 2 s. The power applied can be calculated using the following equation:
P a p p = R P M a p p R P M m a x 3 P m a x
where (RPMmax) is the maximum propeller revolutions per minute and (Pmax) is the maximum installed engine power. The propeller jet velocity is normally calculated using Equations (9) and (10), but in some cases the thrust coefficient is not known, in which case the efflux velocity (V0) can be obtained using the following expression [12]:
V 0 = C 3 P a p p . ρ w D p 2 1 3
Coefficient (C3) is equal to:
  • 1.17 for ducted propellers (propellers with nozzle);
  • 1.48 for free propellers.
The equation above is used for ship bollard pull conditions (delivered by a ship when pulling at zero speed). During an actual manoeuvre, a ship has a certain speed (inertia moment) that must be considered:
V 0 ; v 0 = V 0 1 V D p n
Parameter (V) is vessel speed in (m/s). Jet efflux velocity for a vessel with a speed of more than zero is an important parameter for further analyses of maximum propeller jet bottom velocities (Vb,max), which determine critical sediment resuspension velocities of specific sediment structure [13]:
V b , m a x = 0 = E V 0 ; v 0 D p h t b
where E = 0.71, b = 1.0 for seagoing ships with rudder; E = 0.42, b = 0.275 for seagoing ships without a rudder; E = 0.52, b = 0.275 for seagoing ships with a twin propeller and double rudder.
The distance between the propeller axis and the bottom of the sea is represented by the (ht) value.
h t = C + D p 2
The clearance (C) is the distance between the propeller tip and the seabed. The equation gives the maximum velocity of the propeller jet at the seafloor and provides information about the manoeuvrability of the vessel in terms of a lower sediment resuspension rate. For a more accurate assessment of how much sediment is resuspended during the ship’s manoeuvre, another method is needed that provides a distribution of jet velocity in the vertical and horizontal directions with respect to the propeller axis. There are several authors who have addressed this problem, such as [5,14]. Models are used for the same purposes.
A Schematic view separates the propeller jet velocity into three zones (Figure 4):
  • Efflux zone (0 < X/Dp < 0.35);
  • Zone of flow establishment (0.35 < X/Dp < 3.25), determined by [15];
  • Zone of established flow (3.25 < X/Dp < 50).
where parameter (X) presents horizontal distance behind the propeller.
Figure 4. Propeller jet velocity distribution [16].
Figure 4. Propeller jet velocity distribution [16].
Jmse 11 01619 g004
Propeller jet efflux velocity necessitates two equations. The first is based on axial momentum theory [17], taking into account thrust coefficient:
V 0 = 1.33 n D p C T
The second equation [18] takes into account propeller thrust (T), water density (ρw) and propeller jet contraction (D0). The equation is particularly useful when the FMBS simulation tool is coupled with MIKE 3 FM, which records dynamic ship data, among which is thrust force. It is used to determine the propeller outflow velocity, which is used for further jet velocity distribution.
V 0 = 1.13 1 D 0 T ρ w
The jet velocity field in Zone 1 is determined to be the same (V0 = Vx,r) in longitudinal (X) and vertical (r) directions, up to the axial distance (0 < X/Dp < 0.35) from the propeller.
The maximal propeller jet efflux velocity decay in Zone 2 is determined using the equation proposed by [4]:
V x , m a x V 0 = 1.51 0.175 X D p 0.46 p D P
The propeller jet velocity field (Vx,r_2) in Zone 2 was first determined by the Gaussian normal probability function [19]; later [20], this equation was improved and presented the jet velocity field in Zone 2. The equation is based on pitch to diameter ratio (p/DP).
V x , r 2 V m a x = E X P 1 2 r R m 0 2 σ 2
Parameter (r) determines the radial distance from the propeller axis. The author [20] measured standard deviation (σ) as a constant equal to 0.5 Rm0 up to the downstream distance of X/Dp = 0.5.
σ = 1 2 R m 0 f o r X D p < 0.5
Beyond X/Dp = 0.5 to the end of the zone of flow establishment, the standard deviation was defined as:
σ = 1 2 R m 0 + 0.075 X D p 2 f o r X D p > 0.5
The parameter (Rm0) is the radial distance of maximum velocity from the propeller shaft axis [21]. It is a function of the propeller radius (Rp) and propeller hub radius (Rh):
R m o = 0.67 R p R h
The maximum propeller jet velocity in Zone 3 was determined by [20], who performed experimental measurements based on the propeller geometry characteristics, and proposed the following equation:
V m a x V 0 = A X D p B
Parameters (A′) and (B′) are determined using the following equations:
A = 11.4 C t + 6.65 β + 2.16 P
B = C t 0.216 β 1.024 P 1.87
where the propeller blade area ratio is described with β.
The propeller jet velocity field (Vx,r_3) in Zone 3 was determined according to [17], with the next equation being
V x , r _ 3 V 0 = E X P 22.2 r X 2
The equations presented above are used to determine the velocity field of the propeller jet in the longitudinal (X) and vertical (r) directions corresponding to the propeller axis. The disadvantage of the presented semiempirical equations is that they do not take into account the distance between the propeller axis and the seabed, the bathymetry of the seabed, the deflection of the propeller jet by pier walls, the different depth zones in the ports and the shape of the harbour basins. A possible solution to these drawbacks is to recreate a ship’s manoeuvre with the FMBS and find a suitable method to calculate the propeller jet efflux velocity. There are several computer programs based on Computational Fluid Dynamics (CFD): Environmental Fluid Dynamics Code plus (EFDC+), which is an open source three-dimensional (3D) finite-difference surface water modelling program. The CFD is a branch of fluid mechanics that deals with data structures and numerical analysis to solve flow problems; such programs have already been used by authors to present fluid and sediment motion caused by ship propeller rotation [22,23,24,25,26,27]. Maynord’s model [28] and the Finite Analytical Navier-Stokes Solver (FANS model) [29] were interlinked with Curvilinear Hydrodynamics in three dimensions plus Tableau Input Coupled Kinetic Equilibrium Transport (CH3D + TICKET) [30]. It is a simulation tool to present contamination fate transport and redeposition of the sediment particles from the ship propeller jet. The CFD method is relatively inexpensive compared to experimental investigations. In this work, the hydrodynamic model MIKE 3 FM is coupled with MIKE 3 MT to represent the velocity distribution of the propeller jet corresponding to the bathymetry of the Koper Basin 1 harbour and mud transport after the propeller jet impedes the port bottom. The MIKE 3 FM is presented in the next subsection.

2.3. The Hydrodynamic Model MIKE 3

The three-dimensional, baroclinic model MIKE 3 FM is a commercial nonhydrostatic numerical modelling system designed for a wide range of applications in areas such as oceans, coastal regions, estuaries and lakes. It simulates unsteady three-dimensional flows, taking into account density variations, bathymetry and external influences such as metocean information, tides, currents and other hydrographic conditions. Features of relevance to this study’s hydrodynamic model include: bed resistance, density variations, turbulence modelling, isolated sources and sinks, connected source/sink pairs, pier resistance, particle tracking and discharge calculations [3].
The hydrodynamic model MIKE 3 FM in general is based on a continuity equation, using Cartesian co-ordinates.
u x + v y + w z = S
Parameters (u, v, w) present flow velocity components in the spatial domain (x, y, z).
Simulation MIKE 3 FM uses time integration of the shallow water and the transport equation of fluid/particles is performed using semi-implicit scheme. Using an explicit scheme, the time step interval must be selected, determining the Courant–Friedrich–Levy (CFL) number, which must be less than 1. The free surface is presented with the sigma-co-ordinate transformation approach. The parameters temperature of fluid and salinity are based on a general transport–diffusion equation [31].
Ship propeller jet velocity distribution is based on integral jet model equations, presenting steady state solutions of jet by solving a conservation equation for flux and momentum, temperature and salinity [32].

MIKE 3 Input Data

High-resolution bathymetry within the study area is established as a first step (Figure 5). The hydrodynamic model MIKE 3 FM includes subprogram MIKE Zero with additional features such as: Time Series and Mesh Generator, which are used to represent the sediment resuspension process in this work. The Mesh Generator is used to define the bathymetry; the data input was done using the cartesian co-ordinate system (x, y, z file). The bathymetry was defined as follows: maximum element area 8 m2, minimum allowed angle 26° and maximum number of nodes 55,000.
The second step is to specify bathymetry boundary conditions. Figure (Figure 6) shows the red zone which is specified as the fluid free outflow condition and green zone specifying land as containment zone (zero normal velocity).
The departure manoeuvre on FMBS was performed; ship motion data were recorded. The vessel (RPMapp) is used in Equation (4) to obtain (Papp), which was needed to calculate (Equation (5)) the propeller jet efflux velocity for zero ship speed, and later used in Equation (6) to obtain the propeller jet efflux velocity for nonzero ship speed. MIKE Zero contains time series’; the input data are: Discharge (m3/s), U-velocity component (m/s) and V-velocity component (m/s). The ship’s propeller is determined as a nozzle; the discharge rate (depending on the propeller diameter and efflux velocity when the ship velocity is non-zero), the U-velocity, V-velocity depend on the ship’s heading, the rudder angle and are calculated using the following equation.
D i s c h a r g e = A p r o p . V 0 ; v 0
U v e l . = s i n α h e a d ; r u d d . π 180 V 0 ; v 0
V v e l . = c o s α h e a d ; r u d d . π 180 V 0 ; v 0
Propeller surface (Aprop.) and angle (αhead;rudd) from geographic north, taking into account the ship’s heading and rudder angle (αhead;rudd). The results are inserted in MIKE 3 FM as a standard source (varying in time). The ship’s departure route was determined using 17 sources/locations, as shown in the following figure (Figure 7).
Each source is determined by Equations (19)–(21) in the variant time domain. The final source is the position of the ship when the simulation is stopped on FMBS (8.5 min into the ship’s manoeuvring).
The simulation took 30 min to find out how the velocity field of the propeller jet evolves after a ship manoeuvre due to the variable bathymetry and characteristics of the basin. The sediment structure in the Koper port was determined using the following parameters: density of sediment 2650 (kg/m3), bottom roughness 0.1 m and critical shear stress 0.1 (N/m2). The overall modelling methodology is described in figure below (Figure 8).

3. Results and Discussion

The results of the hydrodynamics and sediment transport modules are presented in this chapter. The numerical domain was vertically calculated at ten different depths. The kinematics of the vessel and the efflux velocities are illustrated in Figure 9. The most important parameters are propeller rates per minute (RPM), propeller efflux velocity (V0), propeller efflux velocity for a nonzero ship speed (V0;v≠0) and ship speed (SOG) in relation to time (8.5 min). The ship speed increases constantly from 0 to 6 m/s over 8.5 min. The departure manoeuvre starts with the command “Dead Slow Ahead”, resulting in the increased propeller flux velocity(V0), which is extremely invasive as the vessel is stationary, and as the vessel begins to move the efflux velocity slowly declines calculated by the nonzero ship speed (V0;v≠0). Initially, efflux velocity is 6 m/s with the discharge rate of 350 m3/s. When the ship’s speed is increasing, efflux velocity is decreasing (black line).
After about 5 min of manoeuvring, the command “slow ahead” is given; a significant increase of (V0) 8.5 m/s and a decrease of (V0;v≠0) 4.7 m/s are observed. The latter is the main parameter influencing mud steering-up process. At the beginning of the manoeuvre, the propeller slip is considerable; resulting in a high (V0;v≠0), with increasing (SOG) and steady (RPM). The conclusion is that the decrease in (V0;v≠0) depends mainly on the relationship between (SOG) and (RPM).
Previous research regarding bottom wash assessment in the Port of Koper has been presented by [33]. In this study, the authors determined the critical propeller jet bottom velocity of 0.35 m/s, at which the jet velocity overcomes the friction of the mud bottom and causes the detachment of sediment particles and their circulation from the seafloor. Figure 10 shows the jet velocity (m/s) generated by the ship’s propeller in bottom layer 1. Figure 11 shows the resulting suspended sediment concentration (kg/m3) in the same layer. The caption is from ship positions 1 and 2, respectively (Figure 7).
Ship propeller generated jet in bottom layer (m/s) Figure 10 shows the corresponding velocity fields of the propeller jet at 134 s when the sediment particles coloured green, yellow and red detach from the bottom. The vectors show the direction and velocity of the propeller flow. The flow is deflected to the right about 100 m ahead of the propeller front due to the starboard rudder angle (STBD) of 11 degrees, the proximity of the pier wall, the horseshoe-shaped port depression and the decrease in port depth from 15 m to 13 m at position long. 401,300. A nearly negligible contribution to the current deflection is due to the Coriolis force, which forces the fluids in the northern hemisphere to turn in the clockwise direction.
Figure 11 shows the total specific sediment concentration (SSC) in (kg/m3) in bottom layer 1; during investigations, it can reach a maximum value of 0.07 (kg/m3). The resuspended sediment particles are further advected by the flow field induced by the propeller jet (Figure 10).
The propeller jet-induced flow and total (SSC) are shown in Figure 12. The duration of the manoeuvre is 510 s. The figures were taken after 1000 s; they show the flow velocity/path and the total (SSC) in bottom layer 1.
As expected, there is a significant decrease in flow velocity (by a factor of three: see Figure 12). From the colour palette, the ship’s departure route can be clearly seen. The flow is deflected clockwise due to the horseshoe-shaped basin. The south-eastern part of the basin is subject to significant erosion of the sediment material, which is influenced by the high velocity flow field.
The low velocity flow eddy surrounded by high velocity fields causes sediment material deposition, marked on Figure 12. The second part of Figure 12 shows the SSC in the same bottom layer and the time step. The colour palette shows the content of resuspended material. The largest amount of SSC is expected at the position shown in the first part of the figure (sediment deposition). Sediment accumulation is expected at the same position.
Validation of the MIKE 3 MT simulation program was performed using a satellite image of Basin 1 from a similar vessel. Figure 13 shows the SSC at the surface (layer 10). It can be seen that the largest sediment concentration is formed in a vortex position at the same location as in the real ship scenario (Google Earth image).
Validation of results is also based on measurements performed 25 July 2013 in the bay of Koper at 38 measurement locations after the departure manoeuvre of the container ship Ital Universo (IMO no.: 9196993). The determination tool of total resuspended particles (SSC) has been water samples, taken at three depths (0.5, 6.0 and 12.0 m). Measurements have been performed after 3 min., when the tugs’ lines have been released. The paper [34] considers measurements from the manoeuvre zone (Figure 14), where it was displayed as a maximum value of (TSS~0.139 kg/m3) and sampled a concentration of (SSC~0.037 kg/m3) in the bottom layer.
Comparing the numerical modelling with the experimental investigations above, there was clearly impressive similarity in the results. Both modelling and sampling were performed in the bottom layer (hydrodynamic model at 1.0 m; and in situ sampling at 0.5 m). The hydrodynamic model showed (Figure 11) noticeable (SSC~0.07 kg/m3) after 2 min 14 s from the start of the manoeuvre. The container ship was first towed astern from berth and after being turned, the vessel accelerated with power similar to that in our case. Their measurements of total resuspended particles were similar (SSC~0.037 kg/m3).
The difference in (SSC) between the real case and the modelled case is mainly due to the greater depth of the former and the presence of two tugs.
The significance here is that now that we are reassured that our model successfully mimics actual circumstances, we can proceed to examine any circumstances, any ship size, UKC, propeller type, and so on, and determine the optimal manoeuvring approach. We are confident as well that adding tugs and thrusters will not have the slightest deleterious effect on the modelling. Perhaps it should also be mentioned that as was implied early on, the effect of an abstract (having nothing to do with nature) decision-making motive and mindset is liable to continue to have a degrading effect on port areas and beyond, and we will have a relatively easy means of displaying these effects. Of course, some protective measures have been introduced, such as the rip-rap technique, and we will be able to monitor their effectiveness [35,36]. Based on our calculation, we can recommend a rip-rap protection area in the northern part of Basin 1 along the entire container terminal and 40 m in lateral distance (Figure 15).

4. Conclusions

The impact of ship traffic on the seabed of Koper Harbour Basin I was analysed using hydrodynamic numerical models based on a high-resolution non-hydrostatic circulation model (MIKE 3 HD FM) with a sediment transport model (MIKE 3 MT FM). The simulation results show that the ship propeller jet can reach velocities of more than 1.5 m/s. The intensity of the ship propeller jet depends on the following factors: propeller dimensions and RPM, distance between the propeller and bottom and whether the manoeuvre is a departure or an arrival. The critical jet bottom velocity at which sediment particles detach from the seafloor is reported to be 0.35 m/s. The jet zone of the propeller can reach distances up to 30–40 times the propeller diameter in the horizontal direction. The simulation allowed validation of the sediment dynamics in the harbour using satellite images of the sea surface showing the generated sediment eddy (Figure 13). Further simulation analysis shows why certain port areas are subject to sediment erosion and others to deposition. A nonhydrostatic circulation model (MIKE 3 HD FM and MT) was tested to determine its ability to represent ship-induced propeller-generated jets. It was determined that the tested programme is suitable for further research, replicating alternative ship manoeuvres using the Full Mission Bridge Simulator and later hydrostatic ship motion data implemented in the MIKE 3 HD FM fluid motion program.
The limitation of this paper aligns well with our intent: now that our model is utile, we intend to use this model to conduct a number of applications and find the optimal means of manoeuvring in the port under all variety of circumstances.

Author Contributions

Conceptualization, Principles and methods, Conducting surveys, Data analysis, Software and Conclusions done by J.S., Principles and methods and Conducting surveys, Validation done additionally by M.P. and A.G. All authors have read and agreed to the published version of the manuscript.

Funding

The publication of the paper is partially financed by the research group at the Faculty of Maritime Studies and Transport, financed by the Slovenian National Research Agency.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data sets used and/or analyzed in this study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors acknowledge the Port of Koper, the Slovenian Maritime Administration, marine pilot Matjaž Felicjan, and Antonio Guarnieri.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Srše, J.; Perkovič, M.; Grm, A.; Androjna, A.; Brcko, T. Sediment resuspension evaluation on two manoeuvres and two different pilots. In Proceedings of the ICTS 2022: 20. Mednarodno Posvetovanje o Prometni Znanosti, Portorož, Slovenia, 23–24 May 2022; Faculty of Maritime Studies and Transport: Portorose, Slovenia, 2022; pp. 332–336. [Google Scholar]
  2. Širca, A.; Rajar, R.; Harris, R.C.; Horvat, M. Mercury transport and fate in the Gulf of Trieste (Northern Adriatic)—A two-dimensional modelling approach. Environ. Model. Softw. 1999, 14, 645–655. [Google Scholar] [CrossRef]
  3. DHI. MIKE 3 Flow Model FM, Hydrodynamic Module, User Guide, 2021; DHI: Singapore, 2021. [Google Scholar]
  4. Hamill, G.A.; Kee, C. Predicting axial velocity profiles within a diffusing marine propeller jet. Ocean. Eng. 2016, 124, 104–112. [Google Scholar] [CrossRef]
  5. Hayes, D.F.; Chintamaneni, R.; Bommareddy, P.; Cherukuri, B. Propwash Impacts on Water Quality Around Dredging and Other Marine Construction Activities. In Proceedings of the 30th Annual WEDA Conference, San Juan, PR, USA, 6–9 June 2010. [Google Scholar]
  6. Yew, W.T. Seabed Scour Induced by Twin-Propeller Ships/Yew Wan Tian. Ph.D. Thesis, University of Malaya, Kuala Lumpur, Malaysia, 2017. [Google Scholar]
  7. Penna, N.; D’Alessandro, F.; Gaudio, R.; Tomasicchio, G. Three-dimensional analysis of local scouring induced by a rotating ship propeller. Ocean Eng. 2019, 188, 106294. [Google Scholar] [CrossRef]
  8. Wei, M.; Chiew, Y.-M. Characteristics of Propeller Jet Flow within Developing Scour Holes around an Open Quay. J. Hydraul. Eng. 2018, 144, 04018040. [Google Scholar] [CrossRef]
  9. Svanberg, M.; Santén, V.; Hörteborn, A.; Holm, H.; Finnsgård, C. AIS in maritime research. Mar. Policy 2019, 106, 103520. [Google Scholar] [CrossRef]
  10. Transas. Description of Transas Ship Motion Mathematical Modeling; Transas: Singapore, 2016. [Google Scholar]
  11. Roh, M.-I.; Lee, K.-Y. Computational Ship Design; Springer: Singapore, 2018. [Google Scholar]
  12. Blaauw, H.G.; Van de Kaa, E.J. Erosion of Bottom and Sloping Banks Caused by the Screw Race of Manoeuvring Ships; Delft Hydraulics Lab: Delft, the Netherlands, 1978; Volume 202, pp. 1–12. [Google Scholar]
  13. International Navigation, A. Guidelines for Protecting Berthing Structures from Scour Caused by Ships; PIANC: Bruxelles, Belgique, 2015. [Google Scholar]
  14. Verhey, H.J. The stability of bottom and banks subjected to the velocities in the propeller jet behind ships. In Proceedings of the 8th International Harbour Congress, Antwerp, Belgium, 13–17 June 1983; Voume 303, pp. 1–11. [Google Scholar]
  15. Stewart, D.P.J. Characteristics of a Ship’s Screw Wash and the Influence of Quay Wall Proximity; Queen’s University of Belfast: Belfast, UK, 1992. [Google Scholar]
  16. Paper, W. EFDC+ Propeller Wash Module White Paper; DSI, LLC Edmonds: Edmonds, WA, USA, 2021. [Google Scholar]
  17. Fuehrer, M.; Römisch, K. Effects of modern ship traffic on inland and ocean waterways and their structures. In Proceedings of the 24th International Navigation Congress, Leningrad, Russia, September 1977; Section I-3. [Google Scholar]
  18. Maynord, S.T. Physical Forces Near Commercial Tows; US Army Corps of Engineers: Vicksburg, MS, USA, 2000. [Google Scholar]
  19. Albertson, M.L.P.; Dai, Y.B.; Jensen, R.A.; Rouse, H. Diffusion of Submerged Jets. Trans. Am. Soc. Civ. Eng. 1948, 115, 1571–1596. [Google Scholar] [CrossRef]
  20. Hamill, G.A. The Scouring Action of the Propeller Jet Produced by a Slowly Manoeuvring Ship. Bull. Perm. Int. Assoc. Navig. Congr. 1988, 85–110. [Google Scholar]
  21. Bergh, H.; Cederwall, K. Propeller Erosion in Harbours; Hydraulics Laboratory, Royal Institute of Technology Stockholm: Stockholm, Sweden, 1981. [Google Scholar]
  22. Guarnieri, A.; Saremi, S.; Pedroncini, A.; Jensen, J.H.; Torretta, S.; Vaccari, M.; Vincenzi, C. Effects of marine traffic on sediment erosion and accumulation in ports: A new model-based methodology. Ocean Sci. 2021, 17, 411–430. [Google Scholar] [CrossRef]
  23. Lam, W.H.; Hamill, G.A.; Robinson, D.J. Initial wash profiles from a ship propeller using: CFD method. Ocean Eng. 2013, 72, 257–266. [Google Scholar] [CrossRef]
  24. Muscari, R.; Dubbioso, G.; Viviani, M.; Di Mascio, A. Analysis of the asymmetric behavior of propeller–rudder system of twin screw ships by CFD. Ocean Eng. 2017, 143, 269–281. [Google Scholar] [CrossRef]
  25. Craig, P.; Jung, J.; Mausolff, Z.; Bastidas, L.; Mathis, T.; Wang, P.-F. Modeling Sediment Resuspension and Transport Processes Induced by Propeller Wash from Ship Traffic. J. Hydraul. Eng. 2023, 149, 04023009. [Google Scholar] [CrossRef]
  26. Gavrilin, S. Validation of Ship Manoeuvring Simulation Models; ITTC: Zürich, Switzerland, 2018. [Google Scholar]
  27. Absalonsen, L. Pacific Northwest LNG; LNG Jetty; Propeller Scour Analysis; Hatch: Los Angeles, CA, USA, 2014. [Google Scholar]
  28. Maynord, S.T. Riprap Protection on Navigable Waterways [Microform]/by Stephen T. Maynord; Technical report; HL-84-3; Hydraulics, L., Ed.; U.S. Army Engineer Waterways Experiment Station: Vicksburg, MS, USA, 1984. [Google Scholar]
  29. Chen, H.-C.; Huang, E.T. Time-Domain Simulation of Floating Pier/Ship Interactions And Harbor Resonance. In Proceedings of the Fourteenth International Offshore and Polar Engineering Conference, Toulon, France, May 2004. [Google Scholar]
  30. Wang, P.-F.; Duarte, I.R.; Richter, K. Evaluation of Resuspension from Propeller Wash in DoD Harbors; ER-201031; Environmental Security Technology Certification Program: San Diego, CA, USA, 2016. [Google Scholar]
  31. DHI. MIKE 21 & 3 Flow Model FM, Hydrodynamic Module, Short Description, 2013; DHI: Denmark, 2013. [Google Scholar]
  32. Jirka, G.H. Integral Model for Turbulent Buoyant Jets in Unbounded Stratified Flows. Part I: Single Round Jet. Environ. Fluid Mech. 2004, 4, 1–56. [Google Scholar] [CrossRef]
  33. Perkovic, M.; Batista, M.; Malacic, V.; Zagar, D.; Jankowski, S.; Gucma, L.; Gucma, M.; Rostopshin, D. Bottom Wash Assessment Using a Full Mission Ship Handling Simulator. In Proceedings of the MARSIM 2012, Singapore, 23–27 April 2012. [Google Scholar]
  34. Soczka Mandac, R.; Covelli, S.; Faganeli, J. Sediment resuspension during vessel manoeuvres in port areas: Evidence from field observations. J. Soils Sediments 2023. [Google Scholar] [CrossRef]
  35. Belyaev, N.D.; Klimovich, V.I. State of the Art in Protection of Bottom Against Scour Caused by Marine Propulsors. Power Technol. Eng. 2019, 53, 159–166. [Google Scholar] [CrossRef]
  36. CIRIA; CUR; CETMEF. Rock Manual: The Use of Rock in Hydraulic Engineer; CIRIA: London, UK, 2007. [Google Scholar]
  37. Srše, J.; Perkovič, M.; Androjna, A.; Brcko, T. Solutions to minimize sediment resuspension in ports. Časopis Pomor. Fak. Kotor J. Marit. Sci. 2022, 23, 47–56. [Google Scholar] [CrossRef]
Figure 1. Black line presents studied area; Basin I of the harbour of Koper, Slovenia [Google Earth image].
Figure 1. Black line presents studied area; Basin I of the harbour of Koper, Slovenia [Google Earth image].
Jmse 11 01619 g001
Figure 2. Study flow chart for the acquisition of dynamic data of a container ship departure manoeuvre and numerical modelling.
Figure 2. Study flow chart for the acquisition of dynamic data of a container ship departure manoeuvre and numerical modelling.
Jmse 11 01619 g002
Figure 3. The investigated area shows the position of the outgoing container ship with a bathymetric depth of 15 m.
Figure 3. The investigated area shows the position of the outgoing container ship with a bathymetric depth of 15 m.
Jmse 11 01619 g003
Figure 5. MIKE Zero Mesh Generator; Port of Koper, Basin I.
Figure 5. MIKE Zero Mesh Generator; Port of Koper, Basin I.
Jmse 11 01619 g005
Figure 6. Bathymetry boundary conditions specification.
Figure 6. Bathymetry boundary conditions specification.
Jmse 11 01619 g006
Figure 7. Propeller jet source determined by 17 locations.
Figure 7. Propeller jet source determined by 17 locations.
Jmse 11 01619 g007
Figure 8. Research methodology review.
Figure 8. Research methodology review.
Jmse 11 01619 g008
Figure 9. Ship propulsion jet velocity characteristics.
Figure 9. Ship propulsion jet velocity characteristics.
Jmse 11 01619 g009
Figure 10. Ship propeller generated jet in bottom layer (m/s).
Figure 10. Ship propeller generated jet in bottom layer (m/s).
Jmse 11 01619 g010
Figure 11. Resulting suspended sediment concentration (kg/m3).
Figure 11. Resulting suspended sediment concentration (kg/m3).
Jmse 11 01619 g011
Figure 12. Current velocity/path and total SSC in bottom layer 1.
Figure 12. Current velocity/path and total SSC in bottom layer 1.
Jmse 11 01619 g012
Figure 13. Sediment concentration (SSC) in surface layer and validation with satellite picture.
Figure 13. Sediment concentration (SSC) in surface layer and validation with satellite picture.
Jmse 11 01619 g013
Figure 14. Bay of Koper. The yellow marks presenting 38 measurements locations and purple line presents departure route from container terminal Basin 1, letters A, B, C presenting prior measurement to verify content of resuspended particles prior to the vessel manoeuvre [34].
Figure 14. Bay of Koper. The yellow marks presenting 38 measurements locations and purple line presents departure route from container terminal Basin 1, letters A, B, C presenting prior measurement to verify content of resuspended particles prior to the vessel manoeuvre [34].
Jmse 11 01619 g014
Figure 15. Port of Koper Basin 1; recommended bottom protection area [37].
Figure 15. Port of Koper Basin 1; recommended bottom protection area [37].
Jmse 11 01619 g015
Table 1. Main particulars of subject vessel.
Table 1. Main particulars of subject vessel.
Ship ParticularsValue/Type
Length Overall (LOA)203.6 m
Moulded Breadth (B)25.4 m
Summer Draft (Ts)9.8 m
Displacement (D)32,025 t
Max engine power (P)15,890 kW
Ship propeller typeFixed Pitch Propeller (FPP)
Propeller diameter (Dp)6.0
Propeller immersion6.67 m
Bow thruster capacity763 kW
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Srše, J.; Perkovič, M.; Grm, A. Sediment Resuspension Distribution Modelling Using a Ship Handling Simulation along with the MIKE 3 Application. J. Mar. Sci. Eng. 2023, 11, 1619. https://doi.org/10.3390/jmse11081619

AMA Style

Srše J, Perkovič M, Grm A. Sediment Resuspension Distribution Modelling Using a Ship Handling Simulation along with the MIKE 3 Application. Journal of Marine Science and Engineering. 2023; 11(8):1619. https://doi.org/10.3390/jmse11081619

Chicago/Turabian Style

Srše, Jure, Marko Perkovič, and Aleksander Grm. 2023. "Sediment Resuspension Distribution Modelling Using a Ship Handling Simulation along with the MIKE 3 Application" Journal of Marine Science and Engineering 11, no. 8: 1619. https://doi.org/10.3390/jmse11081619

APA Style

Srše, J., Perkovič, M., & Grm, A. (2023). Sediment Resuspension Distribution Modelling Using a Ship Handling Simulation along with the MIKE 3 Application. Journal of Marine Science and Engineering, 11(8), 1619. https://doi.org/10.3390/jmse11081619

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop