Решение модифицированным методом Фурье




 

Для решения ДУЧП теплопроводности: в данной краевой задаче используем модификацию метода Фурье на случай полубесконечной прямоугольной области, на двух противоположных сторонах которой нет нулевых граничных условий для функции U (x, t).

Представим искомую функцию U (x, t) в виде суммы двух функций:

Для одной из этих функций - V (x, t) поставим аналогичную краевую задачу с нулевыми условиями на двух противоположных сторонах прямоугольной области D, которую можно решить классическим методом Фурье, вторую функцию - W (x, t) найдем подбором.

V (x, t):
V|t=0 = 0 - W|t=0
V|x=l = 0
V|x=0 = 0
 
l
t
x
 
D  
W(x, t):
W|x=l = 0
W|x=0 = At
 
l
t
x
 
D  

 

I. Найдем функцию W (x, t):

Находим подбором функцию W (x, t), удовлетворяющую условиям:

Этому ДУЧП можно попробовать удовлетворить, представив функцию W (x, t) в виде многочлена второй степени от аргументов x и t:

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

Следовательно, нужно ввести еще коэффициенты, то есть взять многочлен третьей степени по x и t:

Сначала удовлетворяем ДУЧП:

Подставляем в ДУЧП и требуем, чтобы оно тождественно выполнялось по всем x и t:

коэффициенты при

Удовлетворяем остальным условиям:

коэффициенты при

коэффициенты при

Выбрав функцию W (x, t) в виде многочлена третьей степени с неопределенными коэффициентами получили ровно 10 различных условий на эти коэффициенты. Из этих условий определим все коэффициенты и запишем окончательный вид W (x, t).

Проверка:

 

ДУЧП: , Ɐx, Ɐt

Следовательно, функция W (x, t) подобрана верно.

 

II. Найдем функцию V (x, t):

 

Классическим методом Фурье находим функцию V (x, t), как единственное решение краевой задачи, у которой скорректированы нулевые граничные условия.

 

1. Искомую функцию V (x, t) ищем в виде произведения двух функций, каждая из которых зависит только от одного аргумента:

Эта функция позволяет в исходном ДУ реализовать процедуру разделения переменных:

ДУЧП:

2. Перебрасываем нулевые условия при x=0 и x=l на функцию V (x, t):

3. Решаем задачу Штурма Лиувилля для второго уравнения системы (7) вместе с условиями (8)

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

ДУ из системы (9): – ДУ второго порядка относительно функции X (x) линейное, однородное с постоянными коэффициентами, поэтому

где – константы, – ФСЧР, которая находится с помощью корней характеристического уравнения.

Характеристическое уравнение имеет вид:

1 случай: λ > 0 – действительные различные корни k2=-k1

Таким образом, , то есть при λ > 0 задача Штурма-Лиувилля не имеет решения.

2 случай: λ = 0 – действительные одинаковые корни

Таким образом, , то есть при λ = 0 задача Штурма-Лиувилля не имеет решения.

3 случай: λ < 0 – комплексные сопряженные корни

Проверка собственных значений функций на удовлетворение условиям (9):

- верно.

Таким образом, в результате решения задачи Штурма-Лиувилля получено:

4. Решаем второе уравнение системы системы (7):

Возможен только один случай – когда λ<0: при λ<0 имеем линейное однородное ДУ с постоянными коэффициентами:

– ДУ первого порядка относительно функции T(t), относится к типу ДУ с разделяющимися переменными.

5. В соответствии с формулой (6) перемножаем функции Xn(x) и Tn(t) при Ɐn=1,2,3.. и получаем счетное множество Vn (x, t):

Каждая из функций (12) удовлетворяет ДУЧП и нулевым граничным условиям. По исходной задаче остается неудовлетворенным начальное условие при t=0. Для того, чтобы ему удовлетворить, функции (12) имеют произвольные константы Cn. Чтобы все эти константы использовать составим функциональный ряд из функций Vn (x, t):

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

Сделав эти общие предположения относительно составленного ряда (13), обозначим его сумму через V (x, t) и получим ее получим ее выражение в следующем виде:

В силу сделанных предположений относительно ряда, имеем, что составленная функция V (x, t) удовлетворяет ДУЧП и нулевым граничным условиям, так как этому удовлетворяет каждый член ряда.

6. Удовлетворение начальному условию:

Обозначим – известная функция, заданная для

- это равенство представляет собой тригонометрический ряд Фурье по синусам, составленный для функции f (x), где , дополненный нечетным образом на промежуток (-l; 0) и продолженный на всю числовую ось с периодом T=2l.

Следовательно, коэффициенты этого разложения вычисляются по известным формулам Фурье:

Вычисляем Сn:

Запишем функцию V (x, t), используя формулу (14) и значения коэффициентов Сn:

III. Найдем функцию U (x, t):

Складывая найденные функции V (x, t) и W (x, t), запишем функцию U (x, t):

Таким образом, используя модификацию метода Фурье, получили решение краевой задачи Дирихле для одномерного уравнения теплопроводности в полубесконечной прямоугольной области. Это решение имеет вид (16).

 

Управляемые параметры задачи:

A – некоторая постоянная;

l – длина стержня;

a2 – постоянная, зависящая от теплофизических характеристик материала.



Поделиться:




Поиск по сайту

©2015-2024 poisk-ru.ru
Все права принадлежать их авторам. Данный сайт не претендует на авторства, а предоставляет бесплатное использование.
Дата создания страницы: 2019-12-18 Нарушение авторских прав и Нарушение персональных данных


Поиск по сайту: