龔 歡

(溫州大學(xué)數(shù)理與電子信息工程學(xué)院，浙江溫州 325035)

在本文中，考慮以下的流體動(dòng)力學(xué)模型：

其中，t∈(0,T]，Ω?R2是一個(gè)有界凸區(qū)域，且邊界光滑.u(x,t):Ω×[0,T]→?2和p(x,t):Ω×[0,T]→?分別代表速度和壓力.b(x,t):Ω×[0,T]→S，其中S是?2中的單位圓，代表著宏觀分子取向.f:Ω×[0,T]→?2為外力.常數(shù)μ,λ,γ為三個(gè)正常數(shù)分別表示著粘性系數(shù)、彈性系數(shù)和弛豫時(shí)間系數(shù).?b⊙?b為2×2矩陣.如果b是一個(gè)常映射，則上述系統(tǒng)簡(jiǎn)化為不可壓縮Navier-Stokes方程[1]，如果u=0，則上述系統(tǒng)簡(jiǎn)化為熱調(diào)和映照[2].

為了研究方程（1）-（3），賦予如下的初始條件和邊界條件：

其中n表示在?Ω中單位外法向量，這里要求初始向量函數(shù)u0和b0滿足條件divu0=0和|b0|=1.

方程（1）-（3）首先是由林芳華教授在文獻(xiàn)[3]中作為Ericksen-Leslie模型的簡(jiǎn)化模型導(dǎo)出的，用于描述由Ericksen[4-5]和Leslie[6]所提出的向列型液晶流動(dòng)模型，它是在流動(dòng)速度u和類(lèi)似于旋轉(zhuǎn)液晶流的微觀取向b的影響下對(duì)材料時(shí)間演化的宏觀連續(xù)描述.但研究系統(tǒng)（1）-（5）的困難之一是非線性約束|b|=1，通常的方法是采用Ginzburg-Landau加罰函數(shù)去代替（2）式中的2|?b|b，而不考慮非線性約束|b|=1.此時(shí)考慮如下的加罰問(wèn)題：

其中ε＞0為加罰參數(shù).（6）-（8）式首先是由Lin和Liu在文[7]中引入的，文[7]證明了三維問(wèn)題全局弱解和局部強(qiáng)解的存在性.然而由于解的估計(jì)很大程度上取決于ε，所以在（6）-（8）式中當(dāng)ε趨于0時(shí)，加罰問(wèn)題方程（6）-（8）的解是否趨于原問(wèn)題方程（1）-（3）的解仍不清楚.

近幾年許多學(xué)者對(duì)向列液晶流動(dòng)問(wèn)題的數(shù)值方法也進(jìn)行了研究，Liu和Walkington在文獻(xiàn)[8]中對(duì)加罰問(wèn)題方程（6）-（8）的數(shù)值方法首先進(jìn)行了研究，并采用3?Hermite有限元來(lái)逼近bε.為了避免使用Hermite有限元，隨后Liu和Walkington又在文[9]中研究了混合有限元方法，并構(gòu)造了全隱格式，雖然該格式是無(wú)條件穩(wěn)定的，但每個(gè)時(shí)間步必須求解非線性問(wèn)題而需消耗大量的數(shù)值求解時(shí)間，正如文[10]中所指出的，在文[8-9]中導(dǎo)出的誤差估計(jì)依賴(lài)于ε，這將導(dǎo)致非常小的時(shí)間步長(zhǎng)為避免時(shí)間步長(zhǎng)的強(qiáng)約束條件，Becker等在文[10]中通過(guò)引入研究了一種新的混合有限元方法，所提出的完全離散格式僅在求解wε時(shí)是非線性的，并且滿足離散的能量不等式.在文[10]中，Becker等證明了當(dāng)h,τ→0時(shí)數(shù)值解的適定性和收斂性，并沒(méi)有給出誤差估計(jì).最近An和Su在文[11]中研究了求解（1）-（5）式的線性化半隱全離散有限元算法，在構(gòu)造該算法時(shí)沒(méi)有考慮非線性約束條件|b|=1，并證明了該全離散有限元算法具有最優(yōu)階的有限元時(shí)空誤差估計(jì).

