求解对流扩散方的四种差分格式的比较.PDF

第 53 卷 第 10 期 2004 年 10 月 物 　理 　学 　报 Vol . 53 ,No . 10 ,October ,2004 10003290200453 ( 10) 325807 ACTA PHYSICA SINICA 2004 Chin . Phys. Soc . 求解对流扩散方程的四种差分格式的比较 ) ) ) ) ) 何文平1 　封国林1 2 　董文杰3 　李建平2 1) (扬州大学物理科学与技术学院 ,扬州 　225009) 2) ( 中国科学院大气物理研究所 ,北京 　100029) 3) ( 国家气候中心 ,北京 　100081) (2003 年 12 月 19 日收到 ;2004 年 2 月 19 日收到修改稿) 　　利用对流扩散方程 ,在边界和参数存在随机扰动的情况下 ,考察四种差分格式的优劣 ,为求解对流扩散方程提 供一种可靠的差分格式 ,并得到通过空间加密网格的方法可以控制边界 、参数随机影响的结论. 关键词 : 对流扩散方程 , 差分格式 , 随机扰动 PACC : 0340 G , 9260X 参数自忆回溯时间积分格式 ,并应用差分理论讨论 1 引 言 了其稳定性. 数值方法是一种离散化方法 ,常用于解决初值 对流扩散方程是一类基本的运动方程 ,它可描 问题的微分方程(组) ,利用这一方法可以在一系列 述质量 、热量的输运以及反应扩散等众多物理现象 , 离散化点 x 1 , x2 , …, xN 上求出未知函数 u ( x ) 在各 许多物理问题的数学模型都可以归结为对流扩散方 点的近似值 u1 , u2 , …, uN , 自变量的离散值是事先 程(组) . 对流扩散方程被广泛应用于流体力学、化 确定的 , xj 称为节点 ,通常取为等间距 , 即 xj = j h ,j 学 、生物 、气象等领域 ,如化学反应过程 、液体在生物 = 1 ,2 , …,N ,其中 N 为空间网格点总数 , h 称为空 膜两边的渗透等问题都涉及到扩散过程 ,因此 ,为对 间步长 ,必要时可以改变它的大小 , u1 , u2 , …, uN 通 流扩散方程寻找一种快速 、稳定 、实用的数值方法有 常称为初值问题的一个数值解 ,它们还与积分时间 着重要的理论和实际意义[ 1 ,2 ] . 这方面的研究已取得 长度 t 有关 ,因此可用 un 来表示未知函数 u 在离 n j [3 ] 了不少进展 ,刘星等 利用数值方法对采用均分网 ( ) 散点 xj , tn 的近似值. 差分法的核心是将微分方程 格的一维线性无源对流扩散方程在各种边界条件下 化成差分方程