Что сравниваем
Варианты обводов, транца, реданов и подводных элементов — в согласованных состояниях загрузки и при одинаковых заранее заданных ограничениях.
DESIGN & ENGINEERING · CFD
Сравниваем варианты корпуса, оцениваем сопротивление и ходовую посадку, определяем эффективную мощность и по результатам расчёта обоснованно изменяем корпус и движительный комплекс.
ОГЛАВЛЕНИЕ
Статья последовательно разбирает постановку задачи, уравнения, свободную поверхность, подобие режимов, расчётную область, сетку и временной шаг. Далее рассматриваются движения корпуса, интерпретация полей, расчёт мощности, проверка достоверности и состав инженерного отчёта.
Короткий маршрут. Для общего понимания достаточно разделов 01, 05, 06, 10–12 и 17. Для подготовки исходных данных и проверки отчёта дополнительно полезны разделы 07–09, 13, 15, 18 и 19.
ЗАДАЧА РАСЧЁТА
Ценность расчёта — в сопоставлении вариантов и в обоснованном изменении обводов, компоновки или движительного комплекса.
Варианты обводов, транца, реданов и подводных элементов — в согласованных состояниях загрузки и при одинаковых заранее заданных ограничениях.
Сопротивление, давление и трение, волнообразование, просадку, ходовой дифферент, след и изменение характеристик по скорости.
До расчёта определяем целевую величину, допустимые компромиссы и ограничения. По результатам расчёта изменяем только те элементы проекта, влияние которых подтверждено сопоставлением одинаковых режимов.
РЕЗУЛЬТАТ ИССЛЕДОВАНИЯ
Расчёт показывает, где формируется избыточное давление, как развивается волновая система, как меняются смоченная поверхность и ходовая посадка. Поля скорости и давления помогают понять причину результата; интегральные силы и моменты позволяют сравнить варианты количественно.
Выбранная геометрия возвращается в общую 3D-модель. Если меняются масса, центровка, выступающие элементы или режим движения, условия расчёта также уточняются.
СОСТАВ РЕЗУЛЬТАТА
Состав отчётных материалов согласуют до начала расчёта. Отчёт должен содержать данные, необходимые для выбора проектного решения, а не только визуализации полей.
Зафиксированные исходные данные: геометрия, масса, центровка, режимы, физические и граничные условия, способы проверки и границы применимости результата.
Кривые полного сопротивления и, если выбранная методика позволяет, разделение вкладов давления и вязкого трения.
Динамическое погружение, ходовой дифферент, вертикальная сила и продольный момент.
Давление, скорость, волновая система, вихри и след в тех сечениях, которые важны для проектного вывода.
Переход от сопротивления к мощности, необходимой для движения корпуса, и данные для последующего движительного расчёта.
Варианты, сопоставленные по заранее выбранному критерию, сопутствующие потери и ограничения и перечень обоснованных изменений проекта.
ИСХОДНЫЕ ДАННЫЕ
Точный перечень зависит от вопроса. Чем достовернее геометрия, масса, центровка и режимы эксплуатации, тем полезнее выводы.
В КОНТЕКСТЕ ПРОЕКТА
Гидродинамический результат нельзя рассматривать отдельно от массы, центровки, прочности и компоновки. Геометрия с меньшим сопротивлением может забрать нужный внутренний объём, ухудшить производимость или потребовать иной силовой схемы.
Поэтому вариант сначала рассматривают как отдельную гидродинамическую задачу, а затем проверяют в составе всей яхты. Цель расчёта — не оправдать исходные обводы, а определить, какие изменения действительно улучшают проект.

РАЗБОР ПРОЕКТА
Для моторной яхты MLR60 расчёт связывает режим движения, обводы, сопротивление и требуемую мощность. Варианты сопоставляют для одного состояния загрузки и одинаковых ограничений; если разрешены вертикальное перемещение и дифферент, сопоставляют установившуюся ходовую посадку каждого варианта, а не удерживают корпуса на одной искусственно заданной ватерлинии.

