2. Трехточечная схема с весом.
Все эти схемы сводятся к стандартному виду (3.4) и решаются методом прогонки
(3.4)
Коэффициенты Ai, Bi, Ci должны удовлетворять условиям:
(3.5)
Коэффициенты B0 , C0 , F0, AN ,CN ,FN находятся из граничных условий. В данной задаче в зависимости от знака функции p(x,t) ставятся граничные условия и тем самым находятся наши коэффициенты.
1) Когда р>0 задается правое граничное условие:
(3.3?)
Используя уравнения (3.3?) находим коэффициенты AN ,CN ,FN . Коэффициенты B0 , C0 , F0 находятся из дополнительного условия, которое ставится на левом конце.
2) Когда р<0 задается граничное условие на левом конце
(3.3?)
Используя уравнения (3.3?) находим коэффициенты B0 , C0 , F0
Коэффициенты AN ,CN ,FN находятся из дополнительного условия, которое ставится на правом конце.
3.3.1 Центрально-разностная схема
Разностная схема задачи (3.1)-(3.3) имеет следующий вид:
1) р>0. В этом случае граничное условие задается на правом конце:
(3.6)
Используя уравнение (3.6) находим коэффициенты AN =0, CN=1,
Дополнительное условие на левом конце имеет вид:
(3.7)
Приведем уравнение (3.7) к виду :
(3.7?)
Отсюда находим коэффициенты:
2) В случае, когда р<0, граничное условие ставится на левом конце
(3.8)
Используя уравнение (3.8) находим коэффициенты B0,=0, C0=1,
Дополнительное условие на правом конце имеет вид:
(3.9)
Приводим уравнение (3.9) к виду :
(3.9?)
отсюда находим коэффициенты:
Таблица 13. Численное решение уравнения переноса с постоянными коэффициентами центральная разностная схема метод прогонки
-------------kogda p>0--------------50sloy |
||||
N priblijennoe tochnoe pogreshnosti |
||||
0 |
0.03544452 |
0.03678794 |
0.00134342 |
|
1 |
0.03541069 |
0.03558189 |
0.00017120 |
|
2 |
0.03306824 |
0.03441538 |
0.00134714 |
|
3 |
0.03313883 |
0.03328711 |
0.00014828 |
|
4 |
0.03084494 |
0.03219583 |
0.00135089 |
|
5 |
0.03101552 |
0.03114032 |
0.00012480 |
|
6 |
0.02876471 |
0.03011942 |
0.00135472 |
|
7 |
0.02903119 |
0.02913199 |
0.00010080 |
|
8 |
0.02681828 |
0.02817693 |
0.00135865 |
|
9 |
0.02717688 |
0.02725318 |
0.00007630 |
|
10 |
0.02499699 |
0.02635971 |
0.00136272 |
|
11 |
0.02544422 |
0.02549554 |
0.00005132 |
|
12 |
0.02329272 |
0.02465970 |
0.00136698 |
|
13 |
0.02382538 |
0.02385126 |
0.00002588 |
|
14 |
0.02169787 |
0.02306932 |
0.00137145 |
|
15 |
0.02231302 |
0.02231302 |
0.00000000 |
Таблица 14. Численное решение уравнения переноса с постоянными коэффициентами центральная разностная схема метод прогонки
-------------kogda p<0--------------50sloy |
||||
N priblijennoe tochnoe pogreshnosti |
||||
0 |
0.03678794 |
0.03678794 |
0.00000000 |
|
1 |
0.03475182 |
0.03558189 |
0.00083008 |
|
2 |
0.03440516 |
0.03441538 |
0.00001021 |
|
3 |
0.03246493 |
0.03328711 |
0.00082218 |
|
4 |
0.03217504 |
0.03219583 |
0.00002079 |
|
5 |
0.03032529 |
0.03114032 |
0.00081503 |
|
6 |
0.03008771 |
0.03011942 |
0.00003171 |
|
7 |
0.02832337 |
0.02913199 |
0.00080861 |
|
8 |
0.02813396 |
0.02817693 |
0.00004297 |
|
9 |
0.02645027 |
0.02725318 |
0.00080290 |
|
10 |
0.02630518 |
0.02635971 |
0.00005453 |
|
11 |
0.02469766 |
0.02549554 |
0.00079788 |
|
12 |
0.02459330 |
0.02465970 |
0.00006639 |
|
13 |
0.02305773 |
0.02385126 |
0.00079352 |
|
14 |
0.02299077 |
0.02306932 |
0.00007855 |
|
15 |
0.02152320 |
0.02231302 |
0.00078982 |
Текст программы смотри в приложении 5
- Введение
- 1.1 Сеточная область
- 1.2 Сеточная функция. Пространство сеточных функций. Нормы сеточных функций
- 1.3 Аппроксимация дифференциальных операторов
- 1.4 Разностная схема
- 1.5 Корректность разностной схемы
- 1.6 Аппроксимация и сходимость
- 1.7 Неравномерная сетка
- 1.7.1 Построение сеточной области
- 1.7.2 Формирование сетки
- Глава II. Одномерное уравнение переноса с переменными коэффициентами
- 2.1 Постановка задачи
- 2.2 “Явные ” схемы
- 2.3 Неявные схемы
- 2.3.1 Центрально разностная схема
- 2.3.2 Трехточечная схема с весом
- 2. Трехточечная схема с весом.
- 3.3.2 Трехточечная схема с весом
- 3.3.3 Схема “прямоугольник”
- 3.3.4 Схема со сглаживанием
- 3.3.5 Схема прямоугольник со сглаживанием
- 3.3.6 “Шахматная” схема
- Заключение
- Основные положения метода сеток для решения задачи Коши
- 7.2. Порядок аппроксимации разностной схемы. Оценка порядка аппроксимации разностной схемы с весами для нестационарного уравнения теплопроводности.
- Разностные схемы для линейного одномерного уравнения переноса
- 4. Особенности построения конечно-разностных аналогов уравнений переноса
- 6.2. Разностные уравнения, явная и неявная схемы
- Построение разностных схем для уравнений в частных производных
- Устойчивость, консервативность разностных схем. Разностные сетки и преобразование основных уравнений
- Однородные разностные схемы для уравнения теплопроводности