Заглавная страница Избранные статьи Случайная статья Познавательные статьи Новые добавления Обратная связь FAQ Написать работу КАТЕГОРИИ: АрхеологияБиология Генетика География Информатика История Логика Маркетинг Математика Менеджмент Механика Педагогика Религия Социология Технологии Физика Философия Финансы Химия Экология ТОП 10 на сайте Приготовление дезинфицирующих растворов различной концентрацииТехника нижней прямой подачи мяча. Франко-прусская война (причины и последствия) Организация работы процедурного кабинета Смысловое и механическое запоминание, их место и роль в усвоении знаний Коммуникативные барьеры и пути их преодоления Обработка изделий медицинского назначения многократного применения Образцы текста публицистического стиля Четыре типа изменения баланса Задачи с ответами для Всероссийской олимпиады по праву Мы поможем в написании ваших работ! ЗНАЕТЕ ЛИ ВЫ?
Влияние общества на человека
Приготовление дезинфицирующих растворов различной концентрации Практические работы по географии для 6 класса Организация работы процедурного кабинета Изменения в неживой природе осенью Уборка процедурного кабинета Сольфеджио. Все правила по сольфеджио Балочные системы. Определение реакций опор и моментов защемления |
Систему линейных алгебраических уравнений перепишем в виде↑ ⇐ ПредыдущаяСтр 7 из 7 Содержание книги
Поиск на нашем сайте
(7.11) i = 1, 2, …, N-1 при каждом j=0,1,…,M-1. где
Метод прогонки для решения разностной схемы. Система (7.9) – (7.10) является системой линейных алгебраических уравнений с N-1 неизвестными. Полученная задача решается методом прогонки. Пусть . (7.12) После подстановки в (7.11) получится рекуррентная формула
(7.13) Для определения преобразуя (7.10) приводим его к виду (7.14) где Сравнивая (7.14) с (7.12) получаем, что
(7.15)
Теперь из (7.12), (7.15) определяются все
После этого рассматривая совместно (7.12) и (7.9) вычисляются все
В данном случае все условия теоремы 2 из §4 выполняются, поэтому метод прогонки для решения задачи (7.9) – (7.11) является устойчивой.
Расчетная схема
1) Используя заданные функций вычисляются
(7.16)
2) Из рекуррентного соотношения определяются все 3) После этого используя формулу правой прогонки определяются все .
Переменные и блок – схема
В данном случае искомая функция зависит от двух переменных t и х. Поэтому соответствующая сеточная функция зависит от двух дискретных переменных i и j. При программировании мы должны резервировать место в оперативной памяти компьютера для двухмерного массива. БЛОК-СХЕМА
описание термодинамических характеристик грунта. - - - - - Вычисление параметров разностной схемы и начальной функций
J = 0, M - 1, 1
коэффициента прогонки - - - - - - I = N-1, 1, -1
I = 0, N - 2, 1
рис. 3.
Если M и N достаточно большие величины, то в оперативной памяти компьютера может не хватит места для массива . Чтобы избежать этого, вводятся одномерные массивы . Вместо массивов Ai, B i, C i используются идентификаторы A, B, C. Для отводятся одномерные массивы . В формуле (7.16) при определений Ai,Ci,B i используется отношение . Если это выражение очень большое, то вычислительный процесс будет не устойчивой. В этом случае не выполняется теорема 2. Теорема 3. Пусть решение дифференциальной задачи (7.4)-(7.7) θ обладает непрерывными производными до четвертого порядкаи если , то решение разностной схемы (7.8) – (7.10) сходится к решению дифференциальной задачи (7.4) – (7.7). Алгоритм реализации схемы расчета приведен на рис. 3. Задания для лабораторной работы. Составить программу и произвести расчет при следующих начальных данных распределения температуры и теплопроводных характеристиках грунта. 1) Температура воздуха изменяется по закону . Варианты значений θmax и θmin:
Глубина заложения трубопровода Н = 1,5 м; tmax =24 час, считаем, что начальное распределение температуры от трубопровода до поверхности изменяется по закону Температура на поверхности трубопровода равен θ1 = 450.
Теплопроводные характеристики грунта приведены в таблице 1.
Обратная задача для уравнения теплопроводности Постановка задачи
Пусть - коэффициент теплопроводности, - температура в точке z в момент времени t. Тогда процесс распространения тепла на отрезке описывается уравнением
В прямой задаче надо найти по известной функции . Для единственности решения прямой задачи необходимо задать начальные условия
и по одному условию на каждой из границ, например, поток
Если данные прямой задачи (1)-(4) достаточно гладкие и , то решение прямой задачи существует и единственно.
Обратная задача Пусть о решении прямой задачи (1)-(4) известна дополнительная информация Требуется определить один из коэффициентов уравнения (1) (или какую-либо их комбинацию) из соотношений (1)-(5). Численное решение обратной задачи (1)-(5) будем искать при , минимизируя целевой функционал Зададим начальное приближение . Приближение будем вычислять методом простой итерации
Здесь - достаточно малое число, - градиент функционала .
Обратная задача 1: найти коэффициент из соотношений:
Найдем приращение функционала (6):
Здесь . является решением следующей задачи
Рассмотрим сопряженную задачу:
Умножим обе части равенства (12) на функцию и проинтегрируем по области :
Проинтегрируем по частям выражение
Имеем
Тогда учитывая условия (13)-(15) и (17)-(19), получим, что
Последние два слагаемых в правой части равенства (21) имеют второй порядок малости. Тогда получаем следующий градиент функционала Здесь - решение сопряженной задачи.
|
|||||||||||||||||||||||||||||||||||||||||
Последнее изменение этой страницы: 2016-04-23; просмотров: 227; Нарушение авторского права страницы; Мы поможем в написании вашей работы! infopedia.su Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав. Обратная связь - 18.118.10.75 (0.006 с.) |