Главная » Просмотр файлов » Роуч П. Вычислительная гидродинамика

Роуч П. Вычислительная гидродинамика (1185913), страница 90

Файл №1185913 Роуч П. Вычислительная гидродинамика (Роуч П. Вычислительная гидродинамика.djvu) 90 страницаРоуч П. Вычислительная гидродинамика (1185913) страница 902020-08-25СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 90)

Кеицер [1970а] определял Р .. по выражению (5.118) (в котором легко узнать равенство центробежной и центростремительной сил), представляя его односторонними конечными разностями. С учетом условия постоянства энтропии 5 (в невязком газе) вдоль линии тока, совпадающей со стенкой, оно дает возможность найти плотность р., = о(Р, х) (может быть, этот способ можно применить и к случаю падения ударной волны на стенку, когда энтропия вдоль стенки не постоянна). Кенцер дополнительно рассматривал характеристическое соотношение, как это ранее делали Моретти [1968а] и Ьастианон [1969] (см. также Кокли и Портер [1969]).

Хотя такой способ правилен, он не приспособлен к расчету течений вязкого газа, не является консервативным ') и довольно сложен в применении. Другой способ заключается в представлении односторонними разностями членов с потоком в направлении у во всех уравнениях. Например, для уравнения непрерывности будем иметь "+0(э )= Р +' +0(Л ). (5.119) ау ) ау " ду у . Если о 1. О, то, согласно формуле (5.1!9), потоковые величины в направлении у на стенке представляются копечнымп разностямн против потока. Однако при о ы ( 0 будем иметь конечную разность по потоку, что физически абсурдно и приводит к неустойчивости (см.

равд, 3.1.8). Хотя прн помощи этого способа и были получены некоторые сходящиеся решения, в общем случае его рекомендовать нельзя. Лапидус [1967] разработал конечно-разиостную формулировку граничных условий в узловых точках на стенке, дающую конечио-разностную задачу, аппроксимирующую исходную. Однако его способ не является строго консервативяым (хотя и основывается на концепции консервативности) из-за вынужденной перестройки и перекрытия границ пристеночных ячеек, что не согласуется (в обычном смысле слова) с подходом, применяемым во внутренних узловых точках. Кроме того, этот способ также довольно сложен.

') Ыи Моретти, нп Кенпер не применяли консервативных уравнений и, таким образом, совсем пс рассматоивалн пх в своих работах. 9.7 1. Стенка с условием скольекенил Вогачевский и Костофф [1971[, а также Эмери н Ашерст [197! [ в случае течения невязкого газа сформулировали дискретные граничные условия на стенке при наличии вдува, используя расчетную сетку первого типа. В общем случае для стенки со скольжением мы рекомендуем расчетную сетку второго типа. 6.7.1.

