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

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

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

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

следующими из дифференциального уравнения неразрывности и имеющимн вид ди/ду) = 0 и о =О. Действительно, с точностью первого порядка из условия дп/ду~ = 0 следует, что и„+~ = О. Однако из рис. 5.3 интуитивно ясно, что подобное условие не может быть достигнуто при расчете; в результате решение расходится '). Ситуация еще более усложняется при расчете течения вдоль уступа (граница В 5 на рнс. 3.22). Основной причиной неудач схемы с односторонними разностями вблизи точки отрыва потока является отсутствие переноса массы вдоль стенки, поскольку всегда 6(ри)/бхгом = О. В то же время применение метода контрольного объема показывает,что ячейка с центром (й ю) на рис.

5.3 может терять или приобретать массу только через границу гп + '/з. Эту ситуацию можно ') В подобной сптуаиии сходимость может быть достигнута за счет вве. денни члена с искусственной диффузной массы, как н схеме Русзпова (равд. 5.4 3). тогда конаективный поток массы, уносимый нз узловая точки на стенке, при отрыве может быть сбалансирован притоком ее в эту узловую точку за счет яскусственной диффузии Но точность такого способа и даже его аппроксимирующие свойства соанительны, 400 б.7.

Грани сные условия для течений сжил~асмой жидкости изменить, переопределив поток массы, параллельный стенке, для прилежащих к стенке ячеек. В обозначениях, указанных на О О О ~'-1, се О й' Рис. 5.4. Возможная интерпретация скоростей у стенки с условием прилн- пания в расчетной сетке первого типа.

