. 作如下:是單變數函數, 其形式待定. 設, 再令如下:滿足,設處的外法向量為,..." />

如何使一個R^3中徑向函數的Laplacian儘可能地小?

現在問題歸結為: 設fin C^2[0,1] 滿足f(0)=1, f(1)=0,且前兩階導數在0和1處為0. 證明存在xi in[0,1]使|2f

該問題來源於: 在mathbb{R}^3中, 作球A, 再作同心大球B, 試構造徑向函數g使gA內取1, B外取0, 且|Delta g|Bsetminus A中的最大值盡量小.

構造g是為了解決: 設Omegamathbb{R}^3中邊界光滑的有界開區域, 設varphipartialOmega上取1, 在Omega外調和且在無窮遠處為0(即孤立導體電勢). 設m{x}_0inpartialOmega, 且有含於Omega半徑為
ho的內球(切於m{x}_0), 則|
ablavarphi(m{x}_0)|lefrac{C}{
ho}C是一區域無關的常數. 方法是構造函數來估計, 構造得越好, C就能減小. 其實題主是想知道導體表面尖鈍和電荷密度的關係

============================

補充前面問題的聯繫.

為方便設內球球心為原點. 作B_1=B_{
ho_1}(m{0})supsetoverline{Omega}, B_2=B_{
ho_2}(m{0}), 
ho_2><br />
ho_1. 作g:overline{B_2}setminusOmega
ightarrowmathbb{R}如下:

g(m{x})=egin{cases}
1,    m{x}in B_1setminusOmega\
f_{
ho_2-
ho_1}(|x|-
ho_1),    m{x}in B_2setminus B_1\
0,    m{x}in overline{B_2}setminus B_1
end{cases}

f_a(t)是單變數函數, 其形式待定. 設|Delta g|le A_{
ho_2-
ho_1}, 再令z:overline{B_2}setminusOmega
ightarrowmathbb{R}如下:

z(m{x})=g(m{x})-frac{A_{
ho_2-
ho_1}}{6}left[(
ho_2^2-
ho_0^2)frac{
ho_0^{-1}-|m{x}|^{-1}}{
ho_0^{-1}-
ho_2^{-1}}-|m{x}|^2+
ho_0^2
ight]

滿足Delta zge0,zle g,設m{x}_0partialOmega的外法向量為m{n}, 則

-|
ablavarphi(m{x}_0)|=frac{partialvarphi(m{x}_0)}{partialm{n}}gefrac{partial z(m{x}_0)}{partialm{n}}

|
ablavarphi(m{x}_0)|lefrac{A_{
ho_2-
ho_1}}{6}left(frac{
ho_2^2}{
ho}+
ho_2-2
ho
ight)

f_a(t)=frac{6}{a^5}(a-t)^5-frac{15}{a^4}(a-t)^4+frac{10}{a^3}(a-t)^3

A_{
ho_2-
ho_1}
=frac{400+280sqrt{7}}{81(
ho_2-
ho_1)^2}

代入有

|
ablavarphi(m{x}_0)|lefrac{200+140sqrt{7}}{243(
ho_2-
ho_1)^2}left(frac{
ho_2^2}{
ho}+
ho_2-2
ho
ight)


ho_2
ightarrow+infty|
ablavarphi(m{x}_0)|lefrac{200+140sqrt{7}}{243
ho_0}<frac{2.35}{
ho}

可見常數依賴於f_a(t)的構造, 後用24次多項式降到了1.42, 但該方法似無法降到4/3以下. 但已有答主用其他方法得到了常數1!


@岸島柔奈子 現在已經給出了對原始問題的正確的證明。我在這裡只是把他的證明轉寫成「物理」一點的語言,並稍微推廣一下。

命題表述為

考慮任意一導體,電勢為1,無窮遠處電勢為0,導體內可容下一半徑為
ho的球體,並於導體表面相切於mathbf x_0點。那麼mathbf x_0處的電場強度小於1/
ho

