Fletcher-2-rus (1185919), страница 11

Файл №1185919 Fletcher-2-rus (Флетчер К. Вычислительные методы в динамике жидкостей) 11 страницаFletcher-2-rus (1185919) страница 112020-08-25СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

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

Использование обобщенных координат приводит к определенным трудностям. Во-первых, необходимо определить вид уравнений в обобщенных координатах. Данный аспект будет рассмотрен в 9 12.3. Численное решение типичного уравнения в частных производных будет рассмотрено в $ !2.4.

В уравнениях, записанных в криволинейных координатах, появляются дополнительные члены, определяющие отображение физической области на пространство обобщенных координат. Эти дополнительные члены (параметры преобразования) имеют форму производных, например дх/дв, и обычно для них также приходится проводить дискретизацию (9 12.2). Это является дополнительным источником ошибок в решении (п. !2.2.3). Вывод параметров преобразования приведен в $ 12.1. й 12.1. Преобразование координат В данном параграфе будет установлена связь между физическими (х,у, г) и расчетными Я,т1, ~) координатами. Обобщение на преобразования, зависящие от времени (х, у, г, 1-з-~, т1, ~, т), приведено в работах (8!едег, 1978; ТМтпрзоп е! а1., 1985] .

12.1.1. Обобщенные координаты Предполагается, что между физическими и обобщенными координатами существует взаимно однозначное соответствие, которое может быть записано в виде $ = $ (х, у, г), т! = т! (х, у, г), ~ = ~ (х, у, г) (12.1) и соответственно х=х(к, т1, ~) и т. д. Конкретные соотношения определяются при построении сетки в физической области (гл. 13). При заданных функциональных зависимостях $ = $(х, у, г) уравнения могут быть преобразованы к виду, содержащему частные производные по $, т1, ь.

62 Гл. !2. Обобщенные криволинейные коорлинаты В качестве примера можно привести первые производные от компонент скорости и, о и в по х, у и г. Они имеют вид ди дл д$ дк ди ди дЧ дь до до дп дь дв дв дп дь (12.2) дл дв да Здесь фигурирует матрица Якоби преобразования 3, которая равна да да дч дл (12.3) В принципе, если известны аналитические зависимости = й(х,у, г), элементы матрицы ) могут быть определены непосредственно. На практике явные аналитические связи не всегда известны, и более удобно работать с обратной матрнцей Якоби, определяемой выражением дк дк д$ дт! ду ду дй дт! дл дл дй дт! (!2.4) Элементы матрицы )-! можно выразить через элементы мат- рицы 1, если заметить, что -! Транспонироаанная к алгебраяческоиу пополнению Л (у '! Определитель обратного якобиана 13-! ( равен ) 3 ' ) = хй (учгс — уегч) — хч(утгт — усгй) + хс (уагч — учгт), (12.6) где хй — дх/д$ и т. д.

В случае двух переменных Д,т1) выражение (12.6) упрощается (гс —— 1, ус — — хс — — 0) и принимает вид ~ 3 ' ~ = хауч — хчуа. ди ди дх ду до до дк ду дв дв дл ду ди дй до дй дв д$ да дк дч дк д$ ду дч ду дь ду дт! дк дй дк дк дС ду дС дг дь д$ да ду дл дт! дч ду дл дС дь ду дл 63 6 12.1. Преобразование координат Используя (12.5) и (12.6), элементы матрицы ! в (12.3) можно выразить в форме кпуе — к уп $,= !з ! кгкп — кп'1 $у= Упкс — Уек, !з ! кСУ1 "1УС к1Уп — кпУФ !з ! каке — к1кп ~и )1 — 1! Уекп Упк1 !з '! где !1-'! задается выражением (12.6). После построения сетки в физической области (гл.

!3) можно определить дискретную форму элементов, например кь обратной матрицы Якоби (9 12.2). После этого для определения элементов, например $,, матрицы Якоби (12.3) используются уравнения (12.7). Это облегчает дискретизацию уравнений в обобщенных координатах, где они имеют более компактную структуру при записи их с использованием $„ и т. д., а не хе и т.

д. (12.8) Изменения физических координат Лка могут быть связаны с соответствующими изменениями обобщенных координат Л$': ЬХ~ = — Лв~ доз (12.9) (подразумевается суммирование по 1). Следовательно, малое расстояние Лз связано с обобщенными координатами соотноше- 12.1.2. Метрический тензор и физические свойства преобразования Для установления связи между обобщенными, ортогональными и конформными координатами удобно ввести в рассмотрение метрический тензор а», который связан с матрицей Якоби Л (12.3).

