Processing math: 17%

数学物理学报, 2021, 41(1): 126-141 doi:

论文

不可压双曲型液晶的Ericksen-Leslie方程的大Ericksen数的极限的形式分析

程峰,

A Formal Analysis on the Large Ericksen Number Limit for the Incompressible Hyperbolic Ericksen-Leslie System of Liquid Crystals

Cheng Feng,

收稿日期: 2020-03-15  

Received: 2020-03-15  

作者简介 About authors

程峰,E-mail:fengcheng@hubu.edu.cn , E-mail:fengcheng@hubu.edu.cn

Abstract

In this paper, we consider the parameterized incompressible hyperbolic Ericksen-Leslie system of liquid crystals. Formally, letting the parameter vanish, we prove that the limit system admits a local classical solution. Moreover, we formally obtained an estimate on the difference between the parameterized Ericksen-Leslie's incompressible hyperbolic liquid crystal model and the corresponding limit system, which corresponds to a formal energy estimate of the difference of the classical solutions in L2 space.

Keywords: Liquid crystals ; Ericksen-Leslie system ; Classical solution ; Large Ericksen number limit

PDF (428KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

程峰. 不可压双曲型液晶的Ericksen-Leslie方程的大Ericksen数的极限的形式分析. 数学物理学报[J], 2021, 41(1): 126-141 doi:

Cheng Feng. A Formal Analysis on the Large Ericksen Number Limit for the Incompressible Hyperbolic Ericksen-Leslie System of Liquid Crystals. Acta Mathematica Scientia[J], 2021, 41(1): 126-141 doi:

1 引言

不可压液晶的流体动力学理论最早是由Ericksen[2-4]和Leslie[8-9]在1960年代建立的.后来, Hardt, Kinderlehrer和林芳华在文献[13]中最初开始抛物型Ericksen-Leslie方程的数学研究.自此以后, 有大量关于Ericksen-Leslie方程及其逼近形式的存在性、正则性以及部分正则性的研究, 例如文献[22-25].对于更一般的抛物型Ericksen-Leslie方程, 关于二维空间R2中的局部适定性和整体适定性请参见文献[15, 26], 而关于弱解的存在性和唯一性参见文献[14, 21, 27].

对于液晶的双曲型的Ericksen-Leslie方程, 高维情况下的数学研究结果非常少. De Anna和Zarnescu在文献[1]中研究了一个相关的惯性的液晶的Qian-Sheng模型.他们得出了能量守恒律, 同时也证明了有界初值下的局部适定性和小初值的同时系数满足某种衰减性质下的整体适定性.对于无粘性的Qian-Sheng模型, Feireisl等在文献[5]中证明了大初值下的整体耗散解的存在性.最近, 江宁和罗益龙在文献[19]中证明了双曲型不可压Ericksen-Leslie方程的有限初值下的局部经典解的存在性和唯一性.而在Leslie系数的一些约束条件和小初始能量的假设下, 他们也可以证明整体经典解的存在性和唯一性.

最近, Wu, Xu和Zarnescu在文献[29]中研究了不可压缩向列型液晶流的Beris-Edwards系统在控制流动时特征值范围的保留效应和缺陷的动态出现.特别地, 他们为Beris-Edwards系统在大Ericksen数的极限下的部分解耦提供了定性的理解.本文的目标是为了探讨双曲型Ericksen-Leslie方程的解的大Ericksen数的极限.在本文中, 我们将要研究了在大Ericksen数和Reynold数极限下, 如下双曲型Ericksen-Leslie不可压液晶方程的部分解耦问题

{t(ρu)+div(ρuu)+p=div(Σ1+Σ2+Σ3),divu=0,ρ¨d=κΔd+Γd+λ1(˙d+Bd)+λ2Ad,
(1.1)

其中ρ>0表示流体的密度常数, u(t,x)Rn表示流体的速度向量场, d(t,x)Rn表示液晶分子在几何约束|d|=下的方向场.这里, 为棒状向列相液晶分子的平均长度. p(t,x)R表示标量压力场.这里, ˙d¨d分别表示d的一阶物质导数和二阶物质导数, 即

˙d=td+ud,¨d=t˙d+u˙d .
(1.2)

记号

A=12(u+u),B=12(uu)
(1.3)

分别表示表示应变张量的速率和应变速率的反对称部分.记号Σi(i=1,2,3)定义如下

Σ1=12μ4(u+u),Σ2=12κ|d|2Eκ(dd),Σ3=˜σμ(u,d,˙d),
(1.4)

其中En×n单位矩阵, 矩阵(dd)的元素为(dd)ij=idjd, ˜σ(u,d,˙d)是应力张量具有如下分量形式

(˜σ(u,d,˙d))ij=μ1dkdpAkpdidj+μ2di(˙dj+Bjkdk)+μ3dj(˙di+Bikdk)+μ5diAjkdk+μ6djAikdk,
(1.5)

这里参数μi (i=1,2,3,5,6)是Leslie系数, μ4就是通常的粘性系数, κ是测量液晶分子的弹性的Frank常数.通常, λ1,λ2满足λ1=μ2μ3, λ2=μ5μ6.并且, μ2,μ3,μ5,μ6满足关系μ2+μ3=μ6μ5,这一关系被称为Parodi's关系, 它是由Onsager互反关系导出的, Onsager互反关系表示了在非平衡的热力学系统中流体和力之间某些关系的相等关系, 详细内容请参见文献[28].拉格朗日乘子Γ具有如下表示

ΓΓ(u,d,˙d)=12(|˙d|2+κ|d|2λ2dAd).
(1.6)

从原始Ericksen-Leslie模型到不可压缩液晶系统(1.1)的严格推导归功于文献[19-20].由于我们要研究方程(1.1)的渐近行为, 因此研究它的无量纲形式是十分必要的.为此, 我们需要开展如下的无量纲程序.

L,TU分别表示宏观尺度下的单位长度、时间和流体速度, 其中U=LT.ρ,cμ分别表示单位的密度, 单位的声速以及单位的粘性系数.令κ=K2表示单位的Frank常数, 其中K是弹性系数.

为了无量纲化方程(1.1), 我们令

x=Lˆx,t=Tˆt,u=Uˆu,ρ=ρˆρ,d=ˆd,κ=κˆκ,
(1.7)

μ4=μˆμ4,ξ=μˆξ,μ1=μ4ˆμ1,μi=μ2ˆμi,i=2,3,5,6.
(1.8)

分别记马赫数、Reynolds数、Ericksen数以及惯性系数如下

M=Uc,R=ρLUμ,E=μLUκ2,χ=ρU2κ.
(1.9)

现在将表达式(1.7)–(1.8)代入方程(1.1), 借助(1.9)式中的定义, 我们将方程(1.1)改写为如下无量纲形式(去掉所有ˆ)

{ρtu+ρ(uu)+1M2p=1RdivΣ1+1R1EdivΣ2+1RdivΣ3,divu=0,χEρ¨d=κEΔd+Γd+λ1(˙d+Bd)+λ2Ad,
(1.10)

其中Γ=χEρ|˙d|2+κE|d|2λ2dAd.

从上面无量纲形式中, 通过选取不同的尺度, 关于相关的渐进模型已有丰富的研究成果.例如, 让惯性系数χ趋于0, 方程(1.10)形式上收敛到更广泛研究的抛物型Ericksen-Leslie方程, 参见文献[6, 12, 18].考虑Q -张量形式, 例如Beris-Edwards方程, Wu-Xu-Zarnescu在文献[29]中通过让Reynolds数和Ericksen数趋于无穷研究了Beris-Edwards方程的所谓的大Ericksen数极限(在χ=0的情形).

在本文中, 我们对双曲型液晶方程(1.10)的大Ericksen数极限进行了初步的研究.为了表述的简单, 我们令参数ρ,M, χ以及κ1, 并且令Reynolds数R=1ε2和Ericksen数E=1ε (这里选取不同的次数是为了计算的简单).于是, 方程(1.10)可写为

{tuε+uεuε+pε=ε2divΣε1+ε3divΣε2+ε2divΣε3,divuε=0,ε¨dε=εΔdε+Γεdε+λ1(˙dε+Bεdε)+λ2Aεdε,
(1.11)

其中, dε满足几何约束|dε|=1,

Γε=ε|˙dε|2+ε|dε|2λ2dεAεdε.

方程(1.11)有如下初始条件

uε|t=0=u0(x),dε|t=0=d0(x),˙dε|t=0=˜d0(x),
(1.12)

其中d0(x)˜d0(x)满足约束条件和相容性条件

|d0(x)|=1,˜d0(x)d0(x)=0.
(1.13)

对任意的ε>0, 江宁和罗益龙在[19]中证明了, 如果(u0,d0,˜d0)Hs(Rn),s>n2+1满足条件(1.13), 存在Tε>0使得方程(1.11)–(1.12)存在局部经典解

(uε,dε,˙dε)L([0,Tε];Hs(Rn)).

在大初值下, 方程(1.11)–(1.12)的整体适定性仍是个公开的问题.

方程(1.11)–(1.12)的所谓的大Ericksen数极限问题对应的是该方程的经典解(uε,dε,˙dε)随着ε0的渐近行为.形式上, 在方程(1.11)中令ε=0, 我们得到如下的极限方程

{tu+uu+p=0,divu=0,γd+λ1(˙d+Bd)+λ2Ad=0,
(1.14)

这里d满足约束条件|d|=1, 拉格朗日乘子具有如下形式

γ=λ2dAd.

注意到方程组(1.14)的第三个方程对时间变量t只有一阶导数, 因此方程(1.14)的初始条件变为

u|t=0=u0(x),d|t=0=d0(x),
(1.15)

其中d0(x)满足约束|d0(x)|=1.

我们可以猜测到随着ε0, 方程(1.11)–(1.12)的解(uε,dε,˙dε)应该会在某种意义下收敛到方程(1.14)–(1.15)的解(u,d).而本文的一个主要目的就是在经典解的框架下给出他们的误差(uεu,dεd)的一个能量估计.

本文的余下部分是如下组织的.在第2节, 我们将给出一些记号和预备知识.在第3节, 我们将证明极限方程的经典解的局部存在性.在第4节, 我们将陈述我们的主要定理, 同时给出主要定理的证明.

2 记号和预备知识

为了陈述我们的主要结果以及表述的简单, 我们引进一些记号和预备知识.

空间Lp(Tn), 1p表示通常的标量函数或者向量函数的Lebesgue空间, 它的范数记为. p = 2 , L^2({\mathbb T}^n) 是Hilbert空间, 我们记 \langle\cdot, \cdot\rangle 为相应的标量函数或向量函数的 L^2 -内积.我们记 W^{k, p}({\mathbb T}^n) 为通常的Sobolev空间, 相应的范数记为 \|\cdot\|_{W^{k, p}} . p = 2 , 我们用 H^k({\mathbb T}^n) 来表示 W^{k, 2}({\mathbb T}^n) , 相应的内积记为 \langle\cdot, \cdot\rangle_k .给定一个Banach空间 X 以及相应的范数 \|\cdot\|_{X} , 我们用 L^p(0, T;X) , 1\leq p\leq\infty 来表示定义在 (0, T) 上取值在 X 上且满足 \int_0^T \|f(t)\|_{X}^p{\rm d}t<\infty 的所有函数的集合.在本文中, 我们用 C 来表示一个常数, 它可能在不同地方会有不同的取值.

由于我们在Sobolev空间中研究, 我们会很频繁地使用如下引理.

引理2.1  Sobolev空间中的微积分不等式.

1) \forall\ m\in{\mathbb Z} 存在 c>0 使得, \forall\ u, v\in L^\infty\cap H^m({\mathbb R}^n) , 有

\|uv\|_{H^m}\leq c\{\|u\|_{L^\infty}\|v\|_{H^m}+\|u\|_{H^m}\|v\|_{L^\infty} \},

\sum\limits_{0\leq|\alpha|\leq m}\|\partial^\alpha(uv)-u\partial^\alpha v\|_{L^2}\leq c\{\|\nabla u\|_{L^\infty}\|v\|_{H^{m-1}}+\|u\|_{H^m}\|v\|_{L^\infty} \};

2) \forall\ s>\frac{n}{2} , 存在 c>0 使得 \forall\ u, v\in H^s({\mathbb R}^n) , \|uv\|_{H^s}\leq c\|u\|_{H^s}\|v\|_{H^s}.

该引理的证明可参考文献[10].

3 极限方程的经典解的局部存在性

在这一节, 我们考虑极限方程(1.14)–(1.15).注意到 \dot{\rm d} = \partial_t{\rm d}+u\cdot\nabla{\rm d}, 同时如果 \lambda_1<0 , 我们可以将极限方程(1.14)改写为如下形式

\begin{equation} \left\{ \begin{array}{ll} \partial_tu+u\cdot\nabla u+\nabla p = 0, \\ {\rm{divu}} = 0, \\ { } \partial_t {\rm d}+u\cdot\nabla{\rm d}+{\rm Bd}+\frac{1}{\lambda_1}\gamma {\rm d}+\frac{\lambda_2}{\lambda_1} {\rm Ad} = 0, \end{array} \right. \end{equation}
(3.1)

其中 \gamma = -\lambda_2{\rm d}^\top{\rm A}{\rm d} .初始条件为

\begin{equation} u\big|_{t = 0} = u_0(x), \ \ \ {\rm d}\big|_{t = 0} = {\rm d}_0(x), \ \ \ x\in{\mathbb T}^n, \end{equation}
(3.2)

这里 u, {\rm d} 满足周期性边界条件.

方程组(3.1)–(3.2)是一个弱的耦合系统, 它包含流体速度场 u 的Euler方程和方向场 d 的输运方程.在下面的定理中, 我们将给出方程组(3.1)–(3.2)的经典解的存在性的证明.

定理3.1  令 s>\frac{n}{2}+1, \lambda_1<0, u_0\in H^{s+1}({\mathbb T}^n), {\rm d}_0\in H^s({\mathbb T}^n) .则存在 T>0 , 依赖于 s, n, \|u_0\|_{H^s}, \|{\rm d}_0\|_{H^s} , 使得极限方程存在唯一解

u\in L^\infty\big([0, T];H^{s+1}({\mathbb T}^n) \big), \quad {\rm d}\in L^\infty\big([0, T];H^s({\mathbb T}^n)\big).

并且, 对所有的 0\leq t\leq T , x\in{\mathbb T}^n , {\rm d} 满足几何约束 |{\rm d}(t, x)| = 1 .

  由于方程(3.1)中的流体速度场 u 的Euler方程是自适定的.我们也都知道, 如果 n = 2, s>2 , 对于初值速度场 u_0\in H^s({\mathbb T}^2) , Euler方程是整体适定的, 而如果 n = 3, s> \frac{5}{2} , 对于初始速度 u_0\in H^s({\mathbb T}^3) , 存在 T_1>0 依赖于 s, \|u_0\|_{H^s} , 使得Euler方程存在解 u 满足

u\in L^\infty\big([0, T_1];H^s({\mathbb T}^n)\big),

更详细的讨论, 请参见文献[16-17].因此, 我们只需要考虑给定流体速度场 u 下的液晶分子的方向场 {\rm d} 的方程 {(3.1)}_3 .

下面我们将证明, 如果 {\rm d}_0\in H^s({\mathbb T}^n) 满足 |{\rm d}_0| = 1 , 存在 T_2>0 使得方程 (3.1)_3 存在唯一的经典解 {\rm d}\in L^\infty\big([0, T_2];H^s({\mathbb T}^n)\big) , 并且对 t\in [0, T], x\in{\mathbb T}^n , 都有 |d(t, x)| = 1 .

证明分为两个部分.在第一部分, 我们将证明如果 {\rm d}_0\in H^s({\mathbb T}^n) , 线性的输运方程

\begin{equation} \left\{ \begin{array}{ll} { } \partial_t {\rm d}+u\cdot\nabla{\rm d}+{\rm Bd}+\frac{1}{\lambda_1}\gamma {\rm d}+\frac{\lambda_2}{\lambda_1} {\rm Ad} = 0, \\ {\rm d}\big|_{t = 0} = {\rm d}_0 \end{array} \right. \end{equation}
(3.3)

存在唯一的 {\rm d}\in L^\infty\big([0, T_2];H^s({\mathbb T}^n) \big) .这一部分的证明包含下面三个步骤.

第一步  先验估计.

假设对给定的光滑的 (u, {\rm d}) 是方程(3.3)的一个光滑解.对满足 0\leq |\alpha|\leq s 的所有多重指标 \alpha\in{\mathbb Z}^n , 在方程 {(3.3)}_1 的两边应用微分算子 \partial^\alpha , 然后和 \partial^\alpha {\rm d} L^2 -内积, 我们可以得到

\begin{eqnarray} \frac{1}{2}\frac{\rm d}{{\rm d}t}\|\partial^\alpha {\rm d}\|_{L^2}^2 & = &-\langle\partial^{\alpha}(u\cdot\nabla{\rm d}), \partial^\alpha{\rm d}\rangle-\langle \partial^\alpha({\rm Bd}), \partial^\alpha{\rm d}\rangle-\frac{1}{\lambda_1}\langle\partial^\alpha(\gamma{\rm d}), \partial^\alpha{\rm d}\rangle {}\\ &&-\frac{\lambda_2}{\lambda_1}\langle\partial^\alpha({\rm Ad}), \partial^\alpha{\rm d}\rangle\ . \end{eqnarray}
(3.4)

注意到, 如果 s>\frac{n}{2} , H^s 是一个代数.在方程(3.4)两边对多重指标 0\leq|\alpha|\leq s 求和, 结合引理, 我们可以得到

\begin{eqnarray} \frac{1}{2}\frac{\rm d}{{\rm d}t}\|{\rm d}\|_{H^s}^2 \leq C\|u\|_{H^s}\|{\rm d}\|_{H^s}^2+C\|\nabla u\|_{H^s}\|{\rm d}\|_{H^s}^2+C\|\nabla u\|_{H^s}\|{\rm d}\|_{H^s}^4\ , \end{eqnarray}
(3.5)

其中 C 是一个常数, 它依赖于 s, n, \lambda_1, \lambda_2 .注意到 |{\rm d}| = 1 蕴含着 \|{\rm d}\|_{H^s}>0 , 因此由(3.5)式我们可以得到

\frac{\rm d}{{\rm d}t}\|{\rm d}(t, \cdot)\|_{H^s} \leq C\|u\|_{H^{s+1}}\|{\rm d}(t, \cdot) \|_{H^s}+C\|u\|_{H^{s+1}}\|{\rm d}(t, \cdot)\|_{H^s}^3.

对于给定的 u\in L^\infty\big([0, T_1];H^{s+1}({\mathbb T}^n) \big) , 结合Grownwall's不等式, 我们可以得到

\begin{equation} \|{\rm d}(t, \cdot)\|_{H^s} \leq \frac{\|{\rm d}_0\|_{H^s}e^{\int_0^t C\|u(\tau, \cdot)\|_{H^{s+1}} {\rm d}\tau}}{\sqrt{1-\|{\rm d}_0\|^2_{H^s}\int_0^t 2C\|u(\tau, \cdot)\|_{H^{s+1}}e^{\int_0^ \tau 2C\|u(\eta, \cdot)\|_{H^s}{\rm d}\eta} {\rm d}\tau}}\triangleq M(t)\ , \end{equation}
(3.6)

这里 0<t<T_2\leq T_1 , T_2 依赖于 \|{\rm d}_0\|_{H^s}, \|u\|_{H^{s+1}} 和常数 C .

第二步  逼近解的构造.

由于对给定的 u 而言, 方程(3.3)是拟线性的, 因此我们可以用如下的抛物型方程来逼近它

\begin{equation} \left\{ \begin{array}{ll} { } \partial_t {\rm d}^\nu +u\cdot\nabla{\rm d}^\nu +{\rm B}{\rm d}^\nu+\frac{1}{\lambda_1}\gamma^\nu{\rm d}^ \nu+\frac{\lambda_2}{\lambda_1}{\rm A}{\rm d}^\nu-\nu\Delta{\rm d}^\nu = 0, \\ {\rm d}^\nu\big|_{t = 0} = {\rm d}_0\ , \end{array} \right. \end{equation}
(3.7)

其中 \gamma^\nu = -\lambda_2{{\rm d}^\nu}^\top{\rm A}{\rm d}^\nu.

根据抛物方程的经典理论(参见文献[7]), 对任意的 \nu>0 , 方程(3.7)存在唯一的经典解

{\rm d}^\nu\in L^\infty\big([0, T_2];H^s({\mathbb T}^n) \big)\cap L^2\big([0, T_2];H^{s+1}({\mathbb T}^n)\big),

其中 T_2 依赖于 \|{\rm d}_0\|_{H^s}, \|u\|_{H^{s+1}} .并且, 由方程(3.4)–(3.5)的推导过程, 我们可以得出如下能量估计

\begin{equation} \frac{1}{2}\frac{\rm d}{{\rm d}t}\|{\rm d}^\nu(t, \cdot)\|_{H^s}^2+\nu\int_0^t\|\nabla {\rm d}^\nu(s, \cdot)\|_{H^s}^2{\rm d}s\leq C\|u\|_{H^{s+1}}\|{\rm d}^\nu\|_{H^s}^2+C\|u\|_{H^{s+1}}\|{\rm d}^\nu\|_{H^s}^4. \end{equation}
(3.8)

根据方程(3.6), 对任意 \nu>0 , 对任意的 t\in [0, T_2] , 我们有 \|{\rm d}^\nu(t, \cdot)\|_{H^s}\leq M(t) .对方程(3.8)两边关于 t 0 T_2 积分, 我们可以得到

\begin{eqnarray} &&\nu \int_0^{T_2} \|\nabla {\rm d}^\nu(s, \cdot)\|_{H^s}^2{\rm d}s {}\\ &\leq &\int_0^{T_2} \bigg[C\|u(t, \cdot)\|_{H^{s+1}}M(t){\rm d}t+C\|u(t, \cdot)\|_{H^{s+1}}M^3(t)\bigg]{\rm d}t+\frac{1}{2}\|{\rm d}_0\|_{H^s}^2 \triangleq \tilde{M}_{T_2}. \end{eqnarray}
(3.9)

并且, 我们还能推导出解序列 \{{\rm d}^\nu\}_{\nu>0} 在空间 C\big([0, T_2];H^s({\mathbb T}^n) \big) 中是一致有界的.

第三步  收敛性.

在这一步中, 我们将证明解序列 \{{\rm d}^\nu\}_{\nu>0} 是空间 C\big([0, T];H^{s-1}({\mathbb T}^n) \big) 中的Cauchy列.令 0<\nu_1<\nu_2 . {\rm d}^{\nu_1} {\rm d}^{\nu_2} 分别为方程

\begin{equation} \left\{ \begin{array}{ll} { } \partial_t{\rm d}^{\nu_1} +u\cdot\nabla{\rm d}^{\nu_1}+{\rm B}{\rm d}^{\nu_1}+\frac{1}{\lambda_1}\gamma^{\nu_1}{\rm d}^{\nu_1}+\frac{\lambda_2}{\lambda_1}{\rm A}{\rm d}^{\nu_1}-\nu_1\Delta {\rm d}^{\nu_1} = 0, \\ {\rm d}^{\nu_1}\big|_{t = 0} = {\rm d}_0 \end{array} \right. \end{equation}
(3.10)

\begin{equation} \left\{ \begin{array}{ll} { } \partial_t{\rm d}^{\nu_2} +u\cdot\nabla{\rm d}^{\nu_2}+{\rm B}{\rm d}^{\nu_2}+\frac{1}{\lambda_1}\gamma^{\nu_2}{\rm d}^{\nu_2}+\frac{\lambda_2}{\lambda_1}{\rm A}{\rm d}^{\nu_2}-\nu_2\Delta {\rm d}^{\nu_2} = 0, \\ {\rm d}^{\nu_2}\big|_{t = 0} = {\rm d}_0 \end{array} \right. \end{equation}
(3.11)

的经典解, 其中 \gamma^{\nu_1} \gamma^{\nu_2} 分别表示为

\gamma^{\nu_1} = -\lambda_2 {{\rm d}^{\nu_1}}^\top {\rm A} {\rm d}^{\nu_1}, \quad \gamma^{\nu_2} = -\lambda_2 {{\rm d}^{\nu_2}}^\top {\rm A}{\rm d}^{\nu_2}.

由方程(3.10)和(3.11), 可知 {\rm d}^{\nu_2}-{\rm d}^{\nu_1} 满足如下方程

\begin{eqnarray} &&\partial_t({\rm d}^{\nu_2}-{\rm d}^{\nu_1})+u\cdot\nabla({\rm d}^{\nu_2}-{\rm d}^{\nu_1})+{\rm B}({\rm d}^{\nu_2}-{\rm d}^{\nu_1})+\frac{\lambda_2}{\lambda_1}{\rm A}({\rm d}^{\nu_2}-{\rm d}^{\nu_1}){}\\ &&+\frac{1}{\lambda_1}(\gamma^{\nu_2}{\rm d}^{\nu_2}-\gamma^{\nu_1}{\rm d}^{\nu_1})-\nu_2\Delta({\rm d}^{\nu_2}-{\rm d}^{\nu_1}) = (\nu_2-\nu_1)\Delta {\rm d}^{\nu_1}, \end{eqnarray}
(3.12)

相应的初始条件为

\begin{equation} {\rm d}^{\nu_2}-{\rm d}^{\nu_1}\big|_{t = 0} = 0. \end{equation}
(3.13)

这里, \gamma^{\nu_2}{\rm d}^{\nu_2}-\gamma^{\nu_1}{\rm d}^{\nu_1} 可以表示成

\begin{eqnarray} \gamma^{\nu_2}{\rm d}^{\nu_2}-\gamma^{\nu_1}{\rm d}^{\nu_1} & = &-\lambda_2\big\{[ ({\rm d}^{\nu_2})^\top {\rm A} {\rm d}^{\nu_2}]{\rm d}^{\nu_2}-\big[({\rm d}^{\nu_1})^\top {\rm A} {\rm d}^{\nu_1} \big]{\rm d}^{\nu_1} \big\} {}\\ & = &-\lambda_2 \big\{ [({\rm d}^{\nu_2})^\top {\rm A} {\rm d}^{\nu_2}]({\rm d}^{\nu_2}-{\rm d}^{\nu_1})+({\rm d}^{\nu_2})^\top {\rm A}({\rm d}^{\nu_2}-{\rm d}^{\nu_1}){\rm d}^{\nu_1}{} \\ &&+[({\rm d}^{\nu_2}-{\rm d}^{\nu_1})^\top {\rm A} {\rm d}^{\nu_1}]{\rm d}^{\nu_1} \big\}.{} \end{eqnarray}

将方程(3.12)的两边同时与 {\rm d}^{\nu_2}-{\rm d}^{\nu_1} H^{s-1} -内积, 我们可以得到

\begin{eqnarray} &&\frac{1}{2}\frac{\rm d}{{\rm d}t}\|{\rm d}^{\nu_2}-{\rm d}^{\nu_1}\|_{H^{s-1}}^2 +\nu_2\|\nabla({\rm d}^{\nu_2}-{\rm d}^{\nu_1})\|_{H^{s-1}}^2 {}\\ &\leq& C\|u\|_{H^{s-1}}\|{\rm d}^{\nu_2}-{\rm d}^{\nu_1}\|_{H^{s-1}}^2+C\|\nabla u\|_{H^{s-1}}\|{\rm d}^{\nu_2}-{\rm d}^{\nu_1}\|_{H^{s-1}}^2 {}\\ && +\frac{\lambda_2}{\lambda_1}C\|\nabla u\|_{H^{s-1}}\|{\rm d}^{\nu_2}-{\rm d}^{\nu_1}\|_{H^{s-1}}^2+\frac{\lambda_2}{\lambda_1}C\|{\rm d}^{\nu_2}\|_{H^{s-1}}^2\|\nabla u\|_{H^{s-1}}\|{\rm d}^{\nu_2}-{\rm d}^{\nu_1}\|_{H^{s-1}}^2 \\ && +\frac{\lambda_2}{\lambda_1}C\|{\rm d}^{\nu_2}\|_{H^{s-1}}\|{\rm d}^{\nu_1}\|_{H^{s-1}}\|\nabla u\|_{H^{s-1}}\|{\rm d}^{\nu_2}-{\rm d}^{\nu_1}\|_{H^{s-1}}^2 {}\\ && +\frac{\lambda_2}{\lambda_1}C\|{\rm d}^{\nu_1}\|_{H^{s-1}}^2\|\nabla u\|_{H^{s-1}}\|{\rm d}^{\nu_2}-{\rm d}^{\nu_1}\|_{H^{s-1}}^2 +(\nu_2-\nu_1)\|\Delta {\rm d}^{\nu_1}\|_{H^{s-1}}\|{\rm d}^{\nu_2}-{\rm d}^{\nu_1}\|_{H^{s-1}}\\ &\leq& C\|u\|_{H^s}\|{\rm d}^{\nu_2}-{\rm d}^{\nu_1}\|_{H^{s-1}}^2+C\|u\|_{H^s}M(t)^2\|{\rm d}^{\nu_2}-{\rm d}^{\nu_1}\|_{H^{s-1}}^2 {}\\ && +(\nu_2-\nu_1)\|\Delta {\rm d}^{\nu_1}\|_{H^{s-1}}\|{\rm d}^{\nu_2}-{\rm d}^{\nu_1}\|_{H^{s-1}}. \end{eqnarray}
(3.14)

由方程(3.14), 我们可以得到

\frac{\rm d}{{\rm d}t}\|{\rm d}^{\nu_2}-{\rm d}^{\nu_1}\|_{H^{s-1}} \leq C\|u\|_{H^{s-1}}(1+M(t)^2)\|{\rm d}^{\nu_2}-{\rm d}^{\nu_1}\|_{H^{s-1}}+(\nu_2-\nu_1)\|\Delta {\rm d}^{\nu_1}\|_{H^{s-1}}.

由Grönwall's不等式以及方程(3.9), 我们有

\begin{eqnarray*} \|{\rm d}^{\nu_2}-{\rm d}^{\nu_1}(t, \cdot)\|_{H^{s-1}} &\leq& e^{\int_0^t C\|u(s, \cdot)\|_{H^s}(1+M(s)^2){\rm d}s}\int_0^t (\nu_2-\nu_1)\|\nabla {\rm d}^{\nu_1}(s, \cdot)\|_{H^{s}} {\rm d}s {\nonumber}\\ & \leq &e^{\int_0^t C\|u(s, \cdot)\|_{H^s}(1+M(s)^2){\rm d}s} \bigg(\int_0^t\nu_2\|\nabla {\rm d}^{\nu_1}\|_{H^s}^2{\rm d}s \bigg)^{1\over2}\sqrt{\nu_2-\nu_1}t^{1\over2}.{\nonumber}\\ &\leq &e^{\int_0^t C\|u(s, \cdot)\|_{H^s}(1+M(s)^2){\rm d}s}\sqrt{\tilde{M}_{T_2}}\sqrt{\nu_2-\nu_1}t^{\frac{1}{2}}.{\nonumber} \end{eqnarray*}

这意味着 \{{\rm d}^\nu\}_{\nu>0}\; \mbox{为}\; C\big([0, T_2];H^{s-1} \big)\; \mbox{中的Cauchy序列}. 因此, 随着 \nu\to0 在方程(3.8)中取极限, 我们得到对于给定的光滑 u 以及初始条件 {\rm d}\big|_{t = 0} = {\rm d}_0 下方程(3.3)的解 {\rm d} .

现在我们证明了, 如果 {\rm d}_0\in H^s({\mathbb T}^n), u\in L^\infty\big([0, T_1];H^{s+1}({\mathbb T}^n)\big)\cap L^1\big([0, T_1];H^{s+1}({\mathbb T}^n)\big) , 方程(3.3)存在一个解 {\rm d}\in L^\infty\big([0, T_2];H^s({\mathbb T}^n)\big). 然而我们并不知道 {\rm d} 是否满足几何约束条件 |{\rm d}(t, x)| = 1 .

在第二个部分, 我们将证明, 如果初始方向场 {\rm d}_0 满足几何约束条件 |{\rm d}_0| = 1 , 则当 0<t<T_2 时, 方程的解 {\rm d} 也满足几何约束条件.由于 {\rm d} 是方程

\begin{equation} \partial_t {\rm d}+u\cdot\nabla {\rm d}+{\rm Bd}+\frac{1}{\lambda_1}\gamma {\rm d}+\frac{\lambda_2}{\lambda_1}{\rm Ad} = 0 \end{equation}
(3.15)

的经典解, 将方程(3.15)两边同时乘以 {\rm d} , 我们得到

\begin{equation} \frac{1}{2}\partial_t|{\rm d}|^2+\frac{1}{2}u\cdot\nabla |{\rm d}|^2+\frac{1}{\lambda_1}\gamma|{\rm d}|^2+\frac{\lambda_2}{\lambda_1}{\rm d}^\top {\rm Ad} = 0. \end{equation}
(3.16)

\gamma 的表达式代入(3.16)式, 我们得到 |{\rm d}|^2-1 满足如下方程

\begin{equation} \left\{ \begin{array}{ll} { } \partial_t(|{\rm d}|^2-1)+u\cdot\nabla(|{\rm d}|^2-1)+\frac{1}{\lambda_1}\gamma(|{\rm d}|^2-1) = 0, \\ |{\rm d}|^2-1\big|_{t = 0} = 0. \end{array} \right. \end{equation}
(3.17)

解这个关于 |{\rm d}^2|-1 的输运方程(3.17), 容易验证如果 |{\rm d}_0| = 1 , 则当 0<t<T_2 时, |{\rm d}| = 1 也成立.如果令 T = \min (T_1, T_2) , 则定理得证.

4 大Ericksen数和Reynolds数极限

在这一节, 我们回到方程(1.11), 即

\begin{equation} \left\{ \begin{array}{ll} \partial_t u^ \varepsilon+u^ \varepsilon\cdot\nabla u^ \varepsilon+\nabla p^ \varepsilon = \varepsilon^2 {\rm{div}}\Sigma^ \varepsilon_1+ \varepsilon^3 {\rm{div}}\Sigma^ \varepsilon_2+ \varepsilon^2 {\rm{div}}\Sigma^ \varepsilon_3, \\ {\rm{divu}} ^ \varepsilon = 0, \\ \varepsilon \ddot{\rm d}^ \varepsilon = \varepsilon\Delta {\rm d}^ \varepsilon+\Gamma^ \varepsilon {\rm d}^ \varepsilon+\lambda_1(\dot{\rm d}^ \varepsilon+{\rm B}^ \varepsilon {\rm d}^ \varepsilon)+\lambda_2 {\rm A}^ \varepsilon{\rm d}^ \varepsilon, \end{array} \right. \end{equation}
(4.1)

其中 {\rm A}^ \varepsilon, {\rm B}^ \varepsilon, \Sigma_1^ \varepsilon, \Sigma_2^ \varepsilon, \Sigma_3^ \varepsilon, \dot{\rm d}^ \varepsilon, \ddot{\rm d}^ \varepsilon 形式如同在(1.2)–(1.5)式将 u, {\rm d} 换成 u^ \varepsilon, {\rm d}^ \varepsilon .此时, \Gamma^ \varepsilon 具有如下形式

\begin{equation} \Gamma^ \varepsilon = - \varepsilon |\dot{\rm d}^ \varepsilon|^2+ \varepsilon|\nabla {\rm d}^ \varepsilon|^2-\lambda_2{{\rm d}^ \varepsilon}^\top {\rm A}^ \varepsilon {\rm d}^ \varepsilon. \end{equation}
(4.2)

在前面一节, 我们已经证明了, 当 s>\frac{n}{2}+1 , 对于 u_0\in H^{s+1}({\mathbb T}^n) {\rm d}_0\in H^s({\mathbb T}^n) , 极限方程存在唯一解 (u, {\rm d}) 使得

u\in L^\infty([0, T];H^{s+1}), \quad {\rm d}\in L^\infty([0, T];H^s),

其中 T>0 依赖于 n, s, \lambda_1, \lambda_2 和初始值 (u_0, {\rm d}_0) .现在我们将要研究方程(4.1)的经典解 (u^ \varepsilon, {\rm d}^ \varepsilon, \dot{\rm d}^ \varepsilon) 的渐近行为.第一步我们需要给出一个合适的初始条件.事实上, 由方程(1.14)的第三个方程可知

\dot{\rm d} = -\frac{1}{\lambda_1}(\gamma {\rm d}+\lambda_2 {\rm Ad})-{\rm Bd}.

因此, 方程(4.1)可以给出如下初始条件

\begin{equation} u^ \varepsilon\big|_{t = 0} = u_0, \quad {\rm d}^ \varepsilon\big|_{t = 0} = {\rm d}_0, \quad \dot{\rm d}^ \varepsilon\big|_{t = 0} = \dot{\rm d}\big|_{t = 0}, \end{equation}
(4.3)

其中 \dot{\rm d}\big|_{t = 0} = \big\{-\frac{1}{\lambda_1}(\gamma {\rm d}+\lambda_2 {\rm Ad})-{\rm Bd}\big\}\big|_{t = 0} .

我们将验证方程(4.1)–(4.3)的大Ericksen数和Reynolds数极限, 也即方程的解 (u^ \varepsilon, {\rm d}^ \varepsilon, \dot{\rm d}^ \varepsilon) 随着 \varepsilon\to 0 的极限.为此, 我们需要得出方程解 ({u}^ \varepsilon, {\rm d}^ \varepsilon, \dot{\rm d}^ \varepsilon) 的一个一致估计.由文献[19]中所推导的基本能量守恒律, 我们有如下引理.

引理4.1  假设 (u^ \varepsilon, {\rm d}^ \varepsilon, \dot{\rm d}^ \varepsilon) 是方程(4.1)–(4.3)的一个光滑解, 则下面的关系成立

\begin{eqnarray} && {1\over2} \frac{\rm d}{{\rm d}t}\big(\|u^ \varepsilon\|_{L^2}^2+ \varepsilon^3\|\dot{\rm d}^ \varepsilon\|_{L^2}^2+ \varepsilon^3\|\nabla {\rm d}^ \varepsilon\|_{L^2}^2 \big)+{1\over2} \varepsilon^2\mu_4\|\nabla u^ \varepsilon\|_{L^2}^2+\mu_1 \varepsilon^2\|{{\rm d}^ \varepsilon}^\top {\rm A}^ \varepsilon {\rm d}^ \varepsilon\|_{L^2}^2{}\\ &&-\lambda_1 \varepsilon^2\|\dot{\rm d}^ \varepsilon+{\rm B}^ \varepsilon {\rm d}^ \varepsilon\|_{L^2}^2-2\lambda_2 \varepsilon^2\langle \dot{\rm d}^ \varepsilon+{\rm B}^ \varepsilon {\rm d}^ \varepsilon, {\rm A}^ \varepsilon {\rm d}^ \varepsilon\rangle+(\mu_5+\mu_6) \varepsilon^2\|{\rm A}^ \varepsilon {\rm d}^ \varepsilon\|_{L^2}^2 = 0. \end{eqnarray}
(4.4)

如果, 更进一步, \varepsilon>0, \lambda_1<0, \mu_1\geq0, \mu_4>0 \mu_5+\mu_6+\frac{\lambda_2^2}{\lambda_1}\geq0 , 则基本能量守恒律(4.4)是耗散的.

  对于任意的 \varepsilon>0 . (u^ \varepsilon, {\rm d}^ \varepsilon, \dot{\rm d}^ \varepsilon) 为方程(4.1)–(4.3)的经典解.将方程(4.1)的第一个方程的两边同时与 u^ \varepsilon L^2 -内积, 我们得到

\begin{equation} {1\over 2}\frac{\rm d}{{\rm d}t}\|u^ \varepsilon\|_{L^2}^2+{1\over2} \varepsilon^2\mu_4\|\nabla u^ \varepsilon\|_{L^2}^2 = \varepsilon^3\langle- {\rm{div}}(\nabla {\rm d}^ \varepsilon\odot\nabla {\rm d}^ \varepsilon), u^ \varepsilon\rangle+ \varepsilon^2\langle \partial_j\tilde{\sigma}_{ji}, u^ \varepsilon_i\rangle. \end{equation}
(4.5)

将方程(4.1)的第三个方程的两边同时与 \varepsilon^2\dot{\rm d}^ \varepsilon L^2 -内积, 我们得到

\begin{eqnarray} && {1\over2}\frac{\rm d}{{\rm d}t}\big( \varepsilon^3\|\dot{\rm d}^ \varepsilon\|_{L^2}^2 + \varepsilon^3\|\nabla {\rm d}^ \varepsilon\|_{L^2}^2 \big)-\lambda_1 \varepsilon^2\|\dot{\rm d}^ \varepsilon\|_{L^2}^2 {}\\ & = & \varepsilon^3\langle\Delta {\rm d}^ \varepsilon, u^ \varepsilon\cdot\nabla {\rm d}^ \varepsilon\rangle+\lambda_1 \varepsilon^2\langle {\rm B}^ \varepsilon {\rm d}^ \varepsilon, \dot{\rm d}^ \varepsilon\rangle+\lambda_2 \varepsilon^2\langle {\rm A}^ \varepsilon {\rm d}^ \varepsilon, \dot{\rm d}^ \varepsilon\rangle, \end{eqnarray}
(4.6)

其中我们用到关系 {\rm d}^ \varepsilon\cdot \dot{\rm d}^ \varepsilon = 0 .将(4.5)式和(4.6)式加起来, 结合文献[19, Lemma 2.1]中的计算步骤, 我们可以得到下面的基本能量守恒律

\begin{eqnarray*} && {1\over2}\frac{\rm d}{{\rm d}t}\big(\|u^ \varepsilon\|_{L^2}^2+ \varepsilon^3\|\dot{\rm d}^ \varepsilon\|_{L^2}^2+ \varepsilon^3\|\nabla {\rm d}^ \varepsilon\|_{L^2}^2 \big)+{1\over2} \varepsilon^2\mu_4\|\nabla u^ \varepsilon\|_{L^2}^2+\mu_1 \varepsilon^2\|{{\rm d}^ \varepsilon}^\top {\rm A}^ \varepsilon {\rm d}^ \varepsilon\|_{L^2}^2\\ &&-\lambda_1 \varepsilon^2\|\dot{\rm d}^ \varepsilon+{\rm B}^ \varepsilon {\rm d}^ \varepsilon\|_{L^2}^2-2\lambda_2 \varepsilon^2\langle \dot{\rm d}^ \varepsilon+{\rm B}^ \varepsilon {\rm d}^ \varepsilon, {\rm A}^ \varepsilon {\rm d}^ \varepsilon\rangle+(\mu_5+\mu_6) \varepsilon^2\|{\rm A}^ \varepsilon {\rm d}^ \varepsilon\|_{L^2}^2 = 0, \end{eqnarray*}

其中用到抵消关系 \varepsilon^3\langle - {\rm{div}}(\nabla {\rm d}^ \varepsilon\odot\nabla {\rm d}^ \varepsilon), u^ \varepsilon\rangle+ \varepsilon^3\langle\Delta {\rm d}^ \varepsilon, u^ \varepsilon\cdot\nabla {\rm d}^ \varepsilon\rangle = 0. 引理4.1得证.

现在我们可以阐述我们关于双曲型液晶方程的大Ericksen数极限的主要结果.

定理4.1  令 s>\frac{n}{2}+1 . u_0\in H^{s+2}({\mathbb T}^n) 满足 \nabla \cdot u_0 = 0 , {\rm d}_0\in H^{s+1}({\mathbb T}^n) 满足 |{\rm d}_0| = 1 . u\in C\big([0, T_*];H^{s+2}({\mathbb T}^n) \big), {\rm d}\in C\big([0, T_*];H^{s+1}({\mathbb T}^n) \big) 为极限方程(3.1)在初始值为 (u_0, {\rm d}_0) 时的局部经典解.

对任意的 \varepsilon>0 . (u^ \varepsilon, {\rm d}^ \varepsilon, \dot{\rm d}^ \varepsilon) 为方程(4.1)在初始条件(4.3)下的局部经典解, 使得

u^ \varepsilon\in L^\infty\big([0, T_ \varepsilon];H^s({\mathbb T}^n)), \ \ \nabla {\rm d}^ \varepsilon\in L^\infty\big([0, T_ \varepsilon];H^s({\mathbb T}^n)\big), \ \ \dot{\rm d}^ \varepsilon\in L^\infty\big([0, T_ \varepsilon];H^s({\mathbb T}^n)\big),

对某个 T_ \varepsilon>0 .

如果我们记 \omega^ \varepsilon = u^ \varepsilon-u, \theta^ \varepsilon = {\rm d}^ \varepsilon-{\rm d} 以及 \dot{\theta}^ \varepsilon = \dot{\rm d}^ \varepsilon-\dot{\rm d} , 则我们有

\begin{equation} \begin{array}{ll} { } \|\omega^ \varepsilon(t, \cdot)\|_{L^2}^2+ \varepsilon^3\|\nabla\theta^ \varepsilon(t, \cdot)\|_{L^2}^2+ \varepsilon^3\|\dot\theta^ \varepsilon(t, \cdot)\|_{L^2}^2\leq C_T \varepsilon, \\ { } \int_0^T \big( \varepsilon^2\|\nabla\omega^ \varepsilon(t, \cdot)\|_{L^2}^2+ \varepsilon^2\|{{\rm d}^ \varepsilon}^\top({\rm A}^ \varepsilon-{\rm A}){\rm d}^ \varepsilon\|_{L^2}^2 \big){\rm d}t\leq C_T \varepsilon, \ \ \forall\ t\in [0, T], \end{array} \end{equation}
(4.7)

其中 T = \min{(T_*, T_ \varepsilon)} C_T 依赖于 s, n , 和Leslie系数以及初始值 (u_0, {\rm d}_0) .

  形式上, 对方程(4.1)和方程(3.1)作差之后, 我们得到

\begin{equation} \left\{ \begin{array}{ll} { } \partial_t (u^ \varepsilon-u) +u^ \varepsilon\cdot\nabla(u^ \varepsilon-u) +(u^ \varepsilon-u)\cdot\nabla u+\nabla(p^ \varepsilon-p)-\frac{1}{2}\mu_4 \varepsilon^2\Delta(u^ \varepsilon-u)\\ { } = \frac{1}{2}\mu_4 \varepsilon^2\Delta u- \varepsilon^3 {\rm{div}}(\nabla {\rm d}^ \varepsilon\odot\nabla {\rm d}^ \varepsilon)+ \varepsilon^2 {\rm{div}}\tilde{\sigma}, \\ { } \varepsilon\partial_t \dot{\rm d}^ \varepsilon+ \varepsilon u^ \varepsilon\cdot\nabla\dot{\rm d}^ \varepsilon = \varepsilon\Delta {\rm d}^ \varepsilon +(\Gamma^ \varepsilon {\rm d}^ \varepsilon-\gamma {\rm d})+\lambda_1(\dot{\rm d}^ \varepsilon-\dot{\rm d})+\lambda_1({\rm B}^ \varepsilon {\rm d}^ \varepsilon-{\rm Bd})\\ { } +\lambda_2({\rm A}^ \varepsilon {\rm d}^ \varepsilon-{\rm Ad}). \end{array} \right. \end{equation}
(4.8)

将方程(4.8)的两边同时与 u^ \varepsilon-u 相乘, 并在 {\mathbb T}^n 上积分, 我们得到

\begin{eqnarray} &&\frac{1}{2}\frac{\rm d}{{\rm d}t} \|\omega^ \varepsilon\|_{L^2}^2 +\langle\omega^ \varepsilon\cdot\nabla u, \omega^ \varepsilon\rangle +\frac{1}{2}\mu_4 \varepsilon^2\|\nabla\omega^ \varepsilon\|_{L^2}^2 {}\\ & = &\frac{1}{2}\mu_4 \varepsilon^2\langle\Delta u, \omega^ \varepsilon\rangle- \varepsilon^3\langle {\rm{div}}(\nabla {\rm d}^ \varepsilon\odot\nabla {\rm d}^ \varepsilon), \omega^ \varepsilon\rangle+ \varepsilon^2\langle {\rm{div}}\tilde{\sigma}, \omega^ \varepsilon\rangle, \end{eqnarray}
(4.9)

其中我们用到这样的事实:由 {\rm{div}}(u^ \varepsilon-u) = 0 可得 \langle u^ \varepsilon\cdot\nabla\omega^ \varepsilon+\nabla(p^ \varepsilon-p), \omega^ \varepsilon\rangle = 0 . {\rm \theta}^ \varepsilon = {\rm d}^ \varepsilon-{\rm d} \dot{\rm \theta}^ \varepsilon = \dot{\rm d}^ \varepsilon-\dot{\rm d} .注意到

\dot{\theta}^ \varepsilon = \partial_t({\rm d}^ \varepsilon-{\rm d})+(u^ \varepsilon\cdot\nabla {\rm d}^ \varepsilon-u\cdot\nabla {\rm d}) = \partial_t\theta^ \varepsilon+u^ \varepsilon\cdot\nabla\theta^ \varepsilon+\omega^ \varepsilon\cdot\nabla {\rm d}.

\ddot{\rm d} = \partial_t\dot{\rm d}+u\cdot\nabla\dot{\rm d} .将方程(4.8)的第二个方程的两边同时减去 \varepsilon\ddot{\rm d} , 将它改写成如下形式

\begin{eqnarray} \varepsilon\partial_t \dot{\theta}^ \varepsilon+ \varepsilon u^ \varepsilon\cdot\nabla\dot{\theta}^ \varepsilon+ \varepsilon \omega^ \varepsilon\cdot\nabla\dot{\rm d} & = &- \varepsilon\partial_t\dot{\rm d}- \varepsilon u\cdot\nabla \dot{\rm d}+ \varepsilon\Delta {\rm d}^ \varepsilon+(\Gamma^ \varepsilon {\rm d}^ \varepsilon-\gamma {\rm d}) \\ &&+\lambda_1(\dot{\rm d}^ \varepsilon -\dot{\rm d})+\lambda_1({\rm B}^ \varepsilon {\rm d}^ \varepsilon-{\rm Bd})+\lambda_2({\rm A}^ \varepsilon {\rm d}^ \varepsilon-{\rm Ad}). \end{eqnarray}
(4.10)

将方程(4.10)两边同时乘以 \varepsilon^2\dot{\theta}^ \varepsilon 并在 {\mathbb T}^n 上积分, 我们得到

\begin{eqnarray} &&\frac{1}{2}\frac{\rm d}{{\rm d}t}\big( \varepsilon^3\|\dot{\theta}^ \varepsilon\|_{L^2}^2+ \varepsilon^3\|\nabla\theta^ \varepsilon\|_{L^2}^2 \big) + \varepsilon^3\langle\omega^ \varepsilon\cdot\nabla\dot{\rm d}, \dot{\theta}^ \varepsilon\rangle{}\\ & = &- \varepsilon^3\langle \partial_t\dot{\rm d}+u\cdot\nabla\dot{\rm d}, \dot{\theta}^ \varepsilon\rangle + \varepsilon^3\langle \Delta {\rm d}^ \varepsilon, \omega^ \varepsilon\cdot\nabla {\rm d}^ \varepsilon+u\cdot\nabla\theta^ \varepsilon\rangle + \varepsilon^2\langle \Gamma^ \varepsilon {\rm d}^ \varepsilon-\gamma {\rm d}, \dot{\theta}^ \varepsilon\rangle \\ &&+\lambda_1 \varepsilon^2\|\dot{\theta}^ \varepsilon\|_{L^2}^2+\lambda_1 \varepsilon^2\langle {\rm B}^ \varepsilon {\rm d}^ \varepsilon -{\rm Bd}, \dot{\theta}^ \varepsilon\rangle+\lambda_2 \varepsilon^2\langle {\rm A}^ \varepsilon {\rm d}^ \varepsilon-{\rm Ad}, \dot{\theta}^ \varepsilon\rangle + \varepsilon^3\langle \Delta {\rm d}, \partial_t\theta^\varepsilon \rangle . \end{eqnarray}
(4.11)

将(4.9)式和(4.11)式加起来, 我们得到

\begin{eqnarray} &&\frac{1}{2}\frac{\rm d}{{\rm d}t}\big(\|\omega^ \varepsilon\|_{L^2}^2+ \varepsilon^3\|\dot{\theta}^ \varepsilon\|_{L^2}^2+ \varepsilon^3\|\nabla\theta^ \varepsilon\|_{L^2}^2 \big)+\frac{1}{2}\mu_4 \varepsilon^2 \|\nabla\omega^ \varepsilon\|_{L^2}^2{}\\ & = &\frac{1}{2}\mu_4 \varepsilon^2\langle \Delta u, \omega^ \varepsilon\rangle + \varepsilon^2\langle {\rm{div}}\tilde{\sigma}, \omega^ \varepsilon\rangle-\langle\omega^ \varepsilon\cdot\nabla u, \omega^ \varepsilon\rangle- \varepsilon^3\langle\omega^ \varepsilon\cdot\nabla\dot{\rm d}, \dot{\theta}^ \varepsilon\rangle- \varepsilon^3 \langle \partial_t\dot{\rm d}+u\cdot\nabla\dot{\rm d}, \dot{\theta}^ \varepsilon\rangle {}\\ &&+ \varepsilon^3\langle \Delta {\rm d}^ \varepsilon, u\cdot\nabla\theta^ \varepsilon\rangle+ \varepsilon^2\langle \Gamma^ \varepsilon {\rm d}^ \varepsilon-\gamma {\rm d}, \dot{\theta}^ \varepsilon\rangle+\lambda_1 \varepsilon^2\|\dot{\theta}^ \varepsilon\|_{L^2}^2+\lambda_1 \varepsilon^2\langle {\rm B}^ \varepsilon {\rm d}^ \varepsilon-{\rm Bd}, \dot{\theta}^ \varepsilon\rangle {}\\ &&+\lambda_2 \varepsilon^2\langle {\rm A}^ \varepsilon {\rm d}^ \varepsilon-{\rm Ad}, \dot{\theta}^ \varepsilon\rangle + \varepsilon^3\langle \Delta {\rm d}, \partial_t\theta^\varepsilon \rangle, \end{eqnarray}
(4.12)

这里我们用到事实: - \varepsilon^3\langle {\rm{div}}(\nabla {\rm d}^ \varepsilon\odot\nabla {\rm d}^ \varepsilon), \omega^ \varepsilon\rangle+ \varepsilon^3\langle\Delta {\rm d}^ \varepsilon, \omega^ \varepsilon\cdot\nabla {\rm d}^ \varepsilon\rangle = 0 .剩下的工作就是去估计(4.12)式右边的各项.

首先, \varepsilon^3\langle \Delta {\rm d}^ \varepsilon, u\cdot\nabla \theta^ \varepsilon\rangle 这项可如下估计

\begin{eqnarray} \varepsilon^3 \langle \Delta {\rm d}^ \varepsilon, u\cdot\nabla \theta^ \varepsilon\rangle & = & \varepsilon^3 \langle\Delta\theta^ \varepsilon, u\cdot\nabla\theta^ \varepsilon\rangle + \varepsilon^3\langle\Delta {\rm d}, u\cdot\nabla\theta^ \varepsilon\rangle {}\\ & = &- \varepsilon^3\langle \nabla \theta^ \varepsilon, \nabla u\cdot\nabla \theta^ \varepsilon\rangle + \varepsilon^3\langle\Delta {\rm d}, u\cdot\nabla \theta^ \varepsilon\rangle \end{eqnarray}
(4.13)

其次, 我们将 \lambda_1 \varepsilon^2\|\dot{\theta}^ \varepsilon\|_{L^2}^2 这项分拆为如下形式

\begin{eqnarray} \lambda_1 \varepsilon^2\|\dot{\theta}^ \varepsilon\|_{L^2}^2 = \lambda_1 \varepsilon^2\|\dot{\rm d}^ \varepsilon\|_{L^2}^2+\lambda_1 \varepsilon^2\|\dot{\rm d}\|_{L^2}^2-2\lambda_1 \varepsilon^2\langle\dot{\rm d}^ \varepsilon, \dot{\rm d}\rangle. \end{eqnarray}
(4.14)

于是, 在估计 \varepsilon^2\langle {\rm{div}}\tilde{\sigma}, \omega^ \varepsilon\rangle 这项之前, 我们先回顾 \tilde{\sigma} 具有如下形式

\begin{eqnarray} \big(\tilde\sigma(u^ \varepsilon, {\rm d}^ \varepsilon, \dot{\rm d}^ \varepsilon) \big)_{ij} & = &\mu_1{\rm d}^ \varepsilon_k{\rm d}^ \varepsilon_p{\rm A}^ \varepsilon_{kp}{\rm d}^ \varepsilon_i{\rm d}^ \varepsilon_j+\mu_2{\rm d}^ \varepsilon_i(\dot {\rm d}^ \varepsilon_j+{\rm B}^ \varepsilon_{jk}{\rm d}^ \varepsilon_k)+\mu_3{\rm d}^ \varepsilon_j(\dot {\rm d}^ \varepsilon_i+{\rm B}^ \varepsilon_{ik}{\rm d}^ \varepsilon_k) {}\\ &&+\mu_5 {\rm d}^ \varepsilon_i {\rm A}^ \varepsilon_{jk}{\rm d}^ \varepsilon_k+\mu_6{\rm d}^ \varepsilon_j {\rm A}^ \varepsilon_{ik}{\rm d}^ \varepsilon_k\ . \end{eqnarray}
(4.15)

为了陈述的简单, 记

\begin{equation} {\rm I} = \varepsilon^2\langle {\rm{div}}\tilde{\sigma}, \omega^ \varepsilon\rangle = {\rm I}_1+{\rm I}_2+\cdots +{\rm I}_4, \end{equation}
(4.16)

其中 {\rm I}_1, {\rm I}_2, \cdots , {\rm I}_4 是(4.15)式的右边部分表达式.例如, {\rm I}_1 具有如下形式

\begin{eqnarray} {\rm I}_1 & = & \varepsilon^2\langle \partial_j(\mu_1 {\rm d}_k^ \varepsilon {\rm d}_p^ \varepsilon {\rm A}^ \varepsilon_{kp} {\rm d}_i^ \varepsilon {\rm d}_j^ \varepsilon), \omega^ \varepsilon_i\rangle {}\\ & = & \varepsilon^2 \langle \partial_j(\mu_1 {\rm d}_k^ \varepsilon {\rm d}_p^ \varepsilon ({\rm A}^ \varepsilon_{kp}-{\rm A}_{kp}) {\rm d}_i^ \varepsilon {\rm d}_j^ \varepsilon), \omega^ \varepsilon_i\rangle + \varepsilon^2\langle \partial_j(\mu_1 {\rm d}_k^ \varepsilon {\rm d}_p^ \varepsilon {\rm A}_{kp} {\rm d}_i^ \varepsilon {\rm d}_j^ \varepsilon), \omega^ \varepsilon_i\rangle {}\\ & = &-\mu_1 \varepsilon^2\|{{\rm d}^ \varepsilon}^\top({\rm A}^ \varepsilon-{\rm A}){\rm d}^ \varepsilon\|_{L^2}^2+ \varepsilon^2\langle \partial_j(\mu_1 {\rm d}_k^ \varepsilon {\rm d}_p^ \varepsilon {\rm A}_{kp} {\rm d}_i^ \varepsilon {\rm d}_j^ \varepsilon), \omega^ \varepsilon_i\rangle. \end{eqnarray}
(4.17)

{\rm I}_2 具有如下形式

\begin{eqnarray} {\rm I}_2 & = & \varepsilon^2\langle \partial_j(\mu_2 {\rm d}_j^ \varepsilon {\rm B}_{ki}^ \varepsilon {\rm d}_k^ \varepsilon+\mu_3 {\rm d}_i^ \varepsilon {\rm B}_{kj}^ \varepsilon {\rm d}_k^ \varepsilon), \omega^ \varepsilon_i\rangle {}\\ & = & \varepsilon^2 \langle \partial_j(\mu_2 {\rm d}_j^ \varepsilon ({\rm B}_{ki}^ \varepsilon-{\rm B}_{ki}) {\rm d}_k^ \varepsilon+\mu_3 {\rm d}_i^ \varepsilon ({\rm B}_{kj}^ \varepsilon-{\rm B}_{kj}) {\rm d}_k^ \varepsilon), \omega^ \varepsilon_i\rangle {}\\ &&+ \varepsilon^2\langle \partial_j(\mu_2 {\rm d}_j^ \varepsilon {\rm B}_{ki} {\rm d}_k^ \varepsilon+\mu_3 {\rm d}_i^ \varepsilon {\rm B}_{kj} {\rm d}_k^ \varepsilon), \omega^ \varepsilon_i\rangle {} \\ & = &\lambda_1 \varepsilon^2\|({\rm B}^ \varepsilon-{\rm B})d^ \varepsilon\|_{L^2}^2+\lambda_2 \varepsilon^2\langle ({\rm B}^ \varepsilon-{\rm B}){\rm d}^ \varepsilon, ({\rm A}^ \varepsilon-{\rm A})d^ \varepsilon\rangle {}\\ &&+ \varepsilon^2\langle \partial_j(\mu_2 {\rm d}_j^ \varepsilon {\rm B}_{ki} {\rm d}_k^ \varepsilon+\mu_3 {\rm d}_i^ \varepsilon {\rm B}_{kj} {\rm d}_k^ \varepsilon), \omega^ \varepsilon_i\rangle. \end{eqnarray}
(4.18)

{\rm I}_3 具有如下形式

\begin{eqnarray} {\rm I}_3 & = & \varepsilon^2\langle \partial_j(\mu_5 {\rm d}_j^ \varepsilon {\rm d}_k^ \varepsilon {\rm A}_{ki}^ \varepsilon+\mu_6 {\rm d}_i^ \varepsilon {\rm d}_k^ \varepsilon {\rm A}_{kj}^ \varepsilon ), \omega^ \varepsilon_i\rangle {}\\ & = & \varepsilon^2\langle \partial_j(\mu_5 {\rm d}_j^ \varepsilon {\rm d}_k^ \varepsilon ({\rm A}_{ki}^ \varepsilon-{\rm A}_{ki})+\mu_6 {\rm d}_i^ \varepsilon {\rm d}_k^ \varepsilon ({\rm A}_{kj}^ \varepsilon-{\rm A}_{kj}) ), \omega^ \varepsilon_i\rangle {}\\ &&+ \varepsilon^2\langle \partial_j(\mu_5 {\rm d}_j^ \varepsilon {\rm d}_k^ \varepsilon {\rm A}_{ki}+\mu_6 {\rm d}_i^ \varepsilon {\rm d}_k^ \varepsilon {\rm A}_{kj} ), \omega^ \varepsilon_i\rangle {}\\ & = &-(\mu_5+\mu_6) \varepsilon^2\|({\rm A}^ \varepsilon-{\rm A}){\rm d}^ \varepsilon\|_{L^2}^2+\lambda_2 \varepsilon^2\langle ({\rm A}^ \varepsilon-{\rm A}){\rm d}^ \varepsilon, ({\rm B}^ \varepsilon-{\rm B}){\rm d}^ \varepsilon\rangle {}\\ &&+ \varepsilon^2\langle \partial_j(\mu_5 {\rm d}_j^ \varepsilon {\rm d}_k^ \varepsilon {\rm A}_{ki}+\mu_6 {\rm d}_i^ \varepsilon {\rm d}_k^ \varepsilon {\rm A}_{kj} ), \omega^ \varepsilon_i\rangle. \end{eqnarray}
(4.19)

最后, {\rm I}_4 具有如下形式

\begin{eqnarray} {\rm I}_4 & = & \varepsilon^2 \langle \partial_j(\mu_2 {\rm d}_j^ \varepsilon \dot{\rm d}_i^ \varepsilon+\mu_3{\rm d}_j^ \varepsilon\dot{\rm d}_i^ \varepsilon), \omega^ \varepsilon_i\rangle {}\\ & = &\lambda_2 \varepsilon^2\langle \dot{\rm d}^ \varepsilon, ({\rm A}^ \varepsilon-{\rm A}){\rm d}^ \varepsilon\rangle+\lambda_1 \varepsilon^2\langle \dot{\rm d}^ \varepsilon, ({\rm B}^ \varepsilon-{\rm B}){\rm d}^ \varepsilon\rangle. \end{eqnarray}
(4.20)

然后我们记 {\rm II} = \lambda_1 \varepsilon^2\langle {\rm B}^ \varepsilon {\rm d}^ \varepsilon-{\rm Bd}, \dot{\theta}^ \varepsilon\rangle , 它可以改写为如下形式

\begin{eqnarray} {\rm II} & = &\lambda_1 \varepsilon^2\langle {\rm B}^ \varepsilon {\rm d}^ \varepsilon-{\rm B d} , \dot{\rm d}^ \varepsilon-\dot{\rm d} \rangle {}\\ & = &\lambda_1 \varepsilon^2\langle ({\rm B}^ \varepsilon-{\rm B}){\rm d}^ \varepsilon, \dot{\rm d}^ \varepsilon \rangle +\lambda_1 \varepsilon^2\langle {\rm B}\theta^ \varepsilon, \dot{\theta}^ \varepsilon \rangle-\lambda_1 \varepsilon^2 \langle ({\rm B}^ \varepsilon-{\rm B}){\rm d}^ \varepsilon, \dot{\rm d}\rangle. \end{eqnarray}
(4.21)

{\rm III} = \lambda_2 \varepsilon^2 \langle {\rm A}^ \varepsilon {\rm d}^ \varepsilon-{\rm A d}, \dot{\theta}^ \varepsilon\rangle , 并将它改写为

\begin{eqnarray} {\rm III} & = &\lambda_2 \varepsilon^2 \langle {\rm A}^ \varepsilon {\rm d}^ \varepsilon-{\rm A d}, \dot{\rm d}^ \varepsilon-\dot{\rm d} \rangle {}\\ & = & \lambda_2 \varepsilon^2 \langle ({\rm A}^ \varepsilon-{\rm A}){\rm d}^ \varepsilon+{\rm A}({\rm d}^ \varepsilon-{\rm d}), \dot{\rm d}^ \varepsilon-\dot{\rm d}\rangle {}\\ & = &\lambda_2 \varepsilon^2 \langle ({\rm A}^ \varepsilon-{\rm A}){\rm d}^ \varepsilon, \dot{\rm d}^ \varepsilon\rangle -\lambda_2 \varepsilon^2 \langle ({\rm A}^ \varepsilon-{\rm A}){\rm d}^ \varepsilon, \dot{\rm d}\rangle+\lambda_2 \varepsilon^2\langle {\rm A}({\rm d}^ \varepsilon-{\rm d}), \dot{\theta}^ \varepsilon\rangle. \end{eqnarray}
(4.22)

{\rm IV} = \varepsilon^2\langle \Gamma^ \varepsilon {\rm d}^ \varepsilon-\gamma {\rm d}, \dot{\rm d}^ \varepsilon-\dot{\rm d} \rangle , 回顾 \Gamma^ \varepsilon = - \varepsilon|\dot{\rm d}^ \varepsilon|^2+ \varepsilon|\nabla {\rm d}^ \varepsilon|^2-\lambda_2 {{\rm d}^ \varepsilon}^\top {\rm A}^ \varepsilon {\rm d}^ \varepsilon \gamma = -\lambda_2 {\rm d}^\top{\rm Ad} .注意到 |{\rm d}^ \varepsilon| = 1 , 于是 \dot{\rm d}^ \varepsilon\cdot {\rm d}^ \varepsilon = 0. 因此, {\rm IV} 可以改写为

\begin{eqnarray} {\rm IV} & = & \varepsilon^2 \langle \Gamma^ \varepsilon {\rm d}^ \varepsilon, -\dot{\rm d} \rangle - \varepsilon^2 \langle \Gamma {\rm d}, \dot{\rm d}^ \varepsilon-\dot{\rm d}\rangle {}\\ & = & \varepsilon^3\langle |\dot{\rm d}^ \varepsilon|^2 {\rm d}^ \varepsilon, \dot{\rm d}\rangle - \varepsilon^3\langle |\nabla {\rm d}^ \varepsilon|^2 {\rm d}^ \varepsilon, \dot{\rm d}\rangle +\lambda_2 \varepsilon^2\langle ({{\rm d}^ \varepsilon}^\top {\rm A}^ \varepsilon {\rm d}^ \varepsilon){\rm d}^ \varepsilon, \dot{\rm d}\rangle- \varepsilon^2 \langle \gamma {\rm d}, \dot{\rm d}^ \varepsilon-\dot{\rm d}\rangle. \end{eqnarray}
(4.23)

现在将(4.13)–(4.23)式代入(4.12)式, 我们得到

\begin{eqnarray} &&\frac{1}{2} \frac{\rm d}{{\rm d}t}\big(\|\omega^ \varepsilon\|_{L^2}^2+ \varepsilon^3\|\dot{\theta}^ \varepsilon\|_{L^2}^2+ \varepsilon^3\|\nabla \theta^ \varepsilon\|_{L^2}^2 \big)+\frac{1}{2}\mu_4 \varepsilon^2\|\nabla\omega^ \varepsilon\|_{L^2}^2+\mu_1 \varepsilon^2\|{{\rm d}^ \varepsilon}^\top({\rm A}^ \varepsilon-{\rm A}){\rm d}^ \varepsilon\|_{L^2}^2 {}\\ &&-\lambda_1 \varepsilon^2\|\dot{\rm d}^ \varepsilon+({\rm B}^ \varepsilon-{\rm B}){\rm d}^ \varepsilon \|_{L^2}^2 -2\lambda_2 \varepsilon^2 \langle \dot{\rm d}^ \varepsilon+({\rm B}^ \varepsilon-{\rm B}){\rm d}^ \varepsilon, ({\rm A}^ \varepsilon-{\rm A}){\rm d}^ \varepsilon\rangle {}\\ && +(\mu_5+\mu_6) \varepsilon^2\|({\rm A}^ \varepsilon-{\rm A}){\rm d}^ \varepsilon\|_{L^2}^2{} \\ & = &\frac{1}{2}\mu_4 \varepsilon^2\langle\Delta u, \omega^ \varepsilon\rangle-\langle\omega^ \varepsilon\cdot\nabla u, \omega^ \varepsilon\rangle- \varepsilon^3\langle\omega^ \varepsilon\cdot\nabla\dot{\rm d}, \dot{\theta}^ \varepsilon\rangle- \varepsilon^3\langle\ddot{\rm d}, \dot{\theta}^ \varepsilon\rangle+ \varepsilon^3\langle\Delta {\rm d}, u\cdot\nabla\theta^ \varepsilon\rangle {}\\ &&+ \varepsilon^3\langle\Delta\theta^ \varepsilon, u\cdot\nabla\theta^ \varepsilon\rangle+ \varepsilon^2\langle \partial_j(\mu_1 {\rm d}_k^ \varepsilon {\rm d}_p^ \varepsilon {\rm A}_{kp} {\rm d}_i^ \varepsilon {\rm d}_j^ \varepsilon), \omega^ \varepsilon_i\rangle {}\\ &&+ \varepsilon^2\langle \partial_j(\mu_2 {\rm d}_j^ \varepsilon {\rm B}_{ki} {\rm d}_k^ \varepsilon +\mu_3 {\rm d}_i^ \varepsilon {\rm B}_{kj} {\rm d}_k^ \varepsilon), \omega^ \varepsilon_i\rangle+ \varepsilon^2\langle \partial_j(\mu_5 {\rm d}_j^ \varepsilon {\rm d}_k^ \varepsilon {\rm A}_{ki}+\mu_6{\rm d}_i^ \varepsilon {\rm d}_k^ \varepsilon {\rm A}_{kj} ), \omega^ \varepsilon_i\rangle {}\\ &&+\lambda_1 \varepsilon^2\langle {\rm B}\theta^ \varepsilon, \dot{\theta}^ \varepsilon \rangle-\lambda_1 \varepsilon^2 \langle ({\rm B}^ \varepsilon-{\rm B}){\rm d}^ \varepsilon, \dot{\rm d}\rangle-\lambda_2 \varepsilon^2 \langle ({\rm A}^ \varepsilon-{\rm A}){\rm d}^ \varepsilon, \dot{\rm d}\rangle+\lambda_2 \varepsilon^2\langle {\rm A}({\rm d}^ \varepsilon-{\rm d}), \dot{\theta}^ \varepsilon\rangle {}\\ &&+ \varepsilon^3\langle |\dot{\rm d}^ \varepsilon|^2 {\rm d}^ \varepsilon, \dot{\rm d}\rangle - \varepsilon^3\langle |\nabla {\rm d}^ \varepsilon|^2 {\rm d}^ \varepsilon, \dot{\rm d}\rangle +\lambda_2 \varepsilon^2\langle ({{\rm d}^ \varepsilon}^\top {\rm A}^ \varepsilon {\rm d}^ \varepsilon){\rm d}^ \varepsilon, \dot{\rm d}\rangle{}\\ &&- \varepsilon^2 \langle \gamma {\rm d}, \dot{\rm d}^ \varepsilon-\dot{\rm d}\rangle + \varepsilon^3\langle \Delta {\rm d}, \partial_t\theta^\varepsilon \rangle. \end{eqnarray}
(4.24)

剩下的工作是估计(4.24)式的右边各项.为了陈述的简单, 我们记 {\rm R}_1 为如下形式

\begin{eqnarray} {\rm R}_1& = &\frac{1}{2}\mu_4 \varepsilon^2\langle\Delta u, \omega^ \varepsilon\rangle-\langle\omega^ \varepsilon\cdot\nabla u, \omega^ \varepsilon\rangle- \varepsilon^3\langle\omega^ \varepsilon\cdot\nabla\dot{\rm d}, \dot{\theta}^ \varepsilon\rangle- \varepsilon^3\langle\ddot{\rm d}, \dot{\theta}^ \varepsilon\rangle{}\\ &&+ \varepsilon^3\langle\Delta {\rm d}, u\cdot\nabla\theta^ \varepsilon\rangle + \varepsilon^3\langle \Delta {\rm d}, \partial_t\theta^\varepsilon \rangle, \end{eqnarray}
(4.25)

其中 \ddot{\rm d} = \partial_t\dot{\rm d}+u\cdot\nabla\dot{\rm d} .注意到 (u, {\rm d}) 是极限方程(1.14)的一个经典解, 于是

\dot{\rm d} = -\frac{1}{\lambda_1}\gamma {\rm d}-\frac{\lambda_2}{\lambda_1}{\rm Ad}-{\rm Bd}.

因此, \ddot{\rm d} = \partial_t(-\frac{1}{\lambda_1}\gamma {\rm d}-\frac{\lambda_2}{\lambda_1}{\rm Ad}-{\rm Bd})+u\cdot\nabla(-\frac{1}{\lambda_1}\gamma {\rm d}-\frac{\lambda_2}{\lambda_1}{\rm Ad}-{\rm Bd}), 同时注意到 \partial_t u\in L^\infty([0, T];H^{s+1}), \partial_t{\rm d}\in L^\infty([0, T];H^{s-1}) , 我们可以得到

\ddot{\rm d}\in L^\infty([0, T];H^{s-1}).

于是, 由Hölder不等式和Cauchy不等式, 我们有

\begin{eqnarray} |{\rm R}_1|& \leq &\frac{1}{2}\mu_4 \varepsilon^2\|\Delta u\|_{L^2}\|\omega^ \varepsilon\|_{L^2}+\|\nabla u\|_{L^\infty}\|\omega^ \varepsilon\|_{L^2}^2+ \varepsilon^3\|\nabla\dot{\rm d}\|_{L^\infty}\|\omega^ \varepsilon\|_{L^2}\|\dot{\theta}^ \varepsilon\|_{L^2}\\ &&+ \varepsilon^3\|\ddot{\rm d}\|_{L^2}\|\dot{\theta}^ \varepsilon\|_{L^2}+ \varepsilon^3\|\Delta {\rm d}\|_{L^2}\|u\|_{L^\infty} \|\nabla\theta^ \varepsilon\|_{L^2} + \varepsilon^3\|\Delta {\rm d}\|_{L^2}\|\dot{\theta}^ \varepsilon\|_{L^2} \\ &&+ \varepsilon^3\|\Delta {\rm d}\|_{L^\infty} \|\omega^ \varepsilon\|_{L^2}\|\nabla \theta^ \varepsilon\|_{L^\infty} + \varepsilon^3\|\Delta {\rm d}\|_{L^\infty} \|\omega^ \varepsilon\|_{L^2}\|\nabla {\rm d}\|_{L^2} \\ &\leq& \mu_4 \varepsilon^2\big(\|\Delta u\|_{L^2}^2+\|\omega^ \varepsilon\|_{L^2}^2 \big)+\|\nabla u\|_{L^\infty}\|\omega^ \varepsilon\|_{L^2}^2+\frac{1}{2} \varepsilon^3\|\nabla\dot{\rm d}\|_{L^\infty}\big(\|\omega^ \varepsilon\|_{L^2}^2+\|\dot{\theta}^ \varepsilon\|_{L^2}^2 \big)\\ &&+\frac{1}{2} \varepsilon^3\big(\|\ddot{\rm d}\|_{L^2}^2+\|\dot{\theta}^ \varepsilon\|_{L^2}^2 \big)+\frac{1}{2} \varepsilon^3\big(\|\Delta {\rm d}\|_{L^2}^2\|u\|_{L^\infty}^2+\|\nabla\theta^ \varepsilon\|_{L^2}^2 \big) +\frac{1}{2} \varepsilon^3\big(\|\Delta{\rm d}\|_{L^2}^2+\|\nabla{\rm d}\|_{L^2}^2){}\\ && +\frac{1}{2} \varepsilon^3\big(\|\dot{\theta}^ \varepsilon\|_{L^2}^2+\|\nabla\theta^ \varepsilon\|_{L^2}^2 \big) +\frac{1}{2}(1+ \varepsilon^2) \|\Delta{\rm d}\|_{L^\infty}^2\|\omega^ \varepsilon\|_{L^2}^2. \end{eqnarray}
(4.26)

{\rm R}_2

\begin{eqnarray} {\rm R}_2 & = & \varepsilon^3\langle \Delta\theta^ \varepsilon, u\cdot\nabla\theta^ \varepsilon\rangle = \varepsilon^3\int_{{\mathbb T}^n}\sum\limits_{i, j, k = 1}^n\partial_j^2 \theta^ \varepsilon_i (u_k\partial_k\theta^ \varepsilon_i) {\rm d}x \\ & = &- \varepsilon^3\int_{{\mathbb T}^n}\sum\limits_{i, j, k = 1}^n\partial_j\theta^ \varepsilon_i (\partial_j u_k\partial_k\theta^ \varepsilon_i) {\rm d}x. \end{eqnarray}
(4.27)

于是, 由Hölder不等式, 我们有

\begin{equation} |{\rm R}_2| \leq \varepsilon^3\|\nabla u\|_{L^\infty}\|\nabla \theta^ \varepsilon\|_{L^2}^2. \end{equation}
(4.28)

{\rm R}_3

\begin{eqnarray} {\rm R}_3 & = & \varepsilon^2\langle \partial_j(\mu_1 {\rm d}_k^ \varepsilon {\rm d}_p^ \varepsilon {\rm A}_{kp} {\rm d}_i^ \varepsilon {\rm d}_j^ \varepsilon), \omega^ \varepsilon_i\rangle+ \varepsilon^2\langle \partial_j(\mu_2 {\rm d}_j^ \varepsilon {\rm B}_{ki} {\rm d}_k^ \varepsilon +\mu_3 {\rm d}_i^ \varepsilon {\rm B}_{kj} {\rm d}_k^ \varepsilon), \omega^ \varepsilon_i\rangle \\ &&+ \varepsilon^2\langle \partial_j(\mu_5 {\rm d}_j^ \varepsilon {\rm d}_k^ \varepsilon {\rm A}_{ki}+\mu_6{\rm d}_i^ \varepsilon {\rm d}_k^ \varepsilon {\rm A}_{kj} ), \omega^ \varepsilon_i\rangle \\ & = &- \varepsilon^2\langle (\mu_1 {\rm d}_k^ \varepsilon {\rm d}_p^ \varepsilon {\rm A}_{kp} {\rm d}_i^ \varepsilon {\rm d}_j^ \varepsilon), \partial_j\omega^ \varepsilon_i\rangle- \varepsilon^2\langle (\mu_2 {\rm d}_j^ \varepsilon {\rm B}_{ki} {\rm d}_k^ \varepsilon +\mu_3 {\rm d}_i^ \varepsilon {\rm B}_{kj} {\rm d}_k^ \varepsilon), \partial_j\omega^ \varepsilon_i\rangle \\ &&- \varepsilon^2\langle (\mu_5 {\rm d}_j^ \varepsilon {\rm d}_k^ \varepsilon {\rm A}_{ki}+\mu_6{\rm d}_i^ \varepsilon {\rm d}_k^ \varepsilon {\rm A}_{kj} ), \partial_j\omega^ \varepsilon_i\rangle . \end{eqnarray}
(4.29)

由事实 |{\rm d}^ \varepsilon| = 1 , 我们有

\begin{eqnarray} |{\rm R}_3| &\leq &(|\mu_1|+|\mu_2|+|\mu_3|+|\mu_5|+|\mu_6|) \varepsilon^2\|\nabla u\|_{L^2}\|\nabla\omega^ \varepsilon\|_{L^2} \\ &\leq &\frac{C}{\mu_4} \varepsilon^2\|\nabla u\|_{L^2}^2+\frac{\mu_4}{8} \varepsilon^2\|\nabla\omega^ \varepsilon\|_{L^2}^2, \end{eqnarray}
(4.30)

其中 C 依赖于 |\mu_1|, |\mu_2|, |\mu_3|, |\mu_5| and |\mu_6| . {\rm R}_4 为如下形式

\begin{equation} {\rm R}_4 = \lambda_1 \varepsilon^2\langle {\rm B}\theta^ \varepsilon, \dot{\theta}^ \varepsilon \rangle-\lambda_1 \varepsilon^2 \langle ({\rm B}^ \varepsilon-{\rm B}){\rm d}^ \varepsilon, \dot{\rm d}\rangle-\lambda_2 \varepsilon^2 \langle ({\rm A}^ \varepsilon-{\rm A}){\rm d}^ \varepsilon, \dot{\rm d}\rangle +\lambda_2 \varepsilon^2\langle {\rm A}({\rm d}^ \varepsilon-{\rm d}), \dot{\theta}^ \varepsilon\rangle. \end{equation}
(4.31)

注意到 |{\rm d}| = 1, |{\rm d}^ \varepsilon| = 1 , 于是 |\theta^ \varepsilon| = |{\rm d}^ \varepsilon-{\rm d}|\leq2 .因此, 我们有

\begin{eqnarray} |{\rm R}_4| &\leq& C|\lambda_1| \varepsilon^2\|\nabla u\|_{L^2}\|\dot{\theta}^ \varepsilon\|_{L^2}+|\lambda_1|C \varepsilon^2\|\nabla\omega^ \varepsilon\|_{L^2}\|\dot{\rm d}\|_{L^2}+|\lambda_2|C \varepsilon^2\|\nabla \omega^ \varepsilon\|_{L^2}\|\dot{\rm d}\|_{L^2}\\ &&+|\lambda_2|C \varepsilon^2\|\nabla u\|_{L^2}\|\dot{\theta}^ \varepsilon\|_{L^2} \\ &\leq& C \varepsilon^3\|\dot{\theta}^ \varepsilon\|_{L^2}^2+C \varepsilon\|\nabla u\|_{L^2}^2+\frac{\mu_4}{8} \varepsilon^2\|\nabla\omega^ \varepsilon\|_{L^2}^2+\frac{C \varepsilon^2}{\mu_4}\|\dot{\rm d}\|_{L^2}^2, \end{eqnarray}
(4.32)

其中 C 依赖于 |\lambda_1| and |\lambda_2| . {\rm R}_5

\begin{eqnarray} {\rm R}_5 & = & \varepsilon^3\langle |\dot{\rm d}^ \varepsilon|^2 {\rm d}^ \varepsilon, \dot{\rm d}\rangle - \varepsilon^3\langle |\nabla {\rm d}^ \varepsilon|^2 {\rm d}^ \varepsilon, \dot{\rm d}\rangle +\lambda_2 \varepsilon^2\langle ({{\rm d}^ \varepsilon}^\top {\rm A}^ \varepsilon {\rm d}^ \varepsilon){\rm d}^ \varepsilon, \dot{\rm d}\rangle- \varepsilon^2 \langle \Gamma {\rm d}, \dot{\rm d}^ \varepsilon-\dot{\rm d}\rangle \\ &\triangleq &R_{51}+R_{52}+R_{53}+R_{54}. \end{eqnarray}
(4.33)

于是

\begin{equation} {\rm R}_{51} = \varepsilon^3 \langle |\dot{\rm d}^ \varepsilon-\dot{\rm d}+\dot{\rm d}|^2 {\rm d}^ \varepsilon, \dot{\rm d}\rangle = \varepsilon^3 \langle |\dot{\rm d}^ \varepsilon-\dot{\rm d}|^2 {\rm d}^ \varepsilon+2[(\dot{\rm d}^ \varepsilon-\dot{\rm d})\cdot \dot{\rm d}]{\rm d}^ \varepsilon+|\dot{\rm d}|^2 {\rm d}^ \varepsilon , \dot{\rm d}\rangle. \end{equation}
(4.34)

因此

\begin{eqnarray} |{\rm R}_{51}| &\leq& \varepsilon^3\|\dot{\rm d}\|_{L^\infty}\|\dot{\theta}^ \varepsilon\|_{L^2}^2+2 \varepsilon^3\|\dot{\theta}^ \varepsilon\|_{L^2}\|\dot{\rm d}\|_{L^4}^2+ \varepsilon^3\|\dot{\rm d}\|_{L^3}^3 {}\\ &\leq& \varepsilon^3\|\dot{\rm d}\|_{L^\infty}\|\dot{\theta}^ \varepsilon\|_{L^2}^2+ \varepsilon^3\|\dot{\theta}^ \varepsilon\|_{L^2}^2+ \varepsilon^3\|\dot{\rm d}\|_{L^4}^4+ \varepsilon^3\|\dot{\rm d}\|_{L^3}^3. \end{eqnarray}
(4.35)

由类似的方式, {\rm R}_{52} 可以改写为

\begin{eqnarray} {\rm R}_{52} & = &- \varepsilon^3\langle |\nabla {\rm d}^ \varepsilon-\nabla {\rm d}+\nabla {\rm d}|^2 {\rm d}^ \varepsilon, \dot{\rm d}\rangle {}\\ & = &- \varepsilon^3\langle |\nabla {\rm d}^ \varepsilon-\nabla {\rm d}|^2 {\rm d}^ \varepsilon+2[(\nabla {\rm d}^ \varepsilon-\nabla {\rm d})\cdot\nabla {\rm d}]{\rm d}^ \varepsilon+|\nabla {\rm d}|^2 {\rm d}^ \varepsilon, \dot{\rm d}\rangle. \end{eqnarray}
(4.36)

这推出

\begin{eqnarray} |{\rm R}_{52}| &\leq& \varepsilon^3 \|\dot{\rm d}\|_{L^\infty} \|\nabla\theta^ \varepsilon\|_{L^2}^2+2 \varepsilon^3\|\nabla {\rm d}\|_{L^4}\|\dot{\rm d}\|_{L^4}\|\nabla\theta^ \varepsilon\|_{L^2}+ \varepsilon^3\|\nabla {\rm d}\|_{L^4}^2 \|\dot{\rm d}\|_{L^2} {}\\ &\leq & \varepsilon^3 \|\dot{\rm d}\|_{L^\infty} \|\nabla\theta^ \varepsilon\|_{L^2}^2+ \varepsilon^3\|\nabla\theta^ \varepsilon\|_{L^2}^2+ \varepsilon^3\|\nabla {\rm d}\|_{L^4}^2\|\dot{\rm d}\|_{L^4}^2+ \varepsilon^3\|\nabla {\rm d}\|_{L^4}^2\|\dot{\rm d}\|_{L^2}. \end{eqnarray}
(4.37)

因此 {\rm R}_{53} 可改写为

\begin{eqnarray} {\rm R}_{53} & = &\lambda_2 \varepsilon^2 \langle [{{\rm d}^ \varepsilon}^\top({\rm A}^ \varepsilon-{\rm A+A}){\rm d}^ \varepsilon]{\rm d}^ \varepsilon, \dot {\rm d}\rangle {}\\ & = &\lambda_2 \varepsilon^2\langle [{{\rm d}^ \varepsilon}^\top({\rm A}^ \varepsilon-{\rm A}){\rm d}^ \varepsilon]{\rm d}^ \varepsilon, \dot{\rm d}\rangle+\lambda_2 \varepsilon^2\langle ({{\rm d}^ \varepsilon}^\top {\rm A}{\rm d}^ \varepsilon){\rm d}^ \varepsilon, \dot{\rm d}\rangle. \end{eqnarray}
(4.38)

于是

\begin{eqnarray} |{\rm R}_{53}| &\leq &C\lambda_2 \varepsilon^2\|\nabla\omega^ \varepsilon\|_{L^2}\|\dot{\rm d}\|_{L^2}+C\lambda_2 \varepsilon^2\|\nabla u\|_{L^2}\|\dot{\rm d}\|_{L^2} {}\\ &\leq &\frac{\mu_4}{8} \varepsilon^2\|\nabla\omega^ \varepsilon\|_{L^2}^2+\frac{C}{\mu_4}\lambda_2^2 \varepsilon^2\|\dot{\rm d}\|_{L^2}^2+C\lambda_2 \varepsilon^2\|\nabla u\|_{L^2}\|\dot{\rm d}\|_{L^2}. \end{eqnarray}
(4.39)

回顾 \gamma = -\lambda_2 {\rm d}^\top {\rm A}{\rm d} , 于是

\begin{eqnarray} {\rm R}_{54} = \lambda_2 \varepsilon^2\langle ({\rm d}^\top {\rm A}{\rm d}){\rm d}, \dot{\rm d}^ \varepsilon-\dot{\rm d}\rangle. \end{eqnarray}
(4.40)

|{\rm d}| = 1 和Cauchy-Schwartz's不等式, 我们有

\begin{eqnarray} |{\rm R}_{54}| \leq C\lambda_2 \varepsilon^2\|\nabla u\|_{L^2}\|\dot{\theta}^ \varepsilon\|_{L^2} \leq C\lambda_2^2 \varepsilon\|\nabla u\|_{L^2}^2+ \varepsilon^3\|\dot{\theta}^ \varepsilon\|_{L^2}^2. \end{eqnarray}
(4.41)

将估计(4.35), (4.37), (4.39)以及(4.41)代入(4.33)式, 我们有

\begin{eqnarray} |{\rm R}_5| &\leq &(2+\|\dot{\rm d}\|_{L^\infty})\big( \varepsilon^3\|\dot{\theta}^ \varepsilon\|_{L^2}^2+ \varepsilon^3\|\nabla \theta^ \varepsilon\|_{L^2}^2\big)+ \varepsilon^3\big( \|\dot{\rm d}\|_{L^4}^4+\|\dot{\rm d}\|_{L^3}^3+\|\nabla{\rm d}\|_{L^4}^2\|\dot{\rm d}\|_{L^4}^2 \\ &&+\|\nabla {\rm d}\|_{L^4}^2\|\dot{\rm d}\|_{L^2}\big)+\frac{\mu_4}{8} \varepsilon^2\|\nabla\omega^ \varepsilon\|_{L^2}^2+ \varepsilon^2\big(\frac{C}{\mu_4}\lambda_2^2\|\dot{\rm d}\|_{L^2}^2+C\lambda_2\|\nabla u\|_{L^2}\|\dot{\rm d}\|_{L^2} \big) \\ &&+C\lambda_2^2 \varepsilon\|\nabla u\|_{L^2}^2, \end{eqnarray}
(4.42)

其中 C 依赖于 \lambda_2 .将估计(4.26), (4.28), (4.30), (4.32)以及(4.42)代入(4.24)式, 可得

\begin{equation} \frac{\rm d}{{\rm d}t}{\cal Y}(t)+{\cal D}(t) \leq C{\cal H}(t){\cal Y}(t)+C \varepsilon{\cal M}_1(t)+C \varepsilon^2{\cal M}_2(t)+C \varepsilon^3{\cal M}_3(t), \end{equation}
(4.43)

其中

\begin{eqnarray*} {\cal Y}(t) & = &\|\omega^ \varepsilon\|_{L^2}^2+ \varepsilon^3\|\nabla\theta^ \varepsilon\|_{L^2}^2+ \varepsilon^3\|\dot{\theta}^ \varepsilon\|_{L^2}^2, {\nonumber}\\ {\cal H}(t) & = &1+ \varepsilon^2+\|\nabla u\|_{L^\infty}+\|\dot{\rm d}\|_{L^\infty}+\|\nabla\dot{\rm d}\|_{L^\infty} + \varepsilon^3\|\nabla\dot{\rm d}\|_{L^\infty} +(1+\varepsilon^3)\|\Delta{\rm d}\|_{L^\infty}^2 , {\nonumber}\\ {\cal D}(t) & = &\frac{\mu_4}{8} \varepsilon^2\|\nabla\omega^ \varepsilon\|_{L^2}^2+\mu_1 \varepsilon^2\|{{\rm d}^ \varepsilon}^\top({\rm A}^ \varepsilon-{\rm A}){\rm d}^ \varepsilon\|_{L^2}^2-\lambda_1 \varepsilon^2\|\dot{\rm d}^ \varepsilon+({\rm B}^ \varepsilon-{\rm B}){\rm d}^ \varepsilon \|_{L^2}^2 \nonumber\\ && -2\lambda_2 \varepsilon^2 \langle \dot{\rm d}^ \varepsilon+({\rm B}^ \varepsilon-{\rm B}){\rm d}^ \varepsilon, ({\rm A}^ \varepsilon-{\rm A}){\rm d}^ \varepsilon\rangle+(\mu_5+\mu_6) \varepsilon^2\|({\rm A}^ \varepsilon-{\rm A}){\rm d}^ \varepsilon\|_{L^2}^2, \nonumber\\ {\cal M}_1(t) & = &(1+\lambda_2^2)\|\nabla u\|_{L^2}^2, {\nonumber}\\ {\cal M}_2(t) & = &\frac{1+\lambda_2^2}{\mu_4}\|\dot{\rm d}\|_{L^2}^2+\frac{1}{\mu_4}\|\nabla u\|_{L^2}^2+\lambda_2\|\nabla u\|_{L^2}\|\dot{\rm d}\|_{L^2}+\mu_4\|\Delta u\|_{L^2}^2, {\nonumber}\\ {\cal M}_3(t) & = &\|\Delta {\rm d}\|_{L^2}^2\|u\|_{L^\infty}^2+\|\ddot{\rm d}\|_{L^2}^2 +\|\dot{\rm d}\|_{L^4}^4+\|\dot{\rm d}\|_{L^3}^3+\|\nabla {\rm d}\|_{L^4}^2\|\dot{\rm d}\|_{L^4}^2 \\ && +\|\nabla{\rm d}\|_{L^4}^2\|\dot{\rm d}\|_{L^2} +\|\Delta{\rm d}\|^2_{L^2}+\|\nabla{\rm d}\|^2_{L^2}, {\nonumber} \end{eqnarray*}

其中 C 是一个常数, 它依赖于 s, n 和Leslie系数, 但是它不依赖于 \varepsilon .由于 {\cal Y}(0) = 0 , 由Gronwall's不等式, 可得对所有的 0< \varepsilon<1

\begin{equation} {\cal Y}(t)\leq C_T \varepsilon, \quad \forall\ 0<t<T, \end{equation}
(4.44)

其中 C_T 依赖于 s, n , Leslie系数以及极限方程的经典解 (u, {\rm d}) .由(4.44)式的上界, 将它对 t 0 T 积分, 我们得到(4.7)式.定理4.1得证.

注4.1  根据文献[19], 双曲液晶方程的局部经典解 (u^ \varepsilon, {\rm d}^ \varepsilon, \dot{\rm d}^ \varepsilon) 的存在时间是依赖于初始能量和Leslie系数的.然而, 对于我们这里的参数化的双曲液晶方程(4.1)–(4.3), 它的初始能量和Leslie系数都依赖于参数 \varepsilon .特别地, 参数化的双曲液晶方程的Leslie系数参数化为 \mu_1 \varepsilon, \mu_2 \varepsilon, \cdots , \mu_6 \varepsilon .随着 \varepsilon 趋向于0, 它们变得非常小.因此我们并不知道, 对于所有的 0< \varepsilon\leq1 , 方程(4.1)–(4.3)的局部经典解 (u^ \varepsilon, {\rm d}^ \varepsilon, \dot{\rm d}^ \varepsilon) 的存在时间是否有一个一致的下界.因此, 定理4.1仅仅是关于双曲液晶方程的大Ericksen数的极限的一个形式分析的结果.

注4.2  通过标准的能量方法, 似乎我们不能得到关于经典解 (u^ \varepsilon, {\rm d}^ \varepsilon, \dot{\rm d}^ \varepsilon) 的存在时间的关于 \varepsilon 的一致下界.其中一个原因是我们不能得到参数化双曲方程(4.1)–(4.3)的解的一致 H^s -范数估计.主要困难在于能量中参数 \varepsilon 的次数和Leslie张量中 \varepsilon 的次数不匹配造成的.如果我们对参数化的方程(4.1)–(4.3)做 H^s -能量估计, 可以观察到 \|\nabla {\rm d}^ \varepsilon\|_{H^s} \|\dot{\rm d}^ \varepsilon\|_{H^s} 关于 \varepsilon \frac{1}{ \varepsilon} 次的.但是方程的右边会出现一项 \varepsilon^3\|{\rm d}^ \varepsilon\|_{H^s}\|\dot{\rm d}^ \varepsilon\|_{H^s}^3 , 随着 \varepsilon 趋向于 0 , 这一项会变得非常大.

参考文献

De Anna F , Zarnescu A .

Global well-posedness and twist-wave solutions for the inertial Qian-Sheng model of liquid crystals

J Differential Equations, 2018, 264: 1080- 1118

DOI:10.1016/j.jde.2017.09.031      [本文引用: 1]

Ericksen J L .

Conservation laws for liquid crystals

Trans Soc Rheology, 1961, 5: 23- 34

DOI:10.1122/1.548883      [本文引用: 1]

Ericksen J L .

Continuum theory of nematic liquid crystals

Res Mechanica, 1987, 21: 381- 392

URL    

Ericksen J L .

Liquid crystals with variable degree of orientation

Arch Rational Mech Anal, 1990, 113: 97- 120

URL     [本文引用: 1]

Feireisl E , Rocca E , Schimperna G , Zarnescu A .

On a hyperbolic system arising in liquid crystals modeling

J Hyperbolic Differ Equa, 2018, 15: 15- 35

DOI:10.1142/S0219891618500029      [本文引用: 1]

Feireisl E , Rocca E , Schimperna G , Zarnescu A .

Nonisothermal nematic liquid crystal flows with the Ball-Majumdar free energy

Annali di Mat Pura ed App, 2015, 194 (5): 1269- 1299

DOI:10.1007/s10231-014-0419-1      [本文引用: 1]

Ladyzhenskaya O A, Solonnikov N A, Uraltseva N N. Linear and Quasilinear Equations of Parabolic Type. Providence, RI: American Mathematical Society, 1968

[本文引用: 1]

Leslie F . Theory of Flow Phenomena in Liquid Crystals. New York: Academic Press, 1979: 1- 81

[本文引用: 1]

Leslie F .

Some constitutive equations for liquid crystals

Arch Rational Mech Anal, 1968, 28: 265- 283

DOI:10.1007/BF00251810      [本文引用: 1]

Majda A J , Bertozzi A L . Vorticity and Incompressible Flow. Cambrige: Cambridge University Press, 2002

[本文引用: 1]

Gennes P G. The Physics of Liquid Crystals. London: Oxford, 1974

Guillén-González F , Rodríguez-Bellido M A .

Weak solutions for an initial boundary Q-tensor problem related to liquid crystals

Nonlinear Anal, 2015, 112: 84- 104

DOI:10.1016/j.na.2014.09.011      [本文引用: 1]

Hardt R , Kinderlehrer D , Lin F H .

Existence and partial regularity of static liquid crystal configurations

Comm Math Phys, 1986, 105: 547- 570

DOI:10.1007/BF01238933      [本文引用: 1]

Hong M C , Xin Z .

Global existence of solutions of the liquid crystal fow for the Oseen-Frank model in {\mathbb R}^2

Adv Math, 2012, 231: 1364- 1400

DOI:10.1016/j.aim.2012.06.009      [本文引用: 1]

Huang J , Lin F H , Wang C .

Regularity and existence of global solutions to the Ericksen-Leslie system in {\mathbb R}^2

Comm Math Phys, 2014, 331: 805- 850

DOI:10.1007/s00220-014-2079-9      [本文引用: 1]

Kato T .

Nonstationary flows of viscous and ideal fluids in {\mathbb R}^3

J Funct Anal, 1972, 9: 296- 305

DOI:10.1016/0022-1236(72)90003-1      [本文引用: 1]

Kato T .

On classical solutions of the two-dimensional nonstationary Euler equation

Arch Rational Mech Anal, 1967, 25: 188- 200

DOI:10.1007/BF00251588      [本文引用: 1]

Ionescu A D, Kemig C E. Local and Global Well-posedness of Periodic KP-I Equations//Bourgain J, Kenig C E, Klainerman S. Mathematical Aspects of Nonlinear Dispersive Equations. Princeton: Princeton University Press, 2009: 181-212

[本文引用: 1]

Jiang N , Luo Y L .

On well-posedness of Ericksen-Leslie's hyperbolic incompressible liquid crystal model

SIAM J Math Anal, 2019, 51 (1): 403- 434

DOI:10.1137/18M1167310      [本文引用: 5]

Jiang N , Luo Y L , Tang S J .

On well-posedness of Ericksen-Leslie's parabolic-hyperbolic liquid crystal model in compressible flow

Math Models Methods Appl Sci, 2019, 29: 121- 183

DOI:10.1142/S0218202519500052      [本文引用: 1]

Li J , Titi E S , Xin Z .

On the uniqueness of weak solutions to the Ericksen-Leslie liquid crystal model in {\mathbb R}^2

Math Models Methods Appl Sci, 2016, 26: 803- 822

DOI:10.1142/S0218202516500184      [本文引用: 1]

Lin F H , Liu C .

Nonparabolic dissipative systems modeling the flow of liquid crystals

Comm Pure Appl Math, 1995, 48: 501- 537

DOI:10.1002/cpa.3160480503      [本文引用: 1]

Lin F H , Liu C .

Partial regularity of the dynamic system modeling the flow of liquid crystals

Discrete Contin Dynam Systems, 1996, 2: 1- 22

DOI:10.3934/dcds.1996.2.1     

Lin F H , Liu C .

Existence of solutions for the Ericksen-Leslie system

Arch Ration Mech Anal, 2000, 154: 135- 156

DOI:10.1007/s002050000102     

李强.

分数阶扩散的三维液晶方程的整体正则性

数学物理学报, 2020, 40A (4): 918- 924

DOI:10.3969/j.issn.1003-3998.2020.04.009      [本文引用: 1]

Li Q .

Global regularity for the 3D liquid crystals equations with fractional diffusion

Acta Math Sci, 2020, 40A (4): 918- 924

DOI:10.3969/j.issn.1003-3998.2020.04.009      [本文引用: 1]

Wang W , Zhang P , Zhang Z .

Well-posedness of the Ericksen-Leslie system

Arch Ration Mech Anal, 2013, 210: 837- 855

DOI:10.1007/s00205-013-0659-z      [本文引用: 1]

Wang M , Wang W .

Global existence of weak solution for the 2-D Ericksen-Leslie system

Calc Var Partial Differential Equations, 2014, 51: 915- 962

DOI:10.1007/s00526-013-0700-y      [本文引用: 1]

Wu H , Xu X , Liu C .

On the general Ericksen-Leslie system: Parodi's relation, well-posedness and stability

Arch Ration Mech Anal, 2013, 208 (1): 59- 107

DOI:10.1007/s00205-012-0588-2      [本文引用: 1]

Wu H , Xu X , Zarnescu A .

Dynamics and flow effects in the Beris-Edwards system modeling nematic liquid crystals

Arch Ration Mech Anal, 2019, 231 (2): 1217- 1267

DOI:10.1007/s00205-018-1297-2      [本文引用: 2]

/