(五)代数方程的求解-LSEC.ppt

(五) 代数方程的求解 5.1 代数方程系统 5.2 直接法 5.3 主要迭代法 5.4 其他迭代方法 5.1 代数方程系统 有限差分(体积)离散格式提供一个网格点(单元）的代数方程, 以线性代数方程为例： P点和周围邻居点构成计算模板(比差分基架还大） 计算模板(计算分子;解元SE) 5.1 代数方程系统: 计算模板 5.1 代数方程系统: 整体方程系统 流场中每一点都有一个方程(小组), 整个计算域就有一个大型稀疏方程系统 5.1 代数方程系统: 系数矩阵的存储 只存储非零的对角元素 2维5点格式： 5 Ni *Nj 3维7点格式： 7 Ni *Nj*Nk Al,l-Nj=W Al,l-1 =S Al,l =P Al,l+1 =N Al,l+Nj=E 5.2 直接法 5.2.1 Gauss elimination 5.2.2 LU decomposition 5.2.3 Tridiagonal system 5.2.4 Cyclic reduction 5.2.1 Gauss Elimination 5.2.2 LU decomposition 5.2.3 Tridiagonal system (TDMA) 5.2.3 Tridiagonal system:块三对角方程组 5.2.3 Tridiagonal system (cont) 计