При этом будут использоваться тензорные обозначения (Аг(з, 1962!. Предполагается, что физическая область рассматривается в декартовых координатах х'(=х, у, з), з = 1,3, а расчетная— в обобщенных координатах кз(= — $, ть ь), 1= 1, 3. Малое расстояние Лз между двумя точками физического пространства может быть выражено через приращения координат з Ьзз = ~, Ах~Ах». а-1 64 Гл. 12. Обобшенные криволинейные координаты нием з Лаз=~~ ( —, Л5')( — ЬВ1)=узы|'Я1 (12.10) а-1 (подразумевается суммирование по 1 и /), где з г~ даа дла а» доз д$1. а-1 (12.11) Метрический тензор д» устанавливает связь малых изменений обобщенных координат Л~' с расстоянием Лз.

Более подробно метрический тензор рассматривается Арисом (Аг(з, 1962]. ззар =(жр+у, ) Ьт) 2 2122 ь лзч = (ха+ У 2 )т~дй Рис. 12.3. Физические свойства расчетной сетки. (х'+ у') (х х + у у„) $= (х х„ + у у„) (х'„ + »2) , Г (и'. + )!) — (В. . + В, ,) ( = 1 ~ - Р,ч, + ~ипа) 6„ + ~и) где, как и ранее, х =дх/д$, а 13 ( — детерминант матрицы Якоби (12.3) . В двумерном случае длины отрезков, связанные с малыми изменениями $ и т1 соответственно, равны Лз =й»2Я и Ьзч= и = д'12~Л21 (рис.

12.3). В случае двух координат соотношение (12.11) удобно представить в виде матрицы 5 12.1. Преобразование координат к1 сова= )и2 ' Отношение сторон сетки АК определяется отношением длины касательных векторов (при Ь$ = Ь21) Локальная деформация сетки определяется углом О между координатными линиями $ и 21: е= (12.11) (гнала) В трехмерном случае связь физических сеточных параметров с компонентами метрического тензора определена в работе ]Кег((ск, К1ор(ег, 1982].

(12.16) 12.1.3. Ограничения на ортогоналоные и конформные координатор Использование обобщенных координат позволяет рассматривать произвольные геометрии. Однако хорошо известно, что точность решения уменьшается при деформации сетки. Для получения высокой точности сетка должна быть ортогональной или почти ортогональной. Для ортогональных систем координат некоторые члены преобразования пропадают и уравнения упрощаются. Если, кроме того, система координат конформная, происходит дальнейшее упрощение уравнений. Использование 5 К. Флетчер. т.

2 В матричном виде метрический тензор может быть выражен через обратный якобиан (12.6): д = (Л-')тЯ-'. Определяя детерминант, можно получить ]й]и'=~3-'!. (12.13) Это легко получить, особенно в двумерном случае, прямой подстановкой. Метрический тензор у~1 и различные параметры преобразования х и т. д. могут быть связаны с физическими свойствами расчетной сетки (рис. 12.3). Здесь будет приведен ряд формул для случая двух координат.

Площадь ячейки сетки (рис. 12.3) определяется соотношением Площадь = ] д ] ~ а|Ать (12. 14) которое, согласно (12.13), дает физическую интерпретацию обратного якобиана. Физическая ориентация расчетной сетки (касательная к координатной линии $) относительно оси х задается косинусом угла наклона Гл. !2. Обобитениые криволинейные координаты полностью ортогональных или конформных систем координат возможно лишь при сравнительно простой форме границ расчетной области. При этом возникают некоторые ограничения на положение сеточных узлов. Для двумерной ортогональной сетки должно выполняться условие О = 90' (рис.

12.3). Тогда, согласно (12.17), дм —— хах„+ уаун = О. (! 2.18) В трехмерном случае условие ортогональности имеет вид ум=О, 1чь!. (12.19) Если система координат ортогональная, т. е. метрический тензор содержит только диагональные члены дн, удобно ввести обозначения Иа=(йп), 1=1, 2, 3 (нет суммирования). (12.20) Члены Ь< могут рассматриваться как скалярные множители, поскольку на ортогональных сетках при малых изменениях координаты й' в физической области имеет место соотношение Ьз = ЬтЛ4' (нет суммирования).

(12.2 1) В обобщенных координатах ($,т1) уравнение (12.23) можно записать в следующем виде: где дат и т. д. определяются из (12.12), а дца — из соотношения це ца !д! д =х уч — хчуа. (12.25) В двумерном случае из условия ортогональности, согласно (! 2.! 2), следует х„= — уй Ай, уч = хаАК, (12. 22) где Ац — отношение сторон сетки (12.16). Прн АК=! условия (12.22) сводятся к условиям Коши— Римана и сетка получается конформная. Если АК постоянно, но ие равно единице, простое изменение масштаба по $ и Ч приведет к соответствующей конформной системе координат. Уровень сложности уравнений в различных системах координат может быть оценен из рассмотрения уравнения Лапласа. В двумерном случае в декартовых координатах оно имеет вид д'Т д'Т (12.23) $12.1. Преобразование коордииат 67 Другими словами, у — детерминант матрицы д. Такое обозначение, а не (й(, как в (12.4), использовано здесь для компактности.

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

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

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

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