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

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

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

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

Расчет давления на поверхности удобно проводить при помощи легко выводимого соотношения (см. задачу 3.32) 6Р 1 е)ь оз ке дю ' (3.511) З.б.2. Уравнение Пуассона для давления Более точное решение можно найти, интегрируя уравнение Пуассона для давления, которос получается следующим образом. ') Такое интегрирование известной функции [(х, д) по двум везависимым переменным Бакпнгем [1962) называет «кубатурой». где п и з — нормаль и касательная к поверхности (см. Пирсон [1965а] ). Полное описание трудностей, возникаю1цих при численной реализации этого способа, можно найти в работе Масляха и Эпштейна [1970].

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

Этот способ особенно чувствителен к ошибкам в задачах, подобных задаче, прсдставлснной на рис. 3.22, когда путь интегрирования проходит вблизи от угловой точки. Дд 7 Уравнение Пдаггана дяя даамния 277 Продифференцируем уравнение (3.509а) по х; д'Р ! д = — — —, + — — (Ч'и). дх' Ке дх (3.512) Продифференцируем уравнение (3.5095) по у: дзР 1 д = — — + — — Ж'о). ддз Ке дд Сложение уравнений (3.512) и (3.513) дает Из уравнения неразрывности (3.509в) имеем (3.513) (3.514) (3.515) кроме того (3.516) (3.

5! 7) (3. 518) прини- (3. 519) Рассмотрим тождество (3. 520) откуда (3.521) Наконец, получаем д а д ., дзи д'и дз„ дза дх ( ду ( ) дх' + дхдаз + дз дхз + диз С учетом равенств (3.515) — (3.518) уравнение (3.514) мает вид д 5 Рат«ет давления 276 Подставляя этот результат в уравнение (3.519), получаем 1 , до ди ди до 2 дх ду дх ду ' Чзр (3.522) (Из этого уравнения видно, что иа стенке с условием прилппания Ч'Р = О.) Для того чтобы записать уравнение (3.522) через функцию тока ф, вспомним, что и = дну, о = — дфудх; (3. 523) тогда — 2ЧзР=( д е')(д з) (д д )( — д д ), (3.524) или (3.525) Это уравнение Пуассона для давления.

Обычные конечно-разностные формулы с центральными разностями по пространственным переменным (см. равд. 3.1.1) приводят к следующему выражению: 5. =2Гф+«1+ф; — ь1 — эфь1 Фьуы+Ф~ 1,— 24ьт о 1=2( д„з' а„а ( ' ' ' ' ) ), (3.526 ) 1 (3 526) «зч 4 ахну имеющему порядок точности 0(2«хз, Луа). Уравнение (3.525) аналогично уравнению для функции тока Чз«р = ~, причем источниковый член о аналогичен источпиковому члену ~.

Все методы решения уравнения Пуассона, которые были рассмотрены в равд. 3.2, применимы и в данном случае, но как будет выяснено в следующем разделе, граничные условия в этих двух случаях различны. З.б.З. Граничные условия второго рода для давления При решении уравнения Чзф = ~ хотя бы некоторые граничные условия являются условиями первого рода, или граничными условиями Дирпхле: ф(х,у) есть заданная функция вдоль границ. При решении уравнения Пуассона для давления ставятся условия второго рода (граничные условия Неймана), т. е, на границе задается дР(х,у)1дл, где и — нормаль к границе. Величина градиента давления находится пз уравнения (3.509). Упражнение. Запасать ионечно-разностное представление для составляюпгих градиента давления бр(бх и бр/бу через функцию тока ф, Упражнение. Показать, что для стенки с условием прилнпання (Пирсон (1966а]) бр 1 бь (3.627) бп Ке бз Е9 5.5.4, Игераиионньге метода решения где л и з — нормаль и касательная н стенке.

Производная по касательной к стенке б~Яз в уравнении (3.527) легко вычисляется с помощью центральных разностей. Заметим, что здесь при расчете можно брать любой метод, поскольку при решении уравнения Пуассона градиент не пересчитывается, а отсутствие обратной связи исключает возможность появления неустойчивости. Однако предпочтение отдается схемам с центральными разностями, согласующимися со вторым порядком точности 0(ЛЯ) метода в целом.

Для нестационарных решений в уравпепнях (3.509) часто пренебрегают членами ди/д1 и до/д1 и используют простейшее приближение дР,'дп ж 0'). Эти приближенные допущения не обязательны, и следует решать уравнения (3.509) точнее. Заметим, что наличие члена с коэффициентом 1/Гсе отнюдь не означает, что в течениях с большими Ке этим членом можно пренебрегать. Для пограничного слоя, толщина которого б мала, известно, что дР(дп = = 0(б), но в более общем случае величина дР(дп может оказаться большой. (На линии симметрии, конечно, ставится условие бР!би = 0.) В окрестности угловой точки С (рис. 3.22) постановка граничных условий Неймана может привести к двузначности давления в этой точке. Такую двузначность, как и двузначность функции с, в способе ! (рнс.

3.30), рекомендуется сохранить, хотя это ие физично. Ошибка же в «однозначности» Р. — Р, может являться мерой ошибки аппроксимации вблизи угловой точки З.б.4. Итерационные методы решения Граничные условия Неймана выдвигают два специальных требования при решении задачи. Первое требование заключается во введении градиентного условия в уравнения метода последовательной верхней релаксации. Очевидный способ решения здесь таков: вычисляются новые значения на (й + 1)-й итерации во всех внутренних точках сетки, а затем по известной величине бРЯп и вновь вычисленным значениям в точках, смежных с границей, рассчитываются значения функции па этой границе. Для точкн (5 /с) границы В 2 (рис.

3.22), используя в соответствии с формулой (3.380) метод последовательной верхней релаксации, получаем следующие уравнения: ') Это приближение известно иа теории иограничного слоя (см., наиример, Шликтннг (!968]). Фактически е рассматринаемом случае ириближенне будет более точным, поскольку постоянстао Р иредиолагается не поперек коего чогранн мого г ~оя, а лишь а его части. б.тижзйшей к стенке д Д Рас1ет давления 280 для расчета внутренних точек »-1-1 ы гРс Рз+! 2 З с2ра Рс (сь! З (! ) й ) ~~!е! (с»1+ Р1-1, (се!+ Р Р1 ! ьз+ н Рс, — лх Я! (, ! — 2(1+ 5 ) Р! ! 1)+ )т~ ),, (3»528а) (здесь Р,' , — значение на границе, взятое с Ф-й итерации); для определения значений функции на границе 6Р Рс „=-Рь„„— —,Лр.

(3.528б) Однако такой правдоподобный способ не сходится. Решение «ползет» медленно, но нео!раниченно долго. Это является причиной того, что вплоть до конца шестидесятых годов большинство опубликованных работ не содержало расчетов давления, а приводившиеся в остальных работах распределения давления оказывались неверными, хотя они и удовлетворяли какому-либо «критерию сходимости», подобному условию (3.501). Метеоролог Миякода 11962] рекомендует подставлять градиентные граничные условия непосредственно в разностную схему метода последовательной верхней релаксации при расчете внутренних точек, смежных с границами'). Таким образом, уравнение в виде (3.528а) бере~ся только во внутренних точках, отстоящих от границ более чем на одну ячейку.

В точках, смежных с границей, уравнение (3.528а) заменяется следуюшнм: ы г Ь -1- ! 2 1, в+1= с(! ( 1)2) ЬР!+1, !!11+ 1 1 1.(с+1 +(с Р1, !с!2+ 2! 2,1 бР( ! 2 + 11 (кР! '(се! — — ~ ЛР ) — Лх 5! )с~.!— ба !1, (с — 2(1+ 8 ) Р1, (с»11+ Р1, !с+ь (3.529) Второе требование, связанное с постановкой условия Неймана, заключается в том, что градиент давления на границе должен согласовываться с источниковым членом в уравнении Пуассона. По теореме Грина для существования решения диф- 1) Автор весьма нризнатеаен н-ру С. Пнанеку, который указзл ему зту ненную работу. которое разрешается алгебраически относительно члена Рзт! входящего н обе части равенства. После того как сходнмость достигнута, окончательные значения функции на границе можно найти по уравнению (3.528б).

Уравнение (3.529) отличается от (3.528а) н (3.528б) только номером слоя по времени в члене 52(Рс".е!1ь! — ЛР)бп 1! ), Лу). 5.5.5 Характерное велнчина длл отсчета давление 281 ференциального уравнения (3.525) в области Е необходимо выполнение условия Е = ~ Е с!Š— ~ (дР7дп) Ж = О (относительно и ел одномерного случая см, приложение А). Из-за ошибок аппроксимации значения функции на границе обычно не удовлетворяют этому условию, что приводит к медленной расходимости метода последовательной верхней релаксации. Миякода 11962) рекомендует задавать ЬР7бп так, чтобы это условие выполнялось.

Способ, предложенный Брили 1! 974), а также Гхиэ и Гхиа (личиое сообщение), заключается в нахождении дискретизированной величины Е и последующем решении видоизмененного уравнения К'Р = 5 — Е!К где !с — площадь рассчитываемой области.

Если интегрирование для величины Е выполнено надлежащим способом (Гхпа и Гхна; см., например, ниже уравнения (3.533) н (3.534)), то итерационная сходимость как в методе последовательной верхней релаксации, так и в неявной схеме метода чередующихся направлений существенно улучшится. Отметим, что градиентное граничное условие второго рода может потребовать специального подхода, чтобы избежать неопределенности в методе прогонки, применяемом в неявной схеме метода чередующихся направлений (см. приложение А).

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

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

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

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