ОБЛАСТЬ ПРИМЕНЕНИЯ
Расчёт оправдан там, где существует выбор, риск или вопрос, который нельзя с достаточной надёжностью решить одной эмпирической зависимостью.
Когда нужно сопоставить несколько вариантов корпуса и понять не только разницу в сопротивлении, но и её физическую причину.
Когда корпус работает между водоизмещением и развитым глиссированием, а область применимости простой регрессии становится узкой.
Когда требуется связать сопротивление с эффективной мощностью, следом у винта и характеристиками движительного комплекса.
Когда просадка и дифферент заметно меняют смоченную поверхность, сопротивление, обзор или работу выступающих элементов.
Когда волновые системы корпусов взаимодействуют, а расстояние между ними становится частью гидродинамической задачи.
Когда отрыв, вентиляция, динамическая разгрузка и баланс сил нельзя отделить от общей геометрии и положения центра масс.
Когда изменяются масса, двигатель, винт, кормовая геометрия или подводные элементы существующего судна.
Когда фактическая скорость или посадка отличаются от ожидаемых и нужно разделить влияние корпуса, массы, поверхности и движителя.
ПЕРЕХОД К МЕТОДИКЕ
Если геометрия подготовлена и расчётная задача сформулирована, исходные данные можно сразу передать на предварительную оценку. Далее подробно разобрано, как численная модель позволяет определить сопротивление, мощность и необходимые изменения корпуса.
ИНЖЕНЕРНАЯ МЕТОДИКА CFD
CFD, или Computational Fluid Dynamics, — численная гидрогазодинамика. Корпус, вода и воздух становятся расчётной моделью: задаются геометрия, свойства сред, масса, центровка, скорость, движения корпуса и границы области, после чего численно решаются уравнения сохранения массы и импульса.
Цветная картина давления или волн — лишь один из способов представить решение. Инженер работает с полями скорости, давления, объёмной доли фаз и турбулентных величин, а также с интегральными силами и моментами, сопротивлением, изменением вертикального положения и ходовым дифферентом.
Расчётную задачу формулируют через изменение, которое предстоит проверить: обводы, положение редана, расстояние между корпусами, режим работы крыла, расчётную мощность или параметры движителя. Численно устойчивое решение без заранее заданного критерия сравнения не позволяет сделать обоснованный проектный вывод.
Верификация оценивает численную погрешность решения, а валидация — соответствие физической модели экспериментальным данным.
ЧИСЛЕННЫЙ ЭКСПЕРИМЕНТ
В модели сплошной среды вода и воздух описываются полями скорости, давления, плотности и других величин. Компьютер не хранит их во всех точках непрерывного пространства: расчётная область разбивается на конечное число контрольных объёмов, а дифференциальные уравнения превращаются в связанную систему алгебраических уравнений.
Решатель согласует давление, скорость и доли фаз так, чтобы дискретные балансы выполнялись в расчётной области. Затем напряжения на смоченной поверхности интегрируются, положение границы вода–воздух определяется по полю объёмной доли, а инженер получает силы, моменты и распределённые величины, объясняющие их происхождение.
Визуализацию формируют при обработке результатов. Она полезна только вместе со шкалой, режимом, графиками изменения целевых величин во времени и сведениями о численной модели.

