- •Выполнил:
- •Алгоритм решения
- •Реализация алгоритма
- •Z ai |ф*ФА = /ф* cosxdx’ к = 0,4,
- •Разрешающие соотношения
- •Алгоритм решения
- •Результаты расчётов
- •Разрешающие соотношения
- •Алгоритм решения
- •Реализация алгоритма
- •2.1. Аппроксимация решения кусочно-линейными функциями
- •Алгоритм решения
- •Результаты расчетов
- •Разрешающие соотношения
- •Алгоритм решения
- •Реализация алгоритма
- •Разрешающие соотношения
- •Алгоритм решения
- •Реализация алгоритма
- •3.1. Явная разностная схема
- •Алгоритм решения
- •Реализация алгоритма
- •Разрешающие соотношения
- •Программа 3.2
- •Реализация алгоритма
- •Алгоритм решения
- •Реализация алгоритма
- 89,1268a0+178,254a, - 89,1268a2+0,0^ + 0,0a4= 52,7393,
«0,0a0-89,1268a, +178,254a289,1268a3 + 0,0a4= 74,5847,
0,0a0+ 0,0a, - 89,1268a2 +178,254a389,1268a4 = 52,7393,
0,0a0+0,0a, +0,0a2-89,1268a3 + l 19,127a4= 609,968
относительно коэффициентов Г, разложения (2.9) решения диффе ренциального уравнения.
Таблица 2.1
Матрица коэффициентов и правая часть системы линейных алгебраических уравнений для определения коэффициентов разложения решения дифференциального уравнения в ряд по кусочно-линейным пробными функциями
1,0 |
0,0 |
0,0 |
0,0 |
0,0 |
100,0 |
-89,1268 |
178,254 |
-89,1268 |
0,0 |
0,0 |
527,393 |
0,0 |
-89,1268 |
178,254 |
-89,1268 |
0,0 |
745,847 |
0,0 |
0,0 |
-89,1268 |
178,254 |
-89,1268 |
527,393 |
0,0 |
0,0 |
0,0 |
-89,1268 |
119,127 |
699,684 |
Результаты расчетов
Использование Программы 2.1 даёт решение этой системы уравнений:
Т0= 100,0, Т}=103,407, Г2=100,897, Тъ=90,018, Г4-73,2219.
Приближённые решения дифференциального уравнения на от резке [0, л], полученные на основе кусочно-линейного представле ния (2.9), приведены на рис. 2.1 с использованием 4 (а) и 64 (б) сегментов.
Полученное решение даёт возможность определить потоки тепла с торцевых поверхностей стержня. Первое и последнее уравнения системы (2.16)
|
|
SHUT, —SHIXi |
|
|
|
COSJtft - - |
^ - Q 0, |
- - T , + |
- г - - и |
' sm(n-h) |
+ cos(rc-/?)-£?*• |
h m_l |
h m |
c |
|
<2 ‘ = - к { s"l(” ~ 'l) + cos(« - л) ] + |
- г . ) . |
Рис. 2.1. Приближённое решение уравнения АТж(лг)+ W{\) =О
с граничными условиями Г(х)| г=0= Т°
и А,Г'(*)| = - |
а |
-7^] при использовании |
4 (а) и 64 |
(б) сегментов на отрезке [0,7г] |
Точное решение поставленной задачи
T{x) = W0synxlX +\_W0- a { T ° - T j) \x /{ \ + an) +T0 |
(2.18) |
позволяет определить погрешность полученного приближённого ре шения и точные значения потоков на торцах стержня:
|
Х(ж0-а (г ° -7 :))1 |
x(w0-a{T°-T„)) |
||
е ° = ^ И = о = |
W0 cosx:-+ —---- |
г------------- |
L |
= ^ 0+ —---- г------------- |
|
А+ ап |
|
А+ ап |
|
|
Х(г0-а(г°-Г _)) |
X(w0 -а(т° - TJ) |
||
|
W0COSJC+ |
Х+ал |
|
=W0- |
|
|
|
\+ ап |
В табл. 2.2 приведены величины тепловых потоков для прибли жённых решений, полученных при различном числе сегментов. Точ ные значения потоков: Q0=403,34 Вт/м2и QL=1596,66 Вт/м2
Таблица 2.2
Приближённые значения тепловых потоков на торцах стержня при различном числе m сегментов
m |
Q° |
& |
m |
2 ° |
of- |
2 |
363,380 |
596,660 |
i 6 |
403,341 |
1577,45 |
4 |
403,340 |
1303,77 |
32 |
403,339 |
1591,85 |
8 |
403,339 |
1520,54 |
64 |
403,340 |
1595,47 |
На рис. 2.2, а показана зависимость от координаты х погрешности
Д(*) = Г(х) - Тл(х) = Т(х) - X |
(*) • |
|
/=0 |
|
|
В табл. 2.3 приведены значения погрешностей |
|
|
К - F - 7 - . | = a g f M - r . M - g a |
nt)-t,T < p,U ) |
|
|
2m |
m |
8«,2ш= \K ~ Tm\ =max|r2m(x)- Tm(x)( = max £ 7>, ( x ) - ^ 7 ’(pi(x),
Рис. 2.2. Погрешность приближённого решения дифференциального уравнения, полученного на четырёх сегментах с использованием пробных кусочно-линейных функций (а); зависимость погрешностей 8т(—о—)и him (-Д-) приближённого решения от длин сегментов (б)
определённых с помощью точного решения 7\х) дифференциаль ного уравнения и сравнением двух последовательных приближён ных решений Тти Tim дифференциального уравнения при различ ных значениях числа т и 2т сегментов с использованием чебышёвской нормы. Эти же данные представлены в графическом виде на рис. 2.2, б.
Таблица 2.3
Погрешность Ътприближённого решения дифференциального уравнения при различном числе т сегментов
т |
h |
5m |
^т.2т |
т |
h |
5m |
&т.2т |
2 |
1,570796 |
3,00733 |
2,95867 |
16 |
0,196350 |
0,068358 |
0,0657 |
4 |
0,785398 |
1,00522 |
1,00466 |
32 |
0,098175 |
0,016445 |
0,0165 |
8 |
0,392699 |
0,268929 |
0,2585 |
64 |
0,049087 |
0,004208 |
0,0041 |
С применением формулы (В.1) приближенно определяется поря док погрешности численного решения дифференциального уравне ния разложением (2.9) по системе кусочно-линейных пробных функ ций. Для 8тпорядок погрешности оценивается значением
Ьт= (In3,00733 - In0,004208)/(ln 1,570796 - In0,049087) = 1,89622.
Для5т,2m эта величина
Погрешность полученного решения дифференциального урав нения методом Галёркина с использованием системы кусочно линейных пробных функций в обоих случаях можно приближённо оценить как величину второго порядка относительно длин h сегмен тов, т.е.
8т = \т -? т Л = о {и 2), |
ътЛт=|2т;.ср( - £ |
г,ф, |» о (а2). |
|
II 1=0 II |
II 1=0 |
/=О |
II |
В силу этого &т->0 при //->0 или |
Это позволяет утвер |
ждать, что последовательность приближённых решений дифферен циального уравнения, полученных аппроксимацией кусочно-посто янными функциями (2.9), сходится равномерно на отрезке [0, тс].
Выполненные расчёты (см. табл. 2.3, рис. 2.2, б) показывают, что погрешности Ъти §2тпрактически одинаково оценивают погрешность приближённого уравнения. Это означает, что при отсутствии точного решения, когда погрешность 8топределить невозможно, для контроля погрешности можно применять оценку ?>2т. На рис. 2.3 приведена за висимость времени t выполнения расчетов от числа т сегментов раз ностной сетки.
Рис. 2.3. Зависимость времени t выполнения расчетов от числа т сегментов разностной сетки
Выводы
1.Процедура метода Галёркина использована для приближённого решения дифференциального уравнения. Сформирована система ли нейных алгебраических уравнений для определения коэффициентов разложения искомого решения по системе кусочно-линейных функций.
2.Разработана вычислительная программа определения коэффи циентов разложения решения дифференциального решения по сис теме кусочно-линейных функций.
3.С использованием разработанной программы определены ко эффициенты и построены приближённые решения дифференциаль ного уравнения для 2, 4, 8, 16, 32 и 64 сегментов постоянной длины (см. рис. 2.2).
4.Для указанной последовательности разложений определены оценки Ът и bim погрешности приближённого решения (см. табл. 2.3) для различного числа т сегментов.
5.Показано, что с уменьшением длины А сегментов погрешность приближённого решения дифференциального уравнения, определяе мая чебышёвской нормой, уменьшается (см. рис. 2.2, б). Установле но, что погрешность аппроксимации имеет первый порядок относи тельно длин сегментов (шага интегрирования).
6.Выполненное исследование показывает, что последователь ность приближённых решений дифференциального уравнения, полу ченных методом Галёркина при аппроксимации кусочно-линейными
функциями, сходится равномерно на отрезке [0, л].
7. Для численного решения дифференциального уравнения ме тодом Галёркина с использованием кусочно-линейных функций на разностной сетке, содержащей 64 сегмента, на компьютере с процес сором Intel® Pentium® 4 (тактовая частота 2,2 ГГц, объем оператив ной памяти 512 Мбайт) требуется 2,5*10-3 с.
2.2. Аппроксимация решения кусочно-квадратичными функциями
Задание. Методом Галёркина с использованием кусочно-квадра тичных пробных функций построить приближённое решение одно мерного дифференциального уравнения стационарной теплопровод
ности [A.r'(jc)]' + W0sin х =0 |
с граничными условиями |
= Т°, |
\Т /(х]х=к= - a [r(jc)|j=7i -7^]. |
Сформировать систему линейных алгеб- |