Хоча робота Туліна ігнорується емпіриками та теоретиками, на її основі побудований метод Воруса для вивчення теорії для вертикального руху симетричних, двовимірних клинів з кривими та прямим сторонами. Головна відмінність роботи Воруса від роботи Туліна полягала у врахуванні точки наведеної поперечної нормалі швидкості на корпус і заміни сингулярної поведінки бризка складною процедурою розкладання. Ворус таким чином отримав інтегральне рівняння другого порядку. Модель Воруса ускладнена. [13]
Багато інших вчених займалися задачами глісування при наявності випуклості. Наприклад - Маріо (1951) проаналізував глісування в довільних числах Фруда. Камбербач (1958) також вивів формули для двовимірних пластин при великих, але кінцевих числах Фруда. Задача усталеного в’язкого опору, який встановлений при постійній довжині хорди була вирішена Ву (1972). В 1967 році Маріо вирішив дану задачу при врахуванні гравітаційності. Він також враховував бризковий опір, але знову оцінював тільки поверхневі шари. [12]
3 Комп’ютерні методи визначення гідродинамічних характеристик глісуючого комплексу
Передбачення створених хвилею рухів і хвильових навантажень - одне з найважливіших питань при конструюванні судна. Рухи з більшою амплітудою створюють задачу про безпечне пересування суден в воді, у той час екстремальне навантаження може привести до пошкодження структури. Загальне застосування методу малих збурень - один підхід до такої нелінійної задачі, де нелінійні ефекти обчислюються за допомогою збереження квадратних позначень у граничних умовах. Однак, в цьому підході залишаються лінійні припущення.
Загальні методи для передбачення характеристик глісуючго корпусу включають емпіричні рівняння і дослідне випробування. Емпіричні рівняння часто можна застосувати тільки до подібних типів корпуса в малому діапазоні параметрів, у той час як випробування моделі часто дуже дорого, особливо для малого судна.
Зараз зростають вимоги до розвитку методу розрахунку, який орієнтується, в принципі, на широкі можливості. Сучасний розвиток в комп’ютерних характеристиках і чисельних методах дозволили вирішити нелінійні задачі набагато легше, ніж раніше.
Об’єднаний метод Ейлер Лагранжа (МЕЛ) вперше був введений Лонгетом-Хігінсом і Скелетом (1976), для моделювання деформації поверхневих хвиль. Метод МЕЛ використовує підхід повністю нелінійної області часу і застосовується до різноманітних нелінійних задач. Метод моделювання для руху пливучого корпуса у хвилях був розроблений Вінжі і Бревігом (1981), Квінті та ін. (1990), Сеном (1993) і Танізава (1995). Двовимірні взаємодії пливучих тіл з вільною поверхнею можна обчислити раціонально, використовуючи повністю нелінійний підхід. З іншого боку, підхід нелінійної часової області був поширений на тривимірну задачу хвилі судна і вивчений багатьма дослідниками. В принципі, обчислені результати повністю нелінійного підходу були отримані Маскевом (1992), Беком та ін. (1994), Скорпіон та ін. (1996) і Шіракура та ін. (2000). Хоча їх формулювання теоретично точні, чисельно стійкі рішення не можна отримати в деяких випадках обчислення.
Для практичного використання потрібно більше досліджень і числових вимірювань. Числовий аналіз тривимірним методом - інший підхід до нелінійної задачі, який можна розглянути як метод малих збурень. Тобто, лінійна або слабо нелінійна умова вільної поверхні для невстановленої області хвилі представлені у більшості випадків, у той час як миттєва геометрія корпусу враховується в обчисленні в умові поверхні корпусу. Оскільки з цими припущеннями очікуються більш стійкі рішення, були представленні успішні результати обчислення рухів судна (наприклад, Лін і Юу 1990, Накос і ін. 1993, Буннік і Германс 1998, Колагросі і ін. 1999, Ясукава 2000 і Катаока і ін. 2001).
Теорія високошвидкісної смуги (ТВШС), вперше представлена Чапманом (1976), застосовувалась багатьма дослідниками, наприклад, Адачі і Охмасу (1980), Енг і Кім (1981), Охмасу і Фалтінсен (1990), Фалтінсен і Жао (l991). Цей метод часто називають "2.3 D " або "2D+T" теорія, у якій задача тривимірної вільної поверхні корабельної хвилі зведена до двовимірної задачі, яка може бути вирішена послідовно в часовій області. Адачі і Масуа (1996) запропонували метод функції Гріна в 2.5D теорії, де було враховано додаткове позначення, яке відповідає поперечним хвилям у задачі усталеного хвилеутворення. Кашігаві (1995) розробив розширену об’єднану теорію, у якій немає обмеження на порядок поступової швидкості чи частоти коливання. Хоча ці два методи, використовуючи двовимірний підхід - дуже практичні з раціональних, теоретичних і обчислювальних точок зору, обидва методи сумісні з лінійними припущеннями.
З’єднання вищезгаданих двох методів приводить до підходу нелінійної часової області, використовуючи ТВШС, яка є нелінійною версією 2.5D теорії. Взаємодія в низу за течією у тривимірному потоці навколо судна пояснюється ефектом запам’ятовування вільної поверхні. Оскільки задачі граничних умов (крайові задачі) можуть бути описані тим же самим формулюванням у випадках і усталеного і неусталеного потоку навколо судна, потенціали швидкості можна вважати однаковими. Крім того, граничні умови повністю нелінійні в структурі теорії тонкого судна, і в обчисленні можна врахувати геометричні, гідростатичні і гідродинамічні нелінійні характеристики. Калісал і Чан (1989) та Тулін і Ву (1996) розробили чисельні моделювання розбіжних головних хвиль. Фонтайн і Квінт (1997) також показали обчислення головних хвиль і запропонували можливість його застосування для прогнозування удару навантаження. Маруо і Сонг (1994) продемонстрували, що головні хвилі розбиваються при русі високошвидкісного судна, яке застосовувалося для аналізу змочення палуби. Обчислення гідродинамічної сили в задачі усталеного коливання були представлені Кіхара і Найто (1998). Крім того, вони досліджували прогнозування додаткового опору моделі Віглея в регулярних головних хвилях. З подальшим розвитком в цій області, були активно вивчені 2D+T методи разом з процедурою обчислювальної гідродинаміки (CFD). Тулін і Ландріні (2000) представили аналіз розбивання хвиль, використовуючи згладжену частину гідродинаміки (SPH), і Андрілон та Алесандріні (2002) показали результати обчислення, використовуючи обчислювальний пристрій Навье-Стокса з об’ємно-кінцевим (VOF) методом. Ці методи дозволяють моделювати гідродинамічні рухи, включаючи комплексні фази розсіяння хвиль, типу повторного сплеску, формування сплеску вверх і завихреності. Для практичного використання в технічних задач, застосування нових CFD процедур, описаних вище для обчислення гідродинамічної сили, є перспективною задачею. [14]
В роботі [15] представлено результати вивчення використання CFD для оцінки характеристик високошвидкісного глісуючого судна, яке рухається зі сталою швидкість по спокійній воді. Для вивчення використовується неструктурований, багатофазний, кінцевий об’ємний код, який використовує метод об’єму рідини (VOF). Характеристика високошвидкісного судна глибоко пов’язана з орієнтацією корпуса у швидкості, що не може бути відома апріорно. Змінюється підйом глісуючих корпусів і кути атаки, як реакція на область тиску, створену потоком. Для врахування цих змін у положенні корпуса метод моделювання повинен гарантувати, що в підйомі була досягнута динамічна рівновага в момент обробки. Це досягнуто за допомогою ітераційної схеми, у якій область потоку була вирішена для дискретних орієнтацій корпуса. Робота складається з набору експериментальних випробувань моделі, для отримання даних, з якими чисельні результати порівнюються. Для цього було виконано три набори моделювань. Перший набір виконувався для прямого порівняння числових результатів з експериментальними. Другий набір моделювань виконаний для задоволеної умови рівноваги. Третій набір - відповідає стану рівноваги при підйомі та обробці. [15]
4 Основні гідродинамічні характеристики
4.1 Додатковий опір
Як зазначено вище, для конструювання суден необхідно вивчення гідродинамічних характеристик та передбачення поведінки судна в умовах природного хвилювання. Одним з параметрів, які необхідно обов’язково враховувати при конструюванні судна - є додатковий опір хвилі. В [14] увага зосереджується на впливі надводної форми носу на додатковий опорі. В роботі до задачі неусталеного хвилеутворення застосовується нелінійний 2D+T метод, в якому потік описується, використовуючи суперпозицію набігаючої хвилі з іншим збуреннями, які викликані корпусом. Не зважаючи на те, що включені процедури апроксимації, цей метод є практичним для дослідження нелінійних характеристик гідродинамічної сили. Для збільшення енерго-економічних суден, потрібно більше зменшення опору хвиль. Частково, це через те, що відношення додаткового опору хвиль до загального опору стає більшим. Хоча покращення форми корпусу дозволяє зменшити кінську силу для руху в стоячій воді, кінська сила, необхідна для руху в незмінних хвилях. Протилежна до цієї задачі, задача покращення форми корпуса через надводну геометрію коректна для конструювання судна, оскільки характеристика хвилеутворення, яка залежить від підводної форми корпуса, може зберігатися. Однак, використовуючи існуючі методи розрахунку, важко правильно розрізнити різницю між варіантами надводних форм корпуса. Зокрема, нелінійні ефекти динамічної сили, які виникають через розсіювання і дифракцію хвиль ще потрібно вивчити. Поступово, зрозуміло, що надводна форма носу впливає на усталені сили хвиль, тобто на додатковий опір, а не на рух хвилеутворення. Наприклад, вплив надводної форми носу на додатковий опір було експериментально досліджено Найто і ін. (1996). Вони прийшли до висновку, що тупоносі судна з різними надводними формами носу показали різні значення додаткового опору. Тому, розробка обчислювального інструмента, який дасть можливість проектувальникам судна обговорювати переваги надводних форм носу є важливим завданням. [14]
У багатьох застосування в морських умовах додатковий опір відіграє важливу роль. Але багато з існуючих методів недооцінюють додатковий опір при низьких частотах довжини хвилі. Відомо, що для опису руху плавання суда у хвилях дуже гарні результати для багатьох практичних форм корпусу показує теорія стрічки. В останні роки програми обчислювальних машин розвинулись таким чином, що можуть обчислити сили і рухи плавання судна у хвилях за допомогою лінійних дифракційних програм. Фактично, метод використовує лінеаризацію навколо незбуреного потоку навколо судна, що може привести до гарної апроксимації для тонкого судна. Для цього класу суден теорія стрічки і її зміни дають гарні результати. Однак, у випадку коротких хвиль ці методи мають тенденцію недооцінювати додатковий опір. Це формує складну задачу у випадку, якщо робляться спроби оптимізувати форму корпусу. Якщо судно має тупу форму корпуса, локальний сталий рух впливає на величину додаткового опору дуже сильно. У роботі [16] представлено часовий метод, який може вирішити задачі з різними видами лінеаризованих формулювань. Як вхід, програма може використати незбурений потік, потік подвійного корпуса або нелінійний сталий рух. Розглядається симетричне судно з рівномірним плаванням з постійною швидкістю U у хвилях, які поширюються в напрямку, що замикається з передовим напрямком судна. Водна глибина h позначається постійною. Рідина прийнята ідеальною. [16]
