Главная » Просмотр файлов » Н.А. Боголюбов - ММФ (методические материалы)

Н.А. Боголюбов - ММФ (методические материалы) (1125160), страница 13

Файл №1125160 Н.А. Боголюбов - ММФ (методические материалы) (Н.А. Боголюбов - ММФ (методические материалы)) 13 страницаН.А. Боголюбов - ММФ (методические материалы) (1125160) страница 132019-05-11СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 13)

Найти электростатический потенциал внутри кругового цилиндрарадиуса a , основание которого z = 0 и z = h заземлены , а на боковойповерхности поддерживается потенциал cos2φ .Решение:Для электростатического потенциала u имеем следующую краевую задачувнутри цилиндра :Δu = 0,u z=0 = uur=az=h= 0,= cos2φ .Общее решение поставленной задачи имеет следующий вид:u=∞∞∑∑In( h r)πkπkn=0 k=1 In( a)hsinπkz{Ank cosnφ + Bnksinnφ} ,hТак как граничное условие содержит только cos2φ :Ank = 0 n ≠ 2Bnk = 0 при всех n и k .Таким образом решение будет представимо в следующем виде:∞I2( h r)πkπku = cos2φA2ksin z .∑πkhk=1 I2( a)h145Коэффициенты A2k определяются из граничного условия:A2k =h1πksin h z2∫1 ⋅ sin0πk2k(−1)zdz =1−.{}hπkТаким образом, решение поставленной задачи имеет следующий вид:u==∞2cos2φ∑πk=1∞k(−1)1−{} I2( h r)I2(πkI2(kπ(2m + 1)r) sin4hcos2φ∑π(2m + 1)aπm=0 I2()hπkah )sinπkz=hπ(2m + 1)zh.2m + 12.

Найти распределение потенциала внутри прямого кругового цилиндрарадиуса a , на торцах которого задано нулевое электрическое поле, а набоковой поверхности поддерживается потенциал, равный U0 z , U0 = const .Решение:Электростатический потенциал u является решением следующей краевойзадачи :Δu = 0,∂u∂zuz=0r=a=∂u∂z= U0 z .z=h= 0,В этом случае решение от переменной φ не зависит и его можно записать вследующем виде:146u=∞∑I0( h r)πkAk0cosπkk=0 I0( a)hπkz.hКоэффициенты Ak0 определяются из граничного условия:Ak0 ==h1πkcos h z2∫U0 zcos0πkzdz =hU0 h2 , k = 0−U02h1π 2k2 [− (−1)k], k = 1,2, .

. . , ∞Следовательно, решение поставленной задачи имеет следующий вид:I0( h r)∞2U0h1k(−1)u = U0h −1−[]∑ π 2k 22Ik=1πkπk0( h a)cosπkz .h3) Решить уравнение Лапласа внутри сектора кругового цилиндраπ0 ≤ r ≤ a, 0 ≤ φ ≤ , 0 ≤ z ≤ h с граничными условиями :4ur=a=uuφ=0z=h=uφ= π4=uz=0=0 ,= r 4sin4φ .Решение:Так как при z = 0 задано нулевое граничное условие, то решение можнопредставить в следующем виде:147∞∞∑∑u=J4nn=0 k=1()λk4n aJ4n()λk4n rshλk4n zshλk4n hλk4n - корни уравнения:Anksin4nφ , где=0.Необходимо учесть, что разложение решения проводится по собственным функциямπзадачи Дирихле для сектора с углом раствора α = , которые выражаются следующим4образом:λk4n r sin4nφ, n = 0,1, . .

