时间周期双稳型反应扩散方程解的长时间行为
Long-Time Behavior of Solutions to the Time-Periodic Bistable Reaction-Diffusion Equation
DOI: 10.12677/pm.2024.145168, PDF, HTML, XML, 下载: 59  浏览: 114 
作者: 陈 卓:上海理工大学理学院,上海
关键词: 时间周期长时间渐近行为比较原理Time-Periodic Long Time Asymptotic Behaviour Comparison Principle
摘要: 本文研究以下一类时间周期反应扩散方程ut=uxxuxf(t,u), x∈R, t>0.的解的长时间渐近行为,其中,f(t,u)是满足双稳型条件且t具有周期性。将通过引入辅助函数,构造适当的上下解,再运用比较原理,可以得到方程解在无穷远处的性质。
Abstract: This paper focuses on a class of time-periodic reaction-diffusion equationsut=uxxuxf(t,u), x∈R, t>0.of solutions with long time asymptotic behaviour, wheref(t,u)satisfies the bistable condition and t is periodic. The properties of the solutions of the equation at infinity will be obtained by introducing auxiliary functions, constructing appropriate upper and lower solutions, and then applying the comparison principle.
文章引用:陈卓. 时间周期双稳型反应扩散方程解的长时间行为[J]. 理论数学, 2024, 14(5): 122-129. https://doi.org/10.12677/pm.2024.145168

1. 引言

反应扩散方程行波解是反应扩散方程中一种特殊的解析解,展示了在空间中持续传播的波动性质。行波解在自然界的各种动态过程中起着关键作用,例如在生态学、化学和物理学领域等。这些解具有一定的空间和时间结构,在空间中以恒定速度传播。行波解的形式通常是 u ( t , x ) = ϕ ( x c t ) ,其中 ϕ ( ξ ) 是关于 ξ = x c t 的特定函数,而c则代表波速。行波解一直是人们研究的重要课题之一,具体可参看文献 [1] [2] [3] 。行波解的特殊性质使其在方程中具有重要的意义,因为它们展示了物质传播的稳定性和动态行为。

就反应扩散方程而言,较为基础性的问题是研究方程行波解的存在性,通常证明解的存在性的方法为上下解结合不动点定理或者相平面分析法,具体可参看文献 [4] [5] [6] [7] [8] 。关于双稳型反应扩散方程解的稳定性,通常运用半群理论、最大值原理与能量方法,相关结果可参看文献 [9] [10] 。

1997年,文献 [11] 研究了一维情形的双稳型反应扩散方程

u t = Δ u + f ( u ) , t R , x R ,

其中, μ ( 0 , 1 ) 且满足 f ( 0 ) = f ( 1 ) = f ( μ ) = 0 ,他们证明了方程存在唯一行波解 ϕ ( x c t ) 满足

