复杂区域强非线性力学问题求解的小波方法

来源 :兰州大学 | 被引量 : 2次 | 上传用户:wjsj123
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
伴随着人类知识范围的扩展,非线性科学的地位不断上升。由于非线性模型并不满足叠加原理,不能通过对问题的简单分解来进行量化分析,因此不存在一般的获取精确解的解析方法。为了求解非线性问题,数值计算方法在大多数情况下是惟一可行的选择,并占据着至关重要的位置。另外,实际问题还要求可在复杂区域上执行的算法,而目前的传统方法难以同时处理发生在复杂区域上的非线性问题。注意到本研究组在先前工作中提出的小波方法具有求解非线性方程的强大潜力,本文将其扩展到不规则区域上,提出一种兼顾非线性处理与复杂区域求解的高精度小波数值方法。为了形成一套普适性的求解初边值问题的总体方案,本文还给出了一种计算初值问题的小波多步方法。此外,在前人工作的基础上,本文进一步提高了小波方法对非线性方程的计算精度。Coiflet族小波具有适合数值计算的优良属性。作为基础工作,在滤波器系数组的设计上,本文通过改变消失矩参数的方式,给出了几种属于3N+2族的Coiflet小波。它们比前人工作中构建的3N族小波具备更好的光滑性,可将基函数展开的收敛速度提高一阶。本文在理论上提出了一种适用于Coiflet型小波的改进的计算支撑区间外多重积分值的方法,提供了一种直接的多项式型的解析表达式。该式能够快速计算任意点上的积分值,且不再依赖于滤波系数的介入。这减少了可能的数值舍入误差并提升效率,也为后文的数值积分格式作出了铺垫工作。在小波逼近格式方面,为了减少边界延拓引入的额外误差,本文构造了高阶的Lagrange型延拓格式,克服了原有方案采用的低阶差分格式与小波方法自身的高精度并不匹配的问题。该式在15节点下对tanh(x)的逼近误差可以低至10-8量级,优于其它算法。将该式扩展到二维区域,未发现边界附近的误差有明显增大的现象,证明了其有效性。由于强非线性问题对逼近精度的要求很高,本文构造了一种引入Richardson外推技术的高精度小波配点方法。通过引入半步长的方式并调整系数,能够抵消掉低阶误差项,从而提高了算法的收敛速度。它保留了原算法的全部优点,拥有插值性与高阶光滑性,能够解耦方程中的低阶项与高阶项,使得误差与逼近格式自身无关,且容易施加边界条件,可以无缝替代原算法。最后介绍了一种积分型的小波逼近方法。考虑到未来的工作要在更一般的区间上求解问题,插值点的数量可能会跟随边界形状而不断变化,本文通过使用Newton形式的单向延拓对原有算法进行修改。该格式移除了原算法的一些限制,现可使用任意数量的插值点并在任意长度的区间内施行,这是本文在复杂区域上进行计算的核心之一。对比了4、5与6点方案,我们发现取6节点的延拓已经几乎抑制了边界附近的误差波动。在几何形状复杂的区域上,经典算法往往精度受限,这对非线性问题的计算十分不利。部分精度较好的算法通常难于处理不规则区域,且施加边界条件遭遇到困难。为了兼顾两方面的需求,作为本文的主要工作,提出了一种可在任意形状区域上执行的小波方法。该方法具有良好的泛用性,对边界形式没有特殊要求。它采用了将复杂区域嵌入直角坐标网格的处理方式,无须去拟合复杂的曲线边界,不需要耗费大量时间的网格生成工作,可配合各种简单的网格划分技术以提高效率。小波基的高度光滑性质使此方法具备快速的收敛速度,能够容许在相对粗糙的网格上进行操作,并仍可给出较高精度的结果。小波基的插值性质允许该方法能以简洁的方式操作非线性算子。作为强形式的配点方法,无需将方程修改为弱形式,可以直接求解,对变分原理不存在的某些非线性问题同样有效。其高度的稳定性与合适的边界延拓相结合,避免了其它方法中的系数矩阵病态与边界振荡的弱点。此法还能以精确形式满足不同种类的边界条件,而不是采取某种近似方式来施加。该方法直接生成适合大规模计算的稀疏矩阵,避免了某些经典方法中先离散然后根据边界条件修改总体矩阵带来的低效率。为了分析随时间发展的动态问题,本文提出了一种求解初值问题的小波多步算法。通过调整小波消失矩的参数,可构造出一种强稳定的隐式多步方案。这种方法的导出过程并未借助于传统理论,而是从Coiflet小波近似格式得到。然而其一致性、收敛性与稳定性却能满足经典理论给定的必备条件。绘制出的稳定区域图像与阶星图也能从侧面证明这些属性。利用一种小波逼近给出的预测方法,可以与上述隐式方法合并,从而建立出一套完整的显式的预测-校正方案。若引入Richarson外推技术,这种算法还能进一步加速。我们将会把这种方法与空间上复杂区域的小波算法结合起来,以形成一套总体的初边值问题求解方案。最后,本文通过对一些典型数学方程的计算来展示上述小波方法的优点。由于p-Laplacian方程蕴含了很多数值计算中的难关,其数值解答具备较高的实用价值。在导出新算法的过程中,本文将其作为非线性方程的典型范例进行研究。求解过程中利用了先前建立的小波Galerkin方法与新型小波方法的基本思想。小波方法展现出高精度的特性。其中一例显示小波方法达到了10-7量级的精度,远远好于有限元方法。另一例表明小波方法使用70%左右的节点数便达到了与有限单元法相近水准的精度。与两种有限体积方法的对比,表明小波方法拥有更快的收敛速度。当利用积分型的小波方法求解此问题时,它给出的解与打靶法和有限差分方法几乎完全一致。然而小波方法仅使用1/32的步长,其精度便与差分方法在1/800步长下输出的解大致相当。表明小波型方法具有极高的精度。通过小波Richardson配点方法,计算了数个具有代表性的非线性方程以及一个稳态流动问题。数值结果表明此算法提高了计算精度,其预期行为与理论完全相符,取得了5阶的收敛性能。其中一例显示此法在16节点下的精度已经接近了原方法在32节点下的精度。另一例的结果表明这种新方法计算出的解比原方法更平稳。在不同形状的几何区域上,本文计算了非线性Poisson方程、直杆扭转问题与薄板弯曲问题。小波方法不仅精度优异,对边界的形式也不敏感。相比于有限单元法,小波方法收敛十分快速,在1000个节点以内便能接近有限元方法超过6000节点才能达到的精度,表明其良好的计算效率。其中一例显示出在有限元方法收敛较慢且精度不佳的情形,小波方法仍然有能力计算出高精度的解。多个非线性初值问题的算例展示了小波隐式与显式多步方法的精度与收敛性能。其中一例显示出,对于一些同阶的其它算法不能很好处理的问题,小波多步方法仍可提供较优的计算精度。
其他文献
近几年来,引力波的探测和相关研究成为了宇宙学中最热门的课题,同时也给早期宇宙的探索带来了希望。一些物理学家认为,随机引力波起源于早期宇宙的暴胀场衰变和宇宙弦,因此随机引力波的探测对研究早期宇宙具有重要意义。与随机引力波类似,基本粒子的起源也与宇宙的开端密切相关。而且,研究粒子生成可以弥补用引力波研究早期宇宙的不足,例如,粒子生成可以研究引力波不能提供的宇宙的热力学性质。此外,粒子之间的转换会影响宇
近年来,反应扩散方程及其在种群动力学、化学、物理学、材料学等学科中的应用引起了人们的广泛关注,其中关于传播动力学的研究主要集中在传播速度、行波解和整解等方面.本文主要研究波动(shifting)环境下两种群Lotka-Volterra随机扩散竞争系统和三种群Lotka-Volterra非局部扩散竞争系统的传播动力学.首先研究了波动环境下两种群随机扩散竞争系统的生存和传播理论.假设每个种群的增长率沿
风沙运动通过在海岸与干旱地区形成沙丘、沙波纹影响地球表面形态,沙尘的输运也影响着气候、生态系统、水文循环以及地球系统的其他部分。沙粒与床面的碰撞过程是风沙运动研究的重要部分。本文通过实验及数值模拟对粒-床碰撞过程进行研究,定量表征了反弹、溅射颗粒在不同冲击条件下起跳物理量的分布情况。给出了适用于沙粒的击溅函数,并将坡度的影响加入其中。此外,在数值模拟中,通过起跳角度的分布形态的变化给出了表面蠕移运
输运过程,特别是扩散过程,在物理和生物系统中都是普遍存在的.在相当长的时间里,布朗运动(正常扩散)在扩散领域中占领着独一无二的主导地位.然而,在最近几十年里,非布朗运动的扩散(反常扩散)在实验中被大量地发现,并越来越强烈地撼动了布朗运动的绝对统治地位.布朗运动具有遍历的性质,即时间平均等于系综平均.而人们发现的大量区别于布朗运动的反常扩散过程,往往是非遍历的.进一步地,人们在研究反常非遍历的扩散现
新童年社会学在全球范围内得到不同程度的发展,但依然处于边缘地位,其研究内容包括:"社会中的儿童"研究;"儿童的社会"研究;新童年社会学的本土化研究以及理论化探索。超越长久萦绕新童年社会学的二元论述是其实现理论突破的重点。反思西方新童年社会学理论框架解决中国童年问题的适切性,与社会学等其他学科进行沟通与交流,逐步形成研究网络,与以西方为中心的新童年社会学研究进行理论对话,进而推动新童年社会学研究的理
唯象强子物理的研究对人们理解QCD非微扰问题具有重要的指导意义。随着近年来实验数据的不断积累,目前实验上在寻找包括重味介子的激发态、重味重子以及四夸克态和五夸克态的候选粒子的课题上都取得了非常重要的进展,这些实验进展为唯象强子物理的理论研究提供了丰富的内容。在本论文中,我们先从SU(3)对称性的角度,利用色磁相互作用模型,系统地研究了构型为QQQq四夸克态的质量谱,这其中包括明确的四夸克态ccbq
拓扑磁性结构由于其热稳定性高、尺寸小、状态特征多样化等优点,有望作为高密度非易失性磁存储器、信息处理器、传感器等新型磁电功能器件的基本单元,从而解决当前磁性器件的性能瓶颈。对拓扑磁孤子的磁电特性和电场可控性的探索研究,不仅可以使我们对磁性结构的拓扑相变有更深入的理解,而且还可以为下一代自旋电子器件的发展提供新的开发思路。得益于拓扑磁性结构中磁序的非共线分布及其衍生的磁电相互作用,基于inverse
本文的主要内容是关于相对析取语言(即r-析取语言)与相对正则语言(即r-正则语言)的研究,包括三个部分.第一个部分讨论了涉及i-析取语言的某些分解,并利用这些分解实现了真i-析取语言类的纵向分类.第二个部分给出了析取辖区和f-析取辖区的几个新的刻画,开发了完全稠密语言和solid码的一些新性质,并证明了自由幺半群上不存在含有完全稠密同余类的同余.第三个部分讨论了与内缀码,前缀码和后缀码分别对应的三
自旋电子学的研究旨在将自旋自由度加入传统电子器件,探索利用电子的自旋属性代替电荷功能实现信息记录、操作与读取的可能性,从根本上解决集成器件的高能耗问题。磁性异质结构是由磁性金属(FM)或磁性绝缘体(FI)与非磁性金属(NM)组成的复合结构,能够同时实现自旋流的产生与探测。它不仅是研究自旋输运相关现象的理想体系,而且是各种自旋电子学效应得以功能化的核心部件。在前期关于磁性异质结构自旋输运研究中,人们
自旋波(磁子)在信息处理和传输方面表现出了巨大的潜力,但其低激发效率和高耗散强度限制了它的实际应用。同时,传统的微波器件也面临着芯片微型化的挑战。而在单晶衬底上传播的声波(声子)具有低耗散特性、短波长(GHz频率的典型波长在微米量级)、不能在真空中传播等优势,更适用于通讯和信息处理。如果利用声子作为信息处理和传递的手段与载体,则有望得到新一代低功耗、高效率、微型化的移动通讯芯片。因此,实现声子与磁