一类抛物型界面问题的正则性分析
Regular Analysis of a Class of Parabolic Interface Problems
DOI: 10.12677/AAM.2023.122068, PDF, HTML, XML, 下载: 195  浏览: 2,528 
作者: 史雪婷, 杨田洁:中国工程物理研究院,北京;北京应用物理与计算数学研究所,北京
关键词: 对流扩散方程非完美界面条件Henry界面条件De Giorgi迭代H?lder连续性Convective Diffusion Equation Imperfect Interface Condition Henry Interface Condition De Giorgi Iteration H?lder Continuity
摘要: 界面问题用于各种工程应用和物理、化学、生物现象的建模,特别是涉及具有不同扩散性、密度、渗透性或导电性的多种不同材料的现象,其在界面上由一定条件耦合。本文考虑具有溶解质输运的线性两相流模型,分别在相界面处耦合非完美界面条件和Henry界面条件。由于解在界面上的跳跃使得解在各自材料区域上具有比在整个区域上更高的正则性,针对这类界面问题的正则性分析,本文给出一个完整的泛函分析过程,采用De Giorgi迭代方法证明该模型弱解的相关性质,进而证明弱解及其梯度的Hölder连续性。此外,对于Henry界面问题,本文给出了梯度的Lq估计(存在q > 2)。
Abstract: Interface problems are used for various engineering applications and modeling of physical, chemi-cal, and biological phenomena, especially those involving a number of different materials with dif-ferent diffusion, density, permeability or conductivity, which are coupled by certain conditions at the interface. In this paper, a linear two-phase flow model with solute transport is considered, cou-pling imperfect interface condition and Henry interface condition at the phase interface, respec-tively. Since the jump of the solution at the interface makes the solution have higher regularity in the respective material region than on the whole region, for the regularity analysis of such interface problems, this paper presents a complete functional analysis process, and uses the De Giorgi itera-tion method to prove the correlation properties of the weak solution of the model, and then prove the Hölder continuity of the weak solution and its gradient. In addition, for the Henry interface problem, this paper gives the Lq estimation of the gradient (q > 2 is present).
文章引用:史雪婷, 杨田洁. 一类抛物型界面问题的正则性分析[J]. 应用数学进展, 2023, 12(2): 645-678. https://doi.org/10.12677/AAM.2023.122068

1. 引言

当偏微分方程中的系数、或解、或法向流在界面处不连续时,就会产生界面问题。界面问题在自然界以及科学与工程中都有应用。在材料科学中,文献 [1] 提出了一种聚合物基体复合材料的耦合扩散行为模型,在由不同成分构成的复合材料的接触面上产生界面问题。在生物技术中被脂膜所包围的流体所形成的生物膜模型在界面处的相互作用 [2] 。文献 [3] 研究了两相流体动力学界面模型的斯托克斯问题,在界面处具有不连续的密度和粘度系数以及压力溶液。文献 [4] 考虑二维静止热传导椭圆界面问题,其传导系数在光滑的内部界面上是不连续的。文献 [5] 提出一个基于区域分解理论的非完美界面问题的保极值迭代方法,从而使得界面条件自然地嵌入到子域的边界条件中。针对非完美界面问题,文献 [6] 提出了一种保持离散极值原理(DMP)和守恒性的有限体积格式。此外,界面问题也被应用于合金凝固、晶体生长以及在生物系统中 [7] 等。

考虑具有溶解质输运的线性两相流问题(见图1)。设 ω ( R 4 ) = μ ( R 4 ) μ ˜ ( R 4 ) μ min { ess inf Q R 4 , 1 u 1 , k ˜ ( ) , ess inf Q R 4 , 2 u 2 , k ˜ ( ) } ( 1 1 2 s ) ω ( R ) + C ( R α 1 ( [ u 0 1 ] C α 1 ( Ω ¯ 1 ) + [ u 0 2 ] C α 1 ( Ω ¯ 2 ) ) + ( M + F 0 ) R 1 n + 2 p ) . 是一个包含两相不可混溶,不可压缩的流动系统(液–液或液–气)的有界区域,具有光滑边界 ω ( R 4 ) = μ ( R 4 ) μ ˜ ( R 4 ) max { ess sup Q R 4 , 1 u 1 , k ( + ) , ess sup Q R 4 , 2 u 2 , k ( + ) } μ ˜ ( 1 1 2 s ) ω ( R ) + C ( R α 1 ( [ u 0 1 ] C α 1 ( Ω ¯ 1 ) + [ u 0 2 ] C α 1 ( Ω ¯ 2 ) ) + ( M + F 0 ) R 1 n + 2 p ) . 0 < R R 0 1 是一个开子区域且有边界 ω ( R 4 ) ( 1 1 2 s ) ω ( R ) + C R α ( [ u 0 1 ] C α 1 ( Ω ¯ 1 ) + [ u 0 2 ] C α 1 ( Ω ¯ 2 ) + M + F 0 ) , ,故有 0 < α min { α 1 , 1 n + 2 p } 。每个子域中分别包含一个相,这些相通过界面 ω ( R ) C ( R R 0 ) α ( ω ( R 0 ) + R 0 α ( [ u 0 1 ] C α 1 ( Ω ¯ 1 ) + [ u 0 2 ] C α 1 ( Ω ¯ 2 ) + M + F 0 ) ) , R ( 0 , R 0 ] , 分隔开。两相中都含有一种溶解的物质,这种物质由于对流和分子扩散而被输运,并不粘附在界面上。在本文中,假设所考虑的溶解质输运的两相流模型是理想模型,在每个相中都会发生对流传质和扩散传质。假设不会发生相变,反应;界面处没有传质阻力;也不会因为传质而引起界面湍动等。

Figure 1. Two-phase flow model

图1. 两相流模型

在相界面处,考虑定常界面情形。同时对于界面处考虑非完美界面条件 [5] [6] [8] 和Henry界面条件。如果给定的数据,界面 u W ( 0 , T ; V ) 和外边界 X 0 = ( x 0 , 0 ) Γ 光滑,则问题的解在各个区域也非常光滑,但由于界面处的跳跃会使得解的全局正则性降低。在文献 [9] 中,只给出了Henry界面问题弱解的适定性。而在本文中我们着重讨论在上述两种界面条件下的线性两相流模型弱解的性质,我们的主要结论是在给定的Sobolev空间中,利用De Giorgi迭代法来估计线性问题的弱解,得到弱解在界面附近的局部性质,在此基础之上可进一步得到线性模型的弱解及其梯度的Hölder连续性。

设溶解质的浓度为 Q R ( X 0 ) = B R ( x 0 ) × ( 0 , R 2 ] ,标量 0 < R 1 ,这个问题可以用浓度为 d = min { 1 , d i s t { x 0 , Ω } } 的对流–扩散方程来建模。 0 < R d 为界面 0 < α min { α 1 , 1 n + 2 p } 上的单位外法向量,由 C 1 指向 max { osc Q R , 1 u 1 , osc Q R , 2 u 2 } C R α ( M + F 0 + [ u 0 1 ] C α 1 ( Ω ¯ 1 ) + [ u 0 2 ] C α 1 ( Ω ¯ 2 ) ) , 。对于 α , C n , λ , Λ , p 为流体速度场, Q R , i = Q R Q i ( i = 1 , 2 ) 为扩散系数矩阵, F 0 = i = 1 2 f i L 2 ( 0 , T ; L n p n + p ( Ω i ) ) < 为标量传输系数。

非完美界面条件是指解在界面上的跳跃与界面两侧连续的法向流成正比。可得非完美界面问题的数学模型公式如下:

Q T (1.1)

其中 Γ × ( 0 , T ] Q 1 Q 2 满足非完美界面条件。

Henry界面条件要求界面在瞬间平衡的情况下,界面两侧的溶质浓度呈恒定比。同时施加另一个界面条件,即要求法向流在界面处是连续的。故有Henry界面问题的数学模型公式如下:

u (1.2)

其中 u 是分片常数,即在 Γ × ( 0 , T ] ( Γ × [ 0 , T ] ) Ω ,一般有 u C α ( Ω ¯ 1 × [ 0 , T ] ) × C α ( Ω ¯ 2 × [ 0 , T ] ) . 0 < α min { α 1 , 1 n + 2 p } p Q T C 1 + α 满足Henry界面条件。

2. 非完美界面模型

在本节中,讨论非完美界面问题(1.1)。

2.1. 函数分析框架

对于 Γ × ( 0 , T ] C 1 + α ,一般Sobolev空间记为 k p g i C α ( Ω ¯ i × [ 0 , T ] ) 。特别的, w i 。首先引进一些合适的空间:

w i j M 1 , n + 1 + α ( Ω i ; d * ) ( j = 1 , , n , i = 1 , 2 ) Λ 1

[ k p g i ] α + w i j 1 , n + 1 + α Λ 1 . u W ( 0 , T ; V ) ,且

f i M 1 , n + 1 + α ( Q i ; δ ˜ )

由文献 [10] 的定理3.13 (p. 175)可得

u 0 i C 1 + α ( Ω ¯ i )

其中 u C 1 + α ( Ω ¯ 1 × [ 0 , T ] ) × C 1 + α ( Ω ¯ 2 × [ 0 , T ] ) , 表示嵌入。

迹算子:

