双曲空间中曲线的逆曲率流研究
Study on Inverse Curvature Flows of Curves in Hyperbolic Space
摘要: 基于Tsai和Wang的线性插值流方法,重点推导双曲空间中凸曲线的逆曲率流下曲线长度和面积的演化方程,引入等周亏格并证明其单调性以刻画流的演化规律;通过构造与内球中心相关的支撑函数,建立曲率的边界估计,最终证明逆曲率流的长时存在性与收敛性。
Abstract: Based on the research method of Tsai and Wang, a linear interpolation flow model is constructed. We focus on deriving the evolution equations of curve length and area under the inverse curvature flow of convex curves in hyperbolic space, introduce the isoperimetric deficit and prove its monotonicity to characterize the evolution law of the flow. By constructing a support function related to the center of the inner sphere, the boundary estimate of curvature is established, and finally the long-time existence and convergence of the inverse curvature flow is proved.
文章引用:方博文. 双曲空间中曲线的逆曲率流研究[J]. 理论数学, 2026, 16(1): 93-104. https://doi.org/10.12677/pm.2026.161012

1. 引言

1.1. 研究背景与动机

曲率流作为几何分析领域的核心研究方向之一,其核心思想是通过曲线或曲面的曲率驱动几何形体演化,最终收敛到具有良好对称性的几何结构(如圆、球面等)。在欧氏空间中,逆曲率流因具有保持凸性、改变几何形态的优良性质,已被广泛应用于等周不等式证明、正则化等问题,详见[1]-[7]。然而,双曲空间作为具有负常曲率的典型非欧空间,其几何结构(如距离公式、曲率定义、等周不等式形式)与欧氏空间存在本质差异,这使得双曲空间中逆曲率流的研究面临独特挑战:一方面,双曲空间的负曲率特性导致几何量演化方程中出现额外的非线性项,增加了分析难度;另一方面,现有欧氏空间中的曲率流方法(如支撑函数技巧、等周亏格单调性分析)无法直接迁移,需针对双曲空间的几何特性进行改进。

线性插值流作为一类融合了保面积流与保长度流优势的非局部曲率流,由Tsai和Wang首次提出并应用于欧氏空间凸曲线的演化研究。该流的核心优势在于通过参数 α[ 0,1 ] 的调节,可灵活平衡长度与面积的演化趋势,进而实现对几何形状的控制。然而,目前关于线性插值流的研究主要集中在欧氏空间,双曲空间中该类流的基本性质(如长时存在性、收敛性)尚未得到系统探讨。此外,双曲空间中凸曲线的等周不等式与欧氏空间存在显著差异,这为等周亏格的单调性分析带来了新的问题。基于上述研究空白,本文旨在将线性插值流的研究推广到双曲空间,通过构建适配双曲几何的支撑函数与辅助工具,揭示该类流的演化规律,为双曲空间中曲率流的研究提供新的思路与方法。

1.2. 与现有研究的关系

现有关于双曲空间曲率流的研究主要集中于两类:一类是保Quermassintegral的曲率流(如Andrews与Wei [8]研究的保高阶几何量的曲率流),其核心是通过调节流的速度函数保持特定几何量不变,进而证明收敛性;另一类是逆曲率流(如Kwong等[9]研究的二维空间形式中的逆曲率流),重点关注曲率的边界估计与流的存在性。与上述研究相比,本文的线性插值流具有以下特点:

1) 速度函数采用线性插值形式,同时融合了保面积与保长度流的特性,可通过参数 α[ 0,1 ] 连续切换两种特殊情形,具有更强的灵活性;

2) 引入的等周亏格定义适配双曲空间的等周不等式,其单调性直接反映了曲线向测地圆收敛的趋势,为收敛性证明提供了关键工具;

3) 支撑函数的构造基于内球中心,充分利用了双曲空间的径向结构,有效简化了曲率估计的难度。

本文的研究成果既是对欧氏空间线性插值流研究的推广,也是对双曲空间曲率流理论的补充,为后续研究更一般的非局部曲率流奠定了基础。

1.3. 问题引入

设常数 α>0 X 0 : S 1 H 2 是光滑凸闭曲线 γ 0 = X 0 ( S 1 ) H 2 的嵌入。

研究如下类型的闭凸曲线的逆曲率流:

