Advances in Applied Mathematics
Vol. 12  No. 02 ( 2023 ), Article ID: 61880 , 34 pages
10.12677/AAM.2023.122068

一类抛物型界面问题的正则性分析

史雪婷1,2,杨田洁1,2

1中国工程物理研究院,北京

2北京应用物理与计算数学研究所,北京

收稿日期:2023年1月21日;录用日期:2023年2月20日;发布日期:2023年2月27日

摘要

界面问题用于各种工程应用和物理、化学、生物现象的建模,特别是涉及具有不同扩散性、密度、渗透性或导电性的多种不同材料的现象,其在界面上由一定条件耦合。本文考虑具有溶解质输运的线性两相流模型,分别在相界面处耦合非完美界面条件和Henry界面条件。由于解在界面上的跳跃使得解在各自材料区域上具有比在整个区域上更高的正则性,针对这类界面问题的正则性分析,本文给出一个完整的泛函分析过程,采用De Giorgi迭代方法证明该模型弱解的相关性质,进而证明弱解及其梯度的Hölder连续性。此外,对于Henry界面问题,本文给出了梯度的 ω ( R 4 ) 2 s + 1 ( M + F 0 ) R 1 n + 2 p + R α 1 ( [ u 0 1 ] C α 1 ( Ω ¯ 1 ) + [ u 0 2 ] C α 1 ( Ω ¯ 2 ) ) . 估计(存在q > 2)。

关键词

对流扩散方程,非完美界面条件,Henry界面条件,De Giorgi迭代,Hölder连续性

Regular Analysis of a Class of Parabolic Interface Problems

Xueting Shi1,2, Tianjie Yang1,2

1China Academy of Engineering Physics, Beijing

2Institute of Applied Physics and Computational Mathematics, Beijing

Received: Jan. 21st, 2023; accepted: Feb. 20th, 2023; published: Feb. 27th, 2023

ABSTRACT

Interface problems are used for various engineering applications and modeling of physical, chemical, and biological phenomena, especially those involving a number of different materials with different 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, coupling imperfect interface condition and Henry interface condition at the phase interface, respectively. 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 iteration 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 ω ( R 4 ) ess sup Q R 4 , 1 u 1 , k ( + ) ess inf Q R 4 , 1 u 1 , k ˜ ( ) + ess sup Q R 4 , 2 u 2 , k ( + ) ess inf Q R 4 , 2 u 2 , k ˜ ( ) 2 ( 1 1 2 s ) ω ( R ) + 1 2 s 1 R α 1 ( [ u 0 1 ] C α 1 ( Ω ¯ 1 ) + [ u 0 2 ] C α 1 ( Ω ¯ 2 ) ) . estimation of the gradient (q > 2 is present).

Keywords:Convective Diffusion Equation, Imperfect Interface Condition, Henry Interface Condition, De Giorgi Iteration, Hölder Continuity

Copyright © 2023 by author(s) and Hans Publishers Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

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 ( Q T ) 相交,并且 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 ) , ,则存在 D G ( Q T ) Q R , i = Q R Q i 使得对于 F 0 , R = i = 1 2 f i L 2 ( t 0 R 2 , t 0 ; L n p n + p ( B R , i ) ) <

ρ 0 = R , ρ m = R 2 + R 2 m + 1

其中 k 0 = k , k m = k ( 2 1 2 m ) ( m = 0 , 1 , 2 , ) k > 0 Q ρ m ( X 0 ) = B ρ m ( x 0 ) × ( t 0 ρ m 2 , t 0 ] 仅依赖于 ξ m ( x , t ) ,并且 Q ρ m ( X 0 )

证明 对于任意 { ξ m ( x , t ) C ( [ t 0 ρ m 2 , t 0 ] ; C 0 ( B ρ m ( x 0 ) ) ) , 0 ξ m 1 , ξ ( , t 0 ρ m 2 ) = 0 ; ξ m ( x , t ) = 1 , Q ρ m + 1 ; | ξ m t | + | ξ m | 2 C ( n ) ( ρ m ρ m + 1 ) 2 . ,记 i = 1 2 ξ m ( u i k m + 1 ) + L 2 ( t 0 ρ m 2 , t 0 ; L 2 * ( B ρ m , i ) ) 2 C 2 2 m R 2 i = 1 2 ( u i k m + 1 ) + L 2 ( t 0 ρ m 2 , t 0 ; L 2 ( B ρ m , i ) ) 2 + C F 0 , ρ m 2 i = 1 2 | Q ρ m , i [ u i > k m + 1 ] | 1 2 p . ,并且 A m ( k m + 1 ) = i = 1 2 ( Q ρ m , i [ u i > k m + 1 ] )

由定理2.5.1可知,对于 k F 0 , R R 1 n p F 0 , ρ m R 1 n p ,我们可得

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

又因为当 φ m = i = 1 2 ( u i k m ) + L 2 ( t 0 ρ m 2 , t 0 ; L 2 ( B ρ m , i ) ) 2 ,可得以下不等式

φ m = i = 1 2 t 0 ρ m 2 t 0 B ρ m , i | ( u i k m ) + | 2 d x d t ( k m + 1 k m ) 2 | A m ( k m + 1 ) | = k 2 2 2 m + 2 | A m ( k m + 1 ) | ,

其中 L 2 * ( B ρ m , i ) L 2 ( B ρ m , i ) ( i = 1 , 2 ) ,则 ξ m ( u i k m + 1 ) + L 2 ( B ρ m , i ) C | B ρ m , i [ u i ( , t ) > k m + 1 ] | 1 n ξ m ( u i k m + 1 ) + L 2 * ( B ρ m , i ) . ,故由定理2.2.1可得 ξ m ( u i k m + 1 ) + L 2 ( t 0 ρ m 2 , t 0 ; L 2 ( B ρ m , i ) ) 2 | Q ρ m , i [ u i > k m + 1 ] | 2 n ξ m ( u i k m + 1 ) + L 2 ( t 0 ρ m 2 , t 0 ; L 2 * ( B ρ m , i ) ) 2 , ,其中 φ 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 ) . 是抛物距离。下面区分两种情形来估计 k 1 | B R | [ ( t 0 R 2 t 0 B R , 1 u 1 2 d x d t ) 1 2 + ( t 0 R 2 t 0 B R , 2 u 2 2 d x d t ) 1 2 ] k 2 1 | B R | i = 1 2 u i L 2 ( t 0 R 2 , t 0 ; L 2 ( B R , i ) ) 2 1 | B R | i = 1 2 u i L 2 ( t 0 ρ m 2 , t 0 ; L 2 ( B ρ m , i ) ) 2 1 | B R | φ m ,

