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

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

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

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

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

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

6.1. Специальные расчетные сетки Самое простое видоизменение прямоугольной расчетной сетки получается при изменении шага сетки в одном направлении в определенной узловой точке. Как правило, это делается для получения более высокой разрешающей способности сетки (н по возможности более точного решения) в той области, где градиенты параметров потока изменяются быстрее, например в пограничном слое. Для иллюстрации сказанного рассмотрим простейший способ перехода от шага Лх1 между узламп сетки к шагу Лх, в некотором узле ~ = т (рис. б.), а). 426 б.б Спеииальяые расчетные сетки Заметим, что если величина Лхз/Лхг С< 1, то точности представления в точке пт ухудшается до первого порядка малости относительно Лх, (см., например, Блоттнер и Роуч 1!9711).

Выражение для второй производной получается умножением равенства (6.2) на з' = (Лхз/Лхг)' и сложения полученного результата с равенством (6.1): зг) ( — Лх, (1 — з) + —, ~ Лхз+ — —.„~ Лх-(Лх — Лх,) + д( ~ + О (Лх'), (6.6) + О ПЛхз бх' Лх, бх ~ Лх — Лх,), Лхз), (6,7) Для того чтобы последнее выражение имело первый порядок точности в точке г' = ьн должно выполняться равенство з = = О(1 — Лх'-,). Сальвадори и Барон [1961] провели через точка т, гп+ 1 и пч — 1 параболу и получили следующее выражение в конеч. йых разностях: бз( ) 2 з)„,з ~ — (1+ з) ( +(м 1т з(з+1) Лх' (6.

8) 2 Соответствующее выражение для первой производной дается формулой (6.5). Причину более высокого порядка ошибки аппроксимации этих формул легко объяснить при помощи рассмотрения контрольного объема, как зто требуется сделать в следующем упражнении. Упражнение. Проведя границы ячеек между узловыми точками рис. 6.1, и, показать, что при з « 1 точка т удалена от центра ячейки. При з -ч 0 точка гп приближается к правой границе ячейки Можно несколько улучшить физическую интерпретацию, однократно меняя не шаг сетки Лх, а размер ячеек Лх, как показано па рис.

6,1, б. При этом нужны выражения для производных как в точке т, так и в точке (т + 1) = и. Эти выражения аналогичны уже полученным. Упражнение Показать, что полученаые выше формулы для б(гбх и бз()бхз консерватннны Из приведенных выше выражений следуе~, что при быстром изменении шага сетки формальный порядок ошибки аппрокси- В д Спсцпальиьге расчетные сетки 427 мации на самом деле не улучшается, а ухудшается ').

Так, Краудер и Дальтон [1969) провели численный эксперимент с пятью различными сетками с локально меняющимся шагом и пришли к выводу, что для рассматриваемой нми частной задачи наиболее точные результаты дает расчетная сетка с постоянным шагом (см., однако, Блоттнер и Роуч [1971]). Напомним также о явлении отражения волн от места изменения шага сетки, упоминавшемся в равд.

5.4,4, Однако при расчетах было обнаружено, что потеря точности всего решения, особенно при изолированном изменении шага сетки, не столь велика, как это следует из формулы для ошибки аппроксимации (см., например, Мак-Кормак [1971], Шаве и Ричардс [!970), Магнус и Иосихара [1970) ). Хотя вопрос окончательно пе исследован, по уже ясно, что обычно предпочтительный метод повышения локальной разрешающей способности сеткц зак.почается в преобразовании координат, которое будет рассматриваться в равд.

6.2. Но прежде мы рассмотрим еще ряд вопросов, связанных с расчетнымп сетками. Геометрические схемы для одновременного изменения величин обоих шагов Лх и Лу рассматривались в работах Синнота [1960], Ранчела с соавторами [1969], Гиллнса п Лиропа [1969», Кейкера и Уайтло [1970), Даусона и Маркуса [!970). Бахвалов [1959] использовал мелкую сетку с формулами второго порядка у границ, стыкующуюся с грубой сеткой с формулами более высокого порядка во внутренних точках. Преимущества однократного изменения размера ячеек по сравнению с непрерывным изменением обсуждались в работе Робертса [1971). Проблемы стыковки сеток для схемы Лакса — Вендроффа (равд.

5.5.5) рассматривал Саймент [1968). Спеньер [1967] рассматривал методы чередующихся направлений (равд. 3.1.16 и 3.2.6) на сетках с меняющимся шагом. Бердсли [1971] применял полярную сетку с постоянным шагом по углу н с параболическим изменением шага по радиусу; в случае невязкой жидкости для решения уравнения вихря около точки г = 0 здесь потребовались некоторые модификации. Изменение пространственной сетки для решения, развивающегося во времени, часто оказывается желательным для более подробного описания областей с резко меняющимися в пространстве градиентами. Преобразование решения уравнений из ') Схемы коиечио-разпоствого представдсзшя коивсктивиых члсиов, ко.

торые имеют первый порядок точности па равиоисрвых расчетвых сетках (коиечиые разности против потока). ае ухудшаются при измеиеиии шага сетки, в то время как схемы представлеиия диффузионных членов ухудшаются. Для получения второго порядка точности требуются четырехточечиые формулы (см. Саусвслд [)946)). Пирсон [!9681 использовал тректочечиые формулы иа сетках с автоматическв изменяющимся шагом дая квазиод. номерного расчета распростраиеиия ударной волны. 428 б.!. Сггеииальные рпсчетпме сетки одной сетки в другую называется перестройкой ячеек.

Перестройка ячеек сама по себе может изменить решение, внося, например, сглаживающий эффект (в некотором роде искусственная диффузия) пли ошибки, связанные с нарушением консервативности. Разработка машинных про~рами, осуществляющих перестройку ячеек в зависимости от развития во времени решения, является важной н интересной проблемой (см. в этой связи работы Месона и Торна [1970], Батлера [1971] и Кроули [1971]). Рис.

82. Точки нерегулярной гранины на прямоугольной сетке. Приведенные выше формулы, применимые в случае переменных стх и Лу, использовались также для описания нерегулярных границ на прямоугольной сетке. Как показывает рпс. 6.2, для того чтобы значения в смежной с границей точке (й /) согласовывались с значениями в точках Ь1 и (гт, требуются формулы с меняющимися величинами гхх и с!у. Эта процедура предложена давно и иногда дает, по-виднмому, вполне удовлетворительные результаты (см., например, работы Сальвадорн и Барона ]!961], Томана и Шевчика [1966], Синглтона ]1968], Техейры [1966] н Даусона и Маркуса ]1970] для сеток в декартовых координатах, а также работу Лизена [1964] для цилиндрической функции тока при нерегулярной границе).

Однако мы настойчиво рекомендуем воздерживаться от этой процедуры по следующим причинам, (!) Формальный порядок ошибки аппроксимации, как показано выше, увеличнваезся. Заметим, что для произвольной искривленной границы типа крылового профиля (Синглтон [1968]) н общем случае в равномерной сетке может обнаружиться не- 429 Дк паениаесные расчетные сесна которая виутреиияя узловая точка, располо>ке>шая очень близко к границе; при этом величины сГхе/с>х> и Луч/Лд> становятся весьма малыми и, слсдователыю, формальиый порядок ошибки аппроксимации сушествеино возрастает. (2) Скорость сходимости уменьшается (Техейра [19661), и часто бывает трудно определить оптпмальиые релаксацпоипые параметры (разд. 3.2.4, 3.2.6, 3.2.7).

(3) Условия устойчивости могут стать весьма жесткими. Заметим, что обычпое ограничение по числу Кураита дает неравенство с» ( Лх/и (или его двумерный аналог). Голи уменьшать величипу с>х возле граиицы, то максимум величины 61 может быть ограипчеп этим локальным условием. Если имеются границы с условиями прплппаипя, то ограничение, накладываемое иа коивективиые члены, может и пе быть чрезмерным, так как величипа и такнсе локально мала; однако для вязких членов и для границ со скольжением эти локальные условия, по-видимому, будут играть решающую роль. (4) Программа становится более сложной, и поэтому увеличивается как затрачиваемое па ее составление время, так и вероятность появления ошибок. Проблемы аппроксимациоииой сходимости решепий эллиптических уравнений при перегуляриых граипцах иа прямоугольных сетках обсуждались в работах Турайсами (1969а, 19б96).

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

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

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

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