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

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

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

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

При удалении от точки M0 это полебудет убывать, поэтому мы должны потребовать выполнение ряд условийаналогичных в рассмотренном случае при отсутствии проводящейзаземленной поверхности.Тогда мы можем записать математическую постановку нашей задачи вследующем виде:Δu = − 4πqδ(x − x0)δ(y − y0)δ(z − z0), (x, y) ∈ D, -∞ < z < ∞,u ∂Ω = 0,u → 0, uz → 0; z → ± ∞ .Будем искать решение u(x, y, z) , которое допускает вместе со своимипроизводными преобразование Фурье по переменной z .Запишем преобразование Фурье по переменной z :u^(x, y, μ) =∞12π ∫u(x, y, z)e −iμz dz .−∞Применим преобразование Фурье к левой части рассматриваемогоуравнения:1∞2π ∫Δu ⋅ e −iμz dz =−∞= Δ21∞1(Δ2u + uzz)e∫2π−iμz−∞∞2π ∫−∞ue −iμz dz +1∞2π ∫uzze −iμz dz .−∞215dz =∞12π ∫ue −iμz dz = u^(x, y, μ) , а Δ2 - оператор Лапласа в поперечном−∞сечении.Теперь вычислим последний интеграл полученного выражения,используя условия на бесконечности для функции u(x, y, z) .

Интегрироватьбудем два раза по частям:∞12π ∫uzze−iμzdz =−∞z=∞12πuze−iμz+iμz=−∞∞2π ∫uze −iμz dz =−∞=0=iμ2πue−iμzz=∞z=−∞−μ∞12π ∫2ue −iμz dz .−∞=0=u^(x, y, μ)Таким образом, применяя преобразование Фурье к нашему уравнению играничному условию, для Фурье - образа мы получаем следующую задачу:Δ2u^ − μ 2u^ = − 2 2πqe −iμz0δ(x − x0)δ(y − y0), (x, y) ∈ D,u^L= 0,где L - граница области D .Решение полученной задачи удобно искать в виде разложения в рядФурье по системе нормированных на единицу собственных функций vn(x, y)задачи Штурма - Лиувилля :Δvn + λn2vn = 0, (x, y) ∈ D,vnL=0216в поперечном сечении D :∞^ nds .Cn(μ) = uv∫u^(x, y, μ) =C μ v x, y ,∑ n( ) n( )n=1DУмножая уравнение Δvn + λn2vn = 0 на vn(x, y) , интегрируя по области Dи применяя вторую формулу Грина , получим:^ ndS − μ 2 uv^ ndS = − 2 2πqe −iμz0vn(x0, y0),Δ2uv∫∫DDгде:^ ndS = uΔ^ 2vndS = − λn2 uv^ ndS = − λn2Cn(μ) ,Δ2uv∫∫∫D∫DD^ ndS = Cn(μ) .uvDCледовательно коэффициенты Cn(μ) удовлетворяют алгебраическомууравнению:−λn2Cn(μ) − μ 2Cn(μ) = − 2 2πqe −iμz0vn(x0, y0) .Значит:Cn(μ) =2 2πqe −iμz0vn(x0, y0)λn2 + μ 2217.Тогда решение примет вид:u(x, y, z) =1∞u^(x, y, μ)e iμz dμ =2π ∫−∞∞= 2q∞∫ ∑−∞e −iμz0vn(x, y)vn(x0, y0)λn2 + μ 2n=1e iμz dμ =∞e iμ(z − z0)= 2qv x , y v x, ydμ .∑ n( 0 0) n( ) ∫ λn2 + μ 2n=1∞−∞Интеграл в последнем выражении можно вычислить с помощью вычетов,применяя лемму Жордана и замыкая контур в верхней полуплоскости при z − z0 > 0 ив нижней полуплоскости при z − z0 < 0 .Таким образом , мы получим:u(x, y, z) = 2πq∞∑n=1vn(x0, y0)vn(x, y)λne −λnz − z0, где vn(x, y) - собственныйфункции задачи Штурма - Лиувилля в поперечном сечении D .Пример 2 .

Найти потенциал поля точечного заряда величины q ,помещенного внутри двугранного угла величины α , α ∈ (0; 2π) . Грани углапредставляют собой проводящие заземленные плоскости.Решение:Будем строить решение рассматриваемой задачи с помощьюпреобразования Фурье. Тогда математическая постановка нашей задачибудет иметь следующий вид:218Δ3u = −uψ =0где:4πqr0=uδ(r − r0)δ(ψ − ψ0)δ(z − z0),ψ =a= 0,0 < r, r0 < + ∞, 0 < ψ, ψ0 < a, − ∞ < z, z0 < + ∞,1 ∂∂u1 ∂2 u∂2 uΔ3u =r+ 2+ 2 .2()r ∂r∂rr ∂ψ∂zЗапишем преобразование Фурье по переменной z :u^ =1∞2π ∫u(r, ψ, z)e −iμz dz .−∞В пространстве Фурье - образов полученное уравнение приметследующий вид:2 2π1 ∂∂ u^1 ∂ 2 u^2^−iμz 0r+ 2−μu=−δr−rδψ−ψe.()()002r ∂r ( ∂r ) r ∂ψr0Решение полученного уравнения будем искать в виде разложения в рядФурье по системе функций c граничными условиями рассматриваемойзадачи:{Φn} =πnsin ψ , n = 1,2, .

