Advances in Applied Mathematics
Vol. 09  No. 04 ( 2020 ), Article ID: 35230 , 9 pages
10.12677/AAM.2020.94066

Properties of Solutions for a Class of Nonlinear Pseudoparabolic Equations

Jiazhuo Cheng1, Lingyu Jin1, Shaomei Fang1,2

1College of Mathematics and Informatics, South China Agricultural University, Guangzhou Guangdong

2School of Mathematics and Statistics, Yangtze Normal University, Chongqing

Received: Apr. 5th, 2020; accepted: Apr. 17th, 2020; published: Apr. 24th, 2020

ABSTRACT

This paper studies the weak solutions of a class of fractional nonlinear pseudo-parabolic equations. First, a priori estimates of the solution of the equation are obtained through the energy method discussion, and then the Galerkin method is used to construct an approximate solution sequence to prove the existence and uniqueness of the weak solution of the equation.

Keywords:A Priori Estimates, Galerkin Method, Weak Solution

一类非线性伪抛物型方程解的性质

程嘉卓1,金玲玉1,房少梅1,2

1华南农业大学,数学与信息学院,广东 广州

2长江师范大学,数学与统计学院,重庆

收稿日期:2020年4月5日;录用日期:2020年4月17日;发布日期:2020年4月24日

摘 要

本文研究了一类分数阶非线性伪抛物方程的弱解问题。首先通过能量方法讨论得到方程解的先验估计,然后利用Galerkin方法构造近似解序列来证明方程弱解的存在唯一性。

关键词 :先验估计,Galerkin方法,弱解

Copyright © 2020 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. 引言

我们考虑分数阶非线性伪抛物方程的以下初始边值问题

u t + ( Δ ) α u t + ( Δ ) α u μ ( u ) Δ u + ( 1 + σ ( u ) ) u = f ( x , t ) , (1)

u ( x , t ) | x Ω = 0 , (2)

u ( x , 0 ) = u ˜ 0 ( x ) , (3)

其中 t > 0 , x Ω , Ω R 是有界光滑域,T和 α 是正常数,是满足后面指定条件的给定函数。

伪抛物线方程

u t u x x t = F ( x , t , u x , u x x ) , (4)

由Coleman等人于20世纪60年代发现,用来模拟不稳定的简单剪切流的特殊运动状态 [1]。伪抛物线方程用于各种领域,例如均匀流体通过裂隙岩石的渗漏 [2] (三阶项的系数表示岩石的裂隙程度,其减小程度对应于增加裂纹的程度),非线性色散长波的单向传播 [3] [4] (其中u是振幅或旋度),种族迁移的描述 [5] (其中u是人口密度)。由于伪抛物线方程的广泛应用,它们引起了数学家的极大关注,例如讨论了解的存在性、渐近行为、正则性和衰减性 [6] [7] [8] [9]。

描述小振幅长波在非线性色散介质中的传播过程时,经典的伪抛物线方程通常必须考虑耗散机制,以便准确反映实际情况。但是,引起波衰减的机制非常复杂,人们对其了解不多。在这种情况下,人们可能被迫依赖耗散的临时模型 [10]。在对平面波的单向传播进行建模时,需要在模型中增加非线性和色散的耗散,如下式

u t + u x + u u x v u x x u x x t = 0 , (5)

其中 v > 0 是固定常数。方程(5)有时被称为Benjamin-Bona-Mahony-Burgers (BBMB)方程,它是Benjamin-Bona-Mahony (BBM)方程 [3] 和Burgers类型耗散项 [8] 的结合。最近,许多科学家开始注意到出现在工程学和物理学中的三级流体流动数学模型框架中的BBMB方程 [11] [12] [13]。T.Aziz等人研究了以下的BBMB运动方程

ρ u t = μ u y y + α 1 u y y t + 6 β 3 ( u y y ) 2 u y y φ k ( μ u + α 1 u t + 2 β 3 ( u y ) 2 u ) , (6)

