对流占优Burgers方程普遍性随流格式计算稳定性分析.pdfVIP

对流占优Burgers 方程普遍性随流格式的计算稳定性 分析1 杨晓忠，包树蕊，侯晓军 华北电力大学数理学院，北京（102206 ） E-mail：xzy1018@ ，baoshurui_love@163.com 摘 要：本文以最简单的非线性对流扩散方程—Burgers 方程为模型，着重讨论了一类较普 遍的差分格式的计算稳定性问题，用稳定性的启发性分析方法给出了这类格式的稳定性判 据，数值试验证明了得到的稳定性判据的确是保证差分格式计算稳定性的必要条件。 关键词：计算稳定性，随流格式，启发性分析法，数值试验 中图分类号：O241 1. 引言 从物理上分析，耗散特征无论是粘性（动量输送的体现）、热传导（热量输送的体现）或 扩散（质量输送的体现），它都使得原有的场均匀化，从波的角度分析也是类似的。在非线 2 2 ∂ u ∂u ∂u ∂ u 性平流方程中加进形如γ 2 耗散项，就得到著名的 Burgers 方程 +u γ 2 ，其特 ∂x ∂t ∂x ∂x 点是有光滑解，在物理或者力学中，表现为具粘性的耗散过程或紊流现象，它和激波的结构 与湍流的研究紧密相关[1] 。众所周知，对非线性偏微分方程，由于非线性的原因，处理线性 问题常用的叠加原理、积分变换等方法就不灵了，至今还没有处理一般非线性方程的统一方 法，所以对特殊类型的非线性方程，只能用相应的特殊方法去处理。 在数值天气预报的实际计算中，人们常常会遇到这样或那样的计算不稳定现象。为了克服 这些不稳定，一般常用的办法是加光滑或引入耗散项。实践证明，这种加耗散措施确实常常 能克服计算不稳定，但有时仍会出现计算不稳定，甚至还会出现耗散系数γ 加大反而是不稳 定更甚的现象，因此我们对 Burgers 方程数值格式的计算稳定性认真地做些分析是很有意义 的。 文献[2]中蒋锦良给出了对流占优 Burgers 方程的一类新方法— 随流格式，本文用稳定性 的启发性分析方法[1,3]，直接分析迎风随流格式，着重讨论了耗散因子、计算初值和时间步 长对计算稳定性的影响，我们发现 Burgers 方程的计算不稳定主要与原 Burgers 方程解的性 质有密切关系，即主要与耗散系数γ 值的大小和初值密切有关，从而进一步加深了我们对非 线性计算不稳定产生机理的认识，推广了我们以前的工作。 2. 方程和随流普遍性差分格式 考虑一维非线性对流占优 Burgers 方程的初值（Cauchy ）问题： 2 ∂u ∂u ∂ u +u γ 1 2 () ∂t ∂x ∂x u x,0 φ x 2 ( ) ( ) ( ) 其中γ 为大于零的常数。 在研究对流占优型一维 Burgers 方程时,文献[2]提出了随流格式，即 1本课题得到河北省自然科学基金（A2007001027 ）的资助。 - 1 - n+1 n n n 2 n+1 u −u u ⎛ ⎞