Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Методы математической физики и задачи гидроаэродинамики

..pdf
Скачиваний:
7
Добавлен:
12.11.2023
Размер:
14.42 Mб
Скачать

точке границы тела Г в любой момент времени известна темпе­ ратура; к краевому условию второго рода — когда известно зна­ чение теплового потока; к краевому условию третьего рода — когда на границе тела происходит теплообмен с внешней сре­ дой по закону Ньютона: тепловойпоток через площадку еди­ ничной площади пропорционален разности температур тела и внешней среды:

да

О(U ^среды) Is*

дп S

Здесь сг>0 — коэффициент теплопередачи или теплоотдачи, из­ вестной величиной является a«CpeAbi|s=p.

В задаче о малых поперечных колебаниях мембраны усло­ вие закрепления края, когда отклонения точек края от положе­ ния равновесия в любой момент времени равно нулю и |5= 0 , является частным случаем краевого условия первого рода.

Условие упругого закрепления края мембраны, когда при смещениях края вдоль оси Ou возникают силы упругости, по величине пропорциональные отклонениям и(М, t) точек грани­ цы от начального положения ио{М, 0):

«ч да

— o(u — u0)\st

ПРOu1 М —1 0 ^

где а > 0 — константа, зависящая от свойств материала, приво­ дит к краевому условию третьего рода.

Случай колебаний мембраны со свободным краем приводит к краевому условию

£ 1 - ° '

что является частным случаем краевого условия второго рода. При рассмотрении конкретных процессов может возникать задача со смешанными краевыми условиями, когда на части границы задается одно краевое условие, а на остальной части —

другое краевое условие.

Задачи для нестационарных процессов в ограниченных об­ ластях, в которых дополнительно к уравнению задаются началь­ ные условия и краевые условия, называются начально-краевы­ ми задачами для данного уравнения. В случае нестационарных процессов во всем пространстве Rn задают начальные усло­ вия, и задачу тогда называют начальной задачей. Для стацио­ нарных процессов имеет смысл ставить лишь краевые задачи.

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

Рассмотрим постановку задач для уравнения теплопровод­ ности, волнового уравнения, уравнений Лапласа и Пуассона. Будем использовать следующие обозначения:

R%+1= {(Ж, t) |М s R \ t > 0}, R ÿ * = { ( Aî, t) | M ЯV > 0},

DczRn— ограниченная область, Г — ее граница; S = {(M , t) |Â fe е Г , /^ 0 }, / е С * ф ) — функция, непрерывная в области D вместе со своими производными до k-ro порядка включительно.

Начальные задачи

Задача Коши для уравнения теплопроводности

 

 

ut—Lu-\-ft

(2.22)

где Lu=a2Au или Lu=

—— div (k grad и), состоит

в нахожде-

 

ср

 

нии решения уравнения

(2.22) в R"+1 непрерывного в Я"+1

и удовлетворяющего начальному условию Коши

 

и(Л4,0)=<р(М)

(2.23)

(Здесь и'всюду далее рассматриваются классические решения уравнений, т. е. имеющие непрерывные производные всех тех порядков, которые входят в уравнение, и обращающие уравне­ ние в тождество.)

Задача Коши для волнового уравнения

 

uit= L u - \ - f ,

 

(2.24)

где L u= a2àu или £«==— dlv(rgrad«),

состоит в нахождении

р

 

 

решения уравнения (2.24), непрерывно

дифференцируемого в

и удовлетворяющего начальным условиям Коши

 

и(ЛГ,0)=<р(АО, йДЛ4,0)=ф(Л4) V M G R*.

(2.25)

Начально-краевые задачи

Первая (соответственно вторая или третья) начально-крае­ вая задача для уравнения теплопроводности (2.22) в области Q состоит в нахождении решения этого уравнения, непрерывно­ го (соответственно непрерывно дифференцируемого) в замкну­

той области Q и удовлетворяющего начальному условию

и(Л1,0)=<р(М) V M œ D

(2.26)

и краевому условию первого рода (2.19) [соответственно вто­ рого рода (2.20) или третьего рода (2.21)].

Первая (вторая или третья) начально-краевая задача для волнового уравнения (2.24) в области Й состоит в нахождении

решения этого уравнения, непрерывно дифференцируемого в Й и удовлетворяющего начальным условиям

и(Л4,0)=<р(ЛГ), иДЛ/,0)=ф(Ж )

