劉奕含， 陽 鶯,2， 覃柳術(shù)

(1.桂林電子科技大學(xué) 數(shù)學(xué)與計算科學(xué)學(xué)院，廣西 桂林 541004；2.桂林電子科技大學(xué) 廣西高校數(shù)據(jù)分析與計算重點實驗室，廣西 桂林 541004)

Poisson-Nernst-Planck(PNP)方程[1]是由Poisson方程和Nernst-Planck(NP)方程耦合而成的非線性偏微分方程，用于描述溶質(zhì)生物分子體系中移動離子的電擴散，常應(yīng)用于生物分子[1]、電化學(xué)[2]和半導(dǎo)體[3]等相關(guān)領(lǐng)域。因PNP方程的非線性性和非對稱性，其解析解只在極少數(shù)情況下存在，人們提出了各種數(shù)值方法求解PNP方程，如有限差分法[4]、有限體積法[5]及有限元法[6]。有限差分法廣泛用于求解描述生物離子通道等跨膜孔的電擴散PNP方程[7-8]，但將其應(yīng)用于表面高度不規(guī)則的PNP方程時，其精度不高；有限體積法可求解不規(guī)則表面的PNP方程，但高精度有限體積控制元的構(gòu)造較困難；有限元方法可處理邊界不規(guī)則區(qū)域問題，且其解的收斂速度依賴解的正則性。

有限元方法的基本思想是將連續(xù)的求解區(qū)域離散成有限個且按某些方式連接在一起的單元組合體，利用每一單元內(nèi)假設(shè)的近似函數(shù)分片表示求解區(qū)域待求的未知函數(shù)，將連續(xù)的無限自由度問題變成離散的有限自由度問題。為此，利用有限元方法的思想對非線性穩(wěn)態(tài)PNP方程進(jìn)行離散，提高其求解效率。

1 預(yù)備知識

設(shè)Ω?R3是一個多面體，且滿足Lipchitz連續(xù)，?Ω為其邊界，使用標(biāo)準(zhǔn)的Sobolev空間記號Ws,p(Ω)及相應(yīng)的范數(shù)和半范數(shù)[9]，當(dāng)p=2時，有

Ws,2(Ω)=Hs(Ω)，

‖·‖s,p,Ω=‖·‖Ws,p(Ω)。

定義相關(guān)的線性有限元空間：

Sh={v∈H1(Ω):v|?Ω=0∩v|e∈p1(e),

?e∈Γh(Ω)}，

SH={v∈H1(Ω):v|?Ω=0∩v|e∈p1(e),

?e∈ΓH(Ω)}，

考慮非線性穩(wěn)態(tài)的PNP方程：

(1)

方程(1)的邊界滿足齊次的Dirichlet邊界條件p1=p2=φ=0，其中pi(x)為第i種帶電荷量qi的粒子濃度，qi為第i種離子的電荷量，φ(x)為靜電勢，β=1/(κβT)為逆Boltzmann能量，ε(x)為電介質(zhì)常量。

c1(pi,v)+d1(pi,φ,v)+f1(pi,v)=

(2)

c2(φ,w)+d2(pi,w)=(F3,w),

(3)

其中，

c1(pi,v)=(pi,v)，

d1(pi,φ,v)=(βqipiφ,v)，

f1(pi,v)=(f(pi),v)，

c2(φ,w)=(εφ,w)，

2 有限元方法

(4)

(5)

引理1假設(shè)pi∈L∞(Ω)，則有

證明

證畢。

其中C為與h無關(guān)的任意常數(shù)。

(6)

f1(e,vh)=0，

(7)

因此，有

f1(e,pi-vh)-f1(e,e)。

(8)

由引理1及引理2可得，

(9)

同理

(10)

其中ε<1。由于f(pi)是Lipchitz連續(xù)函數(shù)，

(11)

(12)

由式(8)～(12)可得，

因此

記(pi)I=∏hpi，取vh=(pi)I，由插值誤差估計有

‖pi-(pi)I‖1,Ω≤Chs。

其中C為與h無關(guān)的任意常數(shù)。證畢。

3 數(shù)值實驗

求解具有光滑解的非線性穩(wěn)態(tài)PNP方程：

設(shè)Ω=[0,1]3，q1=1,q2=-1，邊界條件和右端函數(shù)依賴于解析解，其中解析解定義為：

表1為H1范數(shù)下有限元解和精確解的誤差，收斂階接近1階，表明有限元法的數(shù)值結(jié)果與理論一致。

表1 有限元解和精確解在H1范數(shù)下的誤差

4 結(jié)束語