PMRes 2023 4
PMRes 2023 4
PMRes 2023 4
of
Naval architecture
No 4 (120) 2023
Vol. 30
ADDRESS OF PUBLISHER
& EDITOR’S OFFICE:
GDAŃSK UNIVERSITY
OF TECHNOLOGY
Institute
of Naval Architecture 4 Piotr Gomułka
G. Narutowicza 11/12 Computer-aided system for layout of fire hydrants on boards
80-233 Gdańsk, POLAND designed vessel using the Particle Swarm Optimization algorithm
17 Yordan Denev
RETROFITTING THE BOW OF A GENERAL CARGO VESSEL ANDEVALUATING
ENERGY EFFICIENCY OPERATIONAL INDEX
24 Jan Sierzputowski, Artur Karczewski, Przemysław Krata
USE OF THE AHP METHOD FOR PREFERENCE DETERMINATION IN YACHT
EDITORIAL STAFF: DESIGN
31 Nguyen Dinh Thach, Phan Van Hung
Wiesław Tarełko OPTIMAL UV QUANTITY FOR A BALLAST WATER TREATMENT SYSTEM FOR
| Editor in Chief COMPLIANCE WITH IMO STANDARDS
Janusz Kozak 43 Se-yeol An, Hyeon-won Jeong, Ohyoung Kim, W. Jaewoo Shim
| Deputy Editors-in-Chief EFFECTS OF SWAY AND ROLL EXCITATIONS ON SLOSHING LOADS IN A KC-1
Wojciech Litwin MEMBRANE LNG TANK
| Deputy Editors-in-Chief 54 Jacek Nakielski
ANALYSIS OF THE ENVIRONMENTAL IMPACT OF THE HULL CONSTRUCTION
OF A SMALL VESSEL BASED ON LCA
61 Ali Zinati, Mohammad Javad Ketabdari, Hamid Zeraatgar
EFFECTS OF PROPELLER FOULING ON THE HYDRODYNAMIC PERFORMANCE
Price: OF A MARINE PROPELLER
single issue: 25 PLN 74 Volodimir Leshchev, Igor Maslov, Oleksandr Palagin, Andrii Naydyonov
TRANSFER FUNCTION FOR A CONTROLLABLE PITCH PROPELLER WITH ADDED
Prices for abroad WATER MASS
single issue: 81 Dmytro Minchev, Roman Varbanets, Oleksandr Shumylo, Vitalii Zalozh,
- in Europe EURO 15 Nadiia Aleksandrovska, Pavlo Bratchenko, Thanh Hai Truong
- overseas USD 20 Digital Twin test-bench Performance for marine diesel engine
applications
WEB: 92 Gang Wu, Guodong Jiang, Changsheng Chen, Guohe Jiang, Xigang Pu,
pg.edu.pl/pmr Biwen Chen
AN EXPERIMENTAL STUDY OF THE EFFECTS OF CYLINDER LUBRICATING
e-mail : pmr@pg.edu.pl OILS ON THE VIBRATION CHARACTERISTICS OF A TWO-STROKE LOW-SPEED
MARINE DIESEL ENGINE
ISSN 1233-2585
ADDRESS OF PUBLISHER 102 Volodymyr Korobko, Serhiy Serbin, Huu Cuong Le
& EDITOR’S OFFICE: EXPLORATION OF A MODEL THERMOACOUSTIC TURBOGENERATOR WITH
A BIDIRECTIONAL TURBINE
GDAŃSK UNIVERSITY 110 Jerzy Girtler, Jacek Rudnicki
OF TECHNOLOGY QUANTUMNESS IN DIAGNOSTICS OF MARINE INTERNAL COMBUSTION
ENGINES AND OTHER SHIP POWER PLANT MACHINES
Institute 120 Piotr Bilon, Wojciech Litwin
of Naval Architecture HUMAN RESOURCE MANAGEMENT DIGITALISATION IN MULTIDISCIPLINARY
G. Narutowicza 11/12 SHIP DESIGN COMPANIES
80-233 Gdańsk, POLAND 129 Minh Duc Nguyen, Ko Tae Yeon, Krzysztof Rudzki, Hoang Phuong Nguyen,
Nguyen Dang Khoa Pham
STRATEGIES FOR DEVELOPING LOGISTICS CENTRES: TECHNOLOGICAL
TRENDS AND POLICY IMPLICATIONS
Editorial
POLISH MARITIME RESEARCH is the scientific journal with a worldwide circulation. This journal is published
quarterly (four times a year) by Gdansk University of Technology (GUT). On September, 1994, the first issue of
POLISH MARITIME RESEARCH was published. The main objective of this journal is to present original research,
innovative scientific ideas, and significant findings and application in the field of :
Naval Architecture, Ocean Engineering and Underwater Technology,
The scope of the journal covers selected issues related to all phases of product lifecycle and corresponding technologies
for offshore floating and fixed structures and their components.
All researchers are invited to submit their original papers for peer review and publications related to methods of the
design; production and manufacturing; maintenance and operational processes of such technical items as:
• all types of vessels and their equipment,
• fixed and floating offshore units and their components,
• autonomous underwater vehicle (AUV) and remotely operated vehicle (ROV).
We welcome submissions from these fields in the following technical topics:
• ship hydrodynamics: buoyancy and stability; ship resistance and propulsion, etc.,
• structural integrity of ship and offshore unit structures: materials; welding; fatigue and fracture, etc.,
• marine equipment: ship and offshore unit power plants: overboarding equipment; etc.
Scientific Board
Chairman : Prof. JERZY GIRTLER - Gdańsk University of Technology, Poland
Vice-chairman : Prof. CARLOS GUEDES SOARES, Universidade de Lisboa, Lisbon, Portugal
Vice-chairman : † Prof. MIROSŁAW L. WYSZYŃSKI - University of Birmingham, United Kingdom
Piotr Gomułka
1
Damen Engineering Gdansk, Gdansk
2
Gdansk University of Technology, Gdansk, Poland
Abstract
The functional layout of fire safety equipment in technical spaces of ships is a time-consuming process. When designing
a ship fire protection system, the designer must manually position each system component in such a way as to meet the
requirements of regulations arising from the technical specification, various legal regulations of maritime conventions
and classification societies of the vessel to be designed. Layout of fire hydrants assisted by a computer that is based
on pre-defined criteria and various constraints could significantly support the designer in working easier and faster.
This paper presents a prototype computer-aided design system that enables optimal placement of fire hydrants using
the metaheuristic Particle Swarm Optimization (PSO) algorithm. This algorithm was used in Rhinoceros 3D software
with its Grasshopper plugin for visualizing the arrangement of fire safety equipment. Various solution arrangements
compared with the fire hydrant placement in real ships are illustrated by a case study. Demonstrating how design work
can be facilitated and what potential benefits can be achieved are presented as well.
𝑛𝑛
∑𝑖𝑖=0 𝑤𝑤𝑖𝑖 ∙𝑝𝑝𝑖𝑖 ∙𝐵𝐵𝑖𝑖𝑛𝑛 (𝑡𝑡)
𝑝𝑝(𝑡𝑡) = 𝑛𝑛 𝑓𝑓𝑓𝑓𝑓𝑓 𝑡𝑡 𝜖𝜖 (0,1) (1) (1)
∑𝑖𝑖=0 𝑤𝑤𝑖𝑖 ∙𝐵𝐵𝑖𝑖𝑛𝑛 (𝑡𝑡)
Significance of compartments electrical equipment, etc. – fire hoses with a length of not
The SOLAS Convention has distinguished different less than 10 [m] and not more than 15 [m] must be used,
classes of bulkheads based on their ability to withstand –– areas on the ship where the probability of fire is slight
fire and smoke on the endangered side. These bulkheads are usually considered places where the risk of fire is
are classified as A, B, and C classes and create thermal and small and negligible, usually with limited human such
structural boundaries on the vessel. Spaces are separated by as voids, cofferdams, tanks, chain lockers, fixed gas fire-
the appropriate class of bulkhead (either bulkheads or decks) extinguishing system storage rooms, and others – in this
depending on the level of fire risk. Details specifying the fire given case, there is no need to install a fire hydrant,
sealing standards to be applied to specific bulkheads between –– other spaces and open decks – fire hoses with a length of
neighboring compartments can be found in SOLAS II-2 Part not less than 10 [m] and not more than 20 [m] must be used.
C. The following types of spaces are most relevant in terms In the developed system, each type of space is assigned to
of fire hydrant distribution: a specific drawing layer (Fig. 3). The number of layers and
–– machinery space category A, which is classically considered their assigned properties can be changed as desired. For
the area with the highest fire risk and requires the highest example, machinery spaces of category A can be assigned
safety standards – fire hoses with a length of not less than to layer 1, while other machinery spaces can be assigned to
10 [m] and not more than 15 [m] must be used, layer 2. Spaces with low fire risk, such as voids or cofferdams,
–– machinery spaces, which include both the previously can be assigned to layer 3. Other spaces and open decks,
mentioned machinery space category A and other spaces where fire hoses are required, can be assigned to layer 4. By
containing propulsion machinery, boilers, fuel units, steam using different layers, the design team can easily distinguish
and internal combustion engines, generators, and main between different types of rooms and apply different design
criteria or safety measures accordingly.
Fig. 3. Possible parameterization options for the length of a fire hose according to the SOLAS convention
Fig. 10. Comparison of the actual arrangement of fire hydrants on a fishing vessel with their simulated distribution
based on two selected criteria: maximizing the degree of coverage of extinguished areas and minimizing the distance from evacuation exits
with weights w1 = 0,6 and w4 = 0,4 respectively
Fig. 11 Comparison of the actual arrangement of fire hydrants on a vessel servicing wind platforms
with their simulated distribution based on two selected criteria: maximizing the degree of coverage of extinguished areas and
minimizing the distance from evacuation exits with weights w1 = 0,6 and w4 = 0,4 respectively
Based on the comparative analysis of the individual and details, the larger the file size, and consequently, the
solutions regarding the arrangement of fire hydrants on the computer’s computational power required for processing and
decks of various ships, the following can be concluded: displaying graphics in real-time becomes greater.
–– both the layout of the fire hydrants obtained using the Tab. 1. Simulation duration for the Wind Platform Vessel, used for wind
classical method and the one obtained using the developed platform servicing, based on the criteria of maximizing the coverage
system meet all the imposed requirements (constraints), degree of extinguished surfaces and minimizing the distance to
evacuation exits
–– the layout of the fire hydrants obtained using the developed
system enables faster access by the firefighting team to number of hydrants
the hydrants, thanks to the criterion of minimizing the
9 12 15 18 25 50
distance from evacuation exits as a component of the
substitute objective function, which significantly increases Max. speed = 0.1 11 11 11
12[s]
14 15
the level of protection against fire hazards on the ship. Iterations: 20 [s] [s] [s] [s] [s]
Considering that the main goal of the developed system Max. speed = 0.2 11 11 11 14 15
12[s]
is to reduce the workload of the designer, the above- Iterations: 20 [s] [s] [s] [s] [s]
mentioned benefits are an additional but very significant
Max. speed =` 0.3 11 11 12 14 15
element of the feasibility of the developed system. So, the Iterations: 20 [s] [s] [s]
12[s]
[s] [s]
basic question remains – how quickly can a solution related
to the arrangement of fire hydrants be obtained using the Max. speed = 0.5 11 11 11 13 15
12[s]
Iterations: 20 [s] [s] [s] [s] [s]
developed system?
In order to answer the question above, a series of simulations Max. speed = 1 11 11 11
12[s]
14 15
were conducted using a computer with typical parameters Iterations: 20 [s] [s] [s] [s] [s]
for computers used in ship design offices. The results of Max. speed = 0.1 21 22 22 22 24 29
the simulation duration for the Wind Platform Vessel are Iterations: 40 [s] [s] [s] [s] [s] [s]
shown in Table 1. Similar simulations were conducted for the
Max. speed = 0.2 21 21 22 23 24 28
remaining three ships. Using the PSO algorithm, the number Iterations: 40 [s] [s] [s] [s] [s] [s]
of iterations, speed, and number of particles (individuals)
in the search space were varied. The type of analyzed vessel Max. speed = 0.3 21 22 22 23 24 28
Iterations: 40 [s] [s] [s] [s] [s] [s]
seemingly had an impact on the algorithm’s speed, but this
was solely due to the size of the graphic file that served as Max. speed = 0.5 21 21 22 22 24 29
the source of the analyzed unit’s general plan. The size of the Iterations: 40 [s] [s] [s] [s] [s] [s]
graphic file depends on various factors, such as the number Max. speed = 1 21 22 22 22 24 29
of graphic elements, resolution, number of layers, styles, Iterations: 40 [s] [s] [s] [s] [s] [s]
and blocks contained in the file. The more graphic elements
Yordan Denev
Technical University of Varna, Bulgaria
Abstract
This report examines the feasibility and impact of retrofitting the bulbous bow on a general cargo ship, in terms of the
energy efficiency operational index (EEOI), in the areas of Western Europe and the Eastern Mediterranean. Three
ship forms were developed and analysed: with a bulbous bow, without a bulbous bow, and with a modified bulbous
bow. The goal in developing the ship forms and conducting the analysis was to achieve minimal differences in the ship’s
characteristics with the same volumetric displacement, aided by PolyCAD software. A route was selected between
two ports: Varna and Rotterdam. The labour intensity of the bulbous bow retrofitting process was evaluated and
approximate values of labours costs and cost for the task were determined. The results obtained for resistance during
ship motion, EEOI, and fuel consumption reductions, or increases, were compared against the retrofitting values. The
return cost of retrofitting is evaluated and measured in terms of fuel saved.
Keywords: retrofitting, bulbous bow, energy efficiency design index, return costs
Hull modelling
The hull form was generated by PolyCAD software. To HULL FORM GENERATION
assess the effect of the modification of the ship’s bow, the
resistance was calculated using the Holtrop and Mennen The main dimensions of the ship were L=120.62 m,
method, for speeds ranging from 0-17 knots. The advantage B= 16.00 m, D= 9.03 m, and d= 6.67 m. The ship had a
of the software is the possibility of recalculating ship double bottom and double side, single deck and one hold,
characteristics in the event of some form of change. After with a 116 TEU container capacity. The hold length was
retrofitting, a small difference in mass displacement appeared. about 84.5m, with a double sided width of 1.3 m per side
and a maximum hold breadth of 13.4 m. The service speed
was 15 kn and the main engine type was a ‘5S35ME’ with
MODEL VALIDATION AND VERIFICATION main engine power of 4350 kW. Three different forms with
similar hull coefficients are shown in Table 1. The original
In this type of analysis, it is important for the calculations hull form was without a bulbous bow (VAR1), while the other
to be within a range of 5% tolerance, which is assumed for two had bulbous bows, where the dimensions of the bulb
engineering calculations. Otherwise, if there is more than were different.
5% tolerance, the impact on the characteristics is significant. Tab. 1. Hull form coefficients
The retrofitting process consisted of modernising the forward
ship hull’s form without making changes to the ship’s main VAR1 VAR2 VAR3
dimensions. Cb 0.78 0.76 0.76
The change in the geometry of the bow was achieved by Cp 0.78 0.77 0.76
mounting a bulb with a specific geometry that corresponds Cm 0.99 0.99 0.99
to the original ship’s form. In this case, in order to ensure a Cw 0.90 0.89 0.87
constant displacement of water, changes were made to the
Δ, m3 10518.34 10518.64 10518.80
coefficients of the shape, specifically the prismatic coefficient
(Cp) and, consequently, the block coefficient (Cb). The
coefficient of the mid-ship section remained the same for all The analysed ship was a general cargo ship with one hold
shapes. Differences in Cb were within 2.5%, and differences of 7000 tDW. The location of the collision bulkhead was at
in Cp were in the range 0.7-2.5%; higher differences in Cw 7% of Lpp, while the engine room bulkhead was at 23% Lpp,
coefficients were in the range 1.8-3.5%. The maximum see Fig. 2 to Fig. 4.
difference in mass displacement was 0.46 t.
In the evaluation of the model shape, a mesh with
rectangular and triangular elements was used. Each type
of element was used in different areas of the ship’s hull. For
example, in the bow and stern regions where the hull shape
has complex curvature in two directions, triangular mesh
elements were used, while rectangular elements were used
in the remaining areas. The transition elements between
rectangular and triangular elements were rhomboidal.
The grid spacing was 0.2 m with a key nudge of 0.002 m; the
number of elements is shown in Fig. 1.
Fig. 2. Variant one (original hull form) without bulbous bow
Despite the fact that the specific fuel consumption for CAPITAL
the EXPENDITURE
Ship steel hull weight FOR RETROFITTING
estimation ANDbyRETURN
can be evaluated the COST
second variant of the ship’s hull form is the lowest, it does mathematical equations presented in [14] but they are not
Fig.7. Voyage fuel cost and specific fuel consumption
not result in the lowest overall fuel cost when evaluating The effects of appropriate
the in this case,
forward part retrofitting willbecause hull
be clearer weight
after is calculated
calculating capital expenditure, retu
total fuel expenses. This is due to the fact that the difference in relation to weight displacement. In the case
costs, and time for return costs. To study this effect, it is necessary to calculatestudy, the
the hull steel weig
in specific
Considering the mainfuel consumption
engine-specific between
consumption, shipthe second
service and
speed, thirdbetween
distance
billable hours volumetric
for its fabrication,and
ports, weight
and their displacement
differences are forms.
for different the same
Afterfor all
production calculation
marine diesel prices, the fuel cost for one voyage was calculated using the CEAS engine calculator
variants is only 0.8 g/kWh in favor of the second one,
the which
return cost and forms,
time are and computer model development is used for hull
calculated.
cific for different engine powers) and is presented in Fig.6. Average marine diesel oil costs 0.42
does not make
er, which corresponds to 546.5a $/mt
significant
of very low impact
sulfuron fuelthe
oilend result,toasprices
(according the fromweight evaluation. The results are shown in Table 2.
July 19,estimation
Ship steel hull weight can be evaluated by the mathematical equations presented in [1
resistance that needs to be overcome with the second hull Tab. 2. Ship hull
) [12]. but they are not appropriate inweights
this case, because hull weight is calculated in relation to weig
form is higher than the third one.
Despite the fact that the specific fuel consumption for the second displacement. In the
variant of the ship's case
hull formstudy, the volumetric and weight displacement
VAR1 VAR2 areVAR3
the same for all forms, a
e lowest, it does not result in the lowest overall fuel cost when evaluating the total
computer fueldevelopment
model expenses. is used for hull weight evaluation. The results are shown in Table 2.
WEIGHT HULL, T 835.1 852.0 894.0
is due to the fact that the difference in specific fuel consumption between the second and third
ENERGY
nts is only 0.8 g/kWh EFFICIENCY
in favor of OPERATIONAL
the second one, which does not make a significant impactBulb area, m 2
on the
Table 2. Ship hull weights
0.00 10.30 9.03
INDEX EVALUATION BASED ON VOYAGE
esult, as the resistance that needs to be overcome with the second hull form is higher than the third
Bulb length, m 0.00 VAR1 2.00
VAR2 VAR3 3.00
PARAMETERS 852.0
1.50 894.02.00 WEIGHT
Bulb radius, m HULL, T0.00 835.1
2
Bulb area, m 0.00 10.30 9.03
EEOI is an indicator for evaluating ship energy efficiency Bulb length, m 0.00 2.00 3.00
ENERGY EFFICIENCY OPERATIONAL INDEX EVALUATION BASED Bulb radius, m
and CO emissions to the environment during a ship’s
2
The difference in hull weights,
0.00 after1.50
retrofitting,
2.00 is about
ON VOYAGE PARAMETERS
operation and through her life cycle. Using the equation 17 t of steel construction for variant VAR2 and 59 t of steel
for EEOI offered in [10], for generated ship hulls,The indexes
difference construction
in hull weights,forafter
variant VAR3. is
retrofitting, Differences of such
about 17 t of magnitude for variant VAR
steel construction
are calculated considering voyage parameters and
and 59 actual
t of steel do not
construction
EEOI is an indicator for evaluating ship energy efficiency and CO2 emissions to the environment affect
for the ship’s
variant carrying
VAR3. capacity
Differences of since,
such during
magnitude thenot affect the shi
do
ship conditions. For a voyage from the port of Varna
carrying to the
capacity
ng a ship’s operation and through her life cycle. Using the equation for EEOI offered in [10], for conceptual
since, during design
the stages,
conceptual a 1%
design reserve
stages, a 1%displacement
reserve is
displacement is provide
port of
rated ship hulls, Rotterdam
indexes with considering
are calculated a speed ofvoyage
15 kn,parameters
a distance
theandofactual
about
results are shown
ship provided;
in Fig.9.the results are shown in Fig. 9.
conditions.
3940
a voyage from thenm, a Varna
port of deadweight
to the portofofabout 7000
Rotterdam t, and
with a fuel
a speed of 15 carbon
kn, a distance of about
content of
nm, a deadweight ofabout
0.86 7000
for light
t, andfuel oilcarbon
a fuel [13], the EEOI
content calculated
of 0.86 fuel oil [13],198000
for light by the
I calculated by Eq. (1) is presented in Fig.8. 196000
Eq. (1) is presented in Fig. 8.
Hull fabrication, mh/t
194000
5,22E-06
192000
5,20E-06 190000
5,18E-06 188000
EEOI, tCO2/tnm
5,16E-06 186000
5,14E-06
184000
5,12E-06
182000
5,10E-06
0 1 2 3 4
5,08E-06
Variant number
5,06E-06
0 1 2 3 4
Variant number Fig. Fig.9.
9. Billable
Billablehours forsteel
hours for steelhull
hull fabrication
fabrication
Jan Sierzputowski *
Artur Karczewski
Przemysław Krata
Gdansk University of Technology, Institute of Naval Architecture, Gdansk, Poland
Abstract
A sailing yacht is a human-centred product, the design of which revolves primarily around the wants and desires of
the future owner. In most cases, these preferences are not measurable, such as a personal aesthetic feeling, or a need for
comfort, speed, safety etc. The aims of this paper are to demonstrate that these preferences can be classified and represented
numerically, and to show that they are correlated with the type of yacht owned. As a case study, the owner’s preferences
for deck equipment are considered. These are determined by pairwise comparisons of the importance rankings for features
previously defined by yacht owners, following the analytic hierarchy process (AHP) method. As a result, a quantitative
representation of these preferences is established, and they are shown to be correlated with the type of yacht. The findings
of the current study show that the yacht owners’ preferences can be represented numerically, leading to a utilitarian
conclusion that concerns the support and even some degree of automation of the design process.
Keywords: Yacht Design, AHP Method, Yacht User Preferences, Sailing Yacht, Human-Centred Design
Fig. 2 Hierarchy structure tree of sailing yacht deck equipment applying the AHP method
Fig. 3 Yacht owners’ preferences, determined based on aggregated consistent individual judgements by yacht type
The individual judgements were corrected, and new priorities impact on the final CR, which was greater than 0.10 only for the
and CR values were calculated. The new results are shown in group of racing yachts. The results obtained from this process
the table in Appendix 2. were very similar to those of the improved answers, as shown
in Fig. 4 and Table 5.
Group aggregation results
After aggregating the results according to the values specific
to each yacht (such as the overall length (LOA), breadth (B), DISCUSSION
displacement (D), etc.), according to the factors of slenderness,
comfort, etc., it was concluded that the purpose of the yacht The use of a nine-point rating scale and a pairwise comparison
had the greatest correlation with the priorities. The individual of features seems to be a good method for determining priorities
judgements were divided into groups related to the type of among yacht owners. The results showed that durability was the
yacht owned, and then aggregated using the geometric mean most important feature for expedition and seagoing cruising
method [32]. As a result, the percentage degree of importance yachts, whereas the efficiency and mass of an accessory were the
was obtained for each of the features of the deck equipment, most important aspects for racing and racing/cruising yachts.
depending on the type of yacht. The results are shown in Fig. 3 The owners of inland/coastal cruising yachts had the most
and Table 4. balanced priorities (regardless of the mass of an accessory) and
In order to explore the impact of improving the respondents’ valued aesthetics most highly.
answers on the relevance of the results, inconsistent judgments The use of the AHP method imposes certain limitations
were also aggregated for comparison. It was found that the on the results. The method used to conduct the study, the
aggregation method using the geometric mean had a positive number of respondents, their questionable proficiency in the
Fig. 4 Yacht owners’ preferences determined based on aggregated inconsistent individual judgements by yacht type before CR improvement
15. S. Ö. Felek, “Parametric sailing yacht exterior and interior 28. A. Jozaghi et al., “A comparative study of the AHP and
design,” 2020, doi: 10.14744/tasarimkuram.2019.30085. TOPSIS techniques for dam site selection using GIS: A case
study of Sistan and Baluchestan Province, Iran,” Geosciences
16. D. Harris, S. Mccartan, B. Verheijden, H. Groningen, and (Switzerland), vol. 8, no. 12, Dec. 2018, doi: 10.3390/
M. Lundh, “European Boat Design Innovation Group: geosciences8120494.
The marine design manifesto. Human factors in safety
management systems view project next generation civil flight 29. F. de Felice and A. Petrillo, “Absolute measurement with
deck concepts view project,” 2014. [Online]. Available: https:// analytic hierarchy process: A case study for Italian racecourse,”
www.researchgate.net/publication/280065046. International Journal of Applied Decision Sciences, vol. 6, no. 3,
p. 209, 2013, doi: 10.1504/IJADS.2013.054931.
17. M. Bilski, “Selected Human Factors in Marina Design,”
Procedia Manuf., vol. 3, pp. 1646–1653, 2015, doi: 10.1016/j. 30. F. de Felice and A. Petrillo, “Multicriteria approach for process
promfg.2015.07.482. modelling in strategic environmental management planning,”
International Journal of Simulation and Process Modelling,
18. T. Bosma, “A human factor’s approach to mega yacht concept vol. 8, no. 1, p. 6, 2013, doi: 10.1504/IJSPM.2013.055190.
design,” Master of Philosophy, University of Strathclyde,
Glasgow, 2013. 31. J. Benítez, X. Delgado-Galván, J. Izquierdo, and R. Pérez-García,
“Improving consistency in AHP decision-making processes,”
19. R. Gill, “Approaches to design,” Des. Stud., vol. 1, no. 3, Appl. Math. Comput., vol. 219, no. 5, pp. 2432–2441, Nov. 2012,
pp. 141–145, Jan. 1980, doi: 10.1016/0142-694X(80)90020-4. doi: 10.1016/j.amc.2012.08.079.
20. P. von Buelow, ‘Using Evolutionary algorithms to aid 32. J. Krejčí and J. Stoklasa, “Aggregation in the analytic hierarchy
designer of architectural structures,” in Creative Evolutionary process: Why weighted geometric mean should be used instead
Systems, Elsevier, 2002, pp. 315–336. doi: 10.1016/ of weighted arithmetic mean,” Expert Syst. Appl., vol. 114,
B978-155860673-9/50050-1. pp. 97–106, Dec. 2018, doi: 10.1016/j.eswa.2018.06.060.
Abstract
Ballast water management is an effective measure to ensure that organisms, bacteria and viruses do not migrate
with the ballast water to other areas. In 2004, the International Maritime Organization adopted the International
Convention on the Control and Management of Ballast Water and Ship Sediments, which regulates issues related to
ballast water management. Many technologies have been researched and developed, and of these, the use of UV rays
in combination with filter membranes has been shown to have many advantages and to meet the requirements of
the Convention. However, the use of UV furnaces in ballast water treatment systems requires a very large capacity,
involving the use of many high-power UV lamps. This not only consumes large amounts of electrical energy, but is
also expensive. It is therefore necessary to find an optimal algorithm to enable the UV radiation for the UV controller
in the ballast water sterilisation process to be controlled in a reasonable and effective manner. This controller helps
to prolong the life of the UV lamp, reduce power consumption and ensure effective sterilisation. This paper presents
a UV control algorithm and a controller for a UV furnace for a ballast water treatment system installed on a ship.
The results of tests on vessels illustrate the effect of the proposed UV controller.
Keywords: UV Quantity Controller, Ship, Ballast Water Management, Viable Organisms Threat, Marine Environment.
Nomenclature
UV dose
j
(k) = UV dose
j
(k–1) · (1 – FV· τ ) +
j
UV dose
a
· FV· τ + i · τ (3)
j
Fig. 7. Variables used by the fuzzy controller structure for the UV furnace
The method used to determine the control output value Uʹ Uʹʹ(s) = 1s · (K3 · Uʹ(s) + X(s)) (10)
is illustrated in Figure 9. We apply Eq. (7) to calculate the value
of U1 for segment A1A2 and U2 for segment B1B2 depending on where X(s) = Ks (U(s)−Uʹʹ(s)) is the value of the anti-saturation
the value of eʹ, as follows: integral. We convert Eq. (10) to the Z domain using the Tustin
method:
2
U1 = A 1 – A2
eʹ1 – eʹ2 · (eʹ1 – eʹ2) + A2 s = T · zz –+11 (11)
To implement Eq. (12) in the PLC, we transform it into The basic components of a fuzzy controller are the fuzzy
a differential equation as follows. stage, the composition rule, and the defuzzification stage.
Since basic fuzzy controllers are only capable of processing
T
Uʹʹ(k) – Uʹʹ(k–1) = current signals, these are called static fuzzy controllers. To
2 · (K3 · Uʹ(k) + K3 · Uʹ(k–1) + X(k) + X(k–1)) extend their application to dynamic control problems, the
necessary kinematics are added to the basic fuzzy controller
T
Uʹʹ(k) = Uʹʹ(k–1) + 2 · to provide it with the derivative or integral value of the signal.
(K3 · Uʹ(k) + K3 · Uʹ(k–1) + X(k) + X(k–1)) When used with these kinematics, the basic fuzzy controller
is called a dynamic fuzzy controller.
T
Uʹʹ(k) = Uʹʹ(k–1) + 2 · For the ballast water treatment UV reactor, this is
(K3 · Uʹ(k) + K3 · Uʹ(k–1) + X(k–1)) + T2 · X(k) a nonlinear object. This nonlinearity is expressed in the
(13) relationship between the flow rate F and the amount of UV,
If we set F = Uʹʹ(k–1)+ T2 · (K3 · Uʹ(k) + K3 · Uʹ(k–1) + X(k–1)), then and between the UV lamp control signal (Udk) and the UV
intensity. The relationship between UVdose and the flow F is
T
Uʹʹ(k) = F + 2 · X(k) (14) described in general terms by the system of equations in (13)
and (14). These systems are difficult to control with PIDs, and
To calculate Uʹʹ(k), we consider three cases: experiments are required to determine the accuracy of the
Case 1: Uʹʹ(k) [0,1] controller. In this case, other controllers are often used (for
Then X(k) = ks · (U(k) − Uʹʹ(k)) = 0. From Eq. (14), we have: example based on fuzzy control, neural control, or adaptive
control) for the system. In the present paper, we build a fuzzy
Uʹʹ(k) = F (15) controller, as this has certain outstanding features: (i) the
Case 2: Uʹʹ(k) >1 control is based on the operator’s experience; (ii) there no
Then X(k) = ks · (1 − Uʹʹ(k)), and from Eq. (14) we get: need for an object model to set up the controller; and (iii) it
can be applied to industrial control devices such as PLCs or
T
Uʹʹ(k) = 2 + k2s · T (F + 2 · ks) (16) microprocessors [43].
Uʹʹ(k) = 2 +2F
ks · T (17)
1 khi Uʹʹ(k) > 1 The fuzzy control structure for the UV reactor is shown in
U(k) = Uʹʹ(k) khi 0 ≤ Uʹʹ(k) ≤ 1 Figure 12. In this structure, the controller has two inputs and
0 khi Uʹʹ(k) < 0 one output signal. The output signal of the controller is in the
(18) range [0,1], corresponding to a UV intensity of the lamp of
The algorithm for calculating the value of the fuzzy controller between zero and Imax.
in the PLC for the UV furnace is summarised in Figure 11. The two inputs of the fuzzy controller are: wrong order
between the UV amount set and the UV amount of water
treated in the reactor. Rate of variation of UV content of ballast
water.
To facilitate the construction of the model and the signal
processing step, all of the input values to the fuzzy controller
are converted to standard values between −1 and 1. To do
this, we need to pass additional inputs in the form of the
coefficients K1 and K 2 , as shown in Figure 13.
Fig. 11. Algorithm for calculating the output value of the fuzzy control
IMO standard
Treatment rated capacity
Name of ship Viable organisms ≥50μm Viable organisms 10-50μm Sampling period
(m3/h)
(org/m3): <10 (org/ml): <10
ANBIEN BAY 500 4 6.8 10h50–14h00, 02/06/2022
TTC PIONEER 200 No live organisms detected 7.8 11h10–12h15, 06/06/2022
ROYAL 89 100 6 4.5 09h50–10h40, 15/08/2022
TRACY 150 1 2.9 15h00–15h30, 20/09/2022
AN THINH PHU 08 100 1 0.2 08h00–09h00, 14/10/2022
THAI BINH 35 50 1 1.7 11h30–12h10, 02/11/2022
1. G. Altug, S. Gurun, M. Cardak, P. S. Ciftci, and S. Kalkan, 11. IMO, International Maritime Organization [IMO], 2004.
“The occurrence of pathogenic bacteria in some ships’ ballast International Convention for the Control and Management
water incoming from various marine regions to the Sea of of Ships’ Ballast Water and Sediments. International
Marmara, Turkey,” Mar. Environ. Res., vol. 81, pp. 35–42, Maritime Organization, London. [Online]. Available: http://
Oct. 2012, doi: 10.1016/j.marenvres.2012.08.005. www.imo.org/en/About/Conventions/List ofConventions/
Pages/International- Convention-for-the-Control-and-
2. L. J. Aridgides, M. A. Doblin, T. Berke, F. C. Dobbs, Management-of-Ships%27-Ballast-Water-and-Sediments-
D. O. Matson, and L. A. Drake, “Multiplex PCR allows (BWM).aspx
simultaneous detection of pathogens in ships’ ballast water,”
Mar. Pollut Bull., vol. 48, no. 11, pp. 1096–1101, Jun. 2004, 12. G. Elidolu, S. I. Sezer, E. Akyuz, O. Arslan, and Y. Arslanoglu,
doi: 10.1016/j.marpolbul.2003.12.017. “Operational risk assessment of ballasting and de-ballasting
on-board tanker ship under FMECA extended evidential
3. J. Chen, Y. Lin, J. Zhou Huo, M. Xia Zhang, and Z. Shang Ji, reasoning (ER) and rule-based Bayesian network (RBN)
“Optimization of ship’s subdivision arrangement for offshore approach,” Reliab. Eng. Syst. Saf., vol. 231, p. 108975, Mar. 2023,
sequential ballast water exchange using a non-dominated doi: 10.1016/j.ress.2022.108975.
sorting genetic algorithm,” Ocean Engineering, vol. 37, no. 11,
pp. 978–988, Aug. 2010, doi: 10.1016/j.oceaneng.2010.03.012. 13. D. K. Gray, I. C. Duggan, and H. J. Macisaac, “Can sodium
hypochlorite reduce the risk of species introductions from
4. I. Staffell and P. Balcombe, “How to decarbonise international diapausing invertebrate eggs in non-ballasted ships?,” Mar.
shipping: options for fuels, technologies and policies,” Energy Pollut. Bull., vol. 52, no. 6, pp. 689–695, Jun. 2006, doi:
Convers. Manag., [Online]. Available: https://www.academia. 10.1016/j.marpolbul.2005.11.001.
edu/38669854/How_to_decarbonise_international_
shipping_options_for_fuels_technologies_and_policies 14. P. P. Sangave, A. Mukherjee, and A. Pandit, “Ballast water
treatment using hydrodynamic cavitation,” Nov. 2014.
5. P. Van Hung, K.-S. Kim, L. Q. Tien, and N. M. Cuong, [Online]. Available: https://www.semanticscholar.org/paper/
“Distribution of oil spill response capability through Ballast-Water-Treatment-Using-Hydrodynamic-Sangave-
considering probable incident, environmental sensitivity Mukherjee/ ca3f5edcc2a28ba175b1a4c0093b70b49bb2a334
and geographical weather in Vietnamese waters,” Journal
of International Maritime Safety, Environmental Affairs, 15. S. Rahman, “Implementation of ballast water management
and Shipping, vol. 2, no. 1, pp. 31–41, Nov. 2018, doi: plan in ships through ballast water exchange system,”
10.1080/25725084.2018.1511240. Procedia Eng., vol. 194, pp. 323–329, Jan. 2017, doi: 10.1016/j.
proeng.2017.08.152.
6. G. Drillet et al., “Improvement in compliance of ships’
ballast water discharges during commissioning tests,” Mar. 16. E. Tsolaki and E. Diamadopoulos, “Technologies for ballast
Pollut. Bull., vol. 191, p. 114911, Jun. 2023, doi: 10.1016/j. water treatment: A review,” Journal of Chemical Technology
marpolbul.2023.114911. & Biotechnology, vol. 85, no. 1, pp. 19–32, Jan. 2010, doi:
10.1002/jctb.2276.
18. C. Grob and B. G. Pollet, “Regrowth in ship’s ballast water 30. IMO, “Globallast guidelines for development of a national
tanks: Think again!” Mar. Pollut. Bull., vol. 109, no. 1, ballast water management Strategy, 1997. Guidelines for
pp. 46–48, Aug. 2016, doi: 10.1016/j.marpolbul.2016.04.061. preventing the introduction of unwanted organisms and
pathogens from ships ballast waters and sediment discharges.
19. V. Raţᾰ and L. Rusu, “Ballast water pollution risk assessment in Resolution A.868 (20).” 2004.
the Black Sea,” in Mechanical Testing and Diagnosis, Jan. 2021,
pp. 35–40. doi: 10.35219/mtd.2020.4.05. 31. G. H. Briton, B. Yao, and G. Ado, “Evaluation of the Abidjan
lagoon pollution,” Journal of Applied Sciences and Environmental
20. E. Briski et al., “Combining ballast water exchange and Management, vol. 11, no. 2, 2007, doi: 10.4314/jasem.
treatment to maximize prevention of species introductions to v11i2.55030.
freshwater ecosystems,” Environ. Sci. Technol., vol. 49, no. 16,
pp. 9566–9573, Aug. 2015, doi: 10.1021/acs.est.5b01795. 32. J. Liu, P. Wang, G. Liu, J. Dai, J. Xiao, and H. Liu, “Study of the
characteristics of ballast bed resistance for different temperature
21. A. Travizi et al., “Macrozoobenthos in the Adriatic Sea and humidity conditions,” Constr. Build. Mater., vol. 266,
ports: Soft-bottom communities with an overview of non- p. 121115, Jan. 2021, doi: 10.1016/j.conbuildmat.2020.121115.
indigenous species,” Mar. Pollut. Bull., vol. 147, pp. 159–170,
Oct. 2019, doi: 10.1016/j.marpolbul.2019.01.016. 33. N. A. Salleh et al., “Pathogenic hitchhiker diversity on
international ships’ ballast water at West Malaysia port,” Mar.
22. P. Mozetič et al., “Phytoplankton diversity in Adriatic ports: Pollut. Bull., vol. 172, p. 112850, Nov. 2021, doi: 10.1016/j.
Lessons from the port baseline survey for the management of marpolbul.2021.112850.
harmful algal species,” Mar. Pollut. Bull., vol. 147, pp. 117–132,
Oct. 2019, doi: 10.1016/j.marpolbul.2017.12.029. 34. G. Romanelli et al., “Ballast water management system:
Assessment of chemical quality status of several ports in
23. J.-H. Park, Y.-B. Sim, S.-Y. Kang, and S.-H. Kim, “Inactivation Adriatic Sea,” Mar. Pollut. Bull., vol. 147, pp. 86–97, Oct. 2019,
of indicating microorganisms in ballast water using chlorine doi: 10.1016/j.marpolbul.2017.12.030.
dioxide,” Ecology and Resilient Infrastructure, vol. 5, no. 3,
pp. 111–117, Sep. 2018, doi: 10.17820/ERI.2018.5.3.111. 35. M. R. First and L. A. Drake, “Life after treatment: detecting
living microorganisms following exposure to UV light and
24. T. McCollin, G. Quilez-Badia, K. D. Josefsen, M. E. Gill, chlorine dioxide,” J. Appl. Phycol., vol. 26, no. 1, pp. 227–235,
E. Mesbahi, and C. L. J. Frid, “Ship board testing of Feb. 2014, doi: 10.1007/s10811-013-0049-9.
a deoxygenation ballast water treatment,” Mar. Pollut. Bull.,
vol. 54, no. 8, pp. 1170–1178, Aug. 2007, doi: 10.1016/j. 36. Z. Manxia, L. Shengjie, T. Xiaojia, L. Xiang, and Z. Yimin,
marpolbul.2007.04.013. “Evaluation of micro-pore ceramic filtration and uv radiation
combination on ballast water treatment,” in 2010 International
25. L. A. Drake, M. A. Doblin, and F. C. Dobbs, “Potential Conference on Digital Manufacturing & Automation,
microbial bioinvasions via ships’ ballast water, sediment, pp. 670–673, Dec. 2010, doi: 10.1109/ICDMA.2010.374.
and biofilm,” Mar. Pollut. Bull., vol. 55, no. 7, pp. 333–341,
Jan. 2007, doi: 10.1016/j.marpolbul.2006.11.007. 37. P. P. Stehouwer, A. Buma, and L. Peperzak, “A comparison
of six different ballast water treatment systems based on
26. S. Vodyanitskaya et al., “Methods of decontamination of ship UV radiation, electrochlorination and chlorine dioxide,”
ballast water with polyguanidine disinfectant,” International Environ. Technol., vol. 36, no. 16, pp. 2094–2104, Aug. 2015,
Journal of Infectious Diseases, vol. 79, p. 77, Feb. 2019, doi: doi: 10.1080/09593330.2015.1021858.
10.1016/j.ijid.2018.11.195.
38. I. Rivas-Zaballos, L. Romero-Martínez, I. Moreno-Garrido,
27. E. Lakshmi, M. Priya, and V. S. Achari, “An overview on the J. Moreno-Andrés, A. Acevedo-Merino, and E. Nebot,
treatment of ballast water in ships,” Ocean Coast Manag., vol. 199, “UV-LEDs combined with persulfate salts as a method to
p. 105296, Jan. 2021, doi: 10.1016/j.ocecoaman.2020.105296. inactivate microalgae in ballast water,” Journal of Water Process
Engineering, vol. 51, p. 103361, Feb. 2023, doi: 10.1016/j.
28. M. N. Tamburri, K. Wasson, and M. Matsuda, “Ballast water jwpe.2022.103361.
deoxygenation can prevent aquatic introductions while
reducing ship corrosion,” Biol. Conserv., vol. 103, no. 3, 39. J. Xiao, Y. Xu, L. Hu, and H. Wu, “Evaluating the treatment
pp. 331–341, Mar. 2002, doi: 10.1016/S0006-3207(01)00144-6. performance of filtration & real-time UV irradiation
Se-yeol An1
Hyeon-won Jeong1
Ohyoung Kim2
W. Jaewoo Shim1,*
1
Department of Chemical Engineering, Dankook University, Gyeonggi-do, 16890, South Korea
2
School of Polymer Science and Engineering, Dankook University, Gyeonggi-do, 16890, South Korea
Abstract
This study investigates the effects of sway and roll excitations on sloshing liquid loads in a tank, using Ansys Fluent
software. The model considered in the study is a 1:50 scaled membrane-type tank, based on a KC-1 membrane LNG
tank designed by Korea Gas Corporation (KOGAS). The volume of fluid (VOF) method is used to track the free surface
inside the tank, and the standard k-ε model is applied to express the turbulent flow of the liquid. To explore the motion
of the tank under excitation, a user-defined function (UDF) and a dynamic mesh technique are employed to control
the external forces exerted on the tank through its motion. The results, in the form of time series data on the sloshing
pressures in the tank under pure sway, roll, and coupled sway-roll, are analysed, with specific ranges for the excitation
amplitudes and frequencies. We show that variations in excitation frequency and amplitude significantly influence
the sloshing loads. Sloshing loads are found to intensify when the excitation frequency matches the tank’s primary
natural frequency, 1.0 ω'1. Furthermore, with coupled sway-roll excitations, the sloshing loads are weakened when the
sway and roll are in-phase and are intensified when these are out-of-phase. Fast Fourier transform analysis provides
insights into the frequency domain, showing that the dominant frequency is 0.88 Hz and it is approximately equal to
the tank’s primary natural frequency, 1.0 ω'1.
Keywords: KC-1 membrane LNG tank; Sloshing; Sway; Roll; Coupled sway-roll; Fast Fourier transform (FFT)
where 𝑢𝑢
⃗ denotes the velocity vector, 𝜌𝜌 the density of the mixture, 𝑝𝑝 the static pre
excitation induces a linear transverse (lateral
𝑋𝑋 andor side-to-side) themotion along the tank’s tran
expression
displacement, for the roll
respectively. 𝑋𝑋0 represents
excitation. Here, the 𝜃𝜃amplitude
represent of lateral
the displacement
lateral and
displacement, 𝜃𝜃
th
whereas roll excitation induces
𝑋𝑋 a tilting rotation of the of tank about itsdisplacement,
longitudinal ax𝜃𝜃
example
An of the
example distribution
of the of 𝛼𝛼
distribution Fig.𝛼𝛼 1.
ofvalues Aninexample
values cells [8]of
in cells
displacement,
[8]the distribution ofthe amplitude
𝛼𝛼 values in cells [8] of the rotational0 displacement, and 𝜔𝜔 is the frequency of the excitations. 0
respectively. represents the amplitude the lateral
centre
the of the bottom
amplitude of the tank
of the rotational acting as theand
displacement, centre
𝜔𝜔 isofthe
thefrequency
coordinate of system. If the sw
the excitations.
NING
the distribution of 𝛼𝛼 values
EQUATIONS in An
Fig. 1. cells excitations
[8] of the distribution of 𝛼𝛼 values in cells [8]
example
are periodic, they can be represented in the form of sinusoidal waves, as sho
Fig. 1. An example of the distribution of 𝛼𝛼 values in cells [8]
G
ofvaluesEQUATIONS
𝛼𝛼study, in cellsin a
Fig. 1. An example of the distribution of 𝛼𝛼 values in cells [8]
[8]conventional numerical
𝑋𝑋 = 𝑋𝑋 𝑠𝑠𝑠𝑠𝑠𝑠 (𝜔𝜔𝜔𝜔) 0
numerical approach is isemployed toapproach
tocalculate istheemployed to calculate the sloshing
issloshing
values
numerical cells [8]
approach employed calculate the sloshing
EQUATIONS
non of liquid fraction in a for the
tank. It two
is assumed phases,
Fig. 1. An thatliquid
example the of
and
fluid
the
gas,
flow
distribution
calculated
isof 𝛼𝛼
incompressible,
values in
using
cells [8] and the turbulent
assumed
is assumed that that thethe fluid fluid flow flow is is incompressible,
incompressible, and and thethe turbulent
turbulent 𝜃𝜃 = 𝜃𝜃0 𝑠𝑠𝑠𝑠𝑠𝑠 (𝜔𝜔𝜔𝜔)
he
NING
sloshing thephenomenon
EQUATIONS volume fraction is represented equationby inthe Eq.Reynolds-averaged
(2): Navier-Stokes (RANS) (6)
al
is
dy, approach
represented
represented
a conventional byisby the employed
the Reynolds-averaged
Reynolds-averaged
numerical to calculate approach theis sloshing
Navier-Stokes
Navier-Stokes employed (RANS)
(RANS)to calculate the sloshing
nd
of The
d,GOVERNING
that
theathe
liquid Reynolds
the
turbulent influid
turbulent
conventional a tank. stresses
flow
flow EQUATIONS
flow is field
field
Itnumerical and
incompressible,
assumed in in thethe
the turbulent
RANS
that
approach RANStheand are
fluid flow
are
is the
calculated field
turbulent
calculated
flow
employed inbased the
based
is incompressible,
to is RANS onon areand
calculate calculated
the thesloshing
turbulent based on
sard
me nted study,
fraction k-ε by
fraction model.
the a for
for conventional
The
Reynolds-averaged
the the two volume
two phases, phases, numerical
fraction
liquid liquid for
Navier-Stokes
and approach
and the
gas, gas,two
is
Eq.
is
phases,
(RANS)
calculated
is calculated
(5)
employed liquid
using using
the to mathematical
and calculate
gas, is calculated expression
theEq.sloshing
(5)
using for the swayexpression
is the mathematical excitation, forwhile Eq. (6) is the m
the sway
oshing
iquid
is employed
h isofemployed inphenomenon
a tank. toIt is is
assumed
calculate
to calculate represented that
the the by
sloshing
the sloshing the
fluid Reynolds-averaged
flow is incompressible, Navier-Stokes
and the (RANS)
turbulent
me non
bulent
duid2): fraction
.Reynolds
ing(2):
flowIn
liquid
flow
phenomenon
flow is
this is stresses
in
equation
field in
incompressible,
study,
a
incompressible,
tank.
ain
isand the It
Eq.
represented
is
RANS
the and
conventional
assumed
(2):
turbulent
thethe are
by
that
calculated
the
turbulent flow
numerical
the fluid
fieldexpression
Reynolds-averaged based
in flow
theonRANS
approach
is incompressible,
for are
isNavier-Stokes the
employed roll to
calculated (2) and excitation,
the
excitation.
based
(RANS)
turbulent
calculate on the while
Here, 𝑋𝑋 and 𝜃𝜃 represent the lateral displacement and th
Eq.
sloshing
(6) is the mathematical expression
is and turbulent
for using
the roll excitation. X and θ of
theHere, represent the lateral
he sloshing phenomenon represented by the Reynolds-averaged Navier-Stokes (RANS)
neynolds
-ε formodel. the stresses
eynolds-averaged
phenomenon
Reynolds-averaged twoThe 1 phases,
[
volume
of and
𝜕𝜕 liquid
(𝛼𝛼
liquid
the
Navier-Stokes
𝜌𝜌
fraction
in
Navier-Stokes
) +
a
𝑛𝑛 𝑛𝑛
andfor
turbulent
tank.
𝛻𝛻 ∙ (𝛼𝛼
gas,
It the
flow
(RANS)
is
𝜌𝜌
is
assumed
(RANS)
𝑢𝑢
⃗
calculated
two
)
field
=
phases,
𝑆𝑆
inthat theusing
displacement,
+ ∑
liquid
the RANS fluid
(𝑚𝑚
and
̇
are
flow gas,
𝑛𝑛in the RANS are calculated based
−
respectively.
calculated
is
𝑚𝑚 ̇
is calculated
incompressible,
) ]
based 𝑋𝑋
on represents
and
0 (2) the turbulent amplitude the lateral displacement, 𝜃𝜃
+𝛻𝛻 The
waction
model.
ield ∙
𝛻𝛻(𝛼𝛼∙inin 𝜌𝜌
Reynolds
(𝛼𝛼
𝑞𝑞equation
the 𝑞𝑞𝑞𝑞
The where
𝑢𝑢
⃗
𝜌𝜌 𝑢𝑢
⃗) =
)
𝑞𝑞𝑞𝑞sloshing
𝑞𝑞RANS the
𝑆𝑆
stresses
=
𝑞𝑞𝜕𝜕𝜕𝜕in 𝛼𝛼
volume Eq.𝑆𝑆
are𝑞𝑞variables
𝑞𝑞𝛼𝛼
+ 𝑞𝑞
(2):
fraction ∑
+ and ∑
𝑝𝑝=1
calculated (𝑚𝑚
the
𝑝𝑝=1 ̇
(𝑚𝑚
for𝑝𝑝𝑝𝑝 ̇
𝑞𝑞
the and
−
turbulent
𝑝𝑝𝑝𝑝
based 𝑞𝑞−𝑚𝑚
two ̇
𝑞𝑞 𝑚𝑚
𝑞𝑞𝑞𝑞̇ ) ]
flow
𝑞𝑞𝑞𝑞
onphases, represent
) ]𝛼𝛼 field
𝑞𝑞 liquid 𝑝𝑝=1 the
and mass (2)
𝑝𝑝𝑝𝑝 (2) transfer
𝑞𝑞𝑞𝑞
gas, is calculatedNavier-Stokes displacement
using
on and the rotational displacement, respectively.
flow
field
ard k-ε model. in the
the 𝜌𝜌 RANSThe volume phenomenon
are
𝑞𝑞 calculated
pistostresses fraction is represented
based
forphases on
the two qthe
by the
pamplitude
phases, Reynolds-averaged
liquid and of[9,gas,the rotational
is calculated displacement, (RANS)
using the amplitude and 𝜔𝜔ofisthe thelateral
frequency of the excitations.
ion
phases,
oequation.equation
phases, 1 𝜕𝜕
from
liquid The
liquid phases
inand Eq.
Reynolds
and 𝑛𝑛
(2):
gas, gas, isq(2):and from
calculated
calculated andusing theusing turbulent to ,flow
respectivelyfield in the 10].
RANS are X0 represents
calculated based on displacement,
me𝜌𝜌 𝑞𝑞 𝑢𝑢
⃗ 𝑞𝑞standard
variables
present
represent
the ) =the
fraction [𝑆𝑆the 𝑚𝑚
𝛼𝛼mass +
equation
̇ 𝑝𝑝𝑝𝑝
(𝛼𝛼 ∑
𝜌𝜌
is𝑞𝑞and
mass the
transfer
𝑞𝑞𝑝𝑝=1 𝑚𝑚in(𝑚𝑚
)transfer
+ source 𝛻𝛻̇ 𝑝𝑝𝑝𝑝
̇ 𝑞𝑞𝑞𝑞 Eq. from −term
𝑞𝑞 𝑚𝑚
∙represent
(𝛼𝛼from 𝑞𝑞̇ 𝑞𝑞𝑞𝑞
𝜌𝜌phases 𝑢𝑢 𝑞𝑞)the
⃗phases
for)] =phase 𝑆𝑆to
𝑝𝑝mass 𝑝𝑝𝛼𝛼𝑞𝑞to 𝑞𝑞∑
𝑞𝑞+transfer
qand in
and
𝑛𝑛
athe
from
𝑝𝑝=1 (𝑚𝑚
cell
(2)
from froṁ phases
[9,
𝑝𝑝𝑝𝑝 − 10].
phases
phases 𝑞𝑞)𝑝𝑝
𝑚𝑚̇𝑞𝑞𝑞𝑞𝑞𝑞Eq. to 𝑞𝑞isand
]liquid
(3) andfrom θgas, represents
phases
(2) 𝑞𝑞 the using
amplitude of the rotational displacement,
1 𝜌𝜌𝑞𝑞𝜕𝜕 𝜕𝜕𝜕𝜕 𝑞𝑞 k-ε model. The volume fraction 𝑛𝑛
for two phases, 0 is calculated
ectively
eource [
term [9, (𝛼𝛼
the 1
10].
𝑞𝑞
for 𝜌𝜌 𝜕𝜕
𝑞𝑞 )
governing
𝑆𝑆 + 𝛻𝛻𝑞𝑞 in
is ∙the(𝛼𝛼
𝑞𝑞 in 𝜌𝜌
equation
source
𝑞𝑞 𝑞𝑞 𝑢𝑢
⃗ 𝑞𝑞 ) =for
term 𝑆𝑆 the
𝛼𝛼 for + ∑
momentum,
phase 𝑝𝑝=1 (𝑚𝑚
𝑞𝑞 iṅ 𝑝𝑝𝑝𝑝
𝑛𝑛 a − as
cell 𝑚𝑚 ̇ follows:
[9,
𝑞𝑞𝑞𝑞 ) ]
10]. Eq. (3) is and
the (2) ω is
governing the frequency of the excitations.
the
source volume
𝜌𝜌𝑛𝑛
∑mass 𝑞𝑞 𝜕𝜕𝜕𝜕
term [ phase
fraction
for phase
(𝛼𝛼
𝛼𝛼𝑞𝑞 𝜌𝜌 equation a cell
a cell
𝑞𝑞)]𝑞𝑞 ) + 𝛻𝛻 ∙ (𝛼𝛼𝑞𝑞 𝜌𝜌𝑞𝑞 𝑢𝑢
in[9, Eq. [9,10]. ⃗ 𝑞𝑞𝑞𝑞Eq.
(2):
10]. ) Eq. =(3) 𝑆𝑆𝛼𝛼(3) is+isthe ∑the governing
governing
𝑝𝑝=1(𝑚𝑚̇𝑝𝑝𝑝𝑝 − 𝑚𝑚̇𝑞𝑞𝑞𝑞 )] (2)
+ the
ables
𝛼𝛼for + ∑
𝑝𝑝=1 𝑚𝑚𝑛𝑛̇ (𝑚𝑚
𝑝𝑝𝑝𝑝
̇
(𝑚𝑚 𝜌𝜌𝑞𝑞̇ −
transfer
and
𝑝𝑝𝑝𝑝 𝑚𝑚
𝜕𝜕𝜕𝜕
𝑝𝑝𝑝𝑝
𝑚𝑚
− ̇ ̇
from
𝑞𝑞𝑞𝑞 𝑚𝑚
𝑞𝑞𝑞𝑞 ̇ represent
𝑞𝑞𝑞𝑞 ) ]
phases 𝑝𝑝
the to mass𝑞𝑞 (2)
and from
transfer
(2) 𝑞𝑞phases
from 𝑞𝑞
phases 𝑝𝑝 to 𝑞𝑞 and from Fig.
phases 𝑞𝑞
2 shows the two-dimensional geometric parameters
ws:𝑞𝑞 the𝑝𝑝=1
lows: momentum, as1 follows: 𝜕𝜕 Fig. 2. Parameters for a two-dimensional rectangular tank (left) and a two-dimensional prismatic tank (right).
∑ 𝑛𝑛
ermles
ely 𝑚𝑚
for
[9, ̇ 𝑝𝑝𝑝𝑝10]. and𝑆𝑆𝛼𝛼𝑚𝑚𝑞𝑞̇ 𝑞𝑞𝑞𝑞
phase in
is represent
a
the cell [[9,(𝛼𝛼
source the
10]. 𝑞𝑞
term 𝜌𝜌massEq.
𝑞𝑞 )for +(3) 𝛻𝛻phase∙is(𝛼𝛼
transfer the𝑞𝑞𝑞𝑞𝜌𝜌from 𝑢𝑢
⃗ a𝑞𝑞 )cell
governing
𝑞𝑞
in = 𝑆𝑆[9,
phases 𝛼𝛼 +
𝑝𝑝
10]. to 𝑞𝑞
𝑝𝑝=1
Eq. (𝑚𝑚
and (3) ̇ from
𝑝𝑝𝑝𝑝
is −
the 𝑚𝑚 ̇
phasesof
𝑞𝑞𝑞𝑞
governing)a] 𝑞𝑞
rectangular tank
(2) and a prismatic model tank. prismatic
The tank (right).
e variables 𝑚𝑚̇𝑝𝑝𝑝𝑝𝑞𝑞 and 𝑚𝑚
𝜌𝜌 ̇ 𝑞𝑞 𝜕𝜕𝜕𝜕
represent
𝜕𝜕 to 𝑞𝑞𝑞𝑞𝑞𝑞and from phases 𝑞𝑞 the mass transfer from 𝑞𝑞
phases Fig. 𝑝𝑝
2. to 𝑞𝑞
Parameters and for
from a phases 𝑞𝑞
two-dimensional rectangular tank (left) and a two-dimensional
+sfer
𝛻𝛻ectively from
(𝜌𝜌𝑢𝑢
∙momentum,
𝛻𝛻 (𝜌𝜌𝑢𝑢
∙ ⃗ 𝑢𝑢 )
⃗⃗ [9, 𝑆𝑆
phases
⃗=
𝑢𝑢 )
𝛼𝛼𝑞𝑞= −𝛻𝛻𝛻𝛻−𝛻𝛻𝛻𝛻 𝑝𝑝 ++ 𝑝𝑝
(𝜌𝜌𝑢𝑢 𝛻𝛻 ⃗ ) 𝑞𝑞
+
(𝜏𝜏𝜏𝜏 ) 𝛻𝛻
(𝜏𝜏𝜏𝜏
∙the∙ source )
++ ∙ (𝜌𝜌𝑢𝑢
𝜌𝜌𝑔𝑔𝜌𝜌𝑔𝑔 ⃗ ++𝑢𝑢
⃗ ) = −𝛻𝛻𝛻𝛻𝑞𝑞 𝑞𝑞 +
𝐹𝐹 𝐹𝐹for phase 𝑞𝑞 in a cell (3) 𝛻𝛻 ∙ (𝜏𝜏𝜏𝜏 ) + 𝜌𝜌𝑔𝑔 + 𝐹𝐹 rectangular tank is considered for the validation of modelling
is 𝛻𝛻
yansfer[9, 10]. from phases
is the sourceto and
term from for phases
phase in a cell [9, 10]. Eq. (3) is the
(3) governing (3)
he 10]. as𝑆𝑆𝜕𝜕𝜕𝜕𝛼𝛼 follows: term Fig.transfer 2 shows [9,(3) 10]. Eq. the (3) is𝑝𝑝 the
to 𝑞𝑞governing
two-dimensional geometric parameters
𝑞𝑞where
se 𝑞𝑞 in
in
momentum, a cell athe cell variables
[9, [9,
as 10]. 10].
follows: 𝑞𝑞 𝑚𝑚̇(3) and
Eq. Eq. 𝑝𝑝𝑝𝑝(3) is theis𝑚𝑚 the ̇ governing represent
𝑞𝑞𝑞𝑞governing the mass from phases assumptions
and fromused phases in 𝑞𝑞the computation, andof theamembrane-
rectangular tank and
𝜕𝜕 Fig. 2 shows the two-dimensional geometric parameters of a rectangular
for
𝑢𝑢
⃗,or,)𝜌𝜌=
denotes
to the
𝜌𝜌 −𝛻𝛻𝛻𝛻
𝑝𝑝,
the the momentum,
the where
density
respectively +velocity
density ∙of(𝜏𝜏𝜏𝜏
𝛻𝛻(𝜌𝜌𝑢𝑢 ⃗ of )[9,
the
𝛼𝛼𝑞𝑞
as
)denotes
++ follows:
vector,
the 𝜌𝜌𝑔𝑔
10]. (𝜌𝜌𝑢𝑢
𝛻𝛻mixture,
∙mixture, 𝜌𝜌𝐹𝐹𝑢𝑢
𝑆𝑆+⃗the )𝑝𝑝the
the
⃗ is = 𝑝𝑝density
velocity
the −𝛻𝛻𝛻𝛻
the static
source +vector,
static of
𝛻𝛻 model
(𝜏𝜏𝜏𝜏for
∙the
pressure,
term ρ(3)
)mixture,
pressure, + the
𝜌𝜌𝑔𝑔
phase tank.
and +
anddensity 𝑝𝑝the
𝜏𝜏𝜏𝜏𝑞𝑞𝐹𝐹 𝜏𝜏𝜏𝜏in
theThe
the ofrectangular
static
a cell the 10]. type
[9,pressure, Eq. (3) and
(3) is𝜏𝜏𝜏𝜏the
tank
tank isgoverning
with
the considered
a prismaticfor the isvalidation
shape the modelof tank tank
modelling
LNG assumptand
𝜕𝜕 𝜕𝜕𝜕𝜕
sor
present
epresent
equation[11,the 12].
the
for (𝜌𝜌𝑢𝑢
mixture, 𝜌𝜌𝑔𝑔 +
gravitationalandp𝛻𝛻the
⃗ ) 𝜕𝜕momentum,
gravitational
the ∙𝐹𝐹(𝜌𝜌𝑢𝑢
static
⃗ 𝑢𝑢
represent
body ⃗)=
body pressure,
force−𝛻𝛻𝛻𝛻
theand
force + and
gravitational
and (𝜏𝜏𝜏𝜏model
𝛻𝛻external )the
∙ external
the + 𝜌𝜌𝑔𝑔stress tank.
+
computation,
bodybody
body tensor
𝐹𝐹force,
forceThe
force, [11,
andrectangular
12].
and theused
external (3)tank theisstudy.
inforce,
membrane-type
body considered for the
tank with validationshape
a prismatic of modelling
is the model assumptiLNG
(𝜏𝜏𝜏𝜏 )
𝜕𝜕𝜕𝜕 (𝜌𝜌𝑢𝑢⃗ )𝜌𝜌𝑝𝑝+the 𝛻𝛻 ∙as(𝜌𝜌𝑢𝑢 follows:
⃗ )pressure,
⃗ 𝑢𝑢 = −𝛻𝛻𝛻𝛻 +and 𝛻𝛻 ∙𝜏𝜏𝜏𝜏(𝜏𝜏𝜏𝜏the
)𝑝𝑝+the 𝜌𝜌𝑔𝑔static
+ 𝐹𝐹 pressure, and 𝜏𝜏𝜏𝜏 the (3)
𝛻𝛻
+ ∙
density
es 𝛻𝛻 the
∙ +
(𝜏𝜏𝜏𝜏
or, 𝜏𝜏𝜏𝜏, 𝜏𝜏𝜏𝜏can
ely
nsor, [11, )
, canof 𝜌𝜌𝑔𝑔
velocity
+
12]. the
bebe +
and
𝜌𝜌𝑔𝑔
The 𝐹𝐹
mixture,
+vector,
𝐹𝐹
𝜕𝜕𝜕𝜕
expressedstress
expressed represent
tensor,
as:as: the the
static
𝜏𝜏𝜏𝜏
density
, gravitational
can be of
(3) the
expressed
(3) the body
mixture, force
computation,
as:
in the study. and external
and the To model
membrane-type the sway and
tank roll
with excitations
a prismatic of the
shapetank,is the
the model LNG
𝜌𝜌 𝜕𝜕density of the mixture, 𝑝𝑝 the static pressure, and 𝜏𝜏𝜏𝜏 the
thegravitational
the
[11, velocity
12]. body force,
vector,
𝜌𝜌𝑔𝑔velocity
and 𝐹𝐹body respectively
representthe
force (𝜌𝜌𝑢𝑢
𝜌𝜌and
the⃗ ) gravitational
+[11,
external ∙ (𝜌𝜌𝑢𝑢
𝛻𝛻 12]. ⃗of⃗ )in
The
𝑢𝑢
body stress
=
body−𝛻𝛻𝛻𝛻
force,
the tensor,
+ 𝛻𝛻𝑝𝑝and
force
study. (𝜏𝜏𝜏𝜏),external
∙ the +can𝜌𝜌𝑔𝑔be+ 𝐹𝐹
bodyfirstforce,
natural thefrequency
𝜏𝜏𝜏𝜏(4) (3) of the prismatic tank was used. The
denotes
ethe 𝜏𝜏𝜏𝜏expressed
mixture, = the
𝜏𝜏𝜏𝜏𝜌𝜌𝑔𝑔
=𝜇𝜇(𝛻𝛻𝑢𝑢 ⃗𝑝𝑝 +
𝑝𝑝𝜇𝜇(𝛻𝛻𝑢𝑢
the ⃗the
𝐹𝐹 +𝛻𝛻𝑢𝑢
static 𝑇𝑇vector,
⃗𝛻𝛻𝑢𝑢 𝑇𝑇
⃗) pressure, 𝜕𝜕𝜕𝜕
) pressure, the
𝜏𝜏𝜏𝜏anddensity
= 𝜇𝜇(𝛻𝛻𝑢𝑢 𝜏𝜏𝜏𝜏⃗ the
𝜏𝜏𝜏𝜏 the + 𝛻𝛻𝑢𝑢 ⃗ 𝑇𝑇 )mixture,
the (4)(4) static pressure, and
1,
an
11, mixture,
12].
be
12]. Theexpressed
and stress as: as:
static
represent
tensor, 𝜏𝜏𝜏𝜏,the be and
cangravitational
expressed as: To
body force model and the
external sway body andfirst
force, natural
roll frequencyofofthe
excitations thetank,
rectangular
the firsttank, ω1, canfrequency
natural be of th
nsor [11,𝑢𝑢 ⃗ 12]. 𝜌𝜌𝑔𝑔 and 𝐹𝐹 represent the 𝜌𝜌gravitational body force and 𝑝𝑝external body force, and 𝜏𝜏𝜏𝜏 the
where
nal
12].
tional body The
body denotes
force
stressforce and the
tensor,
and 𝜏𝜏𝜏𝜏
velocity
external ,
externalcan vector,
body
be force,
expressed
body force, the density
as: To
tankas: of
model the mixture,
the
was used. The first natural sway the determined
static
and pressure,
roll using
excitations Eq. (7),
theand its
tank, correlation
the first with the
of the rectangular tank, 𝜔𝜔1 , can be determine
of natural first
frequency of th
(4) frequency
s, 12].
11, the 12]. dynamic 𝑇𝑇 )viscosity [11, 12].
=ed
ely𝜇𝜇(𝛻𝛻𝑢𝑢
[11, ⃗ tensor
+12].𝛻𝛻𝑢𝑢⃗ The [11, 12]. 𝜌𝜌𝑔𝑔𝜏𝜏𝜏𝜏and
stress tensor, = 𝜏𝜏𝜏𝜏𝜇𝜇(𝛻𝛻𝑢𝑢
, can⃗ be+ 𝛻𝛻𝑢𝑢
𝐹𝐹 represent ⃗ 𝑇𝑇 )the
expressed (4)
ssed as:as:
stress gravitational body force andnatural
externalfrequency
body of
force,the prismatic tank, ω' , is expressed
𝜔𝜔1 , can in ′
𝜏𝜏𝜏𝜏 = 𝜇𝜇(𝛻𝛻𝑢𝑢 ⃗ + 𝛻𝛻𝑢𝑢 ⃗ 𝑇𝑇 ) tank (7), andwasits used. The firstwith
correlation natural frequency
(4)the first naturaloffrequency
the rectangular of the tank,
prismatic
1 tank,be𝜔𝜔determine
1 , is expr
𝑇𝑇respectively [11, 12]. The stress𝜏𝜏𝜏𝜏 tensor,
= 𝜇𝜇(𝛻𝛻𝑢𝑢 ⃗ 𝜏𝜏𝜏𝜏,+can 𝛻𝛻𝑢𝑢⃗ 𝑇𝑇 )expressed
be as: (4) Eq. (8) [13–15]:
(4) ′
dynamic
⃗) )
𝑢𝑢 𝑇𝑇 viscosity [11, 12]. (4)(4) (7), and its correlation
(8) [13–15]: Fig. 2. Parameters with for the first natural
a two-dimensional frequency
rectangular ofand
tank (left) thea two-dimensional
prismatic tank, 𝜔𝜔1tank
prismatic , is(right).
expre
ynamic viscosity [11, 12].
s the dynamic viscosity [11, 12]. 𝜏𝜏𝜏𝜏 = (8)𝜇𝜇(𝛻𝛻𝑢𝑢 ⃗ + 𝛻𝛻𝑢𝑢
[13–15]: ⃗ 𝑇𝑇 ) (4)
where μ is the dynamic viscosity [11, 12]. 𝜋𝜋𝜋𝜋 𝜋𝜋ℎ
where 𝜇𝜇 is the dynamic viscosity [11, 12]. Fig. 2 shows the two-dimensional geometric 𝜔𝜔12 = 𝜋𝜋𝜋𝜋 𝑡𝑡𝑡𝑡𝑡𝑡ℎ (𝜋𝜋ℎ ) of a rectangular
parameters (7) tank and
𝐵𝐵 𝐵𝐵
model tank. The rectangular tank is considered 𝜔𝜔12 = for 𝑡𝑡𝑡𝑡𝑡𝑡ℎ
the (validation) of modelling assumpti
𝐵𝐵 𝐵𝐵
EXCITATION FORCES the computation, and the membrane-type 𝜔𝜔1′
2
[𝐵𝐵1tank−1
𝐻𝐻1 𝑠𝑠𝑠𝑠𝑠𝑠ℎ(
𝜋𝜋𝐻𝐻1 2 (𝜋𝜋𝐵𝐵1 )]
)−𝐵𝐵1 𝐻𝐻1−1 𝑠𝑠𝑠𝑠𝑠𝑠shape
with𝜋𝜋𝐻𝐻𝐵𝐵a prismatic 𝐵𝐵 is the model LNG
𝜔𝜔 ′ 22 = 1 − [𝐵𝐵1 𝐻𝐻1−1 𝑠𝑠𝑠𝑠𝑠𝑠ℎ( 1 )−𝐵𝐵2𝜋𝜋ℎ
𝐻𝐻 −1 𝑠𝑠𝑠𝑠𝑠𝑠2 (𝜋𝜋𝐵𝐵1 )]
in the study. 𝜔𝜔11 𝜋𝜋𝜋𝜋𝜋𝜋𝜋𝜋ℎ(
𝐵𝐵 1 1
) 𝐵𝐵
=1− 𝐵𝐵
2𝜋𝜋ℎ (8)
𝜔𝜔1 2 𝜋𝜋𝜋𝜋𝜋𝜋𝜋𝜋ℎ( )
The tank can experience external perturbations
To model including
theEq’ns. 𝐵𝐵
sway and roll excitations of the tank, the first natural frequency of th
sway and roll excitations. Sway excitationIn this study,
induces a linear (5) to (8) are used to apply external excitations to the tank, in ord
FORCES tank
In was
this
transverse (lateral or side-to-side) prescribed
motion
used.
study,
along motionsThe
the Eq’ns.
first(5)natural
tank’saccording to
In(8) tofrequency
this are used
study,
specific ofapply
the rectangular
toformulae.
Eq’ns. (5) to external
(8) To
tank, 𝜔𝜔external
excitations
are accurately
used 1 ,to
to applysimulate
can
thebe determine
tank,
the in ord
motion
′
(7),
transverse axis, whereas roll excitation and
prescribed
under its motions
induces
these correlation
a tilting
excitations, with
according thetofirst
aexcitations
user-defined natural
to the
specifictank, frequency
in order
formulae.
function (UDF) toTo of
induce
and the prismatic
aprescribed
accurately
dynamic motions
simulate
mesh 𝜔𝜔1 , motion
tank,technique
the is expre
ar
rotation of the tank about its longitudinal
experience external perturbations including (8)
under axis, with
[13–15]:
thesesway
to impose the centre
excitations,
the and roll
external according
a user-defined
excitations.
excitation to specific
functionformulae. To accurately
Sway (UDF) and a dynamic mesh technique ar
effectively. simulate the
of the bottom of the tank acting as thetocentre
impose of the
thecoordinate
external motion ofeffectively.
the tank under these excitations, a user-defined
s a linearsystem.
transverse (lateral or side-to-side) motion along
If the sway and roll excitations are periodic, they can
the excitation
tank’s transverse axis,
function (UDF) and 𝜋𝜋𝜋𝜋 mesh technique𝜋𝜋ℎ
itation induces a tilting rotation of the tank about its longitudinal axis, with 𝜔𝜔athe
2dynamic
1 = 𝐵𝐵 𝑡𝑡𝑡𝑡𝑡𝑡ℎ ( 𝐵𝐵 )
are employed
be represented in the form of sinusoidal waves, as shown to impose the external excitation effectively.
om of the tank acting as the centre of the coordinate system. If The
below. the ′investigation
sway
2 and rollfocused on the effects of𝜋𝜋𝐵𝐵different
𝜋𝜋𝐻𝐻
[𝐵𝐵1 𝐻𝐻1−1 𝑠𝑠𝑠𝑠𝑠𝑠ℎ( 𝐵𝐵 1 )−𝐵𝐵1 𝐻𝐻1−1 𝑠𝑠𝑠𝑠𝑠𝑠2 ( 𝐵𝐵 1 )]
𝜔𝜔1 on the sloshing
excitations loads in the tank, including single
2 = 1−
riodic, they can be represented in the form of sinusoidal waves, as shown below. 2𝜋𝜋ℎ
𝜔𝜔1
excitations (pure sway and pure roll) and
𝜋𝜋𝜋𝜋𝜋𝜋𝜋𝜋ℎ( ) multiple coupled
𝐵𝐵
𝑋𝑋 = 𝑋𝑋0 𝑠𝑠𝑠𝑠𝑠𝑠 (𝜔𝜔𝜔𝜔) excitations '(coupled
(5)sway-roll)'.
(5) In particular, the sloshing
pressures
In this study, Eq’ns. (5) in the
to (8) are tanktoare
used analysed
apply when the
external excitationto
excitations of roll
the tank, in ord
𝜃𝜃 = 𝜃𝜃0 𝑠𝑠𝑠𝑠𝑠𝑠 (𝜔𝜔𝜔𝜔)
prescribed motions according to specific(6)formulae. To accurately simulate the motion
under these excitations, a user-defined function (UDF) and a dynamic mesh technique ar
mathematical expression for the sway excitation, while Eq. (6) is the mathematical
to impose the external excitation effectively.
e roll excitation. Here, 𝑋𝑋 and 𝜃𝜃 represent the lateral displacement and the rotational
spectively. 𝑋𝑋0 represents the amplitude of the lateral displacement, 𝜃𝜃0 represents
the rotational displacement, and 𝜔𝜔 is the frequency of the excitations.
Fig. 2. Parameters for a two-dimensional rectangular tank (left) and a two-dimensional prismatic tank (right)
Parameters for a two-dimensional rectangular tank (left) and a two-dimensional prismatic tank (right).
is coupled with the sway’s in the same-phase (0°) and the In the numerical simulation, the fluid motion was
opposite-phase (180°). modelled as an incompressible unsteady flow, with a no-slip
Fig. 3 shows the time-dependent lateral displacements of condition assumed at the walls. To track the free surface,
a tank subjected to a pure sway excitation having an amplitude which represents the interface between the water and air,
of 0.015 m, while others are the time-dependent rotational the VOF method was employed [7]. To express the turbulent
displacements due to pure roll excitations with amplitudes of flow of the liquid, the standard k-ε model was used [17]. To
1°, 3°, and 5° in the same-phase or the opposite-phase relative calculate the pressure on a computational mesh from the
to the sway’s. The frequency of the excitations was set to the velocity components, the semi-implicit method for pressure-
first natural frequency of the tank, denoted as 1.0 ω'1. link equations (SIMPLE) algorithm was used, in which the
momentum equations were
coupled with an iterative
procedure [18]. A compressive
scheme was used to calculate
the volume fraction for the
interface capturing problem [19].
Fig. 4 illustrates the base
case for the computational
mesh used in the simulation
of the rectangular tank, which
was formed of a combination
of triangular and rectangular
grids. When excitation is
applied to a tank, sloshing-
induced liquid forces often
impacts and shocks certain
zones on the upper-left and
upper-right walls of the tank,
and on the roof. Furthermore,
Fig. 3. Displacement of a tank due to pure sway and pure roll excitations. The the motion of the liquid near the point of impact exhibits
roll excitation is in-phase or out-of-phase with the sway excitation complex flow phenomena [20]. Thus, to accurately capture
the complex flow behaviour of the liquid, smaller triangular
grids (0.0025 m here in the base case) were employed near
VALIDATION the anticipated impact zones on the walls [20, 21]. Larger
rectangular grids (0.005 m in the base case) were used in the
remaining regions to reduce the computational time while
To assess the validity of the modelling assumptions used still maintaining reasonable accuracy. The entire mesh for
for the numerical simulation, a comparison was conducted the base case consisted of a total of 64,620 grid points. To
between the simulation results and experimental data obtained assess the numerical accuracy of the results in terms of the
from a well-known prior study conducted by Hinatsu [16]. sloshing loads in the tank with respect to the grid resolution,
The experimental data selected for the comparison consisted two cases were analysed: the base case and a coarse case of
of measurements of sloshing pressures in a rectangular tank 32,132 grid points, half the number of the base case. It is
under pure sway excitation. A schematic representation of important to note that in both cases, the number of grid points
the rectangular tank used in their experimental investigation significantly exceeded those used in a validation process
is presented in Fig. 4. The tank had a base length of 1.2 (7,200) performed by Rhee [12]. Furthermore, our study
m and height of 0.6 m. For the experiment, the tank was investigated the influence of the time step size on accuracy.
filled with water to 60% of its height from the floor, and the Three time steps were compared: a base time step of 0.0005 s,
remaining spaces were filled with air. Subsequently, the tank a large time step of 0.001 s (twice that of the base case), and an
was subjected to pure sway excitation characterised by an extra-large time step of 0.002 s (twice that of the large case).
amplitude of 0.015 m and a period of 1.404 s. During the
validation process, the simulation results were compared
to the experimental data pertaining to the sloshing impact
pressure at a specific location, denoted as P1 at the tank
wall (as indicated in Fig. 4). The primary objective of this
comparison was to evaluate the degree of agreement between
our simulation data and the experimental data regarding the
sloshing impact pressures at the designated point.
Fig. 7. Schematic of the prismatic model tank (left, actual LNG tank size: H =
Tab. 2. Average pressures at P1 for the prismatic model tank under pure sway
29.71 m, H1 = 4 m, H2 = 16.507 m, H3 = 9.203 m, B = 40.31 m, B1 = 4 m, B2 =
or roll excitations [Units: Pa]
9.203 m, θ1 = θ2 = 45°) and the base case of the computational mesh used for
the tank (right) [Units: m] Pure sway Pure roll
Frequency
Amplitude Amplitude
10. A. B. Desamala, V. Vijayan, A. Dasari, A. K. Dasmahapatra, 21. C. Hu and M. M. Kamra, “An unstructured mesh method
and T. K. Mandal, “Prediction of oil-water flow patterns, for numerical simulation of violent sloshing flows,” Journal
radial distribution of volume fraction, pressure and velocity of Hydrodynamics, vol. 32, no. 2, pp. 259-266, April 2020,
during separated flows in horizontal pipe,” Journal of doi: 10.1007/s42241-020-0019-z.
Hydrodynamics, vol. 28, pp. 658-668, August 2016, doi:
10.1016/S1001-6058(16)60670-4. 22. B. Godderidge, S. Turnock, M. Tan, and C. Earl, “An
investigation of multiphase CFD modelling of a lateral
11. M. L. Hosain, U. Sand, R. Bel Fdhila, “Numerical sloshing tank,” Computers & Fluids, vol. 38, no. 2, pp. 183-
investigation of liquid sloshing in carrier ship fuel tanks,” 193, February 2009, doi: 10.1016/j.compfluid.2007.11.007.
IFAC-PapersOnLine, vol. 51, no. 2, pp. 583-588, 2018, doi:
10.1016/j.ifacol.2018.03.098. 23. H. Akyildiz and E. Ünal, “Experimental investigation of
pressure distribution on a rectangular tank due to the liquid
12. S. H. Rhee, “Unstructured grid based Reynolds-averaged sloshing,” Ocean Engineering, vol. 32, no. 11-12, pp. 1503-
Navier-Stokes method for liquid tank sloshing,” Journal of 1516, 2005, doi: 10.1016/j.oceaneng.2004.11.006.
Fluids Engineering, vol. 127, no. 3, pp. 572-582, 2005, doi:
10.1115/1.1906267. 24. J. H. Jung, H. S. Yoon, and C. Y. Lee, “Effect of natural
frequency modes on sloshing phenomenon in a rectangular
13. O. M. Faltinsen and A. N. Timokha, “An adaptive tank,” International Journal of Naval Architecture and
multimodal approach to nonlinear sloshing in a rectangular Ocean Engineering, vol. 7, no. 3, pp. 580-594, May 2015,
tank,” Journal of Fluid Mechanics, vol. 432, pp. 167-200, doi: 10.1515/ijnaoe-2015-0041.
April 2001, doi: 10.1017/S0022112000003311.
25. S. Wu and Y. Ju, “Numerical study of the boil-off gas
14. O. M. Faltinsen and A. N. Timokha, “Analytically (BOG) generation characteristics in a type C independent
approximate natural sloshing modes and frequencies in liquefied natural gas (LNG) tank under sloshing excitation,”
two-dimensional tanks,” European Journal of Mechanics Energy, vol. 223, no. 15, pp. 1-19, May 2021, doi: 10.1016/j.
– B/Fluids, vol. 47, pp. 176-187, September-October 2014, energy.2021.120001.
doi: 10.1016/j.euromechflu.2014.01.005.
26. K. P. Thiagarajan, D. Rakshit, and N. Repalle, “The
15. N. Parthasarathty, H. Kim, Y. H. Choi, and Y. W. Lee, air-water sloshing problem: Fundamental analysis and
“A numerical study on sloshing impact loads in prismatic parametric studies on excitation and fill levels,” Ocean
tanks under forced horizontal motion,” Journal of the Engineering, vol. 38, no. 2-3, pp. 498-508, 2011, doi:
Korean Society of Marine Engineering, vol. 41, no. 2, pp. 10.1016/j.oceaneng.2010.11.019.
150-155, 2017, doi: 10.5916/jkosme.2017.41.2.150.
27. H. Kim, P. Nanjundan, J. Jeon, and Y. W. Lee, “Numerical
16. M. Hinatsu, “Experiments of two-phase flows for the estimation on applying air-trapping mechanism to suppress
joint research,” In: Proceedings of the SRI-TUHH mini- sloshing loads in a prismatic tank,” Journal of Mechanical
Workshop on Numerical Simulation of Two-Phase Flows, Science and Technology, vol. 34, no. 7, pp. 2895-2902, 2020,
Ship Research Institute, Tokyo, Japan, pp. 12-19, 2001. doi: 10.1007/s12206-020-0621-6.
17. Y. H. Chen, Y. F. Yue, Y. Zhang, R. P. Li, and X. Xu, 28. X. Yuan, Y. Su, and P. Xie, “Frequency characteristics
“Numerical investigation of vibration suppression for the of sloshing resonance in a three-dimensional shallow-
combined device of non-Newtonian fluids coupled elastic water rectangular tank,” Journal of Marine Science and
baffle,” Journal of Applied Fluid Mechanics, vol. 16, no. 3, Engineering, vol. 10, no. 11, pp. 1792-1804, November 2022,
pp. 591-602, 2023, doi: 10.47176/jafm.16.03.1311. doi: 10.3390/jmse10111792.
18. H. Jin, Y. Liu, H. Li, and Q. Fu, “Numerical analysis of the 29. H. Kim, M. K. Dey, N. Oshima, and Y. W. Lee, “Numerical
flow field in a sloshing tank with a horizontal perforated study on sloshing characteristics with Reynolds number
plate,” Journal of Ocean University of China, vol. 16, no. 4, variation in a rectangular tank,” Computation, vol. 6, no. 4,
pp. 575-584, 2017, doi: 10.1007/s11802-017-3369-6. pp. 53-63, October 2018, doi: 10.3390/computation6040053.
Jacek Nakielski
Gdansk University of Technology, Institute of Naval Architecture, Poland
Abstract
In recent years, issues related to the impact of human activity on the natural environment have become pressing, and the
challenge of global warming necessitates immediate action. To support environmental protection efforts, it has become
imperative to adopt a broader perspective when evaluating various products and systems. A valuable tool for such assessments
is a life cycle assessment (LCA), which enables a comprehensive analysis of the entire life cycle of a product.
This paper presents a comparative analysis of the hull of a fast patrol craft, fabricated using three different materials: steel,
aluminium, and composite materials. The LCA covers every stage from material production, through the construction
and use of the hull, to its eventual disposal. A specific criterion was established to evaluate the impact of the hull on the
environment, with clearly defined system boundaries.
In the final section, we draw some conclusions that underscore the importance of reusing construction materials.
By emphasising this approach, ecological footprints can be minimised and a sustainable future can be created.
Fig. 1. Phases of an LCA according to ISO 14040:2006 Environmental Fig. 2. Structure of the aluminium hull [8]
management—Life cycle assessment—Principles and framework [2]
Energy
Process Mass Units SEC Units Unit
consumption
Secondary steel
Primary aluminium
Secondary aluminium
DISCUSSION
1. J. Barreiro, S. Zaragoza, and V. Diaz-Casas, ‘Review of ship 13. F. Nicolae, C. Popa, and H. Beizadea, ‘Applications of life cycle
energy efficiency’, Ocean Engineering, Volume 257, 2022. doi: assessment (LCA) in shipping industry’, 14th International
10.1016/j.oceaneng.2022.111594. Multidisciplinary Scientific GeoConference & EXPO SGEM2014,
Section Air Pollution and Climate Change, Albena, Bulgaria,
2. ISO 14040:2006 Environmental management - Life cycle 17-26 Jun 2014. doi: 10.5593/sgem2014/b42/s19.038.
assessment - Principles and framework.
14. Jianqiao Wu, Ying Pan, Ziyi Ruan, Ziji Zhao, Jing Ai, Jinghai
3. S.D. Chatzinikolaou and N.P. Ventikos, ‘Assessing environmental Ban, and Xianghai Jing, ‘Carbon fiber-reinforced epoxy with
impacts of ships from a life cycle perspective’, in 2nd International 100% fiber recycling by transesterification reactions’, Frontiers
Conference on Maritime Technology and Engineering, MARTECH in Materials Sec. Polymeric and Composite Materials, Volume
2014, Lisbon, Portugal, 15-17 October 2014. doi.org/10.1201/ 9, 2022. doi: 10.3389/fmats.2022.1045372.
b17494
15. M. M. Rabby, M. Rahman, P. P. Das, V. Vadlamudi and
4. L. Joosten, ‘The Industrial Metabolism of Plastics. Analysis R. Rassel, ‘Carbon fibers recycling from degraded prepress and
of Material Flows, Energy Consumption and CO2 Emissions mechanical properties of recycled composite’, North America
in the Lifecycle of Plastics’, PhD [Dissertation], Utrecht Society for the Advancement of Material and Process Engineering
University, 2001. (SAMPE), Seattle, United States, 17 April 2023. doi: 10.33599/
nasampe/s.23.0267.
5. C. Favi, M. Germani, F. Campi, M. Mandolini, S. Manieri,
M. Marconi, and A. Vita, ‘Life cycle model and metrics in
shipbuilding: How to use them in the preliminary design
phases’, in 25th CIRP Life Cycle Engineering (LCE) Conference,
Copenhagen, Denmark, 30 April – 2 May 2018. doi: 10.1016/j.
procir.2017.11.071
Ali Zinati
Department of Maritime Engineering, Amirkabir University of Technology, Tehran, Islamic Republic of Iran
Mohammad Javad Ketabdari
Department of Maritime Engineering, Amirkabir University of Technology, Tehran, Islamic Republic of Iran
Hamid Zeraatgar *
Faculty of Marine Technology, Amirkabir University of Technology, Tehran, Islamic Republic of Iran
Abstract
Propeller performance is typically considered under clean conditions, despite the fact that fouling is an inevitable
phenomenon for propellers. The main objective of this study is to investigate the effects of roughness due to fouling on
the performance of a propeller using a CFD simulation in conjunction with the roughness function model. A simulation
of a clean propeller is verified for a five-blade propeller model using existing experimental results. A roughness
function model is then suggested based on existing measured roughness data. The simulations are extended for the
same propeller under varying severities of roughness. Initially, it is concluded that KT and ηo gradually decrease with
increasing fouling roughness, while KQ increases, compared to smooth propeller. For instance, at J=1.2 for medium
calcareous fouling, KT is reduced by about 26%, KQ increases by about 7.0%, and ηo decreases by 30.9%. In addition, for
the rough propeller, the extra power required is defined as the specific sea margin (SSM) to compensate for the power
loss. A slight roughness causes a large decrease in ηo. A propeller painted with foul-release paint and an unpainted
propeller are found to require 2.7% SSM and 57.8% SSM over four years of service, respectively. Finally, the use of
foul-release paints for propeller painting is strongly advised.
Keywords: Propeller performances, Blade roughness, Frictional resistance, CFD simulation, Fouling
3
NUMERICAL SIMULATION the impact of roughness on ΔU+ depends on the type and
coverage of the roughness.
In this section, the governing equations are introduced, In this study, the roughness parameters measured by
the roughness function is presented, and the numerical Schultz and Flack [29] are used to develop a new roughness
simulations of smooth and roughed propellers are explained. function model. The proposed model is fitted to the roughness
function values reported by Schultz and Flack [29]. One
MATHEMATICAL FORMULATION advantage of our roughness function model is that it can be
MATHEMATICAL FORMULATION applied to all types of fouled surfaces and typical antifouling
The governing equations in this context are those of mass coatings. The values of the sand grain roughness height and
MATHEMATICAL FORMULATION
and momentum conservation, which for compressible flows other relevant data are shown in Table 1 for a range of surface
in the Cartesian coordinate system are as follows [23]: conditions.
The governing equations in this context are those of mass and momentum conservation, which for
compressible
MATHEMATICAL the Cartesian coordinate system are as followsTab.
flows inFORMULATION
MATHEMATICAL
1. Representative coating and fouling conditions [23]
[23]:
The governing equationsFORMULATION
in this context are those of mass and momentum conservation, which for NSTM
CAL FORMULATION Description of condition ks (µm) Rt50 (µm)
𝜕𝜕(𝜌𝜌𝑢𝑢̅𝑖𝑖 ) coordinate system are as follows [23]:
compressible flows in the Cartesian rating*
(4)
(4) =0
The governing equations in𝜕𝜕(𝜌𝜌𝑢𝑢 𝜕𝜕𝜕𝜕context
𝑖𝑖 Hydraulically smoothwhich surfacefor 0. 0 0
𝜕𝜕(𝜌𝜌𝑢𝑢 ̅𝑖𝑖are) those 𝜕𝜕 inofthis
this ̅𝑖𝑖 ) are those of mass and momentum
𝜕𝜕𝜏𝜏̅̅̅
𝜕𝜕𝑝𝑝̅conservation,
conservation,
(4)
The
quations governing
in this context
compressible equations flows + in (𝜌𝜌𝑢𝑢̅𝑖𝑖 𝑢𝑢𝜕𝜕𝜕𝜕
the mass
Cartesian context =
̅̅̅̅̅̅
̅𝑗𝑗 + 𝜌𝜌𝑢𝑢𝑖𝑖 𝑢𝑢𝑗𝑗 ) = −
and ′ 0
momentum
coordinate ′ are thosesystem +
of are mass
𝑖𝑖𝑖𝑖
as and
which
follows momentum
for
[23]: conservation, which (5)
for
𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕𝑖𝑖system 𝑖𝑖 𝜕𝜕𝜕𝜕𝑖𝑖 system 𝜕𝜕𝜕𝜕𝑗𝑗 are asDescription Typical asofapplied condition AF coating 0
NSTM rating* 30 150 ks (µm) Rt50
compressible
ws in the Cartesian flows
𝜕𝜕(𝜌𝜌𝑢𝑢 ̅𝑖𝑖 ) in the
𝑖𝑖
coordinate
𝜕𝜕 Cartesian are
𝜕𝜕(𝜌𝜌𝑢𝑢̅𝑖𝑖 )
as followscoordinate [23]:
𝜕𝜕𝑝𝑝̅ 𝜕𝜕𝜏𝜏̅̅̅ follows [23]:
𝑖𝑖𝑖𝑖 HydraulicallyDeteriorated coating or light slime 10–20 (5) 100 300
where 𝜕𝜕(𝜌𝜌𝑢𝑢 is ̅the density, + 𝜌𝜌𝑢𝑢
(𝜌𝜌𝑢𝑢 ̅̅̅̅̅̅
′
̅𝑖𝑖𝑖𝑖𝑢𝑢 ′
𝑢𝑢̅𝑗𝑗𝑗𝑗𝜕𝜕𝜕𝜕+
is𝜌𝜌𝑢𝑢the ̅̅̅̅̅̅′ ′
=𝑖𝑖0𝑢𝑢𝑗𝑗 )Reynolds =− + stress, 𝑢𝑢 ̅ is (4) the
smooth surface
averaged (4)
Cartesian component
0. 0
𝜕𝜕𝜕𝜕𝑖𝑖 𝜕𝜕𝜕𝜕𝑗𝑗 Typical𝑖𝑖 as applied
(5)
𝑖𝑖 ) 𝜕𝜕𝜕𝜕
= 0𝑖𝑖𝜕𝜕(𝜌𝜌𝑢𝑢̅ )𝜕𝜕𝜕𝜕𝑖𝑖 𝜕𝜕 𝜕𝜕(𝜌𝜌𝑢𝑢 Heavy AF coating
slime 30 0 300 600 30 1
̅𝑖𝑖𝑖𝑖 ) 𝜕𝜕𝑝𝑝̅ 𝜕𝜕𝜏𝜏
̅̅̅ (4)
velocity, and 𝜕𝜕𝜕𝜕𝑖𝑖 p is the mean 𝑖𝑖
̅̅̅̅̅̅
pressure. = 𝜏𝜏
̅̅̅
0 are the 𝑖𝑖𝑖𝑖
mean Deteriorated
viscous coating
stress or
tensorlight slime
components, (5) as follows: 10–20 100 3
𝜕𝜕 is the ′ density,
+ 𝜌𝜌𝑢𝑢 ′ 𝑢𝑢 𝑖𝑖′𝑢𝑢
̅̅̅̅̅̅ Small calcareous fouling or weed 40–60 1,000 1,000
where
𝑢𝑢̅𝑖𝑖 ) 𝜕𝜕𝜕𝜕 𝑖𝑖 𝜕𝜕𝑝𝑝̅ 𝜕𝜕𝜕𝜕 𝑖𝑖
(𝜌𝜌𝑢𝑢
𝜕𝜕𝜏𝜏
̅̅̅
𝑖𝑖
𝑖𝑖𝑖𝑖
̅𝜕𝜕𝜕𝜕
𝑗𝑗
̅𝑗𝑗is + 𝜌𝜌𝑢𝑢
𝑖𝑖 the ′ ′
𝑖𝑖 𝑢𝑢𝑗𝑗 ) =−
𝑖𝑖𝑖𝑖Reynolds
𝜕𝜕𝜕𝜕 𝑖𝑖
+
𝜕𝜕𝜕𝜕 𝑗𝑗
stress, Heavy 𝑢𝑢
̅ 𝑖𝑖 is
slime the
(5) averaged Cartesian component 30 300 6
+ (𝜌𝜌𝑢𝑢̅𝑖𝑖 𝑢𝑢̅𝑗𝑗 + 𝜕𝜕(𝜌𝜌𝑢𝑢
𝜌𝜌𝑢𝑢̅̅̅̅̅̅ ′) = −
𝑖𝑖 𝑢𝑢𝑗𝑗̅ Medium calcareous
or weed fouling 70–8040–60 3,000 3,000 1,000
where ρ is the ) density, 𝜕𝜕𝑖𝑖 + 𝜕𝜕𝜕𝜕 𝜌𝜌𝑢𝑢̅̅̅̅̅̅
′ 𝑢𝑢 ′ is is the Reynolds 𝜕𝜕𝑝𝑝̅stress, 𝜕𝜕𝜏𝜏
̅̅̅ Small calcareous
is the
fouling 1,0
𝑖𝑖 𝑖𝑖′̅̅̅
𝜏𝜏𝑢𝑢𝑖𝑖𝑖𝑖
𝜕𝜕
velocity, and
𝑖𝑖 𝜕𝜕𝜕𝜕𝑖𝑖 where is the
p is the+mean(𝜌𝜌𝑢𝑢 𝑖𝑖 density,𝜕𝜕𝜕𝜕 pressure.
𝑗𝑗
̅𝑖𝑖 𝑢𝑢 𝑖𝑖̅ 𝑗𝑗+𝜕𝜕𝑢𝑢 ̅̅̅̅̅̅̅
the
𝜌𝜌𝑢𝑢 ′𝜕𝜕𝑢𝑢 ̅
Reynoldsare
𝑗𝑗 −the mean
stress, 𝑢𝑢
̅
𝑖𝑖𝑖𝑖 is
viscous the averaged
Description
stress Cartesian
of
tensor condition component
components, NSTM
as 90–100 (5)
rating*
follows: k s (µm) Rt50 (µm)
𝑗𝑗 ) = + 𝑖𝑖 Medium calcareous
̅̅̅̅̅̅ 𝜏𝜏stress,
̅̅̅ = 𝑗𝑗𝜇𝜇(𝑢𝑢̅𝑖𝑖 isvelocity, Hydraulically Heavy calcareous
smoothfouling fouling 0. (6) 70–80 10,000 0 10,0003,000 0 3,0
𝜏𝜏+ )the
𝜌𝜌𝑢𝑢 ′ ′ 𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕mean 𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕
𝑖𝑖 𝑢𝑢𝑗𝑗 is and
surface
density, velocity, averaged the CartesianReynolds
𝑖𝑖p is the 𝑖𝑖 component
𝑖𝑖𝑖𝑖 pressure.
𝜕𝜕𝜕𝜕𝑗𝑗 ̅̅̅
the are
𝑖𝑖𝑖𝑖𝜕𝜕𝜕𝜕
averaged and mean 𝑖𝑖pCartesian
isviscous
the mean
𝑗𝑗 Heavy component
stress
Typical astensor
calcareous*
applied NSTM components,
fouling
AF (2002):
coating Navalas follows:
Ships’ Technical Manual
0 90–100 30 10,000 150 10,
𝑖𝑖
̅̅̅̅̅̅
where
the meanpressure.
where pressure.
is the thedynamic 𝜏𝜏density,
̅̅̅𝑖𝑖𝑖𝑖 are arethe themean mean
𝜌𝜌𝑢𝑢 ′ 𝑢𝑢viscous′viscous 𝜕𝜕𝑢𝑢̅𝑖𝑖 stress stress
𝜕𝜕𝑢𝑢 ̅𝑗𝑗tensor tensor components,
components, 𝑢𝑢̅
as follows:
𝑗𝑗𝜇𝜇(is𝜕𝜕𝑢𝑢̅the Reynolds stress, is(2002):
the averaged Cartesian
* Deteriorated
NSTM coatingNaval
or light Ships'
slimeTechnical Manual component 10–20 (6) 100 300
𝜇𝜇asis follows: 𝜏𝜏
̅̅̅
viscosity. 𝑖𝑖=
𝑖𝑖𝑖𝑖̅̅̅
𝜏𝜏𝑖𝑖𝑖𝑖 = 𝜇𝜇(𝜕𝜕𝜕𝜕𝑗𝑗 + 𝜕𝜕𝜕𝜕)𝑖𝑖
𝑖𝑖 + 𝜕𝜕𝑢𝑢 ̅ 𝑗𝑗 ) Heavy 𝑖𝑖 slime
The proposed roughness (6) function 30
model is as follows:
300 600
velocity, and 𝜕𝜕𝑢𝑢
p ̅𝑖𝑖is 𝜕𝜕𝑢𝑢 the ̅𝑗𝑗 mean pressure. 𝜕𝜕𝜕𝜕𝑗𝑗 ̅̅̅ 𝜏𝜏 𝜕𝜕𝜕𝜕
𝑖𝑖𝑖𝑖 𝑖𝑖 are the mean viscous
The
Small calcareous
proposed stress
(6) tensor
fouling
roughness
or weedcomponents, as 40–60
function model is as70–80 follows:
follows:
1,000 1,000
𝜏𝜏𝑖𝑖𝑖𝑖 = 𝜇𝜇(
̅̅̅ + )
where 𝜇𝜇where
The Reynolds-averaged is the 𝜇𝜇𝜕𝜕𝜕𝜕is𝑗𝑗dynamic
the 𝜕𝜕𝜕𝜕dynamic
𝑖𝑖 viscosity.
Navier-Stokes
viscosity. (RANS) method
Medium calcareous
is employed
Heavy calcareous
fouling
fouling to solve the governing 90–100
3,000
10,000
3,000
10,000
+
ynamic viscosity. 𝜕𝜕𝑢𝑢̅𝑖𝑖 𝜕𝜕𝑢𝑢 ̅𝑗𝑗 0 Technical Manual 𝑘𝑘𝑠𝑠 ≤ 2.5
equations The using Reynolds-averaged the commercial = 𝜇𝜇( CFD
𝜏𝜏Navier-Stokes
̅̅̅ + software )(RANS) STAR-CCM+. * NSTM (2002): Naval Ships'
The flow variables are discretised
+ (6)
The Reynolds-averaged 𝑖𝑖𝑖𝑖 Navier-Stokes (RANS) method method is (6)employed
is employedto solve 1 the
to governing
solve 𝜋𝜋 𝑙𝑙𝑙𝑙𝑙𝑙(𝑘𝑘𝑠𝑠 ⁄2.5)
the governing
in space
veraged using using
Navier-Stokes
equations second-order (RANS)
the commercial methodschemes, 𝜕𝜕𝜕𝜕 𝑗𝑗is employed
CFD 𝜕𝜕𝜕𝜕and
software 𝑖𝑖 theSTAR-CCM+.
toconvection
solve the term
governing
The
The proposed Δ𝑈𝑈of
flow + 𝑙𝑙𝑙𝑙
=the𝜅𝜅 turbulence
variables
roughness (0.2667
are discretised
function model
𝑠𝑠𝑠𝑠𝑠𝑠 [
𝑘𝑘𝑠𝑠+ ) isis used2 𝑙𝑙𝑙𝑙𝑙𝑙 (10)
as follows:
]
as a first-2.5 < 𝑘𝑘𝑠𝑠+ < 25 (7)
equations
where 𝜇𝜇in is using
theCFD the
dynamic commercial viscosity. CFD software STAR-CCM+. The flow variables
10 is used as+a first- are discretised
the
order commercial
upwind space using
scheme. software
second-order The STAR-CCM+. shearschemes, stress The the
and flow
transport variables
convection (SST) are discretised
term ofk-the turbulence
turbulence model
𝑙𝑙𝑙𝑙 (0.2667
{𝜅𝜅 turbulence 𝑘𝑘𝑠𝑠 ) is applied to 𝑘𝑘the
+
𝑠𝑠 ≤ 2.5 𝑘𝑘 + ≥ 25
𝑠𝑠
in spaceorder
cond-order using
where schemes, μ second-order
upwind is the and dynamic
scheme.the convection schemes,
Theviscosity. shearterm stress ofand the
transportthe convection
turbulence (SST)is k- used term
turbulence
as a first- of themodel1 is applied 𝜋𝜋is
to used
the
𝑙𝑙𝑙𝑙𝑙𝑙(𝑘𝑘𝑠𝑠 ⁄2.5) as a first-
+
complete
The RANS
Reynolds-averaged equations, as
Navier-Stokes this combines (RANS) the advantages
method A brief of
explanation
is the
employed
Δ𝑈𝑈 +
= k-
of and
𝑙𝑙𝑙𝑙 the
(0.2667 k-ε + 𝑠𝑠𝑠𝑠𝑠𝑠 [ 2 𝑙𝑙𝑙𝑙𝑙𝑙 (10) ]
𝑠𝑠 ) turbulence
𝑘𝑘solve
development
to the of this models.
2.5 < 𝑘𝑘𝑠𝑠+ <
roughness
governing 25
function
(7)
model is presented in Ap
heme. The
order upwind
completeshear stress RANStransport
scheme. equations, The (SST) shear
as this turbulence
k-stress
combines transport model(SST)
the advantages is applied of the k- toturbulence
k- the
and k-ε 𝜅𝜅 turbulence model models.is applied to the
In
equations this Inmodel, The
using the
Reynolds-averaged
the k-the
k-commercial formulation
formulation Navier-Stokes
CFD isisthe applied
softwarek- and (RANS) to the
STAR-CCM+. A.
inner
method part ofA the
The
1 boundary
brief explanation + layer,
ofk-ε and
the development the k-ε of25
𝑘𝑘𝑠𝑠+ ≥ this roughness
equations,
complete as
this
RANSthismodel, combines the
equations, advantages
as this of
combines applied to the
the k-ε turbulence
inner
advantages part ofmodels. the
of boundary
the k-𝜅𝜅flow
{ 𝑙𝑙𝑙𝑙 (0.2667
layer,
and variables
𝑘𝑘and
k-ε 𝑠𝑠 ) the
turbulence are discretised
models.
k-
n space
eformulation is
formulation employed
formulation to
using second-order the
is toapplied to
free
the solve
free stream,
to the the
stream, inner
schemes, governing
which
which part of
givesgives
theequations better
boundary
better
and the convection using
predictions predictions
layer, A the
andof
brief
Fig. 1 term the function
of
k-ε
flow
explanation
shows of the
separation
of
a schematic
the model
flow
the is
and
development presented
separation adverseof this in
and Appendix
roughness adversefunction A.model is presented
functioninmodel
Appendix
In
he thisstream,
free model,
commercial
pressure which thegives
gradients.
k- To
CFD
formulation
software
better achieve predictions STAR-CCM+.
reliable
isofapplied the
results, The
flow an
toflow
thevariables
separation
appropriate
inner
A.and choice
part
adverse of
of the1isturbulence
Fig.
grids
boundary
shows
representation
crucial. a schematic
is used
layer, of the
and asproposed
a first-
the
representation
k-ε roughness
of the proposed
in Eq. (
pressure
order upwind gradients. To achieve reliable results, an appropriate
compares it choice
k- turbulence with the of grids
measurements is crucial. of Schultz and Flack [29]. Good agreement is ob
formulation
ts. To achieve are discretisedtoscheme.
reliable theresults, free in space
The
stream, shear
an appropriate using which stress
second-order
choice gives transport
of gridsbetterschemes,
(SST)
predictions
is crucial. betweenand the of
roughness
two. the flow model
functionseparation is applied
model and
in Eq.
Fig. 1 shows a schematic representation of the proposed roughness function model in Eq. (7), and
to the
adverse
(7), and compares it with
theof
complete PhysicalRANS modellingequations, of the as roughness this combines geometry is the notadvantages
practical in CFD,the duek- to itsand complexity. The
pressure
Physical
ng
the convection
gradients.
modelling
of thewall-function
roughness geometry Toof termachieve
the is
of the
roughness
not practical
turbulence
reliable in results,
geometry
CFD,
isdue used an
is
to
ascomplexity.
appropriate
not
its
first- ofitThe
a practical
compares the
choice measurements
innear
with CFD, grids
due
measurements k-εto ofturbulence
isitsSchultz
crucial.
ofcomplexity.
Schultz and andmodels.
Flack
Flack [29].
The [29].GoodGoodagreement
agreement is observed
approach is therefore applied to solve the flow equations the wall, rather than
In this
wall-function
proach
order
model,
isusingtherefore
upwind the
approach k-
applied to is
scheme. formulation
solve
The
therefore the flow
shear is stress applied
applied transportto to
nearsolve
the (SST) inner k-ϖ
between
theratherflowthan
theis
part ofobserved
two. the between
boundary
equations near the wall, 7,9 the
layer, two. and
rather the k-ε
than
turbulence-model equations upequations to the wall. the wall,
Physical
formulation
-model
turbulence
equations
modellingto the
up
model
to free
the
ofwall. isstream,
the applied
roughness towhich the complete geometry
gives RANS
better is not equations, practical of
predictions in CFD,the flow dueseparation
to its complexity. Schultz and
and
Flack (2007)
adverse The
using turbulence-model
as this combines the equations
advantages ofupthetok-ϖ theand wall. k-ε turbulence
7,9 6,9 Proposed roughness function model
wall-function
pressure gradients. approach To achieve is therefore reliable applied results, to an solve the flow choice
appropriate equations near the
of6,9 grids is Proposedwall, rather than
Schultz
crucial. and Flack (2007)
models. In this model, the k-ϖ formulation is applied to the roughness function model
using turbulence-model
PROPOSED ROUGHNESS FUNCTION equations up to the wall. 5,9
Physical modelling inner part of the boundary layer, and
of the roughness geometry is not practical in CFD,5,9due to its complexity. Thethe k-ε formulation to
OUGHNESS FUNCTION
the free stream, which gives better predictions of the flow 4,9
PROPOSED
wall-function ROUGHNESS
approach is therefore FUNCTION appliedTo toachieve
solve reliable the flow equations 4,9 near the wall, rather than
separation and adverse pressure gradients.
In general, roughness equations functions areup ∆U+is+ no
using
PROPOSED turbulence-model
results, ROUGHNESS an appropriate choice of obtained grids to the experimentally,
wall.
is crucial. since there (m/s) universal
3,9 roughness
hness functions
function aremodel obtained for all kinds ofFUNCTION
experimentally, roughsince surfaces. there It is should
no universal be noted that∆U
roughness the (m/s)
impact 3,9
of roughness on
Physical modeling of the roughness geometry is not
or all kinds ΔU+ofdepends rough surfaces. on the type It shouldand coverage be notedofthat thethe impact of roughness on
roughness. 2,9
In general, practical roughness in CFD,functions due to its complexity. are obtained Theexperimentally,
wall-function since there 2,9
is no universal roughness
the type and coverage of the roughness.
Inapproach
this study, is therefore
the roughness applied parameters to solvemeasured the flow equations by Schultz near
and Flack [29] are 1,9used
1,9 to develop a
function model
PROPOSED for all kinds ofFUNCTION
ROUGHNESS rough surfaces. It should be noted that the impact of roughness on
In general,
roughness the roughness
+ newparameters
wall,
roughness rather than
function
measured functions usingmodel.
by Schultz turbulence-model
are Theobtained
and proposedFlack [29] model equations
experimentally,
are used is fitted toup to the
develop since a there is
roughness no universal
function values roughness
ΔU dependsto the on
wall. the type and coverage of the roughness. 0,90,9
function model for all kinds of rough surfaces. It should be noted that the impactitof
unction reported
model. by
The Schultz
proposed and Flack
model [29].
is fitted One to advantage
the roughness of our roughness
function values function model is that canroughness on
ultz and be
Flack applied[29]. to
One all types
advantage of fouled
of our surfaces
roughness and typical
function antifouling
model is that coatings.
it can The values of the sand
ΔU + -0,1
In
types depends
this
ofgrainstudy,
PROPOSED
fouled on
thethe
roughness
surfaces roughness type
ROUGHNESS
height
and typicaland and coverage
parametersother
antifouling FUNCTION relevant of
measured
coatings. the
dataroughness.
are by
The shown
values Schultz in the
of Table and1 Flack
sand for a range
-0,1
[29]
1
1 of aresurfaceused to develop 10 a
10
100
100
Innew
heightgeneral,
roughness
conditions.
and roughness
other function
relevant functions
datamodel. are shown areThe obtained Table 1experimentally,
in proposed formodel a rangeisoffitted since
surface to there
the roughness is no universal function roughness
+
ks (log scale)
values
+ (log scale)
In this study, thefor roughness parameters measured byofSchultz and Flack [29] are used to k
is develop ona the values reported by Schultz and
s
function
reported model
by In Schultz general, alland kinds
roughness
Flack of [29]. rough functionsOne surfaces. advantage are It should obtained our
Fig.be1.roughness
noted
Comparison thatof the the
function impact
proposed model
roughness of roughness that it can
Fig. 1. Comparison of the proposedfunction roughness model with
function model with the values
new
ΔU + roughness function model. The proposed modelFig.
be applied depends experimentally,
to on all the types since
type ofTablethere
and
fouled isRepresentative
1.coverage nosurfacesuniversal of coatingroughness
the
and roughness.
typical
and foulingfunction is fitted
1. Comparison
antifouling
conditions
to the
coatings.
[23]
roughness
of the proposed Thereported
function
roughness
values Flack
by Schultz [29]
offunction
and
values
theFlack model
sand[29] with the values reported by Schul
reported
grain Table byRepresentative
model
1.
roughness Schultz
for allheight kindsand of
coating Flack
rough
and andother [29].
surfaces.
fouling One
conditions
relevant It advantage
should [23]
data be noted are ofThe our
that
shown roughness
proposedinroughness Table function1 formodel
function
model
a range is of
has a similar
that
Flack it [29]
surface can
form to the built-in wall function of STAR-
In
be this study,
applied to the roughness
all types parameters
of fouled measured
surfaces by Schultz
and 4typical CCM+, andand
antifouling Flack as[29]
coatings.
is employed arefunction
theThe
wall used of
values toofdevelop
the sand
STAR-CCM+. a The proposed model is introduced
conditions. 4 The proposed
toshown
the fitted
roughness
CFD simulation setup
function
via
model hasB a= similar
theacoefficients 0,
form to the built-in wall function of
C = 0.2667, R+smooth = 2.5, and R+rough = 25
new
grainroughness
roughnessfunction model.
height and Therelevant
other proposeddatamodel is
areCCM+, in to the
Table roughness
1 for function
range of values
surface
and is employed as the wall function of STAR-CCM+. The proposed model is intr
reported
conditions.by Schultz and Flack [29]. One advantage ofto our the CFDroughness
simulationfunction
setup via model is that Bit=can
the coefficients 0, C = 0.2667, R+smooth = 2.5, and R+ro
be applied to all typesTable
of fouled surfaces andcoating
1. Representative typicaland
antifouling POLISH
coatings.
fouling conditions MARITIME RESEARCH,5
The values of the sand
[23] No 4/2023 63
grain roughness height and other relevant data are shown in Table 1 for a range of surface 5
Table 1. Representative coating
4 and fouling conditions [23]
conditions.
The proposed roughness function model has a similar form are 2D. The computational domain consists of two parts: the
to the built-in wall function of STAR-CCM+, and is employed inner region, which rotates with the propeller, and the outer
as the wall function of STAR-CCM+. The proposed model is part, which is fixed.
introduced to the CFD simulation setup via the coefficients
B = 0, C = 0.2667, R+smooth = 2.5, and R+rough = 25 (Eq. (7)),
to replace the STAR-CCM+ coefficients B = 0, C = 0.253,
R+smooth = 2.25, and R+rough = 90. Hence, the mathematical
formulation and flow calculations around the roughed
propeller are the same as for the smooth propeller except
for the roughness function model.
PROPELLER GEOMETRY
GRID GENERATION
This section presents the CFD results for the propeller. The
computational results are first compared with open-water
experimental results for a smooth propeller, for validation
purposes. The simulation procedure is then applied to
a range of fouling conditions to investigate the effects of
fouling roughness on the open-water performance. Finally,
the required extra power is introduced and estimated in
terms of the specific sea margin, SSM, for the power loss of
the roughed propeller. Values of the SSM are obtained for
roughed and painted propellers.
0,05
0,06 0,07
0,08
0,04
Error
0,05
0,06
0,07
0,04 0,03
ErrorError
0,05
0,06
0,03 0,02
0,02 0,04
0,05
0,01
0,01 0,03
0,04
0
0
0,03 0
0,02 0,5 1 1,5 2 2,5 3 3,5 4
0 0,5 1 1,5 2 2,5 3 3,5 4 4,5
0,01
0,02 Cell no. (mil)
Cell no. (mil)
Fig. 5. Propeller surface grids 0
0,01
Fig. 7.0 Grid convergence
0,5 1 results1,5
for KT 2 2,5 3 3,5 4
0
Fig. 7. Grid
0 convergence
Fig.
0,5 7. 1Grid 1,5 forCell
results KT2 no. (mil)
convergence results 3for
2,5 K3,5
T 4
𝑎𝑎 1,2 (15)
details
where
The GCI s,ofpthe GCI
a, q(p
is also Φ
are
ext presented
21
a),calculated
21
, ea 𝐺𝐺𝐺𝐺𝐺𝐺 ,for 21
extin
21
K=are
e𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓 Table
T and
4.
𝑝𝑝intermediate parameters. (15)
𝑟𝑟21 −K1Q. Values for the numerical uncertainty of 2.81% and
1
0,8
2.40% s,arepis
a,obtained Φextfor
21 the21discretisation 10
where
The GCI q(p a),calculated
also , ea ,for eext21KTareand Kerrors
intermediate of K
Q. Values
andtheKQnumerical
parameters.
Tfor , respectively.
K 0,6
The required
uncertainty data and
of 2.81%
T and
details of the
where
2.40% are obtained GCI
s, p , are
q(p
forpresented
), Φ 21
, ea 21
, in
e Table
21
are 4.
intermediate parameters. 0,4
a the ext discretisation errors 10
of K and K , respectively. The required data and
The GCI is also acalculated for KextT and KQ. Values Tfor the Qnumerical
O
0,2 uncertainty of 2.81% and
details of the GCI are presented in Table 4. 0
2.40% are obtainedThe GCI is also
for thecalculated
discretisationfor KT and KQ. of
errors Values for the
KT and KQ, respectively.
0,6 The
0,7 required
0,8 data
0,9 and 1 1,1 1,2
numerical uncertainty
details of the GCI are presented in Table 4. of 2.81% and 2.40% 10are obtained J
for the discretisation errors of KT and KQ, respectively. The KT (CFD) 10KQ (CFD) Eta (CFD)
10 in Table 4.
required data and details of the GCI are presented
KT (Exp.) 10KQ (Exp.) Eta (Exp.)
Fig. 8. Experimental [32] and CFD simulation results (current study) for a VP1304 propeller under
Tab. 4. Uncertainty calculations for KT and 10
KQ Fig. 8. Experimental
smooth open-water conditions [32] and CFD simulation results (current study) for
a VP1304 propeller under smooth open-water conditions
KT KQ The velocity contour at the surface of the propeller is shown in Fig. 9.
e 21
a
0.50% 0.22%
EFFECT OF FOULING ON THE OPEN-WATER PERFORMANCE OF THE
e 21
ext
2.20% 1.88% PROPELLER
12
0,7
When ηrange of asJ. aItfunction
O is plotted can be of observed
ks, it shows the sameKtendency
that T
rapidly as Kdecreases
T, as discussed k s (Fig.
asabove
11(c)). increases up to k = 1000 μm; for k in the range 1000–3000 μm,
0,6 s s
a marginal
Fig. 11(b) shows KQ as adecreasefunction of kof s. It K
can be is seen
observed, whilewith
that KQ increases KTkstends to of ks
up to a value
0,5 T
= 1000 μm. For ks in the range 1000–3000 μm, a marginal increase in KQ is observed, and KQ is
0,4 a constant for
constant for ks larger than 3000
k greater
s μm.
than 3000 μm.
KT
0,3 When ηO is plotted as a function of ks, it shows the same
The slopes for KT, KQ, and O versus ks are larger at a value of J = 1.2 than J = 0.6. Fig. 11(c) shows
0,2 tendency as K , as discussed above (Fig. 11(c)).
that for a higher value of J,T there is a more significant reduction in O. Figs. 10(c) and 11(c) also
0,1
show that a Fig.
slight 11(b)
increaseshowsin roughnessKQ as a function
leads of k s.inItthecan
to a large decrease valuebe of seen that
O for the propeller.
0,0 This is K increases
conclusive
Q
evidence with thatk up
the
s
to
initiala value
roughness of k
up s
=
to 1000
small μm. For
calcareous k in
fouling
s
the
has crucial
range 1000–3000 μm, a marginalby Songincrease
et al. [25], in K is observed,
0,6 0,7 0,8 0,9 1 1,1 1,2
effect. This finding supports those of studies Farkas et al. [28], and Owen et
It can be seen that KT and o gradually decrease Jwith increasing fouling roughness, while KQ Q
al. [23], among others. Therefore, the importance of the initial roughness on the propeller
increases. For a value Smooth
of J propeller
= 1.2 for medium ks = 100 𝐾𝐾 is reduced
ks = 30 calcareous fouling,
𝑇𝑇
ks = by
300about 26% and and isKemphasised,
performance Q
is constant and
for k s larger
propeller painting
than
is
3000 μm.
recommended as a solution. In other words,
KQ increases by about ks =1000
7.0% with respect ks = 3000 ks = 10000
to the smooth propeller. Consequently, a 30.9% decrease The slopes for K K and η versus k are larger at a value of
the rate of required power will increase T, Q, as the roughness
O increases.
s The propeller roughness arises
in o is observed. The effects of medium and heavy calcareous fouling on the open-water from theJ accumulation
performance of the considered propeller are almost (a)
(a) the same. The reason for this may be related
= 1.2 than J
of= 0.6.
fouling Fig. 11(c)
as a function shows that
of the time for a
operatinghigher value ofthatJ, a large
at sea, meaning
10KQ 0,8
0,6 13
0,4
0,2
0,0
0,6 0,7 0,8 0,9 1 1,1 1,2
J
Smooth propeller ks = 30 ks = 100 ks = 300
ks = 1000 ks = 3000 ks = 10000
(b)
Propeller surface/fouling
condition
J KT ΔKT (%) 10KQ ΔKQ (%) ηO ΔηO (%)
ΔK (%)
1.2 0.2700
ΔK (%)
−10.20 0.7884 0.37 0.6541 −10.53
Propeller surface/fouling J K T 10K T Q Q O ΔO(%)
condition 0.6 0.5576 −6.47 1.3627 0.25 0.3907 −6.70
0.6 0.5961 0 1.3593 0 0.4188 0
0.8 0.5056 0 0.8 1.1744 0.4637
0 −8.28 0
0.5481 1.1789 0.38 0.5008 −8.63
Smooth propeller 1.0µm 0.4008
k s = 300 0 0.9993 0 0.6383 0
1.2 0.3007 0
1.0 0.7855
0.3523
0
−12.10 0
0.7311
0.9984 −0.10 0.5616 −12.01
0.6 0.5889 −1.211.2 1.3630 0.2412
0.27 −19.78−1.48
0.4126 0.7742 −1.44 0.5951 −18.60
0.8 0.4969 −1.71 1.1764 0.17 0.5378 −1.88
ks = 30 µm 1.0 0.3901 −2.660.6 0.9992 0.5317
−0.01 −10.81−2.65
0.6214 1.3725 0.97 0.3699 −11.66
−4.62 −0.42
1.2 0.2868
0.8 0.7822
0.4419 −12.60−4.21
0.7003
1.2042 2.54 0.4672 −14.77
0.6 µm 0.5763
k s = 1,000 −3.33 1.3699 0.78 0.4017 −4.07
0.8 0.4832 −4.421.0 1.1825 0.3553
0.69 −16.33−5.07
0.5203 1.0410 4.17 0.5127 −19.68
ks = 100 µm 1.0 0.3757 −6.27 1.0063 0.70 0.5942 −6.92
1.2 0.2700 −10.201.2 0.7884 0.2292
0.37 −23.77−10.53
0.6541 0.8380 6.67 0.5224 −28.54
0.6 0.5576 −6.470.6 1.3627 0.25
0.5279 −11.44−6.70
0.3907 1.3764 1.26 0.3663 −12.54
0.8 0.4637 −8.28 1.1789 0.38 0.5008 −8.63
ks = 300 µm 1.0 0.3523 0.8
−12.10 0.9984 0.4370
−0.10 −13.57−12.01
0.5616 1.2067 2.75 0.4611 −15.88
1.2 µm 0.2412
k s = 3,000 −19.78 0.7742 −1.44 0.5951 −18.60
0.6 0.5317 −10.811.0 1.3725 0.3295
0.97 −17.80−11.66
0.3699 1.0432 4.39 0.5026 −21.26
0.8 0.4419 −12.601.2 1.2042 0.2225
2.54 −26.02−14.77
0.4672 0.8408 7.04 0.5053 −30.89
ks = 1,000 µm 1.0 0.3553 −16.33 1.0410 4.17 0.5127 −19.68
1.2 0.2292 0.6
−23.77 0.8380 0.5279
6.67 −11.44−28.54
0.5224 1.3764 1.26 0.3663 −12.54
0.6 0.5279 −11.44 1.3764 1.26 0.3663 −12.54
ks = 3,000 µm 0.8 0.4370 0.8
−13.57 1.2067 0.4370
2.75 −13.57−15.88
0.4611 1.2067 2.75 0.4611 −15.88
1.0 µm0.3295
k s = 10,000 −17.801.0 1.0432 0.3295
4.39 −17.80−21.26
0.5026 1.0432 4.39 0.5026 −21.26
1.2 0.2225 −26.02 0.8408 7.04 0.5053 −30.89
0.6 0.5279 −11.441.2 1.3764 0.2225
1.26 −26.02−12.54
0.3663 0.8408 7.04 0.5053 −30.89
ks = 10,000 µm 0.8 0.4370 −13.57 1.2067 2.75 0.4611 −15.88
1.0 0.3295 −17.80 1.0432 4.39 0.5026 −21.26
1.2 0.2225 −26.02 0.8408 7.04 0.5053 −30.89
0,7 1,6
0,6 1,4
1,2
0,5
1,0
0,4
KT 10KQ 0,8
0,3
0,6
0,2
0,4
0,1 0,2
0,0 0,0
0 500 1000 1500 2000 2500 3000 0 500 1000 1500 2000 2500 3000
ks (µm) ks (µm)
(a) (b)
(a) (b)
0,8
0,7
15
0,6
0,5
O 0,4
0,3
0,2
ks (µm)
(b)
0,8
Roughness height for painted propellers
0,7
The roughness of a painted propeller over a period of
0,6
operation can be considered as the summation of the roughness
0,5
of the painted surface (k s1) and the accumulation of fouling
over a certain service time (k s2).
O 0,4
0,3
A value of k s1 = 0.17Ra can be used for painted propellers,
0,2
according to Schultz [33] (Table 8). It should be noted that
0,1
in this study, we use the roughness height of Schultz [33] but
0,0
not the corresponding roughness function model.
0 500 1000 1500 2000 2500 3000 Tab. 7. Annual roughness increments [34]
ks (µm)
Annual roughness increment
Coating type
J = 0.6 J = 0.8 J = 1.0 J = 1.2 Rt50 (µm/year)
(c) Traditional AF coating 40–60
(c) of (a) thrust coefficient, (b) torque coefficient,
Fig. 11. Results for propeller performance: graphs
Self-polishing coating, SPC 10–30
Fig. 11. Results for
andpropeller performance:
(c) efficiency graphs of height
versus roughness (a) thrust coefficient,
(b) torque coefficient, and (c) efficiency versus roughness height Foul-release paints, FR 5–15
NOMENCLATURE
A Roughness constant Uτ Friction velocity
B Smooth wall log-law intercept U +
Non-dimensional velocity
C, Cs Roughness constant ΔU+ Roughness function
D Diameter of propeller Vs Speed of the ship
GCI Grid convergence index VA Propeller advance speed
J Advance coefficient w Wake parameter
ks Equivalent sand-grain roughness height y Normal distance from the wall
k s+ Roughness Reynolds number based on ks y+ Non- dimensional distance from wall
KT Thrust coefficient ρ Fluid density
KQ Torque coefficient ϕκ KT and KQ on the kth grid
n Rotational speed of propeller µ Dynamic viscosity
PB, PBf Engine power (smooth, fouled) κ von Karman constant
PD Delivered power τij Mean viscous stress tensor components
PE Effective power ηD Quasi-propulsive efficiency coefficient
Q, Q f Propeller torque (smooth, fouled) ηH Hull efficiency
Ra, Rt50 Roughness height parameters ηm Mechanical efficiency
RT Total hull resistance ηO, ηOf Open-water efficiency (smooth, fouled)
SSM Specific sea margin ηR Relative rotative efficiency
t Thrust deduction factor ν Kinematic viscosity
T Propeller thrust ρui 'uj ' Reynolds stresses
ui Averaged Cartesian component velocity Uτ Friction velocity
2. A. Banawan, M. Mosleh, and I. Seddiek, “Prediction of the 6. A. Nadery and H. Ghassemi, “Numerical investigation of
fuel saving and emissions reduction by decreasing speed of the hydrodynamic performance of the propeller behind the
a catamaran,” J. Mar. Eng. Techno., (3), vol. 12, pp. 40-48, ship with and without WED,” Pol. Marit. Res. , 4 (108), vol.
2013. 27, pp. 50-59, doi: 10.2478/pomr-2020-0065, 2020.
3. M. H. Ghaemi and H. Zeraatgar, “Impact of propeller 7. M. Perić, “Prediction of cavitation on ships,” Brodogradnja,
emergence on hull, propeller, engine, and fuel consumption vol. 73, no. 3, pp. 39-58, 2022.
performance in regular waves,” Pol. Marit. Res. , (116), vol.
29, pp. 56-76, doi: 10.2478/pomr-2022-0044, 2022. 8. W. Tarełko, “The effect of hull biofouling on parameters
characterizing ship propulsion system efficiency,” Pol.
4. M. Burak Samsul, “Blade cup method for cavitation Marit. Res. , 4(84), vol. 21, pp. 27-34X, 2014, doi: 10.2478/
reduction in marine propellers,” Pol. Marit. Res. , 2 (110), pomr-2014-0038.
vol. 28, pp. 54-62, 10.2478/pomr-2021-0021, 2021.
12. R. Townsin, and S. Dey, “The correlation of roughness 24. S. Song, Y. Demirel, and M. Atlar, “An investigation into
drag with surface characteristics,” Fluid Mechanics the effect of biofouling on full-scale propeller performance
and Its Applications, vol 6. Springer, Dordrecht, doi. using CFD,” in: Proceedings of the ASME 2019 38th
org/10.1007/978-94-011-3526-9_10, 1991. International Conference on Ocean, Offshore & Arctic
Engineering OMAE2019, June 9-14, 2019, Glasgow,
13. Y. K. Demirel, “A CFD model for the frictional Scotland, UK, 2019.
resistance prediction of antifouling coatings,” Ocean
Eng., vol. 89, pp. 21-31, 2014. 25. S. Song, Y. Demirel, and M. Atlar, “Propeller performance
penalty of biofouling: CFD prediction,” J. Offshore Mech.
14. A. Vargas, and H. Shan, “Modeling of ship resistance Arct. , vol. 142, no. 6, p. 0601901, 2020.
as a function of biofouling type, coverage, and spatial
variation,” in 2nd Hull Performance & Insight 26. A. Farkas, N. Degiuli, and I. Martić, “The impact of
Conference, pp. 264-281, Ulrichshusen, Germany, biofouling on the propeller performance,” Ocean Eng.,
2017. vol. 219, p. 108376, 2021.
15. C. Grigson, “Drag losses of new ships caused by hull finis,” 27. A. Farkas, N. Degiuli, and I. Martić, “Assessment of the
J. Ship Res., vol. 36, no. 02, pp. 182-196, 1992. effect of biofilm on the ship hydrodynamic performance
by performance prediction method,” International Journal
16. S. Song et al., “Investigating the effect of heterogeneous of Naval Architecture and Ocean Eng., vol. 13, pp. 102-
hull roughness on ship resistance using CFD,” JMSE, vol. 114, 2021.
9, no. 2, p. 1, 2021.
28. A. Farkas, S. Song, N. Degiuli, I. Martić, and Y.K. Demirel,
17. T. Cebeci and P. Bradshaw, “Momentum transfer “Impact of biofilm on the ship propulsion characteristics
in boundary layers,” in Series in Thermal and and the speed reduction,” Ocean Eng., vol. 199, p. 107033,
Fluids Engineering, P. Bradshaw, Ed. Washington: 2020.
Hemisphere Pub. Corp., 1977.
29. M. P. Schultz and K. A. Flack, “The rough-wall turbulent
18. A. Farkas, N. Degiuli, and I. Martić, “Towards the boundary layer from the hydraulically smooth to the fully
prediction of the effect of biofilm on the ship resistance rough regime,” J. Fluid Mech., vol. 580, pp. 381-405, 2007.
using CFD,” Ocean Eng., vol. 167, pp. 169-186, 2018.
30. SVA Hydrodynamic Solutions, “Potsdam Propeller Test
19. A. Farkas, N. Degiuli, L. Martic, and R. Dejhalla, “Impact Case (PPTC) Open Water Tests with the Model Propeller
of hard fouling on the ship performance of different ship VP1304,” 2011.
forms,” JMSE, vol. 8, p. 748, 2020.
31. L. F. Richardson, and J. A. Gaunt, “The deferred approach to
20. A. Farkas, N. Degiuli, and I. Martić, “Impact of biofilm the limit,” Philosophical Transactions of the Royal Society
on the resistance characteristics and nominal wake,” of London: Series A, vol. 226, pp. 299-361, 1927.
Proceedings of the Institution of Mechanical Engineers,
Part M: Journal of Engineering for the Maritime 32. U. Barkmann, H. Heinke, and L. Lubke, “Potsdam
Environment, vol. 234, no. 1, pp. 59-75, 2020. propeller test case (PPTC) test case description,” in: Second
International Symposium on Marine Propulsors SMP’11,
21. A. Farkas, N. Degiuli, and I. Martić, “A novel method for Hamburg, Germany, Workshop: Propeller Performance,
the determination of frictional resistance coefficient for 2011.
C, the roughness constant, is dependent on the type of roughness of the surface under
experiments identified three flow regimes: a hydraulically smooth regime, a transitionally
23 on the type
eration. where C, the roughness constant, is dependent 23
gh regime, and a fully rough regime.
of roughness of the surface under consideration.
xperimentsThe identified
experiments three flow identified regimes: threea hydraulicallyf low regimes: smooth regime, a transitionally
regime, aand hydraulically
a fully rough smoothregime. regime, a transitionally 23 rough regime,
and a fully rough regime.
In the rough regime, for a roughness different from the 24
sand-grain, the roughness constant is23different. Hence,
Eq. (A.5) is applicable to any rough surface, while C is
individually determined.
Volodimir Leshchev
Igor Maslov
Oleksandr Palagin
Andrii Naydyonov
Department of Ship Power Plants and Systems, Danube Institute of the National University „Odessa
Maritime Academy”, Ukraine
Abstract
The relevance of this study lies in the fact that it presents a mathematical model of the dynamics of the propulsion
system of a ship that takes into consideration the mass of water added to it. The influence of this phenomenon on the
resonant frequencies of the propeller shaft is examined, and a transfer function for a controllable-pitch propeller is
obtained for various operating modes. The purpose of the study is to improve the calculation of the dynamic operating
modes of a controllable-pitch propeller by examining the features of a visual models. The VisSim software package is
used in the study. A visual model is developed that considers the influence of the rotational speed on the value of the
rotational inertia attached to the variable-pitch screw of the mass of water, and a special transfer function is proposed.
The study shows that a transfer function of this type has a loop enabling negative feedback. An analysis of the operation
of the propeller shaft at its resonant frequency is conducted based on the application of frequency characteristics using
the transfer functions obtained. We show that in the low-frequency region, a consideration of the added rotational
inertia using the proposed transfer function leads to a significant difference compared to the result obtained with the
existing calculation method.
Keywords: propulsion system, torsional vibrations, shaft line, added water mass, transfer function, propeller screw
MATERIALS AND METHODS All values of frequencies and angular velocities in the text
and figures are given in rad∙s-1 due to the limitations of the
This study involves mathematical modelling of the process VisSim software product.
that occurs during the movement of the ship, that is, the Table 1. Designations and technical parameters of the ship’s propulsion system
interaction between a controllable-pitch propeller and the
aquatic environment for a propulsion system represented
Elements of the
in the form of a three-mass mechanical system (TMS) [21]. propulsion system Symbol
Parameter values Parameter values
In this approach, three distinct masses represent integral before casting after casting
of the ship
components of the propulsion pathway. The first mass Diesel rotational
Jd 4.50·103 kg. m2 4.50kg. m2
corresponds to the diesel engine, and includes its rotational inertia
inertia and the torque it produces. The intermediary mass Rotational inertia
Jr 2.40·103 kg. m2 0.24kg. m2
represents the gearbox or reduction gear, and accounts for of the gearbox
its rotational inertia and its role in moderating the engine Estimated
speed for the propeller. The final mass represents the propeller rotational inertia Jsp 3.2·103 kg. m2 0.32kg. m2
and its shaft, and includes the propeller’s rotational inertia of the screw
and the resistance it faces due to water and other factors. Rotational inertia Jv=(0.2-
(0.6-1.2)·103 kg. m2 (0.6-0.12)kg. m2
The VisSim software package is used here to visualise the of the added water 0.4)Jsp
results of the study. VisSim is a visual-based block diagram Stiffness
programming language tailored for simulation and embedded of the intermediate Cis 0.64·106 N.m 0.64N.m
system development. It has been extensively employed in the shaft
design of control systems and digital signal processing, and Stiffness
facilitates multi-domain design and simulation. The language of the propeller C ps 0.24×106 N.m 0.024N.m
shaft
offers blocks for arithmetic, Boolean, and transcendental
Gear ratio
operations, in addition to digital filters, transfer functions, ig 3.19 3.19
of the gearbox
numerical integration, and interactive graph plotting [22].
Nominal angular
Figure 1 shows a model diagram of the TMS of a ship’s
rotation speed ωdn 40.317rad∙s-1 1
propulsion system [17, 22], with a modified transmission of the diesel engine
function for the CPP in accordance with the one proposed
in this study. Since the action of the dissipative forces of the viscous
action is assumed to be equal to zero, the values of the resonant
frequencies in the shaft line obtained from experiment are
close to the calculated values [22]. The study was conducted
using the VisSim experimental setup shown in Figure 2.
In addition to the TMS, the model shown in Figure 2 also
contains an exciting scanning signal generator (SSG) [17].
(9) (17)
(11)
where Mps is the torque of the propeller shaft with the FPP.
Part of the second term in Eq. (13) can be interpreted as It follows from Eq. (18) that the model of the total transfer
the rotational inertia of the added water for one revolution function of the propeller has the form shown in
of the propeller, and is expressed by the formula: Figure 4.
(14)
Then, the equation for the propeller dynamics takes the The transfer function of the FPP found here differs
form: significantly from the one used today in practice, which is
accepted as:
(21)
Figure 5. Critical frequencies of the torque Ωcr1-1, Ωcr1-2 , and Ωcr1-3 for the
and in the second, the FPP has a rotational inertia that is TMS propulsion system for various values of the moment Jv
increased by 40%:
(23) CONCLUSION
where is the moment reduction coefficient. Our results show that for the propulsion system considered
here, the critical frequency Ωcr1-2 obtained using the propeller
By substituting the values of the circuit elements (Table 1) transfer function proposed in this study is within the existing
into Eq. (23), the following estimate can be obtained: range of changes in critical frequencies, which are typically
estimated in practice by increasing the moment Jsp by 20–40%.
(24) It corresponds to a total moment Jsp = 0.3885 kg ∙ m2, which
represents an increase of 21.9% and confirms the validity of
its application. The proposed method for determining the
The transfer function of the added moment Wv then takes rotational inertia of a screw with a displaced water mass
the value: allows us to calculate the ratio directly, without the need for
approximate empirical formulae such as Kutuzov’s formula,
NOMENCLATURE 10. Klendii, M., Logusch, I., Dragan, A., Tsvartazkii, I.,
Grabar, A. 2022. Justification and calculation of design
and strength parameters of screw loaders. Machinery
CPP Controllable-pitch propeller & Energetics, 13(4), 48-59. https://doi.org/10.31548/
machenergy.13(4).2022.48-59
FPP Fixed-pitch propeller
11. Yoon, M. 2016. A Transfer Function Model of Thrust
SSG Scanning signal generator Dynamics for Multi-Rotor Helicopters. International
Journal of Engineering Research & Technology, 5(1), 15-18.
TMS Three-mass mechanical system
12. Boletis, E., de Lange, R., Bulten, N. 2015. Impact of
propulsion system integration and controls on the vessel
REFERENCES DP and maneuvering capability. IFAC-PapersOnLine,
48(16), 160-165.
1. MacPherson, D.M., Puleo, V.R., Packard, M.B. 2007.
Estimation of entrained water added mass properties for 13. Xiros, N.I. 2004. PID marine engine speed regulation
vibration analysis. https://cutt.ly/8GVZdOe. under full load conditions for sensitivity H∞-norm
specifications against propeller disturbance. Journal of
2. Faltinsen, O., Minsaas, K.J., Liapias, N., Skjørdal, S.O. Marine Engineering & Technology, 3(2), 3-11.
1981. Prediction of resistance and propulsion of a ship in
a seaway. In: T. Inui (Ed.). Proceedings of 13th Symposium 14. Smailova, G., Yussupova, S., Uderbaeva, A., Kurmangaliyeva,
on Naval Hydrodynamics (pp. 1-19). Tokio: University of L., Balbayev, G., Zhauyt, A. 2018. Calculation and
Trondheim. construction of the tolling roller table. Vibroengineering
Procedia, 18, 14-19. https://doi.org/10.21595/vp.2018.19908
3. Król, P. 2021. Hydrodynamic state of art review: Rotor
– stator marine propulsor systems design. Polish 15. Koushan, K. 2006. Dynamics of ventilated propeller blade
Maritime Research, 28(1), 72-82. https://doi.org/10.2478/ loading on thrusters. In: World Maritime Technology
pomr-2021-0007 Conference (pp. 18-21). London: Macmillan Education.
4. Koibakov, S.M., Umirkhanov, M.G. 2013. Model research 16. Senjanović, I., Hadžić, N., Murawski, L., Vladimir, N. 2019.
of ice jams. World Applied Sciences Journal, 25(8), 1158- Analytical procedures for torsional vibration analysis of
1160. https://doi.org/10.5829/idosi.wasj.2013.25.08.13382 ship power transmission system. Engineering Structures,
178, 227-244.
5. Koibakov, S.M., Umirkhanov, M.G. 2013. Icebreaker unit.
World Applied Sciences Journal, 25(8), 1251-1254. https:// 17. Leschev, V.A. 2018. Marine diesel ACS with external
doi.org/10.5829/idosi.wasj.2013.25.08.13423 feedback speed sensor. Modern Engineering and Innovative
Technologies, 5, 11-17.
6. Gayen, D., Chakraborty, D., Tiwari, R. 2017. Whirl
frequencies and critical speeds of a rotor-bearing system
Dmytro Minchev 1
Roman Varbanets 1
Oleksandr Shumylo 1
Vitalii Zalozh 2
Nadiia Aleksandrovska 1
Pavlo Bratchenko 1
Abstract
The application of Digital Twins is a promising solution for enhancing the efficiency of marine power plant operation,
particularly their important components – marine internal combustion engines (ICE). This work presents the concept
of applying a Performance Digital Twin for monitoring the technical condition and diagnosing malfunctions of marine
ICE, along with its implementation on an experimental test-bench, based on a marine diesel-generator. The main
principles of implementing this concept involve data transmission technologies, from the sensors installed on the engine
to a server. The Digital Twin, also operating on the server, is used to automatically process the acquired experimental
data, accumulate statistics, determine the current technical state of the engine, identify possible malfunctions, and make
decisions regarding changes in operating programs. The core element of the Digital Twin is a mathematical model of the
marine diesel engine’s operating cycle. In its development, significant attention was devoted to refining the fuel combustion
model, as the combustion processes significantly impact both the engine’s fuel efficiency and the level of toxic emissions of
exhaust gases. The enhanced model differs from the base model, by considering the variable value of the average droplets’
diameter during fuel injection. This influence on fuel vapourisation, combustion, and the formation of toxic components
is substantial, as shown. Using the example of calibrating the model to the test results of a diesel engine under 27 operating
modes, it is demonstrated that the application of the improved combustion model allows better adjustment of the Digital
Twin to experimental data, thus achieving a more accurate correspondence to a real engine.
INTRODUCTION as the global market for Digital Twin products was valued
at USD $3.8 billion in 2019 and is expected to reach USD
Digital twin technology is a promising solution for current $35.8 billion by 2025 [2].
and future marine power plants. It could be used in various The concept of Digital Twins has been experiencing rapid
ways for various purposes: advanced control, monitoring and development in recent years. Today, Digital Twins are classified
diagnostics, management and data analysis. as the following types: Digital Twin Instance (DTI), Digital Twin
The idea of Digital Twin technology was originally Prototype (DTP) and Performance Digital Twin (PDT) [3].
developed for product lifecycle management issues, by M. However, all types of Digital Twins have common characteristics:
Grieves in 2002, and subsequently refined by J. Vickers of high-fidelity, dynamic, self-evolving, identifiable, multiscale,
NASA in 2010 [1]. It proved to be a very promising approach, metaphysical, and hierarchical [4], [5].
– – 0.5
lb = din j.holesu0We0.21M 0.16 / (3√ 2ρ )τ ;
––0.6 ;
lc = 8.85din j.holesWe0.25M 0.4 ρ
(1)
d32 = 106Ecdin j.holes (ρ–M
0.0733
n = 50…250 rpm, two-stroke 5…12 5…10 10…15 0.1…0.2 0.6…0.8 1.5…3 0.30…0.65
n = 400…750 rpm, four-stroke 8…15 4…9 8…13 0.05…0.15 0.5…0.7 1…1.1 0.45…0.7
n = 750…1500 rpm, two-stroke 10…40 3…7 4…8 0.05…0.1 0.5…0.7 1…1.2 0.5…0.75
n > 1500 rpm, two-stroke 15…30 3…6 3…7 0.04…0.08 0.6…0.8 1.2…2 0.5…0.9
√π
where Φz is the relative combustion duration; and ζa.c0 and
Φz0 are the coordinates of the minimum of the function
ζa.c = ζa.c(Φz0).
Table 1 provides typical values of the coefficients in the
fuel combustion model equations for different types of
engines, while Fig. 5 illustrates the influence of various tuning
parameters on the characteristics of heat release.
Fig. 7. Calculated diagrams of fuel evaporation rate. Medium-speed
diesel engine, 40% of MCR.
Considering variable fuel droplets
size during injection
K3p[N][OH]– K3r[NO][H],
where the square brackets ‘[ ]’ express the volumetric and after fuel injection ends:
concentration of the corresponding matter, and K1p, K1r, K2p,
K2r, K3p, and K3r are constants for direct and reverse chemical
reactions.
The NO formation kinetic equation incorporates the first
and second equations, based on Zvonov’s approach:
In the above equations B1s, B’2s, B’’2s, B3s, and B4s are empiric
coefficients; δd is the droplet core size; Kev is the evaporation
constant; τinj.start and τinj.end are moments of time for injection
start and injection end; xinj.end is the burned fuel fraction for
(15) the moment of injection end; ndisp is the distibution constant
where p is the pressure in the cylinder (bar); R = 8.3144 J/ to consider fuel injection uniformity; k02is the oxidation
(mole·K) is the gas constant; and Tburned is the temperature of coefficient; and [C] is the volumetric soot concentration.
the burned gases.
The Arrhenius Law equations are used to calculate reaction
rate constants:
K = ATB exp (– RT
Ea
), (16)
where A and B are empirical coefficients; and Ea is the
activation energy.
For high-speed engines, the final concentration of CO in
exhaust gases is estimated as an equivalent concentration at
the combustion finish point. For medium-speed and low-
speed diesel engines, the following kinetic equation is used:
d[CO]
dτ
= K1C[CO] [OH] , (17)
where K1C = 7.1 ∙ 1012∙ e-32200RT is the reaction constant; and [CO]
and [OH] are the corresponding CO and OH concentrations.
Calculation of the formation of particulate matter in diesel
engines is performed using Razlejtsev’s approach. The equation
for the instantaneous volumetric soot concentration rate is
as follows:
d[C]
dτ
(
= d[C]
dτ
)kin+(d[C] ) –(d[C]
dτ pol
) –(d[C]
dτ burn
) .
dτ vol
(18) Fig. 9. Soot concentration and nitrogen oxide development in the cylinder.
Medium-speed engine, MCR.
Gang Wu 1
*
Guodong Jiang 1
Changsheng Chen 2
Guohe Jiang 1
Xigang Pu 1
Biwen Chen 2
1
Merchant Marine College, Shanghai Maritime University, Shanghai, China
2
Shanghai Marine Equipment Research Institute, Shanghai, China
Abstract
Two-stroke, low-speed diesel engines are widely used in large ships due to their good performance and fuel economy.
However, there have been few studies of the effects of lubricating oils on the vibration of two-stroke, low-speed diesel
engines. In this work, the effects of three different lubricating oils on the vibration characteristics of a low-speed engine are
investigated, using the frequency domain, time-frequency domain, fast Fourier transform (FFT) and short-time Fourier
transform (STFT) methods. The results show that non-invasive condition monitoring of the wear to a cylinder liner in
a low-speed marine engine can be successfully achieved based on vibration signals. Both the FFT and STFT methods
are capable of capturing information about combustion in the cylinder online in real time, and the STFT method also
provides the ability to visualise the results with more comprehensive information. From the online condition monitoring
of vibration signals, cylinder lubricants with medium viscosity and medium alkali content are found to have the best
wear protection properties. This result is consistent with those of an elemental analysis of cylinder lubrication properties
and an analysis of the data measured from a piston lifted from the cylinder after 300 h of engine operation.
Keywords: two-stroke; low-speed marine diesel engine; cylinder lubricating oils; vibration characteristic; condition monitoring
INTRODUCTION can have direct effects on the hull, machinery, crew, and
passengers when it occurs in ships. The main source of engine
The internal combustion (IC) engine is a classical rotating vibration is an imbalance in the force and moment created by
power machine that has been widely employed for transportation the inertial force of its moving parts, although another source
and power production [1]. In particular, the two-stroke low- is the side thrust and overturning moment generated by the
speed diesel engine is a unique type of IC engine that has been in-cylinder pressure of fuel combustion. Further factors include
used in many medium and large ships due to its advantages the influence of lubricating oils [3, 4], fuels [5-8], and engine
such as multi-fuel feasibility, fuel economy, large output power, faults [9, 10], which can all affect engine vibration. The engine
and durability [2]. With the increasing emphasis on ship safety, contains many internal components, such as the piston-cylinder
more and more researchers are focusing on the hazards to assembly, which require a good lubricating oil to run; hence,
reliability posed by engine vibration. Harmful engine vibration the characteristics of the lubricating oil, and its quality and
Item/parameter Details
Model 6S35ME-B9
6-cylinder, water-cooled, direct
Engine type injection, crosshead type
Combustion order 1-5-3-4-2-6
Fig. 1. Schematic diagram of the experimental setup.
Bore & stroke 350 × 1550 mm
The vibration signals were measured with a tri-axial vibration
Compression ratio 17:1
accelerometer (made by Coinv, China). Detailed specifications
Displacement 894 L of the transducers are shown in Table 3. As can be seen from
Rated output power 3750 kW @ 142 rpm Fig. 1, the accelerometer was mounted on the engine cylinder
cover, and the vibrations in three directions (x, y, z) were all
Maximum indicated pressure 180 bar
acquired by the accelerometer simultaneously. The vibration
Valve opening pressure 380 bar signals were measured under the same engine load condition
for each of the three cylinder lubricating oils. The angle of the
Three different formulas of cylinder lubricating oil (recipes A, engine crankshaft was recorded by the proximity sensor during
B, and C) were prepared and used. These three cylinder oils all the vibration measurements. The coordinate system used for
contained the same percentages of additives to the base oil, which both the accelerometer and the engine is shown to the lower
included compound ester, ester base, amide, phosphate ester, right of Fig. 1.
Sensitivity 100 mV/g where xRMS (m/s ) is the value of the RMS for the calculated
2
25% 98 rpm 812.5 kW 91.64 kN·m To compare the total vibration in each cylinder based on
the kurtosis, Kurtosistotal was derived from the values for the
40% 105 rpm 1300 kW 116.13 kN·m
x, y and z axes using Eq. (4) [15]:
50% 113 rpm 1625 kW 136.40 kN·m
KurtosisTotal = √Kurtx2 + Kurty2 + Kurtz2 (4)
ANALYSIS METHODS
where KurtosisTotal is the total vibration acceleration, and Kurtx,
In order to compare the differences between the vibration Kurty and Kurtz are the values of the kurtosis for accelerations
signals in the time domain, the values of the RMS, kurtosis, along the x, y, and z axes, respectively.
mean, skewness, standard deviation, and squared deviation of
the signals were considered. Mean of acceleration signals
For the frequency and time-frequency domain analyses, the The mean is the sum of all the data in a set divided by the
FFT and STFT were used. These analyses were all performed total number of data points. It is used to reflect the general
using MATLAB (2019 b) and Origin software. overall level of data, and can be calculated using Eq. (5) [9]:
N
RMS of acceleration signal x– = N1 Σ xi (5)
i=1
The RMS represents the strength of the signal, and is also an
expression of the average energy of the signal. The RMS along
Standard deviation of acceleration signals The objectives of this study were to compare the effects of
The standard deviation of the arithmetic square root of the three different lubricating oils on a two-stroke low-speed marine
squared deviation, which indicates the degree of dispersion diesel engine, and to explore the vibration of the engine under
between samples of vibration data, is defined by Eq. (7) [10]: low to medium load conditions. The results of the experimental
study are therefore presented in four sections relating to the
ΣNi=1(xi –x–)2
σ= N (7) time domain, the frequency domain, a time-frequency domain
analysis of the vibration signal, and an elemental spectra
Squared deviation of acceleration signals analysis.
The squared deviation of signals represents the degree
of fluctuation of the signal around the mean; it can also be
understood as the degree of vibration signal shock that exists.
It is calculated using Eq. (8):
Σ (xi – x–)2 (8)
N
δ = N1
i=1
where n is the length of signal X, and Wn = e(–2πi)/n. The lubricating oil not only provides lubrication and
reduces friction between the cylinder and piston ring, but
Short-time Fourier transform also participates in the combustion process. The vibration of
In the STFT, a signal is split into multiple segments using an IC engine mainly arises from the gas force generated by
a window function (in our analysis, the Hamming window the combustion in the cylinders, and a prior study [19] was
function was used). The FFT is then applied to calculate each conducted to explore the in-cylinder combustion process of
segment of the signal, and the spectrograms obtained from a low-speed engine using the CFD method. The usage of different
the FFT are arranged along the time axis, to obtain the time- lubricating oils is also a crucial factor affecting the vibration in
frequency diagram. Hence, one of the advantages of the STFT the engine. To investigate the detailed information for the time
over the FFT is that the time and frequency information can domain, we consider the z 3 s segment of the signal in Fig. 3,
both be retained simultaneously, as the time information is lost where a complete engine cycle is marked by red lines. This signal
when the FFT is used. Another advantage is that the STFT can represents the conditions with lubricating oil recipe B, at an
be applied to non-stationary signals, which is not possible for engine speed of 98 rpm under a load of 25%, along the z-axis.
the FFT. The formula for the STFT is given in Eq. (10) [17]: As shown in the figure, the firing order of the cylinders is 1-5-
+∞ 3-4-2-6. The acceleration of the vibration signal in cylinder 4
STIF(τ, f ) = ∫ x(t)h(t – τ)e–j2πft (10)
–∞
has the maximum value, while that of cylinder 1 shows the
minimum. The reason for this is that the accelerometer was
where h is a window function. f and τ denote the frequency mounted on the head of cylinder 4, whereas cylinder 1 was
and time variables, respectively. farthest from the accelerometer. The acceleration values for
However, the STFT also has certain limitations, as the size the other cylinders decrease with increasing distance from
of the window cannot change with frequency. Thus, the time cylinder 4. The recording period for the vibration signal was
and frequency resolution cannot be optimised at the same 30.1056 s, containing 150,528 samples. Hence, for an engine
time in the STFT. The shorter the window, the more obvious speed of 98 rpm, it can be calculated that the duration of
Fig. 5. Plots of feature parameter values for three engine loads for three
recipes of lubricating oil: (a) mean; (b) skewness; (c) standard deviation;
(d) squared deviation.
Fig. 5 shows the values of the mean, skewness, standard
deviation, and squared deviation. For each of the three loads, the
mean values for recipe A are almost all near zero, while the mean
values for recipe B remain near zero only for the 40% and 50%
load conditions. It is not hard to observe from Fig. 5(a) that the
mean values of all three lubricating oils are relatively steady and
Fig. 4. Plots of feature parameter values for three engine loads for three recipes remain near zero under the 40% load condition. Fig. 5(b) shows
of lubricating oil: (a) RMS; (b) ; (c) kurtosis; (d) . that under the 40% and 50% load conditions, recipe B has low
The RMS values for each recipe of lubricating oil are low values for the RMS, kurtosis, and mean, whereas its skewness
under the 40% load condition. The reason for this is related to values are almost all negative. In contrast, the data points for
the characteristics of this diesel engine, as different lubricating recipe A are uniformly distributed on both sides of zero, and
oils have little effect on the intensity of engine vibrations under do not vary greatly between the load conditions. Furthermore,
Fig. 6. FFT spectra for recipe A under three typical engine loads: Fig. 7. FFT spectra for recipe B under three typical engine loads:
(a) 25%, (b) 40%, (c) 50%. (a) 25%, (b) 40%, (c) 50%.
Figs. 6–8 show FFT spectra for vibration signals under In the range 1.3–1.75 kHz, all three recipes for the lubricating
three different engine loads for the three recipes of lubricating oil show different trends with increasing load. However,
oil. It can be observed that the vibration signatures of the under the 50% load condition, their maximum values are all
diesel engine change when different cylinder lubricating at a relatively low level of below 0.34 m/s2. The changes in the
oils are used, from Fig. 6 to Fig. 8. The changed frequency maximum values for recipes A and B are similar in the high-
domain of 0–2.5 kHz can be divided into four subdomains: frequency subdomain (1.9–2.5 kHz), as both increase with
0.12–0.38 kHz, 0.75–1.3 kHz, 1.3–1.75 kHz, and 1.9–2.5 kHz. rising loads. The increase in the maximum value for recipe A is
In the range 0.12–0.38 kHz, which is the low-frequency region 198%, and for recipe B it is 161%. It is notable that the maximum
of engine vibration, the accelerations all decrease to varying value under 25% load is 0.426 m/s2 for recipe C, which is higher
extents for all three lubricating oils as the load increases, and than for its principal frequency domain. However, under the
especially their maximum values. Due to its high viscosity 40% load condition, the maximum value is only 0.1 m/s2,
and alkaline nature, the decrease in the maximum value for while the value is 0.352 m/s2 under the 50% load condition.
recipe A at a higher load is greater than that for recipes B and C, We conclude that under low to medium loads, using lubricating
with a decline of 60.2%. The maximum values for recipes B oil with recipe C causes high-frequency vibration fluctuations
and C are fairly similar, as they have the same SAE 40 grade. in the engine. This also provides some evidence that a TBN
It is known that the low-frequency components of vibration level that is too low may affect the detergency between the
are the main contributors to noise and discomfort for human piston and liner.
Fig. 9. STFT diagrams for lubricating oils under 40% load: (a) recipe A; (b) recipe B; (c) recipe C.
Author Contributions
8. G. Chiatti, O. Chiavola, and F. Palmieri, „Vibration and 19. X. Zhou, T. Li, N. Wang, X. Wang, R. Chen, and S. Li, „Pilot
acoustic characteristics of a city-car engine fueled with diesel-ignited ammonia dual fuel low-speed marine engines:
biodiesel blends,” Appl. Energ., vol. 185, pp. 664-670, 2017, A comparative analysis of ammonia premixed and high-
https://doi.org/10.1016/j.apenergy.2016.10.119. pressure spray combustion modes with CFD simulation,”
Renew. Sust. Energ. Rev., vol. 173, 2023, https://doi.org/10.1016/j.
rser.2022.113108.
Volodymyr Korobko 1
Serhiy Serbin 1
*
Huu Cuong Le 2
1
Admiral Makarov National University of Shipbuilding, Mykolaiv, Ukraine
2
Institute of Maritime, Ho Chi Minh City University of Transport, Viet Nam
Abstract
The utilisation of the thermal emissions of modern ship power plants requires the development and implementation of
essentially new methods of using low-temperature waste heat. Thermoacoustic technologies are able to effectively use low-
temperature and cryogenic heat resources with a potential difference of 500–111 K. Thermoacoustic heat machines (TAHMs)
are characterised by high reliability, simplicity and environmental safety. The wide implementation of thermoacoustic
energy-saving systems is hampered by the low specific power and the difficulties of directly producing mechanical work.
An efficient approach to converting acoustic energy into mechanical work entails the utilisation of axial pulse bidirectional
turbines within thermoacoustic heat engines. These thermoacoustic turbogenerators represent comprehensive systems that
consist of thermoacoustic primary movers with an electric generator actuated by an axial-pulse bidirectional turbine. The
development of such a thermoacoustic turbogenerator requires several fundamental issues to be solved. For this purpose,
a suitable experimental setup and a 3D computational fluid dynamics (CFD) model of a thermoacoustic engine (TAE)
with bidirectional turbines were created. The research program involved conducting physical experiments and the CFD
modelling of processes in a TAE resonator with an installed bidirectional turbine. The boundary and initial conditions for
CFD calculations were based on empirical data. The adequacy of the developed numerical model was substantiated by the
results of physical experiments. The CFD results showed that the most significant energy losses in bidirectional turbines
are manifested in the output grid of the turbine.
Keywords: waste heat recovery; ship power plant; thermoacoustics; thermoacoustic engine; bidirectional turbine
smaller dimensions. In this regard, it is possible to expect the electric grid. The characteristics of bidirectional turbines in
formation of stable oscillating structures in the TATG resonator, a unidirectional flow were studied using the experimental
which will introduce additional excitations and inhomogeneities installation, the scheme of which is shown in Fig. 3.
at the input to the rectifier of the turbine. They will create a novel Tab. 2. Design parameters of bidirectional turbines
energy-consuming mechanism that should be confined on both Element Value
sides of the bidirectional turbine. In view of these circumstances, Internal diameter of TAE resonator, mm 46
it can be considered appropriate to conduct a complex of studies
Outer diameter of turbine stator, mm 76
to obtain an in-depth understanding of these issues.
Outer blade diameter of turbine, mm 75
Purpose of research. Considering the disparities in the Turbine rotor diameter, mm 50
working processes between OWCs and TAEs, it is sensible to Blade width of rotor, mm 20
carry out research focused on understanding the peculiarities of Number of rotor blades 24
hydrodynamic processes within the TATG resonator to identify Width of stator blades, mm 16
the most rational design solutions. Research should involve
Number of stator blades 18
both physical experiments and numerical calculations using
computational fluid dynamics (CFD) procedures, which requires Differential thermostabilised pressure sensors of MPXV NXP
the creation of appropriate equipment and methods of research. semiconductors were used to measure the pressure drop along
the turbine and on the entrance lemniscate. The equipment of
the experimental stands and the available microprocessor-based
Experimental setup data acquisition system (DAS) provided multiple measurements
of sensor readings in each of the modes, which reduced the
The investigation was conducted utilising an experimental random error.
bidirectional turbine. This turbine was designed and The STM32F407VGT6 microcontroller, which has a 32-bit
manufactured in accordance with the guidelines outlined in architecture and an operating frequency of 168 MHz, was used
[24]. A 3D model of an experimental bidirectional turbine to measure and register signals from pressure sensors, current
was developed in the CAD software developed by SolidWorks. sensors, voltage sensors, etc.
The turbine’s components, including the rotor, guide grids and The equipment of the test stands ensured multiple
fairings, were produced using 3D printing techniques and plastic measurements of sensor readings in each of the modes, which
material, as depicted in Fig. 2. Detailed information about the reduced the random error.
dimensions of the elements of the experimental bidirectional
turbine is provided in Table 2.
The rotor of the turbine was rigidly coupled to a three-phase
brushless electric generator of the Sankyo F2JGL type. The
output voltage from the generator was fed to a diode converter
made of Schottky diodes, and a P517-M laboratory rheostat or
a set of precision resistors was used as a load.
The rotation speed of the turbine was monitored by an SDS
1074CFL digital oscilloscope, which was connected to one of
the windings of the generator.
Fig. 3. Diagram of the experimental stand for the study of the hydraulic resistance
Such a scheme made it possible to simultaneously control of the bidirectional turbine: 1 – lemniscate; 2 – inclined differential
the speed of rotation of the turbine rotor, the root-mean-square pressure gauge; 3 – turbine; 4 – fan; 5 – LATR; 6 – DAS oscilloscope
(RMS) voltage of the generator VRMS and the current in the SDS 1074CFL; 7 – diode bridge; 8 – ammeter; 9 – voltmeter; 10 – resistors
Fig. 9. Distribution of static pressure in the resonator channel with bidirectional turbines at different flow rates
16. Z. Yang, V. Korobko, M. Radchenko, and R. Radchenko, 26. ANSYS, Inc., ANSYS Fluent Theory Guide. ANSYS, Inc., 2013.
“Improving thermoacoustic low-temperature heat recovery
systems,” Sustainability (Switzerland), Vol. 14, No. 19, art. 27. S. I. Serbin, I. B. Matveev, and G. B. Mostipanenko, “Plasma-
No. 12306, 2022, doi: 10.3390/su141912306. assisted reforming of natural gas for GTL: Part II - Modeling
of the methane-oxygen reformer,” IEEE Trans. Plasma
17. T. K. Das, P. Halder, and A. Samad, “Optimal design of air Sci., Vol. 43, No. 12, pp. 3964–3968, 2015, doi: 10.1109/
turbines for oscillating water column wave energy systems: TPS.2015.2438174.
A review,” Int. J. Ocean Clim. Syst., Vol. 8, No. 1, pp. 37–49,
2017, doi: 10.1177/1759313117693639. 28. I. Matveev, S. Serbin, T. Butcher, and N. K. Tutu, “Flow structure
investigations in a Tornado combustor,” in 4th International
18. A. F. O. Falcao and J. C. C. Henriques, “Oscillating-water- Energy Conversion Engineering Conference, AIAA2006-4141,
column wave energy converters and air turbines: A review,” Vol. 2, 2006, pp. 1001–1013, doi: 10.2514/6.2006-4141.
Renewable Energy, 2015, doi: 10.1016/j.renene.2015.07.086.
29. O. Cherednichenko, S. Serbin, and M. Dzida, “Investigation
19. A. Thakker and F. Hourigan, “Modeling and scaling of the of the combustion processes in the gas turbine module of
impulse turbine for wave power applications,” Renewable an FPSO operating on associated gas conversion products,”
Energy, Vol. 29, No. 3, pp. 305–317, 2004, doi: 10.1016/ Polish Maritime Research, Vol. 26, No. 4, pp. 149–156, 2020,
S0960-1481(03)00253-2. doi: 10.2478/pomr-2019-0077.
20. D. Liu, Y. Chen, W. Dai, et al., “Acoustic characteristics of 30. S. Serbin, K. Burunsuz, M. Dzida, J. Kowalski, and D. Chen,
bi-directional turbines for thermoacoustic generators,” “Investigation of ecological parameters of a gas turbine
Front. Energy, Vol. 16, pp. 1027–1036, 2022, doi: 10.1007/ combustion chamber with steam injection for the floating
s11708-020-0702-3. production, storage, and offloading vessel,” International
Journal of Energy and Environmental Engineering, Vol. 13,
21. M. A. Elhawary, A. H. Ibrahim, A. S. Sabry, and E. Abdel- No. 3, pp. 873–888, 2022, doi: 10.1007/s40095-021-00433-w.
Rahman, “Experimental study of a small scale bi-directional
axial impulse turbine for acoustic-to-mechanical power 31. I. B. Matveev, N. V. Washchilenko, and S.I. Serbin, “Plasma-
conversion,” Renewable Energy, 2020, doi: 10.1016/j. Assisted Reforming of Natural Gas for GTL: Part III - Gas
renene.2020.05.162. Turbine Integrated GTL,” IEEE Trans. Plasma Sci., Vol. 43,
No. 12, pp. 3969–3973, 2015, doi: /10.1109/TPS.2015.2464236.
22. C. Iniesta, J. L. Olazagoitia, J. Vinolas, and J. Aranceta,
“Review of travelling-wave thermoacoustic electric-generator
technology,” in Proceedings of the Institution of Mechanical
Engineers, Part A: Journal of Power and Energy, 2018, doi:
10.1177/0957650918760627.
Jerzy Girtler
Jacek Rudnicki *
Abstract
The article provides proof that the diagnostics of marine internal combustion engines and other ship power plant
machines should take into account the randomness and unpredictability of certain events, such as wear, damage, the
variations of mechanical and thermal loads, etc., which take place during machine operation. In the article, the energy
E, like the other forms (methods) that it can be converted into (heat and work), is considered the random variable Et;
at time t, this variable has the mean value , which is the observed value of the statistic with an asymptotically
normal distribution , irrespective of the functional form of the random variable Et. A proof is given
that shows that the expected value estimated in the above way, considering the time t of the performance of task Z
by a marine internal combustion engine or other ship power plant machine, can be used to determine the machine’s
possible action (DM). When compared to the required action (DW) needed for task Z to be performed, this possible action
makes it possible to formulate an operating diagnosis concerning whether the engine or machine of concern is able
to perform task Z. It is assumed that an energy device of this type is able to perform a given task when the inequality
DM≥DW holds. Otherwise, when DM < DW, the device cannot perform the task for which it was adopted in the design
and manufacturing phase, which means that it is in the incapability state, although it still can be started and convert
energy into the form of heat or work..
(1)
(2) (5)
which means that in this case, the value WE of the energy E where , D is the engine (machine) operation;
is equal to the value WL of the work L, i.e. WE = WL . E is the converted (obtained) energy that allows the realisation
Similarly, when the energy E is solely used for generating of the task in the time interval [t1, t2], and t is the time of
the heat Q, then, in the evaluative approach, the following conversion of the energy E.
equivalence occurs: Therefore, if we assume that , then
formula (5) can be written as (Fig. 1)
(3)
(6)
which means that in this case, the value WE of the energy E
is equal to the value WQ of the heat Q, i.e. WE = WQ. The use of formula (5) or formula (6) requires the geometric
The operation of a marine internal propulsion engine or application of a definite integral, and the following inequalities
other ship power plant machine can be considered using must be taken into account during the integration:
the following terms: required operation (DW) and possible
operation (DM) [16]. We can conclude that each engine or
machine is in the capability state, i.e. it can perform a given
task, when The integral given by formula (6) is the Riemann definite
integral [7], with the integration interval defined in this case
(4) as equal to [t1, t2] and the integrand . This
function is integrable in the Riemann sense in the above time
Otherwise, when DM < DW, we can conclude than the interval according to the following formula:
engine or machine is in the incapability state or partial
incapability state [15], [16], [32]. The capability of an engine or (7)
machine can be assessed after comparing the area of required
operation (DW) with that of possible operation (DM). This issue Hence, if we can determine the functional relation between
is discussed in [16]. the energy (E) and time (t) that characterises the operation
In a deterministic approach, the operation of a given of an engine or machine, i.e. the function E = f(t), and this
marine engine or other ship power plant machine can be function is continuous, for instance, in a given time interval
described using a general functional relationship describing [t1, t2], then, according to the second fundamental theorem of
the change in the energy E at an arbitrary time t of machine calculus (Newton-Leibniz theorem), we can write
operation. The operation of the marine engine or other ship
power plant machine analysed for E(t) = f(t) in a given time (8)
interval, e.g. [t1, t2], is shown as the area in Fig. 1.
The application of the Newton-Leibniz theorem is necessary
here because it enables the effective calculation of a definite
integral of any continuous function if an antiderivative of
this function is known.
Fig. 3. Interpretation of a) the accumulation of dissipated energy Er and b) the resulting decrease in the useful energy Eu of an internal combustion engine or other
ship power plant machine: e – portion (quantum) of energy by which the energies Er and Eu change, Erg – dissipated energy limit, Eug – useful energy limit [17]
17. J. Girtler, ‘Aspects of quantum machine diagnosis,’ XXXII 33. J. Rudnicki, ‘Quantitative assessment of operation of
National Symposium ‘Machinery Diagnostics’, Węgierska ship main diesel engine,’IV International Conference on
Górka 2012 (in Polish). Computational Methods in Marine Engineering, Lisbon
2011.
18. F. Grabski, ‘Theory of semi-Markov processes of operation
of technical objects’, Science Notebook of Naval College, 34. B. F. Schutz, An introduction to general relativity theory.
vol. 75A, Gdynia 1982 (in Polish). PWN. Warsaw 2002 (in Polish).
19. J. Gribbin, In Search of Schrödinger’s Cat. Zysk i S-ka. 35. M. Tempczyk, Chaos Theory for the Brave. PWN. Warsaw
Poznań 1997 (in Polish). 2002 (in Polish).
20. S. W. Hawking, A Brief History of Time. From the Big Bang 36. M. White and J. Gribbin, Einstein. A Life in Science. WNT.
to Black Holes. ALFA Publisher, Warsaw 1990 (in Polish). Warsaw 1995 (in Polish).
21. R. Michalski, Diagnostics of Working Machinery. Radom. 37. C. M. Wynn and A. W. Wiggins, The Five Biggest ideals in
Institute for Sustainable Technologies 2004 (in Polish). Science. Prószyński S-ka. Warsaw 1998 (in Polish).
22. M. Moskowits, Is science good. CiS Publisher. Warsaw 38. B. Żółtowski, Fundamentals of Machine Diagnostics.
1997 (in Polish). University of Technology and Agriculture. Bydgoszcz
1996 (in Polish).
23. R.F. Nalewajski, Fundamentals and Methods of Quantum
Chemistry. PWN. Warsaw 2001 (in Polish). 39. B. Żółtowski and C. Cempel, Eds., Engineering of
Machine Diagnostics. The Polish Association of Technical
24. S. Niziński, Elements of Diagnosis of Technical Objects. Diagnostics. Radom Institute for Sustainable Technologies.
University of Warmia and Mazury, Olsztyn 2001 (in Polish). Warsaw, Bydgoszcz 2004 (in Polish).
25. S. Pabis, Methodology and Methods of Empirical Sciences. 40. J. Kiciński, Ed., Modeling and Diagnostics of Mechanical,
PWN. Warsaw 1985 (in Polish). Aerodynamic and Magnetic in Power Turbine Units.
Institute of Fluid-Flow Machinery Polish Academy of
26. Z. Pawłowski, Mathematical Statistics. PWN. Warsaw 1976 Sciences. Gdańsk 2005 (in Polish).
(in Polish).
Piotr Bilon
Wojciech Litwin
Gdansk University of Technology, Institute of Naval Architecture, Gdansk, Poland
Abstract
The digitalisation in the ship design sector is currently applied to the design process itself and is well defined, partially
standardised and practically implemented by both independent design companies and the design departments of
shipyards. The situation is similar in other sectors of engineering. However, the requirements for the digitalisation
of other processes in design and engineering companies have not previously been studied, and the limited financial
resources of ship design companies mean that there is a need for research on the digitalisation needs of multidisciplinary
ship design companies. The implementation of building information modelling (BIM) solutions is costly for design
companies and generates benefits mainly for shipyards and shipowners. The lack of benefits for design companies leads
to the hesitation of managers whenever digitalisation is considered; the scope and scale of the digitalisation, as well
as the selected software and BIM level, are defined by the shipyard requirements. The participation and support of
management in the digitalisation process is one of the key success factors; the expected benefits caused by digitalisation
for the organisation will increase the motivation of managers to invest limited resources in digitalisation. There are no
data that indicate the processes with a high potential for digitalisation and the scale of expected improvements in ship
design companies; therefore, research in this area was performed with a group of project managers from design and
engineering companies. The research focused on collecting the opinions and experiences of the managers related to the
manual management of resources and comparing the poll results with the conclusions from the enterprise resource
planning (ERP) system data analyses. The paper analyses if the digital automation of the resource management process
can lead to the substantial improvement of the operations of multi-project, multidisciplinary engineering ship design
companies.
Keywords: Ship design human resource management, digital automation, project portfolio, multidisciplinary engineering,
ship design, marine engineering
RESULTS
with a dedicated ERP software vendor for the design the implementation
Most of the respondents (1242/1283, orof93%)
the planning procedure
define the elimination of manual resource
management as an important aspect of improving the efficiency of a design company and state that
they develop several variants of future plans as part of manual management, while 75% of all
respondents (957/1283) state their interest in a software solution that supports the preparation of
future variants to support the resource management process in a design and engineering company.
POLISH
The combined MARITIME
results are RESEARCH,
presented in Table 1. No 4/2023 123
Table 1. ERP software user questionnaire results
1. Can you describe your company as a design and engineering multidisciplinary corporation?
No. Answer Res. %
1.1. Yes 1187 93%
1.2. No 96 7%
2. Please indicate your sector of engineering
2.1. Marine (ship design) 430 34%
Most of the respondents (1242/1283, or 93%) define the 7. How is the % of the project management financial budget
elimination of manual resource management as an important consumed on resources management?
aspect of improving the efficiency of a design company and 7.1. <5% 17 1%
state that they develop several variants of future plans as 7.2. <10% 697 54%
part of manual management, while 75% of all respondents 7.3. <15% 565 44%
(957/1283) state their interest in a software solution that 7.4. Other 4 0%
supports the preparation of future variants to support the
8. How big is the impact of resources manual management
resource management process in a design and engineering elimination on the efficiency of engineering company?
company. The combined results are presented in Table 1. 8.1. Critical 671 52%
Tab 1. ERP software user questionnaire results 8.2. Significant 571 45%
1. Can you describe your company as a design and engineering 8.3. Minor 32 2%
multidisciplinary corporation?
8.4. No impact at all 9 1%
No. Answer Res. %
9. Do you prepare several variants of resources management
1.1. Yes 1187 93% projection?
1.2. No 96 7% 9.1. Always 25 2%
2. Please indicate your sector of engineering 9.2. When requested and when I feel 654 51%
2.1. Marine (ship design) 430 34% it is needed
2.4. Plant design 257 20% 10. Would you be willing to use the software tool for automated
resources management even if it creates only variants for
2.5. Airspace 7 1% consideration?
2.6. Energy 113 9% 10.1. Yes 957 75%
3. Please indicate the methodology of resources management in your 10.2. No 258 20%
organization
10.3. I do not know 68 5%
3.1. Manual in spreadsheet 893 70%
3.2. In resources matrix in 387 30%
planning/project management For some companies, manual management is just a part
software of general department management; it might be registered
3.3. In dedicated automated 0 0% as a general management task, just like periodical project
resources management software meetings and even quality control activities. This means that
3.4. Other 3 0% in this first iteration of the research, the assignment of tasks
4. How big is the % of the hourly budget for project management to the appropriate category may be imprecise, due to the fact
consumed on resources management? that the users of the ERP system did not use attributes that
No Answer Res. % clearly define tasks for managers for tasks related to manual
management. This caused the need to manually review each
4.1. <25% 0 0%
stage of the selected projects and each task in order to manually
4.2. <50% 178 14%
assign the tasks and the hours allocated to each task to the
4.3. <75% 673 52% proper group. On the other hand, an unambiguous result is
4.4. Other 432 34% that 18–27% of the financial budget and 22–28% of the hourly
5. How big is the % of the financial budget for project management budget were used for activities related to the management
consumed on resources management? of project work and resources; this is applicable to virtually
5.1. <25% 0 0% all engineering sectors. Considering these numbers in the
5.2. <50% 273 21% context of the average margins on design and engineering
5.3. <75% 437 34% services should focus the attention of employees responsible
for business efficiency on this particular area of the company’s
5.4. Other 573 45%
operations.
6. How big is the % of the project hourly budget consumed on
resources management?
6.1. <5% 767 60%
DISCUSSION
6.2. <10% 389 30%
6.3. <15% 124 10% The subject of project portfolio selection and management
6.4. Other 3 0% is of scientific interest, and it is gaining momentum when it
comes to the number of publications on the topic. Currently,
researchers are exploring the area of decision support systems,
and the trend is to develop more management-oriented
12. M. Rauf, Z. Guan, L. Yue, Z. Guo, J. Mumtaz, and S. 22. W. Litwin, D. Piątek, W. Leśniewski, and K. Marszałkowski,
Ullah, “Integrated Planning and Scheduling of Multiple “50’ Sail Catamaran with Hybrid Propulsion, Design,
Manufacturing Projects under Resource Constraints Theoretical and Experimental Studies,” Polish Maritime
Using Raccoon Family Optimization Algorithm,” IEEE Research, vol. 29, no. 2, 2022, doi: 10.2478/pomr-2022-0012.
Access, vol. 8, pp. 151279–151295, 2020, doi: 10.1109/
ACCESS.2020.2971650. 23. B. Branowski, M. Zabłocki, P. Kurczewski, and A. Walczak,
“Selected Issues in Universal Design of Yachts for People
13. R. Moser, G. Narayanamurthy, K. Kusaba, and G. Kaiser, with Disabilities,” Polish Maritime Research, vol. 28, no.
“Performance of low-cost country sourcing projects – 3, 2021, doi: 10.2478/pomr-2021-0030.
Conceptual model & empirical analysis,” Int J Prod Econ,
vol. 204, pp. 30–43, Oct. 2018, doi: 10.1016/j.ijpe.2018.07.011. 24. L. Venz, A. Pundt, and S. Sonnentag, “What matters for
work engagement? A diary study on resources and the
14. H. Abbasianjahromi and S. Hosseini, “A risk-cost benefits of selective optimization with compensation for
optimization model for selecting human resources in state work engagement,” J Organ Behav, vol. 39, no. 1, 2018,
construction projects,” SN Appl Sci, vol. 1, no. 11, 2019, doi: 10.1002/job.2207.
doi: 10.1007/s42452-019-1570-5.
25. J. Bin Yang and H. Y. Chou, “Subjective benefit evaluation
15. B. Ashuri and M. Tavakolan, “Shuffled Frog-Leaping Model model for immature BIM-enabled stakeholders,” Autom
for Solving Time-Cost-Resource Optimization Problems Constr, vol. 106, 2019, doi: 10.1016/j.autcon.2019.102908.
in Construction Project Planning,” Journal of Computing
in Civil Engineering, vol. 29, no. 1, Jan. 2015, doi: 10.1061/ 26. L. Rihar, T. Žužek, and J. Kušar, “How to successfully
(asce)cp.1943-5487.0000315. introduce concurrent engineering into new product
development?,” Concurr Eng Res Appl, vol. 29, no. 2, pp.
16. E. Acar, G. Bayrak, Y. Jung, I. Lee, P. Ramu, and S. S. 87–101, Jun. 2021, doi: 10.1177/1063293X20967929.
Ravichandran, “Modeling, analysis, and optimization
under uncertainties: a review,” Structural and 27. N. Bibi, Z. Anwar, and A. Ahsan, “Comparison of search-
Multidisciplinary Optimization, vol. 64, no. 5. 2021. doi: based software engineering algorithms for resource
10.1007/s00158-021-03026-7. allocation optimization,” Journal of Intelligent Systems,
vol. 2015, pp. 629–642, 2015, doi: 10.1515/jisys-2015-0016.
17. S. Saeed, M. A. Khan, and R. Ahmad, Business strategies
and approaches for effective engineering management. 2013. 28. L. K. Hansen and P. Svejvig, “Seven Decades of Project
doi: 10.4018/978-1-4666-3658-3. Portfolio Management Research (1950–2019) and
Perspectives for the Future,” Project Management
Journal, vol. 53, no. 3, pp. 277–294, 2022, doi:
10.1177/87569728221089537.
Abstract
Logistics centres are currently performing a key function in the development of countries through their ability to regulate
goods, markets, and transport. This is shown by the infrastructure, cost, goods flow, and quality of logistical services
provided by these centres. Nevertheless, in developing nations or regions with antiquated logistics infrastructure,
conventional logistics centres seem to struggle to manage the volume of commodities passing through them, resulting
in persistent congestion and an unsteady flow of goods inside these facilities. This issue poses a challenge to the
progress of any nation. The emergence of new technology offers a potential avenue to solve the problems inherent in
traditional logistics centres. Most prominently, four technologies (the Internet of Things (IoT), Blockchain, Big Data
and Cloud computing) are widely applied in traditional logistics centres. This work has conducted a thorough analysis
and evaluation of these new technologies in relation to their respective functions and roles inside a logistics centre.
Furthermore, this work proposes difficulties in applying new technologies to logistics centres related to issues such
as science, energy, cost, or staff qualifications. Finally, future development directions, related to expanding policies
in technological applications, or combining each country’s policies for the logistics industry, are carefully discussed.
Internet of Things (IoT) The implementation of sensors and devices that are networked through wireless networks and [10][11][12][13]
Internet-based interaction with the objective of enhancing the value of products and processes.
Blockchain The digital platform facilitates the safe storage and distribution of information across a collective of [14][15][16]
users via the creation of time-stamped, tamper-proof, and indefinitely lasting records. The system
comprises decentralised ledgers that store transactions as data blocks, which are interconnected
by a cryptographic pointer. The aforementioned system exhibits attributes such as distributed
consensus, enhanced security measures, traceability, verification, and transparency of information.
Cloud Computing The online service provides users with the ability to do rapid and streamlined calculations, without [17][18][10]
the need for establishing a tangible infrastructure. This technology facilitates the provision of
computer resources, including networks, servers, storage, applications, and services, with the
ability to access a network that is readily available, easily accessible, and beneficial. The outcome
is a more economically efficient and expeditious resolution with regard to operational platforms,
software, and infrastructures.
Big Data The effective facilitation of decision-making processes is achieved via the management of [19][20][5]
a substantial amount of data, defined by its high volume, rapid velocity, and diversity. This is
accomplished by using cutting-edge analytic approaches that are creative in nature.
The logistics industry is showing an increasing dependence of logistics systems. As a result, the whole manufacturing
on new technologies. Therefore, in order to fully capitalise on sector, including its logistics, is transitioning to a paradigm
the potential of this industry, it is imperative to develop the that is more adaptable and agile, making room for Industry
‘Logistics 4.0’ initiative, which aims to maximise the utilisation 4.0. According to Khatib et al. [22], the heart of Industry
of cutting-edge technologies and implement innovative 4.0 consists of four primary enabling technologies that will
advancements in the logistics field [21]. Governments should increase the adaptability of manufacturing and distribution
aggressively push its ‘Logistics Centre 4.0’ strategy, starting processes: robotics, Big Data, wireless networking, and
with the logistics centre, which serves as the beating heart inexpensive sensors. The interdependencies between the
various 4.0 technologies are depicted in Fig. 1.
In the technologies mentioned above, IoT was acknowledged literature, logistics centre management is a crucial component
as one of the most significant fields of future technology and it of the supply chain that has received increased attention,
is one of the key information and communication technologies with the aim of meeting the rising freight demand and the
(ICT). The IoT is quickly gaining ground in the context of increasingly high standards for logistics services [29]. In the
contemporary wireless telecommunications, particularly last few decades, the warehouse and logistics centre business
with the rapid development of wireless communication has experienced a remarkable development in major trading
technologies [23][24][25]. From an initial emphasis on cities. Planners and social scientists have expressed concerns
machine-to-machine communication and applications in over the social and environmental impacts associated
the ‘ubiquitous aggregation’ of data, the concept of IoT is with this increase. Researchers have shown that transport
continually changing. In other words, the IoT has generated infrastructure and links to supply chain partners are two of
vast amounts of data and various mathematical analytic the most important factors in luring developers to an area
methods can be used to continuously investigate the intricate [30][31]. Although the advantages of the effects and the role of
links between the transactions represented by this data. digitalisation in logistics are still not properly recognised, it
Without a doubt, IoT would be crucial to the deployment cannot be denied that digitalisation will be an essential step in
of smart logistics [26][27], which would fundamentally alter the field of logistics [32][33]; this is reflected in transportation
the design of the logistics system and the logistical operation and warehouse activities [34]. As a typical example, in the
mode. These changes are very important in determining past, real-time tracking of vehicles was performed using
a company’s competitiveness because logistics costs are GPS systems. Then, with the advent of 4.0 technology and
seen as a significant component of overall production costs. blockchain, real-time goods tracking systems using blockchain
Numerous businesses are thinking about how to operate their technology have been applied by many companies. From
warehouses more cost-effectively and efficiently, especially there, digital transformation frameworks are continuously
in light of recent advancements in supply chain and logistics being updated, based on emerging logistics activities. The
technology but there are also many businesses that are hesitant research by Junge et al. [35] presented a framework for digital
about applying smart technologies to a previously traditional transformation in logistics, as demonstrated in Fig. 2.
infrastructure [28]. According to the public and academic
The next problem for the necessary development of 4.0 consistently high in underdeveloped nations. Precisely for
technologies is the e-commerce industry. The quantity of this reason, the Warehouse Management System was born to
products that need to be handled in a logistics centre has resolve the backlog of manual operations at warehouses and
grown over the past ten years, to meet the development of the logistics centres. Warehouse Management Systems (WMS)
e-commerce industry, making warehouse operations more evolved as a result of the emergence of increasingly complex
complex. As a result, older manual methods are no longer tools and algorithms to run warehouses effectively in the
sufficient or practicable to handle this enormous amount 2000s [45]. A WMS has an information system that combines
of activity [36]-[37]. From there, freight and logistics firms software programmes to keep track of, regulate, control, and
may encounter large delays between arrival and departure, manage inventory levels, and optimise warehousing choices.
as a result of poor management of the logistics centre [38]. Order processing, order release, and master data are the three
Extra fees (fines for late transportation and additional main WMS functions. In addition to the above functions, the
expenditure for keeping the tractor-trailer driver on overtime) additional capabilities of WMS include receiving (inbound),
and orders being delivered late are the direct result of this putting things away, and warehouse control [46]. However,
[39]. Long lines of idle vehicles contribute to pollution [40]- WMS stopped meeting the needs of traditional logistics
[41], thus it is best to set up transport and storage facilities centres. For the e-commerce industry, specific characteristics
with standardised loading and warehousing efficiency. The related to the volume and quantity of goods still make it
primary objective of transportation and logistics centres difficult for WMS to operate processing and storage.
is to minimise the overall costs associated with product By recognising the importance of technologies in
transportation, while simultaneously increasing storage applying the management of logistics centres to a smart
capacity through the implementation of short-term strategies environment, this work analyses and makes judgements
and the use of technology. These elements play a crucial role in of Industry 4.0 applications in the intelligent management
the efficient delivery of goods and enhancing transportation of flogistics centres. Applications that reduce the number of
and logistics systems [42]-[43]. Moreover, the absence of pointless procedures can lighten the load in a logistics centre.
cohesion in the managerial procedures within logistics centres Efficiency and operational productivity are then improved.
leads to frequent instances of time wastage and mistakes in These applications also pave the way for future studies in the
the tasks performed by personnel teams. Additionally, one field of smart logistics, where systems are combined with
company’s barcode system will differ from their suppliers’, cutting-edge technology. However, there will be some notable
leading to inconsistencies in the preservation of information obstacles, such as the cost of deploying the technology, the
about the items’ features. This has long been a problem in way in which the location of the logistics centre affects the
conventional logistics centres, where it slows down the arrival environment, and the policies for applying these technologies;
of transporting vehicles and drives up waiting times and these will be of interest in future research.
other expenses [44]. Because of this, the cost of logistics is
Fig. 7. The use of blockchain technology inside a supply chain framework [115]
1. L. Atzori, A. Iera, and G. Morabito, “The internet of 12. K. Witkowski, “Internet of Things, Big Data, Industry
things: A survey,” Comput. networks, vol. 54, no. 15, pp. 4.0 – Innovative Solutions in Logistics and Supply Chains
2787–2805, 2010. Management,” Procedia Eng., vol. 182, pp. 763–769, 2017,
doi: 10.1016/j.proeng.2017.03.197.
2. D. Bandyopadhyay and J. Sen, “Internet of Things:
Applications and Challenges in Technology and 13. H. P. Nguyen, P. Q. P. Nguyen, D. K. P. Nguyen, V. D.
Standardization,” Wirel. Pers. Commun., vol. 58, no. 1, Bui, and D. T. Nguyen, “Application of IoT Technologies
pp. 49–69, May 2011, doi: 10.1007/s11277-011-0288-5. in Seaport Management,” JOIV Int. J. Informatics Vis.,
vol. 7, no. 1, p. 228, Mar. 2023, doi: 10.30630/joiv.7.1.1697.
3. C. Bardaki, P. Kourouthanassis, and K. Pramatari,
“Deploying RFID-Enabled Services in the Retail Supply 14. A. Bekrar, A. Ait El Cadi, R. Todosijevic, and J. Sarkis,
Chain: Lessons Learned toward the Internet of Things,” “Digitalizing the Closing-of-the-Loop for Supply
Chains: A Transportation and Blockchain Perspective,”
20. D. Zhu, “IOT and big data based cooperative logistical 32. Q. Yuan, “Environmental Justice in Warehousing
delivery scheduling method and cloud robot system,” Location,” J. Plan. Lit., vol. 33, no. 3, pp. 287–298, Aug.
Futur. Gener. Comput. Syst., vol. 86, pp. 709–715, Sep. 2018, doi: 10.1177/0885412217753841.
2018, doi: 10.1016/j.future.2018.04.081.
33. Q. Yuan, “Location of Warehouses and Environmental
21. T. T. Le et al., “Management strategy for seaports aspiring Justice,” J. Plan. Educ. Res., vol. 41, no. 3, pp. 282–293,
to green logistical goals of IMO: Technology and policy Sep. 2021, doi: 10.1177/0739456X18786392.
solutions,” Polish Marit. Res., vol. 30, no. 2, pp. 165–187,
2023, doi: 10.2478/pomr-2023-0031. 34. Q. Yuan and J. Zhu, “Logistics sprawl in Chinese
metropolises: Evidence from Wuhan,” J. Transp.
22. E. J. Khatib and R. Barco, “Optimization of 5G Networks Geogr., vol. 74, pp. 242–252, Jan. 2019, doi: 10.1016/j.
for Smart Logistics,” Energies, vol. 14, no. 6, p. 1758, Mar. jtrangeo.2018.11.019.
2021, doi: 10.3390/en14061758.
35. F. Straube, A. L. Junge, P. Verhoeven, J. Reipert, and M.
23. J.-Q. Li, F. R. Yu, G. Deng, C. Luo, Z. Ming, and Q. Mansfeld, Pathway of digital transformation in logistics :
Yan, “Industrial Internet: A Survey on the Enabling best practice concepts and future developments, no. July.
Technologies, Applications, and Challenges,” IEEE 2019. doi: 10.14279/depositonce-8502.
Commun. Surv. Tutorials, vol. 19, no. 3, pp. 1504–1526,
2017, doi: 10.1109/COMST.2017.2691349. 36. R. Narasimhan and S. W. Kim, “INFORMATION
SYSTEM UTILIZATION STRATEGY FOR SUPPLY
24. M. Liu, F. R. Yu, Y. Teng, V. C. M. Leung, and M. Song, CHAIN INTEGRATION,” J. Bus. Logist., vol. 22, no. 2, pp.
“Performance Optimization for Blockchain-Enabled 51–75, Sep. 2001, doi: 10.1002/j.2158-1592.2001.tb00003.x.
Industrial Internet of Things (IIoT) Systems: A Deep
Reinforcement Learning Approach,” IEEE Trans. Ind. 37. S. Kot, I. Goldbach, and B. Ślusarczyk, “Supply chain
Informatics, vol. 15, no. 6, pp. 3559–3570, Jun. 2019, doi: management in SMEs – Polish and Romanian approach,”
10.1109/TII.2019.2897805. Econ. Sociol., vol. 11, no. 4, pp. 142–156, Dec. 2018, doi:
10.14254/2071-789X.2018/11-4/9.
25. C. Qiu, F. R. Yu, H. Yao, C. Jiang, F. Xu, and C. Zhao,
“Blockchain-Based Software-Defined Industrial Internet
61. V. D. Bui and H. P. Nguyen, “A Systematized Review on 73. V. Sharma, I. You, G. Pau, M. Collotta, J. Lim, and J.
Rationale and Experience to Develop Advanced Logistics Kim, “LoRaWAN-Based Energy-Efficient Surveillance by
Center System in Vietnam,” Webology, vol. 18, pp. 89–101, Drones for Intelligent Transportation Systems,” Energies,
2021. vol. 11, no. 3, p. 573, Mar. 2018, doi: 10.3390/en11030573.
62. H. Makino, K. Tamada, K. Sakai, and S. Kamijo, “Solutions 74. R. H. Weber and R. Weber, “Governance of the Internet
for urban traffic issues by ITS technologies,” IATSS of Things,” in Internet of Things, Berlin, Heidelberg:
Res., vol. 42, no. 2, pp. 49–60, Jul. 2018, doi: 10.1016/j. Springer Berlin Heidelberg, 2010, pp. 69–100. doi:
iatssr.2018.05.003. 10.1007/978-3-642-11710-7_4.
63. S. Lee, Y. Kang, and V. V. Prabhu, “Smart logistics: 75. X. Feng, L. T. Yang, L. Wang, and A. Vinel, “Internet of
distributed control of green crowdsourced parcel services,” things,” Int. J. Commun. Syst., vol. 25, no. 9, pp. 1101–1102,
Int. J. Prod. Res., vol. 54, no. 23, pp. 6956–6968, Dec. 2016, 2012.
doi: 10.1080/00207543.2015.1132856.
76. Robert Sheldon, “Definition Sensor,” 2022.
64. A. Kawa, “SMART Logistics Chain,” in Intelligent
Information and Database Systems, 2012, pp. 432–438. 77. M. Wang, C. C. Wang, S. Sepasgozar, and S. Zlatanova,
doi: 10.1007/978-3-642-28487-8_45. “A Systematic Review of Digital Technology Adoption
in Off-Site Construction: Current Status and Future
65. H. P. Nguyen, “Sustainable development of logistics in Direction towards Industry 4.0,” Buildings, vol. 10, no.
Vietnam in the period 2020-2025,” Int. J. Innov. Creat. 11, p. 204, Nov. 2020, doi: 10.3390/buildings10110204.
Chang., vol. 11, no. 3, pp. 665–682, 2020.
78. C. Mao, X. Tao, H. Yang, R. Chen, and G. Liu, “Real-
66. J. Mehmann and F. Teuteberg, “Process reengineering by time carbon emissions monitoring tool for prefabricated
using the 4PL approach,” Bus. Process Manag. J., vol. 22, no. construction: An iot-based system framework,” in
4, pp. 879–902, Jul. 2016, doi: 10.1108/BPMJ-12-2014-0119. International Conference on Construction and Real Estate
Management 2018, 2018, pp. 121–127.
67. Y. M. Tang, G. T. S. Ho, Y. Y. Lau, and S. Y. Tsui, “Integrated
Smart Warehouse and Manufacturing Management with 79. G. Liu et al., “Cyber-physical system-based real-time
Demand Forecasting in Small-Scale Cyclical Industries,” monitoring and visualization of greenhouse gas emissions
Machines, vol. 10, no. 6, 2022. of prefabricated construction,” J. Clean. Prod., vol. 246,
p. 119059, Feb. 2020, doi: 10.1016/j.jclepro.2019.119059.
68. V. Yavas and Y. D. Ozkan-Ozen, “Logistics centers in the
new industrial era: A proposed framework for logistics 80. Y. Zhai et al., “An Internet of Things-enabled BIM platform
center 4.0,” Transp. Res. Part E Logist. Transp. Rev., vol. for modular integrated construction: A case study in Hong
135, p. 101864, Mar. 2020, doi: 10.1016/j.tre.2020.101864. Kong,” Adv. Eng. Informatics, vol. 42, p. 100997, Oct. 2019,
doi: 10.1016/j.aei.2019.100997.
69. S. Winkelhaus and E. H. Grosse, “Logistics 4.0:
a systematic review towards a new logistics system,” 81. M. Valinejadshoubi, A. Bagchi, and O. Moselhi,
Int. J. Prod. Res., vol. 58, no. 1, pp. 18–43, Jan. 2020, doi: “Development of a BIM-Based Data Management System
10.1080/00207543.2019.1612964. for Structural Health Monitoring with Application
to Modular Buildings: Case Study,” J. Comput. Civ.
70. G. S. Cho, “A Study on Establishment of Smart Logistics Eng., vol. 33, no. 3, May 2019, doi: 10.1061/(ASCE)
Center based on Logistics 4.0,” J. Multimed. Inf. Syst., vol. CP.1943-5487.0000826.
5, no. 4, pp. 265-272., 2018, doi: http://dx.doi.org/10.9717/
JMIS.2018.5.4.265. 82. my Murata, “Barometric Pressure Sensor Basics.”
71. W. Hamdy, N. Mostafa, and H. Elawady, “Towards a smart 83. X. Li, H. Chi, P. Wu, and G. Q. Shen, “Smart work
warehouse management system,” Proc. Int. Conf. Ind. Eng. packaging-enabled constraint-free path re-planning for
Oper. Manag., vol. 2018, no. SEP, pp. 2555–2563, 2018. tower crane in prefabricated products assembly process,”
Adv. Eng. Informatics, vol. 43, p. 101008, Jan. 2020, doi:
72. H. CHOW, K. CHOY, W. LEE, and K. LAU, “Design 10.1016/j.aei.2019.101008.
of a RFID case-based resource management system for
112. M. Zheng and K. Wu, “Smart spare parts management 123. Accenture, “The state of of supply The state supply chain
systems in semiconductor manufacturing,” Ind. Manag. and the cloud chain and the cloud,” 2022.
Data Syst., vol. 117, no. 4, pp. 754–763, May 2017, doi:
10.1108/IMDS-06-2016-0242. 124.Maersk, “How cloud computing is shaping the future of
logistics,” 2023.
113. S. N. Wahab, M. I. Hamzah, N. M. Sayuti, W. C. Lee,
and S. Y. Tan, “Big data analytics adoption: an empirical 125. Y. Jiang, “Location and Path Planning of Cross-Border
study in the Malaysian warehousing sector,” Int. J. Logist. E-Commerce Logistics Distribution Center in Cloud
Syst. Manag., vol. 40, no. 1, p. 121, 2021, doi: 10.1504/ Computing Environment,” 2021, pp. 30–40. doi:
IJLSM.2021.117703. 10.1007/978-3-030-67871-5_4.
114. R. W. Ahmad, H. Hasan, R. Jayaraman, K. Salah, and 126.X. Zhang, “Optimization design of railway logistics
M. Omar, “Blockchain applications and architectures for center layout based on mobile cloud edge computing,”
port operations and logistics management,” Res. Transp. PeerJ Comput. Sci., vol. 9, p. e1298, Apr. 2023, doi: 10.7717/
Bus. Manag., vol. 41, p. 100620, Dec. 2021, doi: 10.1016/j. peerj-cs.1298.
rtbm.2021.100620.
127. P. Sharma and S. Panda, “Cloud Computing for Supply
115. S. Johar, N. Ahmad, W. Asher, H. Cruickshank, and A. Chain Management and Warehouse Automation: A Case
Durrani, “Research and Applied Perspective to Blockchain Study of Azure Cloud,” Int. J. Smart Sens. Adhoc Network.,
Technology: A Comprehensive Survey,” Appl. Sci., vol. 11, pp. 19–29, Jan. 2023, doi: 10.47893/IJSSAN.2023.1227.
no. 14, p. 6252, Jul. 2021, doi: 10.3390/app11146252.
128.B. Gupta, A. Gupta, and A. Nagpal, “Implementation of
116. V. Aleksieva, H. Valchanov, A. Haka, and D. Dinev, Warehouse Management Through Cloud Computing,”
“Logistics Model Based on Smart Contracts on Blockchain IJIRST-International J. Innov. Res. Sci. Technol., vol. 3, no.
and IoT,” in EEPES’23, Jul. 2023, p. 8. doi: 10.3390/ 10, pp. 189–192, 2017.
engproc2023041008.
129. V. Sivakumar, R. Ruthramathi, and S. Leelapriyadharsini,
117. J. Chen, S. Xu, K. Liu, S. Yao, X. Luo, and H. Wu, “Intelligent “Challenges of Cloud Computing in Warehousing
Transportation Logistics Optimal Warehouse Location Operations with Respect to Chennai Port Trust,” in
Method Based on Internet of Things and Blockchain Proceedings of the 2020 the 3rd International Conference
Technology,” Sensors, vol. 22, no. 4, p. 1544, Feb. 2022, on Computers in Management and Business, Jan. 2020,
doi: 10.3390/s22041544. pp. 162–165. doi: 10.1145/3383845.3383895.