Как вырастить сетку 2D-airfoil правильно?

[zerganalizer 601]

Я сделал работающую версию генератора 2D-сетки для пакетного CFD-анализа профилей. Она полностью на питоне и строит сетку через API GMSH 4.13.1. Пока это гибридная сетка, но с настройками пограничного слоя, загущения по кривизне, в районе хвостика и следа профиля.

 

Более того - заставить ее считаться в SU2 - было для меня сложно. Заработал расчёт только после указания 'MGLEVEL': '2' и создания в сетке ещё одной физического маркера nearfield и назначить 'MARKER_EULER': '( nearfield )' как границы между прямоугольного пограничного слоя и остальной треугольной сеткой.

 

И вот общий вид - вертикальная линия и хвост просто выведены для вопросов, на сетку не влияют.

 

 

Крупнее:

Район носа:

Хвост:

У меня в настройки вынесено - на сколько сдвинуть профиль относительно лобовой полудуги farfield. Сейчас Move_airfoil=1.1 - это значит носик задвинут вперёд на хорду + 10%, и я не знаю, насколько это правильно. Вы видите - airfoil слева от вертикальной линии "нуля" центра дуги. А вдруг он должен быть справа, посредине или ещё как?

 

Вот базовые настройки:

 

# base mesh sizes
fact = 1                               #Мастштаб ячеек
lc1 = 0.01 * fact                  #Минимальная ячейка
lc2 = 1 * fact                       #Максимальная ячейка
LenNET=50                         #Удаление farfield от профиля, хорд
Move_airfoil=1.1                #Сдвиг профиля от нуля полудуги
Загущение_профиля = 5  #Загущение профиля по кривизне - соотношение размеров max/min ячеек
Загущение_сетки = 90     #Загущение профиля и сетки спереди - сколько точек на 360 градусов дуги.
Thickness=0.01*fact          #Толщина погранслоя профиля
Ratio=1.3                           #Коэффициент роста погранслоя
Size=2.7e-4                       #Начальная толщина ячейки погранслоя

 

Есть ещё куча параметров - по дополнительным загущениям, по алгоритму обработки хвостика профиля и прочее.

 

Вот интересный вопрос - Загущение_сетки = 90, если сделать 720, то лобовая часть перехода от полудуги к носику профиля станет гораздо плотнее, от этого есть профит???

Но меня мучит вопрос - я сделал не по образцам обучалок мешер - им вообще можно пользоваться в SU2 и других CFD? Я видел подобные вопросы про гибридные сетки на CFD-форуме, но информации мало.

 

Или вообще лучше сразу делать упорядоченную трансфинитную, как в родном мешере SU2_EDU? Её гораздо сложнее делать - как по мне. А ещё куча вопросов по соотношению настроек продувки и настроек мешера сетки. Я уже упоминал - для пакетной продувки под разные скорости и углы атаки для точности нужно делать сетку "на лету" и под текущие режимы продувки.

 

Edited June 17 by zerganalizer

[a_schelyaev 410]

В разных отчетах могут по разному приводить данные. Где то исходные и обработанные, где-то сразу обработанные.  Идеальных данных нет, как и не бывает двумерных продувок.

[zerganalizer 601]

7 часов назад, a_schelyaev сказал:

Идеальных данных нет, как и не бывает двумерных продувок.

После предварительной калибровки именно сетки получил на NACA 2412 +/- 90 градусов на 0.1 МАХа и RMS_DENSITY 6+, вблизи больших значений 7+. Гонять детальнее буду в понедельник.

[zerganalizer 601]

Тестирую новую калиброванную по 2-м результатам продувок NACA 2412 и конфиг, на котором я калибровался. Не факт, что при сильно другой скорости результаты линейно (корректно) смасштабируются, но цифры, какие я в баттле получил для профилей по вашим ссылкам теперь может чуть ниже. Словом, я пытаюсь начать баттл на Re 100.000 с тем же шагом до 0.1 МАХа, и поточнее.

 

Пока "кларк" более-менее норм считается, хотя точность я чутка занизил для быстроты теста. Если качество расчёта меня устроит во всём этом диапазоне - я приму это решение за текущий стандарт для дальнейших баттлов.

 

Это решение (мешер + CFD) пока работает от 0.001 МАХа до примерно 1.4 МАХа в диапазоне +/-30 градусов, хотя бОльшая часть диапазона считается +/-90. Точность примерно такая, расчёт в процессе.

 

Запуск расчета: Mach=0.029176066127226598, AOA=16.0
CEff = 9.045397418, Итераций: 1225, RMS_DENSITY = -7.226307772, 
время расчета: 15.9 сек
RE_NUMBER = 698094.9592580154
Режим RANS (Mach = 0.03)

 

 

 

А вообще поищу способ подкалиброваться для более гладких и быстрых поляр. И подстроить визуальное оформление - я введу в набор данных имена профилей, чтобы автоматом их в легенду выводить. Есть куда копать, сюда например:

 

Edited Tuesday at 01:37 PM by zerganalizer