Fletcher-2-rus (1185919), страница 50
Текст из файла (страница 50)
Хотя конкретный вид уравнений в этих двух случаях будет различным. Отличительные свойства КИ8-уравнений и описываемых ими течений будут рассмотрены в $ 16.1. Дополните,чьные предположения и формулы для внутренних течений, например в диффузорах и каналах, будут приведены в 9 16.2. Для внешнего течения около тела (или тел) решение невязких уравнений (гл.
14) дает очень хорошее приближение решения на удаленной границе. Для вязкой области вблизи тела полезно при рассмотрении КЬ15-уравнений различать сверхзвуковое и дозвуковое течения (9 16.3). Для полноты картины методов расчета внешних течений будут рассмотрены некоторые традиционные методы расщепления течений на невязкое течение и течение в пограничном слое. При этом будут рассмотрены более сложные случаи, такие, как образование малых отрывных зон и взаимодействие скачка с пограничным слоем. Наиболее развитые методы, основанные на модели невязкого течения и течения в пограничном слое, относятся скорее к ЙИВ-типу, а не к чисто погранслойным течениям.
2 16.1. Введение В гл. 1 отмечалось, что быстрое развитие вычислительной гидродинамики было вызвано ее широким внедрением в процессы проектирования. При проектировании оборудования, связанного с движением жидкости (обычно воздуха или воды), часто необходимо выбрать единственный режим, при котором дан- 6 !6.1. Введение ное оборудование (например, турбина, диффузор или самолет в целом) будет работать наиболее эффективно. При этом снижение эффективности работы оборудования при выборе близких «нерасчетных» режимов не должно быть слишком велико. Основная роль вычислительной гидродинамики в процессе проектирования состоит в детальном определении стационарных свойств течения в расчетном режиме путем решения соответствующей системы уравнений с необходимыми граничными условиями.
Во многих практически важных задачах имеют место течения с доминирующим направлением. Это означает, что при определенном выборе декартовой системы координат выполняется условие и )> и, ш, если положительное направление оси х является доминирующим. Анализ уравнений (п. 16.1.1) показывает, что вязкие диффузионные члены вдоль основного направления значительно меньше вязких диффузионных членов в поперечном направлении и, следовательно, могут быть отброшены. Данное свойство является ключевым для уравнений пограничного слоя (гл.
15). Если можно пренебречь изменением давления поперек пограничного слоя, что справедливо при малой кривизне линий тока, то для решения задач можно воспользоваться весьма эффективными маршевыми в направлении течения методами. Не удивительно поэтому, что проектирование оборудования часто основано на двух этапах: расчет невязкого течения плюс коррекция на основе теории пограничного слоя. На первом этапе течение предполагается невязким и нетеплопроводным. В результате решения обычно получается распределение давления по поверхности тела (см. гл. 14).
На втором этапе в результате учета вязких и тепловых эффектов в пограничном слое (гл. 15) производится коррекция полей скорости и температуры вблизи поверхности тела. Решение уравнений пограничного слоя позволяет определить сдвиговые напряжения и тепловые потоки на поверхности тела и с учетом распределения давления определить сопротивление тела. В случае необходимости распределение давления может быть подправлено в результате учета толщины вытеснения (п.
14.1.4). Основной целью проектирования оборудования, связанного с движением жидкости, является повышение его эффективности. При этом необходимо добиться выполнения тех же целей на оборудовании меньшего размера, а значит, более дешевом, Данная тенденция вскрывает некоторые недостатки, присущие традиционным методам проектирования, основанным на коррекции невязкого решения на основе приближения пограничного слоя.
Сказанное можно пояснить на примере проектирования 294 Гл. 16. Течения, описываемые К1Ч5-уравнениями Навье — Стокса самолета, максимальная эгфективность которого тесно связана с отношением подъемной силы к силе сопротивления отдельного крыла [Кпс)тегпапп, 1978[. Характерные зависимости подъемной силы и силы сопротивления крыла от угла атаки приведены на рис.
16.!. Максимум отношения подъемной силы к силе сопротивления приходится на большой угол атаки и не слишком удален от точки потери скорости. При оптимальном угле атаки толщина вязкого слоя Е о 20 а га Угол *така, гРал Угол атак .гРал Рис. 16.1. Зависимость подъемной силы и силы сопротивления от угла атаки. на подветренной стороне крыла велика, особенно за ударной волной, которая может возникнуть. Вблизи задней кромки в направлении, поперечном набегаюшему потоку, образуется значитсльный градиент давления [)ч!а)сауаша, 1984).
Данный градиент связан со слиянием подветренного и наветренного потоков. На некоторых нерасчетных режимах крыло будет работать еще ближе к точке потери скорости. При этом произойдет дальнейшее увеличение толщины вязкого слоя на подветренной стороне. Приближение пограничного слоя в таких условиях неприменимо [Ногз1гпап, 1984].
Обращение к решению полных уравнений Навье — Стокса возможно, но неэкономично, поскольку решения придется находить в широком диапазоне параметров проектирования: угол атаки, геометрия крыла, число Рейнольдса, число Маха и т. д. Желательно иметь уравнения, с помощью которых можно было бы моделрровать физические явления (почти) так же точно, как на основе уравнений Навье — Стокса. При этом необходимо, чтобы решение этих уравнений было (почти) столь же экономичным, как и решение уравнений пограничного слоя. Данные уравнения в дальнейшем будут называться укорочснньгми уравнениями Навье — Стокса (КИБ-уравнениями). $16.1.
Введение 295 Свобода определения предполагает возможность интерпретации различных систем уравнений как Кг(Ь-уравнений. Однако всем К148-системам присущи три следующих свойства: 1. Существует доминирующее направление течения, примерно совпадающее с одной из координатных линий.
2. Вязкой диффузией и теплопроводностью в маршевом направлении можно пренебречь по сравнению с вязкой диффузией и теплопроводностью в направлении, поперечном маршевому. 3. КМЬ-уравнения сводятся к уравнениям Эйлера, если пренебречь вязкой диффузией и теплопроводностью в поперечном направлении. Вывод основных КМЬ-уравнений основан на сравнении порядков величин (п. 16.1.1), аналогичном использованному при выводе уравнений пограничного слоя (5 11.4).
Однако в приближенном виде сохраняется уравнение поперечной компоненты импульса, поскольку физика явлений такова, что можно ожидать появления заметных градиентов давления в этом направлении. Для внешних течений важно сохранение всех членов системы уравнений Эйлера. Лишь в этом случае К1чЬ-уравнения применимы во всей расчетной области. На практике при расчете не- вязких областей можно применять и другие методы, но использование глобально применимых уравнений облегчает проведение сшивки решений в вязкой и невязкой областях. Требование сравиимости по экономичности с уравнениями пограничного слоя предполагает возможность построения в вязкой области алгоритма решения, маршевого в направлении течения.
В свою очередь это означает, что в направлении маршевой переменной уравнения должны быть неэллиптическими. Следовательно, на выходной границе ие надо определять граничных условий. Для внутренних течений и сверхзвуковых внешних течений точное решение обычно получается при помощи одного прохода по маршевой переменной. При использовании КХЬ-уравнений для расчета внешних дозвуковых течений требуется повторение проходов (итераций) по маршевой переменной, особенно если образуются небольшие области возвратного течения (отрывные зоны).
Для дозвуковых внешних течений система уравнений имеет эллиптический характер и необходимо задание граничных условий на выходной границе. Однако КМЬ-подход остается экономичным, поскольку, если отрывные зоны малы, достаточно лишь нескольких итераций (проходов по маршевой переменной). В данной главе укороченные уравнения Навье — Стокса будут рассмотрены лишь для решения стационарных задач, где эти уравнения имеют больше всего применений. Однако можно 296 Гл. 16. Течения, описываемые К1Ч5-уравнениями Навье — Стокса рассмотреть и определенный класс нестационарных течений, для которых существует доминирующее направление течения.