Проверяемость результата. Иллюстрация объясняет устройство модели. Чтобы результат можно было проверить, в отчёте указывают исходные условия, размер области, принцип построения сетки, физические модели, критерии сходимости и оценку неопределённости.
ИСТОРИЯ МЕТОДА
Современные CFD-методы в судовой гидродинамике стали возможны благодаря развитию механики жидкости, способов дискретизации, моделей турбулентности, вычислительной техники и процедур сопоставления с экспериментом.
Вводит в уравнения движения жидкости вязкое трение.
Выводит уравнения движения вязкой несжимаемой жидкости и связывает напряжения со скоростью деформации.
Вводит осреднение турбулентного движения, разделяя мгновенную скорость на среднюю и пульсационную составляющие.
Публикуют метод маркеров и ячеек (MAC) для нестационарного течения со свободной поверхностью.
Предлагают метод VOF для переноса объёмной доли жидкости и определения межфазной границы.
Формулирует модель турбулентности k–ω SST.
Публикует обновлённые рекомендации по судовой CFD, расчёту сопротивления, самодвижению и оценке неопределённости.
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ
Физическую постановку определяют геометрия, граничные и начальные условия, свойства воды и воздуха, модель турбулентности, способ описания свободной поверхности и разрешённые движения корпуса. Формулы ниже показывают общий смысл двухфазной несжимаемой модели; решатель использует согласованную дискретную форму этих уравнений.
u — вектор скорости потока.
Для несжимаемых фаз дивергенция поля скорости равна нулю: в контрольном объёме нет источника или стока жидкости. Плотность двухфазной смеси изменяется через объёмную долю фаз.
ρ — плотность; u — скорость; p — давление; τeff — эффективные напряжения; g — ускорение свободного падения; fσ — сила поверхностного натяжения.
В консервативной записи учитываются инерция, давление, эффективные вязкие напряжения, тяжесть и, когда это существенно, поверхностное натяжение. Плотность ρ и вязкость зависят от объёмной доли воды.
Для характерных судовых чисел Рейнольдса течение преимущественно рассматривают как турбулентное. Прямое воспроизведение всего спектра турбулентных масштабов потребовало бы сетки и временного шага, непрактичных для повседневного расчёта полноразмерной яхты. При этом ламинарно-турбулентный переход нельзя автоматически исключать из модельных испытаний, малых аппаратов и отдельных обтекаемых элементов.
В разложении Рейнольдса мгновенная скорость представляется суммой средней и пульсационной составляющих: uᵢ = ūᵢ + u′ᵢ. После осреднения появляются рейнольдсовы напряжения; их влияние на среднее течение описывает модель турбулентности.
Модель k–ω SST Ментера — распространённый инженерный выбор, но не универсальный ответ. Пограничный слой без отрыва, развитый отрыв, вентиляция редана и вихри подводного крыла требуют разных проверок чувствительности и иногда разной детализации турбулентной структуры.
Физический смысл уравнения импульса. Левая часть описывает накопление и перенос количества движения. Правая — действие давления, вязких и турбулентных напряжений, тяжести и межфазных сил. Расчётная программа не «подставляет корпус в формулу», а обеспечивает этот баланс во всех контрольных объёмах и на каждом расчётном шаге времени.

ДИСКРЕТИЗАЦИЯ
Большинство прикладных судовых RANS-расчётов выполняется методом конечных объёмов — FVM. Интегральная форма уравнений записывается для каждой ячейки: потоки массы и импульса через её грани, объёмные источники и накопление должны образовать согласованный баланс.
Значения на гранях получают интерполяцией, поэтому порядок и ограничение численных схем влияют на диффузию, устойчивость и резкость межфазной границы. Соседние ячейки связаны общими потоками; давление и скорость согласуются в общей итерационной процедуре. Локальное соблюдение законов сохранения — важное достоинство FVM, но само по себе оно не гарантирует точность.
В типовой судовой RANS-CFD пространственную дискретизацию обычно выполняют методом конечных объёмов. Метод конечных элементов использует другую дискретную постановку и широко применяется в прочностных расчётах, а также в отдельных задачах механики жидкости.
Способ пространственной дискретизации уравнений сохранения.
Другая численная постановка; особенно привычна в расчётах прочности.
Способ отслеживать границу вода–воздух внутри расчётных ячеек.
ВОДА И ВОЗДУХ
Носовая волна, освобождение транца и интерференция волн не задаются заранее. Они возникают из решения двухфазной задачи для воды и воздуха при выбранной кинематике корпуса. Брызги, тонкие плёнки и захват воздуха требуют существенно более высокого пространственного и временного разрешения, чем общая волновая система.
В модели VOF каждой ячейке назначается объёмная доля воды α.
α — объёмная доля воды; u — вектор скорости смеси.
В практическом решателе применяют ограничивающие схемы, удерживающие α в пределах от 0 до 1, и методы сжатия межфазной границы. Плотность и вязкость смеси определяются через α, а переход 0 < α < 1 занимает конечную полосу ячеек.
Если сетка у поверхности слишком груба, возникает избыточная численная диффузия или выбран слишком большой временной шаг, межфазная граница размывается, а высота волн и характер отрыва потока от транца искажаются. Сетку вблизи свободной поверхности отдельно сгущают, а результат проверяют по сопротивлению, балансу массы и форме волны.