{ t X( x,t )=( ϕ( t ) k α ( x,t ) )ν( x,t ) X( ,0 )= X 0 ( ) (1)

其中 k 是曲线 γ t =X( S 1 ,t ) H 2 的曲率, ν γ t 的单位外法向量, ϕ( t ) 是选定的非局部项,使得相应的面积或长度演化保持稳定。特别地,当 ϕ( t ) 满足特定条件时,流(1)保持曲线所围区域的面积 A( t ) 或长度 L( t ) 。更精确地说,如果

ϕ( t )= 1 L( t ) γ t k α ( x,t )ds

流(1)保持 γ t 围成区域 Ω t 的面积 A( t ) 。如果

ϕ( t )= 1 2π γ t k 1α ( x,t )ds

流(1)保持 γ t 的长度 L( t )

在此,借助Tsai和Wang在文[10]中的方法研究这类线性插值流:

t X( x,t )=( λ 2π k 1α ds+ 1λ L( t ) k α ds k α )ν( x,t ) (2)

其中 0λ1,α>0 。为了简化记号,参考文献[11]

φ= λ 2π k 1α ds + 1λ L( t ) k α ds k α

则线性插值流可表示为:

X( x,t ) t =φ(x,t)ν( x,t ) (3)

定理 1.1 沿着流(3),曲线的长度 L 单调递减,面积 A 单调递增,等周亏格 Δ(t) 单调递减

证明 L A 分别表示闭曲线 X 的长度和围成的代数面积:

L= X ds,A= X X,ν ds

dL dt ( t )= X W,kν ds, dA dt ( t )= X W,ν ds

其中 , 表示内积, W= t X X 的速度向量。则 L( t ) A( t ) 的演化方程:

dL dt = φkds , dA dt = φds (4)

进一步推导得:

dL dt =λ k 1α ds 2π( 1λ ) L( t ) k α ds+ k 1α ds =( 1λ )( k 1α ds 2π L k α ds )

dA dt = λL 2π k 1α ds( 1λ ) k α ds + k α ds =λ( k α ds L 2π k 1α ds )

利用Hölder不等式可得 dL dt 0, dA dt 0

等周亏格:等周亏格的本质是衡量曲线偏离“等周不等式最优解”的程度。

参考文献[12],定义等周亏格:

Δ( t )=L ( t ) 2 4πA( t )A ( t ) 2

结合 L( t ) A( t ) 的单调性可知, Δ( t ) 单调递减。

2. 几何量的演化方程

引理 2.1 曲线各个几何量的演化方程如下:

(1) t ds=k( ϕ( t ) k α )ds

(2) t ν= s ( k α )T

(3) k t = s 2 k α +( k 2 1 )( ϕ( t ) k α )

证明

(1)

ds=ϑdx,ϑ=| X x |= X x , X x

t ϑ= 1 ϑ 2 X tx , X x

X t =( ϕ( t ) k α )ν

2 X tx = x ( ϕ(t) k α )ν+( ϕ( t ) k α ) ν x

由于 ϕ( t ) 是全局时间函数, ϕ( t ) x =0 ,且 ν x = 1 ϑ s ν= k ϑ T

2 X tx , X x =( ϕ( t )k α ) k ϑ ϑ 2 =kϑ( ϕ( t ) k α )

t ϑ= 1 ϑ ϑ( ϕ( t )k α ) k ϑ = k ϑ ( ϕ( t ) k α )

t ds= t ϑdx=k( ϕ( t ) k α )ds

(2)

ν,T =0, t ν,T + ν, t T =0

t T= t ( X s )= s ( X T ) = s ( ϕ( t ) k α )ν+( ϕ( t ) k α ) s ν = s ( ϕ( t ) k α )ν+k( ϕ( t ) k α )T

代入正交条件:

t ν,T + ν, s ( ϕ( t ) k α )ν+k( ϕ( t ) k α )T =0

t ν,T + s ( ϕ( t ) k α )=0

t ν,T = s ( ϕ( t ) k α )

t ν 垂直于 ν ,得 t ν 沿 T 方向

ϕ( t ) 是不依赖空间坐标的全局时间函数,即 s ϕ( t )=0

t ν= t ν,T T= s ( ϕ( t ) k α )T= s ( k α )T=α k α1 s kT

(3) 根据曲率定义:

s T=kν

t s T= k t νk t ν

s t T= s [ s ( ϕ( t ) k α )ν+k( ϕ( t ) k α )T ]

展开得:

s 2 ( ϕ( t ) k α )ν+ s ( ϕ( t ) k α ) s ν+ s [ k( ϕ( t ) k α ) ]T+k( ϕ( t ) k α ) s T

代入 s ν=kT s T=kν

s 2 ( ϕ( t ) k α )ν+k s ( ϕ( t ) k α )T+ s [ k( ϕ( t ) k α ) ]T k 2 ( ϕ( t ) k α )ν

比较左右两边 ν T 的分量:

ν 的分量:

s 2 ( ϕ( t ) k α ) k 2 ( ϕ( t ) k α )= k t (5)

T 的分量:

k s ( ϕ( t ) k α )+ s [ k( ϕ( t ) k α ) ]=k s k α T

s t T= k t νk t ν= k t νk s k α T

ϕ( t ) 不依赖于s, s 2 ϕ( t )=0

k t = s 2 k α +( k 2 1 )( ϕ( t ) k α )

3. 支撑函数

支撑函数:支撑函数是刻画凸曲线几何形态的重要工具,其核心思想是通过曲线到固定点的距离在法向量方向的投影描述曲线的位置与形状。在双曲空间中,由于距离公式的非线性,支撑函数的构造需结合径向结构,本文选择内球中心作为固定点,使得支撑函数充分利用曲线的凸性,简化演化方程的推导。

引理 3.1 p 0 γ t 0 内球 B( p 0 , ρ 0 ) 的中心,定义支撑函数:

u( x,t )= sinh r p 0 ( x ) r ,ν

其中 r p 0 ( x ) 是从点x到固定点 p 0 的距离,其中 x γ t ,t[ t 0 ,min{ T, t 0 +τ } ] ,满足: c 1 u( x,t ) c 2

2C:=sinh( ρ 0 2 )u( x,t )sinh( 2 C 2 )

则支撑函数的演化方程

t u=α k α1 u ss +cosh r p 0 ( x )( ϕ( t )( 1α ) k α )α k α2 u

证明

t u= t ( sinh r p 0 ( x ) r ),ν + sinh r p 0 ( x ) r , t ν

1) 时间导数展开