情形1 对于任意 φ m + 1 C 2 2 m ( 1 + 2 n ) φ m 1 2 p + 2 n R 2 ( 1 n p ) k 4 n 4 p . y m = φ m k 2 | B R | ,不妨设 y m + 1 C 2 2 m ( 1 + 2 n ) y m 1 2 p + 2 n .

如果 y 0 = φ 0 k 2 | B R | θ 0 = C 1 α 2 2 α 2 ( 1 + 2 n ) , ,令 α = 2 n 2 p ,故 i = 1 2 ( u i k ) + L 2 ( t 0 R 2 , t 0 ; L 2 ( B R , i ) ) 2 θ 0 k 2 | B R | , ,所以

lim m y m = 0

因此有

Q R 2 , i

如果 u i 2 k ,则

k 2 1 θ 0 | B R | i = 1 2 u i L 2 ( t 0 R 2 , t 0 ; L 2 ( B R , i ) ) 2

所以

k = F 0 , R R 1 n p + 1 θ 0 | B R | [ ( t 0 R 2 t 0 B R , 1 u 1 2 d x d t ) 1 2 + ( t 0 R 2 t 0 B R , 2 u 2 2 d x d t ) 1 2 ] .

情形2 对于任意 Q T u ,证明同情形1类似。

因此结论成立。

u W ( 0 , T ; V ) ,记 Q ρ , τ ( X 0 ) = B ρ ( x 0 ) × ( t 0 , t 0 + τ ] Q T max { u 1 L ( Q 1 ) , u 2 L ( Q 2 ) } M 整体被界面 δ ( 0 , 1 ] 划分成两个子域 0 < max { ess sup Q ρ , τ , 1 ( u 1 k ) ± , ess sup Q ρ , τ , 2 ( u 2 k ) ± } δ M , u 。设 u D G ( Q T ) = D G ( Q T ; M , λ 1 , m 1 , n , p , F 0 , ρ , τ , C * ) 是个常数,对于 D G ± ( Q T ) ,设

u W 1 , 1 ( B R )

A ( k ) = { x B R | u ( x ) > k }

引理2.5.2设问题 的系数满足假设2.1.2。 l > k ( l k ) | A ( l ) | β R n + 1 | B R A ( k ) | A ( k ) A ( l ) | u | d x , β 是问题的弱下解,并且 u D G + ( Q T ) 。如果

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

Q R ( X 0 ) = B R ( x 0 ) × ( t 0 R 2 , t 0 ] Q T , 0 < R 1 ,参考定义2.4.2,其中参数 μ max { ess sup Q R , 1 u 1 , ess sup Q R , 2 u 2 } 与定义2.4.1相同。

如果 θ ( 0 , 1 ) 是问题的弱上解,且

k < μ (2.38)

i = 1 2 | Q R , i [ u i > k ] | θ | Q R | ,

证明 证明类似定理2.4.1。

定理2.5.3设问题 的系数满足假设2.1.2,如果 δ M H : = μ k > ( M + F 0 ) R 1 n + 2 p , 是问题的弱解,令 max { ess sup Q R 2 , 1 u 1 , ess sup Q R 2 , 2 u 2 } μ H 2 , θ D G + ( Q T ) F 0 = i = 1 2 f i L 2 ( 0 , T ; L n p n + p ( Ω i ) ) < ,且满足条件 和 。对于某一常数 R 0 = R , R m = R 2 + R 2 m + 1 ,初值 k 0 = μ H = k , k m = μ H 2 H 2 m + 1 ( m = 0 , 1 , 2 , ) ,则对于任意 Q R m ( X 0 ) = B R m ( x 0 ) × ( t 0 R m 2 , t 0 ] ,存在 ξ m ( x , t ) C ( [ t 0 R m 2 , t 0 ] ; C 0 ( B R m ( x 0 ) ) ) Q R m , ξ m , k m 使得

Q ρ , τ , ξ , k

其中 i = 1 2 ξ m ( u i k m ) + L 2 ( t 0 R m 2 , t 0 ; L 2 * ( B R m , i ) ) 2 C [ 2 2 m R 2 i = 1 2 ( u i k m ) + L 2 ( t 0 R m 2 , t 0 ; L 2 ( B R m , i ) ) 2 + ( M + F 0 , R m ) 2 i = 1 2 | Q R m , i [ u i > k m ] | 1 2 p ] . 仅依赖于 A m = i = 1 2 ( Q R m , i [ u i > k m ] ) i = 1 2 ( u i k m ) + L 2 ( t 0 R m 2 , t 0 ; L 2 ( B R m , i ) ) 2 i = 1 2 t 0 R m 2 t 0 B R m , i [ u i ( , t ) > k m ] | μ k | 2 d x d t = H 2 | A m | . ,并且 i = 1 2 ξ m ( u i k m ) + L 2 ( t 0 R m 2 , t 0 ; L 2 * ( B R m , i ) ) 2 C ( k m + 1 k m ) 2 | A m + 1 | n 2 n = ( H 2 m + 2 ) 2 | A m + 1 | n 2 n ,

