1. 引言
Cahn-Hilliard方程最初是由Cahn和Hilliard在文献 [1] 中介绍,用以描述熔融合金中复杂的相分离和粗化现象,该合金被淬灭到只有两个不同的浓度相可以稳定存在的温度。随着理论的发展,对Cahn-Hilliard方程的研究也受到众多学者追捧,Cahn-Hilliard方程在流体力学中具有重要应用,并被许多学者广泛研究 [2] [3] [4] [5]。Cahn-Hilliard相场方程经常用来描述油水相分离的行为,在工业上有着重要的应用。在不同的研究领域Cahn-Hilliard方程的模型也各有不同,本文要研究的Cahn-Hilliard方程模型如下 [6] [7]:
(1)
在方程中u表示混合物的浓度,
是一个与界面厚有关的参数(
),n为边界
上的外法向量,其中f是光滑函数F的导数。特别地,取
其中f需要满足:存在常数L,使得
Cahn-Hilliard方程它有一个重要的性质就是可以做Liapunov能量函数 [1]
称作自由能泛函,它能准确的表述混合物的自由能。
2. Cahn-Hilliard方程的能量耗散性
2.1. 能量耗散
本文要研究的稳定性Cahn-Hilliard方程能量稳定的 [8],给出下面的引理
引理1 Cahn-Hilliard方程定义的自由能
满足下面的定律
即Cahn-Hilliard方程是满足能量耗散的。
证明 取
与(1)式的第一个式子做内积,有
式中
,
表示内积。取
由格林公式以及边界条件可得
即
证毕。
2.2. 能量不变二化法
为了构造线性数值格式,我们令(1)式中
的,Cahn-Hilliard方程变形为 [9]
(2)
接下来我们采用能量不变二化法来分析变换后是否满足能量耗散规律。
引理2 Cahn-Hilliard方程采用能量不变二化法后定义的自由能
满足下面的定律
即Cahn-Hilliard方程采用能量不变二化法也是满足能量耗散的。
证明 为了采用能量不变二化法,我们作以下标记:
则(2)式变形为
(3)
用(3)式中三个式子分别与
作内积可得
所以我们采用能量不变二分法满足能量耗散,其中
,这里的
和
是等价的。
3. 稳定性分析
利用格林公式和边界条件,很容易得到Cahn-Hilliard方程的变分形式 [10] [11] [12] 为:求
,使得下面的式子成立
针对上面的变分形式,构造时间方向的半离散格式
(4)
其中τ为时间步长,且该格式满足下面的能量估计式成立:
定理 (能量稳定性)方程(4)的解满足下面的能量稳定,即
其中,
证明 在(4)式中取
,
有
(5)
再把方程(3)中的第三式与
作内积可得
(6)
综合(5)、(6)两式有
整理可得该定理的结果。
证毕。
4. 数值实验
为了能够有效观测Cahn-Hilliard方程相场模型的相位变化情况,采用混合有限元方法给出一个数值算例,选取精确解为
计算区域为方形区域
,其中网格尺寸选择1/64,
。
图1描绘了在
时,T = 1。0e−5,T = 1。0e−4,T = 1。0e−3,T = 1。0e−2时的相位变化图。
图2描绘了在
时,T = 1。0e−5,T = 1。0e−4,T = 1。0e−3,T = 1。0e−2时的相位变化图。
Figure 1. When ε = 0.01, the image of the change in the numerical solution u at different moments
图1. 当ε = 0.01时, 数值解u在不同时刻的变化图像
Figure 2. When ε = 0.1, the image of the change in the numerical solution u at different moments
图2. 当ε = 0.01时,数值解u在不同时刻的变化图像
从上面两图可以看到,所构造的数值算例有效的模拟了Cahn-Hilliard方程的相位变化过程,起初两种物质任意分布,随着时间的推移,两种物质逐渐趋于缓和,最后到达一种稳定状态。
比较图1和图2可以发现,当ε取不同值时,对两种物质的反应速率有很大影响,ε越小反应越剧烈,反应速率越快,相当于起到了催化剂的作用。