Определение Cx и Сy для профиля крыла

[dmiantosha 11]

Сетка гексами, высокого качества... 

 

Думаю, что проволока вносит не более 1% в полученный результат..

[kristeen 355]

@@kol, спасибо за наводку на Накаяму.

 

В процессе поисков нашлась такая картинка:

 

 

@@dmiantosha, знакомый тренд? 

 

source

 

PS Чтотт для Cdrag они снова график от Clift вставили...

[dmiantosha 11]

Блина... ну что же так все сложно... неужели придется простейший профиль в 3D дуть.. =((((((

[kol 12]

Блина... ну что же так все сложно... неужели придется простейший профиль в 3D дуть.. =((((((

Почему? что то я в тексте поста проблемы для оценки профильных характеристик (2D) не увидел...

[kol 12]

Пошаманил с сеткой, взял k-e, вот что получил:

 

p2-ke.jpg

 

Не смотрите на первую точку графика при 0 градусов, там это мой косяк, ввел не те начальные условия.

 

По Сх совпадение очень хорошее.... По Cy - нет.... Опустим отрыв потока на больших углах атаки, почему даже на мелких угла (1,2,3 градусов) я получаю большую подъемную силу, чем в эксперименте (площадь то постоянна)... Есть какие-нибудь идеи по этому поводу?

Турбулентность потока... В каждой А. трубе разный уровень турбулентности, поэтому значения  продувок в разных трубах отличаются и не всегда получается выдержать турбулентность по "ГОСТ"  или вообще измерить и потом учесть. Про проволочку в опытах Накаямы - может не 1%, зависит от числа Рейнольдса т.е. при низких числах уровень турбулентности оказавыет значительно большее влияние, также важна скорость потока сама по себе (влияет через "энергию" потока - написал энергия так как не помню как правильно называется) : http://m-selig.ae.illinois.edu/uiuc_lsat.html - посмотрите про их баблы (пузырики) :)

[kristeen 355]

@@kol, плз, дай ссыль на ту работу Накаямы, которую постоянно поминаешь в суе.

[dmiantosha 11]

В общем Алилуя...

 

И слава я..цам...

 

И терпение и труд все перетрут...

 

В общем если считать в 3D (стационар по шагам от пред. решения к след) :

 

Качественно уже очень похоже на правду. Видно, что ошибка растет с ростом угла, почти уверен, что смогу ее сильно уменьшить, пошаманив сеткой...

 

После 16 градусов (критическим углом является 16.7) наблюдается нестационар:

 

Попробую посчитать на этих углах в трансиенте на этой же сетке, так как сама сила ведет себя очень странно, чем дольше считаешь, тем меньше она становится (последний угол прервал, надоело считать, картина понятна):

 

Оказывается концевой эффект влияет и очень сильно, видно, что давление меняется примерно с половины крыла:

 

Печаль, что приходится считать в 3D..

 

Вижу для себя такой путь решения:

 

1. Постараюсь получить как можно близкий результат к эксперименту.

2. Начну последовательно уменьшать и упрощать модель. 

П.С.: Сильно ли влияет вот это:

 

A wall has been placed at portion(s) of an OUTLET                

boundary condition (at   5.0% of the faces,   5.9% of the area)  

to prevent fluid from flowing into the domain.

 

Начинается на 6 градусах и далее нарастает... 5% - это при 20 градусах.. .на 6 градусах 2.5 

[kristeen 355]

П.С.: Сильно ли влияет вот это:

Вихрь ложится на outlet. Придется отодвинуть выходную границу подальше от крыла.

[kol 12]

 

П.С.: Сильно ли влияет вот это:

Вихрь ложится на outlet. Придется отодвинуть выходную границу подальше от крыла.

 

"A wall has been placed at portion(s) of an OUTLET                

boundary condition (at   5.0% of the faces,   5.9% of the area)  

to prevent fluid from flowing into the domain.

 

Начинается на 6 градусах и далее нарастает... 5% - это при 20 градусах.. .на 6 градусах 2.5 " - это на основе этого сообщения ты сделал такое умозаключение (вихрь наехал на выход)?

@@kol, плз, дай ссыль на ту работу Накаямы, которую постоянно поминаешь в суе.

Вообще нет прямой ссылки на эту инфу... Есть много источников где эта работа используется, я перерыл думаю все в своё время и всё равно собрал "не полный вариант". Что конкретно тебя интересует в этой работе? 

 

В общем Алилуя...

 

И слава я..цам...

 

И терпение и труд все перетрут...

 

В общем если считать в 3D (стационар по шагам от пред. решения к след) :

grafik-good.jpg

 

Качественно уже очень похоже на правду. Видно, что ошибка растет с ростом угла, почти уверен, что смогу ее сильно уменьшить, пошаманив сеткой...

 

После 16 градусов (критическим углом является 16.7) наблюдается нестационар:

solver.jpg

 

Попробую посчитать на этих углах в трансиенте на этой же сетке, так как сама сила ведет себя очень странно, чем дольше считаешь, тем меньше она становится (последний угол прервал, надоело считать, картина понятна):

solver-force.jpg

 

Оказывается концевой эффект влияет и очень сильно, видно, что давление меняется примерно с половины крыла:

CFX_004.jpg

 

Печаль, что приходится считать в 3D..

 

Вижу для себя такой путь решения:

 

1. Постараюсь получить как можно близкий результат к эксперименту.

2. Начну последовательно уменьшать и упрощать модель. 

П.С.: Сильно ли влияет вот это:

 

A wall has been placed at portion(s) of an OUTLET                

boundary condition (at   5.0% of the faces,   5.9% of the area)  

to prevent fluid from flowing into the domain.

 

Начинается на 6 градусах и далее нарастает... 5% - это при 20 градусах.. .на 6 градусах 2.5 

 

"Печаль, что приходится считать в 3D.." - почему такая обречённость про три Д? профиль берёте АТ15 или оригинальный naca (по внешнему виду они одинаковые)?

[dmiantosha 11]

Печалюсь что получаю 1 млн и более... очень долгий расчет получается...

 

A wall has been placed at portion(s) of an OUTLET   - получаю из-за того, что немного хитрю с поворотом домена. Создаю цилиндр, торец - симметрия, вырезаю в нем крыло, делю боковую поверхность пополам, левая сторона - вход, правая - выход. Далее в Pre поворачиваю через Transform Mesh домен на нужный угол... таким образом IN и OUT тоже смещается... Видимо из-за этого я получаю это предупреждение...

 

Но вопрос то все тот же... Сильно ли влияет на результат?

[ccc44 53]

Оказывается концевой эффект влияет и очень сильно, видно, что давление меняется примерно с половины крыла:

Согласно классике распределение циркуляции скорости и подъемной силы вдоль размаха носит эллиптический характер. Хорошее совпадение на линейном участке свидетельствует о отм, что экспериментальные данные получены на крыле конечного размаха с удлинением 5. Методика эксперимента изложена, например, в пособии по лабораторному практикуму ХАИ http://girkay.ho.ua/Stud/Exp_Aero.pdf. Похоже, что с 1930 года она не изменилась. Строго говоря называть такие результаты "аэродинамическими характеристиками профиля" некорректно.

[dmiantosha 11]

На той же сетке делал Mesh Adaptation. Как я и думал ранее, дело именно в ней, так как сразу получил результаты гораздо лучше.

Для углов 18 градусов и далее, потом попробую в трансиенте, но так нельзя сделать Mesh Apadptaion, возможно попробую использовать сетку с пред. шагов.

 

[ccc44 53]

На той же сетке делал Mesh Adaptation. Как я и думал ранее, дело именно в ней, так как сразу получил результаты гораздо лучше.

Это профиль РII-10 АТ-1, стр 40-41? Если да, то здесь Re=845000, поэтому графики более гладкие. Отрыв на 16,7° вы должны увидеть.

[dmiantosha 11]

Экспериментальных данных нет на 16.7... Потом в трансиенте все вместе посчитаю, попробую поймать десятые угла...

[dmiantosha 11]

Попробовал посчитать 18 градусов в трансиенте, ожидал какой-то пульсации Cy... ожидания не оправдались... все ложится на горизонтальную прямую...

[kol 12]

@@kol, плз, дай ссыль на ту работу Накаямы, которую постоянно поминаешь в суе.

Вот то чем я пользовался:http://cfd.mace.manchester.ac.uk/cgi-bin/cfddb/prpage.cgi?11&EXP&database/cases/case11/Case_data&database/cases/case11&cas11_head.html&cas11_desc.html&cas11_meth.html&cas11_data.html&cas11_refs.html&cas11_rsol.html&1&0&0&0&0

[a_schelyaev 366]

 

@@kol, плз, дай ссыль на ту работу Накаямы, которую постоянно поминаешь в суе.

Вот то чем я пользовался:http://cfd.mace.manchester.ac.uk/cgi-bin/cfddb/prpage.cgi?11&EXP&database/cases/case11/Case_data&database/cases/case11&cas11_head.html&cas11_desc.html&cas11_meth.html&cas11_data.html&cas11_refs.html&cas11_rsol.html&1&0&0&0&0

 

А как там материалы скачать - пароль требует и логин?

 

P.S. Оригинальная работа.

[kristeen 355]

А как там материалы скачать - пароль требует и логин?

Если еще не нашел, то вот точки для "профиля А":

#          AIRFOIL GEOMETRY  DSMA661(MODEL A)
#    J      X/C      Y0/C    PHI(DEG)      Y0/C     PHI(DEG)
#                    UPPER SURFACE         LOWER SURFACE
variables="index","x/c","y/c_upper","phi_upper","y/c_lower","phi_lower"
     1    0.000     0.00000   90.00       0.00000    90.00
     2    0.002     0.00739   59.29      -0.00919   124.75
     3    0.005     0.01030   45.00      -0.01037   142.77
     4    0.010     0.01501   33.50      -0.01478   151.22
     5    0.025     0.02253   21.92      -0.01935   167.58
     6    0.050     0.03076   15.74      -0.02330   172.74
     7    0.100     0.04212   10.37      -0.02866   174.59
     8    0.150     0.04979    7.37      -0.03310   175.27
     9    0.200     0.05534    5.35      -0.03677   176.41
    10    0.250     0.05929    3.78      -0.09390   177.87
    11    0.300     0.06205    2.58      -0.04037   179.71
    12    0.350     0.06379    1.36      -0.03985   181.44
    13    0.400     0.06436   -0.10      -0.03799   182.72
    14    0.413     0.06378   -1.42      -0.03568   183.36
    15    0.450     0.06360   -1.62      -0.03527   183.42
    16    0.493     0.06195   -2.72      -0.03259   183.66
    17    0.500     0.06161   -2.89      -0.03214   183.68
    18    0.550     0.05860   -3.99      -0.02888   183.77
    19    0.593     0.05524   -4.93      -0.02603   183.80
    20    0.600     0.05463   -5.07      -0.02557   183.80
    21    0.650     0.04973   -6.10      -0.02225   183.78
    22    0.693     0.04483   -6.90      -0.01943   183.74
    23    0.700     0.04397   -7.01      -0.01897   183.73
    24    0.750     0.03745   -7.83      -0.01573   183.69
    25    0.793     0.03132   -8.37      -0.01298   183.65
    26    0.800     0.03028   -8.43      -0.01253   183.65
    27    0.843     0.02383   -8.60      -0.00980   183.62
    28    0.850     0.02277   -8.60      -0.00936   183.62
    29    0.893     0.01626   -8.61      -0.00664   183.61
    30    0.900     0.01520   -8.61      -0.00620   183.60
    31    0.918     0.01248   -8.61      -0.00507   183.57
    32    0.940     0.00916   -8.54      -0.00372   183.45
    33    0.950     0.00766   -8.46      -0.00313   183.36
    34    0.970     0.00473   -8.23      -0.00200   183.10
    35    0.975     0.00400   -8.18      -0.00173   183.04
    36    0.990     0.00188   -7.91      -0.00096   182.73 
    37    1.000     0.00051   -7.68      -0.00051   182.50 

Источник

 

Там даже сетки есть )

 

@@kol, спасибо.

[a_schelyaev 366]

Бладарю!

[dmiantosha 11]

В итоге.. изначально у меня была сетка в 0.7 млн. элементов, теперь пересчитал на 2 млн, если сделать Mesh Adaptaion, лучше не становится, то есть по секте я уже уперся...

 

Стало действительно лучше, но особенно после срыва потока ошибка 50%, что нельзя назвать удачным совпадением... 

 

Остались еще у кого мысли по этому поводу? Что попробовать?