Главная » Просмотр файлов » Fedorenko-RP-Vvedenie-v-vychislitelnuyu-fiziku

Fedorenko-RP-Vvedenie-v-vychislitelnuyu-fiziku (810773), страница 108

Файл №810773 Fedorenko-RP-Vvedenie-v-vychislitelnuyu-fiziku (Fedorenko-RP-Vvedenie-v-vychislitelnuyu-fiziku) 108 страницаFedorenko-RP-Vvedenie-v-vychislitelnuyu-fiziku (810773) страница 1082020-08-18СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

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

62). Теперь уже будет во- 510 ПРИБЛИЖЕННЫЕ МЕТОДЫ ВЫЧИСЛИТЕЛЬНОЙ ФИЗИКИ 1ч. н семь базисных функций, которые определяются специальным выбором значений р на внешней границе ячейки, а внутрь продолжаются, так же как н раньше, решением уравнения Лапласа с учетом условий на скважине (если она есть). Краевые условия, например для р, и <рз, формулируются так (Е = х/Н, т1 = у/Н): р,(х, 0) = 1 — ЗЦ+ 2ЦВ, р,(0, у) = 1 — 311+ 211«, р, ее 0 на остальных гранях, р (х, 0) =4Ц1 — Ц), р =0 на остальных гранях. На рис. б2 показаны граничные значения Р1 и рь. С каждой счетной точкой на границе ячейки нужно связать свою часть ее контура ог Теперь, кроме «целых» точек (х, лг), появляются «полуцелые» точки (/г, лг+1/2), (/г+!/2, лг) и соответствующие двусторонние контуры о, о«+и о«+из (рис.

63). Опыт применения схемы вто- «1.Н1, рого порядка показал заметное повышение точ- «1,„-«л ности расчета. Но это связано и с заметным возрастанием сложности схемы: в целой точке она 21-точечная, в полуцелой — 13-точечная. При Рис. 63 выборе элементарных контуров о и двусторонних контуров о „,..., как показал опыт, следует придерживаться естественного правила; контуры о« м, о«»нз а«пз не должны перекрываться и должны покрывать все координатные линии сетки (не допуская «пустот»).

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

Вернемся к уравнению (2) для давления, причем и(х, у) считается известной функцией. Реально эта функция известна в узлах сетки, т.е. в виде сеточной функции ие , полученной как «грубое решение» (см. з 30) при расчете раскрытия трещины. Кроме того, известно, что около внешней границы и представляется е виде нД т1) Я~) т~~ + О(т1зи) (12) где $ — длина дуги на дс., т! — расстояние от дс' по нормали. При столь сильном вырождении уравнения (2) «краееое условие» на дб ставится в очень неопределенной форме — е виде требования ограниченности решения.

Этим «условием» надо замкнуть систему разностных уравнений. Ниже в общих чертах мы опишем численную реализацию такого замыкания. Начнем с некоторых терминов. Пусть контур дсг задан, как и в з 30, набором точек (Х!, У!), соединяемых отрезками прямых. Плоскость (х, у) покрыта сеткой узлов (й, и!). Точ- с а-узлы ° д-увен ки, попавшие внутрь б, называются счетными.. В них определена сеточная функция р„,„(н Рис.

64 и, используемая для вычисления коэффициента диффузии, который полагался постоянным в пределах ячейки; эта постоянная величина вычислялась усреднением и тю четырем вершинам ячейки). Разделим счетные узлы на два типа (рис, 64). Узлы в-типа («внутренние») — это узлы (х, !и), у которых все восемь соседних узлов явлшотся счетными. В этих узлах используется е стандартная разностная аппроксимация на 9-точечном шаблоне. Остальные счетные узлы суть узлы д-типа («дефектные»). В этих узлах неприменимо стандартное разностиое уравнение и нет привычною краевого условия, с помощью которого е нееырожденных задачах записываются разностные уравнения в узлах Рис. 65 д-типа.

Для того чтобы разобраться в том, как следует поступать в таком случае, рассмотрим решение в окрестности границм, вводя местную систему координат (рис. 65). В этих переменных уравнение можно записать в форме (13) Здесь мы используем только главные члены асимптотики (12), прн этом Д можно считать постоянной величиной. Преобразуем пгиьлижвнныв мноды вычислительной»и»яки 1ч. и 512 выражение (13) к виду; г г А7 Ч ! — + - А7 Ч вЂ” + А7 Ч вЂ” = а+ ... з ггдд 3 з !7гдл з з7гал дчг 2 дя Считая величину Чиг дгР7д$г пренебрежимо малой по сравнению с остальными, рассмотрим асимптотнческое уравнение Ч + Ч =!7= —.

