稳态三维不可压缩流体非等温有限元模拟.pdfVIP

稳态三维不可压缩流体非等温有限元模拟——秦升学 赵 国群 栾贻国 稳态三维不可压缩流体非等温有限元模拟 秦升学 赵 国群 栾贻国 山东大学，济南，250061 摘要 ：针对幂律型流体 ，建立了稳态三维粘性不可压缩流体非等温流动的有限元模型。采 用罚有限元法将连续性方程 引入到动量方程，获得 了速度场分布。能量方程采用SUPG(stre— amlineupwind／Petrov--Galerkinformulation)方法构造非对称权 函数 ，克服 了对流项 占优 时 数值解的失真现 象，获得 了稳定的温度场分布 。以幂律型流体本构关 系计算了剪切速率对粘 度的影响 ，采用Arrhenius模型计算 了温度对粘度的影响。以幂律 型流体矩形截面稳态收敛 流动为例，计算 了区域 内流体的速度和温度分布规律。计算结果表明，经过有限次迭代可获得 合理速度场分布 ，在能量方程 出现对流项 占优的情况下，即便采用较粗的网格划分 ，仍能获得 稳定 的温度场分布 。 关键词 ：数值模拟 ；SUPG方法 ；粘性流体；有限元 中图分类号 ：0357．1 文章编号：1004—132X(2005)21—1957一O5 FiniteElem entAnalysisofSteadyThreeDimensionalNon-isothermalIncompressibleFluid QinShengxue ZhaoGuoqun IuanYiguo ShandongUniversity，Jinan，250061 Abstract：The modeling forthree dimensionalincompressible viscousnon— isothermalsteady flow waspresentedbyusingpenaltyfiniteelementmethod．Thevelocityfieldwasestablishedbyin— troducingthecontinuityequationintomomentum equation．ByadoptingSUPG(StreamlineUpwind／ PetrovGalerkinformulation)forselectionoftheunsymmetricalweightfunctionforenergyequation， thedistortionofthenumericalresultswhen theconvective term iSdom inatewasovercome．and the stabletemperaturedistributionwasobtained．Thepaperalsocoupledtheviscositywithshearrateby power— law model，and coupled theviscositywithtemperatureby Arrheniusmode1．A power— law fluidflowingprocessthoughthree——dimensionalrectangularconvergentchannelwasselected forthe den1onstrationofthe method．Thevelocity and temperaturedistributionsareobtained．Theresults show thatthereasonable velocity and temperaturefieldscan be obtained within severaliterations． W hentheconvectiveterm isdominantinenergyequation，thereasonabletemperaturefieldcanbealso obtainedeven withthecoarsemesh． Keywords：numericalsimulation；SUPG formulation；ViSCOUSfluid；finiteelement 在一定的困难 。 0 引言