. .{α }Тогда:∞πnu^(r, ψ, μ) =Rn(r)sin ψ ,∑αn=1αRn(r) =2 ^πnu(r, ψ, μ)sin ψdψ .α∫α0219Умножим уравнение , полученное в пространстве Фурье - образов , наπnsin ψ и проинтегрируем по переменной ψ от 0 до α :ααα001 ∂∂ ^πnπnu(r, ψ, μ)sin ψdψ − μ 2 u^(r, ψ, μ)sin ψdψ+r∫ααr ∂r∂r ∫α22 2π −iμz1 ∂ u^(r, ψ, μ)πnπn0δ r − r+ 2sinψdψ=−esinψ .(0)r ∫∂ψ 2αr0α 00Отсюда мы получаем уравнение для радиальной части:dRn1 ∂1 πnr− 2Rn(r) − μ 2 Rn(r) =r ∂r ( dr ) r ( α )2=−4 2παr0e −iμz0sinπnψ0δ(r − r0) .αЕсли мы учтем смысл обобщенных функций - r 2δ(r − r0) = r02δ(r − r0) иумножим полученное уравнение на r 2 , то получим:πnr 2 R″n(r) + rR′n(r) −+ μ 2r 2 Rn(r) =[( α )]2=−4 2παr0e −iμz0sinπnψ0δ(r − r0) .αCледовательно, наша задача свелась к построению Функции Грина уравненияБесселя чисто мнимого аргумента.Потребуем выполнение дополнительных условий:Rn(0) < ∞, Rn(r) < ∞ , при r → ∞ .220Будем искать решение поставленной задачи с заявленнымидополнительными условиями в следующем виде:Rn(r) =C1I πnα (μr), r < r0,C2K πnα (μr), r > r0 .Функции I πnα (μr) и K πnα (μr) представляют собой решения однородногоπn+ μ 2r 2 Rn(r) =уравнения r 2 R″n(r) + rR′n(r) −[( α )]2=−4 2παr0eπnsin ψ0δ(r − r0) , ограниченные при r = 0 и при r → ∞ .α−iμz 0Потребуем выполнение следующих условий сопряжения при r = r0 :Rn(r0 + 0) − Rn(r0 − 0) = 0,R′n(r0 + 0) − R′n(r0 − 0) = −4 2παr0sinπn−iμz 0ψe.0αПодставляя Rn(r) в полученную систему , получаем выражения дляопределения коэффициентов C1 и C2 :C2K πnα (μr0) − C1I πnα (μr0) = 0,μC2K′ πnα (μr0) − μC1I′ πnα (μr0) = −4 2παr0sinπn−iμz 0ψe.0αИ находим:C2 = C1I παn (μr0)K παn (μr0),C1μ[I πnα (μr0)K′ πnα (μr0) − I′ πnα (μr0)K πnα (μr0)] ==−4 2παr0sinπn−iμz 0 πnψeK α (μr0) .0α221Необходимо учесть, что определитель Вронского функций Инфельда иМакдональда равен:W[I πnα (μr0), K πnα (μr0)] == [I πnα (μr0)K′ πnα (μr0) − I′ πnα (μr0)K πnα (μr0)] = −1,μr0и получаем наши коэффициенты:C1 =C2 =4 2πα4 2παsinπn−iμz 0 πnψeK α (μr0),0αsinπn−iμz 0 πnψeI α (μr0) .0αСледовательно:Rn(r) =4 2παI πnα (μr)K πnα (μr0), r < r0,πn−iμz 0sin ψ0eαI πnα (μr0)K πnα (μr), r > r0 .Теперь подставляем полученную радиальную функцию в исходноеуравнение и получаем:u^(r, ψ, μ) =∞=4 2π −iμze 0α∑I πnα (μr)K πnα (μr0)sinπnπnψ0sin ψ, r < r0αα∑I πnα (μr0)K πnα (μr)sinπnπnψ0sin ψ, r > r0 .ααn=1∞n=1Проведя обратное преобразование Фурье , мы получаем решениеисходной задачи:u(r, ψ, z) =1∞u^(r, ψ, μ)e iμz dμ .2π ∫−∞2225.Внутренние и внешние двумерные задачи и методы их решения.Пусть D - область на плоскости, ограниченная достаточно гладкойзамкнутой кривой L .

