涂馨予, 蒲志林
(四川師范大學(xué) 數(shù)學(xué)與軟件科學(xué)學(xué)院, 四川 成都 610066)
?
一類帶一般記憶核的Cahn-Hilliard方程解的能量衰減估計
涂馨予,蒲志林*
(四川師范大學(xué) 數(shù)學(xué)與軟件科學(xué)學(xué)院, 四川 成都 610066)
摘要:采用推廣的帶參數(shù)的Gronwall不等式,研究了一類帶一般記憶核的Cahn-Hilliard方程解的能量衰減估計,證明方程的解是呈指數(shù)衰減的,衰減率可被精確地算出.
關(guān)鍵詞:動態(tài)邊界條件; 記憶核; 過去歷史變量; 能量估計
Cahn-Hilliard方程是由J. W. Cahn和J. E. Hilliard[1]引入,用來描述二元合金在分解時的分離物隨溫度的變化情況,它在材料科學(xué)中有著重要意義.為了描述所謂的二元混合物的旋節(jié)分解(如金屬合金的冷卻),取Ω是R3中的有界集,并帶有充分光滑的邊界Γ,考慮粘性Cahn-Hilliard方程[1]
(1)
其中,u是混合物的相對濃度差,μ是帶有粘性的化學(xué)勢,α≥0,非線性項f(u)是一勢函數(shù)的導(dǎo)數(shù),常數(shù)α=0時,方程(1)是Cahn-Hilliard方程.對這類方程在標(biāo)準(zhǔn)邊界條件下(Neumann邊界條件、狄利克雷邊界條件、周期邊界條件)的柯西問題的討論已經(jīng)很充分[2-5].
為了描述特殊材料(如玻璃)的相分離現(xiàn)象,對方程(1)取μ滿足Neumann邊界條件,即
將通常u的Neumann邊界條件換成非線性的動態(tài)邊界條件,即u滿足
其中,β是正常數(shù),ΔΓ是邊界上的Laplace-Beltrami算子,非線性項g(u)是一充分光滑的邊界勢函數(shù)的導(dǎo)數(shù).許多作者對這類問題進(jìn)行了討論[6-13].
為了加速某些材料的旋節(jié)分離(如玻璃),考慮修改的Cahn-Hilliard方程
其中κ(s)是記憶核函數(shù),它是呈指數(shù)衰減的,即滿足條件κ′(s)+δκ(s)≤0,δ>0,這類Cahn-Hilliard方程已有大量研究[14-22].如在標(biāo)準(zhǔn)邊界條件及光滑勢條件下,3維空間上非恒溫非粘彈條件下的弱解的存在[14]、解的適定性和正則性的結(jié)果[15]已得到.在2維空間上,取α=0已經(jīng)證得了能量解在短的松弛時間的全局吸引子的存在[16].但這類問題在非線性動態(tài)邊界條件下的研究還很少.
這里考慮一類在非線性動態(tài)邊界條件下帶記憶核的Cahn-Hilliard方程[23]:
其中,Ω?R3是具有光滑邊界Γ的有界區(qū)域,α,β是正常數(shù).未知量u=u(x,t):Ω×[0,∞)→R,u是混合物的相對濃度差,u|Γ=ψ(t),u滿足非線性的動態(tài)邊界條件,μ滿足Neumann邊界條件.ΔΓ是邊界上的Laplace-Beltrami算子,η=ηt(x,s):Ω×R+×[0,∞)→R是過去歷史變量,且對任意的s>0和t>0,有ηt(s)=∫0s-Δμ(t-y)dy,顯然ηt(x,s)滿足邊界條件:ηt(x,0)=0,Ω×(0,∞).
對方程(3)給定初值條件如下:
(4)
ν(s)是方程的記憶核,它滿足以下假設(shè)[23]:
(H1)ν(s)∈C1(R+)∩L1(R+),?s∈R+;
(H2) ν(s)≥0,ν′(s)≤0,?s∈R+;
(H5) ν′(s)+ρ(s)ν(s)≤0,?s∈R+;ρ(s)>0,ρ′(s)≤0,?s∈R+;
注滿足假設(shè)(H1)~(H6)的記憶核ν(s)是存在的,如取ν(s)=e-λs,?λ>ρ0.
非線性項f、g滿足如下條件:
(K1) f,g∈C1(L2(Ω));
(K2) f,g∈C2(L2(Ω));
(K3) |f′(y)|≤cf(1+|y|2),|g′(y)|≤cg(1+|y|q),?y∈L2(Ω);
由(K5)可得[23]:
(5)
在文獻(xiàn)[23]中,取非線性記憶核ν(s)滿足ν′(s)+δν(s)≤0,δ>0是正常數(shù),證得了變分解的存在;當(dāng)α>0時,全局吸引子和指數(shù)吸引子的存在;以及α=0時,軌道吸引子的存在.受文獻(xiàn)[4]的啟發(fā),本文考慮非線性動態(tài)邊界條件下,滿足更一般條件的非線性記憶核ν(s),它滿足以下微分不等式[24]:
當(dāng)取ρ(s)=δ時,即為文獻(xiàn)[23]的情況,本文的記憶核更具一般性.再采用推廣的帶參數(shù)的Gronwall不等式,研究了帶這類記憶核的Cahn-Hilliard方程解的能量衰減性,證明當(dāng)t→∞時,方程的解的能量是指數(shù)衰減到0的,并精確地算出了衰減率.
本文在預(yù)備知識中介紹了一些定義、記號、不等式及定理;后面通過構(gòu)建泛函L,對某個正常數(shù)C滿足
1預(yù)備知識
Ω?R3是具有光滑邊界Γ的有界區(qū)域,分別記(·,·),‖·‖;(·,·)Γ,‖·‖Γ為L2(Ω)和L2(Γ)的內(nèi)積和范數(shù),定義算子A:=-Δ,A是L2(Ω)上的非負(fù)自共軛的線性算子,取A的分?jǐn)?shù)力Ar,則可定義一族插值空間如下:
當(dāng)r=0時,V0=L2(Ω),在Vr上定義范數(shù):
定義〈·〉為Ω上的空間平均,即
所以有
定義跡算子TrD為
引入空間
取Hilbert空間
M-1上的內(nèi)積定義如下
定義T=-?s是M-1上的線性算子
方程(3)的解滿足以下質(zhì)量守恒[23]:
(6)
最后,引入乘積Hilbert空間H=V1×M-1.下面介紹一些重要的不等式及引理:
2) Poincaré-Wirtinger不等式:‖u-〈u〉‖p≤C7‖u‖p.
3)W1,2(Ω)上的等價范數(shù)為:‖u‖W1,2(Ω)=‖u‖+‖u‖.
那么
注意:取?a∈(0,1),使得a=m?a-?a+1-l=(m-1)?a+1-l>0,b=1-?a-γ?a>0,定義
則
2一致能量估計
關(guān)于方程(3)在t時刻的能量E(t)為
(7)
有如下的結(jié)果:
定理 2.1假設(shè)(H1)~(H5)和(K1)~(K5)成立,則存在正常數(shù)J,ω,2 即方程的解的能量是呈指數(shù)衰減的,且衰減率為ω,其中: 證明在方程(3)的第一項的兩端同乘以μ,在Ω上積分,可以得到 利用方程(3)的第3項有 利用格林公式可得 由方程(3)的第2項可得 即 (8) 又由方程(3)的第4項及(6)式可得 將以上所求的項帶入到方程(8)可以得到 由(H2)和(K4)可得 (9) 故 (10) 構(gòu)建泛函: 利用方程(3)的第1項,在兩端同時乘以A-1ut,再在區(qū)域Ω上積分可得 由(6)式可得〈ut〉=0,所以有 又由(H3)、(6)式以及H?lder不等式可得 即 所以 (11) 利用方程(3)的第4、3、2項以及格林公式得 由上式及(11)式有 所以 (12) 構(gòu)建泛函 所以有 即 其中C0是與α相關(guān)的正常數(shù). 由(9)、(12)式可得 又由Φ(t)的構(gòu)造可得 所以 所以 有 其中Cf″、Cf?、C2為正常數(shù).又由(5)式可得 所以 故一定存在只依賴于C1、C2、α的常數(shù)C(C充分大)使得 (14) 在(14)式兩端同時乘以e,再在s∈(0,t)上積分得: 所以有 由(13)式可得 由引理1.1,取m=n,l=2,q(s)=q=C,h=0,R=E(0)(1+C0n)可得 (15) 取?a∈(0,1),使得 由(13)及(15)式有 參考文獻(xiàn) [1] CAHN J W, Hilliard J E. Free energy of nonuniform system. I. interfacial energy[J]. J Chem Phys,1958,28:258-267. [2] NOVICK-COHEN A. The Cahn-Hilliard equation[C]//Handb Differ Equ IV: Evolutionary Equations. Amsterdam:Elsevier/North-Holland,2008:201-228. [3] TEMAM R. Infinite-Dimensional Dynamical Systems in Mechanics and Physics[M]. New York:Springer-Verlag,1997. [4] EFENDIEV M, MIRANVILLE A, ZELIK S. Exponential attractors for a singularly perturbed Cahn-Hilliard system[J]. Math Nachr,2004,272:11-31. [5] RYBKA P, HOFFMANN K H. Convergence of solutions to Cahn-Hilliard equation[J]. Commun PDE,1999,24:1055-1077. [7] GILARDI G, MIRANVILLE A, SCHIMPERNA G. On the Cahn-Hilliard equation with irregular potentials and dynamic boundary conditions[J]. Commun Pure Appl Anal,2009,8:881-912. [8] MIRANVILLE A, ZELIK S. Exponential attractors for the Cahn-Hilliard equation with dynamic boundary conditions[J]. Math Models Appl Sci,2005,28:709-735. [9] MIRANVILLE A, ZELIK S. The Cahn-Hilliard equation with singular potentials and dynamic boundary conditions[J]. Discrete Contin Dyn Syst,2010,28:275-310. [10] PRüSS J, RACKE R, ZHENG S. Maximal regularity and asymptotic behavior of solutions for the Cahn-Hilliard equation with dynamic boundary conditions[J]. Ann Mat Pure Appl,2006,185(4):627-648. [11] RACKE R. The Cahn-Hilliard equation with dynamic boundary conditions[C]//Nonlinear Partial Differential Equations and Their Applications. GAKUTO Internat: Ser Math Sci Appl. Tokyo:Gakkotosho,2004,20:266-276. [12] RACKE R, ZHENG S. The Cahn-Hilliard equation with dynamic boundary conditions[J]. Adv Diff Eqns,2003,8:83-110. [13] WU H, ZHENG S. Convergence to equilibrium for the Cahn-Hilliard equation with dynamic boundary conditions[J]. J Diff Eqns,2004,204:511-531. [14] LORENZI A, ROCCA E. Weak solutions for the fully hyperbolic phase-field system of conserved type[J]. J Evol Eqns,2007,7:59-78. [15] CONTI M, COTI ZELATI M. Attractors for the non-viscous Cahn-Hilliard equation with memory in 2D[J]. Nonlinear Anal,2010,72:1668-1682. [16] VERGARA V. A conserved phase field system with memory and relaxed chemical potential[J]. J Math Anal Appl,2007,328:789-812. [17] GALENKO P, JOU D. Diffuse-interface model for rapid phase transformations in nonequilibrium systems[J]. Phys Rev,2005,E6125(71):1-13. [18] GALENKO P, JOU D. Kinetic contribution to the fast spinodal decomposition controlled by diffusion[J]. Physics,2009,A388:3113-3123. [19] GALENKO P, LEBEDEV V. Analysis of the dispersion relation in spinodal decomposition of a binary system[J]. Philos Mag Lett,2007,87:821-827. [20] GALENKO P, LEBEDEV V. Local nonequilibrium effect on spinodal decomposition in a binary system[J]. Int J Thermodyn,2008,11:21-28. [21] GALENKO P, LEBEDEV V. Nonequilibrium effects in spinodal decomposition of a binary system[J]. Phys Lett,2008,A372:985-989. [22] LECOQ N, ZAPOLSKY H, GALENKO P. Evolution of the structure factor in a hyperbolic model of spinodal decomposition[J]. Eur Phys J(Special Topics),2009,177:165-175. [23] CAVATERRA C, GAL C G, GRASSELLI M. Cahn-hilliard equations with memory and dynamic boundary conditions[J]. Asymptotic Analysis,2011,71:123-162. [24] SAID-HOUARI B, MESSAOUDI S A. General decay estimates for a cauchy viscoelastic wave problem[J]. Commun Pure Appl Anal,2014,13:1541-1551. [25] PATA V. Uniform estimates of Gronwall type[J]. J Math Anal Appl,2011,373:264-270. [26] GIORGI C, MUOZ RIVERA J E, PATA V. Global attractors for a semilinear hyperbolic equation in viscoelasticity[J]. J Math Anal Appl,2001,260:83-99. [27] GATTI S, PATA V, ZELIK S. A Gronwall-type lemma with parameter and dissipative estimates for PDEs[J]. Nonlinear Anal,2009,70:2337-2343. 2010 MSC:35B30; 35B99 (編輯周俊) The Energy Decay Estimates for Cahn-Hilliard Equation with General Memory Kernel TU Xinyu,PU Zhilin (CollegeofMathematicsandSoftwareScience,SichuanNormalUniversity,Chengdu610066,Sichuan) Abstract:In this paper the energy decay estimates for Cahn-Hilliard equation with general memory kernel are studied by a generalized Gronwall-type lemma with . It is proved that the solution of such an equation decays exponentially, and the decay rate can be explicitly calculated. Key words:dynamic boundary conditions; memory kernel; the past history variable; energy decay estimates doi:10.3969/j.issn.1001-8395.2016.01.001 中圖分類號:O213.2; O226 文獻(xiàn)標(biāo)志碼:A 文章編號:1001-8395(2016)01-0001-07 *通信作者簡介:蒲志林(1963—),男,教授,主要從事偏微分方程的研究,E-mail:puzhilin908@sina.com 基金項目:國家自然科學(xué)基金(11401409) 收稿日期:2015-03-11