Определение собственных чисел краевой задачи Штурма-Лиувилля

Тип работы:
Реферат
Предмет:
Физико-математические науки


Узнать стоимость

Детальная информация о работе

Выдержка из работы

УДК 536. 2
В. А. Кудинов, В. В. Дикоп, Р. Ж. Габдушев, С. А. Назаренко
ОПРЕДЕЛЕНИЕ СОБСТВЕННЫХ ЧИСЕЛ КРАЕВОЙ ЗАДАЧИ ШТУРМА — ЛИУВИЛЛЯ
Приводится методика определения собственных чисел краевой задачи Штурма — Лиувилля, основанная на совместном использовании точных (метод разделения переменных) и приближённых аналитических (взвешенных невязок) методов. Показана высокая точность получаемых решений и хорошая сходимость к точному решению при увеличении числа приближений.
В методе разделения переменных искомое решение принимается в таком виде, что оно с самого начала точно удовлетворяет дифференциальному уравнению краевой задачи. Собственные числа и неизвестные коэффициенты решения находятся из точного выполнения граничных и начальных условий. Однако не для всех дифференциальных уравнений можно найти функции, точно удовлетворяющие им. Такие функции известны лишь для простейших линейных дифференциальных уравнений применительно к телам классической формы (пластина, цилиндр, шар) при классических линейных граничных условиях. При усложнении любого из этих факторов трудоемкость применения метода Фурье существенно возрастают вплоть до практической невозможности его использования.
В настоящей работе рассматривается метод, основанный на совместном использовании метода Фурье и методов взвешенных невязок (ортогонального метода Бубнова — Галеркина). Ввиду того, что удовлетворение уравнению в данном случае осуществляется путём составления его невязки и требования ортогональности невязки к собственным функциям, дифференциальное уравнение может быть сколь угодно сложным. Не накладывается также ограничений на вид начальных и граничных условий.
В качестве конкретного примера рассмотрим задачу теплопроводности для бесконечно протяжённой пластины при граничных условиях первого рода:
Э0 Э 20
— = - (№ & gt- 0- 0 & lt- л & lt- 1) — (1)
ЭЕо Эл
0(л, 0) = 1- (2)
Э (c)(0, Ео)/Эл = 0- (3)
0(1, Ео) = 0, (4)
где 0(л, Ео) = (Т — Тст)/(Г0 — Тст) — относительная избыточная температура- Тст — температура
пластины при л = 1- Л = х/8 — безразмерная координата- 8 — толщина пластины- Т0 — начальная
температура- Ео = ат/8 — число Фурье- а — коэффициент температуропроводности- т — время.
Следуя методу Фурье, решение задачи (1) — (4) принимается в виде
0(л, Ео) = р (Ео)^(л). (5)
Подставляя (5) в (1), получим
ёр (Ео)/аЕо + т 2р (Ео) = 0- (6)
ё ?(л)/ёл2 + Л2 ?(л) = 0, (7)
2
где т — некоторая постоянная.
Решение уравнения (6), как известно, имеет вид
р (Ео) = А ехр (- /л2 Ео), (8)
где, А — неизвестный коэффициент.
Уравнение Штурма — Лиувилля (7) представим следующим образом:
?11 (л) + 1?(л) = 0, (9)
где 1 = /л2.
Граничные условия для уравнения (9) согласно (3), (4) будут иметь вид
?1 (0) = 0- (10)
?(1) = 0. (11)
Решение задачи в первом приближении разыскивается в виде следующего ряда:
Y (h) = c0 + Chh + Ch2 + Ch + Ch4 +… = XChг, (12)
i=0
где Ci, (i = 0, n) неизвестные постоянные, определяемые из граничных условий (10), (11). В настоящей работе решается задача определения пяти собственных чисел. Так как неизвестных пять, а граничных условий два, то следует добавить дополнительные граничные условия, которые найдём из уравнения (9) путём выполнения самого этого уравнения, а также производных от него в точках h = 0 и h = 1. Эти дополнительные граничные условия будут иметь вид
Y (0) = const = 1- (13)
Y11 (1) = 0- (14)
Y111 (0) = 0. (15)
Граничное условие (13) следует из условия (10).
Подставляя (12) в (10), (11), (13) — (15), получим
C0 = 1- Cj = 0- C2 =-1,2- C3 = 0- C4 = 0,2.
Соотношение (12) принимает вид
Y (h) = 1 — 1,2h2 + 0,2h4. (16)
Для определения первого собственного числа найдём интеграл взвешенной невязки урав-
нения (9)
J[Y11 (h)+AY (h)]h = 0. (17)
0
Подставляя (16) в (17), получим алгебраическое линейное уравнение относительно Я, решение которого дает Я = 2,6.
Точное значение первого собственного числа Я = 2,46 740 110 027 [1].
Собственная функция находится из (12).
Для получения первых двух собственных функций используются следующие основные и дополнительные граничные условия:
Y (0) = 1- Y1 (0) = 0- Y11 (0) = -Я- Y111 (0) = 0- Y1V (0) = Я2- Y (1) = 0- Y11 (1) = 0. (18)
Подставляя (12) в (18) при n = 8, получим
C0 = 1- C1 = 0- C2 =-0,5Я — C3 = 0- C4 = Я2/24-
39 59
C5 =-Я2 + 1,4Я — 3- С6 = 2 — 0,9Я------Я2.
5 8 6 12
После подстановки коэффициентов С{ в соотношение (12) составляется невязка уравнения (9) и требуется ортогональность невязки к функции (12). Отсюда для определения собственных чисел получается полином шестой степени. Первые два корня этого полинома имеют вид
Я = 2,46 740 110 — Я = 22,26 983.
Ввиду того, что уравнение Штурма — Лиувилля удовлетворяется лишь при некоторых дискретных значениях Я (собственных значениях), то остальные корни полинома отбрасываются, как не удовлетворяющие этому уравнению (в этом можно убедиться непосредственной подстановкой в уравнение 9).
Таким образом, во втором приближении первое собственное число до восьмого знака после запятой совпало с точным его значением. Второе собственное число совпадает с его точным значением до второго знака после запятой (точное значение Я2 = 22,20 660 990).
Для получения трёх собственных чисел в дополнение к условиям (18) воспользуемся граничными условиями вида
YV (0) = -Я3- Ym (0) = 0 — Ym1 (0) = Я4- Y1X (0) = 0 — Yх (0) = -Я5. (19)
При выполнении граничных условий (18), (19) для первых трёх собственных функций получаются следующие значения:
Я = 2,4 674 011 001- Я = 22,2 066 135- Я = 62,55 342.
В данном случае первое и второе собственные числа совпадают с точными их значениями соответственно до 10 — го и 4 — го знаков после запятой. Точное значение третьего собственного числа Я3 = 61,68 502 750.
Для получения пяти собственных чисел добавляем следующие граничные условия:
Vх1 (0) = 0- Vх11 (0) = Я6- Vх111 (0)= 0- УХ1Г (0) = -Я7-
V Х7 (0)= 0- V17/1 (0) = Я8- V17/11 (0) = 0- V ^ (0) = -Я9-
V Х1Х (0) = 0- V хх (0) = Я10.
В этом случае имеем
Я1 = 2,4 674 011 002- Я2 = 22,206 610- Я3 = 61,6 850 235- Я4 = 120,90 249- Я5 = 201,0584. Точные значения четвёртого и пятого собственных чисел
Я4 = 120,90 265- Я5 = 199,8595.
Подставляя (8), (12) в (5), для каждого собственного числа будем иметь частные решения
вида
0, (п рв) = А1(l^, 1) ехр (- Хрв), (I = 0,4).
Каждое частное решение точно удовлетворяет граничным условиям (18), (19), (20) и приближённо (в пятом приближении) удовлетворяет уравнению (1) на отрезке 0? п? 1, включая
точки п = 0 и п = 1, ввиду выполнения условий ?77 (0)=-Х и ?77 (1) = 0 (см. условия (18)). Од-
нако ни одно из этих частных решений, в том числе и их сумма
0(, Ее) =? А, V (, Я,)ехр (- Х, Еа),
(21)
,=0
не удовлетворяют начальному условию.
Для выполнения начального условия составим его невязку и потребуем ортогональность невязки к каждой собственной функции, т. е.
? А, V, (, Я,)-1
=0
(22)
(1-е) 1,0
Определяя интегралы в (22), для нахождения коэффициентов
АК (К = 1,5) получим систему из пяти алгебраических линейных уравнений. Её решение дает А1 = 1,274 366 — А2 = -0,427 128-
А3 = 0,257 304 — А4 = -0. 182 685- А5 = 0,150 535.
Приведённые здесь коэффициенты Аг (, = 0,4) найдены из уточнённых значений собственных чисел.
Результаты расчётов безразмерных температур по формуле (21) в сравнении с точными их значениями [3] представлены графически на рисунке.
Анализ результатов расчётов позволяет заключить, что в пятом приближении значения безразмерных температур, полученные по формуле (21), удовлетворительно согласуются с точными их значениями в диапазоне Фурье 0,005? Рв? ?.
Ф і го г избыточной
0,001 Ц005 0,01 0,05 0,1
Г рафики изменения относительной
температуры от числа Бо:
— - точное решение [3]- о — расчёт по формуле (21)
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Кудинов В. А., Карташов Э. М. и др. Тепломассоперенос и термоупругость в многослойных конструкциях. М.: Энергоатомиздат, 1997. 426 с.
2. Кудинов В. А., Аверин Б. В., Габдушев Р. Ж., Стефанюк С. А. Об одном методе определения собственных значений краевой задачи Штурма — Лиувилля // Вест. СамГТУ. Сер. Физ. -мат. науки. 2001. Вып. 12. С. 51−56.
3. Лыков А. В. Теория теплопроводности. М.: Высш. школа, 1967. 600 с.

ПоказатьСвернуть
Заполнить форму текущей работой