曲面上的分支组合Calabi流
Branched Combinatorial Calabi Flows on Surfaces
DOI: 10.12677/AAM.2022.1110791, PDF, HTML, XML, 下载: 224  浏览: 315  国家自然科学基金支持
作者: 高开城, 林爱津*:国防科技大学理学院,湖南 长沙
关键词: 分支结构组合Calabi流组合Ricci势Branch Structure Combinatorial Calabi Flows Combinatorial Ricci Potential
摘要: 受蓝师义和戴道清关于分支组合Ricci流的研究启发,我们将组合Calabi流与分支结构结合,引入三角剖分曲面上的分支组合Calabi流。进一步地,我们利用分支组合Ricci势,组合Calabi能量等工具,在二维欧氏空间E2和二维双曲空间H2中,证明分支组合Calabi流的解长时间存在并且指数收敛到一个分支圆包装度量。我们的工作将葛化彬等人关于组合Calabi流的结果推广到了分支组合Calabi流情形。
Abstract: Inspired by the study of the branched combinatorial Ricci flows by Shiyi Lan and Daoqing Dai, we introduce the branched combinatorial Calabi flows on a fixed triangulated surface by combining the branch structure with the combinatorial Calabi flows. Furthermore, in 2 dimensional Euclidean space E2 and 2 dimensional hyperbolic space H2 , using the branched combinatorial Ricci poten-tial and combinatorial Calabi energy etc., we prove that the solutions to the branched combinatorial Calabi flows exist for the long time and converge exponentially to branched circle packing metrics. Our work extends results on the combinatorial Calabi flows by Huabin Ge et al. to the branched combinatorial Calabi flows.
文章引用:高开城, 林爱津. 曲面上的分支组合Calabi流[J]. 应用数学进展, 2022, 11(10): 7451-7463. https://doi.org/10.12677/AAM.2022.1110791

1. 引言

圆包装(circle packing)理论是近几十年来在复分析和离散几何交叉学科中一个快速发展的研究领域。Thurston [1] 在1985年首次提出利用圆包装来逼近Riemann映射的猜想。Rodin和Sullivan [2] 证明了六边形情形下的Thurston猜想。此后,Beardon和Stephenson [3] 提出圆包装可用于构造离散解析函数以及将复分析的经典结果离散化的想法。分支圆包装(branched circle packing)的概念由Bowers和Stephenson [4] 提出,他们将Andreev-Thurston定理推广到了分支圆包装情形。Dubejko [5] 对分支圆包装进行了更加系统的研究,得到了曲面上分支圆包装存在的充要条件。

受Hamilton引入的Ricci流启发,Chow和罗锋 [6] 首次引进了组合Ricci流的概念,并且用其给出了Andreev-Thurston定理的新证明。此后,各种组合曲率流(包括组合Calabi流,组合Yamabe流等)相继被引入 [7] - [13],特别地,葛化彬在博士论文 [8] 中引入了组合Calabi流,并且证明了在二维欧氏空间 E 2 (2 dimensional Euclidean space)和二维双曲空间 2 (2 dimensional hyperbolic space)中组合Calabi流解的长时间存在性及指数收敛性 [8] [9] [10] [11]。圆包装、组合曲率流可以与计算机编程结合起来,具有广泛的应用价值,例如德国的Bobenko团队 [14] 利用圆包装构造3维空间中的离散极小曲面,顾险峰团队 [15] 成功地把组合曲率流应用到计算机图形学、医学成像等领域。

基于Chow和罗锋的组合Ricci流,蓝师义和戴道清 [16] 引入了分支结构(branch structure),提出了分支组合Ricci流,并且将Chow和罗锋的经典结果推广到了分支组合Ricci流情形。受蓝师义和戴道清 [16] 工作的启发,本文将引入分支组合Calabi流,并且尝试将葛化彬 [8] [9],葛化彬和华波波 [10] 以及葛化彬和徐旭 [11] 等关于组合Calabi流的经典研究结果推广到分支组合Calabi流情形。

本文结构如下:第二章是预备知识。第三章讨论离散拉普拉斯算子的性质。第四章给出分支组合Calabi流解的长时间存在性证明。第五章证明分支组合Calabi流解的收敛性。

2. 预备知识

2.1. 分支圆包装

设X是一个带边或无边的紧致曲面,T是曲面X的一个三角剖分。设 V = V ( T ) b d ( T ) i n t ( T ) E = { e i j } F = { Δ i j k } 分别表示T的顶点集,边界点集,内点集,边集以及面集。三角剖分的权重是一个作用于边集E上的函数 Φ : E [ 0 , π 2 ] 。我们称 ( T , Φ ) 是曲面X的一个加权三角剖分,记加权三角剖分曲面为 ( X , T , Φ )

设P是常曲率曲面上圆的一个集合,若P满足下述条件,则称之为实现加权三角剖分 ( T , Φ ) 的圆包装。

