Главная » Просмотр файлов » Пасконов В.М., Полежаев В.И., Чудов Л.А. Численное моделирование процессов тепло- и массообмена

Пасконов В.М., Полежаев В.И., Чудов Л.А. Численное моделирование процессов тепло- и массообмена (1185910), страница 6

Файл №1185910 Пасконов В.М., Полежаев В.И., Чудов Л.А. Численное моделирование процессов тепло- и массообмена (Пасконов В.М., Полежаев В.И., Чудов Л.А. Численное моделирование процессов тепло- и массообмена.djvu) 6 страницаПасконов В.М., Полежаев В.И., Чудов Л.А. Численное моделирование процессов тепло- и массообмена (1185910) страница 62020-08-25СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

двухточечными схрмами. Пусть х„— один из узлов сетки, р — соответствующее иначе~а> ние приближенного решения. Обозначим через р ~х, х„, у„)' <Л1, точное решение уравнения (1.4 1), принимающее' при х = х значение ра ° Разность р~,х~+ы хв, рп ) ра+г на" (М (йы эывается ошибкой на шаге. Легко проверить,что для схем (1.4.4) — (1.4.8) ошибка на шаге множителем и отличается от ошибки аппроксимации. Для ошибки в целом, т. е. для разности точного и приближенного решений на всем интервале, можно доказать следующее порядковое соотношение: ! у — уь) ю(А)/Ь,.

где со(й) — характерная величина ошибки на шаге. *) Здесь и в дальнейшем звездочками «выделен материал„ который можно пропустить при первом чтении. Глава 2 НАЧАЛЬНЫЕ СВЕДЕНИЯ О.МЕТОДЕ СЕТОК. СЛУЧАЙ НЕСКОЛЬКИХ НЕЗАВИСИМЫХ ПЕРЕМЕННЫХ 5 2Л. Основные понятия ..2Л. ЛЛ. Модельные уравнения и краевые задачи. В этой главе рассматриваются простейшие уравнения математической физики, в частности такие: (2.1Л) ди да — —:-=0 де дхз (2Л.2) да дза —., + —., = О. дх две (2.1.3 Уравнения (2ЛЛ) — (2.1.3) математически моделируют соответственно конвективный одномерный перенос тепла, нестационарное одномерное распространение тепла и стационарное двумерное распределение тепла.

Уравнения (2.1.1~), (2.1.2) являются вволюционнььми, уравнение (2Л.З) — неэволюционное. Для эволюционных уравнений характерно наличие выделенного независимого . переменного, играющего роль времени, а также корректность задач с начальными условиями (существование и единственность решения, непрерывная зависимость от данных задачи). Всюду в.дальнейшем в этой главе, кроме з 2.6, имеются в виду эволюционные уравнения. В математической физике при рассмотрении уравнений с частными производными. обычно требуется определить решение в какой-то области 6 по условиям, заданным на некоторых частях границы этой области (краевая задача). Простым примером является задача Каши для уравнения (2ЛЛ).

"найти решение и(г, 1) . в области ( — <х(+, т>0), удовлетворяющее начальному условию и(0, х) ~р(х), где ~р(х) — заданная функция. Другой пример — первая краевая задача для модельного уравнения теплопроводности (2Л.2). Здесь С есть прямоуголь- 30 Начальное условие и(0, г)) ~р(х), участвующее в .задаче Коши для уравнения (2.1.1), порождает сеточное качильиое условие иг = «р(п»й). (2Л.бг 31 ник (0<х<Х, 0<1< Т); начальное условие эадается.

в виде и(0, х) <р(х), граничные услоеиз — в виде и(г, 0) = У»(г), и(», Х) У»), где ~р(х), /о(э), /,(») — иввестные функции. 2Л.2. Сетка. Сеточные функции. Сеткой наэывается дискретная совокупность точек на плоскости (Г, х) — углое сетки. ое сетки. Важным частным случаем является раеномерная прямоугольная сетка П, х, ), »" г,+пЬг, х х,+ + п»Ьх, где Ь», Ьх — положительные числа, наэываемые. шагами сетки по» и х соответственно; и, и — полые числа (рис. 2Л). Совокупность узлов, соответствующих какому-либо фиксированному значению и, навываетсяслоем.

