低马赫数流场数恒模拟预修正方洁研究.pdfVIP

低马赫数流场数值模拟的预修正方法研究 艾星钟诚文苏耀西 (两北工业大学114信箱西安7100q2) 摘要分析花数值求解低马赫数下非定常叫J^缔N-S方程的定常解时·蟛响隐』℃汁算 格止收敛速度的主要原因，并通过采取预修正的方法米由¨速定常解的收敛过程。义中算例 表l川，通过引入合适的修正矩阵，该方法效果叫娃井易于实现。 1 引 言 虽然在低马赫数下，可压缩流动的控制方程可近似丁不可压流动的控制方程，但人们 仍多采川可压缩流动的控制方程来数值模拟流场，这主要是基于以下原冈： (1)在人们所研究的很多低速流场问题中，同时存在着可压和不可压流动，例如低 速人攻角问题以及多段翼型的绕流问题。如果忽略流动的可压缩性，则不能准 确地反映流场的实际特征。 (2)K期以来，人们已研究并形成了一整套完善的针对可压缩流场的计算软件。因 此，人们希望能充分利用已有的软件，并使其具有更Jh阔地应川范同。 然而在低马赫数F，数值求解可压缩流动控制方程时会遇到圃难。对1j显式格式，我 仃j钼{道为了保证数值求解的稳定性，时间步欧的选取受到最人特征值的限制。对1：低马赫 数Jjq题，由于系统的特征值差异很大，因而使得有些波在一个时间步长过程中变化很小， 从而影响到整个系统定常解的收敛速度。对丁二AD[等隐式近似分解格式，在低马赫数情况 F、眭门。近似分解所产生的误差也同样会极大地影响刨算法的收敛性。 焉 为了解决上述问题，人们针对不同的数值计算格式提出了各种不同的解决方法，预修 止法就是其中一种被广泛应用的方法【l。l。预修1E法的实质就是通过给控制方程中的时间 翻 』mj：!jl乘一修正矩阵，来改变系统的特征值，从而达到提高收敛速度的目的。当我们只关心 、费 ‘ll：定常方程的定常解时，实践证明这种方法能火大加速定常解的收敛过程而并不改变系统 燕 t“ 定常解的数值。 本文将具体应用这种方法来计算绕多段翼型流场的可压缩NS方程的定常解。 2 控制方程 本文采用的控制方程是守恒形式的-二维1E定常可压缩N．S方程。在一般曲线坐标系中 维无姑纲化的N-S方程为： 229 寺+i+_=_(百+—_) a i{ jq Re、西 6q。 这’h哉J将采埘Briley-Mcdonaldi4)提出的ADI隐式方法来求解上式： (，+At罂)(，地譬)△Q；一Atg Db UIi 站巾 √：堡．口：堡 a2 趔 R；堕+堡一上f丝+～dFv) a{aHRc、j{卸； 本文对(3)式中的粘性项采州中心差分格式，对流项JI|I|采用TVD格式{“求解。 3 影响收敛性的因素 采州上述隐式格式进行数值计算时，人”』会发现在低马赫数条i，J=一F，格式的收敛性会 受剑极人地破坏a在分析其产生原因时，人们首先认为可能是由丁在低马赫数F，系统的 特“Ii值筹异很人所造成的a然而在数值实践中人们注意到，采用隐式格式造成上述问题卉勺 t婴原冈是由丁采朋近似分解所带来的误差161。 将(2)式展开可得： 翻 国 ，翻刃 (I+At琵+△。历+△‘2琵历)△Q=一△嘏 14) 该，℃足数值计算时真止求解的，而乘积项AB是影响计算收敛速度的主要冈素。 为了说明上述结论，将(1)式改写为非