ига'Р . З Пгдг д дяг 2 дя Л!3' Оно легко интегрируется. Его общее решение есть Р(Ч) = С, + 2С»Ч !72+ 2!7Ч!7г где С, С, — произвольные постоянные. Очевидно, ограниченное решение можно получить только при С = О. (Именно зто сокращение числа произвольных постоянных в общем решении уравнения второго порядка и является причиной того, что нельзя ставить классическое краевое условие при Ч = О, т.е. на границе области Сг.) Итак, мы получили асимптотику решения около границы: Р(б, Ч) = Сй) + 2(7й) Ч!и+ о(Ч!и). (14) Используя асимптотику (14) запишем соотношения для р„, Р;! 7г! - — С+ 2!7!IЧд, Р, = С + 2!7г7Ч!. Исключая из них С, получаем связь Ра !«Ра + (1 !«) Рд + ! (~7Чд «Чд).

20 (15) Подставляя в (!3) такое решение (предположим, что СД), !7(г)— гладкие функции), можно оценить отброшенные выше члены и оправдать зги действия. В формуле (14) функция С($), конечно, неизвестна, она определяется некоторой процедурой «сшивки» с решением внутри области, Алгоритмическая реализация асимптотики состоит в следующем. Рассмотрим некоторую точку д-типа (точка !7 на рис. 65). Пусть (Х1, У!) — ближайшая к ней вершина границы.

Через вершины с номерами ! — 1, 1, 7+ 1 проведем окружность, которую будем считать участком границы д0. Через точку д-типа проведем нормаль к дб, пересекающую вертикальные и горизонтальные границы ячеек или их диагонали. Среди зтих отрезков выберем ближайпгий к точке д-типа, такой, что оба его конца суть точки д-типа (на рис. б5 они обозначены а и р). Пересечение нормали с отрезком обозначим буквой 5. Значение р в точке з можно проинтерполировать по значениям в точках а и р: р,= мр + (1 — !г) рд, !г != 10, 1). 5!3 8ЗН метод кОнечных супеРэлементов Формулу (15) в итерационном процессе решения системы разностных уравнений можно использовать двумя способами. Первый способ состоит в следующем, Пусть получено некоторое приближение р во внутренних узлах в-типа. Используя (15), доопределим его в точках д-типа. Таким образом, значения Р известны во всех счетных точках. Теперь по стандартной формуле простейшего метода релаксации новые значения ре можно получить-во всех внутренних точках.

Далее процесс повторяется. Уточним некоторые детали. Представим разностную схему в форме ! 1 (16) ! -!! — ! Метод релаксации с ускорением состоит в следующем. В каждой внутренней точке (к, л!) поочередно пересчитывается значение ! = — (с + «~~' Ок ' Р .)!СО О (17) !,/ — ! (звездочка отмечает пропуск слагаемого !'= !'=О). Это предварительное значение. Окончательное значение есть Ра,м' РКН+!О(Рк, Рдт)' (18) где параметр !О, ускоряющий сходимость, подбирается экспериментально, Теория уравнения Лапласа с постоянными коэффициентами указывает диапазон оптимального значения: О! 1.7 —: 1.8.

Найденное по (18) значение р» сразу записывается в массив р, так что формула (18) является «полунеявнойж часть входящих в сумму значений ре ! Ну относится к ч-й итерации, часть — к (я+1)-й. Вычислительный эксперимент показал слабую расходимость этого итерационного процесса. Второй способ состоит в том, что связь (15) заранее вносится в разностные уравнения. Если точка д-типа !1 входит в шаблон в точке (к, л!), причем соответствующие точки а и р тоже входят в этот шаблон, то точка д-типа исключаегся из схемы. Выражение (15) подставляется в (1б), пересчитываются коэффициенты схемы, соответствующие точкам а и р, коэффициент, соответствующий точке !1, делается нулевым.

Далее итерационный процесс выполняется точно так же, как было описано выше, но значения в точках д-типа в расчете фактически не участвуют. После выполнения итерации формула (15) используется для расчета значений Р в точках д-типа. Редко, но все же встречаются ситуации, когда хотя бы одна из точек а и р не входит в шаблон в точке (й, и!).

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

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

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

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