Для краткости часто обозначают ,и о Ьх=й, Ь»=т. Функция, заданная в узлах сетки,наэывает- 'о о о ся сеточной фуикцйей. г Сеточные функции часто рассматривают как элементы вектор' ного пространства, которое может быть нонечномерным или бесконечномерным. При этом ' р пользуются понятием кормы. Простым примером нормы сеточной функции является верхняя грань ее модуля: (!и1 = эпр!и!. 2.1.3, Идея метода сеток.

Согласно'методу сеток дифференциальное уравнение и краевые условия ваменяются сеточными уравнениями, связывающими аначения искомой функции в узлах сетки (сеточная краееая гадача или схема). Построим сеточное уравнение, соответствующее дифференциальному уравнению (2Л.1).

Воспольэуемся равномерной прямоугольной сеткой 1" пт, и О, 1, 2, ...; х = тй, и = О, ~1, ь2, ... Заменим проиэводную ди/д» в точке (пт, тй) раэностным отношением «впередь ( и+» ч1~ и — и //т, а производную ди/дх в той же точке— раэностным отношением «наэад» (й,„— и~ »)/и. Получим '" " ' =О. (2Л.4) Совокупность узлов сетки, используемых при построении сеточного уравнения, называется шаблоном. Для уравнения (2Л.4) шаблон показан на рис.

