第五章对流扩散问题(假扩散高阶格式差分方程的求解).pdf

第五章 对流扩散问题———假扩散 A5 高阶格式差分方程的求解 高阶格式的格式图案如下图所示： WW W P E EE 1D：5点格式 NN NN BB N N B WW W P E EE WW W P E EE S S T SS TT SS 2D ：9点格式 3D ：13点格式 第五章 对流扩散问题———假扩散 当对流扩散方程采用诸如QUICK、SUD等高阶格式时， 进入P点离散方程的相邻点除了近邻点N、S、E、W 、T 和B 以外，还有NN、SS、EE、WW 、TT和BB六个远邻 点。其差分方程如下： 1D a  a  a  a  a  b P P E E W W EE EE WW WW 2D a  a   a   a   a   P P E E W W N N S S a   a   a   a  b EE EE WW WW NN NN SS SS 3D a  a   a   a   a   a   a   P P E E W W N N S S T T B B a   a   a   a   a   a  EE EE WW WW NN NN SS SS TT TT BB BB b 第五章 对流扩散问题———假扩散 由此可见，对一维而言，所得到的差分方程不再是可以 用TDMA直接求解的三对角矩阵方程，而变成一个五对 角阵方程。对二维或三维而言，逐线联立求解的方程组 也不再是可以用TDMA直接求解的三对角矩阵方程，也 变成一个五对角阵方程。 那么，针对这样一种五对角方程，通常我们如何来求解 它们呢？ 第五章 对流扩散问题———假扩散 A5.1 逐点求解 最容易想到和实现的就是逐点求解。实际上这种求解 方法与差分方程是采用几点格式构造的没有关系，具 有普适性。 Jacobi逐点迭代 n1 n n n n a   a  a  a  a  b P P E E W W EE EE WW WW Ga