t u= t ( sinh r p 0 ) r,ν +sinh r p 0 t r,ν (6)

第一项: t sinh r p 0 ( x )

r p 0 ( x ) t = x t , r = ( ϕ( t ) k α )ν, r

t sinh r p 0 ( x )=cosh r p 0 ( x ) ( ϕ( t ) k α )ν, r =cosh r p 0 ( ϕ( t ) k α ) u sinh r p 0

第二项: t r ,ν

¯ ν r 为径向联结在法向量方向的分量, r 是径向向量

¯ ν r ,ν =0

t r ,ν = ¯ x t r ,ν + r , ¯ x t ν = r , ¯ ν ν( ϕ( t ) k α ) = r ,kT ( ϕ( t ) k α ) =( ϕ( t ) k α ) u s k α sinh r p 0

将两项合并:

u t =cosh r p 0 ( x ) ν, r ( ϕ( t ) k α )+sinh r p 0 ( x )k r ,T ( ϕ( t ) k α )

2) 支撑函数的空间导数

u s = s ( sinh r p 0 r ,ν ) =cosh r p 0 T, r r ,ν +sinh r p 0 r , s ν =cosh r p 0 T, r r ,ν +sinh r p 0 k r ,T =ksinh r p 0 r ,T

3) 支撑函数的二阶空间导数

u ss = s ( u s ) = s [ cosh r p 0 T, r r ,ν +sinh r p 0 k r ,T ] =kcosh r p 0 k 2 u+ k s sinh r p 0 r ,T

联立得:

u t =cosh r p 0 ( ϕ( t ) k α ) u sinh r p 0 +sinh r p 0 ( ϕ( t ) k α ) u s k α sinh r p 0 =coth r p 0 ( ϕ( t ) k α )u+ 1 k α ( ϕ( t ) k α ) u s

利用 u ss 的表达式消去 u s

u t =α k α1 u ss +cosh r p 0 ( x )( ϕ( t )( 1α ) k α )α k α2 u

4. 流的存在性

命题 4.1 γ t 是流 t X=( ϕ( t ) k α )N 的光滑凸解,则存在仅依赖于初始曲线 γ 0 和参数 α 的正整数 k 0 ¯ ,使得对所有 t[ 0, ) 成立:

maxk k 0 ¯

