Файл: Арсенин, В. Я. Методы математической физики и специальные функции учеб. пособие.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 16.10.2024

Просмотров: 175

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Мы докажем соответствующие теоремы для одномерного случая. Для двумерного и трехмерного случаев они форму­ лируются и доказываются совершенно аналогично.

Т е о р е м а 1. Если в задачах Коши

 

а2ихх = ut,

и {х,

0) = фх (л-),

— с х > < л :< с о ,

(59)

U

a2vxx = v(,

v(x,

0) = ф2 (лг),

— о о < л ;< с о ,

(60)

для начальных значений фДх) и ф2 (лг) выполняется нера-

венство | ф1(*) -Ф а(*)1< е (61)

при всех значениях х, то для решений этих задач и (дг, t)

иv (х, t) выполняется неравенство

\и{х, t) v (х, i) | sg е

при всех значениях х и t ^

0 .

Д о к а з а т е л ь с т в о .

Используя формулу (38) для

решения задач (59) и (60), а также неравенство (61), получим

ОО

| V (х, t) -

и (х, t) I

5 G(* ~ О I Ч>1 (I) - ф2 (S) I di ==£

 

 

00

со

 

 

 

 

 

=s£e

$ G{x — l,t)d,l = &,

 

 

СО

 

СО

 

 

так как

lj G(x g,

t) dg = 1 .

 

СО

Справедливость последнего равенства устанавливается непосредственным вычислением интеграла. Производя в нем

замену переменной интегрирования

по формуле £ = л' +

4 - а

аЧ ,

получим dl — Y^cPt da

и

 

СО

 

0 0

 

 

$

G ( x - l , f)di = y ^

$

e-«2da = l.

 

— 0 0

— 0 0

 

Теорема доказана.

Для доказательства устойчивости решения задачи Коши к малым изменениям неоднородности в уравнении, очевидно, достаточно рассмотреть решение с нулевыми начальными значениями.

Те о р е м а 2. Каковы бы ни были положительные числа е

иТ, существует такое 6 = 6 (е, Т), что если в задачах Коши

а2ихх + h (х, 0 = Щ, и {х, 0) = О

(62)

153


и

 

 

 

 

a2vxx + /2 (х,

i) = vt,

v {х,

0) = О

(63)

функции (х, t) и /2 (х, t)

при всех значениях х и 0 sg; t С Т

удовлетворяют неравенству

 

 

 

\h(x,

0 1 < б ( е ,

Т),

(64)

то для решений этих задач и (х, t) и v (х, /) выполняется неравенство

| и(х, t) и (х, t) I

8

при всех значениях х и O ^ t ^ T .

 

Д о к а з а т е л ь с т в о . Используя

формулу (54) для

решения задач (62) и (63), а также неравенство (64), для

O ^ t ^ T

получим

 

 

 

 

 

I и (х, t) V (х, t) |

 

 

 

 

 

t

со

 

 

 

 

 

«£ $

^ G ( x - l , t -

т) j /х (£, т) -

/2 (g,

т) \dldx<:

0 — оо

со

 

 

 

 

 

t

 

 

 

i

 

 

5 G ( x - l , t - T ) d l d x = 8\ 1 d x < 6 - 7 ’.

 

0 — со

 

 

 

0

Положив 6 = г/Т,

для

всех

значений х

и для O ^ t ^ T

будем иметь

 

 

 

е.

 

 

| и (х, t) — v (х, t) |

 

Теорема доказана.

 

 

равенством

 

Мы здесь воспользовались

 

СО

^G ( x - l , t - x ) d l = \ ,

ОО

справедливым при любых значениях т < t . Оно устана­ вливается непосредственным вычислением путем замены

переменной интегрирования g = х + а ]/4 а 2 (^ — т).

З а м е ч а н и е . В трехмерном случае для доказательства соответствующих теорем надо использовать равенство

СО

У ~ т1’ z ~ & t)dldx\dl = \,

СО

154


справедливость которого устанавливается непосредствен­ ным вычислением с использованием замены переменных

£ = * + а ] / 4 аЧ, П = г/+ Р V^cM, £= 2 + у YbF t.

