论文部分内容阅读
伴随着人类知识范围的扩展,非线性科学的地位不断上升。由于非线性模型并不满足叠加原理,不能通过对问题的简单分解来进行量化分析,因此不存在一般的获取精确解的解析方法。为了求解非线性问题,数值计算方法在大多数情况下是惟一可行的选择,并占据着至关重要的位置。另外,实际问题还要求可在复杂区域上执行的算法,而目前的传统方法难以同时处理发生在复杂区域上的非线性问题。注意到本研究组在先前工作中提出的小波方法具有求解非线性方程的强大潜力,本文将其扩展到不规则区域上,提出一种兼顾非线性处理与复杂区域求解的高精度小波数值方法。为了形成一套普适性的求解初边值问题的总体方案,本文还给出了一种计算初值问题的小波多步方法。此外,在前人工作的基础上,本文进一步提高了小波方法对非线性方程的计算精度。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节点才能达到的精度,表明其良好的计算效率。其中一例显示出在有限元方法收敛较慢且精度不佳的情形,小波方法仍然有能力计算出高精度的解。多个非线性初值问题的算例展示了小波隐式与显式多步方法的精度与收敛性能。其中一例显示出,对于一些同阶的其它算法不能很好处理的问题,小波多步方法仍可提供较优的计算精度。