V M <=D

(2.27)

и краевым

условиям первого рода

(2.19)

[соответственно вто­

рого рода

(2 .20) или третьего рода

(2.2 1 )].

 

 

Краевые задачи

 

 

Внутренняя задача Дирихле (или первая краевая задача)

для уравнения Лапласа (Пуассона)

 

 

 

 

Д и = 0 (Д и = /)

 

(2.28)

в области D состоит в нахождении решения этого уравнения, непрерывного в D и удовлетворяющего на границе области краевому, условию первого рода

й|г—р. (2.29)

Внутренняя задача Неймана (соответственно третья крае­ вая задача) состоит в нахождении решений уравнений (2.28), непрерывно дифференцируемых в замкнутой области D и удов­ летворяющих на границе, области краевому условию второго рода (соответственно третьего рода)

 

да

г= Р (соответственно

 

 

 

(2.30)

 

дп

 

 

 

Пусть GczRn— область, содержащая внешность шара доста­

точно большого радиуса.

 

 

 

Неймана

или

Внешняя задача Дирихле (соответственно

третья

краевая задача)

состоит

в нахождении

решения

урав­

нения

(2.28) в области

G, непрерывного

в замкнутой области

G (соответственно непрерывно дифференцируемого в С) и удов­

летворяющего на границе Г области G краевому условию пер­

вого рода

(2.29) [соответственно

второго

или

третьего

рода

(2.30)].

В тепловых терминах внутренняя задача Дирихле состоит в отыскании стационарного поля температуры и в области D по заданному распределению температуры и на границе Г этой области.

В § 2.2 было показано, что потенциал скорости безвихрево­ го движения несжимаемой жидкости удовлетворяет уравнению Лапласа (2.10). Выясним вид граничного условия, соответст­ вующего безотрывному обтеканию твердого тела потоком иде­ альной жидкости. При безотрывном обтекании вектор скорости

Г лежит в касательной плоскости к поверхности 5 обтекаемого тела и, следовательно, (V, n ° ) |s = 0, где п° — орт нормали к

поверхности 5. В силу потенциальности векторного поля имеем V — —gradu, где и— потенциал скорости, и

да

( V , n ° ) |s = (-g r a d tf,tt°) \ s = —

дп

Мы приходим к краевому условию

да

0.

(2.31)

=

дп S

Таким образом, задача обтекания твердого тела потенциаль­ ным потоком идеальной жидкости есть задача Неймана с гра­ ничным условием (2.31).

О корректности постановки задачи математической физики

Говорят, что задача математической физики поставлена корректно, если начальные и краевые условия обеспечивают:

1 ) существование решения в рассматриваемой области;

2 ) единственность решения (два решения, удовлетворяющие

одним и тем же начальным и краевым условиям, совпадают в области) ;

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

Почему существенно преждё чем решать задачу математи­ ческой физики, выяснить, корректно ли она поставлена? Иссле­ дуя, существует ли решение, мы выясняем, не переопределена ли задача, не являются ли наложенные на решение условия противоречивыми. Исследуя вопрос о единственности решения рассматриваемой задачи, мы выясняем, не является ли задача неопределенной, когда наложенные условия оказываются недо­ статочными для обеспечения единственности решения. Наконец, установление непрерывной зависимости решения от начальных и краевых условий важно по следующей причине. Если коэф­ фициенты уравнения, которое, как правило, выводится на ос­ нове физических законов, являются точными константами или функциями, то функции, входящие в начальные или краевые условия, находятся приближенно из эксперимента. Поэтому задача, в которой малая ошибка в экспериментальных данных вызывает большие изменения решения, практически бесполезна.

Решение, непрерывно зависящее от начальных и краевых условий, называется также устойчивым, а задача — устойчи­ вой.

Не всякая задача корректно поставлена. Примером некор­ ректно поставленной задачи является задача Коши для урав­

нения Лапласа относительно и(х, у) в полуплоскости у^>Q.

UXX~\~Uyy=0y

и (х,0)= у(х), иу{х>0)==tU).

Убедимся, что эта задача, неустойчива. Действительно, рассмот­ рим решения двух задач:

1) Ди=0,

и (*,0)= 0, tttf(jc,0)=0;

 

%

 

2) Дих= 0 ,

HX(X; , 0 ) = ~ COSXA:, ,^ и^ х ’ —Q.

 

X

ду

При достаточно больших по модулю X начальные функции в

них мало отличаются:

 

 

\и(ху0)—их(л:,0)| =

— cos

 

 

X

С другой стороны, их решения

I

й = 0 и их= — eXÿcosXjt при

Л

* = 0 , У=Уо’> 0 отличаются на величину

I и (0, у0)- их(О, У)0|= 4 l eXÿ° »

|Л|

превосходящую любое наперед заданное число Л1 ;> 0 при соот­

ветствующим образом выбранном Х^так как lim ~ - e x*'« = ooj.

Приведенный пример некорректно поставленной задачи назы­ вается примером Адамара.

Некорректной также является задача Коши для гиперболи­ ческого уравнения с начальными данными на характеристике, так как в этом случае нарушается единственность решения.

ГЛАВА 3

МЕТОДЫ РЕШЕНИЯ ЛИНЕЙНЫХ НАЧАЛЬНО-КРАЕВЫХ «ЗАДАЧ

$ 3.1. ПОЛНОСТЬЮ НЕОДНОРОДНАЯ НАЧАЛЬНО-КРАЕВАЯ ЗАДАЧА. РЕДУКЦИЯ.

МЕТОД ДЮАМЕЛЯ ДЛЯ НЕОДНОРОДНОГО УРАВНЕНИЯ

Рассмотрим в области

Q ={x=(xi,

..., xn)^D czR n, /> 0 }

уравнение относительно и(х,

t)

 

M u = L u + f

(3.1)

идля него начально-краевую задачу с начальными условиями

вобласти D в момент /= 0

Ltu=(f

(3.2)

35

и краевыми условиями на 5 = { * е Г , t^ . 0}, Г — граница облас­ ти D,

Lst t = H-,

(3.3)

где М, L,, Lt, Ls — линейные дифференциальные

операторы:

М, Lt — операторы дифференцирования по времени

t; L, Ls

операторы дифференцирования по пространственным перемен­ ным Х\, ..., хп. В случае, например, начально-краевой задачи

(2.24), (2.27), (2.21) для волнового уравнения имеемМи= —— ,

Lu=a?àu или Lu— I div (Г grad й), начальное условие £*й=<р

?

представляет собой систему двух условий w l^o— ф, ш |Д«о=ф,

краевое условие Lsti=\i имеет вид

= 1», а > 0 .

Полностью неоднородную задачу

(3.1) — (3.3) сводят к бо­

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

Исходную задачу редуцируют прежде всего к задаче с одно­ родными (нулевыми) краевыми условиями. Для этого подби­ рают достаточно гладкую функцию U, удовлетворяющую крае­ вому условию (3.3), и вычитают ее из решения. Тогда для раз­ ности й—и—U в силу линейности операторов М, L, Lt, Ls по­ лучаем

М и=М и —M U=Lu-\- f —M U =Ltt-\-LU - \- f —M Ut Ltu —Ltu — LtU =<p—LtU t Lyi—Lyi—LsU —(i—p = 0 ,

или, полагая ]= f+ L U —MU, <p=<p—LtU, имеем задачу

 

M u = L u + f, Ltiï= y , Lsu = 0 .

(3 .4)

Далее,

представляя й в виде суммы v+ w , подставляя сум­

му в (3 .4)

и используя линейность операторов М, L, Lt, Ls,

получаем

 

 

M v-j-Mw= bo+ Lw -f / ,

 

Ltv-\-Ltw = ip,

 

Lsv + L sw = 0,

откуда следует, что в качестве v a w можно взять решения

36

следующих задач:

A) M v= L v,