其中 y > 0 , t > 0 ,u表示速度分量, μ 表示粘度系数, α 1 β 3 是材料常数 [11]。方程(6)适用于带有多孔介质的刚性板上非定常流动三级流体的数学模型。后来,在 [14] (可以看作 [11] 的延续)中,作者在如下改进的非定常流动三级流体方程中证明了解的存在性,唯一性和指数衰减

u t + α ( Δ ) u t + ( Δ ) u μ ( u ) Δ u + ( γ + β σ ( u ) ) u = f ( x , t ) , (7)

其中 α = β = γ = 1 , μ C ( R ; R ) , σ C 1 ( R ; R ) y ( 0 y z σ ( z ) d z ) y 2 σ ( y ) ( y R )

分数耗散算子 ( Δ ) α 可以看作是Levy稳定扩散过程的无穷小生成器。与微积分方程相比,它可以更准确地描述某些物理现象 [15] [16] [17]。因此,越来越多的科学家致力于分数阶微分方程的研究 [16] [17] [18]。

基于以上工作,我们研究了分数阶方程(1)-(3)的弱解的存在性和唯一性。

2. 预备知识和引理

现对文中将用到的定义、引理和记号做一些说明。记函数空间 L 2 ( Ω ) 的范数通常表示为 ,其标量积由 , 表示, L p ( Ω ) 的范数记为 L p ( Ω ) ,C表示一个正常数,该常数可能在各行之间变化。

定义1:范数 L p ( 0 , T ; X ) 表示实函数 f : [ 0 , T ] X 的Banach空间,使得

f L p ( 0 , T ; X ) = ( 0 T f X p d t ) 1 p < , 1 p < ,

f L ( 0 , T ; X ) = esssup 0 t T f X , p = .

引理1:令 σ μ 均是 R R 上映射并且属于C,那么

μ ( v ) L ( Ω ) < C , (8)

σ ( v ) L ( Ω ) < C , (9)

其中 v H α + 1 ( Ω ) 且常数C依赖于v。

证明:利用Gagliardo-Nirenberg不等式,有

v L ( Ω ) c 2 v H α ( Ω ) η v 1 η < C . (10)

其中 η = n 2 α + 1 α 。应用这个不等式,易知

μ ( v ) L ( Ω ) < C , (11)

σ ( v ) L ( Ω ) < C . (12)

引理2:对于 α 1 2 ,令初值满足 u ˜ 0 ( x ) H α ( Ω ) ,且

μ C ( R ; R ) , μ ( 0 ) = 0 , μ ( z ) > 0 ( z R , z 0 ) , (13)

σ C 1 ( R ; R ) , σ ( 0 ) = 0 , σ ( z ) > 0 ( z R , z 0 ) , (14)

f L 2 ( Ω ) , (15)

那么问题(1)-(3)有一个弱解u,满足

u ( t ) H α ( Ω ) 2 C . (16)

证明:将(1)乘以 u ( t ) ,并关于空间变量x在 Ω 上进行积分,我们很容易就得出了等式

1 2 d d t u ( t ) H α ( Ω ) 2 + u ( t ) H α ( Ω ) 2 + I 1 + I 2 = f ( t ) , u ( t ) , (17)

其中 I 1 = μ ( u ( t ) ) Δ u ( t ) , u ( t ) I 2 = σ ( u ( t ) ) u ( t ) , u ( t )

估计 I 1 。易知

I 1 = μ ( u ( t ) ) Δ u ( t ) , u ( t ) = Ω ( 0 u ( t ) μ ( z ) d z ) u ( t ) d x = Ω ( 0 u ( t ) μ ( z ) d z ) u ( t ) d x 0. (18)

估计 I 2 。考虑(14),易得

I 2 = σ ( u ( t ) ) u ( t ) , u ( t ) 0. (19)

回顾(17),从(18)-(19)可得

d d t u ( t ) H α ( Ω ) 2 + u ( t ) H α ( Ω ) 2 C , (20)

