MAC算法计算二维方腔顶盖流动 优先出版.pdfVIP

第15卷第6期 宜宾 学院 学 报 2015年6月 Vol.15,No.6 JournalofYibinUniversity June,2015 MAC算法计算二维方腔顶盖流动 李江飞，石兆东，段兴华，李岩芳，张 康，逯国强，陈颖超，任亚东 （承德石油高等专科学校热能工程系，河北承德067000） 摘 要：二维方腔流动是不可压缩黏性的典型流动，可以用来检验各种数值算法计算精度和可靠性，目前尚不能求得它的解 析解.基于Matlab编程，采用交错网格MAC算法求解二维方腔流动，计算采用控制容积积分法离散控制方程，对流项和扩散项 采用中心差分格式，得到流动达到稳定状态时各物理量的分布. 关键词：数值模拟；方腔流动；控制容积积分法；MAC算法；离散 中图分类号：TB126 文献标志码：A 文章编号：1671-5365(2015)06-0028-04 引用格式：李江飞，石兆东，段兴华，等.MAC算法计算二维方腔顶盖流动[J].宜宾学院学报,2015,15(6):28-31. LiJF,ShiZD,DuanXH,etal.CalculationofTwo-dimensionalCavityFlowBasedonMAC[J].JournalofYibinUniver⁃ sity,2015,15(6):28-31. CalculationofTwo-dimensionalCavityFlowBasedonMAC LIJiangfei,SHIZhaodong,DUANXinghua,LIYanfang,ZHANGKang,LUGuoqiang,CHENYingchao,RENYadong (DepartmentofThermalEngineering,ChengdePetroleumCollege,Chengde,Hebei067000,China) Abstract:Two-dimensionalsquarecavityflowisatypicalincompressibleviscousflow,whichcanbeusedtotestavari⁃ etyofnumericalalgorithmsforcomputationalaccuracyandreliability,yetitsanalyticalsolutionstillcannotbeachieved. BasedonMatlabprogramming,thestaggeredgridMACalgorithmwasusedtosolvetwo-dimensionalsquarecavityflow. Controlvolumeintegralmethodwasusedtodiscretethecontrolequationsforcalculation.Centraldifferenceschemewas appliedforconvectionanddiffusionterms.Intheend,theflowsimulationresultsofeachphysicalquantitydistributionin steadystatewereobtained. Keywords:numericalsimulation;cavityflow;controlvolumemethod;MAC;discrete 二维不可压缩黏性流体方腔流动顶盖拖动速度 x方向中垂线（x=H/2 ）上的无量纲速度U，y方向中 为u ，方腔的长度和高度均为H,流体密度为ρ、动 垂线（y=H/2 ）上的无量纲速度V，绘制出速度分布 top 曲线，并求出中垂线上 U、V 的平均值. 力粘度为μ. 边界条件：流动速度u、v采用无滑移 | | | | 边界条件，利用动量方程推导压力 [1] p 的