Кривая L - в нашем случае будет являтьсяповерхностью Ляпунова.Рассмотрим внутреннюю задачу Дирихле:Δu = − F(M ), M ∈ D,{u L = f (P), P ∈ L .Определение: классическим решением рассматриваемой задачи , будемназывать функцию u(M ) , дважды непрерывно дифференцируемую вобласти D , непрерывную в замкнутой области D , удовлетворяющую вклассическом смысле уравнению Δu = − F(M ) в области D и граничномуусловию u L = f (P) .При F ∈ L2(D) ∩ C (1)(D) и f ∈ C(L) наша задача имеет единственноеклассическое решение. Для его построения можно повторить все те жерассуждения, что и в трехмерном случае, взяв во второй формуле Гринафундаментальное решение оператора Лапласа на плоскости:11G(Q, M) =ln+ v , v - гармоническая функция.2π rQMu(M ) =∂u(P)∂G(P, M )G(P, M )− u(P)dlP−)∮(∂nP∂nPL− G(Q, M)Δu(Q)dSQ =∫D=∮(G(P, M )L∂u(P)∂G(P, M )− f (P)dlp + G(Q, M)F(Q)dSQ .)∫∂nP∂nPD223Как и в трехмерном случае, в полученном выражении можно убрать∂u(P)слагаемое, содержащее неизвестное значениена границе L , если∂nPвоспользоваться произвольностью гармонического слагаемого v ипотребовать выполнения условия:G(P, M ) = 0, ∀P ∈ L .Тогда:u(M ) = −∮f (P)L∂G(P, M )dlP + G(Q, M)F(Q)dSQ .∫∂nPDТеперь введем следующее определение: функцией Грина внутренней задачиДирихле для оператора Лапласа в двумерном случае будем называть функцию :G(Q, M) =11ln+ v(Q, M), Q ∈ D, M ∈ D,2π rQMудовлетворяющую условиям:1) v(Q, M) - гармоническая функция координат точки Q ∈ D ,непрерывная на D для каждой точки M ∈ D ;2) G(P, M )P∈L= 0 для каждой точки M ∈ D .Для того, чтобы построить решение рассматриваемой двумернойзадачи Дирихле достаточно найти такую функцию v(Q, M) :ΔQv = 0, Q ∈ D,v=−L11ln, PrPM2π224∈ L.Пусть De - дополнение некоторой ограниченной замкнутой области D сгладкой замкнутой границей L до всей плоскости R 2 .