其中C依赖f。那么运用~Gronwall~引理,可推断出

u ( t ) H α ( Ω ) 2 C . (21)

引理2:假设 α 1 2 u ˜ 0 ( x ) H α + 1 ( Ω ) ,且

μ C ( R ; R ) , μ ( 0 ) = 0 , μ ( z ) > 0 ( z R , z 0 ) , (22)

σ C 1 ( R ; R ) , σ ( 0 ) = 0 , σ ( z ) > 0 ( z R , z 0 ) , (23)

f L 2 ( 0 , ; H 1 ( Ω ) ) . (24)

那么问题(1)-(3)有一个弱解u满足

u ( t ) H α ( Ω ) 2 C . (25)

证明:将(1)乘以 ( Δ ) u ,关于空间变量x在 Ω 上进行积分,我们有

u ( t ) H α ( Ω ) 2 + μ ( u ( t ) ) Δ u ( t ) , Δ u ( t ) + 1 2 d d t u ( t ) H α ( Ω ) 2 + σ ( u ( t ) ) u ( t ) , ( Δ ) u ( t ) = f ( t ) , ( Δ ) u ( t ) . (26)

类似(18),我们推断

σ ( u ( t ) ) u ( t ) , ( Δ ) u ( t ) 0. (27)

考虑(22),我们可以获得

2 μ ( u ( t ) ) Δ u ( t ) , Δ u ( t ) 0. (28)

另一方面,有

f ( t ) , ( Δ ) u ( t ) f ( t ) u ( t ) 1 2 ε f ( t ) 2 + ε 2 u ( t ) 2 , (29)

其中 ε 是一个正常数并且充分小。根据(26)-(29)有

d d t u ( t ) H α ( Ω ) 2 + u ( t ) H α ( Ω ) 2 C , (30)

其中C依赖f。利用Gronwall引理,我们可以推断出

u ( t ) H α ( Ω ) 2 C . (31)

引理3:令 α 1 2 u ˜ 0 ( x ) H α + 1 ( Ω ) ,u是(1)-(3)有一个弱解,且

μ C ( R ; R ) , μ ( 0 ) = 0 , μ ( z ) > 0 ( z R , z 0 ) , (32)

σ C 1 ( R ; R ) , σ ( 0 ) = 0 , σ ( z ) > 0 ( z R , z 0 ) , (33)

f L 2 ( 0 , ; H 1 ( Ω ) ) . (34)

满足

0 t u t ( s ) H α ( Ω ) 2 d s < C . (35)

证明:将(1)乘以 u t 并关于时间变量t和空间变量x在 ( 0 , t ) × Ω 上进行积分,整理可得

u ( t ) H α ( Ω ) 2 + 2 0 t σ ( u ( s ) ) u ( s ) , u t ( s ) d s + 2 0 t u t ( s ) H α ( Ω ) 2 d s 2 0 t μ ( u ( s ) ) Δ u ( s ) , u t ( s ) d s = u 0 H α ( Ω ) 2 + 2 0 t f ( s ) , u t ( s ) d s . (36)

从(32),我们可以得到

2 0 t μ ( u ( s ) ) Δ u ( s ) , u t ( s ) d s = 2 Ω 0 u ( t ) 0 y μ ( z ) d z d y d x 0. (37)

考虑(25)和引理1,可以获得

2 | 0 t σ ( u ( s ) ) u ( s ) , u t ( s ) d s | 2 σ ( u ( s ) ) L ( Ω ) 0 t Ω | u ( s ) u t ( s ) | d x d s C + δ 0 t u t ( s ) 2 d s , (38)

其中 δ 是一个正常数。

因此,选择 δ 充分小,根据(36)-(38)有

0 t u t ( s ) H α ( Ω ) 2 d s < C , (39)

其中C依赖f。

3. 主要结果及其证明

借助上述引理,我们可以证明下面的定理。

定理1:对于 α 1 2 令初值满足 u ˜ 0 ( x ) H α + 1 ( Ω ) ,且

