粘性Cahn-Hilliard方程的可解性
Solvability of Viscous Cahn-Hilliard Equation
DOI: 10.12677/AAM.2020.910200, PDF, HTML, XML, 下载: 744  浏览: 1,918 
作者: 黄 梅, 蒲志林, 段 芳:四川师范大学,四川 成都
关键词: 粘性Cahn-Hilliard方程局部解全部解Viscous Cahn-Hilliard Equation Local Solution Global Solution
摘要: 本文研究粘性Cahn-Hilliard方程的可解性问题。通过构造扇形算子并且利用文献中的方法证明了粘性Cahn-Hilliard方程局部解的存在性和整体解的存在性。同时,考虑带有非线性的粘性Cahn-Hilliard方程,研究其可解性问题。
Abstract: In this paper, we study the solvability of the viscous Cahn-Hilliard equation. By constructing sectorial operators and empolying the method in reference, the existence of local solutions and global solutions for viscous Cahn-Hilliard equation are proved. At the same time, the existence of solutions for the viscous Cahn-Hilliard equation with nonlinear term is studied.
文章引用:黄梅, 蒲志林, 段芳. 粘性Cahn-Hilliard方程的可解性[J]. 应用数学进展, 2020, 9(10): 1734-1742. https://doi.org/10.12677/AAM.2020.910200

1. 引言

考虑如下的粘性Cahn-Hilliard方程及初边值问题:

( 1 ν ) u t = Δ μ ( u ) , u = u ( t , x ) , ( t , x ) R × Ω , (1)

μ ( u ) = Δ u f ( u ) + ν u t , (2)

u ( t , x ) = Δ u ( t , x ) = 0 , x Ω , (3)

u N = ( Δ u ) N = 0 , x Ω , (4)

其中 ν [ 0 , 1 ] Ω R N ( N = 1 , 2 , 3 ) 上的一个有界开集, Δ 为Laplace算子, u = u ( t , x ) 是未知函数,f是最高次为 2 p 1 的一个多项式,且有如下形式:

f ( u ) = j = 1 2 p 1 a j u j , a 2 p 1 > 0 , p N , p 2 (5)

粘性Cahn-Hilliard方程是描述二元合金分解的典型Cahn-Hilliard模型的推广,由A.Novick-Cohn在文献 [1] 中提出用来分析粘性一阶变相动力学。 ν = 0 时,是适当限制条件下的Cahn-Hilliard方程; ν = 1 时,是我们熟悉的线性热方程; ν ( 0 , 1 ) 时,是本文所研究的模型。

粘性Cahn-Hilliard 方程解的存在性已有文献研究,例如文献 [2] 研究了广义粘性Cahn-Hilliard方程在RN中的可解性,这一结果在 [3] 中得到了推广;文献 [4] 考虑粘性Cahn-Hilliard方程的解析解,证明了解的存在唯一性。在Neumann边界条件下的粘性Cahn-Hilliard方程,其动力学性质研究可参考文献 [5] [6];无界区域上的粘性Cahn-Hilliard方程的动力学性质研究可参考文献 [7];关于 ν ( 0 , 1 ) 时方程的渐进行为也有详细研究,例如 [8]。文献 [9] 采用质量守恒的crank-Nicolson型有限差分法,对粘性Cahn-Hilliard方程进行了数值求解。

本文研究粘性Cahn-Hilliard方程可解性问题。首先利用扇形算子及文献 [10] 中解的存在性定理,证明了初边值问题(1.1)~(1.4)局部解的存在性;其次在局部解存在的前提下,借鉴文献 [10] 中全局解的存在性定理,证明了初边值问题(1.1)~(1.4)的全局解在一定条件下的存在性,即任意给定 u 0 L 2 ( Ω ) ,初边值问题(1.1)~(1.4)的解u,满足 Δ u 2 是有界的,从而得到全局解存在。最后我们考虑带有非线性的粘性Cahn-Hilliard方程,即

