一维非定常对流扩散方程非均匀网格上的高阶紧致-西安理工大学学报.PDFVIP

　西安理工大学学报 ＪｏｕｒｎａｌｏｆＸｉ’ａｎＵｎｉｖｅｒｓｉｔｙｏｆＴｅｃｈｎｏｌｏｇｙ（２０１３）Ｖｏｌ．２９Ｎｏ．４　 ４７５ 　　文章编号：１００６４７１０（２０１３）０４０４７５０６ 一维非定常对流扩散方程非均匀 网格上的高阶紧致差分格式 赵飞，陈建华，葛永斌 （宁夏大学 数学计算机学院，宁夏 银川７５００２１） 摘要：针对一维非定常对流扩散方程，在非均匀网格上提出了一种两层高精度紧致差分格式，其时 间具有２阶精度，空间具有３～４阶精度；然后采用Ｆｏｕｒｉｅｒ分析方法给出格式的稳定性条件。最后 通过数值算例验证了本文格式对于求解一维含边界层非定常对流扩散问题具有明显的优势。 关键词：非定常对流扩散方程；非均匀网格；高精度紧致差分格式；边界层 中图分类号：Ｏ２４１．８２　　　文献标志码：Ａ ＡＨｉｇｈＯｒｄｅｒＣｏｍｐａｃｔＤｉｆｆｅｒｅｎｃｅＳｃｈｅｍｅｏｎＮｏｎＵｎｉｆｏｒｍＧｒｉｄｓｆｏｒｔｈｅ １ＤＵｎｓｔｅａｄｙＣｏｎｖｅｃｔｉｏｎＤｉｆｆｕｓｉｏｎＥｑｕａｔｉｏｎ ＺＨＡＯＦｅｉ，ＣＨＥＮＪｉａｎｈｕａ，ＧＥＹｏｎｇｂｉｎ （ＳｃｈｏｏｌｏｆＭａｔｈｅｍａｔｉｃｓａｎｄＣｏｍｐｕｔｅｒＳｃｉｅｎｃｅ，ＮｉｎｇｘｉａＵｎｉｖｅｒｓｉｔｙ，Ｙｉｎｃｈｕａｎ７５００２１，Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ：Ａｔｗｏｌｅｖｅｌｈｉｇｈｏｒｄｅｒｃｏｍｐａｃｔｄｉｆｆｅｒｅｎｃｅｓｃｈｅｍｅｏｎｎｏｎｕｎｉｆｏｒｍｇｒｉｄｓｉｓｐｒｏｐｏｓｅｄｆｏｒｓｏｌｖｉｎｇ ｔｈｅ１Ｄｕｎｓｔｅａｄｙｃｏｎｖｅｃｔｉｏｎｄｉｆｆｕｓｉｏｎｅｑｕａｔｉｏｎ．Ｔｈｅｓｃｈｅｍｅｉｓｔｈｅｓｅｃｏｎｄｏｒｄｅｒａｃｃｕｒａｃｙｆｏｒｔｉｍｅａｎｄｔｈｅ ｔｈｉｒｄｔｏｔｈｅｆｏｕｒｔｈｏｒｄｅｒａｃｃｕｒａｔｅｆｏｒｓｐａｃｅ．Ｔｈｅｎ，Ｆｏｕｒｉｅｒａｎａｌｙｓｉｓｍｅｔｈｏｄｉｓｅｍｐｌｏｙｅｄｔｏｇｅｔｔｈｅｓｔａｂｉｌｉｔｙ ｃｏｎｄｉｔｉｏｎｏｆｔｈｅｓｃｈｅｍｅ．Ｆｉｎａｌｌｙ，ｎｕｍｅｒｉｃａｌｅｘｐｅｒｉｍｅｎｔｓａｒｅｃｏｎｄｕｃｔｅｄｔｏｄｅｍｏｎｓｔｒａｔｅｔｈａｔｔｈｅｐｒｅｓｅｎｔ ｍｅｔｈｏｄｈａｓｏｂｖｉｏｕｓａｄｖａｎｔａｇｅｓｆｏｒｓｏｌｖｉｎｇｔｈｅ１Ｄｕｎｓｔｅａｄｙｃｏｎｖｅｃｔｉｏｎｄｉｆｆｕｓｉｏｎｐｒｏｂｌｅｍｓｗｉｔｈｂｏｕｎｄａｒｙ ｌａｙｅｒｓ． Ｋｅｙｗｏｒｄｓ：ｕｎｓｔｅａｄｙｃｏｎｖｅｃｔｉｏｎｄｉｆｆｕｓｉｏｎｅｑｕａｔｉｏｎｓ；ｎｏｎｕｎｉｆｏｒｍｇｒｉｄ；ｈｉｇｈｏｒｄｅｒｃｏｍｐａｃｔｄｉｆｆｅｒｅｎｃｅ ｓｃｈｅｍｅ；ｂｏｕｎｄａｒｙｌａｙｅｒ 间均具有四阶精度；文献［ 　　有限差分法是求解非定常对流扩散方程最常用 ４］提出了一种求解一维 的数值计算方法之一。例如古典显格式、古典隐格 非定常变系数对流扩散方程的时间二阶空间四阶的 式、ＣｒａｎｋＮｉｃｏｌｓｏｎ（ＣＮ）格式、Ｓａｕｌ’ｅｖ格式、Ｐｅａｃｅ 差分格式；文献［５］提出了一种求解一维变系数对 ｍａｎＲａｃｈａｆｏｒｄ格式、Ｄｏｕｇｌａｓ格式等。由于这些格 流扩散方程的交替分段 ＣｒａｎｋＮｉｃｏｌｓｏｎ方法。这些 式普遍精度较低，因此又有不少学者提出了高精度 格式大多都是基于均匀网格构造的。文献［６８］针 的差分格式。如：文献［１］构造出了一种求解一维 对定常对流扩散方程，采用坐标变换方法建立了非 非定常对流扩散方程的指数型高阶紧致差分格式， 均匀网格上的高精度紧致差分格式；文献［９１１］不 在时间上与空间上都达到了四阶精度；文献［２］给 采用坐标变换，而是直接在物理区域上建立了非均 出了一种求解一维含源项非定常对流扩散方程的高 匀网格的紧致差分格式。这些差分