Также, как и втрехмерном случае, для того, чтобы краевая задача для уравненияПуассона или Лапласа в области De имела единственное решение, следуетпотребовать регулярности решения на бесконечности.Определение: функция u(M ) называется регулярной на бесконечности вдвумерном случае, если она ограничена при r → ∞, где r =x2 + y2 .Внешняя задача Дирихле для уравнения Пуассона в двумерном случаеставится следующим образом:Δu = − F(M ), M ∈ De,uL= f (P), P ∈ L,u < ∞.Если функция F является финитной и непрерывно дифференцируемой, афункция f является непрерывной, то рассматриваемая задача имеетединственное, классическое решение.замечание: так как в двумерном случае от функций требуется толькоограниченность на бесконечности, формулы Грина во внешних областяхостаются справедливыми лишь для регулярных функций, гармонических вненекоторой ограниченной области .Так как функция F(M ) является финитной , то для решениярассматриваемой задачи справедливы формулы Грина.Выберем систему координат таким образом, чтобы начало координат Oнаходилось строго внутри области D .225Применим третью формулу Грина для решения рассматриваемой задачи:Ω(M )u(M ) − 2πu∞ =∂u(P)1∂1ln− u(P)lndlP−}∮ { ∂nPrPM∂nP rPML−∫Δu(Q)lnDe1rQMdSQ,где :Ω(M ) = 2π , если M ∈ De , Ω(M ) = π , если M ∈ L , Ω(M ) = 0 , еслиM ∉ De .Перепишем полученную формулу, взяв в качестве точки M началокоординат, в следующем виде:−2πu∞ =∂u(P)1∂1ln− u(P)lndlp−∮ { ∂nPrOP∂nP rOP }L−∫Δu(Q)lnDe1rOQdSQ, так как точка O не принадлежит области De .Пусть M - произвольная точка области De .

Тогда:u(M ) =12π ∮L−u(P)∂u(P)11ln− ln−∂nP ( rPMrOP )∂11ln− lndlP−rOP )∂nP ( rPM226−111Δu(Q) ln− lndSQ .∫2πrOQ )( rQMDe∂u(P)Значениена границе области неизвестно. Применим∂nPстандартный прием для того, чтобы убрать слагаемое, содержащее этонеизвестное значение.Пусть v2 - произвольная гармоническая в области De и регулярная набесконечности функция.Для решения u(M ) рассматриваемой задачи и функции v2 справедливавторая формула Грина в области De :0=∂v (P)∂u(P)v2(P) − u(P) 2dlP − Δu(Q)v2(Q)dSQ .}∮ { ∂nP∫∂nPLDeСкладывая оба равенства , получим:u(M ) =∂u(P)∂G(P, M )−u(P)G(P, M ) dlP−}∮ { ∂nP∂nPL−∫Δu(Q)G(Q, M)dSQ,De111где: G(Q, M) =ln− ln+ v2 .2π ( rQMrOQ )замечание: при построении функции G(Q, M) вместо точки O можновыбрать любую точку строго внутри области D .Определение: Функцией Грина внешней задачи Дирихле для оператора Лапласа вдвумерном случае будем называть функцию:227G(Q, M) =11ln+ v(Q, M), Q ∈ De, M ∈ De , удовлетворяющую2π rQMусловиям:1) v(Q, M) - гармоническая функция координат точки Q ∈ De ,непрерывная в De для каждой точки M ∈ De , имеющая логарифмическуюособенность на бесконечности;2) G(P, M )P∈L= 0 для каждой точки M ∈ De ;3) G(Q, M) регулярна на бесконечности.Следовательно, функция Грина G(Q, M) является решением следующейзадачи:ΔQG(Q, M) = − δ(Q, M), Q ∈ De, M ∈ De,G(P, M )P∈L= 0, P ∈ L,G(Q, M) < ∞ .Так как функция Грина регулярна на бесконечности , то можно показать, чтоона симметрична относительно перестановки точки наблюдения Q и точкиисточника M .228Теперь перейдем к методам решения двухмерных задач:а) Метод электростатических изображений.Как и в рассмотренном нами трехмерном случае, в ряде областей прирешении двухмерных задач, удобно использовать метод электростатическихизображений.Рассмотрим следующий пример1:Найти функцию Грина задачи Дирихле вне круга радиуса a .Решение:Функция Грина является решением следующей задачи:ΔMG(M, M0) = − δ(M, M0), r > a, r0 > a,Gr=a= 0,G < ∞ при r → ∞ .ΔM - оператор Лапласа, где производные берутся по координатам точкиM , а r и r0 - полярные радиусы точек M и M0 .Будем искать функцию G(M, M0) как:G(M, M0) =11ln+ v(M, M0) ,2π rMM0где:ΔM v(M, M0) = 0, r > a, r0 > a,Gr=a=−11lnrPM02π229.a2Пусть M1, ψ0 - точка , сопряженная M0(r0, ψ0) относительно( r0)окружности радиуса a .Функция v = −1a 1является гармонической вне круга радиусаln2π r0 rMM1a и удовлетворяет нашей задаче.Тогда решение нашей исходной задачи примет следующий вид:G(M, M0) =111a 1ln−ln.2π rMM0 2π r0 rMM1пример2: Для любой непрерывной функции f (ψ) построить решение задачиДирихле для уравнения Лапласа в круге в интегральной форме:Δu = 0, r < a, ψ ∈ [0,2π],{u r=a = f (ψ) .Решение:Найдем решение поставленной задачи с помощью формулы :∂G(P, M )u(M ) = − f (P)dlP + G(Q, M)F(Q)dSQ , в которой∮∫∂nPLDфункция Грина G(M, M0) определяется выражением :G(M, M0) =11a 1ln− ln,rMM02π( r0 rMM1 )230r 2 + r02 − 2rr0cos(ψ − ψ0) ,rMM0 =rMM1 =r 2 + r12 − 2rr1cos(ψ − ψ0) ,a2.r1 =r0Подсчитаем производную:∂G∂n=r=aa−lnr01 ∂ln2π ∂r1r 2 + r02 − 2rr0cos(ψ − ψ0)1=r 2 + r12 − 2rr1cos(ψ − ψ0)r=a1 ∂=ln2π ∂r−ln1r 2 + r02 − 2rr0cos(ψ − ψ0)−1r 2r02a2=+ a 2 − 2rr0cos(ψ − ψ0)r=aa 2 − r021=−⋅.2πa a 2 + r02 − 2ar0cos(ψ − ψ0)231−Подставляя найденное выражение в нашу формулу получаем, чторешение в любой точке M0(r0, ψ0) :2πa 2 − r021u(r0, ψ0) =f (ψ)dψ2π ∫ a 2 + r02 − 2ar0cos(ψ − ψ0)0Полученная в задача формула называется интеграл Пуассона .пример3: найти потенциал поля бесконечной заряженной нити с линейнойπплотностью заряда q , помещенной внутри двугранного угла величины параллельноnребру этого угла.

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

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

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

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