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

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

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

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

Дифференцируя выражение, определяющее вихрь, получаем дй д'и д'о д'ф д до ду ду' ду дх ду' дх ду ' (3.436) ') Б самом деле, для профиля Блазнуса в пограничном слое на стенке дзф/дуз = дти/ду' = О; численные расчеты показывают, что в этом случае разностные формулы первого порядка прнаодят к попсе точным результатам, чем разностпые формулы второго порядка. З.З Граннннме деловая 218 Из уравнения неразрывности (2.3) имеем до/ду = — ди/дк. Учитывая это и записывая (3.436) на стенке, находим (3.437) Второй член даи/дка~„равен нулю в силу условия прилппания (и = сопз( = 0).

Член д7/ду[„вычисляется по схеме с разностями назад, имеющей первый порядок точности: дС 1 +0(йу)= —., ~ . де 9 ду] ау аув м (3.438) Подставив это выражение для дв1Р/дуя[„в уравнение (ЗАЗ4) н разрешив относительно с„, получим условие Вудса для вихря на стенке: — — ю + 0(йуа). (3.439) При г,-ь со получается условие, соответствующее обтеканию плоской пластинки и отвечающее условию Вудса (3.439) с добавлением нестационарного члена. Даусон и Маркус [1970] установили, что нестационарный член д~ /д/ в условии (3,440) «весьма сильно снижает устойчивость н должен быть опущен». Условие Вудса с успехом использовал Расселл [1963]; см. также работы Лестера [1961[ н Майкла [1966].

Хын и Макано [1966] независимо вывели условие (3.439) и применили его для расчета течений при умеренных числах Ке (Ке = 300); аналогичное условие было получено и автором этой книги. Теоретические расчеты П. Дж. Тейлора [1968] указали на возможную неустойчивость, связанную с этим условием, прн больших сеточных числах Рейнольдса. В ходе численного эксперимента Ранчел с соавторами [1969] снова независимо пришли к условию (3.439)..Используя итсрационные методы для расчета стационарных течений, они (Ранчел с соавторами [1969]) обнаружили, что данная форма граничного условия второго порядка, вероятно, приводит к расходнмости, особенно при больших Ке и при переменных Лу.

Формулу второго порядка для с в осесимметричном случае оценивали Люгт н Рнмон [1970], сравнив ее с точным решением задачи об обтекании сплюснутого эллипсонда вращения; к сожалению, это решение пригодно только прн Рс = О. Онн рассмотрели выражение, включающее нестационарный член д(. /дб Записанное через локальный радиус кривизны г„это выраже. ние имеет следующий вид: з4 +,/лр — '/,ь +, + '/„лр'йе д~~/дг + 0(Луя). (3,440) 1 — Зу/~г,) + 3 ЗЗ/1(Вг) М.2. Стенка а расчетной сетке лераого тола 219 Наиболее часто используемая формула второго порядка точности для граничного условия для вихря была впервые предложена !венсеном [1959], а позднее Пирсоном [1965а]. Эта формула получается аппроксимацией значений зр в окрестности стенки полиномом третьего порядка с условием дф/ду[„= = и = 0 и имеет вид (см. также Брили [1970)) — тфм+ай -ь~ — Ф+а ! О(Ь а) 2 аут (3.441) Если эту формулу применять вместе с условием (3.441), то расчеты, проведенные по неявной схеме метода чередующихся ') Бердслп (1969) установил, что наилучшей формулой, совместимой с глобальным сохранением физических величин в осесимметричном случае, является выражение Г~ = (2/ — 1)фг,/((/ — '/,) Ьг'], тле г = /Ьг — максимальный радиус.

Прн больших / (г -ь со) это условие переходит в условие (3.436), Эту формулу позже применяли Уилкс и Черчилл [1966] для решения задачи о свободной конвекции при малом числе Грасгофа, но когда Сэмюелс и Черчилл [1967) продолжили эти расчеты, онн обнаружили, что решение неустойчиво, и вернулись к условию первого порядка точности (3.435). Вероятно, при больших числах Рейнольдса расчеты йенсена [1959] также были бы неустойчивы. Расчеты Торранса [1968] н расчеты Бао и Догерти [1969), в которых использовалось условие (3.441), тоже ограничивались умеренными числами Рейнольдса; см.

также работы Саусвелла [1946], Париса и Уитекера [!965). В дальнейшем при сравнении решения осеснмметричной задачи о вращающемся цилиндре с точным решением в одномерном случае Бердсли [1969) обнаружил, что формула йенсена (3.441) даже в случае, когда она устойчива, приводит в действительности к менее точным результатам, чем условие первого порядка '). До последнего времени было широко распространено мнение (предрассудок?), что неустойчивость условия (3.441) как-то связана с тем фактом, что в нем используется информация из точки и + 2, а такие схемы считались обреченными на неудачу. В работе Брили [1970], внесшей в этот вопрос ясность, было доказано, что это мнение неверно. Брили отметил, что представление ф в виде многочлена при выводе условия (3.441) не согласуется с определением и = дф/ду в точке в + 1, если пользоваться центральными разностями; требование согласованности здесь приводит к необходимости использовать для определения и следующую специальную формулу (только в точке и+ 1): и,= — Р[ = ф" ф ~' +' + 0(Л/т) (3.442) м+' ЬЛ (,эз 4 ау г/ дЗ.

ГраниЧнЫе услОВия 220 направлений, оказываются устойчивыми даже при больших числах Рейнольдса. (Конечно, действительно разумной связи между устойчивостью и «согласованностью» нет, причем наше объяснение неустойчивости по-прежнему не более чем упоминавшийся выше предрассудок. По мнению Брили, мы, возможно, просто заменяем несогласованный предрассудок согласованным.) Брили даже вывел формулу для граничного условия, основанную на интерполяционном многочлене Лагранжа четвертого порядка: — звз[з + ~984зы+з 27Иыч.з+4зу +з + 0(я з). (3.443) ьы 18 Дуз у Для согласованности здесь требуется, чтобы скорости в двух точках вблизи стенки вычислялись по особым формулам — 17з)зы+ 9з!зы+~ + 94зщ+з — з)э+з з О сз зз +!— с у), !4зуы — 86зуы4~+ !84зыч-з+4зуыч-з з Ося зз пы+з + ( су) (ЗА44а) (3.444б) и, кроме того, чтобы член д'зР/дуя в уравнении Пуассона был представлен в виде (3.445а) 29Чзы — 84зУы+з + 274зыч.з — 2зуыч.з — (3.4456) 6у' 1ы+з !8 Ьуз Эти формулы также устойчивы при больших Ке.

Значение 9 . полученное по формуле (3.443), очень мало отличается от значения, полученного по условию йенсена (3,441), так как полная ошибка аппроксимации при атом не меняется (Брили [19?О[). Но в неявной схеме метода чередующихся направлений итерации для апач ! (см. равд, 3.1,16) сходятся быстрее, можно проводить расчеты с большими величинами шагов по времени, а суммарное машинное время сокращается вдвое. Однако программа становится сложнее, так как формулы (3.445) при решении уравнения Пуассона с помощью неявной схемы метода чередуюшнхся направлений приводят к появлению членов в узлах, отстоящих от границы более чем на один н!аг; таким образом, расщепленные по времени неявные разностпые уравнения вдоль у в неявной схеме метода чередующихся направлений ие будут трехдиагональнымн.

Для того чтобы исключить эту дополнительную неявность в точках из+ 1 и со+ 2, Брили [1970[ применял простую модификацию метода исключения Гаусса. Может оказаться, что условия (3.445) препятствуют сходнмости при решении уравнения Пуассона методом после- 221 довательной верхней релаксации и ограничивают возможность использования многих прямых методов. (Метод расчета распространения вектора ошибки, рассмотренный в разд.

3.2.8, мог использоваться при обходе узловых точек в направлении х.) Рис. 3.23. Определение значения вихря на стенке, не параллельной оси х. а — пРонзпольный Угол наклона стенки Вп, о — Угол наклона стенки Вп = 45'. Для случая наклонной стенки, проходящей по диагонали через ячейки сетки, как показано на рис.

3.23, значение вихря на стенке можно найти при помощи формулы первого порядка точности (3.435) следующим образом. Значения на стенке в точках тр„ шз, со„ не являющихся узлами сетки, находятся так: где р = Ах~Ау — отношение шагов сетки. Аналогично вычисляются ~, и Е . Далее интерполяцией вдоль стенки определяются значения ~ на стенке в узлах шь гоз, тоз. Этот способ устойчив и, по-видимому, точен (Кемпбелл и Мюллер [1968) ). З.8.2. Стенка в расчетной сетке первого типа 2 (та зрие) Ьге Лпз = (Ьу соз О )т = Лхз 'е 1+ йз (3. 446а) (3. 4466) З.З.

Граничные условия 222 В частном случае, когда угол наклона стенки равен 45' (8 = 1), для точки (1ге,/гс), изображенной на рис. 3.23,б, получается — (3.447) Ьгаь гм — ла где Л = Лх = Лу (О'Лири и Мюллер [1969])'). Такой способ оказывается точнее, чем расчет по формуле 9Г -Ьгмм — А~м,гю л' (3. 448) так как значения зрх входящие в выражение (3.448), берутся в узлах, расположенных дальше от стенки, чем узлы, в которых берутся значения тР, входящие в (3.447). При рассмотрении уравнений движения вязкой жидкости может быть полезно также моделирование условия скольжения на стенке. В этом случае предполагается, что толщина пограничного слоя меньше, чем величина Лу. Это условие эвристически можно моделировать с помощью условия скольжения на стенке, как это делается и в равд, 3.3.5.

На стенке задается значение тр, а вихрь получается из условия Неймана (3. 449) '] Лналогнчную пропедуру, нспользующую оператор, построенный на патпточечном диагональном шаблоне (см. равд. 3.2.10), можно найти а книге Тома н Апельта (1961, с.

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

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

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

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