<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>DIRECT COMPUTATIONAL EXPERIMENT IN STORM HYDRODYNAMICS OF MARINE OBJECTS</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Alexander Degtyarev</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Ivan Gankevich</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Vasily Khramushin</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Saint Petersburg State University</institution>
          ,
          <addr-line>13B Universitetskaya Emb., Saint Petersburg, 199034</addr-line>
          ,
          <country country="RU">Russia</country>
        </aff>
      </contrib-group>
      <pub-date>
        <year>2021</year>
      </pub-date>
      <fpage>5</fpage>
      <lpage>9</lpage>
      <abstract>
        <p>The paper presents and discusses a new computer toolkit for assessing the seaworthiness of a ship in stormy sailing conditions, intended for testing new design solutions for promising ocean-going ships and ships of unlimited ocean navigation, as well as organizing full-function simulators. The presented toolkit can be used by captains to select effective and safe sailing modes, as well as to train personnel. A modern computational experiment allows direct modeling of intense sea waves, which fully recreates the hydromechanics of external influences on the ship, and ensures its full positioning relative to the crests of storm waves. Such experiments require the use of high-performance computing technology based on GPGPU, as well as distributed data preparation and processing. Parameters of three structures of real sea waves, including extremely high heights, are determined as the initial conditions for modeling. The ship also possesses real dynamic characteristics, and she is capable to go by arbitrary courses in relation to waves. The possibility of dynamically changing the parameters of the ship during the experiment is provided. It is required to assess the motion state in various cases of ship loading, including emergency. The most important feature of the software is a full-fledged threedimensional visualization of all storm waves, as well as the spatial position and trajectory of the ship and its parameters. A special experimental environment for engineering surveys for projected ships is created using the graphical tools of OpеnGL. This computer toolkit is considered within the framework of the concept of a virtual testbed.</p>
      </abstract>
      <kwd-group>
        <kwd>hydromechanics</kwd>
        <kwd>direct computational experiment</kwd>
        <kwd>ship design</kwd>
        <kwd>efficiency of storm navigation of a ship</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>-</title>
      <p>ПРЯМОЙ ВЫЧИСЛИТЕЛЬНЫЙ ЭКСПЕРИМЕНТ В
ШТОРМОВОЙ ГИДРОДИНАМИКЕ КОРАБЛЯ,
ОСОБЕННОСТИ ПРОЯВЛЕНИЯ ОБЩЕКОРАБЕЛЬНОЙ
АРХИТЕКТУРЫ
А. Б. Дегтярев 1, И. В. Ганкевич 1, В. Н. Храмушин 1,2
1) Кафедра компьютерного моделирования и многопроцессорных систем, факультет
Прикладной математики – процессов управления СПбГУ.
2) Подсекция мореходных качеств судов в штормовых условиях, Российское
Научнотехническое общество судостроителей им. А.Н.Крылова</p>
      <p>E-mail: i.gankevich@spbu.ru
В работе представляется и обсуждается новый вычислительный инструментарий для
экспериментальной апробации мореходных качеств корабля в условиях штормового плавания,
предназначенный для использования на этапах проектирования перспективного океанского
флота, и который может быть полезен в навигационной практике для обоснования
эффективных и безопасных режимов хода в условиях интенсивного волнения, а также для
обучения штурманского состава безопасным и эффективным методам штормового
судовождения. Эффективность прямого вычислительного эксперимента определяется
включением главной части гидродинамического воздействия трохоидальных волн на корпус
корабля, без ограничений высоты морских волн и с полной динамикой корабля в
поступательных и вращательных движениях, включая дельфинирование, брочинг или
опрокидывание. Начальными и граничными условиями моделирования определяются
характеристики трёх групповых структур морского волнения, без ограничения крутизны волн.
Кинематическая модель перемещений и качки корпуса настроена на управляемое движение
корабля произвольными курсами относительно штормового волнения. Числовая модель
корабля обладает реальными динамическими характеристиками, с возможностью их изменения
в процессе эксперимента, что требуется для выявления аварийных режимов штормового хода, и
может востребоваться в предотвращении морских катастроф. Подобные вычислительные
эксперименты востребуют высокопроизводительные технологии на основе GPGPU, с
распределенной подготовкой и динамической обработкой массивов данных. Важной
особенностью является трехмерная визуализация средствами OpenGL в режиме реального
времени для всех штормовых волн, включая пространственное положение и траекторию
движения корабля, с числовым и графическим представлением процессов штормовой качки во
времени. Данный компьютерный инструментарий рассматривается в рамках концепции
виртуального морского полигона.
Ключевые слова: гидромеханика, прямой вычислительный эксперимент, проектирование
корабля, эффективность штормового плавания корабля.
1. Введение</p>
      <p>
        В работе представляется и обсуждается вычислительный инструментарий для оценки
