波形松弛算法及其在计算流体力学中的应用.pdfVIP

2013年2月 计算数学 第35卷第1期 Feb．，2013 MATHEMATICANUMERICASINICA V01．35，No．1 波形松弛算法及其在计算流体力学中 的应用木1) 杨 熙 (南京航空航天大学数学系，南京210016) 摘 要 本文介绍求解线性常系数微分代数方程组的波形松弛算法，基于Laplace积分变换得到该算 法新的收敛理论．进一步将波形松弛算法应用于求解非定常Stokes方程，介绍并讨论了连续时间 波形松弛算法CABSOR算法和离教时间波形松弛算法DABSOR算法。 关键词：微分代数方程组；鞍点结构；非定常Stokes方程；波形松弛算法 MR(2000)主题分类：65F10，65L80，65N06，65N22，65N40 1．简 介 微分代数方程组的一般形式为 F(t，z，圣)=0， 其中F：Ⅱx耽×现_豫”，这里Ⅱ∈R是一个开区间，眈，D毒∈R”是开集，m，n∈N，N为 自然数集合．方程组(1．1)中圣的含义需要特别注意，它主要包含两种含义：一方面，它表示可 微函数z：ⅡoR”关于其自变量t∈Ⅱ的导数函数；另一方面，它也表示圣是向量值函数F 的独立变量．这样做的意义在于，我们希望能够通过向量值函数F来确定可微函数z使它在 下述意义下成为(1．1)的解函数，即 t∈Ⅱ F(t，z(t)，圣(t))=0，V 对于问题(1．1)，我们需要讨论其解的存在性和唯一性．所谓的唯一性是指：满足如下初值 条件的解函数z是唯一的，即 x(to)=XO， (1．2) 其中to和XO∈R”是给定的初始时间和初始状态；以及满足如下边值条件的解函数z是唯一 的，即 6(z(至)，z(i))=0， (1．3) 其中b：D。×现。Rd，Ⅱ=【±，t-]，对于某些问题其唯一性也会依赖于正整数d．由于微分代数 方程组是由微分方程和代数方程组成的，所以其形态就比单独的微分方程和代数方程都要复 ‘2012年9月4日收到． 1)基金项目：国家自然科学基金资助 万方数据 计算数学 2013正 杂．一方面，它能够反映微分方程的特殊性质；另一方面，它又具有代数方程的某些特点；更重 要的是，还有一些其他现象是由于微分方程和代数方程的组合才出现的．因此，对于微分代数 方程组理论的研究是非常具有挑战性的． 对微分代数方程组的理论分析一直以来有着两种主要思路，一种是几何分析方法，由 簇理论方面的奠基性工作，随后CampbellIs】等建立了关于线性变系数微分代数方程组的理 论…，并将这些理论推广到了非线性微分代数系统[21-24】． 并且取得了很好的数值效果．波形松弛算法现已成为求解大规模微分代数方程组的重要数值 方法，其主要思想是在原始复杂方程组的层面上将问题简单化，即把原来规模较大的问题分解 成许多易于求解的子问题，然后对每个子问题分别单独求解．波形松弛算法是一种迭代算法， 在新的一轮迭代中对各个子问题进行数值积分时，只需要用到其它子问题在前一轮迭代中的 解函数．因此，波形松弛算法非常适合并行求解【15，】，即对每个子问题的求解过程可以由不 同的计算机节点来完成，如果能够适当地分配各个计算机节点的计算任务，则可以获得比较理 想的并行效率．波形松弛算法也是一种连续时间迭代算法，即给定一个解函数的近似函数，在 整个时间区间上计算一个新的近似函数去逼近解函数．在实际求解过程中，需要对连续方程 进行数值离散，因此我们实际