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

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

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

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

Эти способы перечислены в подписи к рнс. З.ЗО. Первые четыре способа были опробованы как со схемой <чехарда», так и со схемой с разностями против потока, последние же три способа были опробованы только со схемой с разностями против потока. В качестве тестовой задачи была выбрана задача об обтекании обратного уступа при Ке = 1О, когда на входной границе задавался профиль Польгаузена,. соответствующий течению в пограничном слое с параметрами б//г = 1 и Л = О, а на твердой стенке задавалось условие прилипання.

(При ббльшнх Ке результаты мало зависели от выбранного способа расчета.) Способ 1 уже был описан. Здесь рассматриваются разрывные значения ~„полученные из граничного условия (ЗАЗ5) на твердой стенке. Такую постановку граничных 'условий предло- 260 ЯЗ.

Граничные условия жил Ричардсон [1910[, а для расчета вихря ее реализовали Том и Апельт [1961), Роуч и Мюллер [!970[, Кейкер и Уайтло [!970). Способ 2 был разработан Кавагути [1965[ в предположении, что функция тока»р симметрична относительно точки С. Здесь в двух точках на стенке определяются фиктивные значения функции тока»[»', а именно»р'„» /,=ф»..,., и»ры и вихрь в точке С находится из уравнения Пуассона (3.365) с учетом этих фиктивных значений тр". Фактически (»с„/с+1) Р ° ° ° °,и ° 'зч ° (»'с+1 „/с) ~7ТЯБ ь у~ а Рис. З.ЗО. Семь различных способов определения значения вихря в угловой точке.

1) разрывные значения: ьа — — 2»Р»с,/с+,/(ЛУ)', ЬЬ вЂ”вЂ” 2»Р»с».», /,/(йх)'1 2) симметрия»р относительно угловой точки С: ьс = 2»р»с, /с+»/(Ьу)'+ 2»р»се»,/с/(ох)~' 3) среднее от значений на стенке: йс»р»с,/с+»/(Лу)'+»р»с+»,/с/(Ьх')1 4) стенка, наклоненная под углол» 46'» ь 2»рр/(Лр)'1 6) отрыв потока в точке С: йс О; 6) отрыв потока в точке ГП 4ь= в все=2»р»с,/се»/(су) ' 7) значение на с~вике перед точкой С: ьс = 2»р»с /с~-»/(ду)5 в данном способе величина »„ в угловой точке определяется так, как будто эта точка является внутренней точкой, причем здесь згз»р = ьс вычисляется с учетом условий прилипания и = = и = О. В том случае, когда Ах = Лу = Л, данный способ эквивалентен другому способу, предложенному Ричардсоном [1910) для закругленного угла, в котором трс и Лр находятся с помощью интерполяции в точке Р на прямой, проходящей через точки (/с,/с+ 1) и (»с+1,/с).

Способы 3 и 4 также испытывались для скругленного каким- либо образом угла. В способе 3 рассматривается одно значение для вихря в угловой точке, равное среднему арифметическому значений вихря на стенках, т. е. в точках А и В, Это дает величину Гс, равную половине значения ~„ полученного в способе 2. То же самое значение может быть получено из уравнения Пуассона (3.365), записанного для угловой точки, при условии, что тр = 0 в точках (»с — 1,/с) и (»с,/с — 1).

Такая интерпретация предлагалась рядом авторов (см., например, Гри»»спэн [1969а]), причем считалось, что условие для величины ~, имеет второй порядок точности. Такая интерпретация неправомочна; поскольку конечно-разностпая форма лапласиана, на которой основано уравнение (3.365), несправедлива в угловой точке, нет оснований предполагать, что в точке С вторая производная З.Ви2.

угловые точки даф/дх' непрерывна по х илн дтф/дуа непрерывна по у, поэтому вывод уравнения (3.365), проведенный с помощью разложения в ряды Тейлора, неприемлем для угловой точки С. Отметим, что пятиточечная схема для лапласиана )/ятр справедлива в точке С только при условии непрерывности производных даф/дха и датр/дуа при переходе через точку С, или, что эквивалентно, при условии непрерывности дп/дх и ди/ду. Но односторонний предел при стремлении к С вдоль границы В 2, когда х : х„ в силу условия прилнпания на стенке В 2 приводит к равенству ди/дх = 0; аналогично ди/ду = 0 вдоль В 5 '). Поскольку ", = ди/ду — ди/дх, ясно, что из непрерывности производных датр'дха и д'тр/дуа в точке С следует, что !., = О.

Отсюда можно прийти к заключению, что если и можно вычислить т)атр в точке С, то этого делать не надо, а можно сразу однозначно записать с, = О. В способе 4 условно считается, что стенка в узле (тс,/с) имеет наклон 45' и формула (3.435) применяется в точке Р (рис. 3.30, б), где значение трл находится с помощью интерполяции по значениям тр',)сч.! и ф„е!,;с.„, '). При р =Лх/Лу= 1 эта формула сводится к условию ьс = 2тргсч.!, тс к!/Лра, или ьс —— = трчс.„!, )г.„!/ЛЯ, где Л = Лх = Лу.