мореходных качеств корабля в условиях штормового плавания, предназначенный к апробации
новых проектных решений для перспективных океанских кораблей и судов неограниченного
океанского плавания [
        <xref ref-type="bibr" rid="ref1">1</xref>
        ], и который может использоваться капитанами для подбора
эффективных и безопасных режимов плавания действующего флота, в том числе для снижения
заливаемости палуб, уменьшения размаха штормовой качки, и создания в целом комфортных
условий обитания и сохранности грузов на борту корабля.
      </p>
      <p>
        Основное внимание в работе обращается на построение и оптимизацию алгоритмов для
проведения интерактивных вычислительных экспериментов с использованием
пространственной графики современных многопроцессорных компьютеров.
1.1 Список основных обозначений
2. Вычислительный эксперимент в гидромеханике корабля
Современный вычислительный эксперимент позволяет проводить прямое трехмерное
моделирование штормового волнения на море [
        <xref ref-type="bibr" rid="ref2">2</xref>
        ], полноценно воссоздающего гидромеханику
внешних силовых воздействий на поверхности моря, что крайне необходимо для воссоздания
реальных сценариев позиционирования надводного относительно гребней штормовых волн. В
качестве исходных и начальных условий в работе задействуется моделирование трех структур
реального ветрового волнения и зыби, в том числе предельно большой высоты с крутыми
гребнями обрушающихся девятых валов. Виртуальная модель корабля строится в
полномасштабных цифровых теоретических чертежах с охватом надводной общекорабельной
архитектуры, со свойственными массовыми и динамическими характеристиками, включая тягу
условных движителей и рулевых устройств для свободного маневрирования в условиях
интенсивного воздействия штормового волнения. Исходными этапами разработки
предусматривается возможность динамического изменения аппликаты центра тяжести
непосредственно в процессе эксперимента, что необходимо для оценок характера и
интенсивности качки корабля в различных режимах штормового плавания, и может
использоваться для возврата корпуса на киль в случае аварийного опрокидывания.
2.1 Обоснование проектных решений по комплексу программ
      </p>
      <p>
        В проектировании вычислительного эксперимента изначально не было определённости
с выбором математических моделей гидромеханики в условиях интенсивной качки и
возможного опрокидывания корабля на предельно высоком штормовом волнении. Требовалась
поэтапная практическая отработка прямого численного моделирования, и потому для поверки
результатов вычислительных экспериментов были созданы два дополнительных программных
комплекса для изучения особенностей штормового трохоидального волнения «Wave» [
        <xref ref-type="bibr" rid="ref2">2</xref>
        ], и для
реализации всего комплекса традиционных гидростатических расчетов в теории корабля
«Hull»: диаграмм остойчивости и кривых элементов теоретического чертежа, с включением
доступных алгоритмов для расчетов корабельного волнообразования и волнового
сопротивления на различных скоростях хода. Предполагается, что корабельное
волнообразование и воздействие штормовых волн на корпус корабля имеют схожую
физическую природу с одинаковым математическим описанием [
        <xref ref-type="bibr" rid="ref1">1</xref>
        ] в форме интеграла
Мичелла. Программа «Hull» также ориентирована на ведение базы данных цифровых
теоретических чертежей [
        <xref ref-type="bibr" rid="ref3">3</xref>
        ] с возможностью их трансформации, а также для построения иных
вариантов аналитических обводов и формы корпуса гипотетических судов.
      </p>
      <p>Эти интерактивные процедуры позволили отработать геометрические построения и