证明 取与定理2.5.1的证明中相同的记号 H 2 2 2 m + 4 | A m + 1 | n 2 n C ( 2 2 m R 2 H 2 | A m | + ( M + F 0 , R m ) 2 | A m | 1 2 p ) .

因为 H > ( M + F 0 ) R 1 n + 2 p > ( M + F 0 , R m ) R 1 n + 2 p | A m | | Q R | ,可得

| A m + 1 | C ( 2 4 m R 2 | A m | + 2 2 m R 2 ( 1 n + 2 p ) | A m | 1 2 p ) n n 2 C ( 2 4 m R 2 ( 1 n + 2 p ) | A m | 1 2 p ) n n 2 .

由定理2.4.4可知,成立

y m = | A m | | Q R | (2.39)

或者

R ( 0 , 1 ] (2.40)

同理,由于 y m + 1 C 2 4 m n n 2 R 2 ( 1 n + 2 p ) n n 2 + n + 2 | A m | ( 1 2 p ) n n 2 C 2 4 m n n 2 y m ( 1 2 p ) n n 2 . ,由定理2.4.5可知,成立

( 1 2 p ) n n 2 > 1 (2.41)

或者

θ ( 0 , 1 ) (2.42)

所以有以下四种情形:

情形1 若(2.39)和(2.41)成立,则

y 0 = i = 1 2 | Q R , i [ u i > k ] | | Q R | θ

由于

lim m y m = 0

所以可得

Q R 2 , i

情形2 若(2.40)和(2.42)成立,则

u i μ H 2

情形3 若(2.39)和(2.42)成立,则

u D G + ( Q T )

情形4 若(2.40)和(2.41)成立,则

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

结合以上四种情形,对于任意 Q ^ 2 R ( X 0 ) = B 2 R ( x 0 ) × ( t 0 R 2 , t 0 ] Q T

0 < R 1 2

其中 μ max { ess sup Q ^ 2 R , 1 u 1 , ess sup Q ^ 2 R , 2 u 2 } 。由引理2.5.1可得,有

0 < μ k δ M

因此结论成立。

推论2.5.1设问题(1.1)的系数满足假设2.1.2。如果 0 < σ < 1 是问题的弱解, u ,记 | B R , i [ u i ( , t 0 R 2 ) > k ] | ( 1 σ ) | B R , i | , Q ^ 2 R , i = Q ^ 2 R Q i ,令 B R , i = B R Ω i ( i = 1 , 2 ) ,则在定理2.5.3的条件下,对于任意 s = s ( σ ) 1 ,存在 H : = μ k 2 s ( M + F 0 ) R 1 n + 2 p , max { ess sup Q R 2 , 1 u 1 , ess sup Q R 2 , 2 u 2 } μ H 2 s , 使得

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

其中 n , λ 0 , p , C * , σ 仅依赖于 F 0 = i = 1 2 f i L 2 ( 0 , T ; L n p n + p ( Ω i ) ) < u D G + ( Q T ) ,并且 X 0 = ( x 0 , t 0 ) Γ × ( 0 , T ]

2.6. 弱解的全局Hölder连续性

因为 Q ^ 2 R a ( X 0 ) = B 2 R ( x 0 ) × ( t 0 , t 0 + a R 2 ] Q T 被界面 0 < R 1 2 分成了两个区域 0 < a 1 μ max { ess sup Q ^ 2 R , 1 a u 1 , ess sup Q ^ 2 R , 2 a u 2 } ,并且陈亚浙 [14] 和GaryMLieberman [15] 都已经证明了 0 < μ k δ M 在单一区域上的内部Hölder连续性和边界Hölder连续性,同时在上一节中已经给出了 0 < σ < 1 在界面上 u 和界面与初值层的交界处附近 | B R , i [ u i ( , t ) > k ] | ( 1 σ ) | B R , i | , t ( t 0 , t 0 + a R 2 ] , 的Hölder连续性。综合以上情况可得

Q ^ 2 R , i a = Q ^ 2 R a Q i

因此由 [15] 可知有以下结论成立。

定理2.6.1设问题(1.1)的系数满足假设2.1.2。对于 B R , i = B R Ω i ( i = 1 , 2 ) ,令 s H : = μ k 2 s ( M + F 0 ) R 1 n + 2 p , ,假设扩散系数矩阵分量 i = 1 2 | Q R , i a [ u i > μ H 2 s ] | C σ a s | Q R a | , ,并且速度场 Q R a = B R ( x 0 ) × ( t 0 , t 0 + a R 2 ] 的分量 Q R , i a = Q R a Q i 。同时对于某个非负常数 n , p , λ 1 , C * ,有下式成立

i = 1 , 2

如果 A R , i ( k , t ) = B R , i [ u i ( , t ) > k ] 是问题(1.1)的弱解,并且 A R , i ( k ) = Q R , i a [ u i > k ] A R , i ( k ) = t 0 t 0 + a R 2 A R , i ( k , t ) d t ,则

k 0 = μ H = k , k l = μ H 2 l ( l = 0 , 1 , 2 , )

并且

A R , i ( k l + 1 , t ) A R , i ( k l , t ) A R , i ( k , t )

3. Henry界面模型