μ C ( R ; R ) , μ ( 0 ) = 0 , μ ( z ) > 0 ( z R , z 0 ) , (40)

σ C 1 ( R ; R ) , σ ( 0 ) = 0 , σ ( z ) > 0 ( z R , z 0 ) , (41)

f L 2 ( 0 , ; H 1 ( Ω ) ) . (42)

那么问题(1)-(3)有一个弱解u,满足

u L ( 0 , ; H α + 1 ( Ω ) ) , u t L 2 ( 0 , T ; H α ( Ω ) ) . (43)

证明:第一步:Galerkin逼近

H α 上取特殊的正交基 { w j } 由拉普拉斯特征函数形成:

Δ w j = λ j w j , w j C ( Ω ) , 0 < λ 1 λ 2 λ j . (44)

定义

u m ( t ) = j = 1 m c m j ( t ) w j , (45)

其中 c m j ( t ) ( j = 1 , , m ) 满足下面逼近方程

( Δ ) α u m t ( t ) , w j + ( 1 + σ ( u m ( t ) ) ) u m ( t ) , w j + ( Δ ) α u m ( t ) , w j + u m t ( t ) , w j μ ( u m ( t ) ) Δ u m ( t ) , w j = f ( t ) , w j , 1 j m , (46)

初始条件

u m ( 0 ) = u 0 m s p a n { w j , 1 j m } ,

u 0 m u ˜ 0 ( m ) H α + 1 ( Ω ) 中强收敛。(47)

显然,此时的方程(46)-(47)是一组常微分方程。应用标准常微分方程的理论,易知当 0 t t m = T 时,方程(46)-(47)具有唯一解。

第二步:取极限

根据以上估算,存在 { u m ( t ) } 的子序列,仍用 { u m ( t ) } 表示,满足

u m u L ( 0 , ; H α + 1 ( Ω ) ) 中弱 收敛, (48)

u m t u t L 2 ( 0 , T ; H α ( Ω ) ) 中弱 收敛。 (49)

应用Riesz定理,有

u m u L 2 ( 0 , ; H α ( Ω ) ) 中强收敛且几乎处处收敛于 Q T (50)

u m u L 2 ( Q T ) 中强收敛且几乎处处收敛于 Q T (51)

Δ u m Δ u L 2 ( Q T ) 中强收敛且几乎处处收敛于 Q T (52)

考虑引理1、引理2和引理3,我们可以得到

0 u m μ ( z ) d z L ( Ω ) < C , σ ( u m ) u m L ( Ω ) < C . (53)

应用控制收敛定理,我们发现

0 u m μ ( z ) d z 0 u μ ( z ) d z L 2 ( Q T ) 中强收敛, (54)

σ ( u m ) u m σ ( u ) u L 2 ( Q T ) 中强收敛。 (55)

通过(47)-(52)和(54)-(55),在(46)中取极限,我们获得

u t ( t ) , w + ( Δ ) α 2 u ( t ) , ( Δ ) α 2 w + 0 u ( t ) μ ( z ) d z , w + ( Δ ) α u t ( t ) , ( Δ ) α 2 w + ( 1 + σ ( u ( t ) ) ) u ( t ) , w = f ( t ) , w , w H α ( Ω ) , u ( 0 ) = u ˜ 0 . (56)

另外

u L ( 0 , ; H α + 1 ( Ω ) ) , u t L 2 ( 0 , T ; H α ( Ω ) ) . (57)

第三步:解的唯一性

令u和v是(1)-(3)的弱解。定义 w = u v ,我们有

w t + ( Δ ) α w t + ( Δ ) α w ( μ ( u ) Δ u μ ( v ) Δ v ) + w + ( σ ( u ) u σ ( v ) v ) = 0 , w ( 0 ) = 0. (58)

将(58)与w相乘,然后关于时间变量t和空间变量x从 ( 0 , t ) × Ω 进行积分,我们得到

