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

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

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

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

Оказывается, свойство транспортивностп имеет такой же физический смысл, как и свойство консервативности. Схемы с разностями против потока, обладающие свойством транспортивности, точнее, чем схемы с центральными разностями для первых производных именно в этом смысле, а не в смысле порядка ошибки аппроксимации. Чтобы подчеркнуть значение свойства транспортивности в противоположность схеме с разностями против потока, рассмотрим схему с разностями по потоку или наветренные разностные схемы (Франкел [1956)). Такая схема неустойчива, но предполагается, что ее можно сделать устойчивой при помощи некоторого конечно-разностпого представления производной по времени.

С точки зрения точности представления только производных эта схема и схема с разностями против потока одинаково приемлемы. Однако в схеме с разностямн по потоку возмущение будет переноситься за счет конвекции только вверх по потоку, а вовсе не в направлении скорости! Это физический абсурд '), и стоит еще раз напомнить то, что было сказано относительно свойства консервативности: точность конечно-разностного представления производных не эквивалентна точности представления дифференциального уравнения.

') Сполдинг предложил называть схему с разностямн против потока «схемой свинарника». Его мысль. заключается з том, что если Ь рассматривать как конпентранию некоторого вещества, то мы должны почунстзозать запах свинарника, когда находимся на его подаетренной, а не на назсгрен. ной стороне (если исключить влияние диффузии).

а1, !т!етодьт решения уравнения яеренота вихря !!О Робертс и Вейс [1966) называют схему Лелевье с разностями против потока «неадекватной», однако в подстрочном примечании признают: «Тем не менее эта односторонняя схема сохраняет знак положительно определенных величин, так же как и лагранжевы схемы; эйлеровы же схемы с центральными разностями не обладают таким свойством». Интересно отметить, что в рассмотренной конечно-разностной схеме каждая узловая точка сетки аналогична конечному элементу в расчетах конечно-элементной модели реактора, иногда используемой инженерами-химиками (Кридер и Фосс 119661). Поэтому, очевидно, оправдано представление этой схемы как схемы «моделирования».

3.1.10. Транспортнвные н нонсврвативные равностные схемы Первая схема с разностями против потока (3.!76) является консервативной, а также транспортивной до тех пор, пока составляющие скорости не меняют свой знак. Покажем для одномерного течения, когда все ит ) О, что эта схема обладает свойством консервативности. Проводя такие же выкладки, как в равд.

3.1.3, получаем вы итьт — и т!и а! ая (3.199) тт т ! т Л (1= У!) (2'=!!+ Ц 11=7!+2Д = иь 11 ! — и1.1 + + иь 11 — и~1 +и~1л — иь1 „,+ . + + +и~1, — иь" 1,, + +» 11,-2 ь 11,-1 + ит» 11 Г 2 22 ("ю = 1, — 11 [1 = 7,~ (3.200) — иь 1,, или 1т ) — „' = иь 1л, — и~ 1, . (3.201) Правая часть представляет собой вычисленную против потока разность двух величин, а пменно потока, втекающего в об- З.С,СО. Траненартивные и кансервативные схемы 111 пасть ст через ее границу 1с за единицу времени, и потока, вытекающего из области сс через ее границу се за единицу времени.

Таким образом, свойство консервативности здесь имеет место. Однако в области, где скорость меняет знак, направление разностей против потока также меняется, н тогда свойство консервативности утрачивается. Покажем зто для одномерного случая при и, ) О для с ( ( и и, ( 0 для с ) й Применяя схему с разностями против потока, получаем Лй иС),. — иЬ(с — — — при с(~с, где и, > О, (3.202) ЛЬ ит~,ьс — иЬсс — — при с > с, где и; < 0; (3.203) тогда с=с, с-с, =и~ 1 с, )с с — иьу!с Л С усу~с — сс~1сн.с Лх х'. Лх Лх+ Лх= с-с+с + — и~ сс с+ + — и~) +иЛ!, + +иьсс 2 — сто!с с + +и~~с с — ю4!с + + ссь ~с е с — иь!с ьа + +и~)„, — иР„в+ -(- + + и~ 11, — и~ ~с. с + + и~ ~с с — нс".1 + и~1с — се~~с вн (3.204) или + = и~ сс, — ис,'1 +, — ис, 1с + и~ 1с с (3.205) с-с~ Первые два члена представляют собой вычисленную против потока разность потоков в единицу времени на двух границах области )г; третий и четвертый члены представляют собой осинки, обусловленные неконсервативностью.

Поскольку ис ) 0 н ис+, ( О, зги члены можно переписать как — иь ~с + иь ~сн.с = — [) ис ~ьс+ ~ ис,с ~ ~с.н) (3 206) 3.1. !негода! решения уравнения переноса вихря !!2 где 7с = — (ис !сс, !ис с!сс ! при и;,)О, 1~-! 0 при ис с(0 о (3.209) (3. 210) с 0 при и;,)О, с с -!- ! ! ись! (~сь! прн и;е! ( О. (3. 21! ) Это всегда имеет место в том случае, когда в качестве ~ь рассматривается плотность (положительно определенная величина), входящая в уравнение неразрывности для сжимаемой жидкости.

х) В плоской задаче такая ошибка появляется только при условии, что составляющая скорости меняет знак при движении в напранлении этой составляющей. Таким образом, если о,, ) О, а о! оы ( О, то ошибка, обусловленная неконсервативностью, появляется, но если ос,, ) О, а осс,, с < О, то такой ошибки нет. Первая ситуапия встречается в отрывных течения., вторая — в случае, когда прп обтекании составляющая скорости о меняет знак при переходе от лобовой частя к кормовой без отрыва потока. Таким образом, разностные решевия Скалы и Гордона (1967) консервативны.

