Точность и погрешность вычислений

[ac100500 0]

Добрый день. Хотел уточнить у специалистов.

Есть простая задача продувки профиля крыла с построением графиков Cl(Cy) от Альфа, Сd(Cx) от Альфа.

Тип задачи: 2D

Провожу расчеты с шагом угла - 1 градус. Например от 0 до 22. Скорость как граничное условие высчитываем.

Данные заносим в таблицу. 

Координаты профиля берем отсюда - http://airfoiltools.com/airfoil/details?airfoil=n2414-il

Результаты имеют отличие - а именно коэффициент лобового сопротивления до 200% на углах до 10 градусов.

Результаты на сайте вышеуказанном получены из софтины и принимаются якобы как достаточно достоверные.

В чем может быть причина - если не учитывать варианты сетки и граничных условий.

Решения при стационарной постановке сходятся на 1600-2100 иттерации.

Получается что много есть примеров обтеканий тел и т.д. но во всех роликах и примерах значение Cd при нулевом угле атаки не меньше чем 1.5е-02

что как бы уже очень много. 

 

 

 

 

Изменено 1 февраля 2017 пользователем ac100500

[piden 2 874]

1 hour ago, ac100500 said:

В чем может быть причина  - если не учитывать варианты сетки и граничных условий.

Все равно лучше начать с сетки и ГУ. О качестве сетки люди иногда имеют весьма своеобразное суждение.

Про Re вы не упоминаете. Вдруг таки скорость неверно вычислили. Нет ведь ни файла проекта, ни детального описания настроек / полученных результатов. Как вычисляются Cx и Cy - встроенными средствами флуента, или в посте по своим формулам? 

 

Ну и наконец результаты из xfloil не эталон. Была когда-то тема по профилям, можете от туда взять какой-нибудь профиль (ссыль в первом посте) и сравниться с экспериментом из книги и с тем, что другие насчитали.

1 hour ago, ac100500 said:

Получается что много есть примеров обтеканий тел и т.д. но во всех роликах и примерах значение Cd при нулевом угле атаки не меньше чем 1.5е-02

что как бы уже очень много. 

А это наводит на мысль, что у вас какое-то свое представление о коэффициенте сопротивления)

Изменено 1 февраля 2017 пользователем piden

[ac100500 0]

По вопросу коэффициента лобового сопротивления в большинстве ламинарных профилей на нулевом угле атаки не должен превышать 1е-02 при Re=1e+06.

В задачах с малым Re (медленные маленькие радиомодели) значения Cx больше. 

 

[ac100500 0]

очень жду помощи специалистов

 

Все вышеописанное касается тестового профиля PW75.... текущие расчеты для угла 5 градусов

Подчеркнул - коэффициент который должен получится - получается - всегда в раза 2-3 больше.

Re - 1.000.000

[ac100500 0]

Пересчитал с моделью SST

 

[ac100500 0]

Продолжаю дальше разбираться в возможности получения подобных данных используя численный метод.

На настоящий момент есть непонимание какие данные необходимы в Reference values. А именно: AREA, LENGTH, DEPTH.

Можно предположить что LENGTH - это длина хорды, AREA - площадь крыла (при 2D постановке вопрос что за площадь), а вот с глубиной большой вопрос.

Во многих видео - полученные коэффициенты к реальным не имеют отношения.

И именно от этих данных зависит и расчет сил и соответсвенно - Cy, Cx.

[Yuri69 1]

DEPTH=1м для 2D-задач

[soklakov 1 475]

В 01.02.2017 в 20:35, ac100500 сказал:

Результаты на сайте вышеуказанном получены из софтины и принимаются якобы как достаточно достоверные.

Как-то у нас был на стажировке студент из Германии. Перед ним стояла задача оптимизации композитного крыла и, соответственно, на одном из этапов нужно было определить сопротивление этого крыла.

Он тоже сравнивался с какой-то программой, пытаясь получить похожие результаты в CFX или Fluent.

Результаты в CFX и Fluent получались одинаковые, но все никак не сходились с этой программкой. Многочисленные проверки его решения не находили никаких ошибок. Ему понадобилось еще пару недель безрезультатных исследований, пока он не решился как следует проверить референсную программку - таки да, несмотря на многолетний опыт её применения, в ней были ошибки.

Так что я бы отнесся повнимательнее к ранее высказанному утверждению:

В 01.02.2017 в 22:14, piden сказал:

 

Ну и наконец результаты из xfloil не эталон.

 

[soklakov 1 475]

В 11.02.2017 в 19:02, ac100500 сказал:

И именно от этих данных зависит и расчет сил и соответсвенно - Cy, Cx.

Строго говоря, в результате расчета Вы получаете давления. Их интеграл по площади - это уже силы. Интеграл программа посчитаем сама без проблем. А чтобы эти силы перевести в коэффициенты нужны как раз все эти AREA, LENGTH и DEPTH. Понимаете? То есть это коэффициенты в формуле вычисления Cy и Cx. Можно выбрать их единичными, тогда значения вычисленных коэффициентов будут равны самим силам.

А пост-процессор, меж тем, может и непосредственно силы выводить.