This paper examines the costs and benefits of reduction measures for the shipping industry to comply with the forthcoming sulphur emission regulations. Sulphur scrubbers and marine gas oil are two promising alternatives for ship owners. However, their economic comparisons are primarily based on a private perspective. This paper provides a wider viewpoint by integrating the private abatement costs of ship owners and the social environmental benefits from emission reduction. The results showed that the price spread between marine gas oil and heavy fuel oil is a determining factor in making this choice. Marine gas oil tends to have higher net present values than scrubbers when the price spread of fuel is less than 231 Euros per tonne. Furthermore, it is more beneficial to install a scrubber on new ships than retrofits. An old ship is not suitable for a scrubber installation when its remaining lifespan is less than 4 years.
In tramp shipping, a preliminary route is required for voyage planning at the pre-fixture stage (before a chartering contract is agreed). Such routes are conventionally designated by using pilot charts or software considering long-term statistical weather. However, it has been experienced by tramp operators that such route solutions often poorly estimated sailing distances for long journeys and thereby cause inappropriate cost estimation and bad voyage plan. To fill this gap, a data-driven methodology is proposed in this paper to establish a practical route library with the consideration of ship sizes, load conditions and seasonality. In this method, it first requires a dividing of ship trajectories into local sea passage and open sea passage. The voyage trajectories made of AIS points are then simplified to pattern nodes based on a speed-weighted geolocation method. Afterwards, the KMeans algorithm is deployed to properly classify these pattern nodes, identifying the most representative nodes (routes) in open sea passages. Simultaneously, the connection points are identified by DBSCAN algorithm, representing local sea passages. Combining the representative routes in open sea passages and the connection points in local sea passages, the most navigated routes between two ports are obtained. Finally, case studies are conducted for the Pacific Ocean and the Atlantic Ocean respectively using global AIS data from tanker vessels to demonstrate the feasibility and effectiveness of this methodology. The proposed route library is capable of providing reliable route references to support the decision-making at the pre-fixture stage.
Sensing data from vessel operations are of great importance in reflecting operational performance and facilitating proper decision-making. In this paper, statistical analyses of vessel operational data are first conducted to compare manual noon reports and autolog data from sensors. Then, new indicators to identify data aberrations are proposed, which are the errors between the reported values from operational data and the expected values of different parameters based on baseline models and relevant sailing conditions. A method to detect aberrations based on the new indicators in terms of the reported power is then investigated, as there are two independent measured power values. In this method, a sliding window that moves forward along time is implemented, and the coefficient of variation (CV) is calculated for comparison. Case studies are carried out to detect aberrations in autolog and noon data from a commercial vessel using the new indicator. An analysis to explore the source of the deviation is also conducted, aiming to find the most reliable value in operations. The method is shown to be effective for practical use in detecting aberrations, having been initially tested on both autolog and noon report from four different commercial vessels in 14 vessel years. Approximately one triggered period per vessel per year with a conclusive deviation source is diagnosed by the proposed method. The investigation of this research will facilitate a better evaluation of operational performance, which is beneficial to both the vessel operators and crew.
Ship air emissions are recognized as one of the key concerns of the maritime industry. Competent authorities have issued various regulations to manage air emissions from ships. Although the authorities are policy makers, the effectiveness of policies is up to the shipping industry who operates the vessels and terminals to fulfill maritime transportation works. Given this characteristic, bi-level optimization model has been widely adopted in studies that optimize policy design or evaluate its effectiveness. The framework of a typical bi-level optimization model for ship emission management problem is given to show the basic structure of similar issues. A series of applications of bi-level optimization model in managing ship emissions is reviewed, including cases of Energy Efficiency Design Index, Emissions Control Area, Market Based Measure, Carbon Intensity Indicator, and Vessel Speed Reduction Incentive Program. We hope this paper can enlighten scholars interested in this area and provide help for them.
A moving static pressure distribution is commonly used to simulate a travelling ship. However, the ship movement changes the fluid velocity around the hull, inducing pressures on the hull surface that are no longer equal to the static pressure. Therefore, we introduce a dynamic pressure correction strategy, which can accurately simulate the impact of the ship movement on the hull-surface pressure and preserve the desired hull shape under both stationary and transient conditions. The strategy is applied to a high-order spectral model and used to investigate ship-induced waves and wave resistance over a both flat and variable topography. We explore various parameters in our study, including the average water depth to ship draft ratio (h(0)/d), the channel width to ship width ratio (W/B), the Froude number (Fr-0 = U/root gh(0)) and variations in bathymetric slope. Compared with experiments on a flat bottom, the numerical results with dynamic correction show better accuracy in the simulation of ship-induced waves and wave resistance than those obtained using a static pressure distribution. The correlation coefficient for wake waves between the numerical and experimental results is improved by approximately 0.25 with the dynamic correction strategy. The amplitude and wavelength of ship-induced mini-tsunamis over a variable topography are found to be reduced when employing a dynamic correction compared with a static pressure distribution, and this effect becomes more pronounced with higher Froude number. The static pressure approach is shown to allow large deformations of the desired hull shape and changes in ship volume which are responsible for the different wave patterns from the two approaches.
The levelized costs of energy (LCoE) of wave power is still not fully competitive with other sources of renewable energy. However, wave energy is partly in a different phase than other renewable energy types and could thus contribute to a better predictability and smoothed power output. This work focuses on co-location of wave and wind power by investigating the intermittency of wind and waves power based on measured historical data from several hundreds of locations worldwide. Employing wind power curves and wave power matrices, the sites are evaluated based on several different metrics. The results indicate that there are several spots where wave power has a much lower intermittency than wind power providing reliable energy supply. Best sites for co-location in terms of energy yield were found in North-Western Europe. However, both wind and wave production have the same seasonal variability in these sites. Only a handful of sites found in California showed the possibility of seasonal power smoothing using the combination of wind and wave.
Mooring systems for floating wave energy converters often rely on floaters to allow for minimum restraints of the body motion in heavy. However, the inclusion of floaters also introduce possible slack-taut scenarios induced by the dynamic response of the floater in relation to the fair-lead point of the mooring. This can increase the occurrence of snap loads. The present study outlines the work to include floaters and sinks into a high-order discontinuous Galerkin model for mooring cable dynamics. Numerical simulations of a mooring leg adapted from the Waves4Power full-scale device are performed, and the results from varying the floater geometry are analyzed.
For this case the floater influence on the occurrence of snap loads was clearly evident. There is a strong correlation between floater pitch response and cable slack in the upper mooring cable. For a floater with constant buoyancy, increasing the floater height and thereby increasing the pitch inertia of the floater is shown to decrease the range of frequencies where cable slack occurs. It is illustrated that for some cases, changing floater geometry can avoid slack altogether. A careful design of the floater geometry can thus make a large difference for the dynamic load factor of the mooring system.
This work presents the verification and validation of the freely available simulation tool MoodyMarine, developed to help meet some of the demands for early stage development of MRE devices. MoodyMarine extends the previously released mooring module MoodyCore (Discontinuous Galerkin Finite Elements) with linear radiation-diffraction bodies, integrated pre-processing workflows and a graphical user interface. It is a C++ implementation of finite element mooring dynamics and Cummins equations for floating bodies with weak nonlinear corrections. A newly developed nonlinear Froude-Krylov implementation is verified in the paper, and MoodyMarine is compared to CFD simulations for two complex structures: a slack-moored floating offshore wind turbine and a self-reacting point-absorber with hybrid mooring.
Marine cables are primarily designed to support axial loads. The effect of bending stiffness on the cable response is therefore often neglected in numerical analysis. However, in low-tension applications such as umbilical modeling of ROVs or during slack events, the bending forces may affect the slack regime dynamics of the cable. In this paper, we present the implementation of bending stiffness as a rotation-free, nested local Discontinuous Galerkin (DG) method into an existing Lax–Friedrichs-type solver for cable dynamics based on an hp-adaptive DG method. Numerical verification shows exponential convergence of order P and P + 1 for odd and even polynomial orders, respectively. Validation of a swinging cable shows good comparison with experimental data, and the importance of bending stiffness is demonstrated. Snap load events in a deep water tether are compared with field-test data. The bending forces affect the low-tension response for shorter lengths of tether (200–500 m), which results in an increasing snap load magnitude for increasing bending stiffness. It is shown that the nested LDG method works well for computing bending effects in marine cables.
The power output from many wave energy converters (WECs) is limited by a finite stroke length in the power take-off (PTO) mechanism. As the PTO approaches its maximum stroke length, an end-stop system needs to be engaged to avoid damage to the machinery. Still the on-set of the end-stop is a nonlinear trigger force, a stiff point in the system. In this respect it is similar to how snap loads in the mooring cables affect the system after a period of cable slack. This paper presents a detailed study into the dynamics of end-stop events and snap loads for a WEC. The WEC is a bottom-mounted linear generator connected to a surface buoy via a steel wire. By comparing a linear spring model with three dynamic mooring line models we conclude that large differences are observed in the low-tension and slack regions of the cable during moderate wave loads, while minor differences are seen in the estimated peak tension. By further varying end-stop parameters we observe that the peak tension in the line changes mildly with the axial stiffness for moderate wave heights. The peak tension is surprisingly unaffected by the introduction of a critical damping level to the end-stop system, despite the significant increase in end-stop force which causes the translator to come to a sudden stop. We discuss how the connection between maximum line force and end-stop parameters is highly dependent on the buoy position in the wave at the instant of end-stop onset.