кинематические алгоритмы для основного вычислительного эксперимента в дружественной
графической среде с использованием традиционных и хорошо поверенных теоретических
методик и математических моделей.
2.2 Особенности реализации прямого вычислительного эксперимента</p>
      <p>Ключевым параметром вычислительного эксперимента является расчетный шаг во
времени, из расчётов с которым строятся геометрические и кинематические аппроксимации в
гидромеханике силового взаимодействия корпуса корабля с предельно интенсивным
штормовым волнением. По аналогии с волновым критерием Куранта в расчётных алгоритмах
контролируется, насколько это возможно, ограничение величин смещений граней числовых
ячеек и малых фрагментов судовой обшивки строго в пределах шага пространственных сеток
для каждого расчётного цикла во времени.</p>
      <p>Другим важным критерием является требование малости трансформации всех
пространственных числовых объектов, внешних источников волнения и, что наиболее
чувствительно, соблюдение плавности изменения форм и устойчивой пространственной
дифференцируемости для сопряженных расчетных ячеек в оконечностях поворачивающегося
корпуса корабля. Для достижения необходимой гладкости пространственных аппроксимаций в
качестве внешнего возмущения задействуются только три групповые структуры с заданными
периодами волн, для которых в вычислительном эксперименте не должно происходить
самопроизвольного повышения частоты, в то время как доплеровские эффекты при отражении
волн от подвижного корпуса корабля могут согласоваться простым уменьшением допустимых
скоростей поворота корабля с установлением меньшего расчётного шага во времени.</p>
      <p>Сформулированные условия проведения вычислительного эксперимента трудно
реализуемы для математических моделей на основе традиционных метацентрических теорий,
также как и в случае интегрирования давлений воды по смоченным поверхностям судовой
обшивки. Собственно такие выводы получались в ходе последовательной реализации
вариантов прямого вычислительного эксперимента по традиционным методикам теории
корабля, в результате чего вычислительные алгоритмы были сведены к исходным
теоретическим предпосылкам с интегралами по объему корпуса и с однократным
интегрированием по спрямленной поверхности действующей ватерлинии для учета влияния
осредненного градиента наклона волновой поверхности.
2.3 О текущей реализации прямого вычислительного эксперимента</p>
      <p>В вычислительном эксперименте представляется вполне естественным использование
только трёх структур морского волнения (рис. 1): активного ветрового; свежей зыби от
недавних ветров под циклоническими вихрями; и пологие или реликтовые волновые
последовательности от далеких штормов или прибрежных волновых резонаторов. В
гидрометархивах хранятся многолетние записи в срочных телеграммах судовых метеопостов из
всех районов Мирового океана о таких волновых структурах.
Рисунок 1. Профили трохоидальных волн с вертикальным смещением среднего уровня, с
учётом действия ветрового напряжения на взволнованной поверхности моря. Сетка отмечает
поверхности равных давлений под гребнями волн, циклоиды – траектории частиц жидкости.
Алгоритм разностной схемы волновых уравнений свободного распространения волн:
   t  
 w  s· s1  s  
 sV  C ·
 
 s  C · t · sV s1V  /
 w  s  
[м / сек] {s  0..ns  1}
[м]
{s 1..ns  1}
,
где: V – скорость и  – местоположение частиц на поверхности воды;  = 2/ – частота на
траекториях трохоидальных волн с периодом  [сек]; t – шаг во времени,  s – условный шаг и
ns – размер сеточной области.</p>
      <p>Движущаяся модель корабля может раздельно взаимодействовать с каждым из
волновых полей, что гидродинамически корректно определяет силовое воздействие волн на
корпус корабля (рис. 2). Результирующее волновое поле получается простым суммированием и
служит лишь для визуализации хода вычислительного эксперимента.
Рисунок 2 Графическое окно с характеристиками корпуса корабля. Справа вверху картушка с
курсом корабля и стрелками характерных длин, высот и направления движения групповых
структур волн. В нижней части экрана графики бортовой, килевой и вертикальной качки.</p>
      <p>На каждом шаге вычислительного эксперимента заново вычисляются геометрические