对于Henry界面问题(1.2),由于界面上的Henry条件,故在弱形式中会出现界面上的积分项,使得在证明相关定理时有一定困难。为了避免界面积分项的出现,下面通过函数变换将问题(1.2)转换成另一种形式。令 l 0 , t 0 t t 0 + a R 2 ,故有 ( k l + 1 k l ) | A R , i ( k l + 1 , t ) | β R n + 1 | B R , i A R , i ( k l , t ) | A R , i ( k l , t ) A R , i ( k l + 1 , t ) | u i | d x . ,此时记Henry界面问题的解为 | A R , i ( k , t ) | ( 1 σ ) | B R , i | ,其中 | B R , i A R , i ( k l , t ) | σ | B R , i | ,故有

| A R , i ( k l + 1 , t ) | C R n + 1 2 l σ | B R , i | H | A R , i ( k l , t ) A R , i ( k l + 1 , t ) | 1 2 ( A R , i ( k l , t ) | u i | 2 d x ) 1 2 . (3.1)

其中 t ( t 0 , t 0 + a R 2 ] 满足界面条件,对于某个 | A R , i ( k l + 1 ) | C R 2 l σ H | A R , i ( k l ) A R , i ( k l + 1 ) | 1 2 ( t 0 t 0 + a R 2 B R , i | ( u i k ) + | 2 d x d t ) 1 2 . ,边界 A R ( k , t ) = i = 1 2 A R , i ( k , t ) 和界面 A R ( k ) = i = 1 2 A R , i ( k ) | A R ( k l + 1 ) | C R 2 l σ H | A R ( k l ) A R ( k l + 1 ) | 1 2 i = 1 2 ( u i k ) + L 2 ( t 0 , t 0 + a R 2 ; L 2 ( B R , i ) ) 的。

3.1. 弱解的存在唯一性

在接下来的讨论中,令问题 的系数满足假设2.1.2。令 ξ ( x ) ,在H上定义标量乘积为

B 2 R ( x 0 )

相应的范数记为 B R ( x 0 ) 。令 ξ ( x ) = 1 ,定义 u D G + ( Q T ) ,所以 k = k l Q ρ , τ = Q ^ 2 R a 等价,并且

i = 1 2 ( u i k l ) + L 2 ( t 0 , t 0 + a R 2 ; L 2 ( B R , i ) ) 2 i = 1 2 ξ ( u i k l ) + L 2 ( t 0 , t 0 + a R 2 ; L 2 ( B 2 R , i ) ) 2 C ( H 2 4 l | B R | + H 2 4 l R 2 | Q R a | + ( M + F ^ 0 , 2 R a ) 2 | Q R a | 1 2 p ) ,

所以 的变分问题为:给定 F ^ 0 , 2 R a = i = 1 2 f i L 2 ( t 0 , t 0 + a R 2 ; L n p n + p ( B 2 R , i ) ) < ,设 H > 2 s ( M + F 0 ) R 1 n + 2 p > 2 s ( M + F ^ 0 , 2 R a ) R 1 n + 2 p i = 1 2 ( u i k l ) + L 2 ( t 0 , t 0 + a R 2 ; L 2 ( B R , i ) ) 2 C H 2 4 l | B R | ,寻找弱下解(弱上解)。

| A R ( k l + 1 ) | 2 C R n + 2 σ 2 | A R ( k l ) A R ( k l + 1 ) | , 使得

0 (3.2)

如果 s 1 既是弱下解,又是弱上解,则称 s | A R ( k s ) | 2 C R n + 2 σ 2 ( | A R ( k 0 ) | | A R ( k s ) | ) C R n + 2 σ 2 | A R ( k 0 ) | C R n + 2 σ 2 | Q R a | C a σ 2 | Q R a | 2 . 是弱解。其中

| A R ( k s ) | = i = 1 2 | Q R , i a [ u i > μ H 2 s ] | C σ a s | Q R a | .

u D G + ( Q T )

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

事实上,根据假设2.1.2可得 Q ^ 2 R ( X 0 ) = B 2 R ( x 0 ) × ( t 0 , t 0 + R 2 ] Q T 0 < R 1 2 上是连续的和强制的,所以由 [10] 可知,线性问题 存在唯一弱解。

3.2. 弱解的全局Hölder连续性

同样对于问题 也可以定义De Giorgi类。

定义3.2.1称定义于 μ max { ess sup Q ^ 2 R , 1 u 1 , ess sup Q ^ 2 R , 2 u 2 } 上的函数 0 < μ k δ M 属于De Giorgi类,如果 0 < σ < 1 u ,且对于 | B R , i [ u i ( , t 0 ) > k ] | ( 1 σ ) | B R , i | , Q ^ 2 R , i = Q ^ 2 R Q i B R , i = B R Ω i ( i = 1 , 2 ) 满足 s 0 = s 0 ( σ ) 1 ,并且 H : = μ k 2 s 0 ( M + F 0 ) R 1 n + 2 p , ,有下式成立:

sup t 0 < t t 0 + R 2 i = 1 2 | B R , i [ u i ( , t ) > μ H 2 s 0 ] | ( 1 σ + 1 2 min { σ , 1 σ } ) | B R | , (3.3)

其中 s 0 n , p , λ 0 , C * , σ ,常数 F 0 = i = 1 2 f i L 2 ( 0 , T ; L n p n + p ( Ω i ) ) < ξ ( x ) B R ( x 0 ) B β R ( x 0 ) 依赖于 ξ ( x ) = 1 ,记De Giorgi类为 0 < β < 1 。如果 0 < a 1 ,且满足 +,则记 Q R a = B R ( x 0 ) × ( t 0 , t 0 + a R 2 ] ;如果 A R a ( k ) = i = 1 2 ( Q R , i a [ u i > k ] ) ,且满足 -,则记 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 ] , 。显然 F 0 , R a = i = 1 2 f i L 2 ( t 0 , t 0 + a R 2 ; L n p n + p ( B R , i ) ) <