V中的点v与P中的圆 C P ( v ) 存在一一映射并且圆 C P ( u ) C P ( v ) 的交角为 Φ = ( e u v )

P保持定向,即若 v 1 v 2 v 3 为T中的面 Δ i j k 按照T的定向排列的三个顶点,则 C P ( v 1 ) C P ( v 2 ) C P ( v 3 ) 在曲面上形成一个正向的三圆组。

下面图1图2是两个圆包装示例(详细信息请参见 [4])。图1图2中均有一个固定的点v被固定在北极点(阴影圆的圆心),就是面向读者的那个点。图1是一个外切圆包装,其中所有权重 Φ 均为0。图2图1具有相同的加权三角剖分,但很明显相较于图1图2更加复杂:图2依旧是外切圆包装,但是它具有两个分支点,分别在北极点和南极点(在球的背面与北极点相对)。沿着阴影圆的外切圆可以发现阴影圆的外切圆共环绕了阴影圆两周,这个圆包装诱导了球面的一个二重覆盖。

Figure 1. Circle packing with weight = 0

图1. 权重均为0的圆包装

Figure 2. Branched circle packing with weight = 0

图2. 权重均为0的分支圆包装

为了更好理解分支圆包装,我们引入圆包装的半径函数概念。若P是实现加权三角剖分 ( T , Φ ) 的一个圆包装,则称给P中的每一个圆 C P ( v ) 的圆心v赋值的半径函数 r : V ( 0 , + ) 为P的圆包装度量。显然,圆包装度量 r = ( r 1 , r 2 , , r N ) 可看作一个N维的向量,其中 N = | V | 是顶点的个数。每一个内点 v i n t ( T ) 的内角和 θ v ( r ) 由圆包装度量决定,说明如下。记点v处的内角和 θ v ( r ) = α v ( r , Δ ) ,其中求和包括所有含有点v的面且 α v ( r , Δ ) = α v ( r , Δ v u 1 u 2 ) 指点v, u 1 u 2 构成的三角形 Δ v u 1 u 2 中点v处的内角。

显然,由圆包装度量的定义可得,任取内点 v i n t ( T ) ,它在圆包装P中的内角和一定为 2 π 的整数倍。若在圆包装P中某个内点u处的内角和等于 2 n π n 2 ,那么我们称u为P中分支度为 n 1 的分支点。若一个圆包装中存在至少一个分支点,则我们称该圆包装为分支圆包装。显然,图2中的南北极点都是P中分支度为1的分支点。

取一固定的加权三角剖分曲面 ( X , T , Φ ) ,考虑作用于边集E上的边长函数为 l : E ( 0 , + ) 。若对于任意三角形 Δ i j k F ,三条边长 l i j , l j k , l i k 都满足三角不等式,则称l为一个分片线性度量。Thurston (见引理13.7.2 [1])指出,任意一个圆包装度量r都可以如下方式决定分片线性度量l。在 E 2 中,令

l i j = r i 2 + r j 2 + 2 r i r j cos ( Φ ( e i j ) ) ,

2 中,令

l i j = cosh 1 ( cosh r i cosh r j + sinh r i sinh r j cos ( Φ ( e i j ) ) ) .

记点 v i 处的内角和为 θ ( v i ) ,则 v i 处的组合高斯曲率定义为

K i = 2 π θ ( v i ) , (1)

K = ( K 1 , K 2 , , K N )

定义2.1取一固定的加权三角剖分曲面 ( X , T , Φ ) ,若点 { b j | j = 1 , , m } 满足 θ ( b j ) = ( 1 + β j ) 2 π ,且任取 v i n t ( T ) \ { b 1 , b 2 , , b m } ,均有 θ ( v ) = 2 π ,则称 b j 是一个分支度为 β j 的分支点。记所有分支点为 { b j | j = 1 , , m } ,且其对应的分支度分别为 β 1 , , β m ,则记P的分支集为

b r ( P ) = { ( b 1 , β 1 ) , ( b 2 , β 2 ) , , ( b m , β m ) } ,

称P为实现三角剖分 ( T , Φ ) 且具有分支集 b r ( P ) 的分支圆包装。设圆包装P实现三角剖分 ( T , Φ ) ,若存在一圆包装度量 r : V ( 0 , + ) 满足在任一点 v i 处组合高斯曲率(1)均满足下述等式

K i ( r ) + 2 β i π = 0 , (2)

则称该圆包装度量为分支圆包装度量,且记为 r b

为简化计算,若内点 v i 不是分支点,可认为它是一个分支度为0的分支点。

定义2.2分支集 b r ( P ) = { ( b 1 , β 1 ) , ( b 2 , β 2 ) , , ( b m , β m ) } V × N 若满足下述条件,则称 b r ( P ) 是一个分支结构(branch structure):存在一个包含至少一个内点的简单闭路径 Γ = { e 1 , e 2 , , e n } 满足不等式

