拉普拉斯方程的
C 2 , α 估计又称Schauder估计, 是椭圆型偏微分方程的基本估计
[1 ] ,首先由Schauder发现并证明。Schauder估计已有多种证明方法, 其中比较常见的有两类。一类借助积分型估计,包括以下几种处理方式: 文献 [
1 ] 给出的经典牛顿位势估计方法;文献[
2 ] 引入的Campanato空间并结合局部积分估计给出的证明;文献[
3 ]利用奇异积分估计证明的解的二阶导数对非齐次项的模连续估计。另一类是扰动方法,如文献[
4 -
7 ]中的证明方法以及文献[
8 ]给出的证明方法。扰动法需要构造一列逼近解,通过逼近解的性质来了解目标解的性质。在这方面,文献[
8 ]的处理方式很巧妙,得到了一种量化形式的模连续估计。用
B 1 = B 1 ( 0 ) 表示
R n ( n ≥ 2 ) 中以原点为中心的单位开球,给定
B 1 上的函数
f , 令
ω r = s u p x - y < r f x - f y , x , y ∈ B 1 。
定理A 设u ∈ C 2 B 1 满足方程Δ u x = f x , x ∈ B 1 , 且∫ 0 1 ω r r d r < ∞ 。 则
s u p x , y ∈ B 1 / 2 , x - y ≤ d D 2 u x - D 2 u y ≤ C n ∫ 0 d ω r r d r + d ∫ d 1 ω r r 2 d r + d u L ∞ B 1 。
定理A量化了非齐次项对方程解的二阶导数的影响。特别地, 定理蕴含了经典 Schauder估计。值得注意的是,文献[
8 ]中的方法也可用于处理完全非线性方程、蒙日安培方程等
[9 -10 ] 。
热方程也有Schauder估计
[11 ] 。其证明方法与拉普拉斯方程类似,一种是借助积分型估计,比较典型的是采用各向异性的Campanato空间理论,另一种是扰动方法
[4 -5 ] 。此外,文献 [
8 ]也给出了与定理A类似的模连续估计。我们将其概述如下。给定点
X = x , t 1 , Y = y , t 2 , x , y ∈ R n , t 1 , t 2 ∈ R ,引入抛物距离
定义抛物柱体
令
ω r = s u p δ X , Y < r f X - f Y , X , Y ∈ Q 1 。
定理B 设u ∈ C x , t 2,1 Q 1 满足方程
∂ t u ( x , t ) - Δ u x , t = f x , t , ( x , t ) ∈ Q 1 ,
且∫ 0 1 ω r r d r < ∞ 。 则
s u p x , y ∈ Q 1 / 2 , δ ( X , Y ) ≤ d D 2 u X - D 2 u Y ≤ C n ∫ 0 d ω r r d r + d ∫ d 1 ω r r 2 d r + d u L ∞ Q 1 。
文献[
8 ]的研究是开创性的,也引起了很多后继研究,如文献[
12 ]用同样方法改进了所谓的偏Schauder估计,即仅针对解的部分导数给出量化估计。目前,偏Schauder估计已经成为一个研究热点,相关研究可参见文献[
13 -
14 ],一个比较清晰的综述可参考文献[
15 ],最近的一些研究进展参见文献[
16 ]以及其中的文献。
本文的目的是给出拉普拉斯方程和热方程的
H 1 弱解的模连续估计。本文采用类似文献[
8 ]的处理方法,但在以下两个方面有所不同: 1) 采用齐次方程的解作为逼近解,而不是非齐次项为常数的方程的解;2) 采用能量估计方法,而非极大值原理。
1 拉普拉斯方程的模连续估计
用B r ( x ) 表示R n n ≥ 2 中以x 为中心,r 为半径的开球,在不影响理解的前提下将其简记为B r 。引入记号
ω t = s u p 0 < s ≤ t , B s ⊂ B 1 s 2 1 B s ∫ B s f 2 1 / 2 , f ∈ L 2 B 1 。
显然,ω 是 t 的单调递增函数。
引理1.1 设函数u 是B r 上的调和函数,即Δ u = 0 , 则对任何α 都有
D α u L ∞ ( B r / 2 ) ≤ C n , k r k 1 B r ∫ B r u 2 1 / 2 ,
此处k = α 。
定理1.2 设u ∈ H 1 B 1 是方程
的一个弱解,且
则
s u p x , y ∈ B 1 / 2 , x - y ≤ d u x - u y ≤ C n ∫ 0 d ω t t d t + d ∫ d 1 ω t t 2 d t + d u L 2 B 1 。
特别地,当f ∈ L p ( B 1 ) 时,
则由定理1.2可直接得到如下的经典结果:
推论1.3 设函数u 是方程(1)的一个H 1 弱解,且存在某个n > p > n 2 , f ∈ L p ( B 1 ) 。则对于α = 2 - n p ∈ 0,1 , 有
u C α ( B 1 2 ) ≤ C n , p f L p B 1 + | u | L 2 B 1 。
推论1.3通常由W 2 , p 估计和Sobolev嵌入定理得到。但W 2 , p -估计的证明是比较困难的。此外, 借鉴下面的证明过程,用扰动方法也可直接证明推论1.3,我们将其留给感兴趣的读者。
定理1.2的证明分为三步。首先,通过解边值问题构造逼近函数列,并用能量方法估计逼近函数列的收敛速度。其次,证明逼近函数列的极限与原方程的解相同。最后,结合积分论证给出所需要的估计。
定理1.2的证明 约定在下面的证明过程中常数C 可以不同但仅依赖于维数。首先,任给x ∈ B 1 2 ( 0 ) ,解边值问题Δ u k = 0 , u k ∂ B k = u , 此处B k = B R k x , R k = 1 2 k , k ≥ 1 。 记v k = u - u k 。则Δ v k = f , v k ∂ B k = 0 。根据能量估计,有
根据引理1.1及ω ( t ) 的单调性,有
u k x - u k + 1 x ≤ C 1 B k + 1 ∫ B k + 1 u k - u k + 1 2 1 / 2 ≤ C 1 B k ∫ B k u - u k 2 1 / 2 + C 1 B k + 1 ∫ B k + 1 u - u k + 1 2 1 / 2 = C 1 B k ∫ B k v k 2 1 / 2 + C 1 B k + 1 ∫ B k + 1 v k + 1 2 1 / 2 ≤ C R k 2 1 B k ∫ B k f 2 1 / 2 ≤ C ω R k , k ≥ 1 。
在下面的证明中,我们经常会用到单调函数积分的性质,不再一一说明。我们有
∑ j = k ∞ u j x - u j + 1 x ≤ C ∑ j = k ∞ ω R j ≤ C ∑ j = k ∞ ∫ R j R j - 1 ω t t ≤
C ∫ 0 R k - 1 ω t t < ∞ 。
于是l i m k → ∞ u k ( x ) 存在,将其记为u ^ x 。 则
其次,对于r < 1 / 8 ,存在唯一的k ≥ 3 ,使得 R k + 1 ≤ r < R k 。因此,
1 B r ( x ) ∫ B r ( x ) u - u ^ ( x ) 2 1 / 2 ≤ C 1 B k x ∫ B k x u - u k - 1 x 2 1 / 2 + C | u k - 1 x - u ^ ( x ) | ≤ C 1 B k - 1 x ∫ B k - 1 x u - u k - 1 2 1 / 2 + C 1 B k x ∫ B k x u - u k - 1 x 2 1 / 2 + C ∫ 0 R k - 2 ω t t ≤ C 1 B k x ∫ B k x u k - 1 - u k - 1 x 2 1 / 2 + C ∫ 0 R k - 2 ω t t 。
为了估计u k - 1 - u k - 1 x , 引入
u k - 1 = ∑ j = 2 k - 1 h j + u 1 , h j = u j - u j - 1 。
当z ∈ B k ( x ) 时,由中值定理和引理1.1可得
h j z - h j x ≤ z - x ∇ h j L ∞ B j + 1 ≤ C z - x R j - 1 1 | B j | ∫ B j h j 2 1 / 2 ≤ C z - x R j - 1 1 | B j - 1 | ∫ B j - 1 f 2 1 / 2 ≤ C z - x ω R j - 1 R j - 1 。
于是
u k - 1 z - u k - 1 x ≤ C z - x ∑ j = 2 k - 1 ω R j - 1 R j - 1 + C z - x u L 2 ( B 1 ) ≤ C R k ∫ R k - 2 1 ω t t 2 + C R k u L 2 B 1 。
因此,
1 B k x ∫ B k x u - u k - 1 x 2 1 / 2 ≤ C R k ∫ R k - 2 1 ω t t 2 + u L 2 B 1 。
于是
1 B r ( x ) ∫ B r ( x ) u - u ^ ( x ) 2 1 / 2 ≤ C ∫ 0 r ω t t + r ∫ r 1 ω t t 2 + r u L 2 B 1 。
由此,结合式(2) 可知l i m r → 0 1 | B r | ∫ B r u = u ^ 。 由勒贝格微分定理有u = u ^ a . e . , x ∈ B 1 / 2 0 。
最后,任取
x , y ∈ B 1 / 2 ( 0 ) 。当
x - y = r < 1 / 8 时,由标准的积分论证(参见文献[
17 ]中定理3.1的证明)得
u x - u y ≤ C 1 B r x ∫ B r x u - u x 2 1 / 2 + C 1 B r ( y ) ∫ B r ( y ) u - u ( y 2 1 / 2 ≤ C ∫ 0 r ω t t + r ∫ r 1 ω t t 2 + r u L 2 B 1 。
当x - y = r ≥ 1 / 8 时,在x , y 中间插入至多7 个点总可以使得相邻两个点的距离小于1 / 8 ,从而使得证明转化为已经得证的情形。证毕。
借助标准奇延拓方法,上述证明还可以处理边界估计,此时要求解在边界为零。
定理1.4 设函数u 是方程
的一个H 1 弱解,且∫ 0 1 ω t t d t < ∞ 。 则
s u p x , y ∈ B 1 / 2 + , x - y ≤ d u x - u y ≤ C n ∫ 0 d ω t t d t + d ∫ d 1 ω t t 2 d t + d u L 2 B 1 + 。
此处B 1 + = B 1 0 ∩ { x | x n > 0 } 。
2 热方程的模连续估计
如前所述,用Q s ( X ) 表示以X 为中心的抛物柱体,也简记为Q s . 约定Q 1 = Q 1 ( 0 ) 并定义
ω r = s u p 0 < s ≤ r , Q s ⊂ Q 1 s 2 1 Q s ∫ Q s f 2 1 / 2 , f ∈ L 2 Q 1 。
显然,ω 是关于r 的单调递增函数。
引理2.1 [18 ] 设函数
u 是
Q r 上齐次热方程的解,即
∂ t u - Δ u = 0 . 则对任何的
α , m 都有
D x α D t m u L ∞ ( Q r 2 ) ≤ C n , k , m r k + 2 m 1 Q r ∫ Q r u 2 1 / 2 ,
此处k = α 。
定理2.2 设u ∈ H 1 Q 1 是热方程
∂ t u ( x , t ) - Δ u x , t = f x , t , ( x , t ) ∈ Q 1 (3)
的一个弱解,且
则
s u p X , Y ∈ Q 1 2 , δ ( X , Y ) ≤ d u X - u Y ≤ C n ∫ 0 d ω r r d r + d ∫ d 1 ω r r 2 d r + d u L 2 Q 1 。
特别地,当f ∈ L p ( Q 1 ) 时,
此时便得到如下的经典结果:
推论2.3 设函数u 是方程(3)的一个H 1 弱解,且存在某个n + 2 > p > n + 2 2 , f ∈ L p B 1 。 则对α = 2 - n + 2 p ∈ 0,1 , 有
u C α , α / 2 ( Q 1 2 ) ≤ C n , p f L p Q 1 + u L 2 Q 1 。
定理2.2的证明 鉴于证明过程与定理1.2的证明类似,简便起见我们仅给出大致过程。首先,任给X ∈ Q 1 / 2 ( 0 ) ,解边值问题Δ u k = 0 , u k ∂ T Q k = u 。 此处∂ T Q k 表示Q k 的抛物边界,
记v k = u - u k 。则∂ t v k - Δ v k = f , v k ∂ T Q k = 0 。 根据能量估计,有
根据引理2.1,有
u k X - u k + 1 X ≤ C 1 Q k + 1 ∫ Q k + 1 u k - u k + 1 2 1 / 2 ≤ C ω R k , k ≥ 1 。
因此,
∑ j = k ∞ u j X - u j + 1 X ≤ C ∑ j = k ∞ ω R j ≤ C ∑ j = k ∞ ∫ R j R j - 1 ω r r ≤
C ∫ 0 R k - 1 ω r r < ∞ 。
于是l i m k → ∞ u k ( X ) 存在,记其为u ^ ( X ) 。从而有
其次,对于r < 1 / 8 ,存在唯一的k ≥ 3 ,使得R k + 1 ≤ r < R k 。因此,
1 Q r ( X ) ∫ Q r ( X ) u - u ^ ( X ) 2 1 / 2 ≤ C 1 Q k X ∫ Q k X u k - 1 - u k - 1 x 2 1 / 2 + C ∫ 0 R k - 2 ω s s 。
为了估计u k - 1 - u k - 1 x , 引入
u k - 1 = ∑ j = 2 k - 1 h j + u 1 , h j = u j - u j - 1 。
当Z ∈ Q k ( X ) 时,由中值定理和引理2.1得
于是
u k - 1 Z - u k - 1 X ≤ C R k ∫ R k - 2 1 ω s s 2 + C R k | u | L 2 Q 1 。
因此,
1 Q k X ∫ Q k X u - u k - 1 X 2 1 / 2 ≤ C R k ∫ R k - 2 1 ω s s 2 + u L 2 Q 1 。
于是
1 Q r ( X ) ∫ Q r ( X ) u - u ^ ( X ) 2 1 / 2 ≤ C ∫ 0 r ω s s + r ∫ r 1 ω s s 2 + r u L 2 Q 1 。
由此,结合条件(2)可知l i m r → 0 1 | Q r | ∫ Q r u = u ^ 。 由勒贝格微分定理,有
最后,任取X , Y ∈ B 1 / 2 ( 0 ) 。当δ ( X , Y ) = r < 1 / 8 时,由标准的积分论证得到
u X - u Y ≤ C 1 Q r X ∫ B r X u - u X 2 1 / 2 + C 1 Q r ( Y ) ∫ Q r ( Y ) u - u ( Y ) 2 1 / 2 ≤ C ∫ 0 r ω s s + r ∫ r 1 ω s s 2 + r u L 2 Q 1 。
当δ ( X , Y ) = r ≥ 1 / 8 时,在X , Y 中间插入至多7 个点,总可以使得相邻两个点的距离小于1 8 ,从而将证明转化为已证明情形。证毕。
国家重点研发计划项目(2024YFA1012803)
国家自然科学基金(12271400)
天津大学极基础项目(2025XJ21-0010)