Некорректное определение расхода моделями RANS

[tik4 0]

очень будет интересно посмотреть

<noindex>http://narod.ru/disk/42826839001.9e9ef9d49...0.4_ke.res.html</noindex>

<noindex>http://narod.ru/disk/42825433001.13fa6323d...0.4_ke.def.html</noindex>

<noindex>http://narod.ru/disk/42825172001.20de5fc60...d/base.prt.html</noindex>

<noindex>http://narod.ru/disk/42827258001.b73795576...495/1.djvu.html</noindex>

<noindex>http://narod.ru/disk/42827321001.bed4ef78e...1f2/2.djvu.html</noindex>

<noindex>http://narod.ru/disk/42827359001.37d8f1a5a...a78/3.djvu.html</noindex>

<noindex>http://narod.ru/disk/42826942001.d4ffec081...5e6/4.djvu.html</noindex>

вот ссылки на def, res, prt и Идельчика (файлы 1-4, перфорированные плиты в файле 3, читать со страницы 401, см. диаграммы 8-3 и 8-5, а после еще 8-7 для корректировки по скорости потока)

сопротивление по Идельчику примерно 41 (это с учетом переходного течения) * 1.2 (коэффициент скорости потока, там Мах примерно 0.08) итого грубо 49...50 (либо 56...57, если считать течение турбулентным). А расчетное примерно 119 (т.е. в 2 раза выше, Re = 29000) по k-e (на остальных моделях примерно столько же, в т.ч. на SST Gamma Theta model (с учетом переходного течения...якобы :)).

Завтра график выложу по 4...5 точкам для данного случая в зав-ти потерь от Рейнольдса по Идельчику и расчету CFX, там все наглядней видно, для Re до 10000 все более менее сходится, а вот выше начинается разбег

Изменено 5 марта 2012 пользователем tik4

[tik4 0]

нашел ошибку в геометрии, там вышеприведенные файлы имеют отношение площадей f = 0.156 отсюда такое сильно завышение, перестроил под 0.2 и картина немного прояснилась (файлы если нужно скину позже, графики тоже, ибо сейчас надо часть точек пересчитать - не успел, работы много :))

пересчитал на f = 0.2 l/d = 0.4 Re=3.3*10^4

пишу примерные результаты... до сотых долей уж не помню)

1. Расчетное сопротивление = 71 (k-e)

2. Идельчик (с учетом сжимаемости и при полностью турб. течении Re > 10^5) = 66

3. Идельчик (с учетом сжимаемости и при переходном течении Re < 10^5) = 58

4. Идельчик (без учета сжимаемости и при переходном течении Re < 10^5) = 41

пытка людей, считавших лопатки Gidra дала информацию, что сжимаемость в ней реализована плохо...

т.е.

1) В CFX переходное течение учитывается неважно (даже спец моделью типа Gamma Theta), но при предположении полностью турбулентного потока с Идельчиком бьется неплохо (для Re 3.3*10^4 погрешность 7...7.5 %, при более низких Re она еще меньше).

2) Но так как течение все-таки переходное то погрешность поболее (CFX завышает сопротивление на 20...25 % относительно Идельчика). Отсюда объясняется расчетное затекание в лопатку, тогда как реальности его не было (возможно перепад на дырке был и не очень большой, но все же был в +)

3) Так как Gidra плохо учитывает сжимаемость, то видно, что разница и впрямь может быть значительной (между 1 и 4 (см. выше) разница 73%), т.е. врут оба расчета, а правда почти посередине :) но к CFX поближе

ну вот вроде все стало и понятно, файлы скину после праздников, если хочется, можно и тот prt, что я сбросил переделать под f = 0.2 и сравнить самим, сетка строится быстро. Можно сравнить и нам том, f что есть, но в Идельчике нет распределения коэффициентов влияния сжимаемости для f = 0.15, а только начиная с 0.2

[IGL 8]

а с реальными газами не пробовали подсчитать?

[tik4 0]

а с реальными газами не пробовали подсчитать?

нет, но можно будет попробовать...не раньше воскресенья :)

мы обычно всегда использовали для воздуха Air Ideal Gas

[tik4 0]

нет, но можно будет попробовать...не раньше воскресенья :)

мы обычно всегда использовали для воздуха Air Ideal Gas

прикинули, принципиально разницы не изменились, только все сдвинулось вниз по Рейнольдсу из-за изменившейся динамической вязкости - у Air Ideal Gas она задана жестко 1.831*10^-5 и не зависит от параметров потока... шаманим дальше :)

[a_schelyaev 366]

Если есть модель турбулентности, то она будет определеять уровень вязкости.

[tik4 0]

Если есть модель турбулентности, то она будет определеять уровень вязкости.

она будет добавлять, насколько я помню, турбулентную вязкость, но заданная вязкость среды также влияет :) а проблему с уровнем потерь пока не решили, выявили, что она особенно актуальна именно для низких F0/F1, где расчетное сопротивление начинает расти резко при меньших Рейнольдсах, чем в Идельчике. Для несильно пережатых дефлекторов данная проблема не так актуальна, там кривые влияния сжимаемости более пологие. А вот расчет несжимаемого газа почти совпал с несжимаемым Идельчиком