j = 1 n [ π Φ ( e j ) ] > 2 ( l ( Γ ) + 1 ) π ,

这里 l ( Γ ) = i = 1 m δ i ( Γ ) β i 且若 Γ 包含点 b i ,则 δ i = 1 ,否则 δ i = 0

Dubejko [5] 证明了当且仅当 b r ( P ) 是加权三角剖分 ( T , Φ ) 的一个分支结构时,实现加权三角剖分 ( T , Φ ) 的分支圆包装存在。

2.2. 分支组合Calabi流

E 2 2 中,Chow-罗锋 [6] 首次提出的组合Ricci流方程分别如下:

d r i d t = K i r i , (3)

d r i d t = K i sinh ( r i ) . (4)

Chow-罗锋证明了组合Ricci流的解长时间存在并以指数速度收敛。在Chow-罗锋研究的基础上,蓝师义和戴道清引入分支结构,提出了分支组合Ricci流,并将Chow-罗锋的结果成功推广至分支情形。在 E 2 2 中的分支组合Ricci流方程分别如下:

d r i d t = ( K i + 2 β i π ) r i , (5)

d r i d t = ( K i + 2 β i π ) sinh r i . (6)

葛化彬 [8] [9] 提出了组合Calabi流的概念,给出 E 2 2 中的组合Calabi流方程分别如下:

d r i d t = Δ K i r i , (7)

d r i d t = Δ K i sinh r i , (8)

这里 Δ 是离散拉普拉斯算子。

受蓝师义和戴道清的分支组合Ricci流启发,可定义分支组合Calabi流。

定义2.3取一固定的加权三角剖分曲面 ( X , T , Φ ) ,其分支结构为 b r ( P ) = { ( b 1 , β 1 ) , ( b 2 , β 2 ) , , ( b m , β m ) } V × N ,则 E 2 中分支组合Calabi流定义为

d r i d t = Δ ( K i + 2 β i π ) r i , (9)

2 中的分支组合Calabi流为

d r i d t = Δ ( K i + 2 β i π ) sinh r i . (10)

当不含分支点即任意 β i = 0 时,则方程(9)和(10)为葛化彬提出的组合Calabi流(7)和(8)。

2.3. 主要结果

下面是本文的主要结果,在 E 2 2 中,我们成功将葛化彬等人的结果 [8] [11] 推广至分支情形下,得到了分支组合Calabi流解的长时间存在性与收敛性。

定理2.4.在 E 2 2 中给定一个固定的加权三角剖分曲面 ( X , T , Φ ) ,其分支结构为 b r ( P ) = { ( b 1 , β 1 ) , ( b 2 , β 2 ) , , ( b m , β m ) } V × N ,任取初值 r ( 0 ) = ( r 1 ( 0 ) , , r N ( 0 ) ) > 0 N 满足 i = 1 N r i ( 0 ) = 1 ,分支组合Calabi流(9)和(10)的解 { r ( t ) } 长时间存在,且当且仅当分支圆包装度量 r b 存在时,解 { r ( t ) } 收敛。若解收敛,则必指数收敛。

3. 离散拉普拉斯算子

E 2 中,令 u i = ln r i ,在 2 中,令 u i = ln tanh r i / 2 ,则分支组合Calabi流方程(9),(10)可统一写为

d u i d t = Δ ( K i + 2 β i π ) . (11)

若点 v i v j 相邻,则记作 i ~ j ,因此任意边 e i j 均可表示为 i ~ j 。在 E 2 中,对任意边 i ~ j ,设

B i j = Δ i j k F θ i j k r j r j , (12)

2 中,对任意边 i ~ j ,设

B i j = Δ i j k F θ i j k r j sinh r j , (13)

A i = sinh r i r i ( Δ i j k F A r e a ( Δ i j k ) ) , (14)

这里 θ i j k 是三角形 Δ i j k 中点 v i 处的内角, A r e a ( Δ i j k ) 为三角形 Δ i j k 的面积, v k v l 是两个顶点且三角形 Δ i j k Δ i j l 具有公共边 e i j

事实上,在 E 2 2 中,均有 B i j = B j i 成立,因为 θ i j k r j r j = θ j i k r i r i 以及 θ i j k r j sinh r j = θ j i k r i sinh r i (见 [6] 中引理2.3)。设 A = d i a g { A 1 , A 2 , , A N } ,定义 L B = ( ( L B ) i j ) N × N 如下:

