半线性椭圆问题的Petrov-Galerkin逼近和双调和方程混合元的一种新格式.pdfVIP

己 l 生 ， J . ‘J 有 限元方法是在古典 助 z一Ga lerkin 变分方法 的基础 上 ，以分 片插值 多项式为 工.具 ，结合计算机 的发展与推广而迅速 发展起来 的一种求解微分方程 的数值 方 法 . co urna t在 1943年提 出在三角形 网格上用逐片线性 函数去逼近 dri ichlet 问题 ， 这是有 限元方法最原始 的思想. 到 了二十世纪五十年代 ，有 限元方法为航 空结 构工程师们所发展 ，随后逐步波及到土木结构工程 . 到六 十年代 ，有限元方法在 许 多·领域开始广泛应用 ，包括船舶 ，巨型建筑 的设计 ，流体力学 以及 电磁学等 . 国 内最早研 究有 限元方 法 的是 冯康 院士 (1920一1993)，他独立 与西方 科学家 奠 定 了 有 限元方法 的数学理论基础 (l[]，[2] ，[3]).从而把有 限元方 法从工程局 限 中解放 出 来给‘出了严 密 的数 学描述 . 二十世纪七+ 年代初 ，丑a吞u*初一刀rezz i([4]，[5][6]，[7]，[s])创立 了混合有 限元方法 的一般数学理论 . 二十世 纪八十年代 初 ， Falk和 osbom 提 出了一种弱化 的混合 有 限元分析方法 ([9])，扩大 了混合有 限元方法 的适用 范 围. 混合有 限元方法 可 以 同时逼近 多个 未 知 函数 ，对 处理 高 阶方 程 和含有两个 或 两个 以上 未 知 函数 的方 程更为便 利 . 利用混合 有 限元方 法还 有很 多优 点 ，例如 在计 算 多孔 介 质 流 时 ， 通 常计算速度 ，如用通 常 的有 限元方法 只能先求 出压力 ，然后求 导得到速度 ， 这，样做将 降低精度 ，而利用混合 有 限元方法求解 ，可 同时求 出压力和速 度 ，提 高 了离散解 的精 度 . 另外混合有 限元解 具有很好 的物理 意 义和物理性 能 但 要 求 所构造 的混合有 限元 空 间必须满足 LB B 相容条件 ，给有 限元 空 间的选 取 带来 了 一 定 的 困难 . 有 限元 方 法 在近似 求 解 椭 圆方 程 边值 问题 及 抛 物 方 程 初 边值 问题 方 面 已经 取.得 了显著 成 就 ，数 学理 论 也相 当完 善 . 然而 ，在 对 于 双 曲 问题 的羊邵良元 方 法 上.，无论是 在 其应 用上还 是理论 分析 上都 存在着 较大 的 困难 . 这 是 因为双 曲 问 题 :具有 许 多不 同于椭 圆 ，抛 物 问题 的复 杂性 质 (例 如 ，解 的正 则 性 很 弱 甚 至 可 能 出现 间断 ，存 在扰 动 的传 播 与反射 ，没有 抛 物 问题 那 样 强 的耗 散 性 质 等 等 .).1975 年 ，G .A .B aker([10])提 出 了一类 非标 准 的 o ale改in 有 限元 方 法 ，也 称 Petrov一o ale南n 乙 有 限元法 . 在 1986 年 ，黄 明游在 (1[ l)] 通过修改证 明方法 ，改进 了 ([10D 的误差估 计 ，使 Petrov一G alerkin 有 限元方法具有 比标准 的 G alerkin 有 限元法更 高 阶 的收敛速 度 . Petrov一G alerkin 型有 限元不 同于 标 准 的 G alerkin有 限元 ，它是 同时利 用 两个 不 同的有 限元 函数 空 间，它 们分别作为 ’，试探 函数 ”和 ’，检验 函数 ’，空 间. 这样做 使得方法具有更大 的灵 活性 . 本文就是将 Petrov一Ga lerkin 逼近模式应用于半线性椭 圆问题 ，并且根据插值 算 子 的性 质进 一步证 明了线性有 限元解 的亏量迭代序列收敛 到 Petrov一Ga lerkin 解 ， 用 上述 同样 的单元 ，对 双调和方程 的混合 有 限元 格式 ，采 用对 流 函数用双 二 次 多项式逼近 ，对 涡 函数用双一次多项式逼近 ，在 自然性 光滑条件下 ，这种逼近 可 以获得与 ciarlet一R aviart 格式 中分别用双二次多项式逼近 同样 的收敛 阶. 本文 的写作安排如下 : 第一章 : 介绍预备知