Краткие сообщения Дифференциальные уравнения
УДК 517.956.3 А.Ю. Сеницкий
ИНТЕГРАЛЬНЫЕ ПРЕДСТАВЛЕНИЯ РЕШЕНИЙ ГИПЕРБОЛИЧЕСКОГО УРАВНЕНИЯ С ПЕРЕМЕННЫМИ КОЭФФИЦИЕНТАМИ
Рассмотрены некоторые случаи интегрируемости гиперболического уравнения общего вида. Показано, что наряду с приемами преобразования уравнений, использующими групповой анализ, методы факторизации, предложенная процедура построения общего решения является весьма эффективной.
В настоящей работе предложена процедура построения общего решения гиперболического уравнения общего вида, которая является весьма эффективной наряду с приемами преобразования уравнений, использующих групповой анализ [1, 2] и факторизационные методы [3, 4].
В области О: (1>0, 1 < г < а} рассмотрим следующее дифференциальное уравнение:
д2и(;, г) д2и(;, г) , ,ди
где А(г), В(г) е С [1, а].
Будем считать, что функция и(гД) удовлетворяет условиям Дирихле, следовательно, ее можно представить интегралом Фурье, записываемым в виде пары формул:
1 +¥
и (г р)=^2р 1и (гґ)
- ірї
;
о
+¥
1 -І-1-*-)
и (г, ї) = -/= 1 и (г, р)вгрёр. у2п о
(2)
(3)
( ) ди (г,0)
Считая, что и(г,0) =------- -----= 0, применим преобразование (2) к (1). Тогда в пространстве изо-
дї
бражений получим следующее уравнение:
д2и (г, р)
+ а (г)дидг; р) + [в(г)+р2 ]г~(;, р) =0.
Представим решение уравнения (4) в такой форме: и(;,р) = Р(;) я) ....,
^ = р^(;),
где ф(г), Т(г) и в^) - дважды непрерывно дифференцируемые функции своих аргументов. Подставляя (5), (6) в (4), приходим к дифференциальному соотношению
(4)
(5)
(6)
& 20
1
■ +
2 (^'):
О
¥"•¥ 2р,х¥ А ¥
(¥')2 + ~р¥7 + ^¥7
&О
+
р • (^')
—[р" + А р' + Вр]о = 0, (7)
которое может быть удовлетворено различными способами. Случай 1. Пусть
& О
1
+
О = 0;
2 (^')2
р( г) • ¥"( г) + 2р'(г) • ¥'( г) + А (г)р( г)¥'( г) = 0;
(8)
(9)
Т"(г) + А (г)р'(г) + В (г)р(г) = 0. Полагая, что Т'(г)=1, т.е.
Т(г)=г,
(10)
(11)
общее решение уравнения (8) записывается в виде
а(*) = с (р)в* + с 2 ( рУ* . (12)
После подстановки равенства (12) в (5) и его обращения по формуле (3), получаем такое представле ние для оригинала [5]:
и (;, г) = +¥ [с, (">"* + С, (р)~р ]"> = <?(;)[/. (г + Т) + /, (г - ')].
Ы2п о
С учетом (11) равенство (9) существенно упрощается. Имеем
ф'(;)+2 а(;)- ф(;) = 0,
(13)
откуда находим
р(г) = е ^ . (14)
В результате подстановки (14) в (10) получаем следующую функциональную зависимость:
В (г) = 4 [2А '(г) + А2 (г)]. (15)
Таким образом, приходим к такому утверждению.
Т е о р е м а 1. Для того чтобы выражение (13) являлось замкнутым решением уравнения (1), достаточно выполнения условий (11), (15).
Случай 2. Пусть в рамках представлений (8) - (10) вместо (11) справедливо равенство
1 1
(16)
Тогда (8) является уравнением Эйлера и его общее решение записывается в элементарных функциях. Имеем
Здесь
[Т'(г)]2 Т2(г) ' го общее решение за
О( *) = С, (р)м, (Т) + С 2 (р)М 2 (Т). М 12(Т) = [Т(г)]
м и(Т) = [Т(г)]
(17)
|2± г1р -4
Р >1 2
1 і М1(Т) = [Т(г)]2, |Р = 2;
1 1 М, (Т) = [Т(г)]Чп Т(г), Р = -
Следовательно,
и(г, і) = рР 1 [С 1 (р)М 1 (Т) + С 2 (р)М 2 (Т)]еір'«Р.
у2р 0
Из (16) следует, что Т(г) = ег. Тогда справедливы такие соотношения:
р(г) = ехр1-1 г I • етрі -11 А(г
(18)
(19)
(гіаг
В(; ) = 4 [2 А'(;)+ А2 (;)-1]
и имеет место следующее предложение.
Т е о р е м а 2. Для того чтобы решить (19), достаточно выполнения второго из равенств (20). Рассмотрим теперь альтернативный (8) - (10) вариант построения общего решения уравнения (1). Случай 3. Пусть имеет место (11), а также следующие равенства:
(20)
+¥
121
Y"(r )• j(r) + 2 j'(r )Y(r) + A(r)Y(r) — 0. [Y'(r)]2 j(r )Y'(r) Y'(r) ’
( Y (/())]2 [A(r)jtr) + B(r)j(r) + j'tr)]— -б; j^Y^r)]
(2l)
(22)
В этом случае равенство (8) принимает вид
S
d 2G(s)
ds2
+ (s2 - 6)g(s) — 0.
Общее решение последнего уравнения можно представить следующим образом [6]:
G(s) = С1 (р)I3cos[ + С2(р)] + ^ sin[ + С2. Выполняя обращение равенств (5), (23), находим
U(r, t) = j= JCi ^.р)-^^3cos[ + С2 (р)] + ^ - ^г) sin[s + С2 (р)]|^Ф. Из соотношения (21) определяем
j(r) = exp - 2 { A (r)dr
Наконец, из уравнения (22) получаем
(23)
(24)
24
2A'(/) + A2 (•)- —
(25)
Т е о р е м а 3. Если коэффициенты уравнения (1) удовлетворяют соотношению (25), то выражение (24) является его общим решением.
Случай 4. Пусть выполняется условие (11) и справедливы следующие соотношения:
'"(;).'(;) аП_Т(;). 2 Р'(;)'(;)
Т'(;) р(;)Т'(;)
[Y'(r)] Y2 (r)
— l;
(2б)
[Y'(r)] j(r) где veR.
Тогда (7) трансформируется в уравнение Бесселя, т.е.
2 d2G(s) dG(s)
[j"(r) + A (r)j'(r) + B (r)j(r)] — -v
(27)
ds2
+ s-
ds
+ (s2 - v2 )G(s) — 0.
Его решение записывается в виде
G(s) = Cl (p)!v (s) + C2 (pYv (s) .
Здесь 1П^), Уп^) - функции Бесселя "п"-го порядка 1-го и 11-го рода.
Выполняя обращение выражения (5), с учетом (28) определяем и(гД). Имеем
Р(; )
U^ t) — I [Сl (p)!v (s) + C2 (ppYv (s)]pdp
Принимая во внимание (ll), из (2б) находим функцию
j(r) — r2 exp
Далее, из (27) следует, что
2A '(r) + A2 (r) +
l - 4v
(28)
(29)
(30)
(3l)
Т е о р е м а 4. При выполнении условия (31) выражение (29) является замкнутым решением уравнения (1).
l22
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
Т.А. Бенгина, М.Ю. Лившиц
КОНЕЧНО- РАЗНОСТНАЯ СХЕМА РЕШЕНИЯ ЗАДАЧИ ТИПА СТЕФАНА С ПРЕОБРАЗОВАНИЕМ КООРДИНАТ
Предлагается численный метод решения задачи типа Стефана, основанный на преобразовании координат, трансформирующий изменяющуюся во времени область в стационарную, и на конечно-разностной схеме решения полученной после этого задачи.
Задачи теплопроводности и диффузии с наличием фазовых переходов классифицируют как задачи с подвижной границей или задачи типа Стефана. Аналитические методы применяются в основном для исследования однофазных или двухфазных линейных моделей различных физических процессов. Получить точное решение нелинейной модели с помощью существующих аналитических методов практически невозможно. В связи с этим при моделировании все чаще используются численные методы решения подобных задач, однако и с их помощью не всегда удается достаточно точно отследить закон движения границы раздела фаз. Подобная проблема возникла при изучении математической модели процесса диффузионного насыщения металла. Многофазная диффузия описывается параболическим уравнением переноса для каждой из фаз, которые в совокупности с краевыми условиями и условиями сохранения баланса на границах фазовых переходов представляют собой задачу Стефана.
Математическая модель газового азотирования имеет вид
дС (х, т) д
дт
дх,
* с €
(1)
где Т> 0; 0 < X < Ь ;С (X, т) — концентрация азота в точке х в момент времени X ; О, (С ) - коэффициент диффузии азота в /'-той фазе.
Краевые условия
дС (х, т)
(С)-
дх
х =Ь
дС (х, т)
дх
= 0,
(2)
(3)
х =0 0
где Ь , — коэффициент массопереноса /-той фазы; Ры — азотный потенциал печной атмосферы; ргм — равновесный азотный потенциал соответствующей фазы.
Начальные условия задачи:
С (х • т)|т.0 = С0(х ) = (4)
X, (0) = Х/о,
где X, (т) — граница /-той фазы.
Условия на границах раздела фаз:
123