定理3.2.1设问题 的系数满足假设2.1.2。如果 l 1 > 1 是问题的弱下解,且对于某常数 i = 1 2 ξ ( u i k ) + ( , t ) L 2 ( B R , i ) 2 ( 1 1 2 l 1 ) 2 H 2 i = 1 2 | B β R , i [ u i ( , t ) > μ H 2 l 1 ] | . sup t 0 < t t 0 + a R 2 i = 1 2 | B β R , i [ u i ( , t ) > μ H 2 l 1 ] | | B R | ( 1 σ ( 1 2 l 1 ) 2 + 4 ε + C γ ( ε ) ( 1 β ) 2 ( | A R a ( k ) | | Q R | ) 1 2 p ) , ,则 Q R ( X 0 ) = B R ( x 0 ) × ( t 0 , t 0 + R 2 ] ;如果 t 0 < t t 0 + a R 2 是问题的弱上解,则 i = 1 2 | B R , i [ u i ( , t ) > μ H 2 l 1 ] | | B R | ( 1 β n + 1 σ ( 1 2 l 1 ) 2 + 4 ε + C γ ( ε ) ( 1 β ) 2 ( | A R a ( k ) | | Q R | ) 1 2 p ) . 。其中 β ( 0 , 1 ) 依赖于 1 β = ( | A R a ( k ) | | Q R | ) 1 3 ( 1 2 p ) ,并且 i = 1 2 | B R , i [ u i ( , t ) > μ H 2 l 1 ] | | B R | ( 1 σ ( 1 2 l 1 ) 2 + 4 ε + C γ ( ε ) ( | A R a ( k ) | | Q R | ) 1 3 ( 1 2 p ) ) .

证明 证明类似定理2.4.1。

注3.2.1如果 q = ε γ ( ε ) ,则可得

ε ( 0 , 1 ] (3.4)

在此基础之上,类似非完美界面模型中的第2.4节,可以得到对应于Henry界面模型的类似定理。进一步,参考第2.5节中的证明,可得Henry界面模型的弱解的全局Hölder连续性,即存在 ε = ψ ( q ) ,有

q 0

3.3. 梯度的Lq估计

对于Henry界面问题(3.1),首先将界面局部拉直,同时假设 ε 0 。对于任意 ψ ( q ) 0 ,存在 ε = ψ ( ( | A R a ( k ) | | Q R | ) 1 3 ( 1 2 p ) ) i = 1 2 | B R , i [ u i ( , t ) > μ H 2 l 1 ] | | B R | [ 1 σ ( 1 2 l 1 ) 2 + C ψ ( ( | A R a ( k ) | | Q R | ) 1 3 ( 1 2 p ) ) ] . 映射 C ψ ( ( | A R a ( k ) | | Q R | ) 1 3 ( 1 2 p ) ) C ψ ( a 1 3 ( 1 2 p ) ) ,使得

a = a ( σ ) > 0

C ψ ( a 1 3 ( 1 2 p ) ) 1 8 min { 1 σ , σ }

a 1

其中 t j = t 0 + j a R 2 ( j = 1 , 2 , , N ) 表示以 s 1 < s 2 < < s N 为球心,r为半径的球。定义

sup t j 1 < t t j i = 1 2 | B R , i [ u i ( , t ) > μ H 2 s j ] | ( 1 σ + j 4 N min { σ , 1 σ } ) | B R | ,

s 1 , s 2 , , s N

反之定义

n , p , λ 1 , γ ( ) , σ

j = 1

所以 i = 1 2 | B R , i [ u i ( , t ) > μ H 2 l 1 ] | | B R | ( 1 σ ( 1 2 l 1 ) 2 + 1 8 min { 1 σ , σ } ) , t ( t 0 , t 0 + a R 2 ] = ( t 0 , t 1 ] . ,并且映射 l 1 = l 1 ( σ ) sup t 0 < t t 1 i = 1 2 | B R , i [ u i ( , t ) > μ H 2 l 1 ] | | B R | ( 1 σ + 1 4 min { 1 σ , σ } ) , 展平。令 s 0 p 1 > l 1 ,所以 | A R a ( μ H 2 p 1 ) | = i = 1 2 | Q R , i a [ u i > μ H 2 p 1 ] | C σ a p 1 | Q R a | = C a σ p 1 l 1 | Q R | , ,且 Q R a ( X 0 ) = B R ( x 0 ) × ( t 0 , t 1 ] 。不妨记 n , λ 0 , p , γ ( ) ,并且选择适当的 H 2 p 1 , s 1 p 1 ,使得 H , l 1 ,并记

sup t 0 < t t 1 i = 1 2 | B R , i [ u i ( , t ) > μ H 2 s 1 ] | | B R | [ 1 σ ( 1 2 ( s 1 p 1 ) ) 2 + C ψ ( ( C a σ p 1 l 1 ) 1 3 ( 1 2 p ) ) ] .

p 1

1 8 N min { 1 σ , σ }

s 1 > p 1 1 σ + 1 8 N min { 1 σ , σ } ,其中 s 1 = s 1 ( σ ) 。定义

sup t 0 < t t 1 i = 1 2 | B R , i [ u i ( , t ) > μ H 2 s 1 ] | | B R | ( 1 σ + 1 4 N min { 1 σ , σ } ) .

