一种新的无矩阵运算的高效差分格式.pdfVIP

中国工程热物理学会 982050 热机气动热力学 学术会议 一种瓤的无矩阵运算的高效差分格式 渤 祁明旭刘斌丰镇平 (讴安交通大学时轮机械研究所，西安。710049) ， 黜 [摘要]本文提出了一种新的显、隐两步差分格式。该格式在构造显式和隐式差分步时， 采用即时修正的思想，充分利用了n+l时刻点的信息值，从而可以加速边界信息在计算 区域的传播速度；同时，在酪式步计算中，结合谱半径恩想，简化了格式，完全避免了 矩阵运算。通过对一维缩放喷管算例的计算表明，新格式可以减少3／4’1／2的计算时间， 是一种高效的差分格式。 关键词：差分格式矩阵运算 ～．概述 随着计算机技术的高速发展，全三维流场的数值计算已成为一种反映某种流动工况 流动特性的非常有效、可靠的手段。然而，通过时间推进法的金三维数值计算来获得计 算流场区域的收敛信息需要耗费大量的计算时间和资源。而在较多的情况下，进行全三 维计算的目的是进行设计或者优化(如透平机械中的时型正设计问题)，这就要求计算 的时间越短越好。进行三维数值模拟的计算时间及计算量在计算条件相同的情况下(如 网格数、出入口参数等)取决于计算所采用的差分格式以及计算方法。近几年来，已经 发展出来许多商效的差分格式【‘】f2Jf”。所有这些差分格式基本上都是采用较大的时间步 以使得边界信息快速传播至整个流场。本文从差分格式上进行了改进．引入了即时修正 和谱半径的思想，在相同的时间步长的条件下，能够加速信息在流场内的传播，并完全 避免了求解隐式步计算时的矩阵运算，大大减小了计算量，从而获得了一种高效的差分 格式。本文给出了新型差分格式的构造，并以准～维的缩放嚷嚷为算铡进行了考核，证 明了本文所发展的格式是正确的、有效的。 二。差分格式的构造 虑 ： 考 矢 量 形 式 的坂 曲型 勰晒～a 十 驴iIIO (1) ～ l 引入Jac。bian系数矩阵：爿=丽OF，式(1)可以写成 型十A罂。0 (2) 一十A_2 执 a‘ 对U”1进行Taylor展开有： (3) u…刈心筹+譬掣O1十o(△t)3m Z 。 匆 ． 6U“=U““～U“ 令 (4) 8Fn“：Fn”一F“ 、 嘞 式(3)的符号右端第二项用米构造显示部分；将式(1)代入可得 AU础型础f一竺1 (5)