Нужна помощь

[AndrewM 0]

Добрый вечер.

Вопрос состоит в следующем:

При использовании пристеночных функции какое граничное условие использовать для давления? Для скорости использую "закон стенки".

Граничные условия для k и e имеют смысл только при y+ [30 , 100]. Как поступать с теми ячейками для которых y+ меньше.

В самом простом случае можно сначала сократить расчетную область на расстояние соответствующее выбранному y+ (перевод в размерное расстояние по приближенной формуле)?

Спасибо за ответы на глупые вопросы

[Lynx 1]

Не знаю, поможет ли, но есть такая вещь как "гибридный закон стенки". Но он применяется несколько наоборот - если используется низкорейнольдсовый вариант модели турбулентности, а у+ местами получается больше, чем нужно (т.е. где-то больше 1). Единственная ссылка, что у меня есть на эту штуку - на немецком: Rung, T. 1999. ‘Formulierung universeller Wandrandbedingungen fur

Transportgleichungsturbulenzmodelle’, Institutsbericht Nr. 02/99,

Hermann-Fottinger-Institut fur Stromungsmechanik, Technische

Universitat Berlin.

Но вот краткое описание из методологии Star-CD

Hybrid_wall.doc

[AndrewM 0]

Спасибо за предоставленную информацию.

Если я правильно понял, то данный метод состоит в следующем:

Первый расчетный узел лежит не на твердой поверхности, а на некотором (любом малом) удалении от нее. И для получения граничных условий в нем (для k и e) мы пользуемся формулами из пердоставленного документа. А для составляющих скорости пользуемся "стандартным законом стенки"

И отличие состоит в том что y+ может лежать в следующем интервале: 0.1 < y+ < 100. Граничные условия меняются на каждой итерации в зависимости от значения средней скорости в первом не граничном расчетном узле.

Но в этом тексте не хватает формул для расчета K_c, y+_m и непонятно что делать с давлением.

[Lynx 1]

На сколько я понял, все не совсем так.

По-моему там предлагается считать интенсивность турбулентности k в пристеночной ячейке обычным образом, с обычными граничными условиями. k_c - это соотвественно полученное таким путем значение k, оно не на стенке, а в центре пристеночной ячейки (на расстоянии y_C от стенки). y+_m считается вот из такого уравнения:

y+_m - (1/kappa)*ln(E*y+_m)=0

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

А касательные напряжения на стенке считаем через взвешенное среднее, с весовым коэффициентом дзета

[AndrewM 0]

Спасибо за ответ.

Попробую реализовать данную схему...

Пристеночная ячейка это та, у которой одна из граней совпадает с твердой поверхностью?

Для epsilon мы пользуемся алгебраическим соотношением только в пристеночной области или во всей расчетной области?

[Lynx 1]

Да, пристеночная это она самая.

Алгебраическое соотношение для epsilon используем только для пристеночных, а дальше как обычно решаем дифференциальное уравнение

[AndrewM 0]

Еще хотелось бы узнать какие граничные условия накладывать на турбулентную вязкость в модели Спаларта-Аллмараса (на входе, выходе и твердой стенке)?