Модель турбулентности описывает перенос импульса пульсациями, которые сетка не воспроизводит непосредственно. Она не моделирует автоматически поступление воздуха со свободной поверхности и не описывает образование пара при падении давления.
Поэтому выражение «двухфазный расчёт» ещё не определяет физику задачи. Нужно назвать сами фазы, механизм их появления и величину, по которой будет проверяться результат.
Определяет структуру пограничного слоя, отрыв и смешение в рамках выбранного турбулентного замыкания.
Возникает, когда воздух получает устойчивый путь от атмосферы к редану, винту, крылу или другой зоне разрежения. Сетка и межфазная модель должны воспроизводить этот канал поступления воздуха.
Связана с падением абсолютного давления до уровня, при котором в воде образуется пар. Нужны давление насыщенного пара и специальная кавитационная модель; обычного VOF вода–воздух недостаточно.
ПОДОБИЕ РЕЖИМОВ
Скорость сама по себе не определяет режим. Для безразмерного сравнения заранее выбирают характерную длину L и последовательно используют её во всех вариантах; для судна это часто длина между перпендикулярами или по ватерлинии, но для специальных задач применяются и другие определения.
V — скорость; g — ускорение свободного падения; L — характерная длина.
Число Фруда сопоставляет инерционные и гравитационные эффекты. Оно особенно важно для волнового подобия и ходовой посадки.
V — скорость; L — характерная длина; ν — кинематическая вязкость.
Число Рейнольдса сопоставляет инерционные и вязкие эффекты. Вязкость зависит от среды и температуры, поэтому её значение фиксируют в исходных данных.
Модельные испытания сохраняют подобие по Фруду. В модельном эксперименте невозможно одновременно сохранить числа Фруда и Рейнольдса. Обычно воспроизводят волновое подобие по Фруду, а вязкую составляющую пересчитывают. Расчёт в натурном масштабе исключает часть масштабных преобразований, но повышает требования к описанию шероховатости, пристеночной области и натурного числа Рейнольдса.
Границы перекрываются и зависят от формы, загрузки и выбранного определения длины. Для глиссирующих корпусов дополнительно применяют объёмное или ширинное число Фруда. Ни одна универсальная граница Fr сама по себе не определяет тип движения.
РАСЧЁТНАЯ ОБЛАСТЬ
Корпус помещается в область со входной, выходной, боковыми, верхней и нижней границами. Их удаляют настолько, чтобы искусственное стеснение, возмущение давления и отражённые волны не меняли искомую величину. При ограниченной глубине дно, напротив, становится физической границей задачи.
На входе задают скорость, доли фаз и параметры турбулентности; на выходе — условия, обеспечивающие свободный выход возмущений без заметного отражения. Начальное давление должно соответствовать гидростатике, а для свободной поверхности применяют зоны гашения или релаксации волн. На корпусе действует условие вязкой непроницаемой стенки.
Половина области допустима лишь для физически симметричного прямолинейного режима без несимметричной геометрии, винтового закручивания или боковых воздействий. Размер области и влияние граничных условий проверяется по чувствительности целевого результата.
ПРОСТРАНСТВЕННОЕ РАЗРЕШЕНИЕ
Качество сетки оценивают по тому, насколько подробно она воспроизводит физически важные изменения потока. Ячейки сгущают у корпуса, свободной поверхности и в других областях с большими градиентами; там, где решение меняется медленно, сетку делают реже.
Призматические слои следуют поверхности и передают изменение скорости от стенки к внешнему потоку.
Сгущение удерживает форму волны, носовую систему, брызговую пелену и освобождение транца.
Острые кромки, зоны отрыва, концевые вихри и вентиляция требуют локального сгущения сетки.
Кормовую часть расчётной области продолжают достаточно далеко, чтобы получить поле скоростей перед винтом и развитие вихревого следа за корпусом.
Безразмерная координата связывает расстояние от стенки до центра первой ячейки, касательное напряжение и вязкость.
y — расстояние до стенки; uτ — динамическая скорость; τw — касательное напряжение на стенке; ρ — плотность; ν — кинематическая вязкость.
Согласно актуальной рекомендации ITTC, при непосредственном разрешении вязкого подслоя обычно стремятся к y+ ≤ 1. Классические пристеночные функции требуют, чтобы первая точка находилась в логарифмическом слое, ориентировочно 30 < y+ < 100.
Выбранный диапазон должен соответствовать конкретной модели и её пристеночной обработке. Смешение разрешённого подслоя и области пристеночных функций на одном корпусе может сделать поле трения несогласованным. Кроме y+, проверяются число слоёв, коэффициент их роста, сохранность слоёв у скул и транца, неортогональность и скачки размера ячеек.
Проверка сеточной сходимости. Если сопротивление заметно меняется при систематическом сгущении сетки, результат ещё зависит от численного разрешения. Разница между двумя красивыми картинками может быть следствием различий в параметрах расчёта, а не геометрии корпуса.
РАСЧЁТ ВО ВРЕМЕНИ
Расчёт свободной поверхности, движения корпуса, вращающегося винта и развитого отрыва обычно ведётся во времени. На каждом физическом шаге решатель согласует связанные поля до принятой внутренней сходимости, однако малые невязки внутри шага ещё не означают выхода расчёта на установившийся или статистически установившийся режим.
u — скорость потока; Δt — шаг по времени; Δx — характерный размер ячейки.
Число Куранта показывает, какую долю характерного размера ячейки проходит поток за один шаг. Для явного алгоритма условие Куранта ограничивает устойчивость; неявный алгоритм может оставаться устойчивым при большем Co, но это не гарантирует временную точность. В VOF дополнительно контролируют число Куранта в зоне межфазной границы.
Шаг выбирают по наиболее быстрому существенному процессу. Для периодического явления ITTC рекомендует как практический ориентир не менее ста шагов на период, после чего проверяется чувствительность результата к Δt. Усредняют только статистически установившийся участок графиков изменения сил и движений.
НОРМИРОВАННЫЙ СХЕМАТИЧЕСКИЙ ПРИМЕР
Сетка и физические модели во всех трёх запусках одинаковы. Сопоставляют среднее значение, амплитуду и фазу колебаний целевой величины; одного внешнего сходства кривых недостаточно.
На узком экране график можно прокрутить по горизонтали.
ДВИЖЕНИЕ КОРПУСА
У быстроходной яхты гидродинамическое давление создаёт сопротивление, вертикальную силу и продольный момент. Поэтому изменение вертикального положения корпуса и дифферента часто решают совместно с потоком — обычно в двух степенях свободы, если остальные движения не входят в поставленный вопрос.

