有限元非结构网格中的AUSM算法研究.pdfVIP

有限元非结构网格中的AUSM算‘法 伍贻兆 王江锋 (南京航空航天大学六系 坏编:210016) 摘要 本文将近年来发展起来的AUSM (AdrectionUpstreamSplittingMethod)算法推广到有限元非结杜网 格中并成功地进行丁三维欧拉方程的数厦计算。利用外插获得空间二阶情度.应用Albabda型限制器以抑制 nAW97N荡，井构.r4J线化逆风二阶情度隐式格式·提高了计算效率此格式具育FVS(FlneVectorSplittiag) 咯式的简单性及FDS(FlaxDeferenceSplitting)格式的高阶情度 最后给出了NACAOM2翼型绕流、平板畸流 及后向台阶喷流流场的计算结果〕 关键词 AUSM ‘育限元 非结构网格 欧拉方程 1.引言 -}FCFD数值分析中.数值求解格式的效率及情度一直受到气们的普遍关注 本文中 将 通量的对流部分与压强部分分离·用适当的方法分别处理线性场及非线性场以提高激波N捉 能力。应用通量外插以获得空间离散的二阶精度。应用Alabda型限制器以抑制数值解的振荡。 并将此法发展成为隐式格式以提高求解效率 非结构网格采用Delaam q法生成 最后给出 了应用木文所述方法在非结构网格上对NACA0012翼型跨声速定常绕流流场进行数值计箕的 结果，并与FVS中VanLeer方法的结果进行了讨比，结果符合的较好.另外对平板喷流乏后 向台阶喷流流场做了数值计算并给出了计算结果 Z理论分析.(=‘1 以一维流的初值问题为例，控制方程为 擎+F二。.，。、一二:::十二 礴I) {V(x,0)=城(x) 式 中 GTR=W,prr,el:=Vrr,pir+p,prrHf 。为单位体积总能.H为总恰 ·ph密度·P为压强。 对理想气体有: ]，二，一‘，)(e一尽，:‘“二{e+p/P (2) 209 摘方N(1)进行有限体积离散后 离散戈式中包含通量F在域交界面 t的边界积分，记F 在域交界而上的值为f 有: 3,;=f.(IIU一，) (3) (3)式为Lipschitz连续函数并满足相容性条件: f(1,.(1，一)=F(UJ) (4) 将通量F写为卢二F+尸的形式，有: F’二AfaP(二MaV,pn,pHfP二二10,p.01 (5) 其中.Al为马赫数 :‘为音速口 交界面七的通量F0及只一可由下式计or: 尽、二m:a了一，叭 P一，二O【,p一;0] (6) 以下以。/代‘替 I‘、告 2.1 ，的定义 应月简单的逆风方法 巾「 电二)巾 (7) 2m,、p,的定义 定义m 为: f九f，L/ 二口一长认一; r1Sa) Aft二iffM )·A4;_,=.切一(A-1.，) i8b) 可 取 叼 (弄f)为: 告(A了二}ml) i盯l?1 刀