б, Стенка со скольжением в расчетной сетке второго типа ! раничные условия на стенке со скольжением проще поставить на расчетной сетке второго типа, которая изображена на рис. 5.2,6. Здесь со стенкой совмещаются не узловые точки сетки (центры ячеек), а стороны ячеек. Поэтому ближайшая к стенке узловая точка ш + ! удалена от стенки на расстояние Лу/2, Внутри стенки можно ввести фиктивную узловую точку ю. Джентри, Мартин и Дали [1966[ рассмотрели в этой расчетной сетке второго типа граничные условия на стенке со скольжением.

Единственным условием, которое здесь необходимо соблюсти, является условие равенства нулю потоков через стенку всех переносимых газодинамических функций. Это условие можно поставить, либо записав специальные уравнения, выражающие равенство этих потоков нулю для пристеночных ячеек с центром в точке ш+1, либо применив к некоторым из этих членов способ отражения, Введя фиктивную ячейку с центром в точке ар внутри стенки, положим (5.!20а) (5.1206) ом омеь Т,„=+/ вь где1=р, и, Е, или Т. Применение способа отражения в расчетных сетках первого и второго типов дает совершенно различные результаты ').

При использовании в точке ш + 1 аппроксимации второго порядка, принятой для стандартных внутренних точек, потоки всех величин 7 на стенке (го+ '/,) обращаются в нуль. Это легко показать при помощи метода контрольного обьема, примененного для уравнений с центральными разностями (см. равд. 3.1.1). Значения потоков величин на стенке (о!) лига определяются следующим образом: (сй „„='/ [(о/) +(с/)„,г)='/ ( — (и[) „+(о/),Д=0 (5.! 2!) Здесь значения 1 на стенке при ш+ '/, отвечают линейной интерполяции; однако необходимо подчеркнуть, что линейная интерполяция используется не для простоты, а потому, что '1 Моретти (1969а, 1969б) рассматривал способ отражения только для соток первого типа.

а 7.2. Стенка с пралааанаел Покажем, что точки внутри стенки являются на самом деле фиктивными. Соотношения (5. 124) (5.125) обеспечивают нулевые потоки величин 1 по нормали черезстенку при ш+ '/,. Однако если за 1 принимать одновременно и р, и рп, то возникает алгебраическое противоречие. Поэтому для получения нулевых потоков через стенку необходимо брать соотношение (5.125) для всех величин 1= р, ри, ро, Е, + Р, игнорируя алгебраические связи между ними. Следует подчеркнуть, что способ отражения должен быть модифицирован для првменения к любой искривленной стенке путем использования лиоо равенства (5.118), либо односторон.

них конечных разностей для бР/бу~ Таким образом, удобный способ отражения после небольшой модификации может давать граничные условия на стенке со скольжением на расчетной сетке второго типа. Эвристические граничные условия на молекулярном уровне для течения около стенки со скольжением на расчетной сетке второго типа разработал Батлер [!967) с использованием уравнений сплошной среды на некотором расстоянии от стенки. б.7.2. Стенка о прилипанием В постановке граничных условий на стенках с прилипанием в течениях вязкого газа к настоящему времени нет полной ясности. уравнение неразрывности здесь не изменяется по сравнению со случаем течения невязкого газа, и плотность лучше рассчитывается в расчетной сетке второго типа, однако другие переменные точнее аппроксимируются в расчетной сетке первого типа. Поэтому невольно напрашивается применение гибридной сетки, и оно, действительно, оказывается успешным.

Однако несколько более простым выходом является расчет значений р около стенки так, как если бы использовалась расчетная сетка второго типа, но найденные значения р приписываются узлам расчетной сетки первого типа. Хотя в ближайшем будущем могут появиться более эффективные способы, однако представляется, что в настоящее время последний способ является наилучшим. Для того чтобы показать преимущества предлагаемого способа, необходимо рассмотреть более илп менее подробно и другие возможности, 396 Б Г )'дпиичиые условия г)ля течений гжилаемой жидкости б.у.й.а. Стенна с прилипанием в расчетной сетке первого типа Для течения вязкого газа в расчетной сетке первого типа три из четырех необходимых граничных условий поставить легко. Составляющие скорости просто равны нулю, и =О, п„=О.

(5.126) В случае заданной температуры стенки Т„ Тв Тст. (5.127) В случае же заданного на стенке числа Нуссельта (безразмерного потока тепла) мы имеем ') )4 в+! (5. 128) Ьу или (5. 129) Т =Т +, — [ч)п Ьу. Если при рассмотрении во внутренних точках ыг+ 1 членов второго порядка, описывающих теплопроводность, надо обеспечить адиабатичность в вычислительном смысле (см. равд. 3.6.4), то следует брать формулу (5.128). Трудность заключается в формулировке граничного условия для плотности. Здесь, как и в случае иевязкого газа, уравнение неразрывности можно аппроксимировать при помощи односторонних конечных разностей.

Если величина и „1 достаточно мала и если в схеме имеется достаточное искусственное затухание, то можно получить устойчивое и сходящееся решение. Так, Скоглунд и Коул [1966) решили задачу о взаимодействии ударной волны с пограничным слоем, используя схему Русанова (равд 5.4.3) ') и односторонние конечные разности для до/дг'[„. Однако когда интенсивность скачка была достаточна для того, чтобы вызвать отрыв пограничного слоя, схема переставила работать. Этот факт подтверждается также работами Роуча и Мюллера [1970[ и Аллена и Чена [1970], посвященными расчету обтекания обратного уступа. Причину отказа схемы легко обьяснить. ') Заметим, что число Нуссельта ие содержит никакой ошибки аппроксимации и задается как парамстр задачи наподобие чисел М или не (см. равд, 3,6.4), Если иа стенке задан ненулевой поток тепла, то иеобхо. димо обратить вккмание иа правильное определение числа Нуссельта на стенке; если размерный коэффициент теплопроводиости й ие является постоян.

иым, то при расчете числа Нуссельта иа стенке надо учитывать, что пря нормировании в качестве характерной величины берется значение й а иевозму. щенном потоке. ') При применении односторонних конечных разностей аии сохраняли члены с явной искусственной диффузией. С другой стороны, при расчете течения с условием скольжения ка стенке Кесслер [19661 использовал способ отражения и полоясил юот члсп равным пулю. 399 5.7.л. Стенка с прилипанием Рассмотрим картину линий тока отрывного течения, представленную на рис. 5.3.

Прн любом сколь угодно малом Лу ясно, что ос +1 ) О. Из условий прилипания для уравнения неразрывности в точке ((,ю) имеем (5. 130а) (5.1306) Если стационарное состояние когда. нибудь будет достигнуто, то бр/61= 0, а из уравнения (5.1306) следует, что и- +1 = О. Это согласуется с граничными условиями прилипания на стенке, (1-1,ю) (й ю) (г'+1, ю) Рис. 5.3. Область отрыва в расчетной сетке первого типа. 1 — линия тока, ограничивающая область отрыва.

Характеристики

Тип файла
DJVU-файл
Размер
12,43 Mb
Тип материала
Высшее учебное заведение

Список файлов книги

Свежие статьи
Популярно сейчас
Как Вы думаете, сколько людей до Вас делали точно такое же задание? 99% студентов выполняют точно такие же задания, как и их предшественники год назад. Найдите нужный учебный материал на СтудИзбе!
Ответы на популярные вопросы
Да! Наши авторы собирают и выкладывают те работы, которые сдаются в Вашем учебном заведении ежегодно и уже проверены преподавателями.
Да! У нас любой человек может выложить любую учебную работу и зарабатывать на её продажах! Но каждый учебный материал публикуется только после тщательной проверки администрацией.
Вернём деньги! А если быть более точными, то автору даётся немного времени на исправление, а если не исправит или выйдет время, то вернём деньги в полном объёме!
Да! На равне с готовыми студенческими работами у нас продаются услуги. Цены на услуги видны сразу, то есть Вам нужно только указать параметры и сразу можно оплачивать.
Отзывы студентов
Ставлю 10/10
Все нравится, очень удобный сайт, помогает в учебе. Кроме этого, можно заработать самому, выставляя готовые учебные материалы на продажу здесь. Рейтинги и отзывы на преподавателей очень помогают сориентироваться в начале нового семестра. Спасибо за такую функцию. Ставлю максимальную оценку.
Лучшая платформа для успешной сдачи сессии
Познакомился со СтудИзбой благодаря своему другу, очень нравится интерфейс, количество доступных файлов, цена, в общем, все прекрасно. Даже сам продаю какие-то свои работы.
Студизба ван лав ❤
Очень офигенный сайт для студентов. Много полезных учебных материалов. Пользуюсь студизбой с октября 2021 года. Серьёзных нареканий нет. Хотелось бы, что бы ввели подписочную модель и сделали материалы дешевле 300 рублей в рамках подписки бесплатными.
Отличный сайт
Лично меня всё устраивает - и покупка, и продажа; и цены, и возможность предпросмотра куска файла, и обилие бесплатных файлов (в подборках по авторам, читай, ВУЗам и факультетам). Есть определённые баги, но всё решаемо, да и администраторы реагируют в течение суток.
Маленький отзыв о большом помощнике!
Студизба спасает в те моменты, когда сроки горят, а работ накопилось достаточно. Довольно удобный сайт с простой навигацией и огромным количеством материалов.
Студ. Изба как крупнейший сборник работ для студентов
Тут дофига бывает всего полезного. Печально, что бывают предметы по которым даже одного бесплатного решения нет, но это скорее вопрос к студентам. В остальном всё здорово.
Спасательный островок
Если уже не успеваешь разобраться или застрял на каком-то задание поможет тебе быстро и недорого решить твою проблему.
Всё и так отлично
Всё очень удобно. Особенно круто, что есть система бонусов и можно выводить остатки денег. Очень много качественных бесплатных файлов.
Отзыв о системе "Студизба"
Отличная платформа для распространения работ, востребованных студентами. Хорошо налаженная и качественная работа сайта, огромная база заданий и аудитория.
Отличный помощник
Отличный сайт с кучей полезных файлов, позволяющий найти много методичек / учебников / отзывов о вузах и преподователях.
Отлично помогает студентам в любой момент для решения трудных и незамедлительных задач
Хотелось бы больше конкретной информации о преподавателях. А так в принципе хороший сайт, всегда им пользуюсь и ни разу не было желания прекратить. Хороший сайт для помощи студентам, удобный и приятный интерфейс. Из недостатков можно выделить только отсутствия небольшого количества файлов.
Спасибо за шикарный сайт
Великолепный сайт на котором студент за не большие деньги может найти помощь с дз, проектами курсовыми, лабораторными, а также узнать отзывы на преподавателей и бесплатно скачать пособия.
Популярные преподаватели
Добавляйте материалы
и зарабатывайте!
Продажи идут автоматически
6374
Авторов
на СтудИзбе
309
Средний доход
с одного платного файла
Обучение Подробнее