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

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

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

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

Добавлен: 16.10.2024

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

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

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

Будем считать, что распределение нейтронов по пространству имеет вид (близко к равновесному) ¥ (х, р) = а (х) -|- b (х) р. Здесь второе слагаемое носит характер поправки к первому (Ь<^.а). В этом приближении

1

П М = § (a-f6p)p*dp = ^ [ l - ( - l ) *

+4 + ^

[1_(_1)Л4*].

- 1

 

 

 

 

 

 

Отсюда

 

 

 

 

 

 

,=

2a,

¥ х = -|& ,

¥ 2=

3 v°*

 

 

 

 

 

 

 

=

0, 1,

получим

 

 

 

 

^

=

( P - a ) ¥ 0,

d¥a =

■ a¥ х.

(29)

dx

 

 

 

 

 

Подставим в (29) ¥ 2 =

¥ 0. Получим

 

 

 

* Л - (а _ а)Цг

-1

^

■«¥,.

(30)

d*

“ (Р а) V°’

3

dx

 

 

Мы получили систему двух уравнений для двух функций ¥ 0 и ¥ х. Удобно перейти от функций ¥„ и ¥ х к концентрации нейтронов п и плотности потока нейтронов j в направлении х. Имеем

1 1

п = ( ¥ (х, р) dp — ¥ 0, / = Г ор¥ (х, р) dp = о¥х.

- 1 - 1

Перепишем систему (30), вводя функции п и / вместо ¥ 0 и ¥ х:

dl.

v($ — a)n,

v dtt

 

dx

 

 

3 Tx = ~ a1'

или

dn

v

d2n

 

v

v ((3 —a) n.

3a d x ’

3a dx2

 

Положим v/(3a) = D. Тогда первое из уравнений дает закон Нернста, где D — коэффициент диффузии, а второе —уравнение диф­ фузии в одномерном стационарном случае:

О щ " + » < Р - с )» = 0.

(31)

Член v (Р — а) п представляет источники.

Заметим, что полученные уравнения справедливы и в анизотроп­

ном случае,

когда а<С Ь, но ¥

=

а + 6р.

 

Если в рассмотренном приближении повторить весь

еывод для

пространственного нестационарного

случая, то мы придем к системе

уравнений,

аналогичной (31):

 

 

 

 

 

/ = — H gradn,

^

=

D A ft+ o (P - a )fi,

(32)

36


§ 8. Типы краевых условий. Постановка краевых задач

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

а) написать уравнение (или систему уравнений), кото­ рому удовлетворяет искомая функция (или система функ­ ций), описывающая исследуемое явление;

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

При решении каждой конкретной физической задачи математическими методами надо ставить вопрос не о реше­ нии соответствующего уравнения, а о решении математи­ ческой задачи в ее полной постановке, вместе с соответ­ ствующими дополнительными условиями. Обычно интере­ сующие нас явления имеют однозначный характер, в то время как описывающие их уравнения имеют множество решений. Поэтому при математической постановке задачи недостаточно написать уравнение (или систему), которому удовлетворяет искомая функция (система функций). Надо также указать дополнительные условия, позволяющие выделить лишь одно интересующее нас решение, описы­ вающее конкретное явление, процесс. Дополнительных условий должно быть «не слишком много», чтобы суще­ ствовало решение, удовлетворяющее им. Таким образом, дополнительные условия должны обеспечить существование и единственность решения.

Характер дополнительных условий мы покажем на примерах задач, рассмотренных в предыдущих параграфах.

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

или стержня

(уравнения (3) и (5)) надо задать начальный

профиль

и (х, 0) = ср (х)

 

и начальную скорость

 

щ (х, 0) = ф (х)

 

точек струны (стержня).

Это начальные условия. Аналогичный вид они имеют

для любого волнового уравнения.

на концах

(краях)

Кроме того, надо

записать режим

струны (стержня).

 

 

и х = 1)

Так, если задан закон движения концов (х = 0

и (0 , /) =

М 0 . и (1>0 =

М-2 (0 *

 

37


то мы будем называть такие дополнительные условия

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

Если задан закон изменения силы, приложенной к концу струны (стержня) и действующей в направлении колеба­ ний, то режим на концах можно записать следующим образом:

Еих \х~о= /1 (0>

Еих U_/ = /2(0

или

их (/, t) = va (t).

их (0 , t) = vx (t),

Это краевые условия второго типа.

Пусть к концу стержня

{х = 1) прикреплена пружина,

действующая вдоль оси х. Тогда сила натяжения Еих на конце будет уравновешиваться силой действия пружины, равной сш. Краевой режим на этом конце можно записать следующим образом:

Eux (l, t) = — сш(/, t),

где. а — коэффициент жесткости пружины, или

 

их (/, 0 +

hu (/, t) =

0 .

 

 

Если

пружина в

свою

очередь

движется

по закону

x = $(t),

то краевой

режим

запишется в виде

 

 

ux (l, t) + h[u{l, t) —р (/)] =

().

 

Это краевое условие

третьего типа.

На

левом конце

(л: = 0) оно запишется

в виде

 

 

 

 

их (0 , /) — h \и (0 , t) р (/)] =

0 .

 

Для дву- и трехмерного случаев рассмотренные типы краевых условий имеют следующий вид:

 

 

и s

= р ( М ,

t)

(первый тип),

(33)

 

 

ди

 

t)

(второй тип),

(34)

 

 

дп 5 = v(M,

 

ди

'Sn + hu) S

= P (

Mt),

(третий тип).

(35)

Здесь

производная • по

внешней нормали к

поверх-

дп

ности

S.

 

 

 

 

 

Такие же краевые

условия

встречаются и в задачах,

приводящих к уравнениям параболического типа. Так, если задается температура на поверхности тела, то имеем краевое условие первого типа. Если задается плотность

38


*

потока тепла — , ди через поверхность тела S, то имеем

краевое условие второго типа. Если же на поверхности тела происходит теплообмен со средой, имеющей темпера­

туру Р (М, t), по закону Ньютона — k

= h [и — Р] \s, -то

имеем краевое условие третьего типа.

условий. Неко­

Встречаются и другие типы краевых

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

ваются однородными,

если правые

части

(u, v, Р) тожде­

ственно равны нулю,

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

Очевидно, такие

же краевые

условия

встречаются и

в задачах, приводящих к уравнению эллиптического типа. Физическое истолкование каждого из них не представляет никаких трудностей.

Краевые условия определяются физической постановкой

задачи

и могут иметь разнообразный характер.

В част­

ности, они могут быть и нелинейными.

 

Теперь мы приведем постановку соответствующих трех

типов

краевых задач для уравнений вида

 

 

 

div (k Vw) — qu + / ( M,

/) = putt

(36)

И

 

div(k Vu) — qu-\-f(M,

t) = pu(,

(37)

где k,

q, p —функции точки M.

 

 

Все рассмотренные нами выше уравнения принадлежали

к этому виду

с q = 0.

 

 

П е р в а я

к р а е в а я задача. Найти функцию и (М , t),

удовлетворяющую в области B =

t> 0)

уравне­

нию (36) (соответственно (37)) и дополнительным условиям:

а) начальным

 

 

 

 

и(М, 0) = ср(М),

щ(М, 0) = ф(М)

для M e D

(соответственно и(М,

0)

= ф(М)),

 

 

б) краевым

 

 

 

 

и(М, t)\s — \a(M, t)

для

^ > 0 .

Вторая (третья) краевая

задача ставится аналогично с за­

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

первого

типа

условием второго

(третьего) типа.

рассмотренные типы краевых усло­

З а м е ч а н и е . Все

вий можно записать одним соотношением

{T i(M )^ + y2 ( M) n } 5 = P(M, t).

39