解线性方程组的TDMA算法在CFD中的应用 ---科学计算选讲结课报告.PPT

小组成员： 邢令利 马奎 学号： 1014207129 1013207126 专业： 工业催化 */16 专业背景 工业化学反应：石油炼制、有机合成、聚合物生产、无机化工、制药、三废治理等 绝大多数工业化学反应是在催化剂作用下进行 研究方向：催化反应的CFD模拟 */16 CFD （计算流体动力学） CFD，软件(Computational Fluid Dynamics），即计算流体动力学, 简称CFD。CFD是近代流体力学，数值数学和计算机科学结合的产物，它以电子计算机为工具，应用各种离散化的数学方法，对流体力学的各类问题进行数值实验、计算机模拟和分析研究，以解决各种实际问题。 计算流体力学和相关的计算传热学，计算燃烧学的原理是用数值方法求解非线性联立的质量、能量、组分、动量和自定义的标量的微分方程组，求解结果能预报流动、传热、传质、燃烧等过程的细节，并成为过程装置优化和放大定量设计的有力工具。 */16 流体力学基本控制方程 以能量守恒定律、质量守恒定律以及牛顿第二定律为基础，结合扩散定律（质量扩散、热量传导和动量传递），可以推出流体力学基本输运方程。  */16 方程具体形式  */16 基本方程的通用形式 将各方程中的物理量如速度、内能等用广义物理量ψ来表示，扩散系数μ、k等用广义扩散系数Γ表示。 方程1-1~1-5可以表示成如下有用的形式：  （1-7） 称式1-7为对流扩散方程，本报告就以对流扩散方程中最简单的形式即略去对流项 和瞬态项 的稳态扩散方程为例，来介绍TDMA算法在扩散方程的数值解法中的应用 */16 稳态扩散问题的有限体积法 一维稳态扩散问题的有限体积法  稳态扩散方程： (1-8) 对于一维问题： (1-9) */16 稳态扩散问题的有限体积法  我们对方程1-9在控制体积CV上进行积分： 对于导数项我们采用中心差分近似： 采用线性化处理： 并令已知项： */16 稳态扩散问题的有限体积法 这样，离散方程可以简化为如下形式： 由于1-10式对于每个网格节点都成立，所以1-10事实上是 一个方程组，赋予恰当的边界条件方程组便具有如下形式： (1-10) (1-11) 式1-11为一个三对角方程组 */16 TDMA算法 三对角矩阵的常用解法为Thomas Algorithm，又叫The tridiagonal matrix algorithm (TDMA)。它是一种基于高斯消元法的算法， 分为两个阶段：向前消元和回代。 CFD中我们常常会遇到形如式1-11的三对角方程组，针对该方程组最常用的方法便是计算机资源占用较小，效率较高的TDMA算法  */16 TDMA算法一般公式  */16 上述1-11中， 是已知的边界值。各方程的通式为： 上述各方程可以重写为：  TDMA算法在CFD中的应用 */16 1-12a 1-12b 1-12c  这些方程可用向前消去和向后代入法求解, 由1-12b解出ψ2代入1-12a，得到： 1-13a 记， 式1-13a可以化为， 1-13b TDMA算法在CFD中的应用 */16  记， 式1-13b可以化为， 1-13c 式1-13c可以从1-12c中消去ψ3，如此进行下去到最后一个方程便构成了向前消去 式1-13c的递推形式为： TDMA算法在CFD中的应用 */16  当消去过程进行到（n+1）处时，由于ψn+1已知，式 1-13c便可以求解得出ψn，依次反推回去便可得到方程组的解ψj。 以上推导仅介绍了三对角方程组的TDMA算法，对于二维问题和三维问题，方程组将不在是三对角方程组，而是五对角和七对角方程组。此时，TDMA算法仍然有效，只是加入了额外的迭代过程。限于篇幅和时间就不再介绍。 TDMA算法在CFD中的应用 */16 */16 LOGO 1/21 1/21 LOGO