close

Вход

Забыли?

вход по аккаунту

?

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

код для вставкиСкачать
ВЕСТН. САМАР. ГОС. ТЕХН. УН-ТА. СЕР. ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ. 1998. № 6
Краткие сообщения
Дифференциальные уравнения
УДК 517.956.3
А.Ю. Сеницкий
ИНТЕГРАЛЬНЫЕ ПРЕДСТАВЛЕНИЯ РЕШЕНИЙ ГИПЕРБОЛИЧЕСКОГО
УРАВНЕНИЯ С ПЕРЕМЕННЫМИ КОЭФФИЦИЕНТАМИ
Рассмотрены некоторые случаи интегрируемости гиперболического уравнения общего вида. Показано, что наряду с приемами преобразования уравнений, использующими групповой анализ, методы
факторизации, предложенная процедура построения общего решения является весьма эффективной.
В настоящей работе предложена процедура построения общего решения гиперболического уравнения общего вида, которая является весьма эффективной наряду с приемами преобразования уравнений,
использующих групповой анализ [1, 2] и факторизационные методы [3, 4].
В области W: {t>0, 1 < r < a} рассмотрим следующее дифференциальное уравнение:
¶ 2U (r, t ) ¶ 2U (r, t )
¶U
+ A (r )
+ B (r ) × U = 0 ,
2
2
¶r
¶r
¶t
(1)
где A(r), B(r) Î C [1, a].
Будем считать, что функция U(r,t) удовлетворяет условиям Дирихле, следовательно, ее можно представить интегралом Фурье, записываемым в виде пары формул:
1
~
U ( r, p) =
2p
U (r, t ) =
Считая, что U ( r,0) =
1
2p
+¥
ò U ( r, t ) e
- ipt
dt ;
(2)
0
+¥
~
ò U (r, p)e
ipt
dp .
(3)
0
¶U ( r,0)
= 0 , применим преобразование (2) к (1). Тогда в пространстве изо¶t
бражений получим следующее уравнение:
~
~
¶ 2U ( r, p)
¶U ( r, p)
~
+ A (r )
+ B (r ) + p2 U (r, p) = 0 .
2
¶r
¶r
[
]
(4)
Представим решение уравнения (4) в такой форме:
~
U ( r, p) = j ( r ) × G( s) ….,
s = pY( r ) ,
где j(r), Y(r) и G(s) - дважды непрерывно дифференцируемые функции своих аргументов.
Подставляя (5), (6) в (4), приходим к дифференциальному соотношению
(5)
(6)
é d 2G
ù
é Y ¢¢ × Y 2j ¢Y A Y ù dG
1
Y2
S ê 2 +
Gú + S ê
+
+
+
[j ¢¢ + A j ¢ + B j ]G = 0 , (7)
ú
( Y ¢) 2 û ë ( Y ¢) 2 jY ¢ Y ¢ û dS j × ( Y ¢) 2
ë dS
2
которое может быть удовлетворено различными способами.
Случай 1. Пусть
d 2G
1
G = 0;
2 +
dS
(Y ¢)2
j ( r) × Y ¢¢( r) + 2j ¢( r ) × Y ¢( r) + A ( r )j ( r)Y ¢( r) = 0 ;
120
(8)
(9)
Y ¢¢( r ) + A (r )j ¢(r ) + B (r )j (r ) = 0 .
(10)
Полагая, что Y¢(r)=1, т.е.
Y(r)=r,
(11)
общее решение уравнения (8) записывается в виде
G(s) = C 1 ( p) e is + C 2 ( p)e - is .
(12)
После подстановки равенства (12) в (5) и его обращения по формуле (3), получаем такое представление для оригинала [5]:
U (r, t ) =
j (r ) +¥
2p
ò [C ( p)e
ipY
1
]
+ C 2 ( p)e - ipY e ipt dp = j (r )[ f1 (t + Y ) + f 2 (t - Y)] .
(13)
0
С учетом (11) равенство (9) существенно упрощается. Имеем
1
j ¢( r ) + A ( r) × j ( r ) = 0 ,
2
откуда находим
j ( r) = e
-
1
A ( r ) dr
2
ò
.
В результате подстановки (14) в (10) получаем следующую функциональную зависимость:
B ( r) =
1
2A ¢( r ) + A 2 ( r) .
4
[
]
(14)
(15)
Таким образом, приходим к такому утверждению.
Т е о р е м а 1. Для того чтобы выражение (13) являлось замкнутым решением уравнения (1), достаточно выполнения условий (11), (15).
Случай 2. Пусть в рамках представлений (8) - (10) вместо (11) справедливо равенство
1
[ Y ¢( r)]
2
=
1
Y ( r)
2
.
(16)
Тогда (8) является уравнением Эйлера и его общее решение записывается в элементарных функциях.
Имеем
G( s) = C 1 ( p)M 1 ( Y) + C 2 ( p) M 2 ( Y ) .
(17)
Здесь
1
M 1,2 (Y) = [ Y( r)]
1 æ 1 2ö
±ç - p ÷
ø
2 è4
M 1,2 (Y) = [ Y(r )]
1 æ 2 1ö
± iç p - ÷
2 è
4ø
1
2
M 1 ( Y) = [ Y(r )] ,
,
1
,
p>
1
;
2
1
;
2
p=
M 2 ( Y) = [ Y(r )] ln Y(r ),
U ( r, t ) =
1
;
2
2
1
2
Следовательно,
p<
2
p=
1
.
2
j ( r ) +¥
2p
ò [C ( p)M (Y) + C ( p)M (Y)]e
1
1
2
2
(18)
ipt
dp .
(19)
0
Из (16) следует, что Y( r) = e . Тогда справедливы такие соотношения:
r
æ 1 ö
æ 1
ö
j (r ) = expç - r ÷ × expç - ò A(r )dr ÷ ,
è 2
ø
è 2 ø
1
B(r ) = 2 A¢(r ) + A 2 (r ) - 1
4
[
]
(20)
и имеет место следующее предложение.
Т е о р е м а 2. Для того чтобы решить (19), достаточно выполнения второго из равенств (20).
Рассмотрим теперь альтернативный (8) - (10) вариант построения общего решения уравнения (1).
Случай 3. Пусть имеет место (11), а также следующие равенства:
121
Y ¢¢( r) × j (r )
[ Y ¢(r)]
2
Y 2 ( r)
j ( r)[ Y ¢( r)]
2
+2
j ¢( r )Y( r) A( r)Y(r )
+
= 0;
Y ¢( r )
j ( r )Y ¢ ( r )
(21)
[ A(r)j ¢(r) + B (r)j (r) + j ¢¢(r)] = -6 ;
(22)
В этом случае равенство (8) принимает вид
d 2 G( s)
S
+ ( s2 - 6)G( s) = 0 .
2
ds
2
Общее решение последнего уравнения можно представить следующим образом [6]:
ì3
ü
3ö
æ
G( s) = C 1 ( p) í cos s + C 2 ( p) + ç1 - 2 ÷ sin s + C 2 ( p) ý .
è s ø
îs
þ
[
]
[
]
(23)
Выполняя обращение равенств (5), (23), находим
j (r ) +¥
ì3
ü
3ö
æ
U (r, t ) =
C 1 ( p)í cos s + C 2 ( p) + ç1 - 2 ÷ sin s + C 2 ( p) ýe ipt dp .
ò
è s ø
2p 0
îs
þ
[
]
[
]
(24)
Из соотношения (21) определяем
é 1
ù
j ( r) = exp ê- ò A ( r) drú .
ë 2
û
Наконец, из уравнения (22) получаем
B ( r) =
1é
24 ù
2A ¢( r) + A 2 ( r ) - 2 ú .
ê
4ë
r û
(25)
Т е о р е м а 3. Если коэффициенты уравнения (1) удовлетворяют соотношению (25), то выражение
(24) является его общим решением.
Случай 4. Пусть выполняется условие (11) и справедливы следующие соотношения:
Y ¢¢(r ) × Y(r )
[ Y ¢ ( r )]
2
Y 2 (r )
[ Y ¢ ( r )] j ( r )
2
+ A( r)
j ¢ ( r ) Y( r )
Y( r )
+2
= 1;
j (r )Y ¢( r )
Y ¢( r )
(26)
[j ¢¢(r ) + A (r)j ¢(r) + B (r)j (r)] = -n
2
,
(27)
где nÎR.
Тогда (7) трансформируется в уравнение Бесселя, т.е.
s2
d 2 G( s)
dG( s)
+s
+ ( s 2 - n 2 )G( s) = 0 .
2
ds
ds
Его решение записывается в виде
G( s) = C 1 ( p)I n ( s) + C 2 ( p)Y n ( s) .
(28)
Здесь In(s), Yn(s) - функции Бесселя "n"-го порядка I-го и II-го рода.
Выполняя обращение выражения (5), с учетом (28) определяем U(r,t). Имеем
U (r, t ) =
j (r ) +¥
2p
ò [C ( p)I (s) + C ( p)Y (s)]e
1
n
2
n
ipt
dp .
(29)
0
Принимая во внимание (11), из (26) находим функцию
j ( r) = r
1
2
é 1
ù
exp ê- ò A( r) dr ú .
ë 2
û
(30)
Далее, из (27) следует, что
B ( r) =
1é
1 - 4n 2 ù
2
(
)
(
)
2
A
¢
r
+
A
r
+
ê
ú.
4ë
r2 û
(31)
Т е о р е м а 4. При выполнении условия (31) выражение (29) является замкнутым решением
уравнения (1).
122
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Зайцев В.Ф., Полянин А.Д. Справочник по нелинейным дифференциальным уравнениям. М.: Наука, 1993. 462 с.
2. Зайцев В.Ф. Дискретно-групповой анализ обыкновенных дифференциальных уравнений // Дифференциальные
уравнения. 1989. Т.25. №3. С. 379-387.
3. Беркович Л.М. Факторизация и преобразования обыкновенных дифференциальных уравнений / Под ред.
Н.Х. Розова. Саратов: Изд-во Сарат. ун-та, 1989. 192 с.
4. Сеницкий Ю.Э., Сеницкий А.Ю. О решении одного дифференциального уравнения с переменными коэффициентами, имеющего приложение в динамической теории упругости // Мат. физика и нелинейн. механика. 1990.
Вып.13. С. 22-25.
5. Бейтмен Г., Эрдейн А. Таблицы интегральных преобразований. Т.1. М.: Наука, 1969. 343 с.
6. Камке Э. Справочник по обыкновенным дифференциальным уравнениям. 3-е изд. М.: Наука, 1965. 703 с.
УДК 519.63
Т.А. Бенгина, М.Ю. Лившиц
КОНЕЧНО- РАЗНОСТНАЯ СХЕМА РЕШЕНИЯ ЗАДАЧИ ТИПА
СТЕФАНА С ПРЕОБРАЗОВАНИЕМ КООРДИНАТ
Предлагается численный метод решения задачи типа Стефана, основанный на преобразовании
координат, трансформирующий изменяющуюся во времени область в стационарную, и на конечно-разностной схеме решения полученной после этого задачи.
Задачи теплопроводности и диффузии с наличием фазовых переходов классифицируют как задачи с
подвижной границей или задачи типа Стефана. Аналитические методы применяются в основном для исследования однофазных или двухфазных линейных моделей различных физических процессов. Получить
точное решение нелинейной модели с помощью существующих аналитических методов практически невозможно. В связи с этим при моделировании все чаще используются численные методы решения подобных задач, однако и с их помощью не всегда удается достаточно точно отследить закон движения границы раздела фаз. Подобная проблема возникла при изучении математической модели процесса диффузионного насыщения металла. Многофазная диффузия описывается параболическим уравнением переноса
для каждой из фаз, которые в совокупности с краевыми условиями и условиями сохранения баланса на
границах фазовых переходов представляют собой задачу Стефана.
Математическая модель газового азотирования имеет вид
¶C ( x , t )
¶ æ
¶C ö
ç Di (C )
÷,
=
¶t
¶x i è
¶x i ø
(1)
где t > 0; 0 < x < L ; C ( x , t ) — концентрация азота в точке x в момент времени t ; Di (C ) – коэффициент диффузии азота в i-той фазе.
Краевые условия
- Di (C )
¶C ( x , t )
¶x
(
= b i p0N - p iN
x =L
)
,
(2)
¶C ( x , t )
=0,
¶x
x =0
(3)
где b i — коэффициент массопереноса i-той фазы; p N — азотный потенциал печной атмосферы; p N —
равновесный азотный потенциал соответствующей фазы.
Начальные условия задачи:
0
C ( x , t) t=0 = C 0 ( x ) = 0 ,
где x i ( t ) — граница i-той фазы.
Условия на границах раздела фаз:
i
(4)
x i (0) = x i 0 ,
123
Документ
Категория
Без категории
Просмотров
4
Размер файла
186 Кб
Теги
решение, уравнения, интегральная, коэффициента, представление, переменных, гиперболическое
1/--страниц
Пожаловаться на содержимое документа