构造对流扩散方程高精度非均匀网格格式的指数变换方法  

AN EXPONENTIAL TRANSFORMATION METHOD FOR CONSTRUCTING HIGH-ORDER SCHEMES OF CONVECTION-DIFFUSION EQUATIONS ON NON-UNIFORM DIFFERENCE GRIDS

作　　者：王烜[1] 杨志峰[1] 华国春[2] 孟占利[1] 

机构地区：[1]北京师范大学环境学院水环境模拟国家重点实验室,北京100875 [2]四川大学水利水电学院,成都610065

出　　处：《北京师范大学学报（自然科学版）》2005年第6期636-640,共5页Journal of Beijing Normal University(Natural Science)

基　　金：教育部博士点基金资助项目(20040027007);国家自然科学基金委员会;水利部黄河水利委员会黄河联合研究基金资助项目(50239010);武汉大学水沙科学教育部重点实验室开放基金资助项目(2003A002)

摘　　要：提出了一种基于非均匀差分网格,构造求解对流扩散方程的高精度格式的指数变换方法.引入指数函数,将对流扩散方程变换为扩散反应方程,消除了数值求解中较难处理的对流项.采用优化差分方法推导出扩散反应方程基于非均匀网格的高精度差分格式,进而通过逆变换得到对流扩散方程的高精度格式.理论分析表明,该方法具有3至4阶精度,当计算区域为均匀网格时取得4阶精度.数值实例表明,在相同的非均匀网格系统中,此方法的计算精度明显优于传统的隐式差分方法.在水环境的实际模拟计算中,根据物理量的变化规律灵活地调整非均匀网格的间距,不仅能增强高精度差分方法的实用性,而且可以取得比均匀网格方法更为精确的计算结果.An exponential transformation method for constructing high-order schemes to solve unsteady convectio-diffusion equations （CDE） on non-uniform difference grids is developed. To eliminate the difficulty in dealing with convection terms, an exponential transformation function is introduced to turn the CDE into their equivalent diffusio-reaction equations （DRE）. Through an extension of optimal difference method to nonuniform grids, the high-order compact scheme for the DRE is obtained, and then the corresponding scheme for the CDE is obtained through converse transformation. Theoretically, the proposed method has the thirdto fourth-order accuracy; its fourth order accuracy is achieved under uniform-grid settings. Two examples are provided to examine the performance of the proposed scheme. Compared with the traditional implicit difference scheme, the proposed scheme can give better results in the same nonuniform grid system. In the actual computation of water environment simulation, high accurate methods using non-uniform grids with proper spacial steps would not only be more practical, but also produce more accurate numerical results than those using uniform grids, with the same number of grid points.

关 键 词：对流扩散方程 差分方法 非均匀网格 指数变换 高精度 

分 类 号：O241.3[理学—计算数学]