i = 1 2 u i 1 + α C ( n , λ , Λ , Λ 1 , m ¯ , α , Q T ) ( i = 1 2 u 0 i 1 + α + i = 1 2 u i α + i = 1 2 f i 1 , n + 1 + α ) .

是有界的。令 u ˜ = β u ,并且 [ u ˜ ] Γ = 0

下面给出问题(1.1)中的系数的有关假设。

假设2.1.1区域 u ˜ u ˜ | Ω i × ( 0 , T ] = u ˜ i 是有界区域, 1 β u ˜ t + 1 β w u ˜ d i v ( 1 β K ( x , t ) u ˜ ) = f ( x , t ) , x Ω 1 Ω 2 , t ( 0 , T ] ; 1 β 1 K 1 ( x , t ) u ˜ 1 n = 1 β 2 K 2 ( x , t ) u ˜ 2 n , x Γ , t ( 0 , T ] ; [ u ˜ ] Γ = 0 , x Γ , t ( 0 , T ] ; u ˜ ( , 0 ) = u ˜ 0 , x Ω 1 Ω 2 ; u ˜ ( , t ) = 0 , x Ω , t ( 0 , T ] . ,界面 u ˜ 0 α 1 ( 0 , 1 ) Ω 上的 Γ 维Hausdroff测度。

假设2.1.2

1) 假设 C 1 + α 1 H = L 2 ( Ω ) 上是一个充分光滑的速度场,并且满足:

i) ( u , v ) H = Ω 1 β u v d x , H

ii) 由于流体是不可压缩的,故有在 V = H 0 1 ( Ω ) u V 2 = u H 2 + u H 2

iii) 在界面 u V 上速度场满足: u H

2) 扩散系数矩阵 V H H V 和标量传输系数 f L 2 ( 0 , T ; L 2 ( Ω ) ) 如下假设:

i) 扩散系数矩阵 u ˜ 0 ( x ) H F L 2 ( 0 , T ; V ) 上是可测的,一致有界的和一致椭圆的,并且 u ˜ W ( 0 , T ; V ) 是对称的,即存在常数 { d u ˜ d t ( v ˜ ) + a ( t ; u ˜ ( t ) , v ˜ ) ( ) F ( t ) ( v ˜ ) , v ˜ V ; u ˜ ( 0 ) = u ˜ 0 ( x ) . 使得

u ˜

并且对于任意 u ˜ ,在 d u ˜ d t ( v ˜ ) = i = 1 2 Ω i 1 β i u ˜ i t v ˜ i d x , 中几乎处处成立

a ( t ; u ˜ ( t ) , v ˜ ) = i = 1 2 Ω i 1 β i w i ( u ˜ i ) v ˜ i + 1 β i ( v ˜ i ) K i u ˜ i d x ,

ii) 传输系数 F ( t ) ( v ˜ ) = i = 1 2 Ω i f i v ˜ i d x . a ( t ; u ˜ ( t ) , v ˜ ) 上是可测的,并且存在常数 V × V 使得

Q T

传输系数矩阵 u ˜ ,则M是半正定矩阵,并且对于任意 u ˜ W ( 0 , T ; V ) ,有

X 0 = ( x 0 , t 0 ) Γ × ( 0 , T ]

在上述假设的基础之上,给出问题 的弱形式。给定 Q ρ , τ ( X 0 ) = B ρ ( x 0 ) × ( t 0 , t 0 + τ ] Q T 。对于 k ,积分可得:

ξ ( x , t ) C ( [ t 0 , t 0 + τ ] ; C 0 ( B ρ ( x 0 ) ) )

0 ξ 1

ξ ( , t 0 ) = 0

sup t 0 < t t 0 + τ i = 1 2 ξ ( u ˜ i k ) ± ( , t ) L 2 ( B ρ , i ) 2 + λ 1 i = 1 2 ( ξ ( u ˜ i k ) ± ) L 2 ( t 0 , t 0 + τ ; L 2 ( B ρ , i ) ) 2 C * [ ( ξ t L ( Q ρ , τ ) + ξ L ( Q ρ , τ ) 2 ) i = 1 2 ( u ˜ i k ) ± L 2 ( t 0 , t 0 + τ ; L 2 ( B ρ , i ) ) 2 + F 0 , ρ , τ 2 i = 1 2 | Q ρ , τ , i [ ( u ˜ i k ) ± > 0 ] | 1 2 p ] ,

所以(1.1)的变分问题为:设 B ρ , i = B ρ Ω i , Q ρ , τ , i = Q ρ , τ Q i 。寻找弱下解(弱上解) 0 < ρ , τ < 1 使得

p > n + 2 (2.1)

如果 λ 1 > 0 既是弱下解,又是弱上解,则称 F 0 , ρ , τ = i = 1 2 f i L 2 ( t 0 , t 0 + τ ; L n p n + p ( B ρ , i ) ) > 0 是弱解。

由文献 [10] (p. 184定理3.16)可知,若要证明线性问题(2.1)弱解的存在唯一性,只需表明 C * n , Λ , p , c 0 , β 1 , β 2 上是连续的和强制的。事实上,根据假设2.1.2和迹算子的有界性可得 D G ( Q T ) = D G ( Q T ; λ 1 , p , n , F 0 , ρ , τ , C * ) 是连续的。对于 u ˜ W ( 0 , T ; V ) 的强制性:任意 u ˜ D G + ( Q T ) ,有

u ˜ W ( 0 , T ; V )

其中

u ˜ D G ( Q T )

故有

D G ( Q T ) = D G + ( Q T ) D G ( Q T )

由于

u ˜ W ( 0 , T ; V )

p > n + 2

f i L 2 ( 0 , T ; L n p n + p ( Ω i ) ) ( i = 1 , 2 )

所以

u ˜ D G + ( Q T )

其中 u ˜ W ( 0 , T ; V ) 是常数。

2.2. 预备知识

以下定义和定理参考文献 [11] 和 [12] :对于任意 u ˜ D G ( Q T ) ,令抛物距离为

C *

假设D是 n , Λ , p , c 0 中的有界区域,对于任意 F 0 , ρ , τ = i = 1 2 f i L 2 ( t 0 , t 0 + τ ; L n p n + p ( B ρ , i ) ) < ,令 u ˜ D G ( Q T ) ,其中 sup t 0 < t t 0 + τ i = 1 2 ξ ( u ˜ i k ) ± ( , t ) L 2 ( B ρ , i ) 2 + λ 2 i = 1 2 ξ ( u ˜ i k ) ± L 2 ( t 0 , t 0 + τ ; L 2 * ( B ρ , i ) ) 2 C * [ ( ξ t L ( Q ρ , τ ) + ξ L ( Q ρ , τ ) 2 ) i = 1 2 ( u ˜ i k ) ± L 2 ( t 0 , t 0 + τ ; L 2 ( B ρ , i ) ) 2 + F 0 , ρ , τ 2 i = 1 2 | Q ρ , τ , i [ ( u ˜ i k ) ± > 0 ] | 1 2 p ] , 0 < α α 1 是D关于抛物距离的直径。

定义2.2.1 (Morrey空间)对于 u ˜ C α ( Ω ¯ × [ 0 , T ] ) . ,令 Γ C 2 表示由 x 0 Γ 中满足

ρ > 0

的所有函数u所组成的赋范线性空间。特别地,若 C 2 是有界区域,对于任意 γ : n 1 ,记 B 1 = Ω 1 B ρ ( x 0 ) = { x B ρ ( x 0 ) | x n > γ ( x 1 , , x n 1 ) } , Π = Γ B ρ ( x 0 ) = { x B ρ ( x 0 ) | x n = γ ( x 1 , , x n 1 ) } , B 2 = Ω 2 B ρ ( x 0 ) = { x B ρ ( x 0 ) | x n < γ ( x 1 , , x n 1 ) } , ,定义 B ρ ( x 0 ) 是由 x 0 中满足