. , k = 1,2, . . .()J4nКоэффициенты ряда определяются из граничного условия при z = h:1Ank =()J4na⊗π4∫∫00=λk4n rsin4nφ2λk4n r sin4φsin4nφrdrdφ =()r 4 J4naδn1J′λk4n a2 ( 4) ()a2⊗22∫5(r J40λk4 r dr)= δn12a 3J5 a(λk4)λk4 (J′4)2λk4 a().Таким образом:u=∞∑k=12a 3J5 a(λk4 (J′4)2 a(λk4)λk4)148shλk4 zshλk4 hλk4 r sin4φ .()J44. Решить уравнение Лапласа внутри тора прямоугольного сечения :a ≤ r ≤ b, 0 ≤ φ ≤ 2π, 0 ≤ z ≤ h с граничными условиями :∂u∂zuz=0r=a=∂u∂z= 0, uz=h= 0,r=b= cos3φ .Решение:Граничные условия при z = 0 и z = h являются однородными граничнымиусловиями второго рода, а при r = a и r = b граничные функции не зависят отz.Следовательно : рассматриваемая задача вырождается в краевую задачу дляуравнения Лапласа внутри кольца a ≤ r ≤ b :Δ2u = 0, a ≤ r ≤ b, 0 ≤ φ ≤ 2π{ur=a= 0, ur=b= cos3φ .Решение данной задачи можно записать в следующем виде:∞A0 R0(a)(r)Rn(a)(r)u=+{Ancosnφ + Bnsinnφ} ,(a)2 R0(a)(b) ∑Rbn=1 n ( )гдеR0(a)(r)r 2n − a 2nr(a).= ln , Rn (r) =nraИз граничных условий находим:Bn = 0, n = 1,2, .

. . ,An = 0, n ≠ 3,A3 = 1 . b3 r6 − a6⇒u= 3 6cos3φr b − a6149.4. Задача Штурма - Лиувилля и уравнение Лапласа : шар и шаровой слой.a) Собственные функции шара.Построим собственные функции шара Ka . Введем сферическую системукоординат (r, θ, φ), 0 ≤ r ≤ a, 0 ≤ θ ≤ π, 0 ≤ φ ≤ 2π с началом в центрешара.Оператор Лапласа в сферической системе координат имеет следующий вид:Δu =1 ∂12 ∂ur+Δθφu ,22()r ∂r∂rrгде Δθφu - сферический оператор Лапласа, который равен:1 ∂∂u1 ∂2 uΔθφu =sinθ+.sinθ ∂θ (∂θ ) sin2θ ∂φ 2Задача Штурма - Лиувилля для шара будет выглядеть следующим образом:Δu + λu = 0, M ∈ Ka,α∂u∂r+ βuu ≠ 0.r=a= 0,которую мы будем решать стандартным методом разделения переменных ,отделяя ее радиальную переменную r :u = R(r)v(θ, φ) ≠ 0 .Подставляем наше решение в исходное уравнение, записанное всферической системе координат и разделив переменные, мы получим:d2 dRrdr (dr )+ λr 2 RR(r)≡−Δθφvv(θφ)=μ .На собственные функции рассматриваемого уравнения мы должныналожить условия ограниченности в Ka и условия периодичности по φ спериодом 2π .150Тогда, задача Штурма - Лиувилля для функции v будет выглядетьследующим образом:Δθφv + λv = 0, 0 < θ < π, 0 < φ < 2π,v(θ, φ) ≡ v(θ, φ + 2π),v(0, φ) < ∞, v(π, φ) < ∞,v(θ, φ) ≠ 0 .Cобственные функции полученной задачи являются сферическимифункциями следующего вида:v = vnm(θ, φ) =Yn(m)(θ, φ)=(m)(cosθ)P(n){sinmφ,cosmφ,Cобственные значения:μ = μn = n(n + 1), n = 0,1, .

. . , ∞, m = 0,1, . . . , n .Для каждого μ = n(n + 1) из полученного уравненияd2 dRrdr (dr )+ λr 2 RR(r)≡−Δθφvv(θφ)= μ , мы можем получить уравнение для радиальнойчасти R(r) :1 dn(n + 1)2 dRr+λ−R=0.22()()r drdrrЕго решение должно удовлетворять граничному условию , при r = a :αdR+ βRdrr=a= 0,и условию ограниченности при r = 0 :151R(0) < ∞ .Сделаем замену: R =y(r). Тогда задача для радиальной части будетrсведена к очередной задаче Штурма - Лиувилля следующего вида:r y″ + r y′ + λr − (n +[22αy′ + (β −αy2a ) r=a12)2y = 0,]= 0,y(0) < ∞ .Общее решение будет иметь вид:y = C1Jn+ 12 ( λr) + C2 Nn+ 12 ( λr) .Если мы учтем поведение функций Неймана в нуле и условиеограниченности, то мы можем определить C2 = 0 .Пусть C1 = 1 . Тогда дисперсионное уравнение для λ примет вид:αa λJ′n+ 12 ( λa) + β −J 1λa) = 0 .(2a ) n+ 2 (Обозначим μ = a λ .Тогда радиальная функция R(r) :R(r) = Rnk(r) =где μk1(n+ 2 )Jn+ 12(μk(n+ 12 )arr), n = 0,1, .

. . , k = 1,2, . . . ,- k - ый корень уравнения :ααμJ′n+ 12 (μ) + βa −J 1 (μ) = 0 , при фиксированном n = 0,1, . . .(2 ) n+ 2152Следовательно, собственная функция шара , имеет следующий вид:unkm(r, v, φ) =11μkn+ 2Jn+ 12rar Yn(m)(θ, φ); n = 0,1, . .

. . , k = 1,2, . . . ,1где μkn+ 2 - корни рассматриваемого уравнения.Собственные значения шара равны:λnk =2n+ 12μk1, где μkn+ 2 - корни рассматриваемого уравнения.aКаждому собственному значению λnk соответствуют 2n + 1 линейно независимыхсобственных функций (rangλnk = 2n + 1) .Вычислим норму собственных функций:unkm2∫=2unkmdV1=rKaгдеYn(m)2Теперь вычислимa2=2= πεmJn+ 122 (n + m)!,2n + 1 (n − m)!1r2Jn+ 12(J′n+ 12 ) (a λ ) + 1 −22Yn(m)а εm =2,2, m = 0,{1, m ≠ 0 .a2= Jn+λr)rdr =∫ 12 (0(n +12)2a2 λ1532Jn+a λ)1(2.Решения 1,2,3 краевых задач в шаре:(задача Дирихле, Неймана и Робена)1.

Задача Дирихле (α = 0, β = 1) :Собственные значения определяются уравнением:Jn+ 12 (μ) = 0, λ =n+ 12μka2,1r2Jn+ 1212a2n+ 12 )(= (J′n+ 12 ) μk.()22. Задача Неймана (α = 1, β = 0) :Собственные значения определяются уравнением:μJ′n+ 12 (μ) −1r1J 1 (μ) = 0 ,2 n+ 22Jn+ 12a2n(n + 1)n+ 12 )2(.=1−Jn+ 1 μk2()221(n+ 2 )2μk[]3. Задача Робена (α = 1, β = h) :Собственные значения определяются уравнением:1μJ′n+ 12 (μ) + ah −Jn+ 12 (μ) = 0, λ =()2154n+ 12μka2Выражение для квадрата нормы можно записать двумя способами:1)2)1r1r2Jn+ 1232Jn+ 123a2=22a=21−n(n + 1) + ah(1 − ah)1(n+ 2 )μk[[]1(n+ 2 )]μk1+4221(n+ 2 )2Jn+ 1 μk2(),(удобна при малых h )− (n +(1 − 2ah)2212)(J′n+ 12 )1(n+ 2 )μk()2(удобна при больших h )б) Собственные функции шарового слоя.Задача Штурма - Лиувилля для шарового слоя D (a < r < b) :Δu + λu = 0,P1[u] = α1 ∂u− β1u∂rP2[u] = α2 ∂u+ β2u∂rr=a= 0,r=b= 0,u(r, θ, φ) ≠ 0, αi + βi ≠ 0, i = 1,2 .Представляем решение в виде:u = R(r)v(θ, φ) .Подставляем его в исходное уравнение, разделяем переменные и получаемзадачу Штурма - Лиувилля для угловой функции v(θ, φ) на сфере r = a :155Δθ,φv + χ v = 0 0 < θ < π, 0 ≤ φ ≤ 2π,1)vθ=0θ=π< ∞,(угловая задача)v(θ, φ) ≡ v(θ, φ + 2π),v(θ, φ) ≠ 0 .1 d2 dRrdr )r 2 dr (2)+ (λ −r2 )P1[R] ≡ α1 dR− β1Rdrr=aP2[R] ≡ α2 dR+ β2 RdrχR = 0,r=b= 0,(радиальная задача на отрезке)=0R(r) ≠ 0, αi + βi ≠ 0, i = 1,2 .Cобственные функции угловой задачи являются сферическими функциями:v = vnm(θ, φ) = Ynm(θ, φ).Собственные значения :χ = χn = n(n + 1), n = 0,1, .