w ( t ) H α ( Ω ) 2 2 0 t μ ( u ( s ) ) Δ u ( s ) μ ( v ( s ) ) Δ v ( s ) , w ( s ) d s + 2 0 t w ( s ) H α ( Ω ) 2 d s + 2 0 t σ ( u ( s ) ) u ( s ) σ ( v ( s ) ) v ( s ) , w ( s ) d s = 0. (59)

利用函数的单调性,我们推断

2 0 t μ ( u ( s ) ) Δ u ( s ) , w ( s ) μ ( v ( s ) ) Δ v ( s ) , w ( s ) d s = 2 0 t 0 u ( s ) μ ( z ) d z 0 v ( s ) μ ( z ) d z , w ( s ) d s 0. (60)

另一方面,我们发现

2 0 t σ ( u ( s ) ) u ( s ) σ ( v ( s ) ) v ( s ) , w ( s ) d s = 2 0 t [ σ ( u ( s ) ) v ( s ) , w ( s ) + σ ( u ( s ) ) w ( s ) , w ( s ) σ ( v ( s ) ) v ( s ) , w ( s ) ] d s . (61)

利用估计,我们有

| σ ( v ( s ) ) σ ( u ( s ) ) | < L 1 | v ( s ) u ( s ) | , (62)

那么,由(62),可知

2 0 t σ ( u ( s ) ) v ( s ) , w ( s ) σ ( v ( s ) ) v ( s ) , w ( s ) d s = 2 0 t ( σ ( v ( s ) ) σ ( u ( s ) ) ) v ( s ) , w ( s ) d s 2 L 1 v ( s ) L ( Ω ) 0 t | | w ( s ) | , w ( s ) | d s L 1 v ( s ) L ( Ω ) 0 t Ω ( w ( s ) 2 + w ( s ) 2 ) d x d s C 0 t w ( s ) H α ( Ω ) 2 d s . (63)

回顾(61),从(63)可推出

2 0 t σ ( u ( s ) ) u ( s ) σ ( v ( s ) ) v ( s ) , w ( s ) d s C 0 t w ( s ) H α ( Ω ) 2 d s . (64)

ρ ( s ) = w ( s ) H α ( Ω ) 2 ,从(58)-(60)和(64),我们有

ρ ( s ) C 0 t ρ ( s ) d s . (65)

此时,应用Gronwall引理,我们有 ρ ( s ) = 0 ,即 w = 0

基金项目

本文受国家自然科学基金项目(No.11271141)和重庆市科学技术委员会(cstc2018jcyjAX0787)资助。

文章引用

程嘉卓,金玲玉,房少梅. 一类非线性伪抛物型方程解的性质
Properties of Solutions for a Class of Nonlinear Pseudoparabolic Equations[J]. 应用数学进展, 2020, 09(04): 551-559. https://doi.org/10.12677/AAM.2020.94066