( 1 ν ) u t + g ( u ) = Δ ( Δ u f ( u ) + ν u t ) , u = u ( t , x ) , ( t , x ) R × Ω ,

其中g是多项式函数。文献 [11] 研究了它的渐进形态,本文研究它的解的存在情况。

2. 预备知识

考虑Banach空间X中的抽象Cauchy问题:

{ d u d t + A u = F ( u ) , t > 0 u ( 0 ) = u 0 (6)

其中, A : D ( A ) X X 是扇形正算子,分数幂算子 A α : X X 定义为:

A α v = 1 Γ ( α ) 0 + t α 1 e A t v d t ,

其中 α ( 0 , + ) 。由 [10] 知道 A α ( α ( 0 , + ) ) 是有界线性算子,并且 A α 是可逆的。记 A α = ( A α ) 1 X α : = D ( A ) ;假设 F : X α X X α 的有界集上Lipschtz连续。

定义2.1 [10] (局部 X α 解):设X是一个Banach空间, α [ 0 , 1 ) u 0 X α 。若存在 τ > 0 和函数 u ( C ( 0 , τ ) , X α ) 满足

1) u ( 0 ) = u 0

2) u C 1 ( ( 0 , τ ) , X )

3) t ( 0 , τ ) , u ( t ) D ( A )

4) t ( 0 , τ ) ,u满足方程(6);

则u称为(2.1)的局部 X α 解。

定义2.2 [10] (全局 X α 解):函数 u = u ( t ) 称为全局 X α 解,如果它满足定义2.1且 τ = +

定理2.1 [10] (局部 X α 解的存在性):对每个 u 0 X α ,存在(6)的唯一 X α u = u ( t , u 0 ) ,它是定义在最大存在区间 [ 0 , τ u 0 ) 上的。即 τ u 0 = + 或者如果 τ u 0 < + ,则 lim t τ u 0 sup u ( t , u 0 ) X α = +

关于方程全局 X α 解的存在性是非常重要的。文献 [10] 引入了如下的一个假设条件H:

1) 存在一个Banach空间Y, D ( A ) Y

2) 存在一个局部有界函数 C : [ 0 , + ) [ 0 , + )

3) 存在一个非减函数 g : [ 0 , + ) [ 0 , + )

4) 存在某个数 θ [ 0 , 1 )

使得 u 0 X α ,满足:

u ( t , u 0 ) Y C ( u 0 X α ) , t ( 0 , τ u 0 ) ,

F ( u ( t , u 0 ) ) X g ( u ( t , u 0 ) Y ) ( 1 + u ( t , u 0 ) X α θ ) , t ( 0 , τ u 0 ) .

定理2.2 [10] (全局 X α 解的存在性):如果假设H成立,则(6)存在全局 X α 解。

引理2.1 [12]:对 η > 0

{ | Δ u | 2 + η | u | 2 } 1 2 , { | Δ u | 2 + η ( Ω u ( x ) d x ) 2 } 1 2 ,

H 2 上的等价范数,同样地,

{ | Δ 2 u | 2 + η | u | 2 } 1 2 , { | Δ 2 u | 2 + η ( Ω u ( x ) d x ) 2 } 1 2 ,

H 4 上的等价范数。

3. 主要结果

本节给出粘性Cahn-Hilliard方程局部解和全局解的存在性定理。令 X = L 2 ( Ω ) ,内积记为 ( , ) ,范数记为 。更一般地,把Banach空间Z上的范数记为 Z

显然,由关于f的表达式(5)知,存在 λ 0 > 0 ,使得对任意 v R

| f ( v ) | < λ 0 . (7)

u L 2 ( Ω ) ,均值定义为:

u = 1 | Ω | Ω u ( x ) d x , x Ω ,

其中 | Ω | Ω 的测度。

对(1)作 Ω 上的积分,利用格林公式可得:

d d t Ω u ( t , x ) d x = 0 ,

因此,

u = c o n s t . = u 0 , (8)