Б) M w=Lw -{-/7

 

I

 

 

Ltv*= ср,

Ltw = 0,

(3.5)

Lsv = 0;

Lsw = О.

 

Задача типа Б в

(3.5)

может быть сведена к задаче типа А.

Действительно,

пусть

М = —

— г *

Рассматри-

 

 

 

dtk

dt}

 

ваемый в задаче

типа

Б

процесс,

например колебательный,

протекает под влиянием непрерывно действующего возмущения J=f(x, t), которое состоит из следующих друг за другом мгно­ венных возмущений ](х, х), действующих в моменты т е [ 0, /].

Процесс, протекающий в последующее время / > т под влияни­

ем мгновенного возмущения J(x, х),

описывается решением

z(x, t, х) следующей задачи типа А:

 

M z= L z,

 

LitXz = y x,

(3.6)

LsZ=0,

в которой начальное условие записывается в развернутом виде следующим образом:

'LUxz = ( z (# ,* ,т), z \(x ,t,x ),...,zj*-1*(л*, ü,т)) |/-т=?:

==s(<Po^»*0» Yi (xy*0»

1 (^»‘^))==(б* 6» •••» 6» f (x,T)).

(3.7)

В силу линейности задачи

(3.6) в момент времени t

процесс,

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

ний во все предшествующие

моменты т е [ 0, /],

описывается

функцией

t

 

 

 

vix , t) =

j1z (x , /, T) dt,

(3.8)

 

0

 

и, следовательно, можно предположить, что эта функция долж­ на быть решением задачи типа Б. Покажем это. Замечая, что

согласно

(3.7) z(x, t, t) —z/(x,

tf /)= ...= z \*~2)(x,

t,

t ) —0, a

z}*-1* (x,

t, t)=J(x, t), для производных от. v(x, t)

no

t полу­

чим следующие выражения:

 

 

 

 

t

t

 

 

 

z t (x,

t) d t = J z\ (x , t , t) dt,

 

*(*-1) (*, ü)=z<*~2) U . *. 0 +

Ç 2Sft-1) (*»

T)dt=

J * ? _1) <■*» *» T)dt»

 

О

 

 

о

 

 

 

 

 

 

(3.9)

 

t

 

 

__

 

■п<*>U , t ) = z < f - v (JC, /, 0 + J 2(/ft) (дг, /, T) d t= 7 U , 0 +

/

0

 

 

 

 

 

T) dt.

 

 

4* J 2}ft) (JC,

 

 

Тогда

 

ft-i

/

 

 

 

 

 

 

Ж®=/(дг, 0 + f */ft>U. *. t) d t+ 2

JJ),

d t=

b] \ z*

<«*»*»

G

 

 

о

 

 