{ y j = x j = : Φ j ( x ) , ( j = 1 , , n 1 ) ; y n = x n γ ( x 1 , , x n 1 ) = : Φ n ( x ) .

的所有函数u所组成的赋范线性空间。

定义2.2.2 (Campanato空间)对于 y = Φ ( x ) . ,以 { x j = y j = : Ψ j ( y ) , ( j = 1 , , n 1 ) ; x n = y n + γ ( y 1 , , y n 1 ) = : Ψ n ( y ) . 表示由 x = Ψ ( y ) . 中满足

Ψ = Φ 1

的所有函数u所组成的赋范线性空间,其上的范数定义为

Φ : x y

Π 表示u在 J = Φ 上的积分平均值,即 J 1 = Ψ

定义2.2.3 (Hölder空间)对于 | det ( J ) | = | det ( J 1 ) | = 1 ,以 Φ ( x 0 ) = y 0 Φ ( Π ) { y n = 0 } 表示满足

R ( 0 , 1 ]

的所有函数u所组成的线性空间,其上的范数定义为

B R ( y 0 ) Φ ( B ρ ( x 0 ) )

显然有

B R , 1 = B R ( y 0 ) { y n > 0 } Φ ( B 1 ) ,

Σ = B R ( y 0 ) { y n = 0 } Φ ( Π ) , 表示由 B R , 2 = B R ( y 0 ) { y n < 0 } Φ ( B 2 ) . 中满足

Q R ( y 0 , t 0 ) = B R ( y 0 ) × ( t 0 R 2 , t 0 ]

所有函数u所组成的赋范线性空间,其上的范数定义为

Q R , i = B R , i × ( t 0 R 2 , t 0 ]

定义2.2.4 称D是(A)型区域,如果存在常数A,使得对于任意 0 < t 0 R 2 < t 0 < T ,都有

v ( y , t ) = u ˜ ( Ψ ( y ) , t ) ,

定理2.2.1 设D是(A)型区域, v | B R , i = v i ( i = 1 , 2 ) ,则当 w ˜ ( y ) = w ( Ψ ( y ) ) 时,

K ˜ ( y , t ) = K ( Ψ ( y ) , t )

其中 f ˜ ( y , t ) = f ( Ψ ( y ) , t ) v ˜ 0 ( y ) = u ˜ 0 ( Ψ ( y ) ) 表示 n ˜ B R , 1 同时成立。

在以后的小节中,记

B R , 2

{ u ˜ t = v t ; x u ˜ = J y v ; d i v x ( K ( x , t ) x u ˜ ) = d i v y ( J K ˜ ( y , t ) J y v ) ; n = J n ˜ | J n ˜ | .

z = J w ˜

A = J K ˜ ( y , t ) J

2.3. 弱解的极值原理

引理2.3.1 (文献( [13] , p. 95)引理5.6)非负序列 1 β v t + 1 β z y v d i v y ( 1 β A y v ) = f ˜ ( y , t ) , y B R , 1 B R , 2 , t [ t 0 R 2 , t 0 ] ; 1 β 1 A 1 y v 1 n ˜ = 1 β 2 A 2 y v 2 n ˜ , y Σ , t [ t 0 R 2 , t 0 ] ; [ v ] Σ = 0 , y Σ , t [ t 0 R 2 , t 0 ] . 满足递推关系式 z , A ,其中 c 0 , λ , Λ r ( 0 , R ] ,则如果 ξ 2 r ( y ) C 0 ( B 2 r ( y 0 ) ) ,必有 τ 2 r ( t ) C ( )

推论2.3.1令 0 ξ 2 r 1 是定义在 B r ( y 0 ) 上的非增非负函数,并且存在 ξ 2 r = 1 使得对于任意 | ξ 2 r | 2 C r 2 ,有 ( t 0 r 2 , + ) ,则当 τ 2 r = 1 ,有 ( , t 0 4 r 2 ]

定理2.3.1 (极值原理)令问题 的系数满足假设2.1.2。如果 τ 2 r = 0 是问题的弱下解,且对于某个常数 0 τ 2 r 1 | d τ 2 r d t | C r 2 ,则

v ( y , t ) L 1 ( Q R )

其中C仅依赖于 v ˜ 2 r ( t ) = ( B 2 r ξ 2 r 2 ( y ) d y ) 1 B 2 r v ( y , t ) ξ 2 r 2 ( y ) d y , t ( t 0 4 r 2 , t 0 ] , B 2 r ( v ( y , t ) v ˜ 2 r ( t ) ) ξ 2 r 2 ( y ) d y = 0. ,且 v ( y , t ) = u ˜ ( x , t ) W ( 0 , T ; V )

证明 令 v ˜ 2 r ( t ) ,for d v ˜ 2 r ( t ) d t L 2 ( ( t 0 4 r 2 , t 0 ] ) ,取测试函数 p > n + 2 ,则 q 1 > n p n + p ,且 f ˜ L q 1 ( Q R ) 。由于 v ( y , t ) 是弱下解,有

sup t 0 r 2 < t t 0 i = 1 2 B r , i | v i v ˜ r ( t ) | 2 d y + λ 0 i = 1 2 Q r , i | v i | 2 d y d t C { 1 r 2 i = 1 2 Q 2 r , i | v i v ˜ 2 r ( t ) | 2 d y d t + r ( n + 2 ) ( 1 2 p ) 4 n i = 1 2 ( Q 2 r , i | f ˜ i | n p n + p d y d t ) 2 ( n + p ) n p } ,

即,

n , β 1 , β 2 , c 0 , Λ , p

由假设2.1.2可得

λ 0

类似强制性的证明可得

λ

其中 Q 2 r Q R

φ ( y , t ) = ( v v ˜ 2 r ( t ) ) ξ 2 r 2 ( y ) τ 2 r 2 ( t ) ,故有

ξ , τ

i = 1 2 B 2 r , i 1 β i ( ( v i v ˜ 2 r ( t ) ) ξ τ ) t ( v i v ˜ 2 r ( t ) ) ξ τ + 1 β i A i ( ξ τ v i ) ( ξ τ v i ) d y = i = 1 2 B 2 r , i f ˜ i ( v i v ˜ 2 r ( t ) ) ξ 2 τ 2 d y + i = 1 2 B 2 r , i 1 β i ( v i v ˜ 2 r ( t ) ) 2 ξ 2 τ d τ d t d y i = 1 2 B 2 r , i 1 β i v ˜ 2 r ( t ) t ( v i v ˜ 2 r ( t ) ) ξ 2 τ 2 d y + i = 1 2 B 2 r , i 1 β i z i ( v i v ˜ 2 r ( t ) ) τ ξ ( ( v i v ˜ 2 r ( t ) ) ξ τ ) d y i = 1 2 B 2 r , i 1 β i A i v i τ 2 2 ξ ( ξ ) ( v i v ˜ 2 r ( t ) ) d y . 充分小,可得

| i = 1 2 B 2 r , i 1 β i v ˜ 2 r ( t ) t ( v i v ˜ 2 r ( t ) ) ξ 2 τ 2 d y | τ 2 | v ˜ 2 r ( t ) t | max { 1 β 1 , 1 β 2 } | B 2 r ( v v ˜ 2 r ( t ) ) ξ 2 d y | = 0 ,

上式两边对t从0到T积分

| i = 1 2 B 2 r , i 1 β i z i ( v i v ˜ 2 r ( t ) ) τ ξ ( ( v i v ˜ 2 r ( t ) ) ξ τ ) d y | 2 ε i = 1 2 1 β i B 2 r , i | ξ τ v i | 2 d y + C i = 1 2 B 2 r , i 1 β i | ξ | 2 | v i v ˜ 2 r ( t ) | 2 d y ,

因此

| i = 1 2 B 2 r , i 1 β i A i v i τ 2 ( 2 ξ ) ( ξ ) ( v i v ˜ 2 r ( t ) ) d y | ε i = 1 2 1 β i B 2 r , i | ξ τ v i | 2 d y + C i = 1 2 B 2 r , i 1 β i | ξ | 2 | v i v ˜ 2 r ( t ) | 2 d y ,

对于 i = 1 2 B 2 r , i f ˜ i ( v i v ˜ 2 r ( t ) ) ξ 2 τ 2 d y 2 ε i = 1 2 B 2 r , i | ξ τ v i | 2 d y + 2 ε i = 1 2 B 2 r , i | ξ | 2 | v i v ˜ 2 r ( t ) | 2 d y + C r n ( 1 2 p ) i = 1 2 ( B 2 r , i | f ˜ i | n p n + p d y ) 2 ( n + p ) n p ,

ε

i = 1 2 B 2 r , i ( ( v i v ˜ 2 r ( t ) ) ξ τ ) t ( v i v ˜ 2 r ( t ) ) ξ τ d y + ( λ 5 ε ) i = 1 2 B 2 r , i | ξ τ v i | 2 d y C { i = 1 2 B 2 r , i ( | d τ d t | + | ξ | 2 ) | v i v ˜ 2 r ( t ) | 2 d y + r n ( 1 2 p ) i = 1 2 ( B 2 r , i | f ˜ i | n p n + p d y ) 2 ( n + p ) n p } .

t ( t 0 4 r 2 , t 0 ] ,故有

sup t 0 4 r 2 < t t 0 i = 1 2 B 2 r , i | ( v i v ˜ 2 r ( t ) ) ξ τ | 2 d y + λ 0 i = 1 2 Q 2 r , i | ξ τ v i | 2 d y d t C { 1 r 2 i = 1 2 Q 2 r , i | v i v ˜ 2 r ( t ) | 2 d y d t + r ( n + 2 ) ( 1 2 p ) 4 n i = 1 2 ( Q 2 r , i | f ˜ i | n p n + p d y d t ) 2 ( n + p ) n p } .

由推论2.3.1可得,对于 ξ r ,有 B r ,即,

B r ξ r 2 d y | B r 2 | = | B r | 2 n

所以在 g ( l ) = B r ( v ( y , t ) l ) 2 ξ r 2 d y 中, g ( l ) ,因此

l = ( B r ξ r 2 ( y ) d y ) 1 B r v ( y , t ) ξ r 2 ( y ) d y = v ˜ r ( t )

推论2.3.2令问题(1.1)的系数满足假设2.1.2。如果 B r | v ˜ 2 r ( t ) v ˜ r ( t ) | 2 ξ r 2 d y 4 B r | v v ˜ 2 r ( t ) | 2 ξ r 2 d y , 是问题的弱解,且对于某个常数 B r | v v ˜ r ( t ) | 2 d y 2 B r | v v ˜ 2 r ( t ) | 2 d y + 2 B r | v ˜ 2 r ( t ) v ˜ r ( t ) | 2 d y 2 ( 1 + 2 n + 2 ) B 2 r | v v ˜ 2 r ( t ) | 2 ξ 2 r 2 d y , w ( y ) W 1 , q ( B r ) ( 1 q < )

w W 1 , q ( B r ) C { w L q ( B r ) + ( B r ξ r 2 ( y ) d y ) 1 B r w ξ r 2 d y } ,

其中C仅依赖于 1 q < n ,且 w L s ( B r ) C r 1 + n s n q w W 1 , q ( B r ) , ( 1 s n q n q ) ,

2.4. 弱解的局部性质

定义2.4.1 称定义于 q = n 上的函数 w L s ( B r ) C r 1 + n s n q w W 1 , q ( B r ) , ( 1 s < ) . 属于De Giorgi类,如果 w w ˜ r L s ( B r ) C r 1 + n s n q w L q ( B r ) , 1 q < n ,且对于 1 s n q n q q = n 1 s < 满足 p > n + 2 ,并且 q 1 > n p n + p ,有下式成立:

f ˜ L q 1 ( Q R ) (2.2)

其中 v ( y , t ) sup t 0 r 2 < t t 0 i = 1 2 B r , i | v i v ˜ r ( t ) | 2 d y + λ 0 i = 1 2 Q r , i | v i | 2 d y d t C { i = 1 2 Q 2 r , i | v i | 2 d y d t + r ( n + 2 ) ( 1 2 p ) 4 n i = 1 2 ( Q 2 r , i | f ˜ i | n p n + p d y d t ) 2 ( n + p ) n p } , ,常数 n , β 1 , β 2 , c 0 , Λ , p λ 0 只与 λ 有关, Q ¯ 2 r Q R Q n + 1 F L l o c m ( Q ) 依赖于 G L l o c m 1 ( Q ) ( 1 < m < m 1 ) ,记De Giorgi类为 F , G 0 。如果 Q r Q ¯ 2 r Q ,且满足(2.2)+,则记 1 | Q r | Q r F m d y d t a { 1 | Q 2 r | Q 2 r G m d y d t + ( 1 | Q 2 r | Q 2 r F d y d t ) m } + θ 1 | Q 2 r | Q 2 r F m d y d t , ;如果 a 1 ,且满足(2.2)-,则记 θ [ 0 , 1 ) 。显然 ε = ε ( a , θ , m , n ) > 0

定理2.4.1 设问题(1.1)的系数满足假设2.1.2。如果 F L l o c m 0 ( Q ) ( m < m 0 < min { m + ε , m 1 } ) , 是问题的弱下解,且对于某个常数 1 | Q r | Q r F m 0 d y d t c { 1 | Q 2 r | Q 2 r G m 0 d y d t + ( 1 | Q 2 r | Q 2 r F m d y d t ) m 0 m } , n , m , a , θ ,则 a ;如果 ε 0 是问题的弱上解,则 Q 4 r 。其中 Q 2 r 依赖于 L q ,并且 p > n + 2

证明 如果 q 1 > n p n + p 是问题的弱下解,令 f ˜ L q 1 ( Q R ) q > 2 上的截断函数, v i L l o c q ( Q R , i ) , ( i = 1 , 2 ) ,并且 Q 4 r Q ¯ 4 r Q R 。记 i = 1 2 Q r , i | v i | q d y d t C { i = 1 2 Q 4 r , i | f ˜ i | q 1 d y d t + r ( n + 2 ) ( 1 q 2 ) ( i = 1 2 Q 4 r , i | v i | 2 d y d t ) q 2 } , 。取测试函数为 i = 1 2 Q 2 r , i | v i v ˜ 2 r ( t ) | 2 d y d t = Q 2 r | v v ˜ 2 r ( t ) | 2 d y d t ( sup t 0 4 r 2 < t t 0 i = 1 2 B 2 r , i | v i v ˜ 2 r ( t ) | 2 d y ) 1 2 t 0 4 r 2 t 0 ( B 2 r | v v ˜ 2 r ( t ) | 2 d y ) 1 2 d t C { i = 1 2 Q 4 r , i | v i | 2 d y d t + r ( n + 2 ) ( 1 2 p ) 4 n i = 1 2 ( Q 4 r , i | f i | n p n + p d y d t ) 2 ( n + p ) n p } 1 2 × t 0 4 r 2 t 0 ( B 2 r | v v ˜ 2 r ( t ) | 2 n n 2 d y ) n 2 4 n ( B 2 r | v v ˜ 2 r ( t ) | 2 n n + 2 d y ) n + 2 4 n d t ,则

s = 2 n n 2 , q = 2

其中

( B 2 r | v v ˜ 2 r ( t ) | 2 n n 2 d y ) n 2 4 n C ( B 2 r | v | 2 d y ) 1 4 .

s = q = 2 n n + 2

( B 2 r | v v ˜ 2 r ( t ) | 2 n n + 2 d y ) n + 2 4 n C r 1 2 ( B 2 r | v | 2 n n + 2 d y ) n + 2 4 n .

由于

r 1 2 t 0 4 r 2 t 0 ( B 2 r | v | 2 d y ) 1 4 ( B 2 r | v | 2 n n + 2 d y ) n + 2 4 n d t C r 1 2 { i = 1 2 Q 4 r , i | v i | 2 d y d t + r ( n + 2 ) ( 1 2 p ) 4 n i = 1 2 ( Q 4 r , i | f ˜ i | n p n + p d y d t ) 2 ( n + p ) n p } 1 4 r 1 1 n ( Q 2 r | v | 2 n n + 2 d y d t ) n + 2 4 n ,

i = 1 2 Q 2 r , i | v i v ˜ 2 r ( t ) | 2 d y d t C r 3 2 1 n { i = 1 2 Q 4 r , i | v i | 2 d y d t + r ( n + 2 ) ( 1 2 p ) 4 n i = 1 2 ( Q 4 r , i | f ˜ i | n p n + p d y d t ) 2 ( n + p ) n p } 3 4 ( Q 2 r | v | 2 n n + 2 d y d t ) n + 2 4 n . (2.3)

i = 1 2 Q r , i | v i | 2 d y d t C r 1 2 1 n { i = 1 2 Q 4 r , i | v i | 2 d y d t + r ( n + 2 ) ( 1 2 p ) 4 n i = 1 2 ( Q 4 r , i | f ˜ i | n p n + p d y d t ) 2 ( n + p ) n p } 3 4 ( Q 2 r | v | 2 n n + 2 d y d t ) n + 2 4 n + C r ( n + 2 ) ( 1 2 p ) 4 n i = 1 2 ( Q 4 r , i | f ˜ i | n p n + p d y d t ) 2 ( n + p ) n p ε i = 1 2 Q 4 r , i | v i | 2 d y d t + C 1 i = 1 2 Q 4 r , i | f ˜ i | n p n + p d y d t + C ( ε ) r 2 ( n + 2 ) n ( i = 1 2 Q 4 r , i | v i | 2 n n + 2 d y d t ) n + 2 n ,

C 1 充分小,对于任意 i = 1 2 f ˜ i L n p n + p ( Q R , i ) ,积分可得

| Q r | (2.4)

对于 ,在 | Q 4 r | = 4 n + 2 | Q r | 上关于t取上确界,有 1 | Q r | i = 1 2 Q r , i | v i | 2 d y d t ε 4 n + 2 1 | Q 4 r | i = 1 2 Q 4 r , i | v i | 2 d y d t + C 1 1 | Q 4 r | i = 1 2 Q 4 r , i | f ˜ i | n p n + p d y d t + C ( ε ) ( 1 | Q 4 r | i = 1 2 Q 4 r , i | v i | 2 n n + 2 d y d t ) n + 2 n . 成立。

如果 ε ( 0 , 1 ] 是问题的弱上解,同理可证 ε 4 n + 2 1 2

注2.4.1 如果 r ( 0 , R ] ,则由 可得

1 | Q r | Q r | v | 2 d y d t 1 2 1 | Q 4 r | Q 4 r | v | 2 d y d t + C 1 1 | Q 4 r | Q 4 r | f ˜ | n p n + p d y d t + C ( 1 | Q 4 r | Q 4 r | v | 2 n n + 2 d y d t ) n + 2 n . (2.5)

定理2.4.2 (局部极值原理)令 m = n + 2 n , m 1 = q 1 ( n + p ) ( n + 2 ) n 2 p , F = | v | 2 n n + 2 , G = ( | f ˜ | n p n + p ) n n + 2 . ,常数 F L l o c m ( Q R ) , G L l o c m 1 ( Q R ) ( 1 < m < m 1 ) . ,则对于任意 1 | Q r | Q r F m d y d t 1 2 1 | Q 4 r | Q 4 r F m d y d t + C { 1 | Q 4 r | Q 4 r G m d y d t + ( 1 | Q 4 r | Q 4 r F d y d t ) m } , ,有

β 1 , β 2 , n , λ , Λ , c 0 , p , i = 1 2 f ˜ i L n p n + p ( Q R , i )

如果 m 0 ( m , m 1 ] ,则有

F L l o c m 0 ( Q R )

其中C仅依赖于 1 | Q r | Q r F m 0 d y d t C { 1 | Q 4 r | Q 4 r G m 1 d y d t + ( 1 | Q 4 r | Q 4 r F m d y d t ) m 0 m } . 的参数, q = 2 n m 0 n + 2 ,并且 q > 2

证明 下面只证明第一种情况。令 Q r | v | q d y d t C { Q 4 r | f ˜ | q 1 d y d t + r ( n + 2 ) ( 1 q 2 ) ( Q 4 r | v | 2 d y d t ) q 2 } . L q ,其中 Q i = Ω i × ( 0 , T ] ( i = 1 , 2 ) 待定。令 L q ,且取 ( i = 1 , 2 ) C μ , k ( Ω ¯ i × [ 0 , T ] ) = { g ( x , t ) : g ( x , t ) 上的截断函数,并且满足

C k

应用公式(2.5)+,有

μ

μ [ 0 , 1 ] , k + } 。首先取 g ( x , t ) C μ , k ( Ω ¯ i × [ 0 , T ] ) = s = 0 k sup Ω ¯ i × [ 0 , T ] | D t s g ( x , t ) | + sup x , y Ω i , t [ 0 , T ] | g ( x , t ) g ( y , t ) | | x y | μ . ,则有

K ( x , t )

K i ( x , t ) C μ , ( Ω ¯ i × [ 0 , T ] ) ,则有

μ ( 0 , 1 )

又因为 C ,则

| K i ( x , t ) K i ( y , t ) | C | x y | μ , ( x , t ) , ( y , t ) Ω i × ( 0 , T ) .

故有

l 1

所以

φ m + 1 i = 1 2 ξ m ( u i k m + 1 ) + L 2 ( t 0 ρ m 2 , t 0 ; L 2 ( B ρ m , i ) ) 2 C | A m ( k m + 1 ) | 2 n i = 1 2 ξ m ( u i k m + 1 ) + L 2 ( t 0 ρ m 2 , t 0 ; L 2 * ( B ρ m , i ) ) 2 C | A m ( k m + 1 ) | 2 n ( C 2 2 m R 2 φ m + C k 2 R 2 ( 1 n p ) | A m ( k m + 1 ) | 1 2 p ) C 2 2 m ( 1 + 2 n ) ( φ m 1 + 2 n R 2 k 4 n + φ m 1 2 p + 2 n R 2 ( 1 n p ) k 4 n 4 p ) .

s = 0 l | D t s K i ( x , t ) | Λ 2 l , Ω i × ( 0 , T ) ; ,则

s = 0 l | D t s K i ( x , t ) D t s K i ( y , t ) | Λ 2 l | x y | μ , Ω i × ( 0 , T ) .

因此

ϵ > 0

Ω ϵ = { x Ω | d i s t ( x , Ω ) > ϵ } . ,则

f ˜ L ( Q T )

由引理2.3.1可得,如果

u ˜ W ( 0 , T ; V ) (2.6)

其中 0 < ϵ < 1 2 ,即

α * < min { μ , α 2 ( 1 + α ) } (2.7)

u ˜ L ( Ω ϵ × ( ϵ T , T ) ) + x u ˜ C α * , 0 ( ( Ω ϵ Ω ¯ i ) × ( ϵ T , T ) ) C ( u ˜ L 2 ( Q T ) + f ˜ L ( Q T ) ) , ,所以在 n , β 1 , β 2 , c 0 , λ , Λ , C , μ , α , ϵ , T , K i C α * , 1 ( Ω ¯ i × [ 0 , T ] ) 中,我们有 Ω i 。现在取k满足 C 1 + α ,则条件(2.7)成立,即(2.6)满足。总结以上关于k的选取,最终取k为

( i = 1 , 2 )

所以结论得证。

有了局部极值原理之后,我们总假设弱解在界面附近是有界的。在此前提之下,将会得到弱解的一些局部性质。为了适应方程的需要,下面对De Giorgi类的定义稍加修改。

定义2.4.2我们称定义于 x u ˜ L ( Ω ϵ × ( ϵ T , T ) ) C ( u ˜ L 2 ( Q T ) + f ˜ L ( Q T ) ) . 上的函数 L q 属于De Giorgi类,如果 q > 2 L q L q ,并且存在 Ω n ( n > 2 ) ,使得对于任意k满足:

Ω (2.8)

Ω 1 Ω 满足 ,记 Γ C 2 ,同样可定义 Ω 2 = Ω Ω 1

引理2.4.1 (De Giorgi lemma) (文献 [13] )设 Γ = Ω 1 ,记 Γ ,则对于 Ω ,有

u = ( u 1 , u 2 )

其中 u 1 : Ω 1 × ( 0 , T ] , u 2 : Ω 2 × ( 0 , T ] 只依赖于n。

定理2.4.3令 u ( x , t ) n Γ Ω 1 ,则存在 Ω 2 ,使得对于 i = 1 , 2 ,如果

w i = w i ( x ) (2.9)

K i = K i ( x , t ) (2.10)

m Γ = m Γ ( x , t )

其中 u i t + w i u i d i v ( K i ( x , t ) u i ) = f i ( x , t ) , x Ω i , i = 1 , 2 , t ( 0 , T ] ; m Γ ( x , t ) [ u ] Γ = K 1 ( x , t ) u 1 n , x Γ , t ( 0 , T ] ; K 1 ( x , t ) u 1 n = K 2 ( x , t ) u 2 n , x Γ , t ( 0 , T ] ; u i ( , 0 ) = u 0 i , x Ω i , i = 1 , 2 ; u 2 ( , t ) = 0 , x Ω , t ( 0 , T ] . 仅依赖于 [ u ] Γ = u 2 | Γ u 1 | Γ 的参数,并且 u 0 ( x ) = ( u 0 1 , u 0 2 )

证明令 u 0 ( x ) u i t + w i u i d i v ( K i ( x , t ) u i ) = f i ( x , t ) , x Ω i , i = 1 , 2 , t ( 0 , T ] ; K 1 ( x , t ) u 1 n = K 2 ( x , t ) u 2 n , x Γ , t ( 0 , T ] ; [ β u ] Γ = 0 , x Γ , t ( 0 , T ] ; u i ( , 0 ) = u 0 i , x Ω i , i = 1 , 2 ; u 2 ( , t ) = 0 , x Ω , t ( 0 , T ] . 。记 β ,取与定理2.4.2中相同的截断函数 Ω i 。在 +中分别取 β = β i > 0 代替 β 1 β 2 ,则

[ β u ] Γ = β 2 u 2 | Γ β 1 u 1 | Γ

u 0 ( x ) ,则

p [ 1 , ]

又因为

W 1 , p ( Ω )

所以

H 1 ( Ω ) : = W 1 , 2 ( Ω )

由条件 可得, { V 1 = H 1 ( Ω 1 ) ; V 2 = { v 2 H 1 ( Ω 2 ) | v 2 | Ω = 0 } ; V = V 1 × V 2 , v V 2 = v 1 H 1 ( Ω 1 ) 2 + v 2 H 1 ( Ω 2 ) 2 ,且 ,则有

H = L 2 ( Ω 1 ) × L 2 ( Ω 2 )

W ( 0 , T ; V ) : = { v L 2 ( 0 , T ; V ) | v t L 2 ( 0 , T ; V ) } ,且 v H 2 = v 1 L 2 ( Ω 1 ) 2 + v 2 L 2 ( Ω 2 ) 2 , v W 2 = v L 2 ( 0 , T ; V ) 2 + v t L 2 ( 0 , T ; V ) 2 . ,则

W ( 0 , T ; V ) C ( [ 0 , T ] ; H ) ,

注意 ,由引理2.3.1可得,存在 tr Γ i : V i L 2 ( Γ ) , ( i = 1 , 2 ) , 使得当 tr Γ u = ( tr Γ 1 u 1 , tr Γ 2 u 2 ) ,有 tr Γ u L 2 ( Γ ) 2 = tr Γ 1 u 1 L 2 ( Γ ) + tr Γ 2 u 2 L 2 ( Γ ) 。所以在 Ω 1 中, Ω 2 。因此结论得证。

定理2.4.4令 Ω Lip Γ = Ω 1 C 2 ,如果 H n 1 Γ n 1 ,对于 w i = w i ( x ) ( i = 1 , 2 ) Ω i w i ( ) [ H 1 ( Ω i ) ] n 满足

w i ( ) L ( Ω i ) c 0 < (2.11)

其中 Ω i d i v w i = 0 ,则存在 Γ ,使得或者

w 1 n = w 2 n = 0 (2.12)

或者

K i = K i ( x , t ) ( i = 1 , 2 ) (2.13)

成立,其中 m Γ = m Γ ( x , t ) ,s仅依赖于 K i ( , t ) = ( k p q i ( , t ) ) p q ,并且 Ω i

在证明定理2.4.4之前,先给出以下两个辅助引理。

引理2.4.2令 K i λ , Λ > 0 ,如果 p , q = 1 n k p q i L ( Ω i ) Λ , a . e . t [ 0 , T ] ; ξ n Ω i λ | ξ | 2 ξ K i ( , t ) ξ Λ | ξ | 2 , a . e . t [ 0 , T ] . ,对于 m Γ = m Γ ( , t ) Γ m ¯ m _ > 0 满足

m _ m Γ ( x , t ) m ¯ , a . e . ( x , t ) Γ × [ 0 , T ] . (2.14)

其中 M = ( m Γ m Γ m Γ m Γ ) r = ( r 1 , r 2 ) 2 ,则对于任意正整数 r M r = a . e . 0 r 1 = r 2 . ,或者

f i L 2 ( 0 , T ; L 2 ( Ω i ) ) (2.15)

或者

v = ( v 1 , v 2 ) V (2.16)

成立,其中 i = 1 2 Ω i u i x v i + w i ( u i ) v i + ( v i ) K i u i d x + Γ ( tr Γ u ) M ( tr Γ v ) d H n 1 = i = 1 2 Ω i f i v i d x . d u d t ( v ) = i = 1 2 Ω i u i t v i d x , ,C依赖于 a ( t ; u ( t ) , v ) = i = 1 2 Ω i w i ( u i ) v i + ( v i ) K i u i d x + Γ ( tr Γ u ) M ( tr Γ v ) d H n 1 ,

证明对于 F ( t ) ( v ) = i = 1 2 Ω i f i v i d x , ,记 u 0 ( x ) H , F L 2 ( 0 , T ; V ) u W ( 0 , T ; V ) ,所以 { d u d t ( v ) + a ( t ; u ( t ) , v ) ( ) F ( t ) ( v ) , v V ; u ( 0 ) = u 0 ( x ) . 。取 u ,则 u 。对于 a ( t ; u ( t ) , v ) ,由引理2.4.1可以得到

V × V (2.17)

应用条件(2.14)可知, a ( t ; u ( t ) , v ) ,则 a ( t ; u ( t ) , v ) ,因此对 应用Hölder不等式,可以得到

v V (2.18)

对 两边关于 a ( t ; v , v ) = i = 1 2 Ω i w i ( v i ) v i + ( v i ) K i v i d x + Γ ( tr Γ v ) M ( tr Γ v ) d H n 1 , 积分,进而可得

i = 1 2 Ω i w i ( v i ) v i d x = 0. (2.19)

a ( t ; v , v ) min { λ , m _ } [ v 1 L 2 ( Ω 1 ) 2 + v 2 L 2 ( Ω 2 ) 2 + Γ ( v 1 v 2 ) 2 d H n 1 ] . v 2 H 1 ( Ω 2 ) 2 c 2 ( v 2 L 2 ( Ω 2 ) 2 + Ω v 2 2 d s ) = c 2 v 2 L 2 ( Ω 2 ) 2 , ,并且在 两边对i求和,有

v 1 H 1 ( Ω 1 ) 2 c 1 ( v 1 L 2 ( Ω 1 ) 2 + Γ v 1 2 d H n 1 ) ,

Γ v 1 2 d H n 1 2 Γ ( v 1 v 2 ) 2 d H n 1 + 2 Γ v 2 2 d H n 1 2 Γ ( v 1 v 2 ) 2 d H n 1 + 2 c 3 2 v 2 H 1 ( Ω 2 ) 2 , a ( t ; v , v ) min { λ , m _ } max { c 1 , 2 c 1 , 2 c 3 2 c 1 ( c 2 + 1 ) } v V 2 = α 0 v V 2 , 上的截断函数,使得在 c 1 , c 2 , c 3 上有 X = ( x , t X ) , Y = ( y , t Y ) Q T 。由于 δ ˜ ( X , Y ) = max { | x y | , | t x t y | 1 2 } . ,类似定理2.4.1的证明,取 n + 1 X D ,则有

D ( X , r ) = D Q r ( X )

其中 Q r ( X ) = B r ( x ) × ( t X r 2 , t X + r 2 ) 。如果 不成立,则 d = d i a m ( D ) ,进而可得 p 1 , θ 0 ,所以

M p , θ ( D ; δ ˜ ) (2.20)

对(2.20)关于l从 L p ( D ) u M p , θ ( D ; δ ˜ ) : = ( sup X D ¯ , d ρ > 0 ρ θ D ( X , ρ ) | u ( Y ) | p d Y ) 1 p < , 求和,可以得到

Ω n

因此

x Ω

引理2.4.3 令 Ω ( x , r ) = Ω B r ( x ) d 1 = d i a m ( Ω ) ,如果 d * ( x , y ) = | x y | M p , θ ( Ω ; d * ) L p ( Ω ) ,对于 u M p , θ ( Ω ; d * ) : = ( sup x Ω ¯ , d 1 ρ > 0 ρ θ Ω ( x , ρ ) | u ( y ) | p d y ) 1 p < , p 1 , θ 0 L p , θ ( D ; δ ˜ ) 满足

L p ( D ) (2.21)

其中 [ u ] L p , θ ( D ; δ ˜ ) : = ( sup X D , d ρ > 0 ρ θ D ( X , ρ ) | u ( Y ) u X , ρ | p d Y ) 1 p < , u L p , θ ( D ; δ ˜ ) = ( u L p ( D ) p + [ u ] L p , θ ( D ; δ ˜ ) p ) 1 p , ,则存在 u X , ρ 使得或者

D ( X , ρ ) (2.22)

或者

u X , ρ = 1 | D ( X , ρ ) | D ( X , ρ ) u ( Y ) d Y (2.23)

成立,其中 0 < α 1 仅依赖于 C α ( D ¯ ; δ ˜ ) ,并且 [ u ] α ; D : = sup X , Y D , X Y | u ( X ) u ( Y ) | δ ˜ ( X , Y ) α < ,

证明 令 u α ; D = sup D | u | + [ u ] α ; D . C α ( D ¯ ; δ ˜ ) C ( D ¯ ) . 上的截断函数,使得 C 1 + α ( D ¯ ; δ ˜ ) 上有 C 1 ( D ¯ ; δ ˜ ) ,其中 [ u ] 1 + α ; D : = i = 1 n sup X , Y D , X Y | D i u ( X ) D i u ( Y ) | δ ˜ ( X , Y ) α + sup ( x , t 1 ) , ( x , t 2 ) D , t 1 t 2 | u ( x , t 1 ) u ( x , t 2 ) | | t 1 t 2 | 1 + α 2 < , 未知,对于 u 1 + α ; D = sup D | u | + i = 1 n sup D | D i u | + [ u ] 1 + α ; D . ,在 X D , 0 < ρ d i a m ( D ) 上类似定理2.4.1的证明,记 | D ( X , ρ ) | A | Q ρ ( X ) | . ,故有

sup t 0 < t t 0 + a R 2 i = 1 2 ξ ( u i k ) + ( , t ) L 2 ( B R , i ) 2 ( 1 + ε ) i = 1 2 ξ ( u i k ) + ( , t 0 ) L 2 ( B R , i ) 2 + γ ( ε ) [ C H 2 ( 1 β ) 2 R 2 | A R a ( k ) | + ( M + F 0 , R a ) 2 | A R a ( k ) | 1 2 p ] ,

其中 n + 2 < θ n + 2 + p 。对于任意整数 L p , θ ( D ; δ ˜ ) C α ( D ¯ ; δ ˜ ) , ,我们有

α = θ ( n + 2 ) p

如果(2.22)不成立,并且应用条件(2.21),故可以得到

A B

其中 A B 。显然对于 B A ,可得

Q T = Ω × ( 0 , T ] , Q i = Ω i × ( 0 , T ] ( i = 1 , 2 ) ;

首先取 p Q T = { ( x , t ) | x Ω ¯ , t = 0 } ( Ω × ( 0 , T ] ) ; 使得 sup p Q T u = max { sup p Q T u 1 , sup p Q T u 2 } ; ,则

v + : = max { v , 0 } , v : = ( v ) + = max { v , 0 } .

由于 y h ( h = 0 , 1 , ) 是关于 y h + 1 C b h y h 1 + ε 的严格递增函数,令 b > 1 是其逆函数。当 ε > 0 ,有 y 0 θ : = C 1 ε b 1 ε 2 ,则 lim h y h = 0 。取 φ ( t ) ,故有

[ k 0 , ) (2.24)

又因为 C > 0 , α > 0 , β > 1 ,取 h > k k 0 使得 φ ( h ) C ( h k ) α [ φ ( k ) ] β 。对于如此确定的常数a,不妨设 d C 1 α [ φ ( k 0 ) ] β 1 α 2 β β 1 是整数N,记 φ ( k 0 + d ) = 0

我们将会通过数学归纳法证明以下结论:存在 u W ( 0 , T ; V ) 使得

p > n + 2 (2.25)

其中 f i L 2 ( 0 , T ; L n p n + p ( Ω i ) ) ( i = 1 , 2 ) 仅依赖于 max { ess sup Q 1 u 1 , ess sup Q 2 u 2 } sup p Q T u + + C F 0 | Ω | 1 n 1 p ,

Step (1) 当 n , λ , m _ , p , T ,已知

sup p Q T u + = max { sup p Q T u 1 + , sup p Q T u 2 + }

F 0 = i = 1 2 f i L 2 ( 0 , T ; L n p n + p ( Ω i ) ) < 充分大,可以得到

k 0 = sup p Q T u +

应用引理2.4.2,对于任意 k k 0 ,有

v 1 = ( u 1 k ) + , v 2 = ( u 2 k ) + (2.26)

其中 v 2 | Ω × ( 0 , T ] = 0 ,C只依赖于 v = ( v 1 , v 2 ) | t = 0 = 0 。在 中取 u 分别代替 i = 1 2 Ω i u i t v i + w i ( u i ) v i + ( v i ) K i u i d x + Γ m Γ ( u 1 u 2 ) ( v 1 v 2 ) d H n 1 i = 1 2 Ω i f i v i d x , ,并且将不等式 代入 的右侧,则有

i = 1 2 Ω i v i t v i + w i ( v i ) v i + ( v i ) K i v i d x + Γ m Γ ( v 1 v 2 ) 2 d H n 1 i = 1 2 Ω i f i v i d x . (2.27)

首先取 i = 1 2 Ω i v i t v i + λ | v i | 2 d x + m _ Γ ( v 1 v 2 ) 2 d H n 1 i = 1 2 Ω i f i v i d x . 足够大使得 的右端方括号中的第二项不大于 λ i = 1 2 v i L 2 ( Ω i ) 2 + m _ Γ ( v 1 v 2 ) 2 d H n 1 C [ v 1 L 2 * ( Ω 1 ) 2 + v 2 L 2 * ( Ω 2 ) 2 ] , ,然后取 2 * = 2 n n 2 使得方括号中的第一项不大于 Ω i [ u i ( , t ) > k ] = { x Ω i | u i ( x , t ) > k } ( i = 1 , 2 ) ,对于选定的 i = 1 2 Ω i v i t v i d x + C i = 1 2 v i L 2 * ( Ω i ) 2 i = 1 2 Ω i f i v i d x i = 1 2 ( ε v i L 2 * ( Ω i ) 2 + C ( ε ) f i L n p n + p ( Ω i ) 2 | Ω i [ u i ( , t ) > k ] | 1 2 p ) , ,我们有

ε

Step (2) 在 i = 1 2 Ω i v i t v i d x + C 0 i = 1 2 v i L 2 * ( Ω i ) 2 C ( ε ) i = 1 2 f i L n p n + p ( Ω i ) 2 | Ω i [ u i ( , t ) > k ] | 1 2 p 上应用 +,然后重复以上步骤,可证明(2.25)对 i = 1 2 v i ( , T ) L 2 ( Ω i ) 2 + C 0 i = 1 2 v i L 2 ( 0 , T ; L 2 * ( Ω i ) ) 2 C F 0 2 i = 1 2 | Q i [ u i > k ] | 1 2 p 成立。

因此有结论成立。

定理2.4.4的证明:

根据定理2.4.4的条件,引理2.4.3成立,令 i = 1 2 v i L 2 ( 0 , T ; L 2 * ( Ω i ) ) 2 C F 0 2 i = 1 2 | Q i [ u i > k ] | 1 2 p . 是由引理2.4.3确定的常数。对于待定的 h > k ,令(2.12)不成立,则(2.15),(2.22)也不成立,而(2.10)成立。所以由引理2.4.3可得,

i = 1 2 v i L 2 ( 0 , T ; L 2 * ( Ω i ) ) 2 = i = 1 2 0 T ( Ω i | v i | 2 * d x ) 2 2 * d t C ( h k ) 2 i = 1 2 | Q i [ u i > h ] | n 2 n ,

因此(2.14)成立。所以可由引理2.4.2得到以下不等式,

( h k ) 2 i = 1 2 | Q i [ u i > h ] | n 2 n C F 0 2 i = 1 2 | Q i [ u i > k ] | 1 2 p .

取s充分大使得 ψ ( k ) = i = 1 2 | Q i [ u i > k ] | ,故(2.9)成立。则由定理2.4.3可得以下结论,

ψ ( h ) ( C F 0 h k ) 2 * [ ψ ( k ) ] n ( p 2 ) p ( n 2 ) .

对于 d = ( C F 0 ) n 2 2 n | ψ ( k 0 ) | 1 n 1 p 2 n ( p 2 ) 2 ( p n ) C F 0 | Q T | 1 n 1 p ,我们有相似的性质。

定理2.4.5令 ψ ( k 0 + d ) = 0 i = 1 2 | Q i [ u i > sup p Q T u + + d | = 0. ,如果 Q i u i sup p Q T u + + C F 0 | Q T | 1 n 1 p max { ess sup Q 1 u 1 , ess sup Q 2 u 2 } sup p Q T u + + C F 0 | Ω | 1 n 1 p . ,对于 u W ( 0 , T ; V ) p > n + 2 f i L 2 ( 0 , T ; L n p n + p ( Ω i ) ) ( i = 1 , 2 ) 满足

max { u 1 L ( Q 1 ) , u 2 L ( Q 2 ) } max { sup p Q T | u 1 | , sup p Q T | u 2 | } + C F 0 | Ω | 1 n 1 p , (2.28)

其中 n , λ , m _ , p , T F 0 = i = 1 2 f i L 2 ( 0 , T ; L n p n + p ( Ω i ) ) < ,则存在 Q T 使得或者

u (2.29)

或者

u W ( 0 , T ; V ) (2.30)

成立。

2.5. 弱解的局部Hölder连续性

下面考虑 X 0 = ( x 0 , t 0 ) Γ × ( 0 , T ] 在界面上 Q ρ , τ ( X 0 ) = B ρ ( x 0 ) × ( t 0 , t 0 + τ ] Q T 和界面与初值层 k , ε ( 0 , 1 ] 的相交处的Hölder连续性。

引理2.5.1 (文献 [14] , p. 140引理4.1)令 ξ ( x , t ) C ( [ t 0 , t 0 + τ ] ; C 0 ( B ρ ( x 0 ) ) ) 是定义于 0 ξ 1 上的非减非负函数,如果它满足

ξ ( , t 0 ) = 0

其中 sup t 0 < t t 0 + τ i = 1 2 ξ ( u i k ) ± ( , t ) L 2 ( B ρ , i ) 2 + λ 1 i = 1 2 ( ξ ( u i k ) ± ) L 2 ( t 0 , t 0 + τ ; L 2 ( B ρ , i ) ) 2 + m 1 t 0 t 0 + τ Γ B ρ ( ξ ( u 1 k ) ± ξ ( u 2 k ) ± ) 2 d H n 1 d t C * [ ( ξ t L ( Q ρ , τ ) + ξ L ( Q ρ , τ ) 2 ) i = 1 2 ( u i k ) ± L 2 ( t 0 , t 0 + τ ; L 2 ( B ρ , i ) ) 2 + F 0 , ρ , τ 2 i = 1 2 | Q ρ , τ , i [ ( u i k ) ± > 0 ] | 1 2 p ] , B ρ , i = B ρ Ω i , Q ρ , τ , i = Q ρ , τ Q i 0 < ρ , τ < 1 是常数,则存在 p > n + 2 λ 1 > 0 使得

λ

其中 m 1 = 2 m _ 仅依赖于 F 0 , ρ , τ = i = 1 2 f i L 2 ( t 0 , t 0 + τ ; L n p n + p ( B ρ , i ) ) > 0

定理2.5.1令 C * n , Λ , p , c 0 D G ( Q T ) = D G ( Q T ; λ 1 , m 1 , p , n , F 0 , ρ , τ , C * ) u W ( 0 , T ; V ) ,则对于 u D G + ( Q T ) ,有

u W ( 0 , T ; V )

其中 u D G ( Q T ) D G ( Q T ) = D G + ( Q T ) D G ( Q T ) 仅依赖于 u W ( 0 , T ; V ) 的参数, p > n + 2 ,并且 f i L 2 ( 0 , T ; L n p n + p ( Ω i ) ) ( i = 1 , 2 ) u D G + ( Q T ) u W ( 0 , T ; V )

证明 记

u D G ( Q T )

C *

对于 n , Λ , p , c 0 F 0 , ρ , τ = i = 1 2 f i L 2 ( t 0 , t 0 + τ ; L n p n + p ( B ρ , i ) ) < ,可以看出以下两种情况中必有一种成立 u W ( 0 , T ; V )

ξ ( x , t ) C ( [ t 0 , t 0 + τ ] ; C 0 ( B ρ ( x 0 ) ) ) (2.31)

Q ρ , τ ( X 0 ) (2.32)

Step (1) 如果(2.31)成立,则由定理2.4.4可得,存在 0 ξ ( x , t ) 1 使得当 ξ ( , t 0 ) = 0 时,

B ρ , i = B ρ Ω i (2.33)

或者

v i = ξ 2 ( u i k ) + 0 ( i = 1 , 2 ) (2.34)

成立,其中 i = 1 2 B ρ , i ( u i k ) + t ξ 2 ( u i k ) + + w i ( ( u i k ) + ) ξ 2 ( u i k ) + + ( ( ξ 2 ( u i k ) + ) ) K i ( u i k ) + d x + m _ Γ B ρ ( ξ ( u 1 k ) + ξ ( u 2 k ) + ) 2 d H n 1 i = 1 2 B ρ , i f i ξ 2 ( u i k ) + d x . ,并且由 可得当 i = 1 2 B ρ , i w i ( ( u i k ) + ) ξ 2 ( u i k ) + d x = i = 1 2 B ρ , i w i ( u i k ) + ( ξ ) ( ξ ( u i k ) + ) d x , 时,有 i = 1 2 B ρ , i ( ( ξ 2 ( u i k ) + ) ) K i ( u i k ) + d x λ i = 1 2 ( ξ ( u i k ) + ) L 2 ( B ρ , i ) 2 i = 1 2 B ρ , i ( ( u i k ) + ) 2 ( ξ ) K i ξ d x .

Step (2) 如果(2.32)成立,则由定理2.4.5可得,存在 i = 1 2 B ρ , i ( ξ ( u i k ) + ) t ξ ( u i k ) + d x + λ i = 1 2 ( ξ ( u i k ) + ) L 2 ( B ρ , i ) 2 + m _ Γ B ρ ( ξ ( u 1 k ) + ξ ( u 2 k ) + ) 2 d H n 1 i = 1 2 B ρ , i ξ t ξ ( ( u i k ) + ) 2 d x + i = 1 2 B ρ , i w i ( u i k ) + ( ξ ) ( ξ ( u i k ) + ) d x + i = 1 2 B ρ , i ( ( u i k ) + ) 2 ( ξ ) K i ξ d x + i = 1 2 B ρ , i f i ξ 2 ( u i k ) + d x . 使得当 i = 1 2 B ρ , i w i ( u i k ) + ( ξ ) ( ξ ( u i k ) + ) d x ε i = 1 2 ( ξ ( u i k ) + ) L 2 ( B ρ , i ) 2 + γ ( ε ) ξ L ( Q ρ , τ ) 2 i = 1 2 ( u i k ) + L 2 ( B ρ , i ) 2 , 时,

ξ ( u i k ) + L 2 * ( B ρ , i ) 2 C ( ( ξ ( u i k ) + ) L 2 ( B ρ , i ) 2 + B ρ , i ( Γ B ρ ) | ξ ( u i k ) + | 2 d H n 1 ) = C ( ξ ( u i k ) + ) L 2 ( B ρ , i ) 2 , (2.35)

或者

i = 1 2 B ρ , i f i ξ 2 ( u i k ) + d x i = 1 2 f i L n p n + p ( B ρ , i ) ξ ( u i k ) + L 2 * ( B ρ , i ) | B ρ , i [ u i ( , t ) > k ] | 1 2 1 p ε i = 1 2 ( ξ ( u i k ) + ) L 2 ( B ρ , i ) 2 + γ ( ε ) i = 1 2 f i 2 L n p n + p ( B ρ , i ) | B ρ , i [ u i ( , t ) > k ] | 1 2 p . (2.36)

成立,并且由(2.36)可得当 ε 时,有 t ( t 0 , t 0 + τ ]

已知 i = 1 2 ξ ( u i k ) ± ( , t ) L 2 ( B ρ , i ) 2 + λ 1 i = 1 2 ( ξ ( u i k ) ± ) L 2 ( t 0 , t ; L 2 ( B ρ , i ) ) 2 + m 1 t 0 t Γ B ρ ( ξ ( u 1 k ) ± ξ ( u 2 k ) ± ) 2 d H n 1 d t C * [ ( ξ t L ( Q ρ , τ ) + ξ L ( Q ρ , τ ) 2 ) i = 1 2 ( u i k ) ± L 2 ( t 0 , t 0 + τ ; L 2 ( B ρ , i ) ) 2 + F 0 , ρ , τ 2 i = 1 2 | Q ρ , τ , i [ ( u i k ) ± > 0 ] | 1 2 p ] . ,则当 ( t 0 , t 0 + τ ] 时, + ,即(2.33), 成立,故 u W ( 0 , T ; V ) ,所以 u D G ( Q T ) 。因此,当 u D G ( Q T ) ,我们有

sup t 0 < t t 0 + τ i = 1 2 ξ ( u i k ) ± ( , t ) L 2 ( B ρ , i ) 2 + λ ¯ i = 1 2 ξ ( u i k ) ± L 2 ( t 0 , t 0 + τ ; L 2 * ( B ρ , i ) ) 2 + m 1 t 0 t 0 + τ Γ B ρ ( ξ ( u 1 k ) ± ξ ( u 2 k ) ± ) 2 d H n 1 d t C * [ ( ξ t L ( Q ρ , τ ) + ξ L ( Q ρ , τ ) 2 ) i = 1 2 ( u i k ) ± L 2 ( t 0 , t 0 + τ ; L 2 ( B ρ , i ) ) 2 + F 0 , ρ , τ 2 i = 1 2 | Q ρ , τ , i [ ( u i k ) ± > 0 ] | 1 2 p ] .

其中 u D G + ( Q T ) 仅依赖于 X 0 = ( x 0 , t 0 ) Γ × ( 0 , T ] 。由引理2.5.1可得,

p > n + 2

其中 Q R ( X 0 ) = B R ( x 0 ) × ( t 0 R 2 , t 0 ] Q T , 0 < R 1 。所以结论成立。

定理2.5.2 令 max { ess sup Q R 2 , 1 u 1 , ess sup Q R 2 , 2 u 2 } C ( 1 R n i = 1 2 u i L 2 ( t 0 R 2 , t 0 ; L 2 ( B R , i ) ) + F 0 , R R 1 n p ) , ,Q与界面 u D G (

参考文献

[1] Gueribiz, D., Jacquemin, F. and Fréour, S. (2013) A Moisture Diffusion Coupled Model for Composite Materials. Eu-ropean Journal of Mechanics A—Solids, 42, 81-89.
https://doi.org/10.1016/j.euromechsol.2013.04.008
[2] Bonito, A., Nochetto, R.H. and Pauletti, M.S. (2011). Dynamics of Biomembranes: Effect of the Bulk Fluid. Mathematical Mod-elling of Natural Phenomena, 6, 25-43.
https://doi.org/10.1051/mmnp/20116502
[3] Voulis, I. and Reusken, A. (2018) A Time Dependent Stokes Interface Problem: Well-Posedness and Space-Time Finite Element Discretization. ESAIM: Mathematical Modelling and Numerical Analysis, 52, 2187-2213.
https://doi.org/10.1051/m2an/2018053
[4] Hansbo, A. and Hansbo, P. (2002) An Unfitted Finite Element Method, Based on Nitsche’s Method, for Elliptic Interface Problems. Computer Methods in Applied Mechanics and Engineering, 191, 5537-5552.
https://doi.org/10.1016/S0045-7825(02)00524-8
[5] Jia, D.X., Sheng, Z.Q. and Yuan, G.W. (2019) An Extre-mum-Preserving Iterative Procedure for the Imperfect Interface Problem. Communications in Computational Physics, 25, 853-870.
https://doi.org/10.4208/cicp.OA-2017-0222
[6] Zhou, H.F., Sheng, Z.Q. and Yuan, G.W. (2020) A Fi-nite Volume Method Preserving Maximum Principle for the Diffusion Equations with Imperfect Interface. Applied Nu-merical Mathematics, 158, 314-335.
https://doi.org/10.1016/j.apnum.2020.08.008
[7] Dillon, R. and Li, Z.L. (2009) An Introduction to the Immersed Boundary and the Immersed Interface Methods. In: Interface Problems and Methods in Biological and Physical Flows, World Scientific Publishing, Singapore, 1-67.
https://doi.org/10.1142/9789812837851_0001
[8] Cao, F.J., Sheng, Z.Q. and Yuan, G.W. (2018) Monotone Fi-nite Volume Schemes for Diffusion Equation with Imperfect Interface on Distorted Meshes. Journal of Scientific Compu-ting, 76, 1055-1077.
https://doi.org/10.1007/s10915-018-0651-8
[9] Gross, S. and Reusken, A. (2011) Numerical Methods for Two-Phase Incompressible Flows. Springer, Berlin.
https://doi.org/10.1007/978-3-642-19686-7
[10] 王耀东. 偏微分方程的 理论[M]. 北京: 北京大学出版社, 1989: 141-202.
[11] Campanato, S. (1963) Proprietà di hölderianità di alcune classi di funzioni. Annali della Scuola Normale Superiore di Pisa—Scienze Fisiche e Matematiche, 17, 175-188.
[12] Campanato, S. (1966) Equazioni para-boliche del secondo ordine e spazi . Annali di Matematica Pura ed Applicata, 73, 55-102.
https://doi.org/10.1007/BF02415082
[13] Ladyzhenskaya, O.A., Solonnikov, V.A. and Ural’Tseva, N.N. (1968) Linear and Quasi-Linear Equations of Parabolic Type. American Mathematical Society, Providence.
https://doi.org/10.1090/mmono/023
[14] 陈亚浙. 二阶抛物型偏微分方程[M]. 北京: 北京大学出版社, 2003: 117-140.
[15] Lieberman, G.M. (1996) Second Order Parabolic Differential Equations. World Scientific Publishing Company, Singapore.
https://doi.org/10.1142/3302
[16] Bögelein, V. (2007) Regularity Results for Weak and Very Weak Solutions of Higher Order Parabolic Systems. Ph.D. Thesis, Salzburg University, Salzburg.
[17] Giaquintav, M. (1983) Multiple Integrals in the Calculus of Variations and Nonlinear Elliptic Systems. (AM-105). Princeton University Press, Princeton.
https://doi.org/10.1515/9781400881628
[18] Li, H.G. and Li, Y.Y. (2017) Gradient Estimates for Parabolic Systems from Composite Material. Science China Mathematics, 60, 2011-2052.
https://doi.org/10.1007/s11425-017-9153-0
[19] Li, Y.Y. and Nirenberg, L. (2003) Estimates for Elliptic Systems from Composite Material. Communications on Pure and Applied Mathematics, 56, 892-925.
https://doi.org/10.1002/cpa.10079