о — стенка, б — внутренняя точка. рнс. 5.4, а, получим (с применением центральных разностей) следующее значение потока массы на границе между ячейками: (рн)с-Нт а, рс-пт,мне-~йьее ~2(РК + рс-Н и> с-иа,ме =',с ( + р )и. (5.131) где и..е и ='7,(и +ис)='/я('/,и,, „+, + '/,ис е,), (5,132) Аналогичная формула имеет место и для нсый, . Этот способ является аппроксимирующим, устойчивым и консервативным да.

же в случае отрывных течений (Роуч и Мюллер [!968, 1970[, Скоглунд и Гей [1969[), Однако вычисления во внутренних н в прнстеночных ячейках проводятся здесь не единообразно. Такой расчет потоков в пристеночных ячейках был бы согласован 401 д7 д Стенка е прилипанием с расчетом потоков во внутренних ячейках (рис. 5 4, б) при применении следующей формулы; и, = ~/,(~/4(и,, +, + 2и,, + и...) 1 + ' 4(и, 7~, + 2и,, + ис ..)). (5.133) Но подобный способ расчета для внутренних точек приводит к тому, что поле скоростей получается неточным и сглаженным. Устойчивые решения дают другие аппроксимации. В частности, широкое распространение получила аппроксимация условия дР(ду~ = О равенством (5. 134) Р =Р„эн Тогда р„находится через Т из уравнения состояния (4.51).

Может показаться, что эта аппроксимация основана на приближении теории пограничного слоя, где поперек пограничного слоя принимается дР7ду ж О (см. Шлихтинг (1968]). В действительности же это гораздо менее жесткое условие, так как постоянство Р предполагается не поперек всего пограничного слоя, а только поперек прилегающего к стенке подслоя толщиной Ьу. Этот способ дает возможность получать устойчивое численное решение как для течения в безотрывном пограничном слое (Курцрок и Мейтс (1956)), так и для течения с отрывом потока, вызванным взаимодействием ударной волны с пограничным слоем (Мак-Кормак [1971) ). Впоследствии Мак-Кормак повторил свои расчеты при более точных граничных условиях и фактически не обнаружил различия в результатах (личное сообщение).

Хотя этот способ привлекает своей простотой и в некоторых случаях дает достаточную точность, в общем случае его рекомендовать нельзя. Во-первых, он не консервативен. Во-вторых— и это более существенно — он дает аппроксимацию решения исходных дифференциальных уравнений в частных производных в некотором смысле, но не аппроксимацию в строго математическом смысле, т. е.

решение конечно-разиостпых уравнений при Лх †« О не стремится к решению исходных дифференциальных уравнений с точными граничными условиями. Более того, при взаимодействии сильной ударной волны с пограничным слоем, при малых числах Рейнольдса, при возникновении отрыва и при наличии сильно искривленных стенок может теряться всякое соответствие между решением конечно-разностных уравнений и исходных уравнений в частных производных.

Нет необходимости применять этот способ, поскольку имеются другие, хотя и несколько более сложные, способы, обеспечивающие аппроксимацию задачи, 402 д К Граничные условия для течений сжимаемой жидкости Очевидно, что такие искусственные приемы, как произвольная экстраполяция значений р или Р на стенку, не обоснованные физически даже в качестве приближенного приема, не обе. спечивают консервативность и аппроксимацию.

Кроме того, в общем случае они неустойчивы при расчете отрывных течений. Мы рекомендуем проводить расчет р около стенки в расчетной сетке второго типа, а затем значения р, полученные на стенке, использовать для нахождения градиента давления около стенки в расчетной сетке первого типа. Полностью эта методика будет изложена в равд. 5.7.2,в.

б.7.2, б, Стенка с прилипанием в расчетной сетке второго типа (5. 135) (5. 136а) (5.! 36б) Рн = Рн»1 и„=и „,, "н пн»1 Линейная интерполяция по значениям (5.136) ист = О. Если задана температура стенки Топ то в принятой линейной интерполяцией имеем Т„='ЫТ е, +Т ), дает и„=О и соответствии с (5.!37) (5.138) откуда Те =2Т, — Таян Условия прилцпания на стенке для скорости могут быть поставлены и в расчетной сетке второго типа, но с ухудшением точности. Практически удобно частично использовать способ отражения. Непосредственное применение способа отражения приводит здесь к серьезным ошибкам, однако из способа отражения мы будем брать лишь методику, которая позволит удобно ставить некоторые нз точных граничных условий; другие же граничные условия будут ставиться явно, не соответствуя способу отражения.

Таким образом, способ отражения здесь будет играть роль лишь некоторого приема программирования, и мы согласны с Моретти [1968а, 1968б] в том, что этот способ не заслуживает названия «принцип»с Определим сначала по способу отражения функции в фиктивных точках внутри границы и исследуем влияниетакогоопределения на члены уравнений сохранения. Это исследование покажет, для каких членов способ отражения дает неправильный результат и, следовательно, какие члены должны быть рассмотрены отдельно. Согласно способу отражения для стенки с прилипанием, показанной на рве.

5.2, б при ш + '/ь имеем 4ОЗ 5.7.2 Стенка с прихипиниел Если же задано число Нуссельта )к(п = (Т „, — Т )/Лу, (5,139) то получим Т„=Т.„— Лу Н . (5. 140) а (Рио).„= а — ~Р.„(~, ) ~ 4 +0(Лу')1, (5.141) что опять является аппроксимацией точных граничных условий, так как член в квадратных скобках стремится к пулю при Лу- -~0. Ошибка аппроксимации для составляющей количества движения в направлении у может даже быть меньше; действительно, в стационарном случае уравнение неразрывности — + — + =0 др д (ри) д (ри) де дх др (5,142а) для стенки с прилипанием сводится к условию (5.1426) Далее, (5.143) или (поскольку о„= 0) (5.144) Таким образом, ошибка аппроксимации первого порядка в выражении (5.123) в стационарном случае будет равна нулю.

Можно показать, что ошибка аппроксимации в стационарном случае имеет вид — (ро )и+, = а '(р "~ в4 ( ~ е ~ ~ + 0 (Лу~) ~, (5,145) ст где в фигурных скобках стоит ошибка в величине потока по нормали к стенке составляющей количества движения в направлении у. Хотя ошибки в определении потоков количества движения, даваемые выражениями (5.14!) и (5.!23) илп (5.!45), являются Комбинация соотношений (5.!35) и (5.!366) дает нулевой поток массы в направлении у при ге + '/, для случая стенки со скольжением.

Но в случае стенки с прилипанием потоки в направлении у составляющих количества движения в направлении х и в направлении у будут получаться с ошибкой. Следуя выводу выражения (5.123), получаем 404 б 7. Граничные условия длл течении сжимаемой жидкости приемлемыми, их можно избежать. Можно рекомендовать изменить выражения в точке то+ 1, явно положив соответствующие потоки равными нулю. Как показывает следующее ниже упражнение, аналогичное положение справедливо и для уравнения энергии. Упражнение, Рассматривая уравнение энергии, показать, что способ отражения дает правильное нулевое значение потока величины Е, + Р в на.

правлении у на стенке только в частном случае адпабатвческой стенки (Л(п = О). (Для более общего аида температурных условий иа стенке этот поток следует явно полагать нулевым.) Применение стандартных конечных разностей с отраженными значениями функций не дает также правильной величины градиента давления дР/ду] +ь В случае расчетной сетки второго типа рекомендуется обратиться к односторонним конечным разностям для дР/ду, что, к сожалению, дает аппроксимацию градиента давления лишь первого порядка: — ' + 0(ду). (5 146) бу )м-~-г йу Была показана устойчивость этого способа даже при отрывных течениях (Аллен (1968], Аллен и Чен (1970], Роуч и Мюллер ]!968, 1970]).

Аллен ]!968], Скоглунд и Гей ]!968] предложили рассчитывать градиент давления по уравнению составляющей количества движения в направлении у, записывая его через односторонние конечные разности; однако не представляется, что этот способ быстро ведет к цели. Аллен 11968] применил улучшенный способ расчета конвек. тивных потоков около степки с прилипанием в расчетной сетке второго типа. Он обнаружил, что в задаче обтекания обратного уступа, показанного на рис. 3.22, на верхней части уступа (граница В 5) иногда могут возникать отрицательные значения плотности. Эта теадеицня усиливается при уменьшении числа Рейнольдса и при измельчении сетки ').

Аллен объясняет это неточностью расчета по линейной интерполяции потока массы в примыкающей к стенке ячейке (ячейка го + з/з на рис. 5.2,б). Величина (ро) )ст =0; кроме того, в стационарном случае пз уравнения неразрывности следует, что д(рп)/др]ст = О. Отсюда видно, что вблизи степки рп изменяется по квадратичному, а не по линейному закону. Поэтому Аллен ввел квадратичную интерполяцию для рп, выбирая не значения в трех узлах сетки, а значение в двух узлах сетки гг + 1, тп + 2 и известное значение ро на стенке ы + '/м т. е, положил (Рп)оезз = /з] (Ро)м„, + 3(ро),„+, — (рп)„еиз1, (5.147) ') Роуч и Л(юппер 1!9701 не обнаружили отрицательным значений плотности, возможно, из.за грубости расчетной сетки.

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

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

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

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