DOI: 10.24937/2542-2324-2022-2-400-59-66 УДК 629.544
А.В. Сверчков1 , M.B. Галушина1, Д. Вастаручас2, П. Полити2
1 ФГУП «Крыловский государственный научный центр», Санкт-Петербург, Россия
2 Danaos Shipping Co. Ltd., Греция
МОДИФИКАЦИЯ КОРПУСОВ КОНТЕЙНЕРОВОЗОВ ПРИ СНИЖЕНИИ СКОРОСТИ ХОДА С ЦЕЛЬЮ СОКРАЩЕНИЯ ОБЪЕМА ПАРНИКОВЫХ ГАЗОВ
Объект и цель научной работы. Объектом исследования являются два судна-контейнеровоза греческой судоходной компании Danaos Shipping Co. Ltd. контейнеровместимостью 8100 и 8500 TEU. Цель работы состоит в экспериментальной проверке модификации формы обводов носовых оконечностей контейнеровозов, направленной на уменьшение гидродинамического сопротивления при снижении эксплуатационной скорости, а также в определении влияния модификации на энергетическую эффективность судна.
Материалы и методы. Экспериментальные данные модельных испытаний, полученные в глубоководном опы-товом бассейне ФГУП «Крыловский государственный научный центр».
Основные результаты. В результате исследовательских работ получены зависимости потребной мощности от скорости хода контейнеровозов с исходной и модифицированной формой обводов носовой оконечности. Заключение. За счет перехода на более низкую скорость хода с одновременной оптимизацией и модернизацией формы обводов носовой оконечности возможное снижение потребной мощности, а соответственно, расхода топлива и объема выбросов углекислого газа, составляет примерно 40 %.
Ключевые слова: контейнеровоз, носовой бульб, модификация носовой оконечности, модельные испытания, индекс энергетической эффективности EEID.
Авторы заявляют об отсутствии возможных конфликтов интересов.
DOI: 10.24937/2542-2324-2022-2-400-59-66 UDC 629.544
A.V. Sverchkov1 , M.V. Galushina1, D. Vastarouchas2, P. Politi2
1 Krylov State Research Centre, St. Petersburg, Russia
2 Danaos Shipping Co. Ltd., Greece
MODIFICATION OF CONTAINERSHIP HULLS AT SPEEDS REDUCED FOR CUTTING DOWN GREENHOUSE GASES
Object and purpose of research. The objects of study are two containerships of 8100 and 8500 TEU capacities operated by Danaos Shipping Co. Ltd, Greece. The purpose is to experimentally verify alteration of bow shapes in the containerships designed to reduce hydrodynamic resistance at a lower operation speed, as well as to examine the impact of these alterations on energy efficiency design index (EEDI).
Materials and methods. Experimental data from model tests were obtained in a deep-water test tank of Krylov State Research Centre.
Main results. Based on the research studies the shaft power was determined versus the containership speed for the hulls with the original and modified bow shapes.
Conclusion. Switchover to lower speeds combined with optimization and modernization of bow shapes will bring a possible reduction of shaft power and, accordingly, less fuel consumption and carbon dioxide emissions by about 40 %. Keywords: containership, bow bulb, bow modification, model tests, energy efficiency design index EEDI. The authors declare no conflicts of interest.
Для цитирования: Сверчков А.В., Галушина М.В., Вастаручас Д., Полити П. Модификация корпусов контейнеровозов при снижении скорости хода с целью сокращения объема парниковых газов. Труды Крыловского государственного научного центра. 2022; 2(400): 59-66.
For citations: Sverchkov A.V., Galushina M.V., Vastarouchas D., Politi P. Modification of containership hulls at speeds reduced for cutting down greenhouse gases. Transactions of the Krylov State Research Centre. 2022; 2(400): 59-66 (in Russian).
Введение
Introduction
История вопроса о выбросах парниковых газов и связанных с этим последствиях насчитывает более 15 лет. Современные требования Международной морской организации (IMO) включают в себя требования по снижению объема парниковых газов в судоходстве к 2050 г. не менее чем на 50 %. При этом главной целью становится достижение, по возможности, минимальных выбросов углекислого газа, зависящих напрямую от количества потребляемого топлива.
Введенный IMO индекс энергетической эффективности (Energy Efficiency Design Index - EEID) определяется как отношение массы выбрасываемого углекислого газа к транспортной работе судна, определяемой мощностью энергетической установки, скоростью хода и грузовместимостью. Возможные пути минимизации индекса энергетической эффективности включают в себя: снижение скорости, снижение потребной мощности за счет оптимизации формы обводов корпуса, установку энергосберегающих устройств, оптимизацию движи-тельного комплекса, смену типа топлива, снижение удельного расхода топлива, увеличение вместимости при сохранении общего водоизмещения судна.
G g, т/сут
Контейнеровместимость, TEU
---- 4000-5000 -- 8000-9000
---- 5000-6000 ------ 9000-10 000
----- 7000-8000 - 10 000+
Рис. 1. График зависимости расхода топлива от скорости судна при различной контейнеровместимости
Fig. 1. Fuel consumption versus ship speed at different container capacity
Потребление топлива зависит от размеров судна и его скорости. Следует отметить, что для контейнеровозов в диапазоне от минимальных скоростей до 14 уз эта зависимость носит характер близкий к линейному, но при скоростях более 14 уз расход топлива начинает увеличиваться экспоненциально. Например, для контейнеровоза вместимостью 8000 TEU с эксплуатационной скоростью 24 уз расход тяжелого топлива составляет 225 т в сутки. При снижении скорости на 3 уз (на 12,5 %) потребление падает до 150 т в сутки (т.е. на 33 %). На рис. 1 показаны зависимости расхода топлива Ge от скорости судна VS при различной контейнеровместимости [1].
Важнейшими параметрами, влияющими на эффективность энергопотребления судами, являются гидродинамическое сопротивление и КПД движителей. Учитывая, что каждое построенное судно на этапе проектирования проходит тщательную проработку в гидродинамических центрах, направленную на численную и экспериментальную оптимизацию формы обводов корпуса и гидродинамических характеристик движительного комплекса, снижение показателей энергоэффективности на 50 % является затруднительным [2]. Единого решения, подходящего для всех типов судов, на данный момент не существует. Каждое судно требует индивидуальной проработки. Поэтому многие судовладельцы рассматривают снижение скорости судна как наиболее простой и надежный способ минимизации выбросов. При этом одновременно встает задача оптимизации формы обводов корпуса под более низкую эксплуатационную скорость. Для крупных судовладельческих компаний эта задача сводится к оптимизации и модернизации формы обводов уже построенных и находящихся в эксплуатации судов. В частности, одна из крупнейших судоходных компаний, греческая корпорация Danaos Shipping Co. Ltd., владеющая современными контейнеровозами вместимостью от 2200 до 13 100 TEU, в 2015 г. начала работы по модернизации своего флота.
Исследуемые суда
Containerships under study
В рамках двух контрактных работ в Крыловском государственном научном центре были выполнены экспериментальные исследования возможной модификации формы обводов носовых оконечностей контейнеровозов CSCL Europe и CMA CGM Attila. Основные технические характеристики исследо-
ванных контейнеровозов приведены в табл. 1. Фотографии судов показаны на рис. 2 и 3.
Оба контейнеровоза имеют бульбообразную носовую оконечность, спроектированную и оптимизированную на проектную скорость 21 уз. При снижении эксплуатационной скорости хода судна
до 18 уз, что соответствует переходу от числа Фру-да Fn = 0,194 к Fn = 0,166, картина волнообразования существенно меняется, поэтому возникает необходимость перепроектирования носового буль-ба, основной задачей которого является снижение волнового сопротивления.
Таблица 1. Основные технические характеристики исследованных контейнеровозов Table 1. Main data of the containerships under study
Наименование CSCL Europe CMA CGM Attila
Год постройки 2004 2011
Судоверфь Корея, Samsung Heavy Industries Co. Ltd Китай, Shanghai Jiangnan Changxing Heavy Industry
Флаг Кипр Мальта
Длина наибольшая, м 333,99 321,46
Ширина B, м 42,8 42,8
Проектная осадка T, м 13,0 13,0 12,0
Длина по ватерлинии LWL, м 316,6 317,6 314,0
Объемное водоизмещение VWL, м3 114 050 114 217 103213
Площадь смоченной поверхности S, м2 16 755 16 569 15 697
Коэффициент общей полноты CB 0,647 0,646 0,640
LWL /B 7,397 7,421 7,336
B/T 3,292 3,292 3,567
Lwl/ о .. . 6,528 6,546 6,694
Контейнеровместимость, EU 8100 8500
Максимальная мощность ЭУ (MCR), кВт 68 470 68 520
Частота вращения гребного винта, об/мин 104 104
Максимальная скорость, уз 22 25,8
Рис. 2. CSCL Europe Fig. 2. CSCL Europe
Рис. 3. CMA CGM Attila Fig. 3. CMA CGM Attila
2. СРЯ = ЬРК /ЬРР - выступающая за носовой перпендикуляр длина бульба ЬРР, отнесенная к длине судна между перпендикулярами ЬРР;
3. С2В = ZB /Тр - высота самой выступающей вперед точки бульба над основной плоскостью 2В, отнесенная к осадке;
4. САВТ = АВТ /Ам - площадь носового бульба на носовом перпендикуляре, отнесенная к площади мидель-шпангоута судна А^;
5. С^ = Аж /AMS - площадь носового бульба в диаметральной плоскости А^, отнесенная к площади мидель-шпангоута судна А^;
6. С,
VPR
■VpR / VWL
Носовой перпендикуляр
Рис. 4. Схема определения параметров бульба Fig. 4. Scheme for defining bulb parameters
Для экспериментальной проверки гидродинамических качеств корпусов контейнеровозов с различными вариантами носового бульба были изготовлены крупномасштабные модели судов длиной около 9,0 м со съемными носовыми оконечностями. На модели судна CSCL Europe были исследованы два варианта носового бульба (исходный и численно оптимизированный); на модели судна CMA CGM Attila экспериментально проверены три варианта (исходный и два варианта, оптимизированные по CFD разными гидродинамическими центрами).
Для описания геометрии носовых бульбов использованы параметры, предложенные в работе Альфреда Крахта [3]:
1. CBB = Bb/Bms - максимальная ширина площади бульба на носовом перпендикуляре BB, отнесенная к ширине судна BMS;
- объем бульба, выступающий за носовой перпендикуляр VPR, отнесенный к объемному водоизмещению судна VWL. Схема определения параметров бульба приведена на рис. 4. В табл. 2 приведены геометрические параметры рассмотренных вариантов носового бульба. На рис. 5 и 6 (см. вклейку) приведены фотографии моделей с различными вариантами носового бульба.
Испытания в опытовом бассейне
Tests in deep-water towing tank
Буксировочные испытания моделей проводились в глубоководном опытовом бассейне ФГУП «Кры-ловский государственный научный центр» в соответствии с методикой, основанной на применении широкоизвестной схемы Фруда [4] и близкой к методике, рекомендованной МКОБ [5].
Во время рабочего пробега измеряются одновременно следующие величины: скорость модели VM, равная скорости буксировочной тележки; полное сопротивление модели RTM; динамический дифферент; динамическое всплытие/погружение.
Таблица 2. Геометрические параметры рассмотренных вариантов носового бульба Table 2. Geometric parameters of the bow bulb options under consideration
CSCL Europe CMA CGM Attila
Параметр Исходный вариант Оптимизированный вариант Исходный вариант Оптимизированный вариант № 1 Оптимизированный вариант № 2
cbb 0,135 0,105 0,129 0,096 0,068
clpr 0,030 0,030 0,029 0,031 0,029
czb 0,679 0,679 0,668 0,743 0,668
cabt 0,076 0,055 0,072 0,065 0,037
cabl 0,132 0,123 0,128 0,144 0,128
c v/'K 0,0025 0,0015 0,0023 0,0019 0,0010
а) б)
Рис. 5. Модель контейнеровоза CSCL Europe: а) исходный вариант; б) оптимизированный вариант Fig. 5. Model of the containership CSCL Europe: a) original version; b) optimized version
Рис. 6. Модель контейнеровоза CM A CGM Attila:
а) исходный вариант;
б) оптимизированный вариант № 1;
в) оптимизированный вариант № 2
Fig. 6. Model of the containership CM A CGM Attila:
a) original version;
b) optimized version No. 1;
c) optimized version No. 2
б)
Рис. 10. Модель контейнеровоза CSCL Europe во время буксировочных испытаний:
а) исходный вариант; б) оптимизированный вариант. Скорость буксировки соответствует скорости натурного судна 18 уз Fig. 10. Containership CSCL Europe model during towing tests:
a) original version; b) optimized version. Towing speed corresponds to 18 knots full-scale speed
Исходный вариант Оптимизированный вариант № 2
б)
Рис. 11. Модель контейнеровоза СМА CGM Attila во время буксировочных испытаний:
а) при осадке Т= 12,0 м; б) при осадке Т= 13,0 м. Скорость буксировки соответствует скорости натурного судна 18 уз Fig. 11. Conta nership СМА CGM Attila model during towing tests:
a) at draft T = 12.0 m; b) at draft T = 13.0 m. Towing speed corresponds to 18 knots full-scale speed
Коэффициент полного сопротивления определяется выражением (1), в котором нижний индекс «М» означает, что данная величина соответствует модели.
Ктлл
C
TM
2 • Р-^
• S,
M
C
F (lgRn)2'58
Cr-10 1,5
Рис. 7. Зависимость коэффициента остаточного сопротивления от числа Фруда, полученная по результатам буксировочных испытаний модели контейнеровоза CSCL Europe
Fig. 7. Residual resistance coefficient versus Froude number obtained from the towing test data for the containership CSCL Europe model
1,0
0,5
Число Рейнольдса определяется как
Rn = Vm'LWLM /vm.
(3)
(1)
В выражении (1) р представляет собой плотность воды в опытовом бассейне, SM - площадь смоченной поверхности корпуса.
Измеренное полное сопротивление модели разделяется на две составляющие: сопротивление трения эквивалентной пластины и остаточное сопротивление. Коэффициент сопротивления трения эквивалентной гладкой пластины рассчитывается по формуле Пранд-тля - Шлихтинга, при этом длина эквивалентной гладкой пластины равна длине корпуса модели по ватерлинии в статическом положении
0,455
В выражении (3) vМ - кинематическая вязкость воды, которая принимается в зависимости от температуры в соответствии с рекомендациями МКОБ [6].
Коэффициент остаточного сопротивления CR и число Фруда вычисляются для каждого рабочего пробега (для каждой скорости буксировки) с использованием следующих выражений:
Fn =
Vm
yjg • LW1
IM
CR = CTM
C
F 0M-
(4)
(5)
(2)
Полученные в ходе буксировочных испытаний зависимости коэффициента остаточного сопротивления от числа Фруда показаны на рис. 7-9. Для модели судна CMA CGM Attila приведены результаты буксировочных испытаний на двух осадках.
: : 1 : : : : 1 : : : 1 :
□ Исходный оульо на судне СÙLL Europe О Оптимизированный вариант бульба
* ...i...
о п о— .... ...
о
0,05
0,10
0,15
0,20
0,25
Fn
Q.10 1,0
Рис. 8. Зависимость коэффициента остаточного сопротивления от числа Фруда, полученная по результатам буксировочных испытаний модели контейнеровоза CMA CGM Attila при осадке T = 12,0 м
Fig. 8. Residual resistance coefficient versus Froude number obtained from the towing test data for the containership CMA CGM Attila model at draft T = 12.0 m
0,5
: i : : : : i : : : 1 i ____ О Исходный бульб на судне CMA CGM Attila А Оптимизированный вариант бульба № 1
—- __ :___ À___
г, о: -«i —- i 15
о О "
.... Л л
П-Äh j
р а ;
.... i
0,10
0,15
0,20
Fn
Рис. 9. Зависимость коэффициента остаточного сопротивления от числа Фруда, полученная по результатам буксировочных испытаний модели контейнеровоза CMA CGM Attila при осадке T = 13,0 м
Fig. 9. Residual resistance coefficient versus Froude number obtained from the towing test data for the containership CMA CGM Attila model at draft T = 13.0 m
Остаточное сопротивление включает в себя так называемое сопротивление формы и волновое сопротивление. Коэффициент сопротивления формы определяется на доволновом режиме, при значениях числа Фруда меньше 0,1. Полученные кривые коэффициента остаточного сопротивления показывают, что оптимизация носового бульба ведет к снижению и сопротивления формы, и волнового сопротивления. На обоих корпусах положительный эффект от снижения объема бульба наблюдается в широком диапазоне скоростей от малого хода до 22 уз. Из двух оптимизированных численным методом вариантов бульба на контейнеровозе CMA CGM Attila лучшим является вариант № 2, отличающийся меньшим объемом и меньшей площадью сечения на носовом перпендикуляре. При этом эффект снижения остаточного сопротивления за счет
АPs, %
10 15 20 25 Vs,y3
Рис. 12. Снижение потребной мощности за счет оптимизации носового бульба на контейнеровозе CSCL Europe
Fig. 12. Reduced shaft power s due to the bow bulb optimization of the containership CSCL Europe
оптимизации носового бульба особенно заметен на контейнеровозе CMA CGM Attila при осадке, равной 12,0 м, при которой исходный бульб с большим объемом и площадью сечения приближен к свободной поверхности воды.
Рис. 10 и 11 (см. вклейку) служат иллюстрацией обтекания носовой оконечности корпусов контейнеровозов при буксировке моделей со скоростью, соответствующей скорости натурного судна 18 уз. На фотографиях отчетливо видно изменение волнообразования в сторону улучшения при установке на модель носового бульба, оптимизированного под меньшую эксплуатационную скорость. Отсутствует неблагоприятная интерференция волн, вызванных корпусом и собственно бульбом, носовая волна выравнивается, уменьшается глубина волновой впадины в районе 19,5 теоретического шпангоута. За счет уменьшения объема бульба существенно снижается величина остаточного сопротивления в диапазоне чисел Фруда Fn от 0,10 до 0,18.
Циклы испытаний со сменными носовыми бульбами включали в себя, кроме буксировочных, самоходные испытания, направленные на получение коэффициентов взаимодействия гребного винта с корпусом. Самоходные испытания проводились со специально изготовленными моделями гребных винтов, геометрия которых соответствует натурным гребным винтам, установленным на судах.
Результаты оптимизации обводов носовой оконечности
Results of bow shape optimization
По окончании испытаний были выполнены оценки ходовых качеств контейнеровозов и получены зависимости потребной мощности от скорости хода. Диаграммы на рис. 12 и 13 иллюстрируют получен-
10 15 20 25 Vs, уз ю 15 20 25 Vs, уз
а) б)
Рис. 13. Снижение потребной мощности за счет оптимизации носового бульба на контейнеровозе CMA CGM Attila: а) при осадке 12,0 м; б) при осадке 13,0 м. Сопоставления исходного варианта с оптимизированным вариантом № 2
Fig. 13. Reduced shaft power due to the bow bulb optimization of the containership CMA CGM Attila: а) at draft 12.0 m; b) at draft 13.0 m. Comparison of the original version with version No. 2
ное за счет оптимизации корпуса снижение потребной мощности, отрицательные значения соответствуют снижению потребной мощности, положительные - возрастанию. Максимальный положительный эффект 16,9 % достигнут при скорости 16 уз для контейнеровоза CMA CGM Attila при осадке, равной 12,0 м.
При скорости 18 уз снижение потребной мощности составило:
■ 5,2 % для контейнеровоза CSCL Europe;
■ 11,9 % для контейнеровоза CMA CGM Attila при осадке, равной 12,0 м;
■ 4,1 % для контейнеровоза CMA CGM Attila при осадке, равной 13,0 м.
В табл. 3 приведены значения удельной потребной мощности PS /D (отношение потребной мощности к весовому водоизмещению судна), полученные для всех испытанных вариантов.
Таблица 3. Удельная потребная мощность Table 3. Specific shaft power
Анализ полученных результатов показал, что при снижении скорости контейнеровозов с 21 до 18 уз потребная мощность снижается на 30-35 %, при одновременной оптимизации формы носового бульба снижение потребной мощности составляет 37-39 %.
Заключение
Conclusion
На сегодняшний день для судовладельцев наиболее доступным и эффективным способом выполнения требований IMO по минимизации выбросов углекислого газа и величины индекса энергоэффективности EEDI является снижение скорости хода судна. Возникающая при этом проблема несоответствия геометрии корпуса судна новому режиму эксплуатации частично может быть реше-
PS ID, кВт/т
Vs, уз CSCL Europe, T = 13,0 м CMA CGM Attila, T = 12,0 м CMA CGM Attila, T = 13,0 м
Исходный вариант Оптимизированный вариант Исходный вариант Оптимизированный вариант № 1 Оптимизированный вариант № 2 Исходный вариант Оптимизированный вариант № 2
18,0 0,148 0,141 0,174 0,163 0,153 0,150 0,144
21,0 0,229 0,222 0,251 0,244 0,238 0,228 0,226
на путем оптимизации формы обводов корпуса, в частности, за счет модификации носовой оконечности и выбора оптимальных параметров носового бульба. За счет перехода на более низкую скорость хода с одновременной оптимизацией и модернизацией формы обводов носовой оконечности у судовладельцев появляется возможность снизить потребную мощность, а соответственно расход топлива и объем выбросов углекислого газа, примерно на 40 %.
Список использованной литературы
1. Notterboom T., Cariou P. Fuel surcharge practices of container shipping lines: Is it about cost recovery or revenue-making // Proceedings of the International Association of Maritime Economists conference (IAME 2009). Basingstoke : Macmillan, 2009. P. 24-26.
2. Пустошный А.В. Проблемы ходкости транспортных судов. Санкт-Петербург : Крыловский государственный научный центр, 2016. 142 c.
3. Kracht A.M. Design of Bulbous Bows // SNAME Transaction. 1978. Vol. 86. P. 197-217.
4. Справочник по теории корабля : В 3 т. Т. 1. Гидромеханика. Сопротивление движению судов. Судовые движители / [ВойткунскийЯ.И. и др.]. Ленинград : Судостроение, 1985. 764 с.
5. ITTC : Recommended Procedures and Guidelines : 7.5-02-02-01. Resistance Test : Rev. 03 / Ed. 26th ITTC Resistance Committee. [S. l.]: ITTC, 2011. 13 p.
6. ITTC : Recommended Procedures and Guidelines : 7.5-02-01-03. Fresh Water and Seawater Properties : Rev. 02 / Ed. 26th ITTC Specialist Committee on Uncertainty Analysis. [S. l.]: ITTC, 2011. 45 p.
References
1. Notterboom T., Cariou P. Fuel surcharge practices of container shipping lines: Is it about cost recovery or revenue-making // Proceedings of the International Association of Maritime Economists conference (IAME 2009). Basingstoke : Macmillan, 2009. P. 24-26.
2. PustoshnyA.V. Problems in propulsion performance of transport vessels. St. Petersburg : Krylov Shipbuilding Research Institute, 2016. 142 p. (in Russian).
3. Kracht A.M. Design of Bulbous Bows // SNAME Transactions. 1978. Vol. 86. P. 197-217.
4. Handbook of ship theory : in 3 volumes. Vol. 1. Hydromechanics. Ship resistance. Ship propellers / [Ya.I. Voit-kunsky et al.]. Leningrad : Sudostroenie, 1985. 764 p. (in Russian).
5. ITTC : Recommended Procedures and Guidelines : 7.5-02-02-01. Resistance Test : Rev. 03 / Ed. 26th ITTC Resistance Committee. [S. l.]: ITTC, 2011. 13 p.
6. ITTC : Recommended Procedures and Guidelines : 7.5-02-01-03. Fresh Water and Seawater Properties : Rev. 02 / Ed. 26th ITTC Specialist Committee on Uncertainty Analysis. [S. l.]: ITTC, 2011. 45 p.
Сведения об авторах
Сверчков Андрей Владимирович, к.т.н., начальник сектора ФГУП «Крыловский государственный научный центр». Адрес: 196158, Россия, Санкт-Петербург, Московское шоссе, д. 44. Тел.: +7 (812) 748-63-26. E-mail: [email protected]. https://orcid.org/0000-0002-9006-8812. Галушина Маргарита Валерьевна, старший научный сотрудник ФГУП «Крыловский государственный научный центр». Адрес: 196158, Россия, Санкт-Петербург, Московское шоссе, д. 44. E-mail: [email protected]. Вастаручас Димитриос, заместитель начальника оперативного отдела и технический директор Danaos Shipping Co. Ltd., Греция.
Полити Параскеви, менеджер по научным исследованиям и развитию Danaos Shipping Co. Ltd., Греция.
About the authors
Andrey V. Sverchkov, Cand. Sci. (Eng.), Head of Sector of Krylov State Research Centre. Address: 44, Moskovskoe sh., St. Petersburg, Russia, post code 196158. Tel.: +7 (812) 748-63-26. E-mail: [email protected]. https://orcid.org/0000-0002-9006-8812. Margarita V. Galushina, Senior Researcher of Krylov State Research Centre. Address: 44, Moskovskoe sh., St. Petersburg, Russia, post code 196158. E-mail: [email protected]. Dimitrios Vastarouchas, DCOO and Technical Director of Danaos Shipping Co. Ltd., Greece.
Paraskevi Politi, R&D Manager of Danaos Shipping Co. Ltd., Greece.
Поступила / Received: 30.03.22 Принята в печать / Accepted: 04.05.22 © Коллектив авторов, 2022