Способы 5 и 6 основаны на введении искусственного отрыва потока в угловой точке, где его можно ожидать. Легко видеть, что для непрерывного потока вихрь равен пулю в точке отрыва (или в точке вторичного присоединения потока). В способе 5 рассматривается одно значение ь„ которое равно нулю, в способе 6 рассматриваются два значения вихря в угловой точке, причем в точке В принимается ~ь = О. В общем случае такой отрыв лучше не вызывать искусственно, а давать возможность ему проявляться в результате решения, поэтому ни один из этих двух способов не считается правильным.

(Как ни странно, этн способы не эффективны даже при искусственно вызванном отрыве, см. Роуч и Мюллер [1970].) Способ 7 предложили Хын и Макано [1966). Идея заключается в следующем; поскольку разделяющая линия тока (предполагается, что она отрывается от угловой точки) почти параллельна степке, расположенной выше по потоку.

от угловой точки, то следует брать одно значение вихря с,„ равное велитпп!е с,. Хын и Макапо применяли для вихря на стенке условие второго порядка точности (3.439) в отличие от рассматриваеь!ого нами условия первого порядка (3.435), Способ 7 приводит по существу к таким же результатам, как и способы 2 и 4. ') Простой апалпя профиля скорости и окрестности точки С прн ди/ду=о покатынает, по ато условие абсурдно. ') Прн р ) 1; селя р ( 1, то фе нахощпся прн попон!я пнтерполянпп по точкап ис + ! !с + П и ис + 1.

ч') дд Граниеные уелоеил 262 Если применяются сетки второго типа (равд. 3.3.3) и пяти- точечная схема для лапласиана ЧаЬ в уравнении переноса вихря, то точка С, расположенная в вершине выпуклого угла, не требует специального рассмотрения. Если же для Чаь берется девятиточечная схема, то для Гс,ус 1 зс»11се1 определения Г» на сетке второго ° и ° типа могли бы использоваться л! способы, аналогичные способам !г л 2, 3 и 4. — ° — е--~- ° ° О'Лири и Мюллер [1969) й,12 ~ также применили постановку разрывных значений вихря Ь, в ° ° угловой точке, подобной точке С а иа рис. 3.22, а также в угловой точке, расположенной в вершине тупого угла в 135' на квадратной сетке, как показано на рис.

3.31, а. В последнем случае значение вихря в точке А (3.499) 2~2 можно брать для составления разности во внутренней точке ((с, )с + 1), а значение вихря в точке В 2»Р~»+ь м»л Ф~» ьь 1»ь~ ьь= але = ле ° ° (3.500) ° ° — для составления разности во внутренней точке ((с + 1, 1с). рие З.ЗП Вихрь на выпуклом угле Такие постановки разрывных в 135' и на задней кромке панской значений ~ в угловой точке рекомендуются в задачах о течении около передней пли задней кромок бесконечно тонкой плоской пластинки; здесь, как показано на рис. 3.31, б, рассматриваются три значения вихря ь», ь» и ь». (Если не предполагается симметрия течения, то, очевидно, на разных сторонах пластинки необходимо задавать различные значения ь.) В связи с втой задачей отметим, что Иосидзава [1970] (а также другие авторы; см. равд.

6.4) численно решал задачу о течении в окрестности передней кромки плоской пластины, используя уравнения Нннье — Стокса в параболиче дЗ.Г2. углосьм тоти ских координатах (Ч,т1), оптимальных для расчета быстрорастущего пограничного слоя. Далее он перешел к новой зависимой переменной О = — Ь Я' + П'), что устранило сннгулярность на передней кромке. Несмотря на то что сравнительное качество каждого из рассмотренных выше семи способов дискутабельно, очевидно, что подходы, основанные иа определении ь, при помощи произвольной экстраполяцпи по значениям во внутренних точках, некорректны п могут привести к неустойчивости, З.ЗЛЯ. б. Сходимость и точность в вершине выпуклого угла Специальные вопросы, связанные с численным решением уравнений эллиптического типа в окрестности выпуклого угла, обсуждали Вудс [1953], Вазов [1957], Лаасонен [1958а, 1958б] и другие авторы.

Для того чтобы продемонстрировать сходимость решения разностного уравнения в случае конечного числа разрывов функции ь, Лаасонен [1958б] рассматривал решение уравнения Пуассона, записанное в интегральной форме при помощи функции Грина (см. Вейнбергер [1985]). Доказанная им теорема требовала, чтобы решение дифференциального уравнения для ь было кусочно непрерывно и чтобы разрывы находились между узловыми точками сетки. Второе из этих условий не удовлетворяется в наших разрывных постановках для ~.

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

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

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

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