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

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

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

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

По определению йо Е= с,(т, — т,) (3.545) (3.546) Но се то се (тките) ср а йо -2 так йо (3.547) "отце где а — скорость звука, /ге — газовая постоянная и у = С /С,. Далее, С,1;а 1 "о тле йе т(с,— с„) м,' (т — 1)м', ' (3.548) 1 ! (3. 549) е (т — 1) или Е= (т — 1) м', (т,/т, — В (3.550) Следовательно, Е- 0 при Мо-э.О для Т~/Те) 1. Но если разность температур, которую можно записать в виде Т~/Та — 1, настолько мала, что Е будет велико прн любом малом, но стационарное поле температур можно получить при помощи прямых неитеративных методов, например при помощи метода расчета распространения вектора ошибки (равд.

3.2.8), Уравнение, аналогичное уравнению (3.543) без диссипативного члена Ф, можно использовать также для решения простейшей задачи о диффузии двухкомпонентной смеси, описываемой элементарным законом Фика; см., например, Берд с соавторами [1960). В этом случае Т заменяется концентрацией, а число Пекле Ре числом Шмидта. 287 3 да Прсдсгаилснис диссинасианаа щунннип отличном от нуля М, и т ~ 1, то роль диссипативного члена будет велика.

Заметим, что при Е = 0 в уравнении (3.543) можно Т заменить на — Т, а зто означает, что решения для убывающей температуры аналогичны решениям для повышавшейся температуры. Однако прн Е '= О такие решения будут различаться из-за положительно определенной диссипативной функции. 3.6.3. Нонечно-разностное представление днсснпатнвной фуннции Поскольку и = дф/ду и о = — дф/дх, уравнение (3.544) можно переписать в виде Э=2[(~ сз ) +(~ н )~+(д и — н с), (3.551) или (3. 552) Применяя обычные конечно-разностные формулы с централь- ными разностями по пространственным переменным (равд.3.1.1), определим Ф со вторым порядком точности 0(Ахи,бу'): 1 ) Ес+ь с+~ — Е;+ь с-1 — тс ь с4~ + 4с 4 1. асад + ~Ес с+1+ Ес с 1 — 2$с с Ес+~ с+ Ес 1 с — 2$с с ~и йхи (3.553) нйи (3.

554) где Ь' — источниковый член в уравнении Пуассона для давления (3.525). Это точное соотношение было бы удобным при уже вычисленных ~ и о', но получаемая таким образом функция Ф не всегда может оказаться положительно определенной, Ошибки аппроксимации могут привести к появлению отрицательных значений функции «и+ 5. Такие ошибки, как правило, малы, но появление отрицательных значений Ф может смутить и ввести в заблуждение.

Однако при вычислении функции Ф на стенках с условием прилипания зто соотношение может быть полезным Поскольку Ф вычисляется только во внутренних точках и никогда не вычисляется на границах, пет необходимости приме. нять где-либо односторонние разности. Выражение (3.553) сохраняет положительную определенность диссипативной функции, Легко показать, что имеет место следующее интересное соотношение: аб. Рис гег гюемрагурм и конггентраигггг при построении графиков на ЭВМ. Поскольку на стенке при условии прилипания 5 = О, мы имеем г1г ~г (3.555) т. е.

здесь Ф вЂ” опять положительно определенная функция. 3.6.4. Граничные условия для температуры и нонцентрации Граничные условия для температуры аналогичны граничным условиям для вихря, за исключением условий на стенках, где онп проше. Вернемся к обозначениям, принятым на рис. 3,22. Температуры стенок могут быть фиксированы, например равны характерным температурам.