证明 构造辅助函数 W( x,t )= k u( x,t )c ,其中 c=sinh( ρ 2 ) 为常数。

由于 u2c ,分母恒正,W定义良好

k=W( uc )Wsinh( 2 c 2 )

目标通过W的上界控制k的上界

曲率的演化方程:

k t = s 2 k α +( k 2 1 )( ϕ( t ) k α ) (7)

支撑函数的演化方程:

u t =α k α1 u ss +cosh r p 0 ( ϕ( t )( 1α ) k α )α k α2 u (8)

k t u t 的表达式(7)和(8)代入 W t = k t ( uc )k u t ( uc ) 2 得:

W t = 1 uc s 2 k α + ( k 2 1 )( ϕ( t ) k α ) uc + α k α1 u ss cosh r p 0 ( ϕ( t )( 1α ) k α )+α k α2 u ( uc ) 2

假设W ( x 0 , t 0 ) 处取得最大值,则在此点处 W s =0 W ss 0

由一阶导数条件 W s = k s ( uc )k u s ( uc ) 2 =0 k s = k u s uc

由二阶导数条件 W ss 0 k ss ( uc )2 k s u s +k u ss 0

将一阶条件代入二阶条件得:

k ss k s 2 k k u ss uc , u ss k s 2 k k( uc ) k ss k 2

s 2 k α = s ( α k α1 k s )=α k α1 k ss +α( α+1 ) k α2 k s 2 α k α1 ( k s 2 k k u ss uc )+α( α+1 ) k α2 k s 2 =α k α u ss +α k α2 k s 2

由于 usinh( 2 C 2 ) cosh r p 0 cosh( 2 C 2 )

W t α k α u ss +α k α2 k s 2 uc + ( k 2 1 )( ϕ( t ) k α ) uc + α k α1 u ss cosh r p 0 ( ϕ( t )( 1α ) k α )+α k α2 u ( uc ) 2 α k α1 W u ss +( k 2 1 uc cosh( 2 C 2 ) uc )ϕ( t )+ c 1 W 2 + c 2 W

其中 C 1 C 2 为正的常数

ϕ( t ) 的非局部估计

(1) 保面积

ϕ( t )= 1 L( t ) γ t k α ds 1 L( t ) γ t ( uc ) Wdssinh( 2 c 2 ) W ˜ , W ˜ = max γ t W (9)

(2) 保长度

ϕ( t )= 1 2π γ t k 1α ds 1 2π γ t ( uc ) Wkdssinh( 2 C 2 ) W ˜ (10)

统一记为: ϕ( t ) C 3 W ˜ , C 3 =sinh( 2 C 2 )

k=W( uc )Wc 代入 ϕ( t ) C 3 W ˜

d W ˜ dt α ( Wc ) α1 W ˜ +( ( Wc ) 2 1 c cosh( 2 C 2 ) c ) C 3 W ˜ + C 1 W ˜ 2 + C 2 W ˜ (11)

W ˜ 足够大时,主导项为 C 1 W ˜ 2 C 3 W ˜ α+2 ,存在阈值 W 0

使得当 W ˜ > W 0 时, d W ˜ dt C 4 W ˜ α+2 <0 ,其中 C 4 >0

结合初始条件 W ˜ ( 0 )= K 0 c 有限,由比较定理知 W ˜ ( t ) 有界

W ˜ (t)max{ W ˜ (0), W 0 }

k W ˜ sinh( 2 C 2 ) k ¯ 0 , k ¯ 0 =sinh( 2 C 2 )max{ k 0 c , W 0 }

命题 4.2 初值为 k 0 ( θ ) 的方程(7)的解 k( θ,t ) 长时存在,如果 k 0 ( θ )>0 ,则 k( θ,t )>0

证明 假设存在最小曲率点 k( x 0 , t 0 )=mink ,则 k t 0

结合演化方程:

k t = s 2 k α +( k 2 1 )( ϕ k α )

在最小值点,

s 2 k α 0 ϕ k α

k t ( k 2 1 )( ϕ k α )0

kmin{ k 0 ,1 }>0 ,即k保持凸性

借助文[12]中的方法,双曲空间内星形曲线的单位外法向量 ν 与曲率 k 可由径向函数 ρ( θ,t ) 表示为:

ν= 1 1+ ρ θ 2 / sinh 2 ( ρ ) ( r ρ θ sinh 2 ( ρ ) θ ) (12)

