logo
Разностные схемы для уравнения переноса на неравномерных сетках

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