Так, можно положить, что в задаче с подогреваемым основанием уступа Т(В5)= — 1, Т(В2)=0, (3 556) а в задаче с охлаждаемым основанием уступа Т(В5)=0, Т(В2)=!. (3. 557) Т~~~ =Т + ~ ~ Лп+ ~ — р —,~ Лггз+ 0(Лгг'), (3.558) дт ~ 1 от или — ~= Ми — (Т „, — Т )(Лп+ 0(ЛП). дТ (3.559) Чтобы поставить условие адиабатичности на стенке, положим Ми=О, илп Т.= Т„„. (3.560) Если для решения уравнения (3.543) применяются явные схемы, то это условие просто используется после вычисления Т гг во внутренних точках, т. е.

принимается, что Т" "' = Т"~.'г. Если же применяются неявные схемы, то условие (3.560) необходимо ввести в схему для расчета внутренних точек так же, как это Возможны, конечно, любые другие распределения температур вдоль стенок, например вдоль границы В 5 можно задать Т(гс, 1') = 1(у). Граничные условия для температуры на стенке можно задать в виде условия Неймана, т. с. задать градиент дТ~дгь В безразмерных переменных дТ)дп = Хц, где Мц — число Нуссельта, которое соответствует безразмерной интенсивности теплопередачи. Наиболее обычным видом граничного условия является случай адиабатической стенки, когда Кц = О.

Если же число гчц не задано, то его интересно вычислить в процессе решения. На сетке первого типа разложение в ряд Тейлора дает З.б.4. Граннянам уалоатая дла температуры и коняентрарии лзэ сделано в схеме Миякоды (3.529) при решении уравнения для давления, При вычислении градиента на стенке особенно важно, чтобы его конечно-разностное представление было согласовано с конечно-разностной схемой, используемой во внутренних точках. Если вычисленные значения температур во внутренних точках рассматривать просто как известные величины, то в случае адиабатичсской стенки для обеспечения условия дТ~дп = 0 можно брать экстраполяцию высокого порядка точности.

Один из возможных вариантов формулы для определения температуры на адиабатической стенке имеет следующий вид: (3.561) Т„= 4Т, ~ — 3Т яа. Однако такое представление не является численно адиабатическим. Для определенности рассмотрим стенку, параллельную оси к при 1 = 1, Тогда в точке (й ш +! ) градиент потока тепла дается членом К'Т~18 а, который определяется равенством (3,562) По схеме с центральными разностями получаем Ьтт т + т., — Зт. (т — т, )ИУ вЂ” (г — т; дРУ рта — й (3.563) Такая запись дает возможность понять смысл члена д'Т!дуа, являющегося градиентом потока тепла в точке (й 2).

Используя терминологию, принятую в методе контрольного объема (разд, 3.1.2), можно сказать, что ина представляет собой поток тепла, вытекающий из узла 1 (на стенке) и втекающий в узел 2. Таким образом, для того чтобы стенка была численно адиабатической, ды должно быть равно нулю, что приводит к условию Ть, = Тьв Значит, если не выполнено условие Т, = Т„,~ь то в вычислительном алгоритме возникает передача энергии между стенкой и жидкостью. Экстраполяции высокого порядка для условия дТ(ду~ = 0 будут численно соответствовать адиабатической стенке только тогда, когда они согласованы с разиостным аналогом высокого порядка, принятым для д'Т(дуа в точке (й рр + 1).

При интерпретации формального порядка величины ошибок аппроксимации этих выражений имеет место широко распространенное заблуждение. Из уравнения (3.559), казалось бы, следует, что число Нуссельта находится только с первым порядком точности. Но это заблуждение. Если, как указано выше, во З.Б, Раенет телаератпры и каннентракии 990 внутренних точках используются соответствующие конечно-разностные представления, то уравнение (3.559) будет давать алгебраически правильную интенсивность теплопередачп в безразмерных переменных, которая фактически получается в результате решения при заданной температуре стенки. Однако если на стенке ставится условие адпабатичностп (3.560), то стенка будет численно адиабатической, т.

е. будет отсутствовать передача энергии от стенки к жидкости. В этом случае число Нуссельта определяется точно и является задаваемым параметром задачи, подобно числу Рейнольдса. Тогда вопрос о формальной оцшбке аппроксимации будет состоять не в том, насколько точно определяется число Нуссельта, а в том, насколько точно рассчитана температура на стенке при заданном числе Нуссельта. Записывая равенство (3.559) в виде Т „, = Т + Хц Лп + О (Лп'), (3. 564) убеждаемся, что температура на стенке определяется со вторым порядком точности. Для нахождения поля температур с граничными условиями Неймана можно бы.чо бы брать и сетку второго типа, изображенпую на рис.

3.24. Задание числа Нуссельта на стенке дает возможность определить температуру Т;„, в узле, находящемся внутри стенки, нри помощи формулы Ты, — — Тра — Иц Лп. (3.565) Тогда численно условие адиабатичности удовлетворяется постановкой условия Т„ , = Тна Но если рассматривается случай, когда на стенке задается температура Т, то линейная интерполяция дает Ты, = 2Т вЂ” Тта', (3.566) полученный при этом результат будет иметь первый порядок точности, н здесь возможно появление ошибки, связанной с нарушением ограниченности решения (см, обсуждение граничного условия на стенке для вихря в разд.

3.3.2). Распределение темнературы на входной границе В 4 потока (рис. 3.22) должно быть задано. Форма этого распределения, конечно, произвольна, но особенно простым и разумным является распределение, соответствующее автомодельному решению в пограничном слое на входной границе (Шлихтинг [19681).

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

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

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

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