参考文献

  1. 1. Coleman, B.D., Duffin, R.J. and Mizel, V.J. (1965) Instability, Uniqueness, and Nonexistence Theorems for the Equa-tion on a Strip. Archive for Rational Mechanics and Analysis, 19, 100-116. https://doi.org/10.1007/BF00282277

  2. 2. Barenblatt, G.I., Zheltov, I.P. and Kochina, I.N. (1960) Basic Concepts in the Theory of Seepage of Homogeneous Liquids in Fissured Rocks [Strata]. Journal of Applied Mathematics and Mechanics, 24, 1286-1303. https://doi.org/10.1016/0021-8928(60)90107-6

  3. 3. Benjamin, T.B., Bona, J.L. and Mahony, J.J. (1972) Model Equations for Long Waves in Nonlinear Dispersive Systems. Philosophical Transactions of the Royal Society A. Mathematical, Physical and Engineering Sciences, 272, 47-78. https://doi.org/10.1098/rsta.1972.0032

  4. 4. Ting, T.W. (1963) Certain Non-Steady Flows of Second-Order Fluids. Archive for Rational Mechanics and Analysis, 14, 1-26. https://doi.org/10.1007/BF00250690

  5. 5. Padron, V. (2004) Effect of Aggregation on Population Recovery Modeled by a Forward-Backward Pseudoparabolic Equation. Transactions of the American Mathematical Society, 356, 2739-2756. https://doi.org/10.1090/S0002-9947-03-03340-3

  6. 6. Bona, J.L. and Dougalis, V.A. (1980) An Initial- and Boundary-Value Problem for a Model Equation for Propagation of Long Waves. Journal of Mathematical Analysis and Applications, 75, 503-522. https://doi.org/10.1016/0022-247X(80)90098-0

  7. 7. Amick, C.J., Bona, J.L. and Schonbek, M.E. (1989) Decay of Solutions of Some Nonlinear Wave Equations (ENG). Journal of Differential Equations, 81, 1-49. https://doi.org/10.1016/0022-0396(89)90176-9

  8. 8. Zhang, L. (1995) Decay of Solution of Generalized Benja-min-Bona-Mahony-Burgers Equations in N-Space Dimensions. Nonlinear Analysis Theory Methods & Applications, 25, 1343-1369. https://doi.org/10.1016/0362-546X(94)00252-D

  9. 9. Medeiros, L.A. and Miranda, M.M. (1977) Weak Solutions for a Nonlinear Dispersive Equation. Journal of Mathematical Analysis and Applications, 59, 432-441. https://doi.org/10.1016/0022-247X(77)90071-3

  10. 10. Bona, J.L., Pritchard, W.G. and Scott, L.R. (1980) Soli-tary-Wave Interaction. Physics of Fluids, 23, 438. https://doi.org/10.1063/1.863011

  11. 11. Taha, A. and Mahomed, F.M. (2014) A Note on the Solutions of Some Nonlinear Equations Arising in Third-Grade Fluid Flows: An Exact Approach. The Scientific World Journal, 2014, Ar-ticle ID: 109128. https://doi.org/10.1155/2014/109128

  12. 12. Hayat, T., Shahzad, F. and Ayub, M. (2007) Analytical Solution for the Steady Flow of the Third Grade Fluid in a Porous Half Space. Applied Mathematical Modelling, 31, 2424-2432. https://doi.org/10.1016/j.apm.2006.09.008

  13. 13. Sajid, M. and Hayat, T. (2008) Series Solution for Steady Flow of a Third Grade Fluid through Porous Space. Transport in Porous Media, 71, 173-183. https://doi.org/10.1007/s11242-007-9118-3

  14. 14. Ngoc, L.T.P., Yen, D.T.H. and Long, N.T. (2018) Existence and Asymptotic Behavior of Solutions of the Dirichlet Problem for a Nonlinear Pseudoparabolic Equation. Electronic Journal of Differential Equations, 2018, 1-20.

  15. 15. Kilbas, A.A.A., Srivastava, H.M. and Trujillo, J.J. (2006) Theory and Applications of Fractinal Differential Equations. North-Holland Mathematics Studies, Vol. 204, Elsevier, Amster-dam.

  16. 16. Pu, X., Guo, B. and Zhang, J. (2012) Global Weak Solutions to the 1-D Fractional Landau-Lifshitz Equation. Discrete and Continuous Dynamical Systems Series B (DCDS-B), 14, 199-207. https://doi.org/10.3934/dcdsb.2010.14.199

  17. 17. Guo, B., Han, Y. and Xin, J. (2008) Existence of the Global Smooth Solution to the Period Boundary Value Problem of Fractional Nonlinear Schrodinger Equation. Applied Mathematics & Computation, 204, 468-477. https://doi.org/10.1016/j.amc.2008.07.003

  18. 18. Zhang, S. (2003) Existence of Positive Solution for Some Class of Nonlinear Fractional Differential Equations. Journal of Mathematical Analysis & Applications, 278, 136-148. https://doi.org/10.1016/S0022-247X(02)00583-8

期刊菜单