其中 r 为径向单位向量, θ 为辐角方向的单位向量。

曲率的表达式可通过双曲空间的第二基本形式推导得:

κ= sinh 2 ( ρ )cosh( ρ )+2 ρ θ 2 cosh( ρ ) ρ θθ sinh( ρ ) ( sinh 2 ( ρ )+ ρ θ 2 ) 3 2 (13)

将法向量表达式与曲率表达式代入 k α 型流方程,消去 ν k ,结合双曲空间内弧长元的计算公式:

ds= sinh 2 ( ρ )+ ρ θ 2 dθ

将流的方程转化为径向函数 ρ( θ,t ) 的标量抛物方程:

{ ρ t =( ϕ( t ) k α ) 1+ ρ θ 2 sinh 2 ( ρ ) ρ( θ,0 )= ρ 0 ( θ )

其中 ϕ( t ) 为全局调项

结合曲率上下界: k 0 _ k( ,t ) k 0 ¯ ,方程为一致抛物型,通过Schauder理论,由曲率k C 2 估计可推导 C 2,α 估计,从而得到 C 估计,由曲率及其导数一致有界,得解在所有 t[ 0, ) 上光滑存在。

定理 4.1 (长时存在性):设 γ 0 是双曲空间 Η 2 上的光滑闭凸曲线,对于任何 α>0 ,保面积/长度的流有光滑凸解 γ t 定义在所有时间 t[ 0, ) ,而且 γ t 的曲率估计k

k 0 _ k( ,t ) k 0 ¯

对所有 t[ 0, ) 成立

证明 假设 k α 型流(1)的最大存在时间 T<+ ,即解 γ t 仅在 t[ 0, ) 上存在。

t T 时,由抛物正则性理论

若解在 t=T 处奇异,则某一几何量趋于极端值

但:1) 曲率k始终被 k _ , k ¯ 控制,不会趋于无穷或零;2) 内半径 ρ c 1 >0 ,不会趋于零;3) 径向函数的导数 ρ θ , ρ θθ 一致有界,不会趋于无穷。

矛盾,故假设不成立,流的最大存在时间 T=+ ,即解 γ t 对所有 t[ 0,T ) 存在。

5. 流的收敛性

指数收敛需证明 t+ 时, γ t C 拓扑下以指数速率逼近测地圆,分“子列收敛”与“全体序列指数收敛”两步。

子列收敛到测地圆:利用等周亏格的单调性与等距变换的紧致性

引理 5.1 存在时间序列 t j + 与等距变换 φ t j : Η 2 Η 2 ,使得 φ t j ( γ t j ) C 拓扑下收敛到与 γ 0 等面积/等长度的测地圆。

证明

1) 等周亏格的极限:

等周亏格 Δ( t )=L ( t ) 2 4πA( t )A ( t ) 2 单调递减有下界,故 Δ(t) Δ 0

Δ >0 ,因 0 + dΔ dt dt<+ ,故 dΔ( t ) dt 0

2) 曲率趋于常数:

保面积流: dΔ dt =2L dL dt 0 ,得 dL dt 0

由长度变化率公式及广义Hölder不等式等号成立条件,得 k( x,t ) k ,其中 k 为常数

保长度流: dΔ dt =( 4π+2A ) dA dt 0 ,得 dA dt 0 ,同理得 k( x,t ) k

3) 高阶导数一致有界:

由Schauder正则性理论,对于抛物方程 t k=a( θ,t ) θθ k+b( θ,t )k+c( θ,t ) ,其中 a( θ,t )=F( t ) k α 有界且非0( F( t ) 有界, k 有界且非0), b( θ,t ) c( θ,t ) 有界,故 k 的各阶导数一致有界。同理,支撑函数 u 及曲线 γ 的各阶导数均一致有界。

4) 等距变换与紧致性:

取内球中心 p 0 ( t j ) ,定义 φ t j p 0 ( t j ) 映射到固定点 p

由Arzela-Ascoli定理, φ t j ( γ t j ) 收敛于光滑闭凸曲线 γ ,且 k( γ )= k

双曲平面内常数曲率闭凸曲线必为测地圆,故 γ 为测地圆。

5) 保面积/长度:

若等距变换保持面积,则 γ 的面积 A( γ )= A 0

若等距变换保持长度,则 γ 的长度 L( γ )= L 0

