强非线性系统的小波自适应封闭算法及其在力学中的应用

来源 :兰州大学 | 被引量 : 0次 | 上传用户:mini_fc
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
非线性系统尤其是强非线性系统的高精度定量求解技术是非线性科学研究中的共性科学问题,同时也是目前亟待解决的关键难题之一。现存的绝大部分数值方法由于不具备封闭性质导致难以针对强非线性问题建立一套统一的定量求解格式,往往需根据具体问题辅以特别的处理技巧。针对这一缺陷,近年来本研究小组原创性地提出了非线性问题的小波封闭解法,实现了弱强非线性问题的统一求解。然而现有小波封闭解法只能高效适用于非常规则的问题域,即一维域、二维矩形域和三维立方体域,这严重限制了该方法的适用范围。同时,这一小波封闭算法只允许采用全域均匀的节点布置,这在具有局部大梯度问题的求解中往往是难以接受。事实上,复杂区域与局部大梯度问题的高效处理也是目前计算力学领域的热点与难点问题之一。典型地如有限单元法,当问题的研究区域高度复杂或解在局部存在陡峭的梯度时,要生成一组可用的高质量网格将会是一件非常繁琐耗时的工作,而且往往难以自动生成,需专业技术人员全程干预。虽然无网格法可以避免网格生成的难题,但其自身将会带来新的问题,如计算效率低下,使用了对计算精度影响显著但又缺乏定量确定理论的人为参数。针对上述问题,本博士学位论文通过第二代小波建立了可逼近定义在任意二维区域上的连续函数的插值格式,并且该小波插值格式允许进行任意局部节点加密。在此基础上结合本研究小组之前所建立的非线性问题小波封闭算法,提出了一套可适用于任意问题域且允许任意细化局部节点分布的非线性问题的小波多分辨封闭解法,并在此基础上进一步建立了对应的自适应求解格式。同时,通过理论分析与数值实验系统地研究了这一方法的求解精度、计算效率、收敛性与稳定性,并应用这一自适应小波封闭算法定量研究了若干典型的力学和物理问题。首先,本文从第二代插值小波出发,基于多分辨分析创造性的给出了任意二维区域上连续函数的逼近格式,该小波逼近格式不仅仍具有插值性还可以进行局部节点加密。随后从理论上严格推导了利用该逼近格式逼近原函数及其导数的误差估计;并根据尺度函数的基本性质,详细地介绍了形函数的导数、积分和连接系数等精确求解过程,这为逼近原函数和后面微分方程的小波求解提供了数值基础;且通过一些数值算例校验了上述的误差分析且验证了程序的正确性。其次,将第二章中介绍地逼小波近格式与Galerkin法相结合,在本小组已有的工作基础上,提出了一套适用于强非线性椭圆型微分方程的局部增强小波封闭算法。我们利用该统一格式研究了对流扩散方程、Bratu方程和L形区域奇异性问题等。采用γ阶的形函数进行求解,若问题的解是光滑的,则无论是线性还是非线性问题,当前小波封闭解法都将具有γ阶收敛速度;特别当γ=6时,Bratu方程数值结果将明显好于现有的诸多数值算法。而若精确解具有奇异性,近似解的收敛速度将受限于精确解的光滑性,值得注意的是,在相同网格布置下,甚至会出现γ越小,计算精度越好的情形。随后,为了讨论局部节点加密的有效性,对若干具有局部剧变或奇异性问题中的大梯度区域进行局部节点加密。数值结果表明,在保证计算精度的前提下,局部细化可极大地减少布点数目。同时为了估计求解最终代数方程组所需花费的计算量,我们对比了不同γ下的系数矩阵的稀疏性,γ越小,矩阵稀疏性将越好。继而我们将算法推广到线弹性力学问题中,并给出了线弹性力学问题的小波统一求解格式,通过中心圆孔板,L形板和中心裂纹等几个经典力学算例再次验证了局部增强小波算法的高效性,且我们得到的应变能将是实际应变能的上限。进一步,基于小波系数反映的是未知函数在该点处的光滑程度,本文构造了未知函数的小波最佳逼近格式。直接逼近L形区域内的具有局部奇异性函数时,取γ=6,仅需371个不规则节点便得到了与28033个均匀节点相同的数值结果,此时最佳逼近格式效果明显。随后基于最佳逼近格式,提出了一种针对椭圆型问题的小波自适应算法,并运用该小波自适应格式求解了几个具有局部大梯度或奇异性问题。对于局部奇异性问题,取γ=6,我们需要均匀布置39042个点方可达到3.98E-03精度,而在自适应格式下,仅需621个点便能达到如此精度。随即,将自适应格式推广到固体力学之中,也得到了令人满意的结果。通过与规则布置节点的数值结果对比发现:自适应算法可以准确的捕捉到问题的局部大梯度或奇异性所在位置;并可实现在该局部区域进行局部增强;且可通过阈值来控制最终求解精度;在保证最大误差范数的前提下,自适应算法将会显著降低整个域内的节点数目。最后,我们针对非线性初边值问题,运用局部增强小波格式进行空间离散,得到了初边值问题的半离散格式,继而利用隐式-显式混合格式或者Runge-Kutta求解初值问题。分别给出了 Burgers方程,KDV方程,非线性Schr?dinger方程和非线性Klein-Gordon方程的统一求解格式;并证明了本算法具有守恒性质。通过数值算例可以发现本文的小波封闭算法可以简单、高效地处理任意非线性问题,空间收敛速度都能保持在6阶左右,计算精度也高于现有多种数值方法。说明了本文的小波封闭算法对于非线性项并不敏感。且能准确的捕捉到激波、孤立波的运动;并模拟了多种孤立波相互作用的整个过程,甚至是爆破解。
其他文献
阶为素数幂的群称为有限p群.有限p群是有限群的一个重要分支,近年来越来越多的群论学者开始关注有限p群,其中计数问题是有限p群的一个重要研究课题.在研究计数问题方面P.Hall做出了重要贡献,他证明了δk(G)三0(mod pd-k+1),其中G是有限p群,δk(G)表示G的指数为pk且不包含Φ(G)的子群个数,1≤k≤d=d(G),d(G)表示G的极小生成系中生成元的个数.本论文是基于该结果来研究
本文主要利用solid码和句法同余研究语言.首先,在l(S)上定义运算*,其中S为字母表A上一 solid码.证明了(l(S),*)为一幺半群,并对其幂等元和正则元进行了刻画.其次,给出简单的新方法证明:由solid码S确定的同余σs为Λ*上的主同余.同时,利用solid码S定义了Λ*上的另外两个同余λS和ρs,并证明它们也是主同余.再次,讨论了稀疏语言与r-析取各层次语言相乘的情况.证明了:若稀
学位
本文主要研究弱耗散波方程解的整体适定性及其长时间动力学行为.首先,本文在局部一致空间中讨论了R3上自治超三次弱阻尼波方程(1)的初值问题,利用有界域上线性波方程的Strichartz估计证明了方程Shatah-Struwe解的整体存在性和唯一性.由于区域的无界性、非线性项的超临界增长以及方程本身的特性带来的困难,我们发展了[32,48,73]中“收缩函数”的思想方法来证明方程(1)的渐近紧性.进一
学位
学位
理论生态学的基本目标是了解生物个体之间以及与环境之间的相互作用如何决定种群的分布和群落的结构.特别地,非均匀环境在维持生物多样性方面发挥重要作用.反应扩散模型为研究非均匀环境影响种群动力学提供了一个良好的框架.本文旨在研究反应扩散模型在生物学中的应用,具体分析三类空间非均匀环境中的反应扩散模型,重点考虑非均匀环境对种群动力学、物种进化、疾病传播等的影响.第二章,我们研究一类非局部反应扩散对流系统,
Dynkin图是数学领域中一类重要的研究对象.众所周知,Dynkin图给出了复半单李代数、(晶体)根系统、Coxeter群、箭图有限表示型、有限型丛代数、反射幺半群等的分类.用Dynkin图作为桥梁,我们主要研究了丛代数理论和量子仿射代数、反射幺半群表示之间的联系.2010年,两位法国数学家Hernandez和Leclerc定义了丛代数幺半群范畴化的概念,他们首先研究了丛代数和量子仿射代数的有限维
学位
学位