半线性对流扩散问题的差分流线扩散法.pdfVIP

第一章 引 言 传统的数值方法无法在现有的计算条件下，给出令人满意的数值模拟:或 者大梯度解被过度磨平，或者在急剧变化区域发生数值振荡.如标准Galerkin有 限元便会遇到上述困难.为此，许多非标准有限元方法被相继提出并发展，如 有限体积法([15],[11])、广义差分法([8])、特征有限元法([4],[9J)和流线扩散方 法[[5]等. 第二节 FDSD方法简要回顾 流线扩散法(称为Streamline-diffusionmethod，简称SD方法)是由 T.J.Hughes和A.Brooks在1980年前后提出的一种沿流线方向附加人工薪性的数 值求解对流占优扩散方程的一种新型有限元方法(参见[5]),随后，Johnson和 Ndvert将SD方法推广到发展型对流扩散(包括一阶纯双曲问题)问题 [(习，[6],[10J).它具有良好的数值稳定性及高阶收敛率，己广泛应用于计算流 体等诸多科学工程计算领域. 然而，传统的SD方法使用时空有限元，导致计算工作量较大，对高维或非 线性问题尤其如此，其编程实现较为复杂，对于非线性问题也不便进行线性化 处理.为此，孙澈[13】提出了差分流线扩散法((FiniteDifference-Streanine Diffusionmethod，简称FDSD方法).该方法仅对空间域作SD方法的有限元离散， 而对时间域作差分离散(如Euler离散或C一离散等). 理论分析表明，FDSD方法保持了SD方法的本质特性，简化了SD方法的算 法结构，方便了SD方法的实际应用，对高维和非线性问题的处理具有明显的优 越性.另一方面，由于FDSD方法的Petrov-Galerkin框架和对时空变元的非一致 性处理，使得对该方法的理论分析较SD方法更加困难，孙澈和沈慧[13]就线性 问题给出了FDSD方法的理论分析，证明了FDSD格式仍具有较标准Galerkin方 法更好的稳定性，且在空间方向具有拟最优阶精度. 1997年，张强和孙澈[16]又进一步推广了1【3]的结论，对一类非线性问题给 出了FDSD方法的误差分析.但文1【6]为避免在每个时间步长求解非线性离散化 方程组，未考虑C-N型格式而仅讨论了Euler型格式，使得在时间方向上的精度 仅为O(At)(At为时间步长).于是张强和孙澈 [17]在[16]的基础上，又进一步 提出了FDSD预测校正格式，证明了在l(若(0))模度量下，FDSD格式关于空间 步长h仍具有最优阶收敛性，关于时间步长At的精度则提高到了。((At)z).之 第一章 引 言 后，张强1【9]在文[16]的基础上，针对线性对流扩散问题的FDSD格式，在仅要 求扩散系数非负的条件下，证明了该格式按l(L2(p)模具有最优阶精度. 第三节 本文研究方法简介 ，.3.1本文研究的基础及研究结果 对流占优扩散现象广泛地存在于现实生活的方方面面，由于该现象本身所 具有的特殊性，使用传统的方法模拟其数值解会导致在局部范围内不能得到满 意的结果.如何设计一种更为有效的数值模拟方法，具有重要的现实意义. 本文在工作(([13],[16],[17],[19])的基础上，针对一类半线性对流占优扩 散问题进行了讨论，给出相应的FDSD格式以及稳定性和误差估计.本文给出人 工粘性参数S的具体选取条件，它仅与空间网格参数h有关:理论分析表明该格 式在空间按}(L(}))模具有所谓的拟最优精度，同时在时间方向具有相应的精 度. 1.3.2本文研究中用到的若干辅助知识 1.3.2.1赋范空间和内积空间的若干性质 1赋范空间的性质 设X是域K上的线性空间，!.}是定义在X上的范数，则Vx,yeX,aeK, 有: (1、正定性:I-IIEXllx112蛆llxl卜0。二=。， (2、齐次性:lax卜lall1xl, (3、三角不等式:卜+YI5:llxll+Iyll- 2内积空间的若干性质 设H是域K上的内积空间，(·，·)是定义在H上的内积，1·}是定义在H 上的范数则Vx,yeH，有((x,y)eH,IxlleH且`dx,y,zeH,aeK有: (1)(x,y)20, (2)(ax,y)=a(x,y),