全体序列的指数收敛:通过线性化扰动方程的特征值分析,证明全体序列指数收敛。

证明

步骤 1:扰动分解

γ 的半径为 ρ ,则保面积流 A 0 =2π( cosh ρ 1 ) ,保长度流 L 0 =2πsinh ρ

对充分大t,径向函数 ρ 可分解为:

ρ( θ,t )= ρ +η( θ,t )

其中 η( θ,t ) 为小扰动, η ( t ) C k 0( t+ )

步骤 2:线性化演化方程

将径向函数 ρ 代入标量抛物型方程,忽略高阶小项,在 ρ= ρ 处进行线性化

k =cothρ ,计算k的线性化导数:

d dε k( ρ +εη )| ε=0 = ρ sinh 2 ρ ( η θθ +η ) (14)

进而 k α 的线性化导数为:

d dε k α | ε=0 = α k α1 ρ sinh 2 ρ ( η θθ +η ) (15)

全局项 ϕ( t ) 的线性化:

d dε ϕ| ε=0 = α k α ρ 2π sinh 2 ρ η dθ (16)

将(15)和(16)代入标量抛物型方程,得 η 的线性化方程:

η=λ( η θθ +η 1 2π ηdθ ) (17)

其中 λ= α coth α1 ρ sinh 2 ρ >0

步骤 3:傅里叶系数衰减

η 展开为傅里叶级数 η( θ,t )= n= + a n ( t ) e inθ ,代入线性化方程(17):

t a n ( t )=λ n 2 a n ( t ) ,解得 a n ( t )= a n ( 0 ) e λ n 2 t ,因此:

η( t, ) L 2 n= + | a n ( 0 ) | e λ n 2 t M 0 e λt

其中 M 0 = n= + | a n ( 0 ) |

由初始曲线的光滑性,傅里叶系数衰减,故 M 0 <+

步骤 4: C k 指数收敛

由抛物方程的正则性估计,对于线性抛物方程 t η=λ θθ η ,有:

η( t, ) C k C k ( η( t, ) L 2 + t η( t, ) C k2 )

η( t, ) C k C k e λt

γ( t, ) C 拓扑下指数收敛到测地圆 γ

参考文献

[1] 潘生亮. 关于凸曲线的一些注记及其对曲率流的应用 [J]. 数学年刊 A 辑 (中文版),2000,(01):53-56.
[2] 潘生亮. 一般平面曲线流的一点注记 (英文)[J]. 数学研究,2000,(01):17-26.
[3] 高来源, 郝瑞霞, 潘生亮. 非局部平面曲线流 [J]. 中国科学: 数学,2024,54(03):407-422.
[4] L. Gao, S. Pan, and D. H. Tsai. On an area-preserving inverse curvature flow of convex closed plane curves[J]. Journal of Functional Analysis, vol. 280, no. 8, p. 108931, 2021.
[5] L. Gao, S. Pan, and D. H. Tsai. On a length-preserving inverse curvature flow of convex closed plane curves[J]. Journal of Differential Equation, vol. 269, no. 7, p. 5802-5831, 2020.
[6] Shengliang Pan, Juanna Yang. On a Non-Local Perimeter-Preserving Curve Evolution Problem for Convex Plane Curves[J], Manuscripta Mathematica, 2008, 127(4): 469-484.
[7] H. Guo and Z. Sun. On a family of inverse curvature flows for closed convex plane curves[J]. Nonlinear Analysis Real World Applications, vol. 50, pp. 1-7, 2019.
[8] Ben Andrews, Yong Wei. Quermassintegral Preserving Curvature Flow in Hyperbolic Space[J], Geometric and Functional Analysis, 2018, 28(5): 1183-1208.
[9] K. K. Kwong, Y. Wei, G. Wheeler, and V. M. Wheeler. On an inverse curvature flow in two-dimensional space forms[J]. Mathematische Annalen, vol. 384, no. 1/2, 2022.
[10] D. H. Tsai and X. L. Wang. On area-preserving and length-preserving nonlocal flow of convex closed plane curves[J]. 2014.
[11] H. Guo and S. Zhu. Linear interpolation on k(a)-type area-preserving and length-preserving curve flows[J]. Monat-shefte fur Mathematik, 2023.
[12] Y. Wei and B. Yang. Volume preserving flows for convex curves and surfaces in the hyperbolic space[J]. Journal of Functional Analysis, 2022.