本文將基于一種新的加罰方法，在不考慮非線性約束條件|b|=1的情況下，構(gòu)造新的一階線性化時(shí)間半離散格式.在時(shí)間步長(zhǎng)和加罰參數(shù)滿足τ=ο(ε)的條件下，我們證明了所構(gòu)造的半離散算法具有一階時(shí)間收斂精度.

1 預(yù)備知識(shí)

本文采用標(biāo)準(zhǔn)的Sobolev空間記號(hào)[12].用(.,.)表示L2(Ω)的內(nèi)積，采用C,C0,C1,C2…表示正常數(shù)，它們可依賴(lài)于u,p,b,f,μ,λ,γ，但不依賴(lài)于加罰參數(shù)ε＞0和時(shí)間步長(zhǎng)τ.

引入空間：

定義雙線性形式：

定義三線性形式：

通過(guò)積分可得到：

記正交投影算子PH:L2(Ω)→H，那么Stokes算子A定義為[1]：

關(guān)于問(wèn)題（1）-（5）的適定性，首先引入如下定理.

定理1 若u0∈D(A)，b0∈H3(Ω)2和|b0|=1在Ω內(nèi)幾乎處處成立，對(duì)于給定的f∈L∞(0,T;H)∩L2(0,T;L4(Ω)2)，則存在T*＜T使得（1）-（5）式存在唯一的局部強(qiáng)解(u,p,b)并滿足[13]：

記0=t0＜t1＜…tN=T是時(shí)間間隔[0,T]的一個(gè)均勻化分割，時(shí)間步長(zhǎng)，τ=T/N和tn=nτ.記un=u(x,tn)，pn=p(x,tn)，bn=b(x,tn)，fn=f(x,tn).對(duì)序列記取初始值B0=b0，U0=u0，構(gòu)造如下的線性化加罰迭代格式：

為了研究上述迭代格式的時(shí)間收斂階，本文將應(yīng)用如下的不等式：

2 主要結(jié)論

從（25）式中減去（18）式，并從（24）式、（26）式中減去（19）式得誤差方程為：

本文的主要結(jié)論如下.

定理2 在定理1的正則性條件下，對(duì)足夠小的ε＞0，若時(shí)間步長(zhǎng)τ滿足τ=ο(ε)，則有：

證明：因?yàn)閡0=U0=u0和b0=B0=b0，因此（30）式和（31）式對(duì)m=0成立.假設(shè)（30）式和（31）式對(duì)m≤n成立，根據(jù)數(shù)學(xué)歸納法，那么只需要證明（30）式和（31）式對(duì)m≤n+1也成立.（28）式乘以然后在Ω上積分可得到：

現(xiàn)在開(kāi)始估計(jì)上式右端.根據(jù)

并使用（15）式，（16）式，（22）式，（30）式（對(duì)m≤n成立），Holder’s不等式和Young’s不等式，可估計(jì)I1如下：

其中C1＞0與C0無(wú)關(guān).這里使用

再由|bn|=1可得：

利用Holder’s不等式和Young’s不等式，可估計(jì)I3如下：

這里使用

把I1,I2,I3和I4的估計(jì)帶入（32）式得到：

利用Holder’s不等式和Young’s不等式有：

其中C2＞0與C0無(wú)關(guān).

再根據(jù)

有：

這里使用了（30）-（31）式（對(duì)m≤n成立）和

把I5,I6,I7帶入（36）式中可得：

其中C3＞0與C0無(wú)關(guān).

把（35）式，（38）式，（39）式加起來(lái)，得：

使用離散的Gronwall’s不等式[14]，并根據(jù)τ=ο(ε)，存在足夠小的τ0＞0使得當(dāng)τ＜τ0時(shí)有：

這樣就證明了定理2的結(jié)論.