不可压缩流两重稳定有限体积算法应用研究.pdfVIP

ComputerEngineeringandApplications计算机工程与应用 2015，51(18) 255 不可压缩流两重稳定有限体积算法应用研究 杨建宏 YANG Jianhong 宝鸡文理学院 数学系，陕西 宝鸡 721013 DepartmentofMathematics，BaojiUniversityofArtsandSciences，Bali，Shaanxi721013，China YANG Jianhong．Two—levelstabilized finitevolumemethod forincompressibleflow．ComputerEngineeringand Applications，2015，51(18)：255—260． Abstract：Inthispaper,two ‘levelstabilizedfinitevolum emethodsareconsideredwhicharebasedon localGausslnte。 graltechniqueandtwo—levelgridalgorithm fortheincompressibleflow．Theerroranalysisshowsthatthetwo—levelstabi- lizedfinitevolumemethodsprovideanapproximatesolutionwiththeconvergencerateofthesameorderastheusualsta— bilizedfinitevolum esolutionsolvingtheincompressibleflow problemsonafinegridofrarelatedchoiceofmeshwidths． Theperformanceofthreekindsiterativeschemeoftwo-—levelstabilizedmethodsarecomparedinefficiencyandprecision aspectsbyaseriesofnumericalexperiments．Itdiscoversthatthesimpleschemeisbetterthantwoothersonaccuracyand efficiency．ThereisthepoornumericalaccuracyfortheNewtonscheme，butthePicardschemeismoresuitabletoincom— pressibleflow withlow viscositycoeffi cient． Keywords：incompressibleflow；finitevolumemethods；two—levelstabilizedmethods；Navier-Stokesequations 摘 要：通过将局部高斯积分稳定化方法和两重网格算法思想紧密结合，提出了粘性不可压缩流体的两重稳定有限 体积算法。将该算法的三种迭代格式进行了效率的分析比较。理论分析和数值 实验发现：当粗、细网格尺度比例选 择适当时，两重算法与传统算法具有相 同精度解的同时，效率大大提高；对不同格式的两重有限体积算法进行比较 分析发现：Simple格式计算效率最高，Picard~式次之 ，Newton格式较低。 关键词：不可压缩流体；有限体积算法；两重稳定方法；Navier．Stokes方程 文献标志码：A 中图分类号：TP39 doi：10．3778~．issn．1002．8331．1309．0292 l 引言 性和非线性特征。不可压缩性要求求解它的有限元对 有限体积方法 (FVM)也叫作控制体积方法或一阶 必须满足LBB条件，即要求对速度逼近的多项式的次 广义差分方法 ，它是求解偏微分方程的一种重要数值工 数高于对压力逼近的多项式的次数，这在实际计算中需 具。该方法已经广泛应用于科学计算和工程领域，尤其 要用到两套不同的计算网格 ，显然计算难度和工作量都 大大增加了。而低次等阶有限元PI—P 1或Q 是在流体动力学 ，热量传导能量传递 ，油藏工程等方面 _ 一 一 1Q·一1