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

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

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

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

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

Теперь возникает следуюшнй вопрос: что будет происходить в том случае, когда механизм диссипации для отбора энергии от коротковолновых компонент отсутствует и когда при расчетах нельзя различить компоненты с длиной волны Л ( 2Ьх? Неожиданный ответ заключается в том, что тогда энергия пере- распределяется и вновь переходит к длннноволновым компонентам, искажая эти представляюшие наибольший интерес компоненты и даже приводя к некоторому роду численной неустойчивости (Филлипс [1959[). Наличие в расчетах какого-либо затухания, физического нли численного, обеспечивает диссипативный механизм н тем самым облегчает решение этого вопроса.

Таким образом, расчеты течений при малых Ке н расчеты течений невязкой жидкости прн помощи численных методов, вводящих искусственное затухание коротковолновых компонент, меньше страдают из-за ошибок, обусловленных неразличимостью, Слово нерпзличцмость используется также для описания того факта, что на конечной сетке просто нельзя добиться разлнчимости некоторых частот. (Именно из-за этого энергия может перераспределяться к длинноволновым компонентам,) То обстоятельство, что два таких распределения, как = соз [п(пт.+ п)1[, имеют в точности одинаковые значения в узловых точках ~', является следствием простых тригонометрических тождеств. Такая же дискретизация или «сеточное» пред- а!. Методы решения уравнения переноса вихря 196 ставление янления во времени вызывает стробоскопические эффекты в движущихся картинах, например колеса вагона кажутся замедляющими движение и вращающимися в обратном направлении.

Такое смещение частот является неизбежным следствием дискретизации; см. Хемминг [1962, с. 276, 303]. В заключение этого раздела рассмотрим обобщение схемы Лейта на случай двух пространственных переменных. Очевидным обобщением уравнения (3.224) будет уравнение ~»+1 ~л С (9» !» ) + Се (~л 2~» + ~л ) С (~л ~л! ) + Сз (Ц 2~~ + ~» ) (3 252) или в очевидных сокращенных обозначениях (ьл) 1 1 Сзбз (~~) С 6 (йл) 1 ! Сзбз щ (3.253) где С, = иИ/Ьх, а С„= об!/Ьу. Эта схема не является схемой второго порядка точности.

Схема второго порядка должна была бы включать член с производной дзь/дхду в разложениях в ряды Тейлора. Более того, данная схема неустойчива (Лейт [1965]). Таким образом, конечно-разностная схема (3.252) является примером двумерной неустойчивой схемы, полученной в результате комбинации двух одномерных схем, каждая из которых устойчива. Лейт [1965] дает устойчивую двумерную схему, которая основывается на понятии дробных шагов по времени (Марчук [1965]) и которую теперь обычно называют схемой расщепления по времени. Ее идея заключается в последовательном применении каждой одномерной схемы по отдельности, причем результаты, получающиеся на первом промежуточном шаге, лишены фязического смысла '). Эта схема записывается так; ьз+,'" = ~! г — 9 С,б„(ь") + — С',б'-„'(ь"), (3.254а) р+! ~л+!н С б влети) + Сзбз (~ле!н) (3.254б) Она является аппроксимирующей, в отличие от схемы (3.253) обеспечивает полный второй порядок точности и устойчива при условиях С, ( 1, Св ( 1.

') Эта конпепния похожа на расшеплеиие, ранее использованное Пасменом и Ракфордом [1955) в неявной схеме метода чередующихся направлений (см. равд 3.1 !6), но не зквнвзлспгна ему. В монографии Яненко [!967) рассматриваются н другие методы дробных шагов по времени.

127 3 l 13 Схема Лейта По-видимому, любую устойчивую одномерную схему можно применять и в случае двух пространственных переменных, когда проводится расщепление по времени, причем условия устойчи. ности для одномерной схемы не меняются. В задачах обтекания тел такое расщепление по времени приводит к трудностям, связанным с граничными условиями на промежуточном шаге Как мы увидим, для вычисления граничных значений гь На поверхности тела используются значения функции тока зр во внутренних точках, но вычислять значения фату на промежуточном шаге, кажется, не имеет смысла. При достаточно малых шагах по времени на втором шаге схемы (3.254), вероятно, можно брать значение ~н с предшествующего слоя по времени. Как влияет такая постановка граничного условия на устойчивость и точность какой-либо схемы расщепления по времени, пока что не установлено.

(При решении метеорологических задач, рассматривавшихся Лейтом, трудностей не возникает.) Фромму [!97!] при помощи введения в одношаговую схему членов со смешанными производными удалось добиться устойчивости в схеме, которая была устойчива в одномерном случае и неустойчива в двумерном. !засщепление по времени может привести помимо улучшения устойчивости и к улучшению точности (см.

следующее упражнение). Упразхнение. Применить схему с разностями против потока в случае двух пространственных переменных и отсутствия вязкости для уравнения переноса с постоянными коэффициентами дь дь дь — = — и — — о —. д! дх ду ' Показать, что схема расщсплснвя по времени в отличие от схемы без расщепления приводит к точному репгсни~о, котла числа Куранта ио1/дх = ) и па1/Ьу = 1. Некоторые схемы, аналогичные схеме Лейта, обсуждаются в работах Касахары [1965) н Фишера [1965а).

Эти схемы и схема Лейта похожи на схему Лакеи — Вендроффа и ее двухшаговые варианты, которые будут рассматриваться в гл, 5, Хотя схемы Лакса — Вендроффа были разработаны для течений сжимаемой жидкости, они представляют интерес и для течений несжимаемой жидкости (Липли [1965) ), несмотря на то что эти схемы приводят к сильному затуханию коротковолновых возмущений.

Другие схемы, построенные на идее расщепления по времени и обладающие лучшими дисперсионными свойствами, будут рассматриваться в следующих разделах, а сейчас мы продолжим обсуждение некоторых более простых схем. 8.1 Метет!ы решения уравнения переноса вихря 128 3.1.14. Неявные схемы Рассмотренные выше схемы являются явными, т.

е. в них для вычисления значений на (и+ 1)-м слое по времени необходимы только известные значения на п-м, (и — 1)-м, ... слоях, Теперь приступим к обсуждению неявных схем, в которых в простраис!венных производных используются значения на п + ! слое по времени и поэтому для продвижения расчета нужно одновременно решать систему уравнений на и+ 1 слое. Запишем общую схему для модельного уравнения, описывающего течение повязкой жидкости, в следующем виде: Г."" = 1" — ( Л!) бих.

(3.255) Здесь будет рассматриваться представление пространственной производной только центральными разностями, так что Ьей !.т ! (3. 256) Ьх 2 ах где слой по времени н пока еще не определен. Аналогично, для уравнения диффузии ноложим Ь" ' ' = Ь" + (а Л1) б 8/бх (3.257) Если пространственные производные в (3.255) или (3.257) выписать на п-м слое по времени, то получится рассмотренная ранее явная схема с разностями вперед по времени и центральными разностями по пространственной переменной и с ошибкоч порядка 0(ЛА Лх').

Если же член 6~4х в уравнении конвекции (3.255) записать на новом (п+ 1)-м слое, то получится так называемая полностью неявная схема ьп-т тп ( А1) был~ /, ° (3.258) Ошибка при этом по-прежнему имеет порядок 0(ЛАЛхе), но эта схема обладает существенным преимуществом в смысле устойчивости. Исследуя устойчивость методом фон Неймана и полагая С = и А!1!хх, получаем (т +! )т (С)2) (тп+! ( те,— те) )т"" ' (1 + С7 гбп 6) = (т", 6— 1 1 — С!е!пв 1+ С! Мп О 1+ С'Мп'О 1+ Сеып'О ! [! + С'Мпевр 1+ Се Мп'О ' (3.262) Таким образом, для полностью неявной схемы имеем ) 6! ( 1 независимо от величины С.

Данная схема абсолютно устойчива, Знв!4. Неяаньге схемы 129 что дает возможность вести расчеты с произвольно большим шагом по времени, а это является большим преимуществом. Полностью неявная схема первого порядка точности абсо- лютно устойчива также и для уравнения диффузии') (Лаасо- нен [1949], Рнхтмайер н Мортон [1967]), При с( = схЛ1/Лхз по- лучаем ~" е~ = ~" + а Л!б ~~~ ('бх, (Ун+! Рч ! )гне!д!( ш 1,-'" 2) (г" е' [1 — 2а' (соз Π— 1)] = 1', 6= 1 !+2а(1 — созз) ' (3.263) (3. 264) (3. 265) (3.266) Так как 1 — сов О ) О для любых О, имеем /6] ( 1 для любых г( или любых Лй Заметим также, что 6 ) О для любых О; как было указано ранее, это условие отвечает случаю, когда осцнлляции, обусловленные чрезмерно большим шагом по времени, отсутствуют.

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

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

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

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