( L B ) i j = { k i B i k , i = j B i j , i ~ j 0 , (15)

K ( r ) + 2 β π 的雅可比矩阵为L,由于 L = ( L i j ) = ( ( K i + 2 β i π ) u j ) = ( K i u j ) ,由 [9] 中引理3.3可得下述引理。

引理3.1在 2 中,L正定;在 E 2 中,L是一个 N × N 维的半正定矩阵,它的秩为 N 1 。此外,L的零空间为 K e r ( L ) = { t ( 1 , , 1 ) T | t }

葛化彬在 [11] 中指出在 E 2 中, L = L B ,在 2 中, L = L B + A ,此外, L B 半正定,A是正定的。

设f是作用在顶点集V上的任意函数,它可看作一个N维向量,则离散拉普拉斯算子Δ定义为

Δ f i = ( Δ f ) i = ( L T f ) i = j = 1 N ( K i + 2 β i π ) u j f j . (16)

葛化彬在 [9] 中指出,在 E 2 中,

Δ f i = j ~ i B i j ( f j f i ) , (17)

此外,葛化彬和徐旭在 [11] 中指出,在 2 中,

Δ f i = j ~ i B i j ( f j f i ) A i f i . (18)

3.1. 对离散拉普拉斯算子的估计

这里对 B i j A i 进行估计,我们得到下述引理。

引理3.2参考引理3.1 [11] 和引理3.1 [9] 在 E 2 中,任取 e i j ,其中 1 i , j N ,可得

0 < B i j < 2 3 , (19)

2 中,可得

0 < B i j < 1 ,

0 < A i < d i cosh 1 d cosh 1 ,

这里 d = max 1 i N { d i } d i 是点 v i 的度,即与 v i 相邻的点的个数。

此外,葛化彬在 [10] 中给出的引理2.2和2.3是证明方程解的长时间存在性和收敛性的关键步骤,我们可推广他的结果,得到下述两个引理。

引理3.3在 2 中,给定一个固定的加权三角剖分曲面 ( X , T , Φ ) 。任选双曲三角形 Δ i j k ,其三个顶点处圆的半径分别为 r i , r j , r k ,对于任意常数 c 1 1 ,均存在一个充分大的常数 C = C ( c 1 ) > 0 使得当 r i C 时,有

r i [ ( 1 + c 1 ) A r e a ( Δ i j k ) + c 1 θ i j k ] 0.

基于引理3.3,我们可以得到下述引理。

引理3.4对于任意常数 c 1 1 ,均存在一个充分大的常数 C = C ( c 1 ) > 0 使得当 r i C 时,有

A i c 1 j i B i j .

c 1 = 1 ,则引理3.3和3.4即为葛化彬 [10] 中的引理2.2和2.3。因此,我们的结果将 [10] 中引理2.2和2.3从 c 1 = 1 推广至任意常数 c 1 1

下面我们将给出上述两个引理的证明。

证明[引理3.3]:取一点 v ¯ i ,使得 r ¯ i r i 充分接近, r ¯ i > r i 且点 v ¯ i 与点 v i 位于边 e i j 的同侧。在 [11] 中,葛化彬指出当圆包装度量 r i 增加时,对应点 v i 的内角 θ i 会增加,显然 θ j i ¯ k > θ j i k ,所以点 v i 处于双曲三角形 Δ i ¯ j k 的内部,过点 v i v i S e i ¯ j ,垂足为点S,详情可参考图3

Figure 3. A hyperbolic triangle for lemma 3.3

图3. 引理3.3的一个双曲三角形

证明该引理等价于证明对于任意 r ¯ i > r i ,下述不等式成立。

( 1 + c 1 ) A r e a ( Δ i ¯ j k ) + c 1 θ i ¯ j k ( 1 + c 1 ) A r e a ( Δ i j k ) c 1 θ i j k 0.

2 中的Gauss-Bonnet公式得,在双曲三角形 Δ i j k 中有 θ i j k + θ j i k + θ k i j = π A r e a ( Δ i j k )

x = θ j i ¯ k θ j i k y = θ k i ¯ j θ k i j ,可得

( 1 + c 1 ) A r e a ( Δ i ¯ j k ) + c 1 θ i ¯ j k ( 1 + c 1 ) A r e a ( Δ i j k ) c 1 θ i j k = A r e a ( Δ i ¯ j k ) A r e a ( Δ i j k ) c 1 x c 1 y = A r e a ( Δ i ¯ i j ) + A r e a ( Δ i ¯ i k ) c 1 x c 1 y = A r e a ( Δ i ¯ i j ) c 1 x + A r e a ( Δ i ¯ i k ) c 1 y . (20)

所以我们仅需证明 A r e a ( i ¯ i j ) c 1 x 0 即可。由于 r ¯ i > r i ,可得 l i ¯ j > l i j θ i i ¯ j > θ i ¯ i j 。考虑到 θ i i ¯ j + θ i ¯ i j + x = π A r e a ( Δ i ¯ i j ) < π ,因此 θ i ¯ i j < π 2 。由于 r ¯ i 充分接近于 r i ,所以 l i ¯ i 1 。由双曲余弦定理得 cos x 的极限存在,即

cos x = cosh l i j cosh l i ¯ j cosh l i ¯ i sinh l i j sinh l i ¯ j 1 , r i + .

因此可以找到一个正常数 C 1 使得当 r i C 1 时,有 x π 4 ( 1 + c 1 ) 成立。

γ = θ v i S v j ,当 γ < π 4 时,显然有 γ + ( 1 + c 1 ) x < π 2 ,然而由 2 中的Gauss-Bonnet公式得

γ + x + π 2 + A r e a ( Δ S v i v j ) = π ,

因此 A r e a ( Δ i ¯ i j ) > A r e a ( Δ S v i v j ) > c 1 x

γ π 4 时, sinh l S v j sinh l i j = sin γ 2 2 。由于 cos γ = sin x cosh l S v j cos x = tanh l S v j tanh l i j 并且 sin sin γ = sinh l S v j sinh l i j ,然后由 [10] 得

sin ( A r e a ( Δ S v i v j ) ) sin x cosh r i 1 2 . (21)

考虑到当 x > 0 时, cosh x 递增且无上界,因此这里我们可以选取一个常数 C 2 使得当 r i C 2 时有 cosh r i 1 2 c 1 。因此可得

sin ( A r e a ( Δ S v i v j ) ) c 1 sin x sin c 1 x ,

因为x和 c 1 x 都属于 ( 0 , π 2 ) c 1 1 ,所以右侧不等式成立。注意到 A r e a ( Δ S v i v j ) ( 0 , π 2 ) ,因此可得

A r e a ( Δ i ¯ i j ) > A r e a ( Δ S v i v j ) > c 1 x ,

C = max { C 1 , C 2 } ,则引理3.3证毕。

下面我们将证明引理3.4。

证明[引理3.4]:

A i c 1 j ~ i B i j = Δ i j k F A r e a ( Δ i j k ) r i sinh r i c 1 j ~ i ( θ i j k r j sinh r j + θ i j l r j sinh r j ) = Δ i j k F A r e a ( Δ i j k ) r i sinh r i c 1 Δ i j k F ( θ j i k r i sinh r i + θ k i j r i sinh r i ) = Δ i j k F ( A r e a ( Δ i j k ) c 1 θ j i k c 1 θ k i j ) r i sinh r i = Δ i j k F [ ( 1 + c 1 ) A r e a ( Δ i j k ) + c 1 θ i j k ] r i sinh r i 0. (22)

证毕。

定义3.5在 E 2 2 中,给定一个固定的加权三角剖分曲面 ( X , T , Φ ) ,其分支结构为 b r ( P ) = { ( b 1 , β 1 ) , ( b 2 , β 2 ) , , ( b m , β m ) } V × N ,定义分支组合Calabi能量E为

E ( r ) = i = 1 N ( K i ( r ) + 2 β i π ) . (23)

3.2. 分支组合Ricci势

E 2 2 中,因为 ( K i + 2 β i π ) u j = ( K j + 2 β j π ) u i (参考引理2.3 [6]),所以光滑1-形式 i = 1 N ( K i ( u ) + 2 β i π ) d u i 是闭的。蓝师义–戴道清 [16] 定义分支组合Ricci势如下。

定义3.6给定一个固定的加权三角剖分曲面 ( X , T , Φ ) ,其分支结构为 b r ( P ) = { ( b 1 , β 1 ) , ( b 2 , β 2 ) , , ( b m , β m ) } V × N 。在 E 2 中,定义分支组合Ricci势为:

F ( u ) = a u i = 1 N ( K i ( u ) + 2 β i π ) d u i , u N , (24)

2 中,定义分支组合Ricci势为

F ( u ) = a u i = 1 N ( K i ( u ) + 2 β i π ) d u i , u < 0 N , (25)

这里 a N 是任意常数。

引理3.7在 E 2 中,若分支组合Calabi流方程的解 { r ( t ) } 存在,则 i = 1 N r i ( t ) = i = 1 N r i ( 0 ) 以及 i = 1 N u i ( t ) = i = 1 N u i ( 0 ) 均是常数。

由 [16] 中定理1,可得分支组合Ricci势的凸性:

引理3.8在 2 中,分支组合Ricci势F在 ( , 0 ) N 上严格凸。在 E 2 中,F在 N 上凸且在超平面 U = { ( u 1 ( t ) , u 2 ( t ) , , u N ( t ) ) T | i = 1 N u i ( t ) = 0 } (或超曲面 P = { ( r 1 ( t ) , r 2 ( t ) , , r N ( t ) ) T | i = 1 N r i ( t ) = 1 } )上严格凸。

由引理3.8可得下述刚性定理,参考文献 [4] [5] [16]。

引理3.9在 E 2 2 中,记圆包装度量 ( r 1 ( t ) , r 2 ( t ) , , r N ( t ) ) T 到组合高斯曲率 ( K 1 ( t ) , K 2 ( t ) , , K N ( t ) ) T 的映射为 Π : > 0 N N 。则

2 中,Π是单射,所以组合高斯曲率唯一决定圆包装度量。

E 2 中,在超平面 U = { ( u 1 ( t ) , u 2 ( t ) , , u N ( t ) ) T | i = 1 N u i ( t ) = 0 } 上,Π是单射,所以在不考虑刚性运动下组合高斯曲率唯一决定圆包装度量,即K与aK决定同一圆包装度量,这里a是任意非零常数。

葛化彬(参考定理B.2 [9])和葛化彬–徐旭(参考引理B.1 [11])证明了组合Ricci势在 E 2 2 中的真性,这里我们可将他们的结果推广到分支情形。

引理3.10在 2 中,分支组合Ricci势F(25)满足

lim u + F ( u ) = + , (26)

E 2 中,分支组合Ricci势F(24)满足

lim u + F ( u ) | u U = + . (27)

由于引理B.1 [11] 得 F * ( u ) = a u i = 1 N K i ( u ) d u i 满足真性,显然有引理3.10成立。

4. 分支组合Calabi流解的长时间存在性

在这章中,我们将证明在 E 2 2 中,分支组合Calabi流的解 { r ( t ) } 长时间存在。

定理4.1在 E 2 2 中,给定一个固定的加权三角剖分曲面 ( X , T , Φ ) ,其分支结构为 b r ( P ) = { ( b 1 , β 1 ) , ( b 2 , β 2 ) , , ( b m , β m ) } V × N 。任给初值圆包装度量 r ( 0 ) ,分支组合Calabi流的解 { r ( t ) } 长时间存在。

为了证明该定理,我们需要下面的引理(参考引理3.5 [6])。

引理4.2对于任意的 ϵ > 0 ,均存在一个充分大的常数l,使得当 r i > l 时,双曲三角形 Δ i j k 中点 v i 处的内角 θ i j k

θ i j k < ϵ .

现在我们开始证明定理4.1。

证明:设 d i 为点 v i 处的度,即与 v i 相邻的点的个数。记 d = max { d 1 , , d N } ,则 ( 2 d ) π < K i < 2 π 。在 E 2 中,

Δ ( K i + 2 β i π ) = [ L T ( K + 2 β π ) ] i = j = 1 N L i j ( K j + 2 β j π ) = j i ( B i j ) ( K j + 2 β j π ) L i i ( K i + 2 β i π ) = j i B i j ( K j K i + 2 β j π 2 β i π ) . (28)

由引理3.2得 0 < B i j < 2 3 ,可得 | Δ ( K i + 2 β i π ) | 上有界,即

| Δ ( K i + 2 β i π ) | < 2 3 N ( d + 2 β max ) π c ,

这里 β max = max { β 1 , β 2 , , β N }

因此 | d u i d t | < c ,进而我们得到解 { r ( t ) } 满足

c 2 e c t r i c 3 e c t ,

这里 c 2 = min { r 1 ( 0 ) , r 2 ( 0 ) , , r N ( 0 ) } c 3 = max { r 1 ( 0 ) , r 2 ( 0 ) , , r N ( 0 ) } ,因此在 E 2 中分支组合Calabi流(9)的解 { r ( t ) } 长时间存在。

2 中,根据引理4.2,可得对于任意点 v i ,存在一个充分大的 l > 0 ,使得当 r i > l 时,三角形 Δ i j k 中点 v i 处的内角 θ i j k 小于 π d ,因此 K i > π 。令(18)中 f = K + 2 β π ,则

1 sinh r i d r i d t = j ~ i B i j ( K j K i + 2 β j π 2 β i π ) A i ( K i + 2 β i π ) < j ~ i B i j ( 2 π K i + 2 β max π 2 β i π ) A i ( K i + 2 β i π ) = ( 2 π + 2 β max π ) j ~ i B i j ( j ~ i B i j + A i ) ( K i + 2 β i π ) ( 2 π + 4 β max π ) j ~ i B i j ( π + 2 β min π ) ( j ~ i B i j + A i ) = ( π + 2 β min π ) [ 1 + 4 β max 2 β min 1 + 2 β min j ~ i B i j A i ] . (29)

令引理3.4中 c 1 = 1 + 4 β max 2 β min 1 + 2 β min 1 Q = max { l , C , r i ( 0 ) + 1 } ,则当 r i > Q 时,有 1 sinh r i d r i d t < 0 ,所以解 { r ( t ) } 有上界。

下面我们将证明解 { r ( t ) } 有下界。这里我们需要考虑分支组合Ricci势 F ( u ) ,对 F ( u ) 求导得 d F ( u ) d t = ( K + 2 β π ) ( L T ) ( K + 2 β π ) 0 ,因此 F ( u ) 沿着分支组合Calabi流是单调不增的。又因为引理3.10指出 lim u + F ( u ) = + ,所以存在一个仅依赖于三角剖分,分支结构与初值 r ( 0 ) 的正常数 δ ,使得对于任意点 v i 和时间t,均有 u i ( t ) δ ,因此解 r i ( t ) 有下界,即

r i ln 1 + e δ 1 e δ .

综上,在 2 中分支组合Calabi流的解 { r ( t ) } 长时间存在,证毕。

5. 分支组合Calabi流解的收敛性

在这一章中,我们将证明在 E 2 2 中分支组合Calabi流(9)和(10)解 { r ( t ) | t [ 0 , + ) } 的收敛性。

定理5.1在 2 E 2 中,给定一个固定的加权三角剖分曲面 ( X , T , Φ ) ,其分支结构为 b r ( P ) = { ( b 1 , β 1 ) , ( b 2 , β 2 ) , , ( b m , β m ) } V × N 。任取一初值圆包装度量 r ( 0 ) = ( r 1 ( 0 ) , , r N ( 0 ) ) > 0 N 满足 i = 1 N r i ( t ) = 1 ,设 { r ( t ) | t [ 0 , + ) } 是分支组合Calabi流长时间存在的解,则当且仅当存在一个分支圆包装度量 r b 时,解 { r ( t ) } 收敛。若解收敛,则必指数收敛。

为了证明定理5.1,我们需要下面的引理。

引理5.2在 E 2 2 中,给定一个固定的加权三角剖分曲面 ( X , T , Φ ) ,其分支结构为 b r ( P ) = { ( b 1 , β 1 ) , ( b 2 , β 2 ) , , ( b m , β m ) } V × N 。设存在一个分支圆包装度量 r b ,分支组合Calabi流的解 { r ( t ) } M ,这里M是 > 0 N 中的一紧集,则解 { r ( t ) } 指数收敛到分支圆包装度量 r b

下面我们将给出引理5.2的证明。

证明:对 K i + 2 β i π 求导,得到

d ( K + 2 β π ) d t = d ( K + 2 β π ) d u d u d t = L 2 ( K + 2 β π ) .

因为L在 2 中是正定的,在 E 2 中半正定,所以我们可取L的最小正特征值为 λ 1 ,然后对分支组合Calabi能量进行求导得到

E ( t ) = 2 ( K + 2 β π ) T ( K + 2 β π ) = 2 ( K + 2 β π ) T L 2 ( K + 2 β π ) 2 λ 1 2 ( K + 2 β π ) T ( K + 2 β π ) = 2 λ 1 2 E ( t ) .

因此我们得到 E ( t ) E 0 e 2 λ 1 2 t | K i ( t ) + 2 β i π | | K ( t ) + 2 β π | = E ( t ) E 0 e λ 1 2 t 。又因为 { r ( t ) } M ,M是 > 0 N 中一紧集,所以沿着分支组合Calabi流, L i j r i sinh r i 是有界的。因此在 E 2 中,我们有

| d r i d t | = | r i Δ ( K i + 2 β i π ) | = | r i j = 1 N L i j ( K j + 2 β j π ) | m e λ 1 2 t ,

2 中,可得

| d r i d t | = | sinh r i j = 1 N L i j ( K j + 2 β j π ) | n e λ 1 2 t ,

这里 m , n 是两个充分大的正常数。这说明,无论是在 E 2 还是 2 中,分支组合Calabi流的解 { r ( t ) } 指数收敛。记 r ( t ) r * , t + 以及 K ( t ) K * , t + ,显然 K * + 2 β π = 0 ,所以由引理3.9得 r * 实际上是一个分支圆包装度量 r b ,证毕。

下面我们将给出定理5.1的证明。

证明:在 E 2 中,根据引理5.2的证明,容易得到,若分支组合Calabi流的解 { r ( t ) | t [ 0 , + ) } 收敛,则分支圆包装度量 r b 存在。

另一方面,假设分支圆包装度量 r b 存在,设 { r ( t ) | t [ 0 , + ) } 是分支组合Calabi流的解,令 u i ( t ) = ln r i ( t ) ,由 i = 1 N r i ( 0 ) = 1 得初值圆包装度量 u 0 在超平面 U = { ( u 1 ( t ) , u 2 ( t ) , , u N ( t ) ) T | i = 1 N u i ( t ) = 0 } 上,由引理3.7可得解 { u ( t ) | t [ 0 , + ) } U 。令分支组合Ricci势定义中的 a = u b ,由引理3.9得分支组合Ricci势的梯度映射 F : u K + 2 β π 是单射,并且 u = u b 是F唯一的临界点。由引理3.1得F的海塞矩阵L是半正定的,所以有 F ( u ) F ( u b ) = 0 。进一步地,我们还可以得到

d F d t = i = 1 N ( K i + 2 β i π ) d u i d t = ( K + 2 β π ) T L T ( K + β π ) 0.

所以F沿着分支组合Calabi流是单调递减的,即 0 = F ( u b ) F ( u ) F ( u ( 0 ) ) ,由引理3.10,可得 lim u + F ( u ) | u U = + ,因此 { u ( t ) } ( F | u U ) 1 [ 0 , F ( u ( 0 ) ) ] ( F | u U ) 1 [ 0 , F ( u ( 0 ) ) ] 是一个紧集,因此 { r ( t ) } M ,这里M是 > 0 N 中的一紧集,所以由引理5.2,可得分支组合Calabi流方程解 { r ( t ) } 指数收敛到分支圆包装度量 r b

2 中,证明与在 E 2 中完全类似,此处省略。证毕。

6. 结论

分支组合Calabi流是四阶曲率流,无法建立类似Chow和罗锋 [6] 在组合Ricci流证明中的极大原理,这是研究的一大难点。我们推广了葛化彬等人 [8] [9] [10] [11] 在组合Calabi流中使用的几何方法,克服了难点,证明了分支组合Calabi流解的长时间存在性与收敛性。

引理3.3和引理3.4中的一致估计是分支组合Calabi流解的长时间存在性与收敛性证明的关键步骤,为此我们将葛化彬和华波波 [10] 中的类似估计从 c 1 = 1 成功推广到了任意常数 c 1 1

综上所述,本文对于组合曲率流的理论研究具有一定的借鉴价值,同时还有潜在的应用价值。

基金项目

国家自然科学基金面上项目(12171480),湖南省杰出青年基金项目(2022JJ10059),湖南省自然科学基金面上项目(2020JJ4658)资助。

NOTES

*通讯作者。

参考文献

[1] Thurston, W. (1976) Geometry and Topology of 3-Mainfolds. Princeton Lecture Notes.
[2] Rodin, B. and Sullivan, D. (1987) The Convergence of Circle Packings to the Riemann Mapping. Journal of Differential Geometry, 26, 349-360.
https://doi.org/10.4310/jdg/1214441375
[3] Beardon, A.F. and Stephenson, K. (1990) The Uniformization Theorem for Circle Packings. Indiana University Mathematics Journal, 39, 1383-1425.
https://doi.org/10.1512/iumj.1990.39.39062
[4] Bowers, P. and Stephenson, K. (1996) A Branched Andreev-Thurston Theorem for Circle Packings of the Spheres. Proceedings of the London Mathematical Society, 73, 185-215.
https://doi.org/10.1112/plms/s3-73.1.185
[5] Dubejko, T. (1995) Branched Circle Packings and Discrete Blaschke Products. Transactions of the AMS, 347, 4073- 4103.
https://doi.org/10.1090/S0002-9947-1995-1308008-1
[6] Chow, B. and Luo, F. (2003) Combinatorial Ricci Flows on Surfaces. Journal of Differential Geometry, 63, 97-129.
https://doi.org/10.4310/jdg/1080835659
[7] Feng, K., Lin, A.J. and Zhang, X.X. (2020) Combinatorial p-th Calabi Flows for Discrete Conformal Factors on Surfaces. Journal of Geometric Analysis, 30, 3979-3994.
https://doi.org/10.1007/s12220-019-00224-0
[8] 葛化彬. 组合方法与几何方程[D]: [博士学位论文]. 北京: 北京大学, 2012.
[9] Ge, H.B. (2018) Combinatorial Calabi Flows on Surfaces. Transactions of the AMS, 370, 1377-1391.
https://doi.org/10.1090/tran/7196
[10] Ge, H.B. and Hua, B.B. (2018) On Combinatorial Calabi Flow with Hyperbolic Circle Patterns. Advances in Mathematics, 333, 523-538.
https://doi.org/10.1016/j.aim.2018.05.039
[11] Ge, H.B. and Xu, X. (2016) 2-Dimensional Combinatorial Calabi Flow in Hyperbolic Background Geometry. Differential Geometry and Its Applications, 47, 86-98.
https://doi.org/10.1016/j.difgeo.2016.03.011
[12] Lin, A.J. and Zhang, X.X. (2019) Combinatorial p-th Calabi Flows on Surfaces. Advances in Mathematics, 346, 1067-1090.
https://doi.org/10.1016/j.aim.2019.02.011
[13] Luo, F. (2004) Combinatorial Yamabe Flow on Surfaces. Communications in Contemporary Mathematics, 6, 765-780.
https://doi.org/10.1142/S0219199704001501
[14] Bobenko, A., Hoffmann, T. and Springborn, B. (2006) Minimal Surfaces from Circle Patterns: Geometry from Combinatorics. Annals of Mathematics, 164, 231-264.
https://doi.org/10.4007/annals.2006.164.231
[15] Gu, X.D. and Yau, X.T. (2008) Computational Conformal Geometry. Advanced Lectures in Mathematics, Vol. 3, Higher Education Press of China, Beijing.
[16] Lan, S.Y. and Dai, D.Q. (2007) Variational Principles for Branched Circle Patterns. Nonlinear Analysis, 67, 498-511.
https://doi.org/10.1016/j.na.2006.06.008