2.2. Уравнения (2.1.4), (2.1.5),вместе с описанием сетки составляют схему, приближающую задачу Коши для уравнения (2.1.1). * Приведем еще один, более сложный пример. Для урав- нения — — — — /((,х) = 0 ди ди дхз (2Л.6) рассматривается краевая задача с дополнительными условиями (2Л.7) (2Л.8) 0<х<Х, 0< 2<Т, и(0, х) = ср(х), — "," ~ + 7,(() = О, (2Л.9) (2Л.10) , В области (2Л.7) введем сетку х„тй, т О, 1, 2, ..., М„Ь Х/ЛХ; (о~Л.11) 1"=пт, п=О, 1, 2, ...,,М, т=Т/)т'. (2.1Л2) Уравнение (2Л.6) заменим сеточным на четырехточечном ~п 1,т~ оьт-о ос~~ Рис.

2.2. Ркс. 2.3. шаблоне (пт, (т — 1)Ы, (пт, тЫ, (пт, (т+1)й), ((и+ + 1)т, тЫ (рис. 2.3): ""+' ", + '"-' /"„- = О. (2'.~.13) Начальное условие (2.1.8) ааменяется очевидным образом: 1 и' = (тЬ) (2ЛЛ4) 22 Граничные условия (2.1.9), (2.1.10) интерпретируем' так: — '„«+ ог (пт) = О, (2.1Л5) "„" '+дз( )=О. (2.1.16) В данном случе схему, соответствующую задаче (2Л.6)— (2Л.1О), составляют сеточные соотношения (2.1.18)— (2Л.16). 2Л.4, Ошибка аппроксимации схемы на решении. Близость схемы н исходной краевой аадачи проще всего оценивается по величине нееяеки, получающейся при подстановке точного решения в уравнения и граничные условия сеточной краевой задачи.

Обозначим (2.1.17) всю совокупность уравнений, входящих в краевую аадачу, т. е. основное дифференциальное уравнение и краевые (начальные и граничные) условия. Запишем аналогично сеточную краевую задачу В,(и,) = О. , (2.1Л8) Для сокращения записей мы предполагаем здесь и в дальнейшем, что т т(Ь), т- О при Ь- О, так что сетка определяется одним параметром Ь. Ошибкой аппроксимации схемы (2Л.18) на точном решении задачи (2ЛЛ7) называется сеточная функция се» =" =Л,(и). Схема называется аппроксимируюи(ей на точном решении, если при Ь- О ошибка аппроксимации стремится к нулю: ем- О при Ь- О.

Если а»=О(Ь«), то говорят, что порядок аппроксимации равен р. Рассмотрим схему (2Л.4), (2Л.5). Предположим, что точное решение и = и(г, х) имеет непрерывные равномерно ограниченные производные второго порядка. Пусть, кроме того, т-ЬЬ, Ь сопэФ. Пользуясь формулой Тейлора, получим где» и «» означают, что вторые производные соответствуют некоторым точкам, близким к точке (пт, тЬ). Так как начальное условие переносится в схему без ошибки, то и,=О(й), порядок аппроксимации первый. В.

М. Пас»оиов и др. 33 х Рассмотрим также схему (2Л.13) — (2ЛЛ6). Подстав- ляя точное решение в левую часть (2.1.13) и применяя формулу Тейлора с начальной. точкой 1=пт, х=иЬ, 1 ди 1 зэки получим невязку -т — —,Ь' —.При т=гЬх, г сопэ$, 2 д1з 12 дх4 невязка уравнения (2.1ЛЗ) есть 0(Ь*). Невязка начально- го условия (2Л.14) равна нулю. Невязки граничных ус-. ловий (2.1.15), (2ЛЛ6), как легко видеть, являются величинами первого порядка относительно Ь. Таким об- разом, ошибка аппроксимации для (2ЛЛЗ) — (2Л.16) есть 0(Ь), т. е.

имеем аппроксимацию первого порядка. Заменив приближенные граничные условия (2.1.13), (2.1Л6) более точными (см. (1.3.5)), получим схему, име- ющую второй порядок аппроксимации: з + о, (пт) = О, 2Л.5. Ошибка аппроксимации иа классе функций. По- нятие ошибки аппроксимации вводят и другим способом. Полагают аь = Яь(и) — Л(и)„где и — произвольная до- статочно гладкая функция из некоторого функциональ- ного класса К Легко видеть, что в этом смысле схема (2Л.4), (2Л.5) также имеет первый 'порядок аппроксима- ции.

Покажем на примере того же уравнения (2Л.4), что порядок, аппроксимации на точном решении может быть выше порядка аппроксимации для класса функций, об- ладающих той же гладкостью, что и точное решение. Пусть Ь=1, т. е. т=Ь, Применяя формулу Тейлора, получим и+1 и и и им — ит и,и — и„, + ь /ди ди 1~ Э 1 ди ди1 =~ — + — ) + — 2~ — — — ~ +0(Ь). \д8 дх)щ ' ~ д1з дхз) и Ясли и(1, х) есть точное решение уравнения (2ЛЛ), то, дифференцируя (2.1.1), легко выводим соотношение д~и ди д~ дхз В этом случае и, = В,(и) = 0(Ч). Для произвольной функции и(С, з)) имеем иь = Лл(и) — Л(и) = 0(Ь), так как д'и/дд и д'и/дх' могут принимать различные значения.

й 2.2. Примеры сеточных аппроксимаций ,2.2Л.' Способ конечных разностей. При построении се' точного уравнения мы заменили частные производные ди/дС, ди/дх разностными отношениями. Замена отдельных частных производных, входящих в уравнение, сеточными выражениями, аппроксимирующими зти производные, представляет наиболее простой и часто употребляемый путь построения аппроксимирующих схем — способ конечных разностей.

Приведем наиболее часто применяемые сеточные аппроксимации для производных ди/дС, ди/дх: и(С+т,х) — и(С,х) ди т дзи дС 2 дС« и(С, +Ь) — и(С,. ) ди Ь дти Ь дх 2 дхт — разностные отношения «вперед»', и(С,х) — и(С вЂ” т,х) ди «д и т дС 2 дС« и(С,х) — и(С,х — Ь) ди Ь д»и Ь дх 2 дх« — — + ... (2.2.4) — разностные отношения «назад»; и(С ст,х) — и(С вЂ” т,х) ди «д и + (225) 2т дС 6 д«з и (С,х + Ь) — и (С, х †' Ь) ди Ь« д и 2Ь вЂ” + — — + ° дх '6 — с(ентрально-разностньсе отношения. Для аппроксимации второй производной дй/дх' обычно применяется симметричное разностное отношение второго порядка и (С, х+ Ь) — 2и (С, х) + и (С, х — Ь) д и Ь» дх«+ 12 дх« (2.2.7) Все эти формулы непосредственно получаются с помощью 3* 35 соответствующих тейлоровских разложений; многоточиями обозначены члены высших порядков.. 2.2.2.

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

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

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