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

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

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

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

В опубликованных работах по вычислительной гидродинамике такой подход не использовался. Важное замечание о аппроксимационной сходнмостн сделал Чен [1968, 1970]. Даже если и удается найти предел 1!т 1;(Л) л-»о или добиться выполнения равенства ь(Лз) = ь(Л!)+и при произвольно малом е, для задач внешнего обтекания может проявиться эффект вычислительных граничных условий, Строго говоря, необходимо также проверять сходимость в интересующей нас области, когда «бесконечно удаленные» границы (верхняя, нижняя, входная и выходная) отодвигаются от этой области. Проверка сходимости только по шагу сетки или при помощи методов более высокого порядка не дает информации для установления сходимости в смысле уменьшения влияния гра ннц.

В этой связи стоит упомянуть работу Гамилеца и Рааля [1969]. Они решали задачу об обтекании кругового цилиндра в полярпой системе координат с разностной сеткой, постираю- щейся до т = гь! все граничные условия, соответствующце «бесконечно удаленной» границе, ставились на окружности такого ') дппрокснмапиопная скодимость с достоверностью монотонна для уравнения Лапласа !7~ф = 0 и лля простого уравнения диффузии дсгщ = = ывтйг, но не обязател~ но для уравнения Пуассона влн для нелинейного уравнения с копвсктпвнымп и диффузионными членамн другие ссылки по приме. нсюпо экстраполяпоонных ьчетодов для дпфференпоальных уравнений в частных провзводных можно найти в работе Берджеса [!97!1 и в равд 59 настоящей книги.

Шенхерр и Черчилл !!9701 рассмотпеля экгтрвполяпию вестапионарных рещеннй к стапионарному пределу для >равнг пня диффузии; Штеттер (!9701 рассмотрел экстраполяпию Ричардсона для жестких уравнений. 8.4. Критерии сеодимости и начальные условие 273 радиуса. Расчеты проводились при двух значениях гь; предполагая, что скорость сходимости подчиняется квадратичному закону, эти авторы оценили коэффициент сопротивления кругового цилиндра при помощи экстраполяционного метода Ричардсона, пРоводЯ экстРаполЯцию до 1ььгь = О. Другую ннформацию об ошибках аппроксимации можно полуььььть 1) вычисляя ошибки, связанные с нарушением свойства консервативности (см.

равд. ЗЛ.З) для неконсервативных схем, 2) вычисляя контурный интеграл от дсььдп по кривой, охватывающей границу тела (см. задачу 3.32), и 3) сравнивая два значения для давления в угловой точке контура тела (см. равд. 3.5.2). Если для нахождения стационарного решения полной задачи о течении несжимаемой жидкости применяются нестационарные двухслойные схемы и если в итерационном процессе для нахож„.ь.ь ач-ь дения т(ьа" ь и (~ ) используются неявные достаточно хорошо сходящиеся ') методы, то начальные условия будут несущественны. (Начальные условия, очевидно, определяют начальное не- стационарное решение, но даже это влияние экспоненциально затухает по времени, см., например, Бердсли [19691.) Далее, начальные условия обычно ие оказывают существенного влияния на требуемое для расчета машинное время, хотя многие исследователи предполагают противоположное. Начальные условия слабо влияют на требуемое машинное время, поскольку ошибки для выбираемого начального приближения обычно ограничены и по величине на много порядков превосходят величину, входящую в критерий сходимости.

Например, если безразмерный расход в канале нормирован так, что т)ь,„= = 0(1), то начальное приближение т)ьо = 0 во всех внутренних точках будет давать ошибку только порядка О (1). Очень хорошее начальное приближение для поля течения с отрывом в расширяющемся канале может давать ошибки в величине т)ьо порядка 0(10 '), Однако если в критерии сходимости взять величину в = 1Π— в, то улучшение будет незначительным. Для лучшего понимания процесса сходимости следует помнить, что при расчете течений с большими числами 11с ошибка в значении вихря ~ в единственной узловой точке вблизи входной границы и вне пограничного слоя должна за счет коивекции уноситься из рассчитываемой области.

Тогда нужное для сходпмости число шагов по времени будет ограничено снизу величиной и = т/Л/, где т — время переноса частицы нз указанного положения через всю рассчитываемую область. Это ') Есльь вто не так, то устойчивость метала в пелом можно повысить ва с ьет умеиьшеьиья ЛЬ ня начальной ставни расчета (Бао и Догерти 1!9691). 274 Д4, Критерии сяодимости и на«альньы услоаия время т не зависит от размера ошибки и поэтому приближенно представляет собой время, нужное для сходимости, В качестве характерного примера несущественности начальных условий рассмотрим задачу об обтекании обратного уступа (рис.

3.22). Автор данной монографии решал эту задачу, принимая в качестве начальных условий т)т = 0 во всех внутренних точках н вдоль границы В 1 — В 5 — В 2, задавая на входной границе значение ар, соответствующее течению в пограничном слое, и считая, что граница В 3 является «крышкой>, т. е. тр(В 3) = = а)т(1,2); для вихря всюду полагалось ь = О. Такое начальное приближение кажется совершенно неразумным. Однако после первой итерации прн решении уравнения Пуассона с граничными условиями на входной границе, заданными по формулам (3.473), всюду появилась отличная от нуля скорость конвекции.