{ ϕ + c ϕ + f ( ϕ ) = 0 , ϕ ( ) = 1 , ϕ ( + ) = 0 ,

并且行波解 ϕ ( x c t ) 具有渐近稳定性。

文献 [12] 研究了在具有变漂移系数的半线性反应扩散方程,漂移系数对解的长时间行为的影响。她的

研究结果表明,当反应项f满足 f ( 0 ) = f ( 1 ) = f ( μ ) = 0 lim sup x + u 0 ( x ) < μ 时,则下列方程

u t = u x x + k ( x ) u x + f ( u ) , t > 0 , x R ,

的解 u ( t , x ) 是增长的,那么对一些辅助函数 m ( t ) 和常数 x 0 ,有以下的渐近行为

sup x > 0 | u ( t , x ) φ ( x c t + m ( t ) + x 0 ) | 0 , t + .

文献 [13] 证明了对于不带有漂移系数的时间周期反应扩散方程

u t = Δ u + f ( u ) , t R , x R ,

解的存在性、唯一性与稳定性,证明了其时间周期行波解 ϕ ( t , x c t ) 满足

{ ϕ t c ϕ ξ ϕ ξ ξ f ( t , ϕ ) = 0 , t R , ξ R , ϕ ( t + T , ξ ) = ϕ ( t , ξ ) , ϕ ( t , ) = 1 , ϕ ( t , + ) = 0 , t R .

其中,行波解 ϕ 和速度c仅依赖于反应项f,并且c的符号与 0 T 0 1 f ( t , s ) d s d t 的符号相同。且当 t + 时,若初值 u 0 ( x ) 满足 lim sup x + u 0 ( x ) < μ lim inf x + u 0 ( x ) > ρ ,则对于一些常数 x 0 ,解 u ( t , x ) 满足

sup x R | u ( t , x ) ϕ ( x c t + x 0 ) | 0 , t + .

受到文献 [12] 和 [13] 的启发,本文主要研究如下时间周期反应扩散方程的初值问题

{ u t = u x x + u x + f ( t , u ) , x R , t > 0 , u ( 0 , x ) = u 0 ( x ) , x R , (1)

其中, 0 u 0 ( x ) 1 。假设在 f ( t , u ) 关于t是周期的,也就是说,存在 T R ,使得 f ( t + T , u ) = f ( t , u ) 。接下来,假设 f ( t , u ) [ 0 , 1 ] 上是双稳型的,也就是说存在一个周期函数 T t ,使得

f ( t , 0 ) = f ( t , 1 ) = f ( t , T t ) = 0 ,

并且对于任意的 t R ,函数f满足,在 u ( 0 , T t ) 上, f ( t , u ) < 0 ;在 u ( T t , 1 ) f ( t , u ) > 0 且对于任意 t R ,有

f u ( t , 0 ) < 0 , f u ( t , 1 ) > 0 ,

即,0和1是 f ( t , u ) 的稳定零点。也就是说,存在 ρ ( 0 , 1 2 ) μ > 0 ,使得,在任意 t R ,有

f u ( t , u ) τ , u [ 0 , ρ ] [ 1 ρ , 1 ] , (2)

现在,我们来做一些技术性的假设

0 T 0 1 f ( t , s ) d s d t > 0

也就是说, c > 0 。并且,我们假设方程(1)的解 u ( t , x ) 是增长的,即,当 t + 时,解 u ( t , x ) 局部一致趋近于1。

本文的工作是研究方程(1)的解的长时间行为。

2. 主要结论

定理1.1 如果 lim sup x + u 0 ( x ) < ρ ,并且方程(1)的解 u ( t , x ) 是增长的,那么对于一些常数 x 0 ,则有

sup x > 0 | u ( t , x ) ϕ ( t , x c t + x 0 ) | 0 , t + ,

其中 ϕ 是行波解,c为传播速度。运用于−x轴方向上也有类似的结论。

为了更好的证明定理1.1,本文将会引用文献 [9] 中的一个辅助函数,并且运用文献 [13] 的一个估计。

现在给出

η ( t ) = 2 3 ln 2 3 ( t + 1 ) , t > 0 ,

可知, η ( t ) 是递增函数,并且满足

η ( 0 ) = 0 , sup t 0 η ( t ) 1 ,

0 + e a η ( t ) d t < + , a > 0.

因为 ϕ ( t , ) = 1 ϕ ( t , + ) = 0 ,那么存在 M > 0 ,使得

0 < ϕ ( t , ξ ) < ρ 2 , t R , ξ M (3)

1 ρ 2 ϕ < 1 , t R , ξ M . (4)

因为 ϕ ξ ( t , ξ ) < 0 ,则存在 a > 0 ,使得

ϕ ξ ( t , ξ ) > a , t R , M ξ M . (5)

我们从文献 [13] 中可知,存在常数 C 2 > 0 β > 0 ,使得

| ϕ ξ ( t , ξ ) | C 2 e β | ξ | , t R , ξ R .

q = min { τ 2 , β } .

现在,我们给出解 u ( t , x ) 在初值条件下的上解与下解。

u + ( t , x ) = ϕ ( t , x c t r e ω t α ( t ) ) + r e q η ( t ) , t > T 1 , x 1 ,

u ( t , x ) = ϕ ( t , x c t + r e ω t + α ( t ) ) r e q η ( t ) , t > T 1 , x 1 ,

其中,

α ( t ) = 0 t r C e q η ( s ) d s

并且,对于 T 1 > 0 ω 为常数时,有 ω e ω t 1 , t > T 1

引理2.1 对于任意 r ( 0 , ρ 2 )

(i) 如果有

u 0 ( x ) ϕ ( 0 , x + 1 ) r , x > 0

u ( t , 1 ) 1 r e q η ( t ) , t > 0 ,

则有

u ( t , x ) ϕ ( t , x c t + e ω t + α ( t ) ) r e q η ( t ) , t > T 1 , x 1. (6)

(ii) 如果有

u 0 ( x ) ϕ ( 0 , x 1 ) + r , x > 0

则有

u ( t , x ) ϕ ( t , x c t e ω t α ( t ) ) + r e q η ( t ) , t > T 1 , x 1. (7)

证明:(i) 当 t > 0 x R 时,我们定义

u ( t , x ) = max { ϕ ( t , x c t + e ω t + α ( t ) ) r e q η ( t ) , 0 }

现在,我们首先验证初值条件。将 t = 0 代入 u ( t , x ) ,可以求出

u ( 0 , x ) = ϕ ( 0 , x + 1 ) r u 0 ( x ) , x > 0.

再验证边界条件,当 x = 1 t > 0 的时候,可以得出

u ( t , 1 ) = ϕ ( t , 1 c t + e ω t + α ( t ) ) e q η ( t ) 1 r e q η ( t ) u ( t , 1 )

接下来,仅需要证明

L u = u t u x x u x f ( t , u ) 0.

t > 0 x 1 时候,使得 u ( t , x ) > 0

可以得到,

L u = α ( t ) ϕ ξ ω e ω t ϕ ξ ϕ ξ + r q η ( t ) e q η ( t ) + f ( t , ϕ ) f ( t , u ) = r C e q η ( t ) ϕ ξ ω e ω t φ ξ φ ξ + r q η ( t ) e q η ( t ) + f ( t , ϕ ) f ( t , u )

其中, ϕ ϕ ξ 的取值为 ( t , x c t + e ω t + α ( t ) )

t > T 1 x 1 时,使得 x c t + e ω t + α ( t ) M ,可以得到

1 ρ 2 ϕ ( t , x c t + e ω t + α ( t ) ) < 1 ,

1 ρ u ( t , x ) < 1.

通过 r < ρ 2 ,且根据(1.2)可得,

f ( t , ϕ ) f ( t , u ) τ r e q η ( t ) .

根据 ϕ ξ < 0 , q τ 2 , η ( t ) 1 q e q t 1 ,可以得到

L u τ 2 r e q η ( t ) 0.

对于 t > T 1 x 1 ,使得 M x c t + e q t + α ( t ) M ,那么可以得到 ϕ ξ a

因此,现在取C充分大,使得

L u = r C e q η ( t ) ϕ ξ ω e ω t ϕ ξ ϕ ξ + r q η ( t ) e q η ( t ) + f ( t , ϕ ) f ( t , u ) a r C e q η ( t ) + r q η ( t ) e q η ( t ) ω e ω t ϕ ξ ϕ ξ + f u ( t , u ) L r e q η ( t ) 0.

对于 t > T 1 x 1 ,使得 x c t + e q t + α ( t ) M ,有

0 < ϕ ( t , x c t + e q t + α ( t ) ) < ρ 0 u ( t , x ) ρ ,

现在,可以从(2)得到,

f ( t , ϕ ) f ( t , u ) τ r e q η ( t ) .

根据 ϕ ξ < 0 , q τ 2 , η ( t ) 1 r q e q t 1 ,可以得到

L u r C e q η ( t ) ϕ ξ τ 2 r e q γ ( t ) 0.

因此,综上所述, L u 0

至此,完成了引理2.1 (i)的证明。然后,可以从比较原理得出(6)。

对于(ii)的证明,与证明(i)的方法类似,可以证得 L u + 0

引理2.2 设 u ( t , x ) 是(1)的解,且它是增长的。那么,对于任意 0 < c 1 < c ,存在常数 λ ,使得

sup 0 < x c 1 t ( 1 u ( t , x ) ) = ο ( e λ t ) , t + .

证明:令 0 < c 1 < c ,并且 c = ( c 1 + c ) / 2 。可以根据比较论证和常数 λ > 0 ,证得

1 u ( t , c 1 t ) = ο ( e λ t ) , t + . (8)

现在取

δ : = c c 1 > 0.

然后取 r ( 0 , ρ ) ,使得

r a δ f u ( t , u ) L + τ ,

其中, ρ τ 和a分别由(2)和(5)定义。现在取 γ ( 0 , η ) ,使得

γ 2 γ τ .

t > 0 x R 时,定义

u ( t , x ) : = max { ϕ ( t , x c t ) r e γ x , 0 } .

下面,将要证明 u ( t , x ) 是(1)的一个下解。

根据(6),并且由于 ϕ ( 0 , x ) = Ο ( e μ ξ ) , x + ,所以存在 T 1 > 0 ,使得

u ( T 1 , x ) ϕ ( 0 , x ) r e γ x , x 0

u ( t , 0 ) 1 r , t T 1 .

可知,当 x 0 时,则有 u ( T , x ) u ( 0 , x ) 和当 t T 1 时,有 u ( t , 0 ) u ( t , 0 ) 。那么,只需要证明当 t T 1 x 0 时,有 u > 0

L u : = ( u ) t ( u ) x x ( u ) x f ( t , u ) 0.

通过计算,可以得出

L u = ( c c 1 ) ϕ ξ + ( r β 2 r β ) e β x + f ( t , ϕ ) f ( t , u ) α ϕ ξ + τ r e β x + f ( t , ϕ ) f ( t , u ) ,

其中, ϕ ϕ ξ 取值于 ( t , x c t ) 。接下来,运用(3)中所定义的M,可以得到,当 t T 1 x 0 时,使得 x c t M x c t M ,有 1 ρ u ( t , x ) 1 0 u ( t , x ) ρ

因此,

f ( t , ϕ ) f ( t , u ) τ r e γ x

由于 ϕ ξ < 0 ,则可得 L u 0 。当 r T 1 x 0 时,使得 M x c t M ,则有, ϕ ξ > a 。那么,可有(8)得出

L [ u ] a α + ( τ + f u ( t , u ) L ) r e γ x 0.

综上所述,由比较原理可知,当 t T 1 x 0 时,

u ( t , x ) u ( t , x ) ϕ ( t , x c t ) r e γ x .

由此可得,(8)成立。这就完成了证明。

现在准备证明定理2.1。

证明:设 u ( t , x ) 是(1)的解,固定 r ( 0 , ρ ) ,因为 ϕ ( t , ) = 1 ϕ ( t , + ) = 0 ,使得当 T 2 = k 1 T k 1 Z ,其中T为周期且 x 0 时,

u ( T 2 , x ) ϕ ( 0 , x + 1 ) + r .

通过引理2.1,可知,当 t > T 1 x > 0 时,有

u ( T 1 + t , x ) ϕ ( t , x c t + r e q t α ( t ) ) + r e q η ( t ) .

通过引理2.2,可知,当 T 3 = k 2 T > 0 k 2 Z ,使得

u ( T 3 , x ) ϕ ( 0 , x 1 ) r , x 0

再根据引理2.1,可知,当 t > T 1 x 1 时,

u ( T 3 + t , x ) ϕ ( t , x c t + e q t + α ( t ) ) r e q η ( t ) .

又因为 lim t + α ( t ) < + ,那么存在常数 α 1 α 2 ,使得当 t max { T 2 , T 3 } x 1 时,有

ϕ ( t , x c t + e q t + α 1 ) r e q η ( t T 2 ) u ( t , x ) ϕ ( t , x c t e q t + α 2 ) + r e q η ( t T 1 ) . (9)

取序列 { t n : = n T } n N ,使得 t n + , n + 。再令

u n ( t , x ) = u ( t + t n , x + c t n t n ) .

通过(9),可得,当 t t n + max { T 2 , T 3 } x 1 c t n + e q t n 时,

ϕ ( t , x c t + e q ( t + t n ) + α 1 e q t n ) r e q η ( t + t n T 2 ) u n ( t , x ) ϕ ( t , x c t e q ( t + t n ) + α 2 e q t n ) + r e q η ( t + t n T 1 ) ,

再通过抛物线估计,当 t R x R 时,序列 u n ( t , x ) 局部一致收敛到方程

( u ) t = ( u ) x x + f ( t , u )

的整体解 u ( t , x ) ,可以得到,当 t R x R 时,有

ϕ ( t , x c t + α 1 ) u ( t , x ) ϕ ( t , x c t + α 2 ) .

那么,从文献 [13] 所得的关于 ϕ ( t , x c t ) 的稳定性结果可以知道,存在 x 0 R

使得

u ( t , x ) ϕ ( t , x c t + x 0 ) .

所以,对任意 r > 0 ,存在 N > 0 ,有

| u ( t N , x + c t N e q t N ) ϕ ( 0 , x + x 0 ) | < r ,

通过引理2.1,可以得到对 t > t N x 1 ,有

ϕ ( t , x c t + e q t N + α N ( t ) + x 0 ) r e q γ ( t t N ) u ( t , x ) ϕ ( t , x c t e q t N α N ( t ) + x 0 ) + r e q γ ( t t N ) ,

其中, α N ( t ) = t N t r C e q γ ( s ) d s 。那么,当 t N 充分大时, α N ( t ) 可以任意小,并且由于r可以任意小和 | ϕ ξ |

界,可得,当 t > t N t N 充分大,有

sup x 1 | u ( t , x ) ϕ ( t , x c t + x 0 ) |

足够小。故,可得

sup x 1 | u ( t , x ) ϕ ( t , x c t + x 0 ) | 0 , t + .

这就完成了证明。

参考文献

[1] Britton, N.F. (1986) Reaction-Diffusion Equations and Their Applications to Biology. Academic Press, London.
[2] Fife, P.C. (2013) Mathematical Aspects of Reacting and Diffusing Systems. Vol. 28. Springer Science & Business Media, Berlin, Heidelberg.
[3] Ma, S. and Zou, X. (2005) Existence, Uniqueness and Stability of Travelling Waves in a Discrete Reaction-Diffusion Monostable Equation with Delay. Journal of Differential Equations, 217, 54-87.
https://doi.org/10.1016/j.jde.2005.05.004
[4] Berestycki, H. and Nirenberg, L. (1992) Travelling Fronts in Cylinders. Annales de lInstitut Henri Poincaré C, Analyse non linéaire, 9, 497-572.
https://doi.org/10.1016/s0294-1449(16)30229-3
[5] 叶其孝, 李正元. 反应扩散方程引论[M]. 北京: 科学出版社, 1990.
[6] 王明新. 非线性抛物型方程[M]. 北京: 科学出版社, 1997.
[7] Guo, H. and Monobe, H. (2021) V-Shaped Fronts around an Obstacle. Mathematische Annalen, 379, 661-689.
https://doi.org/10.1007/s00208-019-01944-y
[8] Hamel, F., Monneau, R. and Roquejoffre, J.M. (2005) Existence and Qualitative Properties of Multidimensional Conical Bistable Fronts. Discrete and Continuous Dynamical Systems-Series A, 13, 1069-1096.
https://doi.org/10.3934/dcds.2005.13.1069
[9] Xin, J.X. (1992) Multidimensional Stability of Traveling Waves in a Bistable Reaction-Diffusion Equation, I. Communications in Partial Differential Equations, 17, 1889-1899.
https://doi.org/10.1080/03605309208820907
[10] Levermore, C.D. and Xin, J.X. (1992) Multidimensional Stability of Traveling Waves in a Bistable Reaction-Diffusion Equation, II. Communications in Partial Differential Equations, 17, 1901-1924.
https://doi.org/10.1080/03605309208820908
[11] Fife, P.C. and McLeod, J.B. (1977) The Approach of Solutions of Nonlinear Diffusion Equations to Travelling Front Solutions. Archive for Rational Mechanics and Analysis, 65, 335-361.
https://doi.org/10.1007/BF00250432
[12] Uchiyama, K. (1985) Asymptotic Behavior of Solutions of Reaction-Diffusion Equations with Varying Drift Coefficients. Archive for Rational Mechanics and Analysis, 90, 291-311.
https://doi.org/10.1007/BF00276293
[13] Alikakos, N., Bates, P. and Chen, X. (1999) Periodic Traveling Waves and Locating Oscillating Patterns in Multidimensional Domains. Transactions of the American Mathematical Society, 351, 2777-2805.
https://doi.org/10.1090/S0002-9947-99-02134-0