故粘性Cahn-Hilliard方程是满足质量守恒的,这对后面范数估计起着至关重要的作用。

令算子 A = ( Δ ) B ν + ρ I ( ρ > 0 ) B ν = Δ [ ( 1 ν ) I ν Δ ] 1 ,其定义域为:

D ( A ) = { ϕ H 4 : u N = ( Δ u ) N | Ω = 0 } ,

首先讨论算子A的性质,有如下的命题:

命题 3.1当 ν ( 0 , 1 ) 时,算子 B ν 是X到X上的有界算子,更进一步, Δ B ν : H 2 ( Ω ) X 是一个扇形算子。

证明:算子 B ν 还可写为

B ν = 1 ν I 1 ν ν [ ( 1 ν ) I ν Δ ] 1 ,

由扇形算子定义可知: M 0 > 0 ,使得

( ( 1 ν ) I ν Δ ) 1 M 0 , (9)

故可证得 B ν : L 2 ( Ω ) L 2 ( Ω ) 是一个有界算子。

注意到,算子 Δ B ν 有以下形式:

Δ B ν = 1 ν ( Δ ) 1 ν ν 2 I + ( 1 ν ) 2 ν 2 [ ( 1 ν ) I ν Δ ] 1 ,

它是扇形算子 1 ν ( Δ ) 的一个有界扰动,因此 Δ B ν 是一个扇形算子(参考 [10] )。

定理 3.1 当 α [ 1 2 , 1 ) ν ( 0 , 1 ) 时,初边值问题(1)~(4)存在局部 X α 解。

证明:由命题3.1可知,要证明局部 X α 解的存在性,还需要验证非线性项

B ν f : X α ( Ω ) X ,

满足局部Lipschtiz连续。由(5)可知, f ( u ) 关于 u Ω 局部Lipschitz连续。

更进一步,有如下的Sobolev空间嵌入关系:

{ H 2 ( Ω ) L ( Ω ) n 3 H 2 ( Ω ) W 1 , 4 ( Ω ) n 3 X 1 2 ( Ω ) H 2 ( Ω ) , (10)

令U是 X α 上的有界子集,则 B ν f : X α X 在有界集U上是Lipschitz连续的,即: ϕ , ψ U

B ν ( f ( ϕ ) ) B ν ( f ( ψ ) ) ( ( 1 ν ) I ν Δ ) 1 L ( X , X ) Δ f ( ϕ ) Δ f ( ψ ) M 0 [ f ( ϕ ) Δ ϕ f ( ψ ) Δ ψ ] + [ f ( ϕ ) | ϕ | 2 f ( ψ ) | ψ | 2 ] C U ϕ ψ X α ,

其中 C U 是关于U的常数。因此初边值问题(1)~(4)的局部 X α 解存在。

在局部 X α 解存在的基础上进一步讨论全局解的存在性。

定理 3.2当 α [ 1 2 , 1 ) ν ( 0 , 1 ) 时,在条件(5)和(7)成立的情况下,初边值问题(1)~(4)存在唯一的全局 X α 解。

证明:用 v = u t + ϵ 0 u 与(1)作内积有

1 2 d d t Δ u 2 + ϵ 0 Δ u 2 + ( 1 ν ) v 2 ( 1 ν ) ϵ 0 ( u , v ) + ν ϵ 0 ( Δ u , v ) + ν v = ( Δ f ( u ) , v ) , (11)

由Poincare不等式得:

u 2 λ u 2 , (12)

由(12)以及Hölder不等式和Young不等式可得

u 2 = Ω u Δ u d x 1 2 u 2 + λ 2 Δ u 2 ,

u 2 λ Δ u 2 , (13)

更进一步,由(12)和(13)得

u 2 λ 2 Δ u 2 ,

又由Hölder不等式和Young不等式可得

( 1 ν ) ϵ 0 ( u , v ) ( 1 ν ) ϵ 0 λ 2 2 Δ u 2 ( 1 ν ) ϵ 0 2 v 2 , (14)

ν ϵ 0 ( Δ u , v ) ϵ 0 2 ν λ 4 Δ u 2 ν v 2 , (15)

将(14)和(15)代入(11)可得

1 2 d d t Δ u 2 + c 1 Δ u 2 + c 2 v 2 ( Δ f ( u ) , v ) , (16)

其中 c 1 = ϵ 0 ( 1 ν ) ϵ 0 λ 2 2 ϵ 0 2 ν λ 4 , c 2 = ( 1 ν ) ( 1 ϵ 0 2 )

最后由Young不等式可得

( Δ f ( u ) , v ) = ( f ( u ) , v ) λ 0 2 4 c 2 Δ u 2 + c 2 v 2 , (17)

将(17)代入(16)得

1 2 d d t Δ u 2 + c 3 Δ u 2 1 ,

其中 c 3 = c 1 λ 0 2 4 c 2 ,取适当的 ϵ 0 ,使得 c 3 > 0 。则由Gronwall引理得:

Δ u 2 M ,

其中M是一个大于0的常数。再由(8)得

Δ u 2 + | u | 2 M + | u 0 | 2 .

Y = { ϕ H 2 ( Ω ) : ϕ N = ( Δ ϕ ) N | Ω = 0 } ,由引理2.1可知Y上的等价范数为 { Δ u 2 + | u | 2 } 1 2 。为了利用定理2.2完善证明,还需近一步处理。

将(1)改写为:

u t = ( B ν Δ + ρ I ) u B ν f ( u ) ρ u ,

其中 ρ > 0 且有

0 < ρ R e σ ( B ν Δ + ρ I ) ,

由空间嵌入定理(10),在Y空间上有下列先验估计:

B ν f ( u ( t ) ) ρ u ( t ) B ν L ( X , X ) f ( u ( t ) ) + ρ u ( t ) B ν L ( X , X ) sup | s | c 4 u ( t ) H 2 ( Ω ) ( | f ( s ) | + ρ ) ( 1 + u ( t ) H 2 ( Ω ) ) = g ( u ( t ) H 2 ( Ω ) ) ,

其中 g : [ 0 , ) [ 0 , ) 是一个增函数,证毕。

以上的讨论过程是粘性Cahn-Hilliard方程的解的存在情况研究的新思路,拓展了粘性Cahn-Hilliard方程研究的方法,用新的方法证明了粘性Cahn-Hilliard方程局部解和全局解的存在性。

4. 推论

这一节我们考虑粘性Cach-Hilliard方程带有非线性的形式,即考虑下面的初边值问题:

( 1 ν ) u t + g ( u ) = Δ μ ( u ) , u = u ( t , x ) , ( t , x ) R × Ω , (18)

μ ( u ) = Δ u f ( u ) + ν u t , (19)

u ( t , x ) = Δ u ( t , x ) = 0 , x Ω , (20)

u N = ( Δ u ) N = 0 , x Ω , (21)

其中f的定义如前面的(5)而g的定义如下:

g ( u ) = i = 1 2 q 1 b i u i , b 2 q 1 > 0 , q N , q 2 , (22)

有不等式:

g ( s ) s q b 2 q 1 s 2 q c , s R , (23)

由关于g的表达式(22)知

λ 1 > 0 , v R , | g ( v ) | < λ 1 , (24)

其中 λ 1 是一个常数。显然该初边值问题不满足质量守恒,下面讨论该初边值问题的局部 X α 解和全局 X α 解。

根据上一节中局部解的讨论过程,(10)同样成立。很容易可以得到当 α [ 1 2 , 1 ) 时该初边值问题的局部 X α ( α [ 1 2 , 1 ) ) 解存在。令V是 X α 空间上的有界集, ϕ , ψ V ,则该方程的非线性项 B ν f ( ( 1 ν ) I ν Δ ) 1 g : X α L 2 ( Ω ) 满足Lipschitz连续,即

[ B ν ( f ( ϕ ) ) ( ( 1 ν ) I ν Δ ) 1 g ( ϕ ) ] [ B ν ( f ( ψ ) ) ( ( 1 ν ) I ν Δ ) 1 g ( ψ ) ] M 0 Δ f ( ϕ ) Δ f ( ψ ) + M 0 g ( ϕ ) g ( ψ ) M 0 [ f ( ϕ ) Δ ϕ f ( ψ ) Δ ψ ] + [ f ( ϕ ) | ϕ | 2 f ( ψ ) | ψ | 2 ] + M 0 ( g ( ϕ ) L ( Ω ) + g ( ψ ) L ( Ω ) ) C V ϕ ψ X α ,

其中 C V 是关于V的常数。故可以知道该初边值问题的局部 X α ( α [ 1 2 , 1 ) ) 解存在。

定理4.1 当 α [ 1 2 , 1 ) ν ( 0 , 1 ) 时,在条件(5),(7)和(22)~(24)成立的情况下,初边值问题(18)~(21)存在唯一的全局 X α 解。

证明:用 w = u t + ϵ 1 u 与(18)作内积有u

1 2 d d t Δ u 2 + ϵ 1 Δ u 2 + ( 1 ν ) w 2 ( 1 ν ) ϵ 1 ( u , w ) + ν ϵ 1 ( Δ u , w ) + ν v + ( g ( u ) , w ) = ( Δ f ( u ) , w )

根据Hölder不等式和Young不等式可得

( 1 ν ) ϵ 1 ( u , w ) ( 1 ν ) ϵ 1 λ 2 2 Δ u 2 ( 1 ν ) ϵ 1 2 w 2 , (25)

ν ϵ 1 ( Δ u , w ) ϵ 1 2 ν λ 4 Δ u 2 ν w 2 , (26)

由(25)和(26)可得

1 2 d d t Δ u 2 + c 5 Δ u 2 + c 6 w 2 + ( g ( u ) , w ) ( Δ f ( u ) , w ) ,

其中 c 5 = ϵ 1 ( 1 ν ) ϵ 1 λ 2 2 ϵ 1 2 ν λ 4 , c 6 = ( 1 ν ) ( 1 ϵ 1 2 )

又由

( Δ f ( u ) , w ) = ( f ( u ) , w ) λ 0 2 4 c 6 Δ u 2 + c 6 w 2 ,

1 2 d d t Δ u 2 + ( c 5 λ 0 2 4 c 6 ) Δ u 2 + ϵ 1 Ω g ( u ) u d x Ω g ( u ) u t d x ,

由(23)以及Young不等式和格林公式有

1 2 d d t Δ u 2 + ( c 5 λ 0 2 4 c 6 λ 1 2 λ 2 2 ) Δ u 2 + ϵ 1 q b 2 q 1 u L 2 q ( Ω ) 2 q c | Ω | + 1 2 u t 2 , (27)

再将 u t 与(18)作内积得

u t 2 + u t 2 + 1 2 d d t Δ u 2 + Ω g ( u ) u t d x = ( Δ f ( u ) , u t ) ,

由(7),(24)以及格林公式和Yong不等式得到

1 2 u t 2 1 2 d d t Δ u 2 + ( λ 0 2 λ 2 + λ 1 2 λ 2 2 ) Δ u 2 , (28)

将(28)带入(27),再由空间嵌入定理有

d d t Δ u 2 + c 7 Δ u 2 + ϵ 1 q b 2 q 1 u L 2 q ( Ω ) 2 q c | Ω | ,

其中 c 7 = c 5 λ 0 2 4 c 6 λ 0 2 λ 2 λ 1 2 λ 2 。取适当 ϵ 1 ,使得 c 7 > 0 ,则由Gronwall引理得: M 1 > 0

Δ u 2 M 1 ,

Δ u 2 + u 2 ( 1 + λ 2 ) M 1 .

将(18)改写为:

u t = ( B ν Δ + ρ I ) u B ν f ( u ) ( ( 1 ν ) I ν Δ ) 1 g ( u ) ρ u ,

由空间嵌入定理(10),在Y空间上有下列先验估计:

B ν f ( u ( t ) ) ( ( 1 ν ) I ν Δ ) 1 g ( u ( t ) ) ρ u ( t ) B ν L ( X , X ) f ( u ( t ) ) + M 0 g ( u ( t ) ) + ρ u ( t ) max ( B ν L ( X , X ) , M 0 ) sup | s | c 8 u ( t ) H 2 ( Ω ) ( | f ( s ) | + | g ( s ) | + ρ ) ( 1 + u ( t ) H 2 ( Ω ) ) = g ( u ( t ) H 2 ( Ω ) ) ,

其中 g : [ 0 , ) [ 0 , ) 是一个增函数,证毕。

以上的结果,拓展了粘性Cach-Hilliard方程的讨论范围,得出了在一定条件下带有非线性的粘性Cach-Hilliard方程局部解和全局解的存在性。当然,我们所作的假设条件还可以得到更一般的结论,这将是我们后面所要努力的。

参考文献

[1] Novick-Cohen, A. (1988) On the Viscous Cahn-Hilliard Equation, Material Instabilities in Continuum Mechanics and Related Mathematical Problems. Clarendon Press, Oxford, 329-342.
[2] Dlotko, T., Kania, M.B. and Sun, C.Y. (2012) Analysis of the Viscous Cahn-Hilliard Equation in RN. Journal of Differential Equation, 252, 2771-2791.
https://doi.org/10.1016/j.jde.2011.08.052
[3] Bui, L.T.T., Dao, A.N. and Díaz, J.L. (2017) Critical Case for the Viscous Cahn-Hilliard Equation. Electronic Journal of Differential Equations, 176, 1-8.
[4] Choo, S.M. and Chung, S.K. (2003) Asymtotic Behaviour of the Viscous Cahn-Hilliard Equation. Journal of Applied Mathematics and Computing, 11, 143-154.
https://doi.org/10.1007/BF02935727
[5] Rossi, R. (2007) Global Attractor for the Weak Solutions of a Class of Viscous Cahn-Hilliard Equations, Dissipative Phase Transitions. Series on Advances in Mathematics for Applied Sciences, 71, 247-268.
https://doi.org/10.1142/9789812774293_0013
[6] Carvalho, A.N. and Dlotko, T. (2008) Dynamics of the Viscous Cahn-Hilliard Equation. Journal of Mathematical Analysis and Applications, 344, 703-725.
https://doi.org/10.1016/j.jmaa.2008.03.020
[7] Cholewa, J.W. and Dlotko, T. (2010) Dynamics of the Modified Viscous Cahn-HillIiard Equation in RN. Topological Methods in Nonlinear Analysis, 35, 277-294.
[8] Elliott, C.M. and Stuart, A.M. (1996) Viscous Cahn-Hilliard Equation II. Analysis. Journal of Differential Equations, 128, 387-414.
https://doi.org/10.1006/jdeq.1996.0101
[9] Choo, S.M. and Chung, S.K. (2004) A Conservative Nonlinear Difference Scheme for the Viscous Cahn-Hilliard Equation. Journal of Applied Mathematics and Computing, 16, 53-68.
https://doi.org/10.1007/BF02936150
[10] Cholewa, J.W. and Dlotko, T. (2000) Global Attractors in Abstract Parabolic Problems. Cambridge University Press, Cambridge.
https://doi.org/10.1017/CBO9780511526404
[11] 任永华. 粘性Cahn-Hilliard方程的整体吸引子[J]. 应用数学进展, 2016, 5(4): 818-822.
[12] Pazy, A. (1983) Semigroups of Linear Operator and Applications to Partial Differential Equation. Spring Verlag, New York.
https://doi.org/10.1007/978-1-4612-5561-1