характеристики и гидродинамические параметры корпуса c учётом кинематики на
взволнованной поверхности моря. На первом цикле вычислений волновое поле не создается, и
результаты расчётов переприсваиваются исходным числовым объектам для фиксации
массовых, инерционных и гидростатических параметров корабля на тихой воде без качки.
Вычислительный эксперимент образуется последовательностью алгоритмов с
 
векторными и тензорными величинами: B  vW·iV – вектор архимедовых сил всплытия;

W  {0,0,Volume} – вертикальный вектор массы (объема) корпуса;
   
V  g·(B W )·t /(Volume·r) – скорость поступательных перемещений корпуса в корабельной
системе координат. После применения коэффициентов демпфирования, привносящих влияние
вязкости и присоединенных масс жидкости, такая вызванная скорость просто суммируется с
   
текущей скоростью хода судна. R  r ·(vB Gravity) – плечо момента архимедовых сил</p>
      <p>Важным элементом вычислительного эксперимента служит ограничение ускорения при
выходе на заданную скорость хода и установка двух режимов циркуляции как: «руль на борт» и
«полборта». Это позволяет моделировать возникновение крена на циркуляции и потерю хода на
штормовом волнении, дельфинирование при захвате корпуса попутной волной и другие
эффекты качки и ходкости по всем шести степеням свободного маневрирования судна.</p>
      <p>Ускорению вычислений служат упорядоченные по шпациям числовые структуры для
описания поверхностей судовой обшивки, что даёт возможность включения реентерабельных
процедур для распараллеливания процесса всего вычислений. На каждом шаге эксперимента
заново интерполируется действующая ватерлиния и ставятся бимсы для динамического
разделения подводного и надводного объемов корпуса. По подводной части вычисляются
измененные объемы, смещённые центры и новые моменты инерции массы корпуса, с
использованием градиентов наклона поверхности волн определяется интегральный объём и
осреднённое направление действия архимедовых сил всплытия. Метацентрические
характеристики корпуса используются только для формального контроля динамики корабля в
рамках традиционных представлений теории корабля, для чего дополнительно вычисляются
площадь, положение центра и моменты инерции площади действующей ватерлинии.
3. Заключение</p>
      <p>Важнейшей особенностью разрабатываемого программного обеспечения является
полноценная трехмерная визуализация всех полей штормового волнения, также как и
пространственного положения и траектории движения корпуса, его гидродинамических
центров и метацентрических параметров мгновенной остойчивости. С использованием
графических средств OpеnGL для изменения прозрачности и выбора вариантов контурных
прорисовок судовых обводов и надводной архитектуры корабля, создаётся особая
экспериментальная среда инженерных изысканий в проектировании перспективного флота
повышенной штормовой мореходности. Новый инструментарий полезен в практическом
судовождении для оценки вариантов штормового маневрирования корабля, для выбора
эффективного курса и скорости хода при оптимальной комфортности и безусловной
безопасности реального штормового плавания.
Литература</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1]
          <string-name>
            <surname>Храмушин</surname>
            <given-names>В. Н.</given-names>
          </string-name>
          <article-title>Поисковые исследования штормовой мореходности корабля. История эволюционного становления корабельного дела, о единении морских инженерных наук и хорошей морской практики</article-title>
          . Lambert Academic Publishing.
          <year>2018</year>
          .
          <volume>470</volume>
          стр.
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2]
          <string-name>
            <surname>Degtyarev</surname>
            <given-names>A</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Gankevich</surname>
            <given-names>I</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kulabukhova</surname>
            <given-names>N</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Khramushin</surname>
            <given-names>V.</given-names>
          </string-name>
          (
          <year>2020</year>
          ).
          <article-title>Computational Model of Unsteady Hydromechanics of Large Amplitude Gerstner Waves</article-title>
          .
          <source>EPJ Web of Conferences. 226</source>
          .
          <year>02009</year>
          .
          <volume>10</volume>
          .1051/epjconf/202022602009
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <surname>Богданов</surname>
            <given-names>А. В.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Храмушин</surname>
            <given-names>В</given-names>
          </string-name>
          . Н. База данных: «Vessel»
          <article-title>- Проектные чертежи для анализа гидростатических характеристик, остойчивости и ходкости корабля</article-title>
          .
          <source>Роспатент. База данных СПбГУ № 2015621368 от 8 сентября 2015 г.</source>
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>