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

 


Информация о работе «Разностные схемы для уравнения переноса на неравномерных сетках»
Раздел: Информатика, программирование
Количество знаков с пробелами: 74851
Количество таблиц: 18
Количество изображений: 9

Похожие работы

Скачать
59485
4
20

... на первой  и последующих  итерациях равна: ; (3.22) . (3.23) Критерием завершения итерационного процесса является условие: ,(3.24) где  - заданная точность расчета [4]. 4. Методы оценки термонапряженного состояния 4.1 Физические основы возникновения термических напряжений При изменении температуры происходит объемное расширение или сжатие твердого тела. Неравномерный нагрев ...

Скачать
36871
3
34

... диаметрах критического сечения представлены на рисунке 2.24 Рисунок 2.24 - Зависимость оптимальной высоты поднятия фурмы от давления при различных диаметрах критического сечения сопла Лаваля 3. Численное исследование движения жидкости Приведены уравнения Навье - Стокса установившегося осесимметричного движения несжимаемой вязкой жидкости в переменных функция тока - вихрь. Проведено ...

Скачать
11306
2
0

... системы на ЭВМ, а так же требование его экономичности обуславливают применение регулярных сеток, расположение узлов в которых подчиняется определённым закономерностям. В практике численного моделирования микроэлектронных структур примеяются как непрерывные прямоугольные (неравномерные), так и треугольные сетки (рис.2.). Треугольная сетка позволяет с меньшим количеством дополнительных узлов сгущать ...

Скачать
243425
1
0

... . Реакции узлов более высокого уровня менее зависят от позиции и более устойчивы к искажениям. Структура Неокогнитрон имеет иерархическую структуру, ориен­тированную на моделирование зрительной системы челове­ка. Он состоит из последовательности обрабатывающих слоев, организованных в иерархическую структуру (рис. 10.8). Входной образ подается на первый слой и передается через плоскости, ...

0 комментариев


Наверх