記該導體在空間中產生的電勢分布為w(mathbf r)。記該球球心位置為mathbf y_0,考慮這樣一個導體球在電勢為1時在空間中產生的電勢分布為v(mathbf r)=
ho/|mathbf r - mathbf y_0|。考慮將原導體上的電荷分布固定,球形導體上電荷分布固定並反號,則由疊加原理,它們共同在導體外空間中產生的電勢分布為w(mathbf r)-v(mathbf r)。注意到這個函數在原導體表面處處大於等於零,以及在無窮遠處趨於零,我們將證明這個函數在導體外全空間也處處大於等於零。

若存在一點小於零,則存在其他地方流入該點的電場線,並必要流向電勢更低點,但無窮遠處電勢為零,故在有限遠處必存在一個電勢最低點,與自由空間中電勢無極值矛盾。故該電勢在導體外的全空間大於等於零。

那麼考慮mathbf x_0處,此處電勢為零,則電場線必指嚮導體內。

故由原導體電荷分布產生的電場小於等於由新加球形導體電荷分布產生的電場,即1/
ho

證畢。

這個命題可以簡單推廣成

考慮任意一導體,電勢為1,無窮遠處電勢為0,導體可被一半徑為
ho的球體容下,並於導體表面相切於mathbf x_0點。那麼mathbf x_0處的電場強度大於1/
ho

證明完全是類似的,這裡只要證明函數w(mathbf r)-v(mathbf r)在球體外處處小於等於零即可。

可以看到,原命題是對導體表面電場強度和表面曲率之間很強的一個估計,而第二個命題的意義就不是那麼大了,僅對全凸的導體才有比較好的效果。


First version: 或許是輔助函數取得不夠好, 目前我只能證明到梯度模長不超過
ho^{-2}.

Second version: 在原來所作的輔助函數前乘上係數 
ho, 就可以得到最終想要的估計.

======================================================================

Omegasubsetmathbf R^3 為有界光滑區域, 它在 x_0inpartialOmega 處滿足內球條件, 即存在 y_0inOmega 以及 
ho>0 使得 B_
ho(y_0)subsetOmega,~overline{B_
ho(y_0)}cappartialOmega={x_0}.mathbf R^3setminusOmega 上的調和函數 u 滿足u|_{partialOmega}=1,~lim_{x	oinfty}u(x)=0. 則有 |
abla u(x_0)|leq
ho^{-1}.

證明. 由橢圓正則性理論, 
abla u(x_0) 有意義. 讓 v(x)=1-u(x),~forall xinmathbf R^3setminusOmega. 極值原理表明 vOmega 上的正函數, 並且 v|_{partialOmega}=0. 我們的目標是證明 |
abla v(x_0)|leq
ho^{-1}.

定義 w:mathbf R^3setminus{y_0}	omathbf Rw(x)=1-
hocdot|x-y_0|^{-1}, 容易驗證

Delta w(x)=0,~forall xinmathbf R^3setminus{y_0}.

w(x_0)=v(x_0)=0,~w(x)>0, ~forall xin(mathbf R^3setminusoverlineOmega)cup(partialOmegasetminus{x_0}).

由 Kelvin 變換以及比較原理可知, |v(x)|leq w(x),~forall xinmathbf R^3setminusOmega.


u=(x_0-y_0)/
ho, 注意到 v|_{partialOmega}=0, 從而有 |
abla v(x_0)|leqlim_{t	o0+}frac{w(x_0+t
u)}{t}=
ho^{-1}.


另一位答主好像已經解決了原問題,我只算了一下題主歸結到的那個問題

首先f二階導在端點的限制是沒有意義的,因為二階導可以做一個很陡的跳變。這時候Laplacian最小的解是當x小於某個值時取某個負值,x大於某個值的時候取相反數。算得的最大值下界為2/(1-2^(1/3))≈7.7,比8小一點


如果只是想知道曲率與電荷密度的關係很簡單吧,我猜電荷密度反比與這一點平均曲率


推薦閱讀:

如何理解焓和焓變?
物理學和信息學的信息概念有不同嗎?怎麼理解?
如何理解量子霍爾效應?
如何評價《萬有引力新論》(謝兆帥著)這本書?
美國大學的物理系中,在「量子信息與量子計算」領域成果最突出、實力最強的是哪所學校?

TAG:數學 | 物理學 | 數學分析 | 實變與泛函 | 偏微分方程 |