. . , m = 0,1, . . . , n .Тогда общее решение радиальной задачи при χ = n(n + 1) примет вид:R = C1Jn+ 12 ( λr)r+ C2Nn+ 12 ( λr)r.Подставляем общее решение в граничные условия и получаем следующуюсистему:C1 p1(λ, a) + C2q1(λ, a) = 0,{C1 p2(λ, b) + C2q2(λ, b) = 0,156где:p1(λ, a) = α1 λJ′n+ 12 ( λa) − (β1 +q1(λ, a) = α1 λN′n+ 12 ( λa) − (β1 +J 12a ) n+ 2 (α1N 12a ) n+ 2 (λa),α1p2(λ, b) = α2 λJ′n+ 12 ( λb) + (β2 −J 12b ) n+ 2 (α2q2(λ, b) = α2 λN′n+ 12 ( λb) + (β2 −N 12b ) n+ 2 (α2λa),λb),λb),Приравниваем к нулю определитель полученной системы и получаемдисперсионное уравнение для определения собственного значения λ :p1(λ, a)p2(λ, b)=q1(λ, a)q2(λ, b)Теперь мы можем выразить из нашей системы C2 :p1(λ, a).C2 = − C1q1(λ, a)Пусть C1 = q1(λ, a) , тогда согласно полученному общему решению , мыможем записать собственную функцию в следующем виде:R=Jn+ 12 ( λr)rq1(λ, a) −Nn+ 12 ( λr)rp1(λ, a).Также , мы можем учесть полученное дисперсионное уравнение ипереписать собственную функцию в виде:R=p1(λ, a)p2(λ, b)Jn+ 12 ( λr)rq2(λ, b) −157Nn+ 12 ( λr)rp2(λ, b).Таким образом, собственные функции шарового слоя выражаются:Jn+ 12 ( λr)unkm(r, θ, φ) =rq1(λ, a) −Nn+ 12 ( λr)rp1(λ, a) Yn(m)(θ, φ),угловаярадиальнаяn = 0,1, .

. . , m = 0,1, . . . , n, k = 1,2, . . . ,где λ =1(n+ 2 )λk- k - ый корень рассматриваемого уравнения при каждомфиксированном n .Теперь вычислим квадрат нормы:b π 2π2unkm=∫∫ ∫2(r)Y n(m) (θ, φ)r 2sinθdθdφdr =Rnk2a0 0b2∫{aJn+ 12 ( λr)q1(λ, a) − Nn+ 12 ( λr)p1(λ, a) rdr ⊗}2π π⊗∫∫Y n(m) (θ, φ)sinθdθdφ = = Rnk22Yn(m)2.0 0Квадрат нормы радиальной функции Rnk вычисляется аналогично.Ее окончательная формула имеет следующий вид:Rnk22= 2π2[ p2(λ, b) ]p1(λ, a)12)n+(α2 2β2 ++ α22 1 −(2b )λb 212)n+(α1 2− β1 ++ α12 1 −(2a )λa 22,158где λ =2−n+ 2λk( )1.в) Частные решения уравнения Лапласа в сферической системе координат.Для построения решения решения краевой задачи для уравнения Лапласа вшаре, вне шара и в шаровом слое, необходимо получить специальные решения,которые мы будем называть шаровыми функциями .Давайте сначала построим эти решения: введем сферическую системукоординат (r, θ, φ) и найдем решения уравнения Лапласа в следующем виде:u(r, θ, φ) = R(r)v(θ, φ) .Подставляем наше решение в уравнение Лапласа , разделяем переменные иполучаем:d2 dRrdr (dr )R≡−Δθφvv=λ ,1 ∂∂v1 ∂2 vгде : Δθφv =- cферический операторsinθ+ 22()sinθ ∂θ∂θsin θ ∂φЛапласа .Таким образом, мы получаем краевую задачу для определения угловой частиv(θ, φ):Δθφv + λv = 0, 0 < θ < π, 0 ≤ φ ≤ 2π,v(θ, φ) = v(θ, φ + 2π),v(θ, φ)θ=0,πзадача Штурма - Лиувиллядля< ∞, v ≠ 0сферического оператора Лапласа.и уравнение для радиальной части R(r) :r 2 R″ + 2rR′ − λR = 0 .159Рассмотрим угловую задачу .

