Решение краевых задач конечно-разностными методами (Лабораторная работа № 2), страница 4

Тест 3:

x

Y(h=0.1)

Y(h=h/2)

Y(h=h/4)

Y*

0

 1.02557640571455e-003

 2.97503302121971e-004

 7.94865844440280e-005

0

0.1

 1.21148612403335e-001

 1.20327897941378e-001

 1.20087038928366e-001

 1.20000000000000e-001

0.2

 2.81242366935388e-001

 2.80351068213686e-001

 2.80092797301028e-001

 2.80000000000000e-001

0.3

 4.81308426818694e-001

 4.80367512308628e-001

 4.80096898910928e-001

 4.80000000000000e-001

0.4

 7.21348667447361e-001

 7.20377802614601e-001

 7.20099499762230e-001

 7.20000000000000e-001

0.5

 1.00136533317407e+000

 1.00038262094259e+000

 1.00010078539091e+000

 1.00000000000000e+000

0.6

 1.32136105313010e+000

 1.32038277175794e+000

 1.32010097531264e+000

 1.32000000000000e+000

0.7

 1.68133883823151e+000

 1.68037919234678e+000

 1.68010032692246e+000

 1.68000000000000e+000

0.8

 2.08130207329876e+000

 2.08037296572215e+000

 2.08009914057894e+000

 2.08000000000000e+000

0.9

 2.52125451455745e+000

 2.52036534120336e+000

 2.52009776740043e+000

 2.52000000000000e+000

1

 3.00120030488225e+000

 3.00035776906816e+000

 3.00009662178809e+000

3

Выводы:

Чтобы построить консервативную схему для стационарного одномерного уравнения теплопроводности, мы воспользовались методом баланса, т.е. интегрируя на интервалах сетки нашли соответствующие выражения  для нахождения коэффициентов схемы:

, .

Учитывая некоторые упрощения данных формул нахождения коэффициентов: ,а также то, что погрешность аппроксимации первых краевых условий равна 0, можно сделать вывод, что погрешность аппроксимации краевой задачи Дирихле  растет от  погрешности исходной схемы .  Для задачи Ньютона погрешность возникает также при аппроксимации самих краевых условий.

Интерполируя подынтегральные функции, выбирая разностную схему для аппроксимации потока, мы можем получить тот или иной порядок аппроксимации исходной задачи.

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

Что касается практических тестов, то величина порядка аппроксимации с использованием 3-х вложенных сеток варьируется от 1.3 и выше в зависимости от вида приближаемого решения и линейности коэффициента k . На первом и третьем тесте мы получили желаемые результаты, на втором тесте порядок аппрокисмации не подтвердился.

Из результатов численных исследований следует, что построенная схема с выбранными коэффициентами является устойчивой. Погрешность непрерывна на всем интервале, на границах приближается к нулю, следовательно, метод сходиться.