原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
解线性方程组的TDMA算法在CFD中的应用----科学计算选
解线性方程组的TDMA算法在CFD中的应用 -- -- 科学计算选讲结课报告 报告人：刘波 学号： 1009207047 专业： 化学工程 研究方向：乙烯裂解炉的CFD模拟 流体力学基本控制方程 流体力学基本控制方程 方程具体形式见下表： 基本方程的通用形式 将各方程中的物理量如速度、内能等用广义物理量ψ来表示，扩散系数μ、k等用广义扩散系数Γ表示。方程 1-1~1-5可以表示成如下有用的形式： 稳态扩散问题的有限体积法 一维稳态扩散问题的有限体积法 稳态扩散问题的有限体积法 我们对方程1-9在控制体积CV上进行积分： 对于导数项我们采用中心差分近似： 源项采用线性化处理： 并令已知项： 稳态扩散问题的有限体积法 这样，离散方程可以简化为如下形式： 由于1-10式对于每个网格节点都成立，所以1-10事实上是 一个方程组，赋予恰当的边界条件方程组便具有如下形式： 解三对角方程组的TDMA算法 CFD中我们常常会遇到形如式1-11的三对角方程组，针对该方程组最常用的方法便是计算机资源占用较小，效率较高的TDMA算法，下面就给出这种算法的推导。 上述1-11中， 是已知的边界值。各方程的通式为： 上述各方程可以重写为： 解三对角方程组的TDMA算法 解三对角方程组的TDMA算法 解三对角方程组的TDMA算法 结束 * * 以能量守恒定律、质量守恒定律以及牛顿第二定律为基础，结合扩散定律（质量扩散、热量传导和动量传递）我们可以推出流体力学基本输运方程。 （1-7） 称式1-7为对流扩散方程，本报告就以对流扩散方程中最简单的形式即略去对流项 和瞬态项 的稳态扩散方程为例，来介绍TDMA算法在扩散方程的数值解法中的应用 稳态扩散方程： 对于一维问题： (1-8) (1-9) (1-10) (1-11) 式1-11为一个三对角方程组 这些方程可用向前消去和向后代入法求解, 由1-12b解出ψ2代入1-12a，得到： 1-12a 1-12b 1-12c 1-13a 记， 式1-13a可以化为， 1-13b 记， 式1-13b可以化为， 1-13c 式1-13c可以从1-12c中消去ψ3，如此进行下去到最后一个方程便构成了向前消去 当消去过程进行到（n+1）处时，由于ψn+1已知，式 1-13c便可以求解得出ψn，依次反推回去便可得到方程组的解ψj。 以上推导仅介绍了三对角方程组的TDMA算法，对于二维问题和三维问题，方程组将不在是三对角方程组，而是五对角和七对角方程组。此时，TDMA算法仍然有效，只是加入了额外的迭代过程。限于篇幅和时间就一一不在介绍。 式1-13c的递归形式为：
文档评论（0）