1. 引言
奇异摄动边值问题因来源于物理、力学、生物等领域而受到广泛的关注和研究 [1] [2] [3] [4] ,奇异摄动理论和方法已经成为研究非线性问题的有力工具之一。近几年,在流固耦合 [5] 、对流扩散 [6] 、界面交互 [7] 等问题中出现了大量具有不连续系数和源项的奇摄动边值问题。例如,2012年,丁海云等在文 [8] 中考虑了一类带有不连续项的二阶半线性奇摄动边值问题,运用边界函数法和缝接法获得了一致有效的渐近解。2017年,包立平在文 [9] 中研究了具有不连续源的半线性微分方程组边值问题,用 Aumann 介值定理对解在间断点处进行了光滑地连接,并证明了解的一致有效性。
最近,Rao和Chawla讨论了带有边界层和内部层的奇摄动半线性反应-扩散方程 [10] :
其中,
,且
,
,
,
他们建立了一种基于Shishkin网格的数值计算方法。
受上述工作的启发,本文在更一般的条件下讨论如下具有不连续源项的二阶半线性系统的奇摄动边值问题:
(1.1)
其中
为
上的充分光滑的函数,
为
的参数,
为给定常数,函数
在
处间断。问题(1.1)可以看作纯量情形 [8] [11] 的自然推广。
令
,则(1.1)式等价于如下的四维吉洪诺夫系统:
(1.2)
作如下假设。
(H1)函数
定义如下:
其中,函数
在
,
上充分光滑,且
。
(H2)方程组
,
,分别存在唯一解
。
(H3)存在
和
使得对任意
成立
由假设(H3)可知,4阶矩阵
的特征根有2个正实部和2个负实部,式中
为2 × 2单位矩阵。因此,辅助系统
(1.3)
和
(1.4)
的平衡点(0,0)皆具有两维的稳定流形和两维的不稳定流形。记系统(1.3)的稳定流形为
;记系统(1.4)的不稳定流形为
。
存在可逆矩阵
使得
其中,
均为2阶矩阵,且
的特征值的实部为负值,
的特征值的实部为正值。
(H4)
且
。
注:若(1.1)为线性系统,即若
是
和
的线性函数,则条件(H3)自然蕴含(H4)成立。
2. 构造形式渐近解
考虑到源项
在
处间断,把问题(1.1)看成如下左、右问题的耦合。
左问题(L)
(2.1)
右问题(R)
(2.2)
其中
,
为依赖于
的待定参数。
按照缝接法思想,我们首先构造出左问题(L)和右问题(R)的形式渐近解,然后在间断点
处对形式解进行光滑地连接,由此确定待定参数
和
先考虑左问题(L)。由边界函数法 [4] 可设其形式解为:
(2.3)
式中
(2.3a)
(2.3b)
(2.3c)
其中
。
将式(2.3)代入式(2.1)中,再根据不同的尺度分离可得:
(2.4)
(2.5)
(2.6)
对
进行
幂级数展开:
(2.7)
其中,
为
的函数。
将式(2.3a)和(2.7)代入式(2.4)中,得
(2.8)
(2.9)
其中
为
的函数,由假设(H1)和(H2)可知(2.8)式有唯一解,又由假设(H3)知,从式(2.9)可依次求出
和
。
将(2.3b)代入(2.5)中,得
(2.10)
(2.11)
其中,
为
的函数。
由边界层函数的性质可知:
。
利用文 [4] 中引理3.5.6,由假设(H3)和(H4)可知,存在
使得
同理,由式(2.11)知:存在
,且
,使得
将(2.3c)代入(2.6)中,可得
(2.12)
其中,
为
的函数。
由边界层的性质可知:
。
利用文 [4] 中引理3.5.6,由假设(H3)和(H4)可知:
其中
。
类似地,可求出右问题的形式渐近解
(2.13)
其中
,
分别是左、右边界层校正项,且满足
其中,
。
为左、右问题的解在
处光滑连接,须有
(2.14)
将式(2.3),(2.13)代入(2.14),考虑到
在
处的指数衰减性,从而有
(2.15)
令
,对(2.12)变形为:
两边对
在
上积分得:
即
类似地,有
由连接条件(2.15)知,
即
考虑第二个分量,类似可得
其中,以上诸式中
,
(H5)假设方程关于
的方程组
有唯一解
。
至此,我们构造出原问题(1.1)的形式渐近解
(2.16)
3. 主要结论及证明
定理1 假设(H1)-(H5)成立,对于充分小的
,问题(1.1)存在解
,满足
(3.1)
其中
为式(2.16)所定义。
证明:构造辅助函数
其中,
显然,有
。由边界条件易得,
对
,由条件(H3)及中值定理得
由形式解的构造过程知
因此,存在
,当
时,有
那么,取
时,有
对
,存在
,
,当
时,
同理可证,在
上有
类似可证,在
上有
因此,
和
分别是问题(1.1)的一个下解和上解。由微分不等式定理知 [2] ,问题(1.1)存在解
且满足(3.1)。证毕。
4. 例子
考虑如下的边值问题:
(4.1)
其中,
显然,由
可知,
,且由
知
。
可将问题(4.1)看成如下左、右问题的耦合。
左问题(L):
(4.2)
右问题(R):
(4.3)
其中
,
为
有关的待定常数。
由渐近解的构造方法可求出左、右问题的零阶近似:
左问题(4.2)的边界层函数
:
右问题(4.3)的边界层函数
:
其中,
。
由式子(2.24),可得
,因此边值问题(4.1)的零阶近似解为: