"/>
  • 
    

    
    

      99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

      粘性Cahn-Hilliard方程的半線性Crank-Nicolson格式

      2018-04-28 02:33:49
      關(guān)鍵詞:粘性先驗(yàn)邊界條件

      李 娟

      (南京審計(jì)大學(xué) 金審學(xué)院 基礎(chǔ)部, 江蘇 南京 210023)

      ut-αuxxt+βuxxxx-(u3-u)xx=0,
      x∈Ω, t>0,

      (1)

      邊界條件

      u(x,t)=uxx(x,t)=0,
      x∈?Ω, t>0,

      (2)

      初值條件

      u(x,0)=u0(x), x∈Ω,

      (3)

      其中,Ω為有界區(qū)間,α≥0,β>0,u表示兩相流的濃度,u3-u是內(nèi)部固有化學(xué)勢,β為界面能量參數(shù),α表示粘性系數(shù).當(dāng)α=0時(shí),方程(1)為標(biāo)準(zhǔn)Cahn-Hilliard方程.本文考慮α>0的情形.

      在研究玻璃和聚合物兩相分離的連續(xù)模型中,將分子間的摩擦力考慮進(jìn)來,提出了粘性Cahn-Hilliard方程.作為重要且應(yīng)用更廣泛的相場模型方程,學(xué)者已對(duì)其有了廣泛的研究.文獻(xiàn)[1]指出文獻(xiàn)[2]忽略了反映粘性影響的項(xiàng)αuxxt,提出了粘性Cahn-Hilliard方程;文獻(xiàn)[3]指出粘性Cahn-Hilliard可被看做相變的奇異極限;文獻(xiàn)[4]研究了該方程的亞穩(wěn)內(nèi)部層動(dòng)態(tài);文獻(xiàn)[5]對(duì)1、2、3維粘性Cahn-Hilliard方程進(jìn)行了理論分析;文獻(xiàn)[6]建立了該方程穩(wěn)態(tài)解的Morse分解和全局吸引子的結(jié)構(gòu);文獻(xiàn)[7]討論了帶有時(shí)間周期勢的粘性Cahn-Hilliard方程的解的存在性、唯一性和漸進(jìn)性;文獻(xiàn)[8]給出了基于邊界條件的最優(yōu)控制,并證明了方程最優(yōu)控制解的存在性;文獻(xiàn)[9]研究了方程中通道對(duì)小參數(shù)的限制;文獻(xiàn)[10]研究了帶有動(dòng)態(tài)邊界條件和雙障礙源的粘性Cahn-Hilliard方程的最優(yōu)控制問題;文獻(xiàn)[11]研究了非標(biāo)準(zhǔn)的粘性Cahn-Hilliard方程系統(tǒng)解的適定性和長時(shí)間行為;文獻(xiàn)[12-13]討論了黏性Cahn-Hilliard方程弱解的存在性.

      2 差分格式的建立

      下面對(duì)問題(1)~(3)建立差分格式.令v=u2ux,則方程(1)等價(jià)于

      ut-αuxxt+βuxxxx+uxx-3vx=0,
      x∈[0,L], t>0,

      (4)

      v=u2ux, x∈[0,L], t>0.

      (5)

      (6)

      (7)

      (8)

      其中

      為截?cái)嗾`差.將(7)和(8)式分別代入(6)式可得

      (9)

      (10)

      其中

      假定問題(1)~(3)的解是適當(dāng)光滑的,則存在正常數(shù)c0,使得

      (11)

      由邊界條件和初值條件可得

      (12)

      (13)

      (14)

      (15)

      (16)

      (17)

      在上面差分格式中,可利用(14)式求得第一、二時(shí)間層的數(shù)值解.由于其為非線性方程組,故需迭代求解.當(dāng)un、un-1、un-2(n≥2)為已知時(shí),(15)式為關(guān)于un+1的線性方程組,可通過解線性方程組求得其余時(shí)間層的數(shù)值解.

      3 差分格式解的先驗(yàn)估計(jì)

      下面將給出差分格式解的先驗(yàn)估計(jì)式.設(shè)Vh={u=(u0,u1,…,uM),u0=uM=0}.對(duì)任意的u,v∈Vh,引入如下內(nèi)積和范數(shù):

      為便于先驗(yàn)估計(jì)式的推導(dǎo),引入如下引理.

      引理3.1[20]對(duì)任意的u,v∈Vh,存在

      (u,δxv)=-(δxu,v), (u,xv)=-(xu,v),
      ‖xu‖≤‖δxu‖.

      引理3.2[21]對(duì)任意的v∈Vh和p≥2,存在

      ‖v‖p≤μ(‖v‖1-γ‖δxv‖γ+‖v‖),

      下面,對(duì)于差分格式的數(shù)值解un,給出‖un‖,‖δxun‖的先驗(yàn)估計(jì)式,進(jìn)一步得出‖un‖∞的有界性.

      (18)

      利用引理3.1可得

      (19)

      由內(nèi)積定義知

      (20)

      由引理3.1和柯西不等式可得

      (21)

      將(20)和(21)式代入(19)式可得

      (22)

      Fn=‖un‖2+α‖δxun‖2.

      由(22)式可得

      (23)

      (24)

      (25)

      從而有

      (26)

      (27)

      結(jié)合(24)和(27)式可知,對(duì)于0≤n≤N-1,(27)式均成立.

      由離散的Gronwall不等式可得

      從而有

      (28)

      (29)

      根據(jù)上面定理,易證差分格式的解是有界的.在引理3.2中,取p=∞有

      (30)

      由定理3.1中的估計(jì)式(28)~(30)可得,存在不依賴于h、τ的正常數(shù)B,使得

      ‖un‖∞≤B, 0≤n≤N,

      (31)

      即差分格式的解是有界的.

      4 差分格式的收斂性

      (32)

      (33)

      (34)

      (35)

      定理4.1假設(shè)問題(1)~(3)的解適當(dāng)光滑.差分格式(14)~(17)式的解收斂于問題的精確解,收斂階為時(shí)間、空間方向二階收斂,即存在正常數(shù)c,使得

      ‖en‖∞≤c(τ2+h2), 0≤n≤N.

      (36)

      證明采用數(shù)學(xué)歸納法證明定理成立.由(35)式知,當(dāng)n=0時(shí),(36)式成立.

      I) 當(dāng)n=1,2時(shí),‖en‖∞≤c(τ2+h2)成立.

      (37)

      由引理3.1可得

      (38)

      下面估計(jì)(38)式右端的第二項(xiàng)和第三項(xiàng).對(duì)于第二項(xiàng)有

      (39)

      從而有

      (40)

      對(duì)于第三項(xiàng)有

      (41)

      將(40)和(41)式代入(38)式,可得

      (42)

      由引理3.1和柯西不等式可得

      將上式代入(42)式可得

      (43)

      令En=‖en‖2+α‖δxen‖2.由(43)式可得

      其中

      從而有

      (44)

      由(44)式可得

      從而有

      (45)

      (46)

      由引理3.2可得

      即當(dāng)n=1,2時(shí),(36)式成立.

      (II) 假設(shè)定理結(jié)論對(duì)0≤n≤m均成立,證明‖em+1‖≤c(τ2+h2)成立.

      (47)

      從而有

      (48)

      下面估計(jì)上式右端的第二、三項(xiàng).對(duì)于第二項(xiàng)有

      成立,從而有

      (49)

      對(duì)于第三項(xiàng)有

      (50)

      將(49)和(50)式代入(48)式可得

      (51)

      注意到

      并將其代入到(51)式可得

      (52)

      其中

      在(52)式中,用l代替n,并對(duì)l從2到n求和可得

      從而有

      (53)

      (54)

      其中

      由(53)式可得

      由離散的Gronwall不等式可得

      En+1≤c4exp(c5T)(τ2+h2)2, 2≤n≤m,

      從而有

      應(yīng)用引理3.2可得

      即當(dāng)n=m+1時(shí),(36)式成立.由數(shù)學(xué)歸納法,定理結(jié)論成立.

      5 數(shù)值算例

      1) 驗(yàn)證差分格式的精度.由于無法計(jì)算該問題的精確解,為驗(yàn)證數(shù)值收斂階,定義如下最大模誤差

      對(duì)于充分小的空間步長h,定義時(shí)間收斂階

      對(duì)于充分小的時(shí)間步長τ,定義空間收斂階

      表 1 差分格式(14)~(17)的最大模誤差和時(shí)間收斂階

      表 2 差分格式(14)~(17)的最大模誤差和空間收斂階

      6 小結(jié)

      圖 1 當(dāng)β=0.5,α=0.5,0.05,0.005,0時(shí),(1)~(3)式的數(shù)值解曲面

      [1] NOVICK-COHEN A. On the viscous Cahn-Hilliard equation[C]//Material Instabilities in Continuum and Related Mathematical Problems. Oxford:Oxford Univ Press,1988.

      [2] CAHN J W, HILLIARD J E. Free energy of a non-uniform system I. interface free energy[J]. J Chem Phys,1958,28(2):258-267.

      [3] BAI F, ELLIOTT C M, GARDINER A, et al. The viscous Cahn-Hilliard equation. part I:computations[J]. Nonlinearity,1995,8(2):131-160.

      [4] REYNA L G, WARD M J. Metastable internal layer dynamics for the viscous Cahn-Hilliard equation[J]. Methods Appl Anal,1995,2(3):285-306.

      [5] ELLIOTT C M. Viscous Cahn-Hilliard equation II. analysis[J]. J Diff Eqns,1996,128(2):387-414.

      [6] GRINFELD M, NOVICK-COHEN A. The viscous Cahn-Hilliard equation:morse decomposition and structure of the global attractor[J]. Trans Am Math Soc,1999,351(6):2375-2406.

      [7] LI Y, YIN J. The viscous Cahn-Hilliard equation with periodic potentials and sources[J]. J Fixed Point Theory Appl,2011,9(1):63-84.

      [8] ZHAO X, LIU C. Optimal control problem for viscous Cahn-Hilliard equation[J]. Nonlinear Analysis,2011,74(17):6348-6357.

      [9] THANH B L T, SMARRAZZO F, TESEI A. Passage to the limit over small parameters in the viscous Cahn-Hilliard equation[J]. J Math Anal Appl,2014,420(2):1265-1300.

      [10] COLLI P, FARSHBAF-SHAKERM H, GILARDI G, et al. Optimal boundary control of a viscous Cahn-Hilliard system with dynamic boundary condition and double obstacle potentials[J]. SIAM J Control Optim,2015,53(4):2696-2721.

      [11] COLLI P, GILARDI G, PODIO-GUIDUGLI P, et al. Well-posedness and long-time behavior for a nonstandard viscous Cahn-Hilliard system[J]. SIAM J Appl Math,2011,71:1849-1870.

      [12] 喻朝陽. 黏性Cahn-Hilliard方程弱解的存在唯一性[J]. 四川師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2012,35(2):227-230.

      [13] 姜金平,董超雨. 粘性Cahn-Hilliard方程在H1中的弱解存在性[J]. 貴州大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,32(6):1-3.

      [14] MOHYUD-DIN S T, YLDRM A, SARAYDN S. Approximate series solutions of the viscous Cahn-Hilliard equation via the homotopy perturbation method[J]. World Appl Sci,2010,11:813-818.

      [15] CHOO S M, CHUNG S K, LEE Y J. A conservative difference scheme for the viscous Cahn-Hilliard equation with a nonconstant gradient energy coefficient[J]. Appl Numer Math,2004,51(2):207-219.

      [16] CHOO S M, KIM Y H. Finite element scheme for the viscous Cahn-Hilliard equation with a nonconstant gradient energy coefficient[J]. Appl Math Comput,2005,19(1):385-395.

      [17] MOMANI S, ERTURK V S. A numerical scheme for the solution of viscous Cahn-Hilliard equation[J]. Numer Methods Partial Diff Eqns,2008,24(2):663-669.

      [18] SHIN J, CHOI Y, KIM J. An unconditional stable numerical method for the viscous Cahn-Hilliard equation[J]. Discrete Contin Dyn Sys,2014,B19(6):1737-1747.

      [19] CHOO S M, CHUNG S K. A conservative nonlinear difference scheme for the viscous Cahn-Hilliard equation[J]. Appl Math Comput,2004,16(1):53-68.

      [20] 孫志忠. 偏微分方程數(shù)值解法[M]. 2版. 北京:科學(xué)出版社,2012.

      [21] ZHOU Y. Applications of Discrete Functional Analysis of Finite Difference Method[M]. New York:International Academic Publishers,1990.

      [22] LI J, SUN Z Z, ZHAO X. A three level linearized compact difference scheme for the Cahn-Hilliard equation[J]. Sci China Math,2012,55(4):805-826.

      [23] 李娟,高廣花. 求解Fisher-Kolomogorov方程的三層線性化緊差分格式[J]. 西南民族大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,41(5):634-639.

      [24] 李娟,高廣花. 二維擴(kuò)展Fisher-Kolmogorov方程的線性化緊差分格式的最大模誤差分析[J]. 西南師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2017,42(3):12-21.

      猜你喜歡
      粘性先驗(yàn)邊界條件
      一類具有粘性項(xiàng)的擬線性拋物型方程組
      一類帶有Stieltjes積分邊界條件的分?jǐn)?shù)階微分方程邊值問題正解
      帶有積分邊界條件的奇異攝動(dòng)邊值問題的漸近解
      帶粘性的波動(dòng)方程組解的逐點(diǎn)估計(jì)
      基于無噪圖像塊先驗(yàn)的MRI低秩分解去噪算法研究
      基于自適應(yīng)塊組割先驗(yàn)的噪聲圖像超分辨率重建
      粘性非等熵流體方程平衡解的穩(wěn)定性
      基于平滑先驗(yàn)法的被動(dòng)聲信號(hào)趨勢項(xiàng)消除
      先驗(yàn)的廢話與功能的進(jìn)路
      帶Robin邊界條件的2維隨機(jī)Ginzburg-Landau方程的吸引子
      濮阳市| 普格县| 靖安县| 凤城市| 巴楚县| 宁河县| 河间市| 阿克陶县| 密云县| 乌拉特后旗| 托克逊县| 乌鲁木齐市| 浑源县| 信阳市| 青阳县| 额尔古纳市| 成武县| 连城县| 巴彦县| 汽车| 镇原县| 河池市| 洞头县| 临朐县| 临城县| 盐亭县| 灵武市| 老河口市| 延安市| 准格尔旗| 武汉市| 汉沽区| 丘北县| 仪征市| 共和县| 井研县| 门源| 南漳县| 彭阳县| 正阳县| 景泰县|