з) Такой частный случай в действительности эквввалентен применению схемы (3.!76), когда скорость не меняет знак. Эта схема при расчетах используется вместо схемы (3 !76) не во всех точках только потому, по в ней для вычислений и логических операторов 1Р ири реалнзапии программы на Фортране требуются несколько большие затраты машинного времени, чем в схеме (3.!76). и тогда 1 + = ий 1, — и9 ),, — ( ! и, ( йс + ! и,+, ~~се Д.

(3.207) с=с, Рассмотрим случай, когда и не меняет знака ') при переходе от точки ! к точке !+1. Тогда член, соответствующий ошибке, можно рассматривать как искусственный сток функции с. Если же и, ( О, а ссьм > О, то этот член можно рассматривать как искусственный источник " з). Так как на участке между точками ! и ! + 1 существует точка, где и = О, то и, н ись! обычно малы и при уменьшении шага пространственной сетки ошибка, обусловленная неконсервативностью, стремится к нулю. В двух модификациях схемы с разностями против потока можно устранить появление искусственного источника. Первая модификация основана на первой схеме с разностями против потока (3.176), когда скорость не меняет знак ни между точками ! — 1 и !', нн между точками ! и !+ 1.

Если же скорость меняет здесь знак, то конечно-разностная схема строится при помощи метода контрольного объема, охватывающего точку с'). Итак, мы имеем Л". с. +! + (3. 208) ЛС Лх а!.!!. Вторая схема с разностями против патока !!3 Легко проверить, что такая модификация делает первую схему с разностями против потока консервативной и транспоргивной.

При этой модификации требуется меньше арифметических операций, но по точности она уступает второй схеме с разностями против потока, к обсуждению которой мы теперь приступим. 3.1.11. Вторая схема с разностями против потока Во второй схеме с разностями против потока, называемой также схемой с донорными ячейками (Джентри, Мартин и Дали [1966) ), по каждую сторону от узловой точки пространственной сетки находятся некоторые средние значения скоростей на границах ячейки; знак этих скоростей определяет, из какого именно узла сетки надо взять значения ь для написания разностей против потока.

В одномерном случае будем иметь Вь, аяз — и й (3. 2! 2) Ьх где иа = '!а(и,+! + и), ис = !!а (и! + и! !) (3.213) (возможны и некоторые другие способы осреднения). Значения ь берутся так: при ил~О, 1а= ! ьт+! при ия ( О, при иь) О, при иь < О. Легко проверить, что эта схема является консервативной и транспортивной. Эту схему просто интерпретировать с точки зрения метода контрольного объема, если величины скоростей на границах ячеек находятся как средние значения, а соответствующие величины ~ определяются направлением потока.

(Замечание. Если величины с на сторонах ячеек определять тоже как средние, то получится схема с разностями вперед по времени н центральными разностями по пространственным переменным, не обладающая свойством транспортнвности.) По сравнению с первой схемой с разностями против потока в рассматриваемой схеме требуется проведение дополнительных вычислений для скоростей: дополнительное численное дифференцирование функции тока ф для получения и и дополнительные расчеты средних значений (3.213).

Однако эта схема точнее первой схемы при аппроксимации производной д(ив)/дх, так ЗЛ. Методы решения уравнения переноса вихря 114 как сохраняет кое.что от второго порядка точности, которым обладают схемы с центральными разностямн. Рассмотрим случай, когда ~ постоянно, т. е. = ~;е~ = ~, а и является функцией пространственной переменной.

Тогда уравнение (3.212) принимает вид (3 215 ) Ьх (3.215 а) инь — и й бс Ьх или Ль, и,,— и, сн. с бс 2дх (3.215б) би,. Зà — 4Ь,, + Ьш (3. 216) бх 2дх что дает второй порядок точности для члена, описывающего конвекцию. Таким образом, вторая схема с разностями против потока обладает как свойством консервативности, так и свойством транспортивности и сохраняет кое-что от второго порядка точности, присущего схемам с центральными разностями для пространственных производных. Превосходство второй схемы над первой схемой с разностями против потока было фактически продемонстрировано Торрансом 11968) при расчете плоских течений.

Теперь можно объяснить удивительное согласование результатов (при Рее = 100 отклонение менее 5осв) решения задачи о течении вязкой жидкости внутри замкнутой прямоугольной области с одной подвижной границей, полученных прн помощи этой второй схемы и при помощи схемы второго порядка точности, использованной Торрансом. Вдоль средней части стенок, где применимо приближение пограничного слоя, влияние члена со схемной искусственной вязкостью ае мало (см.

обсуждение в разд. 3.1.8). Вблизи же углов скорости малы, поэтому здесь коэффициент физической вязкости сс» ссе. А во вращающейся центральной части области течения вихрь ~ меняется слабо, и поэтому рассматриваемая вторая схема имеет почти второй порядок точности в соответствии с уравнением (3.2!5б), Схема Ранчела с соавторами 11969) алгебрапчески эквивалентна второй схеме с разностями против потока.

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

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

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

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