Все остальные выкладки и оценки остаются прежними.

ЗАДАЧИ

1. Начальный ток и начальное напряжение в полубесконечном однородном проводе C X xsC co равны нулю. Самоиндукция единицы длины провода пренебрежимо мала. С момента ( = 0 к концу провода приложена постоянная э. д. с. Е0. Найти напряжение в проводе.

2.Найти распределение температуры в неограниченном простран­ стве, вызванное тем, что в начальный момент i 0 на сферической поверхности радиуса г0 выделилось мгновенно Q равномерно распре­ деленных единиц тепла. (Построение функции влияния мгновенного сферического источника тепла.)

3.Найти концентрацию диффундирующего вещества в неогра­

ниченном

пространстве,

начальная

концентрация

которого

равна

и \t_„ = u0 = const

для

0 ^ r < R

и и;е_0=

0 для

r > R .

предполагая,

4.

Решить

задачу 3

для

полупространства

z >

0,

что Zj <

R, (0, 0,

z0)—координаты центра сферы,

в которой

началь­

ная концентрация равна и0.

Рассмотреть случаи, когда:

а)

плоскость

z = 0 непроницаема для диффундирующего вещества;

б)

на плоскости

г= 0 поддерживается концентрация, равная нулю.

 

 

ОйС д: < со

5.

Найти

температуру

полубесконечного

стержня

с теплоизолированными

концом

и боковой

поверхностью,

обуслов­

ленную действием

тепловых

источников плотности

Q(t) на отрезке

(а, Ь)

(0 < а < Ь )

начиная с момента t = t0.

 

 

 

 

 

 

решить

6.

С помощью функции источника,

найденной в задаче 2,

краевую задачу

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а2 [игг -\-

u^j + / ( г ,

0 =

“ь

0 < r ,

t < c o ,

 

 

 

 

и (г, 0) = ф (г),

\ и \ < с о ,

г2 = х2-[-г/2-|-г2.

 

 

7.

Найти

распределение температуры в неограниченном простран­

стве,

вызванное тем,

что в начальный

момент t — 0

на

каждой еди­

нице длины бесконечной цилиндрической поверхности радиуса г„ выделилось Q равномерно распределенных единиц тепла. (Построение

функции

влияния

мгновенного

цилиндрического

источника

тепла.)

8. С помощью функции источника, найденной

в задаче 7,

решить

краевую

задачу

 

 

 

 

 

 

a2 (urr+

y u r}-b-f(r,

t) = U(,

0 < r ,

t < c o ,

 

 

и (г, 0) = ф (г), ] и ! < со , г2 = дс2 + у2.

 

9. Построить функцию влияния мгновенного точечного источника на бесконечной прямой для уравнения

а?ихх —-hu — щ.


Г л а в а VII

МЕТОД ФУНКЦИЙ ГРИНА РЕШЕНИЯ КРАЕВЫХ ЗАДАЧ ДЛЯ УРАВНЕНИЙ ЭЛЛИПТИЧЕСКОГО ТИПА

Метод функций Грина, описанный в гл. VI, приме­ няется к решению весьма широкого круга линейных задач. В настоящей главе мы рассмотрим применение его к решению краевых задач для уравнений эллиптического типа.

Простейшим и наиболее часто встречающимся в прило­ жениях уравнением этого типа является уравнение Лап­ ласа Ды = 0. Его решения — гармонические функции — обладают рядом специфических свойств, которые пона­ добятся нам при изучении свойств функций Грина. Поэтому мы рассмотрим прежде всего некоторые свойства гармони­ ческих функций.

§ 1. Вторая формула Грина. Простейшие свойства гармонических функций

Все результаты, к которым мы придем в этой главе, будут получены из небольшого числа формул и соотно­ шений.

Выводом этих формул мы и займемся прежде всего.

1. Пусть функции и(М) и v (М) обладают следующим свойствами:

1) непрерывны вместе с частными производными пер­

вого порядка всюду в замкнутой области D, ограниченной поверхностью S, кроме, может быть, конечного числа точек;

2 ) интегрируемы вместе с частными производными пер­ вого порядка в области D;

3) имеют интегрируемые в области D частные произ­ водные второго порядка.