其中 Q R a ( X 0 ) = B R ( x 0 ) × ( t j 1 , t j ] ( j = 2 , , N ) 。记 j = 2 , , N s 0 s > s 0 sup t 0 R 2 < t t 0 i = 1 2 | B R , i [ u i ( , t ) > μ H 2 s 0 ] | ( 1 σ + 1 2 min { σ , 1 σ } ) | B R | , i = 1 2 | Q R , i [ u i > μ H 2 s ] | C σ s | Q R | C σ s s 0 | Q R | . 为部分拉直的界面上的单位外法向量,由 C σ s s 0 θ ( 0 , 1 ) 指向 max { ess sup Q R 2 , 1 u 1 , ess sup Q R 2 , 2 u 2 } μ H 2 s . 。所以

D G ( Q T )

u D G ( Q T ) X 0 = ( x 0 , t 0 ) Γ × ( 0 , T ] ,故 可转换为

Q ^ 2 R ( X 0 ) = B 2 R ( x 0 ) × ( t 0 R 2 , t 0 ] Q T (3.5)

并且 0 < R 1 2 满足假设2.1.2的(1),(2)中的(i),不妨令其中的常数仍记为 μ ˜ min { ess inf Q ^ 2 R , 1 u 1 , ess inf Q ^ 2 R , 2 u 2 } 。给定 0 < k μ ˜ δ M ,设 0 < σ < 1 u 是截断函数,并且满足 | B R , i [ u i ( , t 0 R 2 ) < k ] | ( 1 σ ) | B R , i | , ,在 Q ^ 2 R , i = Q ^ 2 R Q i B R , i = B R Ω i ( i = 1 , 2 ) s = s ( σ ) 1 ;在 H : = k μ ˜ 2 s ( M + F 0 ) R 1 n + 2 p , 上有 min { ess inf Q R 2 , 1 u 1 , ess inf Q R 2 , 2 u 2 } μ ˜ + H 2 s ,在 u 上有 Γ × ( 0 , T ] ( Γ × [ 0 , T ] ) Ω ω ( R ) 。其中C与r无关。对于任意 ( 0 , R 0 ] ,定义

ω ( ν R ) η ω ( R ) + K R α , R ( 0 , R 0 ] ,

则有

0 < ν , η < 1

对于问题(3.5)的弱解 0 < α 1 ,函数 K 0 有弱导数 0 < β α

定理3.3.1 (Coccippoli型不等式)对于问题 ,设常数 C 1 ω ( R ) C ( R R 0 ) β ( ω ( R 0 ) + K R 0 α ) , R ( 0 , R 0 ] , β , C ,则弱解 ν , η , α 满足

u D G ( Q T )

其中C依赖于 X 0 = ( x 0 , t 0 ) Γ × ( 0 , T ] Q R ( X 0 ) = B R ( x 0 ) × ( t 0 R 2 , t 0 ] Q T 0 < R 0 1 有关, 0 < R R 0

证明 对于问题 ,取测试函数为 max { osc Q R , 1 u 1 , osc Q R , 2 u 2 } C ( R R 0 ) α ( osc Q R 0 , 1 u 1 + osc Q R 0 , 2 u 2 + ( M + F 0 ) R 0 1 n + 2 p ) , ,以下简记截断函数为 α ( 0 , 1 n + 2 p ) ,类似定理2.4.1的证明,有

C 1

由于

D G ( Q T )

Q R , i = Q R Q i

osc Q R , i u i = ess sup Q R , i u i ess inf Q R , i u i

( i = 1 , 2 )

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

μ ( R ) = max { ess sup Q R , 1 u 1 , ess sup Q R , 2 u 2 } , μ ˜ ( R ) = min { ess inf Q R , 1 u 1 , ess inf Q R , 2 u 2 } , (3.6)

(3.6)两边对 max { osc Q R , 1 u 1 , osc Q R , 2 u 2 } ω ( R ) = μ ( R ) μ ˜ ( R ) osc Q R , 1 u 1 + osc Q R , 2 u 2 . 积分,故有

0 < μ ( R ) ( μ ( R ) 1 2 ω ( R ) ) < δ M (3.7)

下面记 0 < ( μ ˜ ( R ) + 1 2 ω ( R ) ) μ ˜ ( R ) < δ M ( i = 1 , 2 ) 上的截断函数,故有

| B R 2 , i [ u i ( , t 0 ( R 2 ) 2 ) > μ ( R ) 1 2 ω ( R ) ] | 1 2 | B R 2 , i | , (3.8)

对于函数 | B R 2 , i [ u i ( , t 0 ( R 2 ) 2 ) < μ ˜ ( R ) + 1 2 ω ( R ) ] | 1 2 | B R 2 , i | . s = s ( 1 2 ) 1 R ( 0 , R 0 ] 时取得最小值,则

1 2 ω ( R ) = H 2 s ( M + F 0 ) R 1 n + 2 p ,

所以

μ ( R 4 ) = max { ess sup Q R 4 , 1 u 1 , ess sup Q R 4 , 2 u 2 } μ ( R ) H 2 s = μ ( R ) ω ( R ) 2 s + 1 ,

因此结合(3.7)可得结论。

对于任意 2 s ( M + F 0 ) R 0 1 n + 2 p = δ M ,可得

R ( 0 , R 0 ] (3.9)

其中C与r无关。由Sobolev嵌入可得:

i) 当 ω ( R 4 ) ω ( R ) ( 1 1 2 s + 1 ) 时,

s = s ( 1 2 ) 1

ii) 当 R ( 0 , R 0 ] 时,

1 2 ω ( R ) = H 2 s ( M + F 0 ) R 1 n + 2 p ,

结合(3.9)可得Poincare不等式:

μ ˜ ( R 4 ) = min { ess inf Q R 4 , 1 u 1 , ess inf Q R 4 , 2 u 2 } μ ˜ ( R ) + H 2 s = μ ˜ ( R ) + ω ( R ) 2 s + 1 , (3.10)

其中当 R ( 0 , R 0 ] 时, ω ( R 4 ) ω ( R ) ( 1 1 2 s + 1 ) ;当 δ M H 时, R R 0

所以有以下结论:

定理3.3.2对于问题(3.5),设常数 2 s ( M + F 0 ) R 1 n + 2 p δ M H ω ( R ) 2 δ M ( R R 0 ) 1 n + 2 p ω ( R 4 ) 2 M ( R R 0 ) 1 n + 2 p 2 s + 1 δ ( M + F 0 ) R 1 n + 2 p ,则弱解 0 < R 1 满足

ω ( R 4 ) ω ( R ) ( 1 1 2 s + 1 ) + C 2 s ( M + F 0 ) R 1 n + 2 p ,

其中C依赖于 C > 0 δ ω ( R ) C ( R R 0 ) α ( ω ( R 0 ) + ( M + F 0 ) R 0 1 n + 2 p ) , 有关, α ( 0 , 1 n + 2 p )

引理3.3.1 [16] [17] 令 u D G ( Q T ) 是开集, Γ × ( 0 , T ] Q Q T ,几乎在Q中 α ( 0 , 1 n + 2 p ) 。假设对于 C 1 ,有

i = 1 , 2 (3.11)

其中 [ u i ] C α ( Q ¯ i ) C d α ( M + F 0 d 1 n + 2 p ) , d = min { 1 , d i s t { Q , p Q T } } 是固定常数。则存在 Q i = Q Q i 使得

α , C

并且

n , λ 1 , p , γ ( ) , δ

其中c依赖于 F 0 = i = 1 2 f i L 2 ( 0 , T ; L n p n + p ( Ω i ) ) < ,并且当 X 0 = ( x 0 , t 0 ) ( Γ × ( 0 , T ] ) Q 时, Q R ( X 0 ) = B R ( x 0 ) × ( t 0 R 2 , t 0 ]

注3.3.1在(3.11)的右边用 Q R , i = Q R Q i ( i = 1 , 2 ) 代替 R ( 0 , d ] ,引理3.3.1的结论仍成立。

定理3.3.3 (梯度的 max { osc Q R , 1 u 1 , osc Q R , 2 u 2 } C ( R d ) α ( M + F 0 d 1 n + 2 p ) . 估计)对于问题(3.5),设常数 R d 1 R n + 2 + α Q R , i Q i | u i ( X ) u ˜ i | d X C ( n ) R α osc Q R , 1 u 1 C d α ( M + F 0 d 1 n + 2 p ) , u ˜ i = 1 | Q R , i Q i | Q R , i Q i u i ( X ) d X ,则存在 u i L l o c 1 , n + 2 + α ( Q i ; δ ˜ ) 使得

u i C l o c α ( Q ¯ i ; δ ˜ )

并且对于任意 δ ˜ ,有

[ u 1 ] C α ( Q ¯ 1 )

其中C与r无关。

证明 应用定理3.3.2 (2r代替r)以及Hölder不等式可得

[ u 2 ] C α ( Q ¯ 2 )

应用(3.10) (取 X 0 = ( x 0 , t 0 ) ( Γ × ( 0 , T ] ) Q ),则

Y 1 = ( y 1 , t Y 1 ) Q ¯ 1

应用(3.10) (取 t 0 > t Y 1 ),则

δ ˜ ( X 0 , Y 1 ) d

再一次应用定理3.3.2 (2r代替r)以及Hölder不等式可得

R 1 = δ ˜ ( X 0 , Y 1 )

所以

Y 1 Q R 1 , 1 ( X 0 )

由定理3.3.1可得

| u 1 ( X 0 ) u 1 ( Y 1 ) | osc Q R 1 , 1 u 1 C ( R 1 d ) α ( M + F 0 d 1 n + 2 p ) .

其中 [ u 1 ] C α ( Q ¯ 1 ) C d α ( M + F 0 d 1 n + 2 p ) . δ ˜ ( X 0 , Y 1 ) > d 有关。上式两边同除以 | u 1 ( X 0 ) u 1 ( Y 1 ) | 2 M ( δ ˜ ( X 0 , Y 1 ) d ) α , ,并且 [ u 1 ] C α ( Q ¯ 1 ) 2 M d α . ,故有

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

适当选取 Y 2 = ( y 2 , t Y 2 ) Q ¯ 2 使得 X 0 = ( x 0 , 0 ) Γ 。因此对于 Q R ( X 0 ) = B R ( x 0 ) × ( R 2 , R 2 ) 可得

Q R ( X 0 ) (3.12)

定义

Γ × ( 0 , T ]

Q R , 1

所以

Q R , 2

同时(3.12)可重新表述为

υ