= / ( ^ э О + | M zdx= = J+ ^ L züx = = /+ !* {[

2 d t j = 7 + ^ ,

T . e. о — решение уравнения Mv = Lv+f. Полагая

в (3.8) и'в

(3 .9) / = 0, получаем v(x, 0)=vt(x,

0) = ...=

 

0) = 0, т. е.

L,t>=0. Наконец, так как Lsz 0, то Lsu—0.

 

не меняет

Заметим, что сдвиг по

времени

t'= tг, х'=х,

уравнения в задаче (3.6),

а начальные

условия в

момент t = х

преобразует в начальные условия при t'—0.

Таким образом/

функция

t

 

 

 

 

 

 

 

 

 

w (x ,t)—^ z {x ,t —т, t)d t,

 

(ЗЛО)

о

где z(х, t, х) — решение задачи

M z= L z ,

I , z = (О,..., 0 ,7 (•*»*))>

(3.11)

Lsz = О,

дает решение задачи типа Б в (3.5). Этот метод представления решения задачи типа Б для неоднородного уравнения с нуле­ выми начальными и краевыми условиями в виде интеграла (3.10), где z(x, t, х ) — решение импульсной задачи (ЗЛ1), на­ зывается методом Дюамеля.

§ 3.2. МЕТОД ФУРЬЕ РЕШЕНИЯ НАЧАЛЬНО-КРАЕВОЙ ЗАДАЧИ

ДЛЯ ОДНОРОДНОГО УРАВНЕНИЯ С ОДНОРОДНЫМИ КРАЕВЫМИ УСЛОВИЯМИ

Рассмотрим основные идеи метода.

1.Так как линейная комбинация решений однородного уравнения, удовлетворяющих нулевым краевым условиям, есть также решение данного уравнения, удовлетворяющее нулевым краевым условиям, то сначала находят бесчисленное множест­ во таких решений и из них конструируют решение, удовлетво* ряющее и начальным усло­ виям.

2.Отыскание указанного

выше множества решений сво­ дят к решению системы более простых краевых и начальных задач для функций меньшего числа переменных, представ­ ляя решение в виде произведе­ ния функции только времени t на функции только пространст­ венных переменных Хи

Рассмотрим сначала решение этим методом задачи о сво­ бодных колебаниях струны с закрепленными концами

utt= a 2uxx в 2 = { 0 < ^ < / , / > 0},

я |,-о=<?(*), а/|/-о=Ф (л) V *< =[0,/],

 

я |х- о = а |х«*/=0 V / > 0 .

 

(3.12)

Положим и(х,

t) = T(t)X(x),

другими словами,представим

ре­

шение в виде «стоячей волны»

(рис. 3.1). Подставляя это

про­

изведение в уравнение

и

краевые условия,

получим

 

 

 

 

ХГ'=а?Х"Т;

 

(3.13)

X

( 0 )

Т

X{ l ) T{ t ) = О W

> 0.

(3.14)

В уравнении (3.13) разделим переменные. Для этого разделим его на а2ХТ

Т" _ X"

а*Т ~ X

Левая часть полученного равенства зависит только от t, а пра­ вая— только от х. Это возможно лишь тогда, когда обе части равенства тождественно равны константе. Обозначим эту кон­ станту через —Я/. Тогда равенства

Т" _ _ Х " _ ^

c f l T ~

X

1

 

 

 

Т"-\-а2\Т = 0 ;

 

 

 

 

(3.15)

 

 

 

аг" + ха: = о.

 

 

 

 

(3.16)

 

 

 

 

 

 

 

 

Далее, так

как

мы ищем

ненулевые

решения и(х,

t),

то

Т(У)Ф 0, и

равенства

(3.14)

возможны

лишь

в

том

случае,

когда

 

 

 

 

 

 

 

 

 

 

 

«,

 

 

 

 

 

X (0 )= 0 ,

Х (1 )= 0.

 

 

 

(3.17)

Краевая задача

(3.16), (3.17) не при всяких значениях X

име­

ет ненулевое

решение. Действительно,

при

К О

характеристи­

ческое уравнение А2+Х =0 уравнения

(3.16)

имеет два

различ­

ных действительных

корня

kit2= ±

Y

— X

и

общее

решение

уравнения

(3.16)

имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

X (х)=

 

- f c2z~v:=lx.

 

 

 

 

Согласно краевым условиям

(3.17) получаем систему Ci + C2= 0,

£1е/-л/_|_£2е-^ -Л,= 0 ,

которая

имеет

единственное

нулевое

решение Ci=c2= 0. Следовательно,

при ХфО краевая

задача

(3.16) , (3.1?)

не Имеет ненулевых решений.

 

решение Х(х)=г=

При Х= 0

уравнение (3.16)

имеет

общее

= С\Х+С2 и краевым условиям

(3.17)

можно

удовлетворить

лишь при

Ci=C2= 0,

т. е. и

в

этом

случае

краевая

задача

(3.16) , (3.17) не имеет ненулевых решений.

При Х>0 характеристическое уравнение k2+X имеет мни­

мые корни k\t2= ± i п и общее решение уравнения (3.16) име­

ет вид

X (jr)=£i cos ÿr\x-j-c2 sin Y \x .

Согласно краевым условиям (3.17) получаем систему

Ci= 0 , с2 sin У X/= 0 ,

которая имеет ненулевое решение ci = 0, с2=^0

sin

|^ Х /= 0.

Отсюда

^ |щ 4 ÏS

>* II а

и в этом случае

 

^ n (^ )= sin

п — 1 , 2 ,....

лишь когда

(3.18)

(3.19)

(Можно положить С21, так как решение задачи (3.16), (3.17)