156

Тогда, как известно,

R [и, у] = — ^ vL [и] dx =

D

= jj £ (V«, Vy) dx +

^ quv dx— i k v ^ d x *).

(1)

D

D

S

 

Здесь L [u] = div (k Vu) — qu. Вычитая R [и, v] из R[v, u],

получим вторую формулу Грина:

J {уТ [ы] - ы! [ у]} Л =

J

(2)

D

S

 

Для одномерного случая вторая формула Грина имеет вид

i

 

 

 

 

§ {vL [u]-uL[v]} dx = k ( ^ v ~ - u j x)

(3)

о

 

 

 

 

Сл е д с т в и е . Пусть L [и] =

div (k V«). Дели для такого

оператора Ь[и] в формуле (2 )

положить

у = 1 , а б каче­

стве и(М) взять решение уравнения L\u]=f(M),

непре­

рывное вместе с частными

производными

первого порядка

в области D = D-\-S, то получим

 

 

$ k t l da = jj f(M)dr.

 

(4)

S

D

 

 

 

Для двусвязной области D', ограниченной двумя кон­ центрическими сферами Sr и Sr, (Rt <CR) с центром в точке М0, формула (2) запишется в виде

§ {vL [и] — иЬ [у]} dx =

=

S

S k {VT r - Utr)dG-

W

 

SR

sRi

 

 

 

Минус перед интегралом no Sr, появился

потому,

что на

поверхности

с.

д

= —

д

 

Sr,

выполняется равенство дп

 

 

З а м е ч а н и е .

При произвольных функциях ф (М) и

f(M) (даже

непрерывных) вторая краевая

задача

может

не иметь решения. Действительно, пусть L[n] = div (kVu). Тогда для решения второй краевой задачи и (М), непре­ рывного вместе с частными производными первого порядка

*) См. формулу (8) гл. IV.

157


в D = D + S,

должно

выполняться соотношение (4).

По­

скольку ^ | s =

cp(M ),

то получим

 

 

^ / (М) dx — ^ k д~ do =

jj кц>(М) do.

(6 )

Ь

 

s

s

 

Таким образом, функции f{M) и ф (М) должны быть свя­ заны соотношением

^ f (М) dx = ^ kq> (М) do.

(7)

оs

Вчастности, если f{M) = 0, то функция ф (М) должна удовлетворять условию

§ /еф do = 0 .

(8 )

s

 

Легко понять физический смысл соотношений (7) и (8), если и (Л4) интерпретировать как стационарную темпера­ туру, a k(M) как коэффициент теплопроводности. Тогда соотношение (7) выражает следующий очевидный факт: для того чтобы существовало стационарное решение, необ­ ходимо, чтобы количество тепла, образующееся в области D за промежуток времени At от действия внутренних источ­

ников, было равно суммарному потоку тепла,

уходящего

через границу области S за тот же промежуток времени.

2.

Функция и (Л4) называется гармонической в области D'

если она

непрерывна в D и в

 

каждой точке области D

удовлетворяет уравнению Лапласа Аи = 0.

что в трех­

Непосредственной

проверкой

убеждаемся,

мерном пространстве

функция

1/гМР является

гармониче­

ской всюду, кроме точки, где гмр = 0 .

 

 

В двумерном случае

функция

In (1/гМР) является гар­

монической всюду, кроме точки,

где гМР 0 .

 

1 /f%^2

В m-мерном пространстве

(ш^= 2) функция

является гармонической всюду, кроме точки Р.

также

Функции 1п(1/гд,р),

1/гмр> 1/г™ ~ 2 называются

фундаментальными решениями уравнения Лапласа.

Применяя формулу (4) к гармоническим в области D функциям (непрерывным вместе с их частными производ­ ными первого порядка в D + S), получим

158