Санкт-Петербургский государственный университет
Кафедра вычислительных методов механики деформируемого тела
Макунин Алексей Вячеславович
Магистерская диссертация
Моделирование внутренних задач гидродинамики
на основе f - модели турбулентности
Направление 01.04.02
Прикладная математика и информатика
Магистерская программа надежность и безопасность сложных систем
Научный руководитель,
доктор физ.-мат. наук,
профессор
Павловский В.А.
Санкт-Петербург
2016
Содержание
Аннотация
3
Введение
4
1 Современное состояние теории турбулентности
1.1 Актуальные проблемы современной гидродинамики
1.2 Вычислительная гидродинамика . . . . . . . . . . .
1.3 Турбулентность в ANSYS . . . . . . . . . . . . . . .
1.4 f − модель турбулентности . . . . . . . . . . . . . .
.
.
.
.
8
8
11
14
16
2 Ламинарное течение в круглой трубе
2.1 Течение Пуазейля в трубе . . . . . . . . . . . . . . . . . . . .
2.2 Течение в начальном участке круглой трубы . . . . . . . . .
2.3 Закон сопротивления для течеия в трубах . . . . . . . . . .
19
19
21
21
3 Феноменологическая f − модель турбулентности
3.1 Физический смысл функции f . . . . . . . . . . . . . . . . .
3.2 Система уравнений модели . . . . . . . . . . . . . . . . . . .
3.3 Определение феноменологических констант . . . . . . . . .
27
27
34
37
4 Расчет течения жидкости в круглой цилиндрической трубе на основе f − модели
4.1 Постановка задачи о течении в трубе . . . . . . . . . . . . .
4.2 Численный расчет и сравнение результатов . . . . . . . . . .
4.3 Верификация на основе ламинарного течения . . . . . . . .
4.4 Интегрирование задачи в рамках f − модели . . . . . . . . .
46
46
49
52
53
Заключение
55
Список литературы
56
2
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
Аннотация
В данной магистерской диссертации производится численный расчет
характеристик турбулентного течения жидкости на основе f?модели турбулентности. В качестве объекта исследования взято течение несжимаемой
вязкой жидкости в круглой цилиндрической трубе. В исследовании задействованы методы механики сплошных сред, гидродинамики, численные
методы решения краевых задач, а также аналитические методы. Характеристики течений, рассчитанные по изложенным расчетным методикам,
сравниваются с известными экспериментальными данными и результатами, полученными ранее другими авторами.
3
Введение
К концу позапрошлого столетия развитие гидродинамики подошло
к этапу, когда ее дальнейшее существование как единой науки, использующей в своих построениях строгий математический аппарат, стало невозможным. Это объясняется тем, что ввиду быстрых темпов развития техники и появления новых технических устройств, возник спрос на количественные характеристики реальных величин, которые в силу сложности
рассматриваемых процессов, современная на тот момент гидродинамика
дать не могла. Это привело к тому, что дальнейшее развитие этой науки
пошло по двум различным путям. Первый – это теоретическая гидродинамика, основанная на уравнениях Эйлера движения жидкости при отсутствии сил трения. Иной подход – это принявшая отчетливый эмпирический характер, основывающаяся в большинстве своем на экспериментальных данных и радикально отличающаяся от теоретической методиками и
спектром решаемых задач – гидравлика [1].
Решение обширного спектра гидродинамических задач подразумевает применение строгого математического аппарата интегрирования дифференциальных уравнений, лежащих в основе математической модели, которая служит эквивалентом рассматриваемой задачи, а также граничных
и начальных условий. Суммарные характеристики, описывающие процесс,
находятся из довольно общих теорем механики - законов сохранения количества движения, массы, энергии и других. Большая сложность, однако,
состояла в том, что результаты, полученные методами классической гидродинамики, во многом шли вразрез с результатами опытов. Особенно резко
это коснулось в таких важных с практической точки зрения вопросах, как
расчет сопротивления, оказываемого жидкостью на движущееся в ней тело, нахождение характеристик потока жидкости вблизи твердой стенки и
др.
В начале прошлого века были выведены уравнения движения жидкости с учетом вязкого трения (уравнения Навье-Стокса). Но из-за больших
математических трудностей в нахождении их аналитических решений, применить эти уравнения для решения задач о течении жидкости при наличии
трения удавалось только редких частных случаях.
В идеальной жидкости касательные силы отсутствуют и на поверх4
ности соприкосновения твердого тела с жидкостью имеется разность касательных скоростей, то есть происходит скольжение жидкости. Напротив, в
действительной жидкости на обтекаемую твердую стенку передаются касательные силы, что приводит к прилипанию жидкости к стенке[2]. В данной
работе будет рассматриваться идеальная вязкая жидкость.
Следующим значимым шагом в развитии этой теории стало введение в 1904 году Прандтлем понятия пограничного слоя [2]. Эта гипотеза,
помимо своей физической наглядности в описании важной роли вязкости
в задаче о сопротивлении вязкой жидкости, дала возможность преодолеть
математические трудности и теоретически исследовать жидкости с трением. До сравнительно недавнего времени теоретические исследования динамики жидкости в большинстве своем проводились для идеальной (несжимаемой и не обладающей трением) жидкости. Влияние трения и сжимаемости в заслуживающей того степени стало учитываться гораздо позднее.
При движении жидкости без трения между слоями, в рассмотрение входят
только нормальные силы (давления). Касательные же силы (напряжения
сдвига) в такой модели отсутствуют. Касательные силы связаны с наличием у жидкости вязкости. В идеальной жидкости также отсутствует трение между жидкостью и твердой стенкой и вводится понятие скольжения
жидкости вдоль стенки. В действительных жидкостях действуют как нормальные, так и касательные силы, возникающие между слоями жидкости
и на границе жидкости и поверхности обтекаемого ею тела, а также между
жидкостью и твердой стенкой, что приводит к прилипанию жидкости к
стенке [1].
На сегодняшний день, все теории, описывающие турбулентность, являются полуэмпирическими, то есть для их построения используются данные, полученные в результате экспериментов. Подобно описанной ситуации, когда в рамках одного раздела науки имеется дело с почти не пересекающимися между собой путями исследования, в современной гидродинамике раздельно изучаются две условные группы турбулентных течений.
Первая – это полуэмпирические теории турбулентности (Прандтль, Карман, и др.), в которых устанавливаются алгебраические соотношения между касательными напряжениями и градиентами осредненных скоростей.
Вторая группа – это интенсивно развивающиеся в настоящее время мето5
ды численного интегрирования средствами ЭВМ уравнений Рейнольдса с
помощью различных моделей турбулентности, которые носят чисто эмпирический характер.
Схематически, такой метод можно представить как совокупность уравнений Рейнольдса и модели турбулентности, благодаря которой замыкается
полученная система уравнений. В сравнении с полуэмпирическими теориями из первой группы, численные методы дают возможность исследовать
гораздо более широкий спектр турбулентных течений, в первую очередь
в связи с наличием большего числа эмпирических констант, необходимых
для реализации метода на ЭВМ. Помимо этого, необходимо также задавать
и «пристенные функции» - в том случае, если производится расчет вблизи
границы с твердой стенкой.
Ламинарный и турбулентный режимы течений в современной гидродинамике также рассматриваются раздельно. Следствием данного подхода
является возникновение обособленных друг от друга и разных по своему содержанию теорий этих режимов течения. В рамках каждого из них имеют
место свои реологические соотношения, которые приводят к разным выражениям ряда гидродинамических характеристик, среди которых профили
скоростей и переменные, описывающие рассматриваемое течение в рамках
выбранной модели. Это иллюстрирует тот факт, что задачи о рассмотрении
перехода потока из ламинарного режима течения в турбулентный нетривиальны и по сей день требуют детального рассмотрения [36].
Определяющим фактором при выборе математической модели течения является критическое число Рейнольдса Reкр , которое характеризует
переход от ламинарного режима течения к турбулентному и зависит от
типа геометрии течения. К примеру, для круглых труб Reкр ≈ 2300. В
некоторых случаях, если течение не достаточно исследовано экспериментально (например, при проектировании нового оборудования), становится
невозможным точно указать Reкр и в результате невозможно достоверно утверждать, с каким режимом течения предстоит иметь дело. Между
тем, использование математической модели течения, несоответствующей
реальной картине, то есть являющейся результатом изначально неверного
предположения - естественно - самым критическим образом скажется на
достоверности результатов подобных расчетов.
6
Существующие модели турбулентности нельзя считать универсальными, т.к. точность и достоверность результатов, полученных при их применении к различным видам турбулентных течений, будут различаться.
Это означает, что исследования в области гидродинамики турбулентных
течений далеки от завершения. Сама по себе проблема описания турбулентности является многоплановой и нетривиальной и, как упоминалось,
не решена в рамках какого-то одного подхода.
Решить эту проблему на чисто феноменологическом уровне можно
с помощью феноменологической модели, разработанной В. А. Павловским
[17]. Данная модель является своего рода альтернативой гипотезе длины
пути перемешивания Прандтля и получила название f − модели. Для исследования f − модели на предмет ее работоспособности, в настоящей работе будет рассматриваться внутреннее течение Пуазейля, ввиду широкой
распространенности данного вида течения.
Цели исследования
• Изучение природы турбулентности и актуальных методов ее исследования.
• Рассмотрение течения Пуазейля в гладкой цилиндрической трубе.
• Применение f - модели турбулентности для расчета внутренних течений.
• Формулировка математической модели рассматриваемого течения.
• Разработка и верификация алгоритма численного решения данной
задачи.
• Получение численных результатов для различных чисел Рейнольдса.
• Получение аналитического решения задачи для различных чисел Рейольдса.
• Сопоставление численных результатов с аналитическим решением и
имеющимся экспериментальными данными.
7
1
Современное состояние теории турбулентности
1.1
Актуальные проблемы современной гидродинамики
При рассмотрении течений различных видов, возникает фундаментальный вопрос: при каком условии течения каких-либо жидкостей около
геометрически подобных тел при одинаковых направлениях натекания будут геометрически подобны картины линий тока этих течений? Согласно
[2], таким критерием механического подобия двух течений является число
Рейнольдса Re.
Каждому типу геометрии течения соответствует также свое критическое число Рейнольдса Reкр , которое характеризует переход течения потока жидкости из ламинарного режима в турбулентный. Reкр является определяющим фактором при выборе модели течения жидкости (для круглой
трубы, к примеру, Reкр ≈ 2300). Следует также отметить тот факт, что в
определенных случаях (к примеру, при разработке нового устройства было
получено новое, мало исследованное экспериментально течение), указать
значение Reкр становится проблематично. Как упоминалось во введении,
такая неясность приводит к невозможности дать точный ответ на вопрос
о выборе той или иной модели течения.
Применение моделей турбулентности на практике приводит к значительному упрощению использование специальных вычислительных комплексов, содержащих широкий перечень математических моделей, описывающих течения жидкости и предназначенных для расчетов разнообразных течений [3]. Однако, всякий раз, когда рассматривается определенное
турбулентное течение, возникает также и проблема в выборе модели турбулентности. Эта проблема обусловлена тем, что каждая такая модель ориентирована на некоторый класс течений, для которого её использование является особенно эффективным. На выбор модели турбулентности зачастую
оказывают влияние и другие параметры, такие как необходимая точность
ожидаемого результата и имеющиеся ресурсы: как вычислительные , так и
временные [4].
8
Существующие модели турбулентности, как уже говорилось, нельзя
считать универсальными, поскольку их применение к различным классам
турбулентных течений не позволяет получать результаты одинаково высокой степени точности и достоверности. Это означает, что научные исследования в области гидродинамики турбулентных течений, в том числе и в
рамках феноменологического подхода, по-прежнему остаются незаконченными.
Проблема описания турбулентности является нетривиальной и многоступенчатой и до настоящего времени остается неразрешимой в рамках
какого-то одного подхода. С другой стороны, разносторонний подход к ее
рассмотрению и использование различных методов исследований дает надежду на значительное продвижение. B связи с этим, полезным может оказаться такой новый подход к описанию течений вязкой жидкости, который
был бы основан на использовании реологического соотношения, являющегося независимым от режима течения, и удовлетворительно, в согласии с
опытом, описывал бы поле осредненных скоростей и трение на стенках для
каждого конкретного числа Рейнольдса [37,38]. В настоящее время предпринимаются попытки разработать подобную модель.
Достичь этой цели можно также и на чисто феноменологическом
уровне, проводя анализ обширных экспериментальных данных о течениях
вязкой жидкости, по большей части в турбулентном режиме - для ламинарного режима теория и опыт хорошо согласуются. Основы такой феноменологической модели, являющейся в некотором смысле чисто феноменологической альтернативой гипотезе длины пути перемешивания Л. Прандтля,
содержатся в работе В. А. Павловского [5]. За основу данного исследования
взята именно эта феноменологическая теория течения вязкой жидкости.
Эта модель получила название «f − модель» турбулентности, благодаря
феноменологической величине f , характеризующей меру турбулентности.
Теория нуждается в дальнейшем развитии и распространении ее на различные виды течений.
Турбулентный режим течения чрезвычайно сложен для теоретического описания, поскольку поля физических величин, гидродинамических
характеристик потока, изменяются случайным образом - как в пространстве, так и во времени. При этом частота этих изменений весьма высока
9
[6]. Из-за этого проблема об описании турбулентности до сих пор не нашла
своего окончательного решения, поскольку невозможно описать поведение
рассматриваемой величины (скорость, плотность, давление и др.) во всех
деталях, как функцию, зависящую лишь от координат и, если поток не
стационарен, времени. Тем не менее, достоверная теория для данного типа течения, в силу его распространенности на практике, является весьма
востребованной.
Наблюдения показывают, что вместе с хаотичностью имеет место
некоторая регулярность картины течения во времени и пространстве. А
именно, в некоторой фиксированной точке потока регулярно повторяется с течением времени одна и та же картина течения. С другой стороны,
в некоторой фиксированный момент времени одна и та же картина течения регулярно повторяется в пространстве. Эти факты свидетельствуют
о том, что теоретическое описание турбулентных потоков возможно. При
рассмотрении турбулентности могут быть выделены статистически точные
осредненные значения рассматриваемых величин и их корреляций (осредненных повремени и пространству). Это позволяет изучать физические
закономерности, характеризующие ее. В настоящее время в теории турбулентных течений развивается несколько различных направлений [5]. Поскольку пульсации гидродинамических характеристик представляют собой
случайные поля-процессы (их неупорядоченное изменение происходит как
в пространстве, так и во времени), то наиболее естественным математическим аппаратом для их описания является теория случайных полей в
самой общей постановке. Однако, несмотря на то, что данный подход развивается уже более полувека, в теоретическом плане его успехи ограничиваются в основном исключительно однородной и изотропной турбулентностью, попытки же развития статистического подхода на наиболее интересные с точки зрения практики неоднородные сдвиговые течения оказались
безуспешными. В связи с развитием вычислительной техники и совершенствованием вычислительных алгоритмов и процедур, все большее внимание
уделяется расчету турбулентных течений исходя непосредственно из уравнений Навье - Стокса - метод прямого численного моделирования (Direct
Numerical Simulation, DNS) [8]. Недостатками данного подхода являются
затраты рабочего и машинного времени, и необходимость в каждом кон10
кретном случае тщательной проверки адекватности используемых приближенных уравнений [39-41]. В последнее время особое внимание уделялось
системам нелинейных неконсервативных дифференциальных уравнений с
детерминированными параметрами, приводящими (в некотором интервале значений этих параметров) к столь неупорядоченным решениям задач,
что описываемые ими поля величин уместно рассматривать как случайные.
Такого рода решения весьма неустойчивы, сильно реагируют на малые изменения начальных и граничных условий. Но это касается только самих
неизвестных, а не их осредненных значений, которые, как правило, весьма
устойчивы.
Интересны так же попытки трактовать турбулентные течения через
призму уравнений больцмановского типа. В обоих случаях кинетические
уравнения понимаются не буквально, а как дальнейшее развитие представлений Рейнольдса и Прандтля о том, что турбулентные течения могут
рассматриваться как квазимолекулярные движения элементарных объемов
жидкостей.
1.2
Вычислительная гидродинамика
В связи с ростом мощности компьютеров и совершенствованием вычислительных алгоритмов стремительно развивается вычислительная гидродинамика (ВГД), в которой также существует несколько подходов к вопросу о моделированиия турбулентных течений. Эти подходы базируются
либо на использовании уравнений Навье-Стокса в исходной форме, либо
на их модификациях [9] с аппроксимацией членов, из-за чего из рассмотрения исключаются мелкомасштабные турбулентные пульсации. Основными методами являются прямое численное моделирование (Direct Numerical
Simulation - DNS), метод крупных вихрей (Large Eddy Simulation - LES)
и решение осредненных по Рейнольдсу уравнений Навье-Стокса (Reynolds
Averaged Navier-Stokes - RANS) [8]. Существуют также и различные комбинированные подходы, содержащие в себе черты выше перечисленных
подходов: RANS, LES и DNS. Методы прямого численного моделирования и крупных вихрей позволяют получать более качественную картину
турбулентности, чем RANS, однако их использование связано с большими
вычислительными затратами, особенно это касается DNS. Поэтому самым
11
распространенным на сегодняшний день в вычислительной гадродинамике
является RANS [42].
Наибольшие успехи прикладной теории турбулентности связаны с
феноменологическим подходом [7, 10, 11-12]. Такой подход основан на преобразованных уравнениях Навье-Стокса. Изменение масштаба осреднения
в них (осреднение по Рейнольдсу), превращает их в уравнения течения особой сплошной среды — «жидкости Рейнольдса». Осредненные по Рейнольдсу уравнения Навье-Стокса дают возможность исследовать турбулентные
течения, опираясь не на мгновенные пульсационные значения скорости и
давления, а на их статистически точные осредненные характеристики. Однако система, состоящая из уравнения течения «жидкости Рейнольдса» и
уравнения неразрывности даже для несжимаемой жидкости не является
замкнутой, ввиду чего необходимы дополнительные определяющие соотношения — реологические уравнения состояния. Основополагающими при
построении этих соотношений являются гипотезы турбулентной вязкости
(Буссинеск), пути перемешивания (Прандтль) и самоподобия стационарных турбулентных течений (Карман) [47]. Исходя из них, был выведен
универсальный логарифмический профиль скорости, являющийся ближней асимптотикой для большинства стационарных плоских турбулентных
потоков. Применительно к течению в трубах перечисленные гипотезы привели к универсальному закону сопротивления и формуле для дефекта скорости, охватывающего практически весь экспериментально доступный интервал изменения числа Рейнольдса [48].
Эмпирические модели, использующие понятие турбулентной вязкости, можно классифицировать следующим образом [13]:
1 алгебраические модели
2 модели с одним дифференциальным уравнением переноса характеристик турбулентности (однопараметрические модели)
3 модели с двумя дифференциальными уравнениями переноса характеристик турбулентности (двухпараметрические модели)
4 модели с большим числом уравнений.
12
Алгебраические модели связывают турбулентную вязкость µT с градиентом поля скоростей. Они ведут свое начало от работ Л. Прандтля
[8], который ввел понятие длины пути перемешивания l, что позволило
записать
турбулентную
вязкость для простых сдвиговых течений в виде
. Эта формула позволила решать ряд задач теории турбулентных течений, важных для инженерных приложений. Однако эту формулу
невозможно распространить на трехмерный случай. Также к недостаткам
алгебраических моделей следует отнести их ограниченную область применения, т.е. привязкой к конкретному типу течения, предположение гипотезы пути перемешивания о состоянии локального равновесия турбулентности, т.е. о соблюдении баланса между генерацией турбулентной энергии и ее
диссипацией в каждой точке расчетной сетки. Также, в соответствии с данной теорией при достижении скоростью локального максимума, турбулентные напряжения исчезают, что, вообще говоря, не верно. Учет подобных
фактов для подобных моделей приводит к их усложнению, что лишает их
своей привлекательности. В настоящее время данные модели применяются
в основном для расчета течений в непосредственной близости от стенки
[4]. Попытки устранения недостатков алгебраических моделей привели к
созданию моделей, содержащих одно балансовое дифференциальное уравнение переноса турбулентности, которое выводится либо для турбулентной
вязкости, либо для кинетической энергии турбулентности [14].
Модели турбулентности, содержащие два дифференциальных уравнения переноса, берут свое начало в работах А. Н. Колмогорова [15], они
позволяют рассчитывать двумерные и трехмерные течения и в настоящее
время являются наиболее распространенными. Для данного класса моделей предполагается решение дифференциальных уравнений переноса для
двух характеристик турбулентности, которые позволяют определить турбулентную вязкость посредством алгебраических соотношений. Среди моделей, содержащих два дифференциальных уравнения переноса, в настоящее время наиболее распространены модели k − ε и k − ω, а также их
различные модификации [43]. В этих моделях в качестве первой характеристики турбулентности используется величина кинетической энергии
турбулентных пульсаций k. В k − ε модели второй характеристикой является диссипация энергии турбулентности ε, а в k − ω — удельная скорость
13
диссипации, ω ∼ kε .
В зависимости от модели, для вихревой вязкости справедливы зави2
симости: k − ε − νT ∼ kε , k − ω − νT ∼ ωk
В частности, в моделях группы k − ε турбулентная вязкость выражается как
k2
µT = ρCµ ,
ε
где Cµ - эмпирическая константа; k - кинетическая энергия турбулентности;
ε - диссипация энергии турбулентности.
Для расчета величин k и ε используются соответствующие уравнения
переноса, которые можно получить трансформацией векторного уравнения
Навье-Стокса
1.3
Турбулентность в ANSYS
Прежде чем перейти к расмотрению задачи о турбулентном течении
в рамках представленной В. А. Павловским f − модели, рассмотрим еще
один важный вопрос: возможность численного решениия задач, связанных
с турбулентными течениями жидкостей, при помощи имеющихся на сегодняшний день пакетов прикладного программного обеспечения. Лидером
в этой области является линейка продуктов корпорации ANSYS Inc.
ANSYS является универсальным и наиболее распространенным пакетом программного обеспечения, основывающимся на методе конечно элементного (МКЭ) анализа. С момента его появления около 30 лет назад, он успешно развивающаяся и в настоящее время. ANSYS широко используется специалистами в сфере автоматизированных инженерных расчётов (Computer-Aided Engineering , CAE) и конечно - элементного решения практически всех типов пространственных задач механики деформируемого твёрдого тела и механики конструкций, задач механики жидкости и газа, теплопередачи и теплообмена, электродинамики, акустики, а
также механики связанных полей. Моделирование и анализ, как альтернатива дорогостоящим и порой трудноосуществимым экспериментам типа
«проектирование — изготовление — испытания» делает его по-настоящему
незаменимым инструментом во многих производственных сферах. ANSYS
предоставляет наиболее широкий спектр современных физических моде14
лей. Также имеется возможность проводить мультидисциплинарные расчеты в рамках интегрирующей среды ANSYS Workbench.
ANSYS DesignModeler - программный модуль, разработанный для создания и подготовки расчетной модели. Он прост в использовании, построен на полностью параметрической среде; присутствует возможность прямой двусторонней связи с ведущими CAD пакетами.
ANSYS CFD распологает наиболее полным набором моделей турбулентности, среди которых широко используемые k − ε, k − ω и модель Рейнольдсовых напряжений (RSM) для сильно закрученных или анизотропных течений. Благодаря наблюдаемой в последнее время тенденции увеличения производительности и снижения стоимости ЭВМ [16], модели крупных вихрей (LES) и более экономичные модели неприсоединенного вихря
(DES) становятся все более популярны при решении производственных задач. Новые модели позволяющие рассчитывать ламинарно-турбулентный
переход и современная модель (SAS), автоматически определяющая масштаб турбулентных вихрей, обычно применяются в случаях, когда точности стационарных моделей турбулентности недостаточно. Различные пристеночные функции и метод enhanced wall treatment (EWT) позволяют максимально точно описать потоки, ограниченные стенками.
ANSYS CFX
Программный комплекс, сочетающий возможности анализа гидрогазодинамических процессов, многофазных потоков, химической кинетики,
горения и многих других. CFX обеспечивает принципиально новый уровень решения задач вычислительной гидрогазодинамики за счет уникального сочетания технологий, начиная от прямого интерфейса к большинству CAD систем и заканчивая возможностью проводить сопряженный
жидкостно-структурный анализ совместно с ANSYS Multiphysics. Широкий выбор моделей турбулентности, в сочетании с линейным решателем
с технологией "Algebraic Coupled Multigrid"позволяет добиться высокой
точности результатов при решении различного класса задач. CFX имеет
возможность параллельного счета на любых параллельных вычислительных системах [46] (с распределенной (SMPD) и общей (SIMD) памятью)
под различными ОС (включая Windows, LINUX, UNIX). CFX обладает
15
хорошей масштабируемостью (до 64 процессоров).
FLUENT
FLUENT - программный комплекс, предназначенный для решения
задач механики жидкостей и газов. С 2006 года FLUENT входит в состав
ANSYS.
FLUENT использует неструктурированную сеточную технологию (типы элементов - гексаэдры, тетраэдры, призмы и пирамиды). Адаптация
расчетной сетки позволяет получить точное решение для областей с большими градиентами потока. Это наиболее важно при пограничных слоев.
Возможности использовать модели динамической адаптивной сетки позволяют проводить такие расчеты, как: потоки в цилиндрах, клапаны и
другие. Во FLUENT включены ламинарные и турбулентные модели гидродинамики, теплопередачи, фазовых переходов и радиации, а также модели
для расчета кавитации, сжимаемых жидкостей, теплообмена, теплопроводности, реальных газов, и модуль для расчета влажного пара.
FLUENT включает в себя множество моделей турбулентности: несколько версий модели k − ε, модели k − ω, Reynolds stress модель (RSM), LES
и DES.
Постпроцессор FLUENT обладает возможностью отображать результаты расчета в векторном и контурном видах, а так же отображать траектории движения частиц. Также имеется возможность анимирования результатов расчетов. FLUENT позволяет выполнять параллельные вычисления на Windows, Linux, и Unix платформах. При этом могут быть использованы многопроцессорные машины или кластеры. FLUENT обладает
клиент-серверной архитектурой и довольно хорошей масштабируемостью.
1.4
f − модель турбулентности
Как было отмечено ранее, модель течения жидкости, которую можно
было бы применить в расчетах, не задаваясь вопросом о режиме течения,
востребована во многих исследовательских областях.
Гипотеза длины пути перемешивания Л. Прандтля на полуэмпирическом уровне описывает пристенные турбулентные течения и тем самым
16
удовлетворяет потребностив инженерных расчетах таких течений. «Теория
приводит к логарифмическим профилям осредненных скоростей, хорошо
согласующимися с экспериментальными данными за счет надлежащего выбора двух констант. Следует заметить, что она не применима в непосредственной близости к стенке, поскольку краевое условие прилипания в рамках её невыполнимо. Поэтому для придания формуле Л. Прандтля физического смысла приходится вводить понятие вязкого подслоя. В настоящее
время используются двух- и трехслойные модели течения, позволяющие
сопрягать зону потока с чисто молекулярным трением с областью развитого турбулентного течения, где трение имеет молярную природу, и профиль
скорости является логарифмическим. Естественно, что при использовании
таких многослойных моделей число эмпирических констант увеличивается» [17].
Можно пойти иным путем и построить чисто феноменологическую
теорию, альтернативную гипотезе длины пути перемешивания Прандтля,
которая также приводит к логарифмическим профилям осредненных скоростей, обеспечивает краевое условие прилипания и содержит всего лишь
две феноменологические константы.
Такая теория, а именно, дифференциальная модель течения несжимаемой вязкой жидкости чисто феномнологического характера − f − модель, являющаяся альтернативой теории длины пути перемешивания, предложена в работе [17] В. А. Павловским. В ней содержится безразмерная скалярная величина − весовая функция f , принимающая значения на интервале от нуля до единицы и уравнение переноса этой величины [18]. Данная
функция является отношением турбулентной взкости к суммарной. Модель
особенна тем, что дает предельные перходы: результаты решения совпадают с данными для ламинарного режима при малых числах Рейнольдса и
для турбулентного при больших. Также при больших значениях Re решения дают логарифмические профили скоростей и логарифмические законы
сопротивления, характерные для турбулентного режима течения. Неоспоримое преимущество данной модели в том, что она позволяет производить
расчеты течений не задумываясь, какой при этом режим течения реализуется − ответ на этот вопрос дают результаты расчета. Система уравнений
f − модели позволяет задать градиент трения вблизи стенки, обеспечивая
17
тем самым выполнение условия прилипания.
В следующих главах рассмотрена f − модель и на ее основе построено решение задачи о течении вязкой несжимаемой жидкости в круглой
цилиндрической трубе.
18
2
2.1
Ламинарное течение в круглой трубе
Течение Пуазейля в трубе
Будем рассматривать прямую круглую трубу с постоянным по всей
длине диаметром. Скорость течения на стенках трубы вследствие прилипания равна нулю, в середине трубы она имеет наибольшее значение. В точках цилиндрических поверхностей с осями, совпадающим с осью трубы,
скорость течения постоянна.[19] Отдельные концентрические слои скользят один по другому так, что скорость везде имеет осевое направление.
Движение такого вида называется ламинарным течением. На достаточно
большом расстоянии от входа в трубу распределение скоростей течения
вдоль радиуса не зависит от координаты в продольном направлении.
Движение жидкости в трубе происходит под действием перепада давления в направлении оси трубы, но в каждом поперечном сечении, перпендикулярном к оси трубы, давление можно рассматривать как постоянное
[2]. Движение каждого элемента жидкости ускоряется вследствие перепада
давления и замедляется вследствие напряжения сдвига, вызванного трением. Другие силы на жидкость не действуют.
Рассмотрим цилиндр длиной l и радиусом y внутри жидкости, ось
которого совпадает с осью трубы. В направлении этой оси на него действуют силы давления p1 πy 2 и p2 πy 2 , приложенные к входному и выходному основаниям цилиндра соответственно, а также касательная сила 2πylτ ,
действующая на боковую поверхность цилиндра. Приравнивая разность
сил давления (p1 − p2 )πy 2 касательной силе, получаем в качестве условия
равновесия в направлении движения уравнение:
τ=
p 1 − p2 y
l
2
(1)
В рассматриваемом случае скорость u уменьшается с увеличением координаты y и равна нулю при y = R. Поэтому на основании закона трения
du
τ = µ du
dy следует принять, что τ = −µ dy . Подставляя в (1), получим:
du
p1 − p2 y
=−
.
dy
µl 2
19
Интегрируя с начальным условием u(y) = 0 при y = R, получаем:
u(y) =
p1 − p2 2
(R − y 2 ).
4µl
(2)
Т.о. имеем параболическое распределение скоростей по радиусу трубы.
Наибольшее значение скорость имеет в середине трубы, где она равна
um =
p1 − p2 2
R.
4µl
Полное количество Q жидкости, протекающее сквозь сечение трубы (расход жидкости) определяется как объем параболоида вращения и следовательно,
π 2
πR4
Q = R um =
(p1 − p2 ).
(3)
2
8µl
Введя среднюю скорость течения через поперечное сечение трубы, равную
ū =
Q
,
πR2
запишем (3) в виде:
l
ū.
(4)
R2
Закон, выражаемый формулой (2) - закон Пуазейля ламинарного течения в
трубе. Такое течение может наблюдаться только в трубах со сравнительно
малым диаметром и при сравнительно небольших скоростях. С возрастанием скорости или при увеличении диаметра трубы характер течения в корне
меняется: перепад давления уже не зависит линейно от скорости течения
в соответствии с формулой (4). Кроме того, скорость течения распределяется по поперечному сечению трубы значительно равномернее и вместо
упорядоченного слоистого течения происходит течение, в котором на продольную скорость налагаются беспорядочные поперечные составляющие,
вызывающие сильное перемшивание жидкости в поперечном направлении.
Иными словами, жидкость переходит в турбулентный режим течения.
p1 − p2 = 8µ
20
2.2
Течение в начальном участке круглой трубы
Во входном поперечном сечении профиль скоростей имеет прямоугольную форму, но затем под воздействием трения он постепенно вытягивается и, наконец, на некотором расстоянии от входа в трубу, принимает
форму параболы. Л. Шиллер нашел приближенное теоретическое значение
длины начального участка. Однако измерения И. Никурадзе показали, что
на протяжении первой трети начального участка развитие профиля скоростей хорошо совпадает с расчетом Шиллера, но дальше переход измеренного профиля в параболический профиль Пуазейля переходит медленнее.
Приближенным решением этой задачи занимались и другие ученые.
Рис. 1: Распределение скоростей в начальном участке круглой трубы при
ламинарном течении.
По Толботу при числе Рейнольдса Re = 103 течение становится устоявшимся на расстоянии около сорока диаметров трубы от входного сечения,
что хорошо совпадает с результатами эксперимента [2].
2.3
Закон сопротивления для течеия в трубах
Для ламинарного течения связь между перепадом давления и количеством протекающей жидкости Q = πR2 ū определяется чисто теоретически
и при этом хорошо совпадает с опытом. Для турбулентного режима подобную зависимость можно установить только экспериментально, т. к. чисто
теоретический расчет турбулентных течений в настоящее время невозможен. [22] Связь между перепадом давления и расходом устанавливается
21
законом сопротивления для движения в трубе. В литературе выведено множество формул, определяющих сопротивление в трубе, причем более рание
из нх выведены без учета закона подобия Рейнольдса и зависят от выбора
единиц. В настоящее время таким фомулам придают безразмерный вид,
вводя безразмерный коэффициент сопротивления ξ, определяемый соотношением
ξ ρ ¯2
p1 − p 2
=
u,
L
d2
(5)
где d = 2R - диаметр трубы.
Согласно [2], наибольшее касательное напряжение имеет место у стенки трубы, где оно равно
p1 − p0 R
τ0 =
,
(6)
L 2
причем эта формула верна независимо от режима течения. Т.о., касательное напряжение на стенке может быть определено экспериментально, путем
измерения перепада давления.
Сравнивая равенства (5) и (6), получаем формулу
ξ
τ0 = ρu¯2 .
8
(7)
В 1911 г. Блаузиус получил эмпирическую формулу для коэффициента сопротивления в гладких трубах
ūd
ξ = 0, 3164
ν
− 14
,
известную как закон сопротивления Блаузиуса.
В этой формуле ūd/ν = Re - число Рейнольдса.Согласно этой формуле, безразмерный коэффициент сопротивления есть функция только числа
Рейнольдса. Закон Блаузиуса примним только для Re ≤ 105 . Отсюда видно, что в данном диапазоне чисел Рейнольдса перепад давления пропорционален ū7/4 .
И. Никурадзе произвел измерения сопротивления и распределения
скоростей для широкой области чисел Рейнольдса: 4 · 103 ≤ Re ≤ 3, 2 · 106 .
С увеличением числа Рейнольдса профили становятся все более пол22
ными (Рис.2).
Рис. 2: Распределение скоростей в гладкой трубе при различных Re по
Никурадзе. На оси ординат отложены безразмерные скорости u/U , а на
оси абсцисс - безразмерные расстояния y/R
Уравнения профиля скоростей можно записать в виде:
y 1/n
u
=
,
U
R
(8)
где U - максимальная скорость. В этом соотношении показатель степени
1/n слабо зависит от Re. Это наводит на предположение, что существуют
асимптотические законы, содержащие логарифм как предельное значение
очень малой степени. В сокраценой записи такой закон распределения скоростей имеет вид:
φ = A ln η + D,
(9)
где
A=
1
,
χ
1
D = − ln β
χ
свободные постоянные. При сравнении этого закона с экспериментальыми
23
данными Никурадзе, находим: A = 2, 5; D = 5, 5. Следовательно, χ =
0, 4; β = 0, 111. Таким образом, закон распределения скоростей в гладких
трубах при очень больших числах Рейольдса имеет вид
φ = 2, 5 ln η + 5, 5,
(10)
или
φ = 5, 75 lg η + 5, 5.
Из гипотезы пути перемешивания Прандтля в совокупности с гипотезой l = χy и гипотезы подобия Кармана, получаем [2] универсальный
закон распределения скоростей
y
U −u
=f
v∗
R
Явное выражение для функции f Ry можно получить, воспользовавшись (10):
U −u
R
= 2, 5 ln .
(11)
v∗
y
Логарифмический закон скоростей (17) допускает экстраполирование на произвольно большие числа Рейноьдса. Будем полагать, чт этим же
свойством обладает и подлежащий выводу закон сопротивления. Выполним вывод, следуя изложению Прандтля [26].
Имея распределение скоростей (11), вычислим среднюю скорость по
поперечному сечению:
ū = U − 3, 75v∗ .
(12)
Результаты экспериментов, выполненных Никурадзе, дают вместо
3,75 другое число, а именно
ū = U − 4, 07v∗ .
Из формулы (7) имеем
ξ=8
v 2
∗
ū
.
(13)
Далее, из универсального закона распределения скоростей (10) сле24
дует, что
Rv∗
U = v∗ 2, 5 ln
+ 5, 5 .
ν
Подставляя в (12), получим
ū = v∗
Rv∗
2, 5 ln
+ 1, 75 .
ν
(14)
Величину Rv∗ /ν, входящую в последнее равенство, можно представить как произведение числа Рейнольдса ūd/ν на выражение, зависящее
от ξ:
√
1 ūd v∗
ūd ξ
Rv∗
√ .
=
=
ν
2 ν ū
ν 4 2
(15)
Выразив теперь v∗ через ξ посредством соотношения (13) и подставля
в (14), получим закон сопротивления
1
ūd p
√ = 2, 035 lg
ξ − 0, 91.
ν
ξ
√
√
Отложив значения 1/ ξ как ординаты на абсциссах lg(Re ξ), получим графическое изображение закона сопротивления для гладких труб в
виде прямой.
Рис. 3: Универсальный закон распределения для гладкой трубы. Кривая(1)
соответствует закону Прандтля, (2) - закону Блаузиуса.
25
Как показывает рис.3, этот результат блестяще подтверждается измерениями: все экспериментальные значения укладываются на прямую 1 с
численными коэффициентами, лишь немного отличающимися от теоретических значений. Эта прямая имеет уравнение
1
ūd p
√ = 2, 01 lg
ξ − 0, 8.
ν
ξ
Закон, выражаемый этим уравнением, называется универсальным законом сопротивления Прандтля для гладких труб. Он проверен И. Никурадзе [21] вплоть до Re = 3, 4 · 106 и дает хорошее совпадение с измерениями. В соответствие с условиями вывода, этот закон может быть экстраполирован до произвольно больших чисел Рейнольдса.
26
3
Феноменологическая f − модель турбулентности
3.1
Физический смысл функции f
Будем рассматривать дифференциальную феноменологическую модель течения вязкой несжимаемой жидкости - f -модель, в определенном
смысле альтернативную теории длины пути перемешивания. В ней содержится безразмерная скалярная величина - некоторая весовая функция f ∈
(0, 1) и уравнение переноса этой величины. Эта функция представляет собой отношение турбулентной вязкости к суммарной, включающей в себя и
молекулярную вязкость.
В турбулентном режиме молекулярная вязкость мало влияет на касательные напряжения. Тензор Рейнольдсовых напряжений можно выразить
через тензор суммарных напряжений σΣ , который включает в себя как турбулентные, так и чисто вязкостные напряжения.
Согласно f - модели, тензор Рейнольдсовых напряжений имеет вид:
R = f σΣ ,
где f - некоторая безразмерная величина, равная нулю на стенке и близкая
к единице в области развитого турбулентного течения.
Смысл величины f легче всего продемонстрировать на примере простого сдвигового течения, когда имеется лишь одна компонента скорости
− продольная, которая зависит только от поперечной координаты. Для
такого течения:
τ
f = T,
τΣ
причем
τΣ = τL + τT ,
где τΣ , τL , τT - соответственно, суммарное, вязкостное и турбулентное
касательные напряжения.
Как функция от перемещения и критерия Рейнольдса, f (y, Re) обладает следующими свойствами:
при Re → 0, f → 0, при Re → ∞, f → 1.
27
Так, под функцией f можно понимать скалярную безразмерную меру
режима течения потока в рассматриваемой точке.
В случае простого сдвигового течения:
τL = µ du
dy – вязкое напряжение длёя сдвиговых течений (Ньютон),
τT = µT du
dy турбулентное напряжение (Буссинеск),
где µT - турбулентная вязкость.
Здесь и далее, под скоростью u будем понимать осредненное по Рейнольдсу значение скорости. При ламинарном режиме осредненное и мгновенное значения скорости совпадают.
Имеем:
µT
,
f=
µ + µT
откуда видно, что f показывает соотношение молекулярной и турбулентной
долей вязкости в потоке. Отсюда турбулентная вязкость:
µT = µ
f
1−f
и в результате, после подстановки, получаем формулу для суммарного напряжения
µ du
.
(16)
τΣ =
1 − f dy
На рис.4 представлен характерный вид функции f по сечению канала: на стенке она равна нулю, а на оси достигает максимальное значение.
Рис. 4: Поведение функции f по сечению канала.
а) характерный вид; б) простейшая аппроксимация
28
Далее рассмотрим установившееся турбулентное течение несжимаемой вязкой жидкости в плоском канале шириной 2h. В проекуии на продольную координату, уравение движения сплошной среды в напряжениях
примет вид:
dp dτ
− + Σ = 0,
dx
dy
или
dτ
dp dτ
− + L + T = 0.
dx
dy
dy
Учитывая для τL закон вязкого трения Ньютона, имеем:
dp
d2 u dτT
− + ρν 2 +
.
dx
dy
dy
Здесь ρ− плотность, p− давление, ν− кинематическая вязкость, u− осредненная продольная скорость, τL = −ρu0 v 0 − чисто турбулетное касательное
напряжение, x и y − соответственно, продольная и поперечная координаты, причем далее координату y будем отсчитывать от стенки.
Интегрирование этого уравнения дает:
−
dp
du
(y − h) + ρν
+ τT = 0,
dx
dy
или, если ввести ввести динамическую скорость
s
,
dx ρ
этот интеграл примет вид:
du
y τT
2
ν
= v∗ 1 −
− .
dy
h
ρ
Так как для данного течения τΣ = ρv∗2 (1 − y/h), то, учитывая f = τT /τΣ ,
запишем это выражение в безразмерном виде:
dv
= Re∗ (1 − η)(1 − f ),
dη
(17)
где η = y/h, Re∗ − число Рейнольдса, вычисленное по динамической
29
скорости:
h
u
Re∗ = v∗ , v = .
ν
v∗
Далее надо определить f так, чтобы она наилучшим образом описывала экспериментальные данные, в частности, скоростей и сопротивлений. Прежде всего можно аппроксимировать f простейшим образом, в виде
кусочно-линейной функции так, как это показано на рис. 5б и тем самым
построить самый простой вариант теории. Задаваясь значениями α и β
можно рассчитать профиль и сопротивление. Эти значения подбираются
так, чтобы обеспечивалось наилучшее согласование с опытными данными.
Очевидно, они зависят от Re∗ . Анализ показывает, что для Re∗ = 102 величины α ≈ 0, 95, β ≈ 0, 80; для Re∗ = 103 α ≈ 0, 70, β ≈ 0, 999, для
Re∗ = 104 α ≈ 0, 49, β ≈ 0, 9999, т.е. β быстро стремится к единице.
Для того, чтобы такое отличие было заметным, от величины f полезно перейти к безразмерной турбулентной вязкости νtv = νt /(v∗ · h). Связь f
и νtv определяется формулой:
νtv
Re∗ · νtv
f
,f =
=
Re∗ · (1 − f )
1 + Re∗ · νtv
(18)
В настоящее время накоплено достаточно достоверных данных относительно величины νtv . Схематично они изображены на графике рис.5a
[22].
Рис. 5: Поведение безразмерной турбулентной вязкости по сечению потока:
а) – данные опытов, б) – простейшая аппроксимация.
Её можно аппроксимировать в виде кусочно-линейной функции, представленной на рис.5б следующим образом:
30
(
νtv
=
νtv0 η/η0 , η ≤ η0
νtv0 , η0 ≤ η ≤ 1
(19)
Тогда, учитывая соотношение (18), формулу (17) можно записать в
виде
Re∗ · (1 − η)
dv
=
dη
1 + Re∗ · νtv
(20)
Интегрирование этого уравнения с учётом (19) даёт для пристенной
области, где η ≤ η0 :
η0
v η
+
v = v ln 1 + Re∗ νt0
νt0
η0
1
+
Re∗
η0
νtv0
2
η
η
− Re∗ νtv0
(21)
ln 1 + Re∗ νtv0
η0
η0
а для зоны, где η0 ≤ η ≤ 1:
Re∗
v = v(η0 ) +
1 + Re∗ νtv0
η2
η02
η−
− η0 +
2
2
(22)
Если положить νtv0 = const = 0, 08, что неплохо согласуется с данными [16], то величина η0 будет функцией числа Рейнольдса Re∗ . Анализ обработки экспериментальных данных даёт следующую эмпирическую формулу для η0 , справедливую для турбулентного режима течения (примерно
для η0 > 102 ):
η0 = 0, 25 + 16/[lg(1 + Re∗ )]4
(23)
По сопротивлению этот вариант с νtv0 = 0, 08 и η0 по формуле (23)
приводит к удовлетворительному согласованию с опытом, но профили скоростей вблизи стенки оказываются недостаточно полными. Поэтому с целью устранения этого недостатка, а также для обеспечения предельного
перехода к ламинарному режиму течения величину νtv0 также следует считать функцией от Re∗ . В результате обработки опытных данных можно
получить следующие эмпирические формулы:
31
0, 71Re4∗ − 250Re3∗ + 8, 6 · 107 · Re2∗
η0 = 1 −
(3100 + Re∗ )4,018
νtv0 =
0, 117Re4∗
515Re3∗
4
+
+ 9, 5 · 10 ·
(3100 + Re∗ )4,018
Re2∗
(24)
На графике рис.6 показаны профили скоростей, вычисленные по формулам (21)-(24) на фоне известного логарифмического закона Л. Прандтля.
Соответствующий закон сопротивления приведён на рис.1.4 (число
Рейнольдса здесь Re∗ = ucp h/ν) на фоне асимптотик согласно [34].
Рис. 6: Профили скоростей в универсальных координатах для различных
Re∗ : 1 - Re∗ = 103 ; 2 - Re∗ = 104 ; 3 - Re∗ = 105 ; пунктир соответствует
закону u/v∗ = 2, 5ln(yv∗ /ν) + 5, 5.
Изложенная «скелетная» теория приводит к логарифмическому профилю в зоне развитого турбулентного течения. При этом она даёт возможность описания течения во всём диапазоне чисел Рейнольдса, т.е. она
охватывает как ламинарный, так и турбулентный режимы течения. Переходная зона между ними с помощью формул (24) не описывается.
32
Рис. 7: Закон сопротивления в плоском канале; сплошная линия – расчёт,
штриховые соответствуют ламинарному режиму и закону Блазиусовского
типа.
Анализ «скелетной» феноменологической модели показывает, что даже сравнительно грубое, кусочно-линейное представление функции f (или
связанной с ней турбулентной вязкостью νtv ) позволяет получать логарифмические профили и законы сопротивления, почти удовлетворительно согласующиеся с экспериментом. При этом, в отличие от длины пути перемешивания Л. Прандтля выполняется условие прилипания и, кроме того, нет
необходимости полагать постоянным касательное напряжение по сечению
потока.
Более тонкий вариант изложенной теории с использованием меры
турбулентности f можно получить, если эту величину представить не двумя, а тремя и более участками. В пределе можно получить непрерывную
кривую для функции f , которая подчиняется дифференциальному уравнению второго порядка.[31,32] При этом к тому же возникает возможность
записать уравнение движения жидкости в виде, справедливом и для ламинарного, и для турбулентного режима течения.
33
3.2
Система уравнений модели
Выражение для касательного напряжения (16), которое включает в
себя вязкие и турбулентные напряжения можно, как и формулу Ньютона,
обобщить на трёхмерный случай. Тогда тензор напряжений τ будет линейно связан с тензором скоростей деформаций (осреднённых, если поток
турбулентен):
τ=
где
2µ
S
1−f
(25)
T
1 ~
~ · ~v
S=
∇ · ~v + ∇
2
(26)
Если подставить реологическое соотношение (25) в уравнение движения сплошной среды в напряжениях
ρ
d~v
~ +∇
~ ·τ
= −∇p
dt
то получаем уравнение движения жидкости [29,31]:
ρ
d~v
µ
µ
~ +
~
= −∇p
∆~v +
d · ∇f
2
dt
(1 − f )
(1 − f )
(27)
В скалярном виде, в декартовой системе координат имеем три скалярных уравнения
∂ui
∂ui
ρ
+ uj
∂t
∂xj
∂p
µ
∂ 2 ui
=−
+
+
∂xi 1 − f ∂xk ∂xk
µ
+
(1 − f )2
∂ui ∂uk
+
∂xk
∂xi
∂f
(28)
∂xk
где i, j, k = 1, 2, 3.
Видно, что уравнение (27) представляет собой модификацию уравнения Навье-Стокса. При f = 0, при чисто ламинарном режиме, это уравение
переходит в уравнение Навье-Стокса.
Для замыкания системы уравнений, содержащей уравнение движения (27) и уравнение неразрывности, необходимо ещё одно скалярное урав34
нение – уравнение для функции f , которая образует некоторое поле f (~r, t).
Такое уравнение предложено в работе [24] – уравнение переноса
величины f :
~ × ∇f
~ ) · (∇
~ × ~v )
df
µψ(f ) ~ ~
(∇p
p
p
ρ = µ∆f +
∇f ∇f + ρ(1 − f )
dt
1−f
2|S : S| 2|Ω : Ω|
(29)
Здесь S – тензор скоростей деформаций, определяемый выражением
(26), Ω – антисимметричная часть тензора градиентов деформаций
T
1 ~
~ · ~v
∇ · ~v − ∇
Ω=
2
Двоеточие в выражении (29) обозначает двойное скалярное произведение тензоров, величина ψ(f )- алгебраическая функция f :
ψ(f ) =
2α + β(1 − f )
α + β(1 − f )
(30)
Величины α и β – две феноменологические константы:
α = 2, 5; β = 8, 5,
(31)
которые найдены из обработки опытных данных для широкого класса
пристенных турбулентных течений.
Видно, что уравнение (29) напоминает уравнение Бюргерса. Анализ его структуры показывает, что величина f изменяется в пределах
0 ≤ f < 1, что и следовало ожидать.
Система уравнений, состоящая из уравнения движения (27), уравне~ · ~v = 0 и уравнения (29) содержит три неизвестния неразрывности ∇
ные – скорость ~v , давление p и величину f , которая заключена в пределах
0 ≤ f < 1 согласно определению и структуре уравнения (29). Эта система
является замкнутой и позволяет решать задачи об изотермических течениях несжимаемой вязкой жидкости независимо от режима течения. [33]
Граничными условиями для этой системы уравнений будут условия прилипания и вязкого ньютоновского трения на твёрдой границе S с нормалью
~n, а также равенство нулю функции f на этой границе и условие на производную модуля скорости по этой функции:
35
v∗2
~v |S = 0, |2S · ~n|S = ,
ν
= (α + β)v∗ .
df S
(32)
В этих условиях v∗ – динамическая скорость, которая выражается
через касательное напряжение на стенке τW :
v∗ =
|τW |
ρ
(33)
Структура реологического соотношения (25) показывает, что тензоры
τ и S соосны, но связь между ними нелинейна при f > 0 (при f = 0, т.е. в
ламинарном режиме, эта связь линейна). [25] Тем самым соотношение (25)
может описывать и турбулентные течения, если пренебрегать анизотропией
и памятью турбулентного потока. Если же такие упрощения недопустимы,
что вполне может быть для потоков, в которых несоосность тензоров τ и S
существенна, то реологическое соотношение (25) нужно обобщать на этот
случай.
Однако для большинства турбулентных течений при их расчёте достаточно ограничиваться лишь учётом нелинейности, что подтверждается
успешным применением теории длины пути перемешивания Л. Прандтля.
В дальнейшем для описания течений вязкой несжимаемой жидкости
будем использовать следующую систему уравнений:
d~v
µ
µ
~ +
~
ρ = −∇p
∆~v +
S · ∇f
dt
(1 − f )
(1 − f )2
~ · ~v = 0
∇
~ × ∇f
~ ) · (∇
~ × ~v )
df
µψ(f ) ~ ~
(∇p
p
ρ = µ∆f +
∇f ∇f + ρ(1 − f ) p
dt
1−f
2|S : S| 2|Ω : Ω|
(34)
Здесь ~v – скорость; f = f (~r) – безразмерная скалярная величина,
мера турбулентности, 0 ≤ f < 1;µ – динамическая вязкость,p – давление;
d и w – тензор скоростей деформаций и спин, соответственно, характеризующие симметричную
и антисимметричную части тензора градиентов
~ v d + w = ∇~
~ v ; две точки обозначают двойное скалярное проскоростей ∇~
изведение тензоров; Ψ(f ) – скалярная функция от f :
36
ψ(f ) =
2α + β(1 − f )
α + β(1 − f )
где α и β – феноменологические константы; α = 2.5, β = 8.5.
Все величины в этих уравнениях понимаются как осредненные по
Рейнольдсу (в ламинарном режиме, очевидно, осредненные и мгновенные
значения величин совпадают). Краевые условия для этой системы уравнений – условия прилипания и вязкого ньютоновского трения на твердой
границе S с нормалью ~n, а также равенство нулю функции f на этой стенке
и условие на производную модуля скорости по этой функции:
v∗2
~
~v |S = VT B , |2d · ~n|S = ,
ν
= (α + β)v∗ .
df S
(35)
В этих выражениях ν – кинематическая вязкость, V~T B – скорость
твёрдой границы, v∗ – динамическая скорость, связанная с трением на стенp
ке τw , v∗ = |τw |/ρ, ρ – плотность жидкости.
3.3
Определение феноменологических констант
В любой феноменологической теории константы находятся из условия согласования решения задачи и данных опыта. Как известно, эталонной задачей для поиска значений этих констант является задача об установившемся течении жидкости в круглой цилиндрической трубе. Для этого течения имеются обширные экспериментальные данные, позволяющие
найти достоверные значения феноменологических констант. После нахождения таких значений теория обычно проверяется путём решения других
задач очередного сравнения с опытом. В необходимых случаях происходит
корректировка теории. Такой путь построения теории характерен для всех
феноменологических подходов [6].
Рассмотрим задачу о течении в трубе радиуса R. Здесь удобно ввести в рассмотрение безразмерную координату η, отсчитываемую от стенки
трубы:
η=
R−r
, 0 ≤ η ≤ 1,
R
37
(36)
где r – радиальная координата (отсчитывается от оси). Кроме того
удобно ввести и безразмерную скорость v = u/v∗ , где v∗ – динамическая
скорость, связанная с касательным напряжением на стенке τW , ρv∗2 = |τW |.
Оперируя с этими безразмерными величинами, можно компактно записывать решения задач как для ламинарного, так и для турбулентного режима
течения [34]. После нахождения профиля скорости в виде
v = v(η)
можно найти среднюю скорость
ucp
= vcp = 2
v∗
Z1
v(1 − η)dη
0
и далее коэффициент сопротивления
λ=
8
,
2
vcp
который будет выражаться будет выражаться через число Рейнольдса по
динамической скорости
Re∗ =
v∗ R
ν
С обычным числом Рейнольдса
Re =
2Rucp
ν
оно связано зависимостью
Re = 2Re∗ vcp
Интересно заметить, что из последнего выражения следует vcp = Re/(2Re∗ ),
что позволяет коэффициент сопротивления записать в виде
Re∗
Re
Это соотношение справедливо для любого режима течения, разница для
этих режимов будет в связи между Re∗ и Re, что определяется реологией.
λ = 32
38
В ламинарном режиме течения закон вязкого трения Ньютона
τ =ν
du
dr
и закон распределения касательных напряжений
τ = τW
r
r
= −ρv∗2
R
R
приводят к уравнению
du
r
= −ρv∗2
dr
R
После введения безразмерных величин η и v его можно записать в виде:
µ
dv
= Re∗ (1 − η)
dη
После интегрирования, с учётом условия прилипания η = 0, v = 0 можно
записать выражение для безразмерного профиля скоростей (течение Пуазейля):
v = Re∗ (η − η 2 /2)
(37)
Средняя по сечению скорость
Z1
vcp = 2
η2
Re∗
Re∗ η −
(1 − η)dη =
2
4
0
Число Рейнольдса
Re∗
Re2∗
Re = 2Re∗ vcp = 2Re∗
=
4
2
Коэффициент сопротивления
ξ=
8
128
64
=
=
2
vcp
Re2∗
Re
(38)
Аналогичные выкладки можно сделать для любой модели течения. В случае единой феноменологической ламинарно-турбулентной модели решение
задачи [34] приводит к профилю скорости, неявно зависящему от
39
координаты η
v = βf − αln(1 − f )
(39)
где функция f определяется из трансцендентного уравнения
Re∗ η − η 2 /2 =
αf
− βln(1 − f )
(1 − f )
(40)
При малых числах Рейнольдса это решение сводится к обычному течению
Пуазейля. После вычисления можно построить профиль в универсальных
координатах, (y+ , vu∗ ) где y+ = yvν+ = Re∗ η, а также найти vcp и далее число
Re и коэффициент сопротивления λ.
Далее полученное решение необходимо сравнить с экспериментальными данными. Для турбулентного режима течения в трубе, для течения
при больших числах Рейнольдса такими экспериментальными данными являются данные Никурадзе [27]. Они приводят в части сопротивления к
закону сопротивления Прандтля-Никурадзе:
p
1
√ = 2lg Re ξ − 0, 8
ξ
или, что практически то же самое: λ = 0, 0032 + 0, 221/Re0,237 .
В части профиля скоростей они дают для логарифмической зоны
(достаточно удалённой от стенки) закон
v = 5, 75 lg y+ + 5, 5
Условие прилипания для этого профиля не выполняется, что приводит к необходимости вводить в рассмотрение зону вязкого подслоя.
Однако можно предложить, на основе анализа экспериментальных
данных Никурадзе, следующую эмпирическую однослойную формулу для
турбулентного профиля скоростей, для которой выполняется условие прилипания:
v = 1, 25ln 1 + Re2∗ η 2 + arctg(Re∗ η)+
+
Re2∗ η 2
2
3
3,
93
+
5,
5η
−
4,
5η
(41)
1 + Re2∗ η 2
40
Эта формула справедлива для Re∗ > 300, т.е. для Re ≥ 104 . В таблице
1 приведены расчётные профили скоростей согласно формулы для двух
чисел Рейнольдса, для которых имеются опытные данные Никурадзе.
Таблица 1
η
Re = 2, 3 · 104 Re = 3 · 106
v/vmax
v/vmax
0
0
0
0,1
0,702
0,8
0,2
0,785
0,855
0,3
0,838
0,891
0,4
0,88
0,919
0,5
0,914
0,942
0,6
0,943
0,962
0,7
0,966
0,978
0,8
0,984
0,989
0,9
0,995
0,997
1
1
1
На графике рис. 8 показаны эти расчётные профили на фоне опытных
точек Никурадзе [36].
Видно, что формула (41) хорошо аппроксимирует профили скоростей.
В таблице 1.2 показаны результаты расчётов течений на основе эмпирической формулы (41) [30,35]. Для фиксированных значений числа Рейнольдса, вычисляемого по динамической скорости, приведены соответствующие
значения обычного числа Рейнольдса Re и коэффициента сопротивления.
41
Рис. 8: Профили скоростей для чисел Рейнольдса: 1 - Re = 2, 3·104 , расчёт,
согласно (41) 2 - Re = 3 · 106 , расчёт; точки – экспериментальные данные
Никурадзе [21].
Таблица 2
Re∗
Re
628,5
2, 3 · 104
4,361 0,0238
1,378
1000
3, 89 · 104 4,590 0,0211
1,324
5,041 0,0168
1,226
5, 04 · 105 5,703 0,0125
1,098
2523,7
104
1, 1 · 105
lgRe
ξ
lg (1000ξ)
51150
3, 0 · 106
6,790 0,0083
0,920
105
6, 2 · 106
6,790
0,083
0,920
106
7, 35 · 107 7,866 0,0059
0,772
На графике рис. 9 показан закон сопротивления (сплошная линия) на фоне
формулы Блазиуса (пунктирная линия):
ξ = 0, 3164/Re1/4
(42)
Любопытно заметить, что эта формула может быть аппроксимирована выражением ξ = 0, 1 · lg (100/Re)1/4 . Сплошная линия, соответствующая рас42
чёту по (41) практически точно следует формуле Прандтля-Никурадзе (3.3).
Погрешность составляет не более 4%.
Рис. 9: Закон сопротивления, соответствующий формуле (41); пунктирная
линия соответствует формуле Блазиуса (42).
Видно, что кривая сопротивления практически точно следует формуле Прандтля-Никурадзе.
Кроме аппроксимационной формулы (41) для турбулентных профилей скоростей, обобщающей экспериментальные данные и приводящие к
закону сопротивления Прандтля-Никурадзе, можно предложить и другую,
также хорошо аппроксимирующие опыты:
1
v = 1, 25ln 1 + 81Re2∗ η 2 + arctg(9Re∗ η)+
9
Re2∗ η 2
2
3
+
3,
93
+
5,
5η
−
4,
5η
(43)
1 + Re2∗ η 2
Результаты расчётов по этой формуле приведены в таблицах 3 и 4.
43
Таблица 3
Re = 2, 3 · 104 Re = 3 · 106
v/vmax
v/vmax
η
0
0
0
0,1
0,704
0,801
0,2
0,785
0,855
0,3
0,836
0,889
0,4
0,875
0,915
0,5
0,907
0,934
0,6
0,934
0,956
0,7
0,957
0,971
0,8
0,976
0,984
0,9
0,990
0,993
1,0
1,0
1,0
Таблица 4
Re∗
Re
lg Re
ξ
lg (1000ξ)
628,65
2, 3 · 104
4,361
0,0236
1,374
103
3, 91 · 104 4,592
0,0209
1,320
104
5, 06 · 105 5,704
0,0125
1,096
6,477 0,00926
0,967
105
6, 21 · 106 6,793 0,00829
0,918
106
5, 9 · 107
7,867 0,00590
0,770
51045
3, 0 · 106
Соответствующие профили скоростей и закон сопротивления показаны на
рис. 8 и 9. Видно, что и формула (43) хорошо аппроксимирует опытные
данные для турбулентного режима течения. Однако она, в отличие от (41),
более тщательно описывает ещё и зону течения вблизи оси.
Опираясь на формулу (37) и аппроксимационные формулы (41) и
(43) была реализована итерационная процедура для поиска значений и ,
фигурирующих в выражении (39). Она привела к следующим значениям
феноменологических констант
44
α = 2, 5
β = 8, 5
Эти значения будут использованы далее при решении задач на основе единой феноменологической ламинарно-турбулентной теории.
45
4
Расчет течения жидкости в круглой цилиндрической трубе на основе f − модели
4.1
Постановка задачи о течении в трубе
Итак, подходим к вопросу о стабилизированном изотермическом течении жидкости с постоянными физическими свойствами в цилиндрической круглой трубе.
Рис. 10: Модель течения в круглой цилиндрической трубе.
Согласно f -модели:
ρ
ρ
dv
µ
µ
= −∇p +
∆v +
S · ∇f,
dt
1−f
(1 − f )2
df
ψ(f )
(∇p × ∇f ) · (∇ × v)
p
= µ∆f + µ
∇f · ∇f + ρ(1 − f ) p
.
dt
1−f
2|S : S| 2|Ω : Ω|
(44)
(45)
В уравнении движения и уравнении переноса скалярной (безразмерной) функции f , 0 ≤ f < 1 величины:
)
ψ(f ) = 2α+β(1−f
, α = 2, 5; β = 8, 5;
α+β(1−f
)
h i
kg
ρ, m
- плотность; p, [P a] - давление;
3
µ, [P a · s] - динамическая вязкость;
v̄, ms - вектор скорости;
¯ − (∇v̄)
¯ T ) - тензор вращения;
Ω = 12 (∇v̄
d
d
∂
¯
dt - оператор материальной производной,
dt = ∂t + (v̄ · ∇);
¯ - вектор Гамильтона (набла).
f¯ - вектор интенсивности массовых сил, ∇
Уравнение неразрывности:
¯ · v̄ = 0.
∇
(46)
Задачу будем решать в цилиндрической системе координат (r, ϕ, z)
[29]. Для стабилизированного течения вектор скорости v̄ = v̄(vr , vϕ , vz )
46
имеет лишь оду компоненту - осевую vz , которая зависит лишь от одной
координаты - радиальной. Обозначим эту компоненту vz как u = u(r).
Скалярная величина f также зависит радиальной координаты, f = f (r).
Уравнение движения (44), пренебрегая массовой силой, принимает вид:
∂p 1 d
− +
∂z r dr
µ du
r
1 − f dr
Для стабилизированного течения p = p(z),
перейти к обыкновенным производным:
dp 1 d
−
+
dz r dr
∂p
∂z
µ du
r
1 − f dr
= 0.
=
dp
dz
= const, и тогда можно
= 0.
(47)
Уравнение (45) переноса функции f принимает вид:
q d
µ
r dr
2
dp df
df
ψ
df
dz dr
r
+µ
− (1 − f ) du = 0.
dr
1 − f dr
dr
(48)
Уравнение неразрывности (46) тождественно удовлетворяется:
0 ≡ 0.
Запишем уравнения в безразмерном виде, вводя следующие безразмерные переменные:
• безразмерную радиальную координату, отсчитываемую от стенки трубы.
r
η = 1 − ; dr = −Rdη.
R
На стенке, при r = R величина η = 0, а на оси, при r = 0, η = 1. Тем
самым, 0 ≤ η ≤ 1.
• безразмерную скорость
v=
u
,
v∗
где v∗ - динамическая скорость:
s
или
2
dp
= ρv∗2 .
dz
R
После введения этих переменных уравнения системы запишутся в виде:
d (1 − η) dv
= −2Re∗ (1 − η),
dη (1 − f ) dη
(49)
2
d
df
dv ψ(f )
df
df
dv 1
(1 − η)
+
+ 2Re∗ (1 − f ) = 0. (50)
dη (1 − η) dη
dη
dη (1 − f ) dη
dη
v∗ R
, ν - кинематическая вязкость. Безразмерный комν
ucp d
плекс Re связан с числом Рейнольдса Re =
:
ν
Здесь Re∗ =
Re =
ucp · 2R v∗ R 2ucp
=
·
= Re∗ · 2vcp ,
ν
ν
v∗
где vcp - безразмерная средняя скорость по сечению трубы,
Z1
(1 − η)vdη.
vcp = 2
0
Граничные условия уравнений (49),(50):
η=0:
v = 0,
dv
= Re∗ ,
dη
f = 0,
df
Re∗
=
.
dη
α+β
(51)
Система (49),(50) с граничными условиями (51) решается следующим
образом: задается значение Re∗ (10, 100, 1000, 10000), решается задача и
находятся профили скоростей v(η) и, кроме этого, f (η). Вычисляем vcp и
далее Re. Кроме этого, находим коэффициент сопротивления по формуле:
ξ=
8
2
vcp
.
48
4.2
Численный расчет и сравнение результатов
Рассмотрим результаты решения задачи посредством f − модели при
разных значениях Re∗ , соответствующих ламинарному и турбулентному
режимам. На рис.11-12 представлены профили безразмерной скорости и
меры турбулентности f .
Рис. 11: а. Re∗ = 101 ; б. Re∗ = 102 .
Рис. 12: а. Re∗ = 103 ; б. Re∗ = 104 .
Из графиков видно, что с ростом числа Рейнольдса профили скорости
увеличивают полноту. Это наглядно проиллюстрировано на рис.13:
49
Рис. 13: Профили скорости (а) и меры турбулентности (б) для чисел Рейнольдса.
На рис.14 показан универсальный закон распределения скоростей для
ламинарного и турбулентного режима течений в гладкой трубе. Результаты расчетов согласно f − модели сопоставляются с экспериментальными
данными.
Рис. 14: Универсальный закон распределения скоростей для гладкой трубы. Сплошная линия - расчет, маркерами нанесены данные [2].
50
Приведенные на графике выше данные позволяют сделать вывод о
хорошем согласовании расчетов с опытными данными.
Для буферной зоны результаты расчета и экспериментальные данные
Лауфера [23] показаы на рис.15.
Рис. 15: Буферная зона. Сплошная линия - расчет, пунктир - асимптотики,
маркеры - экспериментальные данные [23].
Результаты расчета коэффициента сопротивлений проиллюстрированы на рис.16. В качестве асимптотик используются законы сопротивления для ламинарного режима, закон Блаузиуса для турбулентного и закон
Прандтля-Никурадзе для гладких труб.
Рис. 16: Закон сопротивления для гладкой трубы.
51
На рисунке видно, что кривая сопротивления из ламинарного участка
переходит на участок, довольно точно соответствующий закоу ПрандтляНикурадзе.
На основании этих фактов можно говорить о хорошей согласованности расчетных данных с аналитическими, и, следовательно, о возможности
использования f − модели для решения задач о турбулентных течениях
жидкости.
4.3
Верификация на основе ламинарного течения
Для верификации выкладок выполним расчет при f = 0, то есть для
ламинарного режима течения. Вместо системы (49), (50) будет
d
dv
(1 − η)
= −2Re∗ (1 − η).
dη
dη
(52)
Неизвестная v(η) находится из решения уравнения с граничными
условиями:
η = 0 : v = 0;
Задается значение Re∗ , после нахождения профиля v(η) вычисляем
vcp и далее находим коэффициент сопротивления.
Итак, сначала интегрируем (52). Первый интеграл:
dv
η2
(1 − η) = −2Re∗ η −
+ C1 .
dη
2
Из условия η = 1,
dv
= 0 находим C1 = Re∗ . Отсюда имеем:
dη
dv
= Re∗ (1 − η).
dη
Интегрируем второй раз:
v = Re∗
η2
η−
2
+ C2 .
Из условия η = 0, v = 0 находим C2 = 0.
52
В итоге безразмерный профиль скоростей равен
v = Re∗
η2
η−
2
,
Средняя по сечению трубы скорость равна
Z1
Z1
(1 − η)v(η)dη = 2
vcp = 2
0
(1 − η)Re∗
dη =
Re∗
.
4
v
= 2(2η − η 2 ).
vcp
Коэффициент сопротивления равен ξ =
4.4
0
Отсюда получаем соотношение:
ξ=
η2
η−
2
8
128
=
, откуда имеем:
2
vcp
Re2∗
64
, что и следовало ожидать.
Re
Интегрирование задачи в рамках f − модели
f − модель позволяет получать для простых сдвиговых течений решения в квадратурах.
Общая схема решения:
1. задается значение Re∗
2. вычисляются профили скоростей v(η) и величина v∗ /U0
3. Находится Re и затем ξ
При построении численной процедуры системы нелинейных дифференциальных уравнений второго порядка, описывающие движение жидкости, сводились к системам уравнений первого порядка [44,45].
Обозначим y1 = v; y2 = v 0 ; y3 = f ; y4 = f 0 . В итоге, решаемая система
уравнений и краевые условия принимают вид:
dy1
dη = y2 ,
dy3 = y ,
4
dη
dy2 1
y2 y4
+
2 = −Re∗ ,
dη
1−y
(1−y
3
3)
dy4 + ψ(y3 ) + y42 + Re∗ y4 (1
dη
1−y3
y2
53
− y3 ) = 0;
η = 0 :
η = 1 :
Re∗
α+β ;
Re∗
0; y4 = − α+β
.
y1 = 0; y2 = Re∗ ; y3 = 0; y4 =
y1 = 1; y2 = −Re∗ ; y3 =
Общая схема численного интегрирования, согласно [28], такова:
1. Задается значение Re∗
2. Задается начальное значение пристрелочного параметра v∗ /U0
3. С заданой точностью в точке η0 стыкуются решения задач Коши
«вперед» и «назад»
4. По полученному решению вычисляются значения Re, Rp и ξ
Методы, задействованные при решении:
1. метод пристрелки (для решения краевой задачи)
2. метод Рунге-Кутта 4го порядка (для решения задачи Коши)
3. метод Ньютона (решение систем нелинейных уравнений)
4. метод Симпсона (вычисление определенного интеграла)
Верификация вычислительной процедуры была реализована посредством сравнения численных решений задачи с решениями в квадратурах,
определенных при одинаковом шаге по η. В результате было установлено, что результаты решений, полученные этими двумя методами, хорошо
согласуются.
54
Заключение
1. Для описания течений жидкости при произвольных числах Рейнольдса можно использовать феноменологическую f − модель течения, которая позволяет выполнить расчеты, не задумываясь о том, какой
режим течения реализуется в конкретном рассматриваемом случае.
Ответ на этот вопрос дадут результаты проведенных расчетов.
2. В рамках используемой модели допускается сравнительно простая
формулировка граничных условий.
3. Результаты расчетов на основе f − модели для пристенных течений
вязкой несжимаемой жидкости показывают, что модель обеспечивает
удовлетворительное их согласование с известными экспериментальными данными.
55
Список литературы
[1] Лойцянский Л. Г. Механика жидкости и газа: Учеб. для вузов. - 7- е
изд., испр. - М.: Дрофа, 2003. - 840с.
[2] Шлихтинг Г. Теория пограничного слоя. М.: Наука. 1974. - 712с.
[3] Хинце И. О. Турбулентность, ее механизм и теория. М.: Физматгиз,
1963. - 630с.
[4] Монин А. С., Яглом А. М. Статистическая гидромеханика, ч. II. М.:
Наука, 1967-720с.
[5] Белов И. А., Исаев С. А., Коробков В. А. Задачи и методы расчета
отрывных течений несжимаемой жидкости. — Л.: Судостроение, 1989.256с.
[6] Павловский В. А. Краткий курс механики сплошных сред. СПб.: изд.
СРбГТУ РП. 1993. - 209с.
[7] Кочин H. Е., Кибель И. А., Розе Н. В. Теоретическая гидромеханика,
ч. I, II. М.: Физматгиз, 1963. - 584, 728с..
[8] Турбулентность. Принципы и применения./Сборник статей/Под ред. У.
Фроста, Т. Моулдена./пер. с англ. М.: Мир, 1980.
[9] Бабкин А. В., Селиванов В. В. Основы механики сплошных сред: Учебник для втузов. - 2-е изд., испр. - М.: Изд-во МГТУ им. Н. Э. Баумана,
2004. - 376с.
[10] Новожилов В. В., Павловский В. А. Установившиеся турбулентные течения несжимаемой жидкости: Монография. СПб.: Изд. центр СПбГМТУ, 1998. - 484с.
[11] Белов И. А., Исаев С. А. Моделирование турбулентных течений, 1998.106с.
[12] Reynolds О. On the dynamical theory of turbulent incompressible viscous
fluids and the determination of the criterion // Phil. Trans Royal Soc.1894.-Vol. 186-pp. 123-161.
56
[13] Смольяков А. В., Ткаченко В. М. Измерение турбулентных пульсаций.
Л. : Энергия, 1980.
[14] Монин А. С., Яглом А. М. Статистическая гидромеханика, ч. I. СПб.:
Гидрометеоиздат, 1992 — 693с.
[15] Короткин А. И., Роговой Ю. А. Метод расчета продольных средних
скоростей в пристенных турбулентных течениях несжимаемой жидкости. СПб.: Мор Вест, 2009. — 121с.
[16] Лукащук С. Ю. Введение в многопоточное программирование и гридтехнологии: Учеб. пособие. Уфимск. гос. авиац. техн. ун-т. – Уфа: УГАТУ, 2008. – 132 с.
[17] Павловский В. А., Шестов К. В. Применение f-модели турбулентности для расчета внутренних задач гидродинамики и тепло-массообмена
// Актуальные проблемы морской энергетики: материалы третьей Всероссийской межотраслевой научно-технической конференции. – СПб.:
Изд-во СПбГМТУ, 2014. 181 с.
[18] Макунин А.В. Применение f-модели для расчета пристенных течений
несжимаемой вязкой жидкости // Процессы управления и устойчивость. Том 2(18). № 1. СПб.: Издательский дом Федоровой Г.В., 2015. с.
213-216
[19] Ламб Г. Гидромеханика. М-Л.: ОГИЗ - Гостехиздат. 1947. - 900с.
[20] Kuethe A. M., Schetzer J. D., Fonduations of aerodynamics. John Wiley,
New York 1959.
[21] Jakob M., Heat transfer, т. 1 и 2, John Wiley, New York 1957.
[22] Новожилов В. В. Теория плоского турбулентного пограничного слоя
несжимаемой жидкости. Л.: Судостроение. 1977. - 165с.
[23] 88. Reichardt H. Gesetzm?ssigkeiten der geradlinigen turbulenten CouetteStr?mung Mitt. № 22 des Max-Planck-Inst. f?r Str?mungsforschung und der
AVA Gottingen. 1959.
57
[24] Глушко Г. С. Некоторые особенности турбулентных течений несжимаемой жидкости с поперечным сдвигом//Изв. АН СССР. Сер. Механика
жидкости и газа. 1971. №4 с. 128-136.
[25] Коловандин Б. А. Моделирование теплопереноса при неоднородной
турбулентности. Минск: Наука и техника. 1980. - 184 с.
[26] Брайтон Д., Джонс П. Полностью развитый турбулентный поток в
каналах кольцевого сечения//Теоретические основы инженерных расчётов. 1964. № 4, с.240-242.
[27] Helsten A. Some improvements in Menter’s к —со SST turbulence model.
Tech, rep., American Institute of Aeronautics and Astronautics (1998).
AIAA-98-2554.
[28] Чистов А.Л. Единая ламинарно?турбулентная дифференциальная модель для течения вязкой несжимаемой жидкости. Вестник СПбГУ.
Сер.10. 2008, вып.4, с. 103?106
[29] Макунин А. В. Расчет стабилизированного изотермического течения
жидкости с постоянными физическими свойствами в круглой цилиндрической трубе на основе f-модели турбулентности // Молодой ученый. — 2016. — №7. — С. 104-111.
[30] Новожилов В. В., Павловский В. А. Установившиеся турбулентные течения несжимаемой жидкости: 2-е изд., испр. - СПб.: Изд-во С.-Петерб.
ун-та,2012. - 453 с.
[31] Миллионщиков М. Д. Турбулентные течения в пограничном слое и
трубах. М.: Наука, 1969. - 52с.
[32] Robertson J.M., Johnson H. F. Turbulence Structure in plane Couette
flow//ASME J. Eng. Mech. Div. 1970. Vol. 96. p. 1171-1182.
[33] Robertson J.M., Johnson H. F. Turbulence Structure in plane Couette
flow//ASME J. Eng. Mech. Div. 1970. Vol. 96. p. 1171-1182.
[34] Bech К. H., Tillmark N., Alfredsson P. H., Andersson H. I. An investigation
of turbulent plane Couette flow at low Reynolds numbers //J. Fluid Mech.,
Mar 10, 1995. p. 286, 291-325.
58
[35] Себеси Т., Брэдшоу П. Конвективный теплообмен. М.: Мир. 1987. 592с.
[36] Павловский В. А., Пастухова Е. В., Шестов К. В. Турбулентное установившееся неизотермическое течение жидкости с постоянными физическими свойствами в прямой круговой трубе//Сб. Кораблестроительное
образование и наука. 2003. Т.П. СПб. Изд-во СПб ГМТУ. с. 37-43.
[37] FLUENT 6.3 User’s Guide Inc., 2005.
[38] Wilcox D. С. Turbulence Modeling for CFD. La Canada, California.: DWC
Industries Inc., 1998. - 540p.
[39] Hanjalic K. One-point closure models for buoyancy-driven turbulent
flows//Annu. Rev. Fluid Mech. 2002.Vol.34. p.321-347.
[40] Spalart P. R., Almaras S. R., A one-equation turbulence model for
aerodynamic flows//La Rech. Aerospatiale. 1994.Vol.l, p.5-21.
[41] Launder B.E., Spalding D.B. The Numerical Computation of Turbulent
Flows // Сотр. Meth. Appl. Mech. Eng. - 1974. - Vol. 3. - P. 269-289.
[42] Самарский А. А. Введение в численные методы. М.: Наука, 1987. —
286с.
[43] Калиткин H. Н. Численные методы. М.: Наука, 1978. - 512с.
[44] Чистов А. Л. О построении моделей турбулентных течений в пакетах
прикладных программ//Процессы управления и устойчивость. Труды
XXXVIII международной научной конференции аспирантов и студентов (CPS’07). СПб.: Изд-во СПбГУ, 2007. с.221-226.
[45] Павловский В. А., Чистов А. Л. Связь между истинной и изотропной
диссипациями через осредненные параметры течения//Проблемы экономии топливно-энергетических ресурсов на промпредприятиях и ТЭС:
Межвуз. сб. науч. трудов/ГОУВПО СПбГТУ РП. 2007. с.3-4.
[46] Газизов Р.К., Юлдашев А.В., Ямилева А.М. Опыт использования гибридных вычислительных систем при моделировании линейной сварки
59
трением // Rational Enterprise Management / Рациональное Управление
Предприятием (REM), 2012, № 04/2012, Автоматизация аэрокосмической отрасли – С. 50-52.
[47] Чистов А. Л. Единая ламинарно-турбулентная дифференциальная модель для течений вязкой несжимаемой жидкости. Вестник СПбГУ. Сер.
10. 2008, вып.4, с. 103-106.
[48] Пантакар С. Численные методы решения задач теплообмена и динамики жидкости. М.: Энергоатомиздат, 1984. - 152 с.
60
Отзывы:
Авторизуйтесь, чтобы оставить отзыв