Godunov方法在复杂流场数值模拟中的应用研究.pdfVIP

Godunov方法在复杂流场 数值模拟中的应用 纪楚群 李 骏 北京空气动力研究所 摘要 本文根据VanLeer提出的方法给出了求解三维非定常Euler方程和三维定常 Euler方程的二阶Godunov方法.方法要点为:在各体积离散单元中，用流动变量的分段 线性分布代替一阶方法中的分段常量分布，在确定线性分布的斜率中引人单调性限制条 件;沿离徽单元界面法向构造一维Riemann问题，由Riemann间断解确定单元侧面的通 量;用预测一修正两步法对离徽方程推进求解.用上述方法结合单区网格或多区重迭网格 给出了复杂外形的流场数值棋拟结果. 1 引言 各种外形复杂的飞行器在超声速飞行时流场将出现复杂的激波相交及反射现象.在各 类计算有激波流场的数值方法中，Godunov方法是一种非常有效的方法.该方法的特点 是将Riemann间断解直接引人流场数值模拟中，由干Riemann问题是可得到分析解的物 理问题。因此该方法是分析解和数值模拟成功的结合. 最早的Godunov格式为一阶精度的单调守恒格式.随后Godunov方法在提高计算 精度及向多维间题扩展方面得到进一步发展和应用R-wl 本文在前人工作基础上，给出了求解三维非定常Euler方程和三维定常Euler方程的 两类二阶Godunov方法.即时间相关法和空间推进法. 针对不同的飞行器外形特点，本文采用单区网格或多区网格.单区网格用代数预测一 抛物化求解Laplace方程修正的方法生成.多区网格用重迭网格方法生成. 用上述方法给出航天飞机、捆绑式运载火箭等外形的超声速复杂流场的数值模拟结 果.计算结果表明二阶Godunov方法能正确地模拟飞行器的复杂流场. 2数值计算方法 2.1Euler方程及有限体积离散 在直角坐标系 (x,y,z)中非定常Eul“方程在体积为△V的单元体上积分得: JJJ子()dVff‘U(，F,G),WdS一0 上述方程中S为单元体表面积，风为表面单位外法向向童.列向量W,U,F,G为 W二 (P,Pu,Pv,Pw,ph。一P)r 209 U=fpu，P+pu21puvIpuw,puh,飞r F二于pv，puv，p+pv2，pvw，pvh,下r G二{pw，puw,pvw，P+pwz，pwh,)r 上式中P.P,h。及 “，，，w分别为压力、密度、总熔及速度沿坐标z.y,z的分且. 其中总焙为h,为 P . 1 h。二一Y 一 宁 育 (u1+v2+,z) (2) v一 p ‘ 将(I)式写作下面的离散方程: ((W甲D△V:))i.}h+,一(WLaV).j.,+Hr-iiz，一Hf+1/3.ak 减 曰 用 +H，一、/z，一Hu+I/zk+W,4k:一H1.ak+，/: (3) 啪 通 其中万为单元体侧面通云，如Hi+1/2，为 祖 与 H)十:/:，，二△t(USy,+FS:二+GS,y)丁+，/:，+(PLV)，十1/alk (4) 书 上式中，S，:，S:二，S=，分别为该侧面积在坐标平面v一:，:一二，二一，上的投影 值.如计算网格取不随时间变化的固定网格，CVOA，其它侧面通蚤H的形式类似.离