论文部分内容阅读
在当代,微分方程无处不在,各个科学领域的研究都围绕着微分方程模型.为了与实际相符,模型形式日趋复杂,比如地震波波动模型.只有经典的原始的微分方程才可以求得解析解,对于大部分地震波波动模型目前只能简化以后进行数值模拟.随着研究的深入,对于更复杂的地震波传播模型,在数值模拟不易进行时,考虑研究解的定性理论,也就是不求解直接研究解的分布和性态,从而探讨地震波的传播特征.方程解的振动性是微分方程定性理论的重要分支.本文的研究内容分为两个部分,第一部分是在常微分方程解的振动性的基础上讨论了中立型时滞脉冲偏微分方程和方程组、分数阶脉冲偏微分方程和分数阶脉冲时滞偏微分方程组解的振动性.在振动性的讨论中,利用平均值方法将偏微分方程转化为常微分方程或者不等式,从而得到偏微分方程解的振动性,并尝试将振动性的研究运用于各向同性声波方程.在分数阶偏微分方程振动性的讨论中,分别利用变量代换以及分数阶导数定义与Γ函数的关系两种不同的方法将分数阶转化为整数阶.第二部分,将Entropy-TVD格式推广至一维浅水波方程,并利用三个数值实例验证了Entropy-TVD格式的有效性,并将这个格式与标准的Godunov格式在分辨率、数值精度阶数和计算成本等方面进行了比较.论文取得的主要成果和结论如下:(1)本文研究了两类时滞脉冲偏微分方程及方程组的振动性.利用平均值法、格林公式和边界条件将所要研究的非线性脉冲时滞双曲方程边值问题解的振动性转化成二阶脉冲微分不等式解的振动问题,接着利用Riccati变换将这个二阶脉冲微分不等式降为一阶,利用辅助函数得到所求边值问题解振动的充分条件.在研究一类中立型脉冲时滞抛物系统在两类边界条件下解的振动性时,首先利用平均值法、格林公式、边界条件以及垂直相加法将脉冲时滞偏微分方程组转化为脉冲时滞常微分不等式组.接着利用变量代换来处理脉冲项,将复杂的分段连续情形转化为连续的状态来考虑,将所研究的问题转化为普通一阶常微分不等式解的振动问题.这样的处理可以极大限度地让已有的大量的一阶常微分方程或者不等式解的振动理论得到推广应用,使得研究空间更为广泛.尝试将微分方程振动理论运用于各向同性声波方程中,并得到结论.(2)基于分数阶微分方程在反常扩散、多孔介质力学、非牛顿流体力学等学科中的广泛应用,本文讨论了一类分数阶脉冲偏微分方程和一类分数阶脉冲时滞偏微分系统解的振动性.基于分数阶导数给研究带来的困难,本文采用两种不同的方法将分数阶偏导数转化为整数阶导数,这样就可以利用已有的整数阶微分方程解的振动理论处理分数阶微分方程解的振动性.本文采用的第一种方法是直接利用Γ函数进行变量代换,第二种方法是利用Modified Riemann-Liouville分数阶导数与Γ函数之间的关系.对于转化之后的微分方程,综合应用Riccati变换和微分不等式,得到了这两类分数阶脉冲偏微分方程在不同边界条件下解的振动准则.(3)本文将Entropy-TVD格式推广至一维浅水波方程.首先详细描述了Entropy-TVD格式,介绍了这个格式的一些性质然后运用于一维浅水波方程.给出了三个数值实例,表明了Entropy-TVD格式的有效性,并研究了Entropy-TVD格式的数值精度阶数和计算成本.Entropy-TVD格式比标准Godunov格式减少了数值耗散,具有更好的分辨率.为了提高线性特征场和非线性特征场的精确度,本文还建立了两个HS重构并将深度和速度作为两片常函数.Entropy-TVD格式包含四个数值实体,即数值熵、数值速度、深度和流量.文中验证了这个格式保留了深度和流量守恒,而且满足熵条件.本论文的创新之处主要表现在以下三个方面:(1)在对偏微分方程解的振动性的讨论中,利用Green公式的推导更好地处理了非线性项,有助于处理非线性地震波波动方程.利用Riccati变换对所研究的二阶常微分方程组进行降阶,使研究更为简便.利用变量代换将分段连续函数转化为连续函数,更有效地处理了脉冲项.这样可以处理更多的存在多种突发扰动的系统.将振动理论运用于声波方程,为研究复杂介质中或者更复杂的比如带有脉冲和时滞的波动模型提供理论基础.(2)在对分数阶微分方程的讨论中,其中分数阶导数的定义采用Modified Riemann-Liouville分数阶导数,修正了原先推导中的漏洞.目前,在对分数阶微分方程解的振动性的讨论中分数阶偏微分方程并不多见,带脉冲时滞的方程少之又少,基本上没有对偏微分方程组进行讨论.本文利用整数阶变量代换的方法处理了所讨论方程中的脉冲项,并利用垂直相加法得到了分数阶脉冲时滞偏微分方程组解的振动性.(3)本文将一阶精确Entropy-TVD格式推广到了一维浅水波方程,为了提高线性特征场和非线性特征场的精确度,建立了两个HS重构并把深度和速度作为两片常函数.这个格式包含四个数值实体,即数值熵、数值速度、深度和流量.Entropy-TVD格式比标准Godunov格式减少了数值耗散,具有更好的分辨率.(4)本文将熵格式推广到地下水溶质运移方程,首先采用分裂方法将地下水溶质运移方程分成对流方程和弥散方程,对流方程是一个双曲型方程,利用熵格式求解,弥散方程的空间离散用二阶中心格式离散时间离散用简单的向前差分.通过数值试验,对不同对流强度的地下水溶质运移方程进行了数值计算,计算结果表明熵格式没有出现过量问题,没有出现非物理振荡,数值弥散小,特别适合强对流问题的数值计算.