Характеристики

Тип файла
PDF-файл
Размер
2,57 Mb
Тип материала
Высшее учебное заведение

Список файлов книги

Свежие статьи
Популярно сейчас
Как Вы думаете, сколько людей до Вас делали точно такое же задание? 99% студентов выполняют точно такие же задания, как и их предшественники год назад. Найдите нужный учебный материал на СтудИзбе!
Ответы на популярные вопросы
Да! Наши авторы собирают и выкладывают те работы, которые сдаются в Вашем учебном заведении ежегодно и уже проверены преподавателями.
Да! У нас любой человек может выложить любую учебную работу и зарабатывать на её продажах! Но каждый учебный материал публикуется только после тщательной проверки администрацией.
Вернём деньги! А если быть более точными, то автору даётся немного времени на исправление, а если не исправит или выйдет время, то вернём деньги в полном объёме!
Да! На равне с готовыми студенческими работами у нас продаются услуги. Цены на услуги видны сразу, то есть Вам нужно только указать параметры и сразу можно оплачивать.
Отзывы студентов
Ставлю 10/10
Все нравится, очень удобный сайт, помогает в учебе. Кроме этого, можно заработать самому, выставляя готовые учебные материалы на продажу здесь. Рейтинги и отзывы на преподавателей очень помогают сориентироваться в начале нового семестра. Спасибо за такую функцию. Ставлю максимальную оценку.
Лучшая платформа для успешной сдачи сессии
Познакомился со СтудИзбой благодаря своему другу, очень нравится интерфейс, количество доступных файлов, цена, в общем, все прекрасно. Даже сам продаю какие-то свои работы.
Студизба ван лав ❤
Очень офигенный сайт для студентов. Много полезных учебных материалов. Пользуюсь студизбой с октября 2021 года. Серьёзных нареканий нет. Хотелось бы, что бы ввели подписочную модель и сделали материалы дешевле 300 рублей в рамках подписки бесплатными.
Отличный сайт
Лично меня всё устраивает - и покупка, и продажа; и цены, и возможность предпросмотра куска файла, и обилие бесплатных файлов (в подборках по авторам, читай, ВУЗам и факультетам). Есть определённые баги, но всё решаемо, да и администраторы реагируют в течение суток.
Маленький отзыв о большом помощнике!
Студизба спасает в те моменты, когда сроки горят, а работ накопилось достаточно. Довольно удобный сайт с простой навигацией и огромным количеством материалов.
Студ. Изба как крупнейший сборник работ для студентов
Тут дофига бывает всего полезного. Печально, что бывают предметы по которым даже одного бесплатного решения нет, но это скорее вопрос к студентам. В остальном всё здорово.
Спасательный островок
Если уже не успеваешь разобраться или застрял на каком-то задание поможет тебе быстро и недорого решить твою проблему.
Всё и так отлично
Всё очень удобно. Особенно круто, что есть система бонусов и можно выводить остатки денег. Очень много качественных бесплатных файлов.
Отзыв о системе "Студизба"
Отличная платформа для распространения работ, востребованных студентами. Хорошо налаженная и качественная работа сайта, огромная база заданий и аудитория.
Отличный помощник
Отличный сайт с кучей полезных файлов, позволяющий найти много методичек / учебников / отзывов о вузах и преподователях.
Отлично помогает студентам в любой момент для решения трудных и незамедлительных задач
Хотелось бы больше конкретной информации о преподавателях. А так в принципе хороший сайт, всегда им пользуюсь и ни разу не было желания прекратить. Хороший сайт для помощи студентам, удобный и приятный интерфейс. Из недостатков можно выделить только отсутствия небольшого количества файлов.
Спасибо за шикарный сайт
Великолепный сайт на котором студент за не большие деньги может найти помощь с дз, проектами курсовыми, лабораторными, а также узнать отзывы на преподавателей и бесплатно скачать пособия.
Популярные преподаватели
Добавляйте материалы
и зарабатывайте!
Продажи идут автоматически
6417
Авторов
на СтудИзбе
307
Средний доход
с одного платного файла
Обучение Подробнее