在本文中,我们创建一个新尝试延长静电学的吉布斯变分框架。我们的方法非常简单,完全简单的概念上。与之前的许多尝试([
12- - - - - -
14),三),我们明确排除电场和电位移独立热力学变量的列表。相反,我们占极化(或磁化),通过添加一个术语的“传统的”热弹性系统的自由能。额外的术语表示静电场中的势能积累在整个空间。不同作者选择这一项的不同:
E
2,
E
·
D,等等。我们选择被积函数在最简单的形式
E
2。我们构建我们的方法在精确的非线性连续介质理论和依赖欧拉坐标作为独立空间变量。
我们可以重写(
17)(见[
22,
23,
27])
(18)
∇
我
ℵ
我
j
=
0
,aleph张量的
ℵ
我
j,由
(19)
ℵ
我
j
≡
ρ
∂
ψ
∂
∇
我
U
k
一个
k
·
·
j
+
1
8
π
E
k
E
k
- - - - - -
1
4
π
E
k
D
k
z
我
j
+
1
4
π
D
我
E
j
,可以被认为是应力张量的可极化的物质。我们可以重写张量
ℵ
我
j作为
(20)
ℵ
我
j
≡
ζ
我
j
+
Γ
我
j
,静电伽马张量在哪里
Γ
米
k是由
(21)
Γ
我
j
≡
1
8
π
E
k
E
k
- - - - - -
1
4
π
E
k
D
k
z
我
j
+
1
4
π
D
我
E
j
。方程(
17)可以写在另一个深刻的形式:
(22)
∇
我
ζ
我
j
=
- - - - - -
∇
我
Γ
我
j
。在可极化的可变形的物质,没有一个张量,
ζ
我
j或
Γ
我
jdivergence-free。
γ张量
Γ
我
j也可以认为是一个许多可能的麦克斯韦应力张量的概括
T
我
j:
(23)
T
我
j
≡
- - - - - -
1
8
π
E
k
E
k
z
我
j
+
1
4
π
E
我
E
j
,自
Γ
我
j正值
T
我
j当极化消失了。其他可能的麦克斯韦应力张量的概括,
(24)
T
1
我
j
≡
- - - - - -
1
8
π
E
l
D
l
z
我
j
+
1
4
π
D
我
E
j
,
(24 b)
T
2
我
j
≡
- - - - - -
1
8
π
E
l
D
l
z
我
j
+
1
4
π
D
j
E
我
,
(24)
T
3
我
j
≡
- - - - - -
1
8
π
E
l
D
l
z
我
j
+
1
8
π
D
我
E
j
+
D
j
E
我
,可能是审美吸引力比γ张量
Γ
我
j。我们相信伽马张量的优势比其他可能的概括是变分的起源和其能力来帮助解决这一问题的稳定性计算的基础上第二能量变化。
一个更有用的材料可极化张量是贝丝张量
ℶ
·
j
我
·,或者是
张量的电化学张量的潜力。它被定义为
(25)
ℶ
·
j
我
·
≡
ρ
ψ
z
我
k
- - - - - -
ℵ
我
k
+
Γ
我
k
B
k
j
,在张量
B
k
j矩阵逆的吗
一个
k
j中定义的(
15)。我们下面显示,贝丝张量
ℶ
·
j
我
·满足条件的散度为零
(26)
∇
我
ℶ
·
j
我
·
=
0
,类似于张量
ℵ
我
j。贝丝张量
ℶ
·
j
我
·可以写成
(27)
ℶ
·
j
我
·
=
ρ
B
k
j
χ
我
k
,在哪里
χ
我
k是博文
对称张量的化学势
(28)
χ
我
j
=
ψ
z
我
j
- - - - - -
1
ρ
ζ
我
j
=
ψ
z
我
j
+
1
ρ
ℶ
我
j
- - - - - -
ℵ
我
j
。
的
对称的张量
χ
我
j应区别于一般
不对称张量的化学张量
μ
我
j:
(29)
μ
我
j
=
z
∘
k
我
z
l
k
χ
j
l
,在哪里
z
∘
米
我逆变的度规张量吗
最初的配置。
一个可以分析模型的极化向量
P
我是固定的(
20.]。然后
ψ
±空间常数,但可能仍然依赖于温度。
8。在静电学Divergence-Free张量
我们最后的三个方程的证明(
18),(
26)和(
37消失的散度。剩下的两个身份同样可以证明。首先,让我们重写gimel张量
ℷ
我
j如下:
(38)
ℷ
我
j
=
ϵ
P
z
我
j
+
1
8
π
E
k
E
k
- - - - - -
1
4
π
E
k
D
k
z
我
j
+
1
4
π
D
我
E
j
。第一学期(
38),我们有
(39)
∇
我
ϵ
P
z
我
j
=
z
我
j
∂
ϵ
P
∂
P
k
∇
我
P
k
。用热力学的身份
(40)
∂
ϵ
P
∂
P
我
≡
E
我
,我们可以重写(
39),
(41)
∇
我
ϵ
P
z
我
j
=
E
我
∇
j
P
我
。第二项的(
38),我们有
(42)
∇
我
z
我
j
1
8
π
E
k
E
k
- - - - - -
1
4
π
E
k
D
k
=
- - - - - -
1
4
π
E
k
∇
j
E
k
- - - - - -
P
k
∇
j
E
k
- - - - - -
E
k
∇
j
P
k
,这可以从以下的身份链:
(43)
第二项
=
∇
我
z
我
j
1
8
π
E
k
E
k
- - - - - -
1
4
π
E
k
D
k
b (43)
=
∇
k
1
8
π
E
k
E
k
- - - - - -
1
4
π
E
k
D
k
(43)
=
1
4
π
E
k
∇
j
E
k
- - - - - -
∇
j
E
k
D
k
- - - - - -
E
k
∇
j
D
k
d (43)
=
1
4
π
E
k
∇
j
E
k
- - - - - -
∇
j
E
k
E
k
+
4
π
P
k
- - - - - -
E
k
∇
j
E
k
+
4
π
P
k
(43)
=
- - - - - -
1
4
π
E
k
∇
j
E
k
- - - - - -
P
k
∇
j
E
k
- - - - - -
E
k
∇
j
P
k
。第三项(
38),我们有
(44)
∇
我
1
4
π
D
我
E
j
=
1
4
π
D
我
∇
我
E
j
=
1
4
π
E
我
∇
我
E
j
+
P
我
∇
我
E
j
。结合(
41)- (
44),我们发现
(45)
∇
我
ℷ
我
j
=
∇
我
ϵ
P
z
我
j
+
z
我
j
∇
我
1
8
π
E
l
E
l
- - - - - -
1
4
π
E
l
D
l
+
1
4
π
∇
我
D
我
E
j
(45 b)
=
E
我
∇
j
P
我
- - - - - -
1
4
π
E
我
∇
j
E
我
- - - - - -
P
我
∇
j
E
我
- - - - - -
E
我
∇
j
P
我
+
1
4
π
E
我
∇
我
E
j
+
P
我
∇
我
E
j
(45度)
=
- - - - - -
1
4
π
E
我
∇
j
E
我
- - - - - -
P
我
∇
j
E
我
+
1
4
π
E
我
∇
我
E
j
+
P
我
∇
我
E
j
。最后,使用对称性质
∇
我
E
j
≡
∇
j
E
我我们到达身份(
37)。
9。准静态演化
可以假定准静态演化类比(
4)。在nondeformable阶段的情况下,它读取
(46)
J
=
- - - - - -
K
N
我
N
j
ℷ
我
j
- - - - - -
+
。同样的方法可以应用到一个孤立的情况下域与固定总量未重排。在这种情况下,应该稍微修改演化方程考虑表面扩散。图
2说明了这种方法在二维情况下的一个实现。最初的准静态演化循环域和固定极化矢量导致伸长极化的方向向量
P
我,并最终形态不稳定。