К моменту п = 30 формировалась вполне правдоподобная зона возвратно-циркуляционного течения, а это указывало на то, что начальное приближение оказалось лучше, чем можно было ожидать. При таком грубом подходе для окончательной сходи- мости при Ее ) 1 потребовалось такое жс машинное время, как н при общепринятом подходе, заключающемся в расчете очередного варианта при начальном приближении, взятом по результатам предыдущего варианта, полученным при ином це или иных условиях на входной границе. Начальные условия могут оказывать влияние в методах расчета стационарных течений при использовании «непрерывных замещений» при задании начального приближения в нестационарных методах, если уравнение Пуассона решается с недостаточной степенью точности, и в неявных методах, если граничное условие на стенке (с".+')ь+' недостаточно проптерировано.

В этих случаях плохие начальные условия могут привести к неустойчивости, связанной с нелинейностью уравнений. (В двух последних случаях неустойчивость можно предотвратить уменьшением М на начальной стадии расчета.) Даже в случае простейших уравнений, когда они решаются при помощи многослойных схем, начальные условия могут вызвать возникновение лишенных смысла осцнлляций. Поскольку данные на предыдущем слое по времени иногда могут быть очень хорошим начальным приближением для последующего слоя, точность начального условия может также оказать существенное влияние на итерируемое решение уравнения Пуассона. Сошедшееся решение аРамт на новом слое по времени отличается от тря только за счет источникового члена ~ н за счет граничных условий для Чт Если Л1 достаточно мало или если сходимость по времени нестационарной задачи в целом почти достигнута, то са+' СЯ и ьР"Ы ж ф" на границах (Для яадаст, подобных задаче о течении внутри замкнутой прямо- У д Д Численное интегрированна угольной области с подвижной стенкой, гр фиксировано вдоль всех границ и в течение всего времени.) Тогда выбор решения грн в качестве начального приближения (!рн+!) будет очень эффективным, Действительно, если сходимость для вихря почти достигнута, то для функции тока ф критерий сходимости (гу'+! — гр'(( е часто удовлетворяется за одну итерацию, так что применяемый метод автоматически сводится к непрерывным замен!синям, как в случае расче~а стационарных течений с помощью итераций.

Линч и Райс (1968) показали, что для неявных схем метода чередующихся направлений скорость сходнмосп! выше в случае гладких ошибок в начальных данных. На практике это обычно выполняется, поскольку гладкими являются и окончательное решение, и начальное приближение (включая случай чг = О). Точность начальных условий может играть более важную роль в случае сверхзвуковых течений, где плохие начальные данные могут привести к распространению паразитных волн.

З.б. Расчет давления ЗЛэ.!. Численное интегрирование для определения давления Уравнения (2,1) — (2.3), записанные в безразмерных пере- менных (см. первое упражнение в равд. 2.4), имеют вид ди ди ди дР ! — +и — + о — = — —. + — охи, дг дх ду дх Ке до до до дР ! — +и — + о — = — — + — т1О, д! дх ду ду Ке ди до — + — =О, дх ду где Р— безразмерное давление; Р = Р/(рОо). (3.510) (3.509а) (3. 509б) (3.509в) Одно из преимуществ работы с уравнениями, описывающими течение несжимаемон жидкости, заключается в том, что здесь число занисимых переменных может быть уменьшено. Давление исключаемся из уравнений количества движения в переменных (и, о, Р) при помощи перекрестного дифференцирования, как в равд.

2.2. Теперь мы будем находить поле давления по известному численному решению для ф и ~. Уравнение для давления представляет собой уравнение Пуассона и аналогично уравнению для функции тока. Однако при его решении методом последовательной верхней релаксации возникают значительные трудности из-за иного типа граничных условий, которые ставятся в этом случае. з Д Расчет дав ~енин 276 Чтобы определить поле давления, расчет можно начинать с произвольной точки, принимая в ней произвольное характерное значение Р (постоянная интегрирования); дискретизированные уравнения (3.509), содержащие дР[дх в дР]ду, интегрируются численно, причем составляющие скорости вычисляются по известному решению для функции тока. (Такой подход применяли многие авторы, в частности Кавагути [1965], Пирсон [1965а], Бургграф [1966], Хын и Макано [1966],Томан и Шевчик [!966], Римон и Чен [1969], Лавап с соавторами ]1969], Сон и Ханратти [1969], Маслях н Эпштейн [1970], Ша. вит и Лаван [1971].) Если интересуются только давлением на поверхности тела, как, например, при определении коэффициента сопротивления, то такой метод может оказаться приемлемым.

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

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

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

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