非定常边界层方程的微求积解法分.pdfVIP

第lo卷第2期 计 算 物 理 v。I．10，N。．2 l993年 6 月 CHINESE JOURNAL OF COM PUTATIONAL PHYSICS Jun．，1993 非定常边界层方程的微分求积解法 昌 L／／7 耳 f南京航空学院，2100I6) 摘 要 对微分求积法的关键问题一加权系数的确定进行 了研究。在 Lagrange插值 逼近理论的基础上发展 了计算备加权系数的方法，克服了原方法Ill对 同格节点选取有严格 限制的局限性，使之更适用于流体力学问题的研究。该方法用于求解非定常边界层方程的 结果表明，仅用较少的网格点就能得到足够精度的数值解。 ^ ． 关键词 挑 常 边 鲁 ， ， 劳奄霹 1 引 言 微分求积法 (Di腩rentialQuadratureMethod)I。。是一种快速求解非线性偏微分方程 的新方法 为保证足够的计算精度，在数值计算中得到广泛应用的有限差分和有限元法需将 求解区域刊分成的网格或节点数 目是相当大的 (即使仅需求解域 中少数点的函数值 )，从而 导致大的计算量 。而微分求积法仅需少量网格点就能得到足够精度的计算结果，这是该方法 的显著特点。本文发展了该方法，推导出从一阶到高阶，从一维 到多维 的加权系数表达式， 并用于非定常边界层方程的计算 2 微分求积法 1)概述 微分求积法是将函数在网格节点处对空间某方 向的偏导数表示成该坐标方 向所有网格点 上的函数值在加权意义下的和，使得偏微分方程离散为常微分方程组 (非定常问题 )或代数 方程组 (定常 问题 )而便于求解 。 以一维非定常问题为例 假设待求函数 “( ，)充分光滑，XEa【，6】，则 “在节点 处 对 的一阶导数可表示为 ： “(X，)=∑4 “(，) (1) J I 式中．日为加权系数， ∈I“，61，i=1，2，．．．，Ⅳ 将 (I)式代人给定的非定常偏微分方程后 ，将得到一组关于时间f的常微分方程，从而 可用经典方法 ，如 Runge—Kutta或 Adams—Moulton法 ，积分求解 。 Bellman等人采用了两种试验函数族 【·I以得到 (1)式中的加权系数。需要指出的是， 由 于高阶 Vandermonde型矩阵求逆的困难以及节点须选 ShiftedLegendre多项式之根的限制 而影响其应用。 2)各阶偏导数的加权系数 991年 12月 18日收到原稿，I992年 9月30 日收到修改稿 计 算 物 理 第 10卷 为克服上述问题，本文在函数的Lagrange插值逼近的基础上，发展了计算加权系数的 方法．仍以一维非定常问题为例。设函数 “( ，)， ∈【口，6】用 (N一1)次 多项式逼近 ， 选 用 Lagrange插值基函数，略去截断误差后得 ： “(，)=∑rs()“(，f) (2) J；1 式中，Lagrange插值基函数 ‘()为