Фиксация корпуса на проектной ватерлинии может дать устойчивое поле, но искажённое значение сопротивления. Особенно чувствительны переходные и глиссирующие корпуса, катамараны, суда с реданами и подводными крыльями.
Свободная посадка не означает произвольные шесть степеней свободы. Выбираются только физически значимые движения; задаются масса, координаты центра масс и необходимые моменты инерции. Переход к квазистационарному или периодическому режиму контролируется по силам, моментам, перемещениям и их средним значениям.
Полученную посадку оценивают по требованиям проекта: обзору, клиренсу, погружению транца, работе движителя и соответствию расчётной загрузке.
Состояние загрузки — не только водоизмещение. Лёгкое судно, полный запас топлива и переходный эксплуатационный вариант могут иметь одинаковую или близкую массу, но разную центровку и посадку. Для активных степеней свободы фиксируются необходимые координаты центра масс и инерционные характеристики; для частично заполненных танков учитывается свободная поверхность жидкости. Результат одного состояния нельзя без проверки линейно переносить на другое.
ИНТЕРПРЕТАЦИЯ
Инженерное заключение строят по совокупности результатов. Сопротивление без картины потока не показывает, какую часть геометрии следует менять; поле давления без интегральных сил не показывает величину эффекта; малые невязки без устойчивых графиков изменения сил не подтверждают завершение расчёта.
Коэффициент давления Cp показывает распределение относительно выбранного опорного давления. Вклад отдельного участка в сопротивление или подъёмную силу определяют с учётом его площади и направления нормали.
Модуль и направление пристеночного напряжения помогают находить области малого трения и возвратного течения. Итоговый вязкий вклад определяется интегрированием по всей принятой смоченной поверхности.
Сопоставляются продольные и поперечные сечения волн, фаза интерференции и освобождение транца. Один высокий гребень не равен автоматически большему волновому сопротивлению.
Поле оценивается в зафиксированных сечениях и системе координат. Для движителя важны не только средняя скорость, но и неоднородность, закрутка и выбранный способ осреднения.
Критерии Q или λ₂ зависят от порога визуализации. Их сопоставляют с полями скорости и давления, а не используют как самостоятельное доказательство потерь.
В отчёте задаётся интервал усреднения и показывается разброс сопротивления, сил и движений. Отдельное мгновенное значение не характеризует установившийся или периодический режим.
Проверка причин различий. Распределения давления и скорости позволяют сформулировать гипотезу, которую затем проверяют контролируемым изменением одной группы параметров и сопоставлением интегральных величин. Если одновременно изменились форма, водоизмещение и центровка, вклад каждого фактора из одной пары расчётов определить нельзя.
ДВИЖИТЕЛЬНЫЙ БАЛАНС
Расчёт буксируемого корпуса даёт сопротивление при заданной скорости и посадке. Для выбора двигателя дополнительно учитывают поле притока, взаимодействие корпуса и движителя, характеристики винта или водомёта, потери в редукторе и валопроводе и эксплуатационный запас. Эти составляющие нельзя заменить одним произвольным процентом.
RT — полное сопротивление; ρ — плотность воды; S — принятая площадь смоченной поверхности; V — скорость.
Безразмерная запись требует явно определить площадь S и одинаково вычислять её для сравниваемых режимов и вариантов.
PE — эффективная мощность; RT — полное сопротивление; V — скорость.
Мощность, непосредственно затрачиваемая на движение корпуса с заданной скоростью.
PD — мощность, переданная движителю; ηD — общий пропульсивный коэффициент.
Общий пропульсивный коэффициент связывает эффективную мощность корпуса с мощностью, переданной движителю.
Переход от PD к тормозной мощности PB дополнительно учитывает механический КПД передачи. Для прогноза скорости нужны характеристики двигателя, редуктора, винта или водомёта, диаметр, обороты, кавитационные ограничения и поле притока. CFD-расчёт самодвижения решает иной, более полный вопрос, чем расчёт сопротивления голого корпуса.
РАСЧЁТНАЯ ПРАКТИКА ASTRATELUM
Для всех вариантов сохраняют одну расчётную задачу, физические условия, модели и критерии проверки. Сетки для разных вариантов геометрии не обязаны совпадать буквально, но должны обеспечивать сопоставимую детализацию ключевых зон, одинаковый способ описания пристеночной области и оценённую численную неопределённость. Иначе различия могут быть вызваны дискретизацией, а не формой корпуса.
Проверяем замкнутость, устраняем нерасчётные детали, фиксируем ватерлинию и выступающие элементы.
Для каждого варианта задаём состояние загрузки, координаты центра масс, необходимые инерционные данные и диапазон скоростей. До запуска определяем целевую величину, ограничения и допустимый компромисс.
Определяем масштаб, турбулентность, свободную поверхность, разрешённые движения и способ представления движителя.
Проверяем область, локальные сгущения, пристеночные слои, временной шаг и зависимость выбранной целевой величины от плотности сетки.
Считаем согласованные режимы, контролируем невязки и балансы, графики изменения сил и моментов, посадку, свободную поверхность и достижение статистически установившегося режима.
Анализируем поля давления, трения, волн и следа и формулируем обоснованное изменение геометрии, которое можно проверить повторно.
Повторяем варианты в сопоставимой постановке, сопоставляем разницу с численной неопределённостью и фиксируем вывод вместе с областью его применимости.
Условия сопоставимости вариантов. Если изменение обводов уменьшило сопротивление одновременно с водоизмещением, полезным объёмом или смоченной поверхностью, это уже другой проектный компромисс. Перед серией фиксируют, что остаётся неизменным — водоизмещение, центровка, объём, габариты, посадка или мощность, — а изменения остальных величин показывают в отчёте вместе с гидродинамическим выигрышем.
ВЫЧИСЛИТЕЛЬНЫЕ ЗАТРАТЫ
В типовых задачах Astratelum окончательная сетка корпуса со свободной поверхностью нередко содержит от нескольких миллионов контрольных объёмов; ориентир около пяти миллионов характеризует обычный для нас объём вычислений, а не норматив качества. Итоговый размер задают геометрия, режим, требуемая детализация пограничного слоя, волн, следа и локальных отрывных явлений.
Рабочий ориентир, который уточняется по чувствительности целевой величины.
Время непрерывного расчёта одного режима — от трёх суток. Расчёт сложного окончательного режима может потребовать больше времени; подготовка геометрии, построение сетки и анализ результата в этот срок не входят.
Скорости, геометрические варианты и проверочные сетки считаются отдельно.
Время зависит не только от количества ячеек. На него влияют временной шаг, число внутренних итераций, модель турбулентности, VOF, разрешённые движения корпуса, движитель, качество сетки и требуемая статистическая устойчивость.
Подготовка геометрии, построение сетки и интерпретация результата выполняются отдельно от непрерывного вычисления режима. Серия вариантов может занимать несколько недель расчётного времени, особенно если каждый вывод проверяется сгущением сетки.
ПРОВЕРКА ДОСТОВЕРНОСТИ
Установившиеся силы и малые невязки показывают итерационную сходимость конкретного запуска. Верификация оценивает численную ошибку при решении выбранной математической задачи. Валидация показывает, насколько математическая модель воспроизводит эксперимент; при сравнении учитывают неопределённость расчёта и измерений.
Исследуются систематически изменяемые сетки и временные шаги, а также чувствительность к области, схемам и итерационным допускам. Результатом является оценка численной неопределённости выбранной величины.
Расчёт сравнивают с результатами соответствующего модельного эксперимента или натурными измерениями. Расчёт контрольного корпуса проверяет процедуру, но не заменяет валидацию новой яхты во всём диапазоне режимов.
Совместно рассматриваются численная неопределённость, погрешность измерений и неопределённость исходных параметров. Разность расчётного и измеренного значений без этих оценок не подтверждает достоверность модели.
НОРМИРОВАННЫЙ СХЕМАТИЧЕСКИЙ ПРИМЕР
Геометрия, режим и численные модели остаются неизменными. Сетки строят по одному принципу и последовательно уменьшают характерный размер ячейки.
На узком экране график можно прокрутить по горизонтали.
Три близких числа сами по себе ещё не гарантируют монотонной сходимости. Проверяются коэффициент сгущения, порядок метода, поведение разностей и применимость выбранной оценки неопределённости.
НЕЗАВИСИМАЯ ИНЖЕНЕРНАЯ ОЦЕНКА
Перед серией вариантов исходный корпус полезно проверить ещё одним способом. Для обычного водоизмещающего или умеренно быстроходного однокорпусного судна применяют метод Холтропа — Меннена; для установившегося глиссирования призматического корпуса — метод Савицкого. Оба расчёта выполняют для одной геометрии, массы, центровки и скорости, а сравниваемые величины приводят к одинаковым единицам и определениям.
Сравнивают форму кривой полного сопротивления и эффективной мощности. Метод используют только в опубликованном диапазоне исходных коэффициентов. Для многокорпусных судов, подводных крыльев и необычных обводов такая оценка неприменима либо требует отдельно подтверждённых поправок.
Исходная работаСопоставляют сопротивление, ходовой дифферент, динамическую подъёмную силу и смоченную длину. Базовая модель относится к установившемуся движению призматического глиссирующего корпуса; реданы, переменная килеватость, интерцепторы и крылья требуют самостоятельных моделей или обоснованных поправок.
Исходная работаЕсли разница превышает ожидаемую неопределённость или меняет проектный вывод, проверяют исходные данные, единицы измерения, область применимости формулы, сетку, граничные условия, режим и степени свободы корпуса. Ошибочную постановку исправляют и расчёт повторяют; обоснованное физическое различие объясняют и фиксируют в отчёте.
Классическая оценка не заменяет сопоставление с экспериментом или натурными данными. Совпадение с методом Холтропа — Меннена или Савицкого повышает доверие к постановке, но само по себе не доказывает точность CFD для нового корпуса. Расхождение также не означает автоматически ошибку: сначала проверяют, описывают ли методы одну физическую задачу и находятся ли исходные параметры в пределах их применимости.
СОПОСТАВЛЕНИЕ С НАТУРОЙ
Для установившейся скорости Astratelum использует модуль относительного отклонения расчётного значения от результата ходовых испытаний. Целевой критерий — не более 10% в подтверждённой области применения методики.
Критерий применяется при достоверной построечной геометрии, известной массе и центровке, согласованном состоянии поверхности, корректных характеристиках движительного комплекса и сопоставимых условиях испытаний.
Он относится к оговорённому режиму и не является универсальной гарантией для любой будущей модели. В отчёте вместе с абсолютным отклонением фиксируют знак ошибки, диапазон скоростей и неопределённость натурных данных.
Vcalc — расчётная скорость; Vtrial — скорость, измеренная на ходовых испытаниях.
Проектная цель: не более 10%Критерий применяют только к оговорённому режиму; неопределённость расчёта и испытаний оценивают отдельно.
Область применимости критерия. Опубликованные валидационные исследования показывают расхождения в несколько процентов для отдельных судов и режимов. Эти значения относятся к конкретной геометрии, эксперименту и расчётной постановке и не переносятся автоматически на другую модель.
ГРАНИЦЫ МЕТОДА
Численный расчёт не компенсирует неизвестную массу, неточную построечную геометрию, неверные данные движителя или отсутствующую проверку чувствительности. Вывод действует только для исследованных состояний, физических моделей и критериев.
Одна расчётная скорость может отвечать на узкий вопрос, но не показывает форму кривой сопротивления, переход между режимами и запас движителя.
Для всех вариантов нужны сопоставимые границы расчётной области, физические модели, детализация ключевых зон и критерии сходимости.
Если корпус должен свободно изменять погружение и дифферент, искусственно заданная ватерлиния способна исказить силы и картину потока.
Эффективная мощность корпуса не учитывает движительный КПД и потери установки.
Число ячеек без проверки сеточной сходимости и контроля y+ не подтверждает результат.
Картинка не заменяет балансов, историй сил, исследования сетки и времени, а при заявленной валидации — экспериментальных данных.
ИНЖЕНЕРНЫЙ ОТЧЁТ
Полная воспроизводимость расчёта требует исходных файлов, версии программы и настроек решателя. Инженерный отчёт должен позволять проследить, для какой геометрии, состояния загрузки, физических моделей, сетки, критериев и области применимости получен каждый вывод.
Каждый расчётный случай получает идентификатор, связанный с именем и версией геометрии, состоянием загрузки и набором настроек. Отдельно фиксируются единицы, начало и направление осей, точка приведения моментов, включённые и исключённые детали, версия расчётной программы и дата запуска.
Поля снабжаются шкалой, единицами, системой координат и положением сечения; графики — интервалом и способом усреднения. Таблица случаев показывает, какие параметры были неизменными, какие менялись и какая целевая величина относилась к каждому сравнению.
В заключении отчёта указывают, что непосредственно вычислено, как поля и балансы объясняют результат, какие изменения геометрии обоснованы расчётом и какую её версию следует передать на следующую стадию проектирования.
Структура инженерного вывода. В отчёте полезно различать: «получено в расчёте» — численное значение или поле; «наблюдается при сравнении» — воспроизводимое различие вариантов; «инженерный вывод» — объяснение и рекомендуемое действие. Такое разделение отделяет непосредственно рассчитанные данные от их интерпретации и проектной рекомендации.
ПЕРЕД ЗАПУСКОМ
Расчёт запускают после фиксации следующих данных и критериев.
КРАТКИЙ СЛОВАРЬ
ИСТОЧНИКИ
Источники задают принятую терминологию, процедуры проверки и границы применимости приведённых утверждений.
Постановка судовой CFD, область, сетка, турбулентность, свободная поверхность, сходимость и представление результата.
Открыть документ ITTCРасчёт сопротивления, эффективной мощности, модели и натуры.
Открыть документ ITTCОтдельно оценивают неопределённость численного решения, модели и эксперимента.
Открыть документ ITTCСостав внутренней методики: версия программы, геометрия, состояния, степени свободы, модели, сетка, критерии и определение выходных величин.
Открыть документ ITTCПрактические примеры оценки сеточной неопределённости для судовых задач.
Открыть пример ITTCМоделирование движителя и равновесия сопротивления с тягой.
Открыть документ ITTCОфициальное описание дискретизации методом конечных объёмов.
Открыть документациюКлассическая работа по VOF с объёмной долей жидкости в ячейке.
DOIОдна из основ нестационарного численного расчёта со свободной поверхностью.
DOIРабота, в которой сформулирована модель SST.
DOIКонтрольные случаи и определения для моделей турбулентности.
Открыть ресурс NASAОценка численной неопределённости по систематическому сгущению сетки.
DOIСистемный подход к V&V судовых расчётов.
DOIИсходная работа, лежащая в основе осреднения Рейнольдса.
DOIРегрессионный метод оценки сопротивления и мощности однокорпусных судов в опубликованной области применимости.
DOIКлассическая инженерная модель установившегося глиссирования и её область применимости.
DOIПример сопоставления полноразмерной CFD-модели с ходовыми испытаниями.
DOIПРАКТИЧЕСКИЕ ВОПРОСЫ
ПРЕДВАРИТЕЛЬНАЯ ПОСТАНОВКА ЗАДАЧИ
Опишите, что требуется сравнить или проверить, какие режимы движения важны и какая геометрия уже подготовлена. Определим достаточность исходных данных и состав расчёта, который даст вывод, пригодный для принятия проектного решения.
Николай Андреевич Малиновский — инженер и основатель Astratelum. До начала работ согласуем состав исходных данных, этапы расчёта и порядок взаимодействия.