论文部分内容阅读
非线性发展方程在力学领域、生物医学工程领域、控制工程领域等理论研究和工程应用中都有着十分重要的意义。材料力学中的粘性振动和生物医学中的神经传播的模型可归结为一类非线性发展方程一非线性拟双曲方程,对该类方程解的存在性和进行数值求解方法的研究,是材料力学和生物医学研究领域的重要课题。但其高维情形的数值方法和理论分析难度较大。波浪是海洋中最常见的物理现象,其由深海向近岸传播过程中所产生的各种非线性效应对近岸人类生产生活和近岸地貌变化有着较大影响。而近岸浅水区波浪非线性较强且形式复杂,导致数值模拟难度增大。Boussinesq类方程是目前常用的一类浅水非线性发展方程,对该类方程进行理论分析和计算方法研究,为进一步保障近海生产作业、进行海岸工程设计和施工等有重要意义。本文对非线性拟双曲方程,以及在海岸工程中得到广泛应用的Boussinesq类浅水波浪模型这两类非线性发展方程进行了理论分析与计算方法的研究。
针对材料力学中的粘性振动及神经传播过程中涉及的一类非线性拟双曲方程,目前计算方法的研究一般局限于一维问题,其高维非线性情形的理论分析和数值求解难度较大。本文基于有限元配置法,对该类非线性发展方程的二维初边值问题,采用分片双三次Hermite插值多项式空间作为求解的逼近函数空间,建立了半离散和全离散格式,并对两种格式证明了数值解的存在唯一性;应用微分方程先验估计的理论和技巧得到L2最佳阶误差估计。数值实验结果表明,在保证整体误差估计和不增加计算量的前提下,本方法比传统有限元方法有更高的逼近精度,进一步扩展了配置法的应用范围。
基于交替方向有限元法和变网格法,本文提出了针对上述非线性发展方程三维形式的交替方向变网格有限元方法。基于对二维方程的分析,对其三维方程构造了交替方向变网格有限元格式;通过在原方程中加入扰动项,将高维问题化为一系列简单的一维问题逐次求解;应用常规恒等变换的技巧,使得方程能够进行交替方向计算,且可以并行进行。对该格式的误差分析表明,该格式形式简便、隐定,易于并行计算,适用于大规模科学工程计算。
现有Boussinesq类方程,在提高其非线性的同时,通常方程的最高阶导数的阶数也会随之增加,很多情况下都保留了高于三阶的导数项,这给数值计算带来很大困难。通过重新假定波高水深比和水深波长比的关系,本文建立了不增加方程最高阶导数前提下提高精度的新型二维Boussinesq方程。
对近岸浅水区波浪的分析和计算求解中,求解区域范围通常比较大,需要消毛较大的计算量和计算时间。对经典的Boussinesq方程,本文提出了基于区域分解的Boussinesq方程并行计算方法。基于区域分解方法和重叠度概念,提出了用于并行求解的重叠区域分解方法;结合子区域校正算法,采用单位分解函数合理分配重叠区域上的校正量,构造了一类针对Boussinesq方程的新型并行有限差分算法;通过对孤立波进行算例分析,讨论了子区域重叠度和迭代次数对算法的影响。实验结果表明,该并行计算方法可以有效地提高计算效率,对于单位分解函数集取线性函数时,计算中只需要较小的重叠度和较少的迭代次数即可满足实际需要。