原创力文档- 1、本文档共5页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
第13卷第13期2013年5月科学技术与工程Vo1．13No．13May2013
1671—1815(2013)13—3686—05ScienceTechnologyandEngineering⑥2013Sci．Tech．Engrg．
数学
污染物对流扩散方程的预测校正紧差分格式
杨录峰
(北方民族大学信息与计算科学学院，银川750021)
摘要结合预报校正线性多步法与高阶紧致差分格式方法的优点，空间导数采用四阶紧致差分格式进行离散之后，对得到
的空间半离散格式采用改进的预报校正的线性多步法进行时间推进，得到一种时空方向均为四阶精度的求解非线性对流扩
散方程的高精度方法。数值试验表明该格式可以有效求解非线性对流扩散方程，验证了格式的良好性能。
关键词紧致差分格式预报校正线性多步法非线性对流扩散方程
中图法分类号0241．82；文献标志码A
对流扩散方程常常用来描述大气、河流等污染
1数值方法
中污染物质的扩散与分布、工业领域的流动与传热
传质、交通流等众多物理现象。在描述湖泊、河流、对于非线性对流扩散方程
大气及地下水中污染物的浓度分布、流体流动、泥r一O^(u))一t,rtZ~：，t)；0b,0t≤T
沙输移等过程时，经常会遇到求解对流扩散方程。M(0，t)=(t)，(b，t)：(t)；
对这类方程特别是非线性对流扩散方程的数值计0t≤T
算方法的研究具有重要的理论和现实意义。主要(，0)=g()；ab
的研究方法有构造特征差分格式和有限元方法及(1)
其各种变形方法等_1]。谢志华等_4给出了时间推式(1)中：(0，T]是时间域，0，分别表示流体
进采用四阶Runge—Kutta方法，空间导数采用的速度和常黏性系数，，t)(t)(t)，g()均
QUICK格式求解此类问题。本文通过分别求解一为已知的足够光滑函数。
阶导数和二阶导数与函数值线性组合构成的线性为构造微分方程问题(1)的有限差分逼近，首
方程组。针对非线性对流扩散方程本文结合高阶先将求解区域[0，b]×[0，T]进行网格剖分，选取
紧致差分格式与显、隐式Adams方法的优点J，对正整数和Ⅳ，并令空间步长：h=(b—a)／M，时
于空间导数上采用三点四阶紧差分格式进行离散，间步长：丁=Ⅳ，网格点为(，t)，其中：=
得到半离散形式的非线性常微分方程，为避免非线0+ih，i=0，l，…，M，t=，凡=0，1，…，Ⅳ，?
性方程组的计算，时间方向上采用改进的四步预测表示在网格节点(i，t)处的速度值。
校正的线性多步法进行推进，从而使截断误差达到1．1紧致差分格式的构造
O(+h)，最后给出了数值算例验证了该方法具对于一阶导数(对流项)构造三点四阶紧致格
有良好的性能。式，设系数Ⅱ，0，，m，m，m…满足：
文档评论（0）