其中C与 i = 1 , 2 有关。由定理3.3.1可知,存在 u i , υ ( + ) = { max { u i , υ } ( x , t ) Q R , i Q i ; υ ( x , t ) Q R , i Q i , 使得 u i , υ ( ) = { min { u i , υ } ( x , t ) Q R , i Q i ; υ ( x , t ) Q R , i Q i . ,并且

X 0 = ( x 0 , 0 ) Γ

Q R ( X 0 ) = B R ( x 0 ) × ( R 2 , R 2 ) ,则 u W ( 0 , T ; V ) ,并且

max { u 1 L ( Q 1 ) , u 2 L ( Q 2 ) } M <

本节讨论的是界面附近梯度的 k : = max { ess sup Q R , 1 Ω 1 u 1 , ess sup Q R , 2 Ω 2 u 2 } , δ M max { ess sup Q R , 1 Q 1 u 1 k , ess sup Q R , 2 Q 2 u 2 k } > 0 , 估计,对于 u k ( + ) = ( u 1 , k ( + ) , u 2 , k ( + ) ) D G + ( Q T ) = D G + ( Q T ; M , λ 1 , m 1 , p , n , F 0 , ρ , τ , C * , δ ) ,采用同样的证明方法,均可得到单一区域上梯度的 λ 1 , p , m 1 , p , n , F 0 , ρ , τ , C * 估计。

3.4. 梯度的内部Hölder连续性

在本节中,对于Henry界面问题(3.1),不考虑对流项,并且其弱形式中没有界面积分项。设问题(3.1)的系数满足假设2.1.2。定义 u W ( 0 , T ; V )

k ˜ : = min { ess inf Q R , 1 Ω 1 u 1 , ess inf Q R , 2 Ω 2 u 2 } , δ M max { k ˜ ess inf Q R , 1 Q 1 u 1 , k ˜ ess inf Q R , 2 Q 2 u 2 } > 0 , 关于t是 u k ˜ ( ) = ( u 1 , k ˜ ( ) , u 2 , k ˜ ( ) ) D G ( Q T ) = D G ( Q T ; M , λ 1 , m 1 , p , n , F 0 , C * , δ ) 连续,关于x是 u W ( 0 , T ; V ) 阶Hölder连续; X 0 = ( x 0 , 0 ) Γ ,

Q R ( X 0 ) = B R ( x 0 ) × ( 0 , R 2 )

进一步, 0 < R 1 还满足以下假设:

假设3.4.1 max { u 1 L ( Q 1 ) , u 2 L ( Q 2 ) } M < ,即存在常数 α 1 ( 0 , 1 ] [ u 0 i ] C α 1 ( Ω ¯ i ) < ( i = 1 , 2 ) 使得

0 < R R 0 1

并且,对于任意整数 0 < α min { α 1 , 1 n + 2 p } ,存在 C 1 (依赖于l),使得

max { osc Q R , 1 u 1 , osc Q R , 2 u 2 } C ( R R 0 ) α ( osc Q R 0 , 1 u 1 + osc Q R 0 , 2 u 2 + R 0 α ( 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

参考文献 [18] [19] ,可得到梯度的内部Hölder连续性。

定理3.4.1设问题(3.1)的系数满足假设2.1.2和假设3.4.1。 F 0 = i = 1 2 f i L 2 ( 0 , T ; L n p n + p ( Ω i ) ) < μ ( R ) , μ ˜ ( R ) , ω ( R ) 是问题的弱解,则对于 u k ( + ) D G + ( Q T ) k ess sup Q R , i Ω i u i ,有

| B R 2 , i [ u i ( , 0 ) > k ] | = 0.

其中C依赖于 μ k 2 s ( M + F 0 ) R 1 n + 2 p , max { ess sup Q R 4 , 1 u 1 , k ( + ) , ess sup Q R 4 , 2 u 2 , k ( + ) } μ H 2 s μ μ k 2 s . u k ˜ ( = ) D G ( Q T ) 范数 k ˜ μ ˜ 2 s ( M + F 0 ) R 1 n + 2 p , 。特别地,

min { ess inf Q R 4 , 1 u 1 , k ˜ ( ) , ess inf Q R 4 , 2 u 2 , k ˜ ( ) } μ ˜ + k ˜ μ ˜ 2 s .

4. 结论

本文分别考虑了耦合非完美界面条件和Henry界面条件的两相流模型。对于界面模型弱解的相关性质利用De Giorgi迭代法给出了详细证明,例如极值原理,局部极值原理等。在此基础之上,给出弱解及梯度的Hölder连续性。对于Henry界面模型,我们也给出了梯度的 ω ( R 4 ) ω ( R ) = μ ( R ) μ ˜ ( R ) 2 s + 1 ( M + F 0 ) R 1 n + 2 p + k k ˜ . 估计(存在 k k ˜ = max { ess sup Q R , 1 Ω 1 u 1 , ess sup Q R , 2 Ω 2 u 2 } min { ess inf Q R , 1 Ω 1 u 1 , ess inf Q R , 2 Ω 2 u 2 } ess sup Q R , 1 Ω 1 u 1 ess inf Q R , 1 Ω 1 u 1 + ess sup Q R , 2 Ω 2 u 2 ess inf Q R , 2 Ω 2 u 2 R α 1 [ u 0 1 ] C α 1 ( Ω ¯ 1 ) + R α 1 [ u 0 2 ] C α 1 ( Ω ¯ 2 ) , )的详细证明。

致谢

感谢导师袁光伟研究员对本文的建议和帮助。

文章引用

史雪婷,杨田洁. 一类抛物型界面问题的正则性分析
Regular Analysis of a Class of Parabolic Interface Problems[J]. 应用数学进展, 2023, 12(02): 645-678. https://doi.org/10.12677/AAM.2023.122068

参考文献

  1. 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. 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. 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. 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. 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. 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. 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. 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. 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. 10. 王耀东. 偏微分方程的 理论[M]. 北京: 北京大学出版社, 1989: 141-202.

  11. 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. 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. 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. 14. 陈亚浙. 二阶抛物型偏微分方程[M]. 北京: 北京大学出版社, 2003: 117-140.

  15. 15. Lieberman, G.M. (1996) Second Order Parabolic Differential Equations. World Scientific Publishing Company, Singapore. https://doi.org/10.1142/3302

  16. 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. 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. 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. 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

期刊菜单