Т.А. Шапошникова - Курс лекций по уравнениям в частных производных (1127861), страница 8
Текст из файла (страница 8)
Лекция 17.17.1. Функция Грина. Задача Дирихле для уравнения Лапласа.Теорема 17.1 (Лиувилль). Пусть u(x) - гармоническая в Rn , неотрицательная функция. Тогда u = const.Доказательство: Зафиксируем точку x0 и шар QxR0 радиуса R с центром в нашей точке.Поскольку производная гармонической функции – также функция гармоническая, по теореме о среднем имеем:ZZ∂u1∂u1(x0 ) = x0(x) dx = x0u(x)νj (x)dS∂xj|QR | QxR0 ∂xj|QR | SRx0Мы использовали формулу Стокса, чтобы перейти к интегрированию по границе шара. А теперь используемеще одну теорему о среднем, на этот раз из курса математического анализа:Z|S x0 |1= x0 νj (ex)u(x)dS = Rx0 νj (ex)|u(x0 )|x|R ||QR |SR0x0Здесь |SR| = wn Rn−1 ,а |QxR0 | =wn nn R .Таким образом,|∂un(x0 )| 6 |u(x0 )| → 0 R → ∞∂xjR∂u(x0 )| = 0 ∀ j = 1, .
. . , n.Ничего не мешает нам выбрать радиус шара сколь угодно большим, а значит, | ∂xjСледовательно, u = const, что и требовалось.Задача: Пусть u(x) - гармоническая в Rn и u(x) > −C(1 + |x|m ), где c, m > 0 u(x). Показать, что в этомслучае u(x) есть полином степени не выше [m].Рассмотрим задачу Дирихле для оператора Лапласа:(∆u(x) = 0 x ∈ Ωu(x) = ϕ(x) x ∈ ∂Ω, f ∈ C(∂Ω)TПод классическим решением понимается решение из класса C 2 (Ω) C(Ω).Для такого решения ранее была получена формулаZZZ∂E(|x − x0 |)∂uu(x0 ) =E(|x − x0 |)∆u(x) dx +u(x)dS −E(|x − x0 |) (x) dS∂ν∂νΩ∂Ω∂ΩПусть существует функция g(x, x0 ) со следующими свойствами:g(x, x0 ) ∈ Cx2 (Ω) , ∆x g(x, x0 ) = 0 ∀ x0 ∈ Ω и при этом g(x, x0 )|x∈∂Ω = −E(|x − x0 |)|x∈∂Ω , ∀ x0 ∈ ΩЗапишем вторую формулу Грина для u(x), g(x, x0 ):ZZ∂g∂u(u∆g − g∆u) dx =u−g∂ν∂νΩ∂ΩПоскольку функции u и g – гармонические в области, то левая часть формулы обращается в ноль.
Теперьприбавим к левой и правой частям уже упоминавшееся соотношение:ZZ∂E(|x − x0 |)∂udS −E(|x − x0 |) (x) dSu(x0 ) =u(x)∂ν∂ν∂Ω∂ΩИ поскольку g(x, x0 ) + E(|x − x0 |) = 0 на ∂Ω, получаем:Z∂(E(|x − x0 |) + g(x, x0 ))u(x0 ) =u(x)dS∂ν∂ΩОпределение: Функцией Грина называется функция G(x, x0 ) = E(|x − x0 |) + g(x, x0 ), где функция g(x, x0 )введена ранее. Точка x0 называется полюсом.Итак, если u ∈ C 2 (Ω) – решение задачи Дирихле, то:Z∂G(x, x0 )u(x0 ) = u(x0 ) =u(x)dS∂ν∂ΩЛемма 17.2. ∀ x1 , x0 ∈ ΩG(x1 , x0 ) = G(x0 , x1 )37Доказательство: Исходя из определения функции Грина, ∆x G(x, x0 ) = 0 при 6= x0 .
Вырежем вокруг точекx и x0 шарики маленького радиуса ε, и то что осталось обозначим за Ωε :Ωε = Ω(Qxε 0 ∪ Qxε 1 )Введем функции u(x) = G(x, x0 ), и v(x) = G(x, x1 ). В области Ωε они гармонические, поэтому по второй формулеГрина:ZZ ∂v∂u −vdS0=(u∆v − v∆u) dx =u∂ν∂νΩε∂ΩεТак как функция g определялась условиемg(x, x0 )x ∈ ∂Ω= −E(|x − x0 |)x ∈ ∂Ω, ∀ x0 ∈ Ω,то G(x, x0 )|∂Ω = 0 и u(x)|∂Ω = v(x)|∂Ω = 0. ПоэтомуZZ∂G(x, x1 )∂G(x, x1 )dS +G(x, x0 )dS0=G(x, x0 )x0x∂ν∂νSε 1SεПерейдем к пределу при ε → 0, воспользовавшись тем, что∂E(|x − x1 |)∂E1=−(r) = −∂ν∂rwn rn−1Получаем,чтоZxSε 1G(x, x0 )∂G(x, x1 )1dS = −∂νwn εn−1ZxSε 1G(x, x0 ) dS → −G(x1 , x0 )Аналогично переходя к пределу в другом слагаемом, имеем0 = −G(x1 , x0 ) + G(x0 , x1 )что и требовалось доказать.18. Лекция 1818.1. Интеграл Пуассонаu(x) ∈ C 2 (Ω), Ω - ограниченная область в Rn , ∂Ω ∈ C 1 .Будем решать следующую задачу (1):(∆u = f (x) в Ωu|∂Ω = ϕ(x)ZZ∂Gu(x0 ) =f (x)G(x, x0 )dx +ϕ(x)(x, x0 )ds∂νxΩ∂ΩG(x, x0 ) - функция Грина, G не более чем единственна.ρΩ = Q0R ⊂ Rn , G(x, x0 ) = E(|x − x0 |) − E R|x − x∗ |Если(∆u = f (x), x ∈ Ω,,то u(x0 ) =R0SRu ∈ C 2 (Q0R )u|SR0 = ϕ(x)0)ϕ(x) ∂G(x,xdsx .∂νx∂G(x, x0 )∂|x − x0 |ρρ∂|x − x∗ ||x∈SR0 ,x0 ∈Q0R = E ′ (r)− E ′ ( r1 ),∂ν∂νRR∂νгде r = |x − x0 |, r1 = |x − x∗ |.nPxj −x0j∂|x−x0 |=∂ν|x−x0 | = cos γ, где γ - угол между ν и x − x0 .
Аналогичноj=1222∂|x−x∗ |∂ν= cos β.−ρПо теореме косинусов ρ2 = R2 + r2 − 2Rr cos γ ⇒ cos γ = R +r.2RrАналогично2 24R2 + Rρ2r − ρr2R2 + r12 − ρ211 ρ2 + r 2 − R 2cos β ===2Rr1ρ2r2R rRρ382−nТ.к. E(r) = − ωnr(n−2) , то E ′ (r) = ωn r1n−1 .Тогда∂G(x, x0 )ρR 2 − ρ2R2 − |x0 |21|x∈SR0 ,x0 ∈Q0R = E ′ (r)(cos γ − cos β) = E ′ (r)=∂νRRrωn R|x − x0 |nСледовательно,R2 − |x0 |2u(x0 ) =ωn RZr2 = R2 + ρ2 − 2Rρ cos α, тогдаu(x0 ) =ϕ(x)dsx|x − x0 |n0SRR2 − |x0 |2ωn RZ0SR- интеграл Пуассонаϕ(x)dsx(R2 + |x0 |2 − 2R|x0 | cos α)n/20Пусть ϕ ∈ C(SR).
Тогда функция u(x), задаваемая интегралом Пуассона, есть решение задачи (1). Необходимопроверить 2 условия:R0)1) u(x0 ) = ∂Ω ϕ(x) ∂G(x,xdsx - гармоническая по x0 , x0 ∈ Q0R .∂νx02) ∀ x̂ ∈ SR∃ lim u(x0 ) = ϕ(x̂)x0 →x̂Докажем это.1) При x 6= x0 ∆x G(x, x0 ) = 0 ⇒R ∆x0 G(x0 , x) = 0 ⇒ ∆x0 G(x, x0 ) = 0 в силу симметричности G.x0 ∈ Ω ⊂⊂ Q0R ⇒ ∆x0 u(x0 ) = ∂Ω ϕ(x) ∂ν∂x ∆x0 G(x, x0 )dsx = 02) В силу решения задачи Дирихле задача(∆u = 0, x ∈ Q0Ru|SR0 = 1имеет единственное решение u ≡ 1, тогдаR2 − |x0 |21=ωn RZdsxrnR2 − |x0 |2ωn RZϕ(x̂)dsxrnϕ(x̂) =0SR0SR R2 − |x |2 Z ϕ(x) − ϕ(x̂) R2 − |x |2 Z |ϕ(x) − ϕ(x̂)|00|u(x0 ) − ϕ(x̂)| = dsdsx =6x ωn R00rnωn RrnSRSRZZR2 − |x0 |2|ϕ(x) − ϕ(x̂)|R2 − |x0 |2|ϕ(x) − ϕ(x̂)|=ds+dsxxn0ωn RrωRrnnσδ(ε)SR rσδ(ε)Первое слагаемое обозначим за I1 , второе за I2 .0Здесь мы пользовались тем, что ∀ ε ∃δ > 0 : x ∈ σδ = Qx̂δ ∩ SR⇒ |ϕ(x) − ϕ(x̂)| < 2ε .I1 6ε R2 − |x0 |22 ωn RZ0SRdsxε=nr20Что можно сказать об I2 ? При x ∈ SRr σε имеем |x − x0 | > a > 0, как только |x0 − x| < δ1 .
ТогдаI2 6 2 max|ϕ(x)|a−n0SRωn Rn−1 2(R − |x0 |2 ) = c1 (R2 − |x0 |2 )ωn RС другой стороны, ∃δe > 0 : |x0 − x̂| < δe ⇒ R2 − |x0 |2 <εc1 .Тогда I2 < 2ε , а это и требовалось.18.2. Неравенство ХарнакаТеорема 18.1 (Неравенство Харнака). Пусть u(x) — гармоническая в Q0R функция, непрерывная вплотьдо границы; u(x) > 0. Тогда ∀ x ∈ Q0RRn−2 (R − ρ)Rn−2 (R + ρ)6 u(x0 ) 6,n−1(R + ρ))(R − ρ)n−1 )39где ρ = |x0 |.ZR2 − |x0 |2u(x)dsx ; R − ρ 6 r 6 R + ρ0ωn RrnSRZZR2 − |x0 |21R2 − |x0 |21u(x)ds6u(x)6u(x)dsxx0ωn R(R + ρ)n SR0ωn R(R − ρ)n SR0u(x0 ) =НоR0SRu(x)dsx = u(0)ωn Rn−1 , откуда и получаем нужное нам утверждение.
18.3. Обратная теорема о среднемПусть u ∈ C(Ω) и для любого шара Qxr 0 ⊂ Ω для u справедлива теорема о среднем по сфере. Тогда u(x) гармоническая в Ω функция.Доказательство.Задача(∆v = 0, x ∈ QxR0 , Qxr 0 ⊂ Ωv|SRx0 = u|SRx0имеет единственное решение, которое задается интегралом Пуассона.Тогда имеем v(x) ≡ u(x) в QxR0 , что следует из равенства u(x) и v(x) на ∂QxR0 , выполнения теоремы о среднем ипринципа максимума.18.4.
Теорема об устранимой особенностиПусть u(x) - гармоническая в Ω r {x0 } функция, m(ρ) = sup |u(x)|.xΩrQρ 0Пусть m(ρ) 6 a(ρ)E(ρ), где a(ρ) → 0 при ρ → 0.Тогда u(x) можно доопределить в x0 таким образом, что полученная функция будет гармонична в Ω.19. Лекция 19Доказательство теоремы об устранимой особенностиРассмотрим Qxρ10 : Qxρ10 ∈ ΩНайдем:(∆v(x) = 0, x ∈ Qxρ10v(x) = u(x), x ∈ Qxρ10Наша цель показать, что u(x) = v(x) в шаре всюду, кроме его центра x0 . Рассмотрим ω(x) = v(x) − u(x), x ∈Qxρ10 /Qxρ10 . ∆ω = 0, x ∈ Qxρ10 /Qxρ10 ω = 0 на Sρx10max |v(x)| 6 max|u(x)| = MxSρ10xQρ10Возьмем произвольное ε > 0, Qxρ10 / {x0 } , |ω(x)| 6 εE(|x − x0 |)max|ω(x)| 6 max|u(x)| + max|v(x)| 6xxxSρ10Sρ10(19.1)Sρ10∀ ε > 0∃ρ 6 ρ0 : M 6 2ε |E(ρ)|С другой стороны по условию теоремыmax|u(x)| 6 a(ρ)|E(ρ)| ∀ ε > 0∃eρ0 : ∀ ρ 6 ρe0 , a(ρ) 6xSρ0⇒ max|u(x)| 6xSρ0ε2ε|E(ρ)|2Тогда (19.1) можно продолжитьПо принципу максимума:6 ε|E(ρ)|, ∀ ρ 6 ρb|ω(x)| 6 ε|E(|x − x0 |)|, ∀ x ∈ Qxρ10 / {x0 }Следовательно, полагая u(x0 ) := v(x0 ), получаем гармоническую функцию во всем шаре4019.1.
Теория потенциалаu(x) ∈ C 2 (Ω), ∂Ω ∈ C 1Вспомним, чтоu(x0 ) =ZE(x, x0 )∆udx +ΩZu∂E(x, x0 )dS −∂ν∂ΩZE(x, x0 )∂udS∂ν∂Ω- представляется в виде трех потенциалов.Введем для каждого из потенциалов обозначения:Z1P0 (x) = ρ(ξ)dξ, n > 3 - объёмный потенциал.|x − ξ|n−2ΩP1 (x) =P2 (x) =ZΓZµ(ξ)Γσ(ξ)1dSξ - потенциал простого слоя.|x − ξ|n−2∂1dSξ - потенциал двойного слоя.∂νξ |x − ξ|n−2cn |x−ξ|1 n−2 .Здесь E|x − ξ| =Теорема 2Пусть µ(x), σ(x) ∈ L1 (Γ).
Тогда P1 (x), P2 (x) - гармонические функции в Rn /ΓДоказательствоМожно дифференцировать 2 раза P1 , P2 под знаком интеграла.∀ x ∈ Ω1 ; Ω, ∂Ω = Γ, Ω0 = Rn /Ω Ω1 ⊂⊂ ΩВозьмем Ω1 , ρ(Ω1 , Γ) > α > 0. В P1 (x) можно дифференцировать под знаком интеграла. То же самое и для Ω0Z1∆x P1 (x) = µ(ξ)∆xdSξ = 0|x − ξ|n−2ΓТеперь посмотрим на потенциал двойного слояnnXX11∂|x − ξ| k1ξk − xk k∂=−(n−2)ν=−(n−2)ν =ξ∂νξ |x − ξ|n−2|x − ξ|n−1∂ξk|x − ξ|n−1|x − ξ| ξk=1k=1= −(n − 2)ТогдаP2 (x) =ZΓcos (r, νξ )rn−1∂1σ(ξ)dSξ =∂νξ |x − ξ|n−2Zσ(ξ)ΓnX∂1ν ξ dSξ∂ξk |x − ξ|n−2 kk=1Теперь по тем же соображениям можно дифференцировать и потенциал двойного слоя.Посмотрим, как ведут себя P1 , P2 на бесконечности.
Если ξ ∈ Γ, |x| >> 1, ⇒ |x − ξ| > |x| − |ξ| >для P1 :Z2n−2M0|P1 (x)| 6|µ(ξ)|dSξ =n−2|x||x|n−2|x|2ТогдаΓСледовательно для n > 3, P1 (x) → 0, x → ∞.По тем же соображениям|P2 (x)| 6|P0 (x)| 6M1|x|n−2M2|x|n−2Теорема 19.1 (3). Пусть ρ(x) ∈ C(Ω). Тогда P0 (x) ∈ C 1 (Rn ) и P0 (x) гармоническая функция при x ∈ Rn /Ω.41x ∈ Ω ⊂⊂ Rn /Ω ⇒ ∆x P0 (x) =x ∈ Ω, Qxε - шар.RΩZP0 (x) 6ρ(ξ)∆x |x−ξ|1 n−2 dξ = 0|ρ(ξ)|Ω∩QxεZ1dξ +|x − ξ|n−2Zx ∈ Rn /Ω ⇒ |P0 (x)| 6 M∀ x0 ∈ Rn|P0 (x0 + h) − P0 (x0 )| 6Z+|ρ(ξ)|Zxx|ρ(ξ)|Ω∩Qxδ|Z1dξ 6 M1,ε|x − ξ|n−21dξ +|x0 + h − ξ|n−2Z|ρ(ξ)|Ω∩Qxδ1dξ гладкая поверхность.|x0 − ξ|n−2|ρ(ξ)| 11 dξ 6 ε−n−2n−2|x0 + h − ξ||x0 − ξ|3xZxQδ 0Ω∩Qxδ1dξ|x − ξ|n−2|ρ(ξ)|1dξ 6 max |ρ(x)||x0 − ξ|n−2Ω|ρ(ξ)|Z1dξ|x − ξ|n−21dξ+|x0 − ξ|n−211 dξ−|ρ(ξ)| n−2n−2|x0 + h − ξ||x0 − ξ|Ω/QδΩ/Qδ|ρ(ξ)|Ω∩QxδΩ/QδZZxΩ/QεQxεΩ/QxεZ|ρ(ξ)|1|ρ(ξ)|dξ 6 max |ρ(x)||x − ξ|n−2ΩΩ∩Qxε∀ ε > 0 ∃h0 > 0, ∀ |h| < h0Z1dξ = C1|x − ξ|n−2Zδrdr = C2 δ 2 <0εпри δ ≪ 131dξ|(|h| < δ/2)(|x0 + h − ξ| 6 |x0 − ξ| + |h| < 3δ/2) 6|x0 + h − ξ|n−26 max |ρ(x)|C2Ω3δ/2Z0rdr = C3 δ 2 6 ε/3 при δ ≪ 1Окончательно получаем, что при малых h|P0 (x0 + h) − P (x0 )| < ε → P0 ∈ C(Rn )Задача 19.1.∂P0∈ C(Rn )∂xk∆P0 (x) = 0, x ∈ Rn /Ω∆P0 (x) = −(n − 2)ωn ρ(x), x ∈ Ω′ , ωn = |S10 |(n = 2 : ∆P0 = −2πρ(x))Теорема 19.2 (4).
Пусть ρ ∈ C 1 ((Ω)). Тогда P0 (x) ∈ C 2 (Ω), при x ∈ Ω(−(n − 2)ωn ρ(x), n > 3∆P0 (x) =−2πρ(x), n = 2n > 3, вычислим∂P0 (x)=∂xkZΩ∂1ρ(ξ)dξ = −∂xk |x − ξ|n−242ZΩρ(ξ)∂1dξ =∂ξk |x − ξ|n−2= lim −ε→0Zρ(ξ)Ω/Ωxε∂1dξ∂ξk |x − ξ|n−2И получим то, что требуется доказать в условии. 20. Лекция 20RОбъёмный потенциал: P0 (x) = Ω ρ(ξ) |x−ξ|1 n−2 dξRПотенциал простого слоя: P1 (x) = Ω µ(ξ) |x−ξ|1 n−2 dξRПотенциал двойного слоя: P2 (x) = Ω σ(ξ) |x−ξ|1 n−2 dξ∆Pi (x) = 0, x ∈ Rn /Γ, i = 1, 2;20.1.