资源描述：
《求解Euler方程的区域分解方法与并行算法》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、计 算 物 理第17卷第4期Vol.17,No.42000年7月CHINESEJOURNALOFCOMPUTATIONALPHYSICSJul.,2000[文章编号]10012246X(2000)0420360207求解Euler方程的区域分解方法与并行算法123吕晓斌, 兰黔章, 朱自强(11中国科学院力学研究所,北京 100080;21上海交通大学动力系,上海 200030;31北京航空航天大学流体所,北京 100083)[摘 要]将复杂形状区域划分成多块子区域,研究发展了一种多块区域之间迎风守恒型的内边界耦合方法,实现相邻子区域解的光滑过渡,使多区耦合
2、得到总体流场的数值解。对二维翼型跨音速流动和圆弧形隆起物超音速流动等进行了分区数值计算,并将计算结果与单区计算结果和实验结果作了比较。并行分区计算引入“先进先出”的同步控制等待机制,实现了高效率并行计算,还分析了影响并行效率的主要因素。[关键词]计算流体力学;区域分解方法;并行处理;Euler方程[中图分类号]O35;V211[文献标识码]A0 引言复杂外形无粘性可压缩流场的数值计算,主要涉及到偏微分方程组的求解。该类问题的特征波速是有限的,如果把复杂几何区域分解成若干块子区域,子区域之间通过一定的耦合条件进行信息交换,以实现相邻子区域解的光滑过渡,复杂区域总体流
3、场可实现并行求解。因而[1]分区方法是一个有效地并行处理复杂区域流场的数值方法基础。本文结合矢通量分裂数值计算方法和信息传播的特征理论,采用完全对接网格,研究发展了一种迎风型内边界通量守恒的耦合条件,应用于二维跨音速翼型绕流、圆弧形隆起物超音速流动激波反射现象的分区并行计算。1 控制方程和数值方法在直角坐标系(x,y)中,二维非定常Euler方程组的微分形式为5U5F(U)5G(U)++=0(1)5t5x5y式中,U是守恒矢变量,F(U)、G(U)是矢通量ρρuρvρuρu2+pρuvU=,F(U)=,G(U)=(2)ρv2ρuvρv+pE(E+p)u(E+p)v
4、[收稿日期]1999204220[基金项目]国家自然科学基金资助项目[作者简介]吕晓斌(1965～),男,山东,博士后,从事计算空气动力学并行计算方面的研究.第4期吕晓斌等:求解Euler方程的区域分解方法与并行算法361理想气体状态方程为122p=(γ-1)[E-ρ(u+v)](3)2其中,ρ为密度,u、v为速度分量,p为压力,γ为比热比,空气γ=114,E为单位体积的总能。用有限体积方法在每一个网格单元控制体上进行积分,计算变量定义在网格中心点上,在边界上计算通量,得到一组空间离散方程Vi,jdUi,j=-dtŽFi+1/2,jSi+1/2,j-ŽFi-1/2
5、,jSi-1/2,j+ŽFi,j+1/2Si,j+1/2-ŽFi,j-1/2Si,j-1/2(4)式中ŽF=nxF(U)+nyG(U)(5)其中,S是控制体Vi,j的边界面,n为边界面的外法向单位矢量。如果网格是足够光滑的,边界物理量由二阶插值公式计算,则方程式(4)是二阶的。[3]本文按VanLeer矢通量分裂方法计算控制体边界面对流矢通量。对二维流动情况,根[4]据边界面的法向流动速度进行矢通量分裂。公式(4)中ŽFi+1/2,j,ŽFi,j+1/2等不能直接应用±±VanLeer通量分裂公式得到ŽFi+1/2,j,ŽFi,j+1/2等,必须先用一个变换T得到
6、局部坐标系中法-1向通量,然后按一维分裂方法对法向通量作分裂,分裂后再用一个反变换T,将通量变换回±来才能得到ŽFi+1/2,j等。亚音速流动(
7、Mn
8、<1)时,边界面法向矢通量分裂公式为ρc2±(1±Mn)4±(γ-1)MHn=h1c[n±2]/γ(6)±h1cMt±(E+p)/ρh1±±(ρc/4)(1±M)2式中h1表示第一项,即h1=±n,Mn为边界面的法向马赫数,Mt表示与界面法向垂直的马赫数分量。超音速流动时,
9、Mn
10、E1,此时法向特征速度(qn+c,qn,qn-c)全为正或全为负。由信息传播理论,边界面中点物理量全部由上游一侧插值计算,固壁物面边界上
11、满足壁面无穿透条件,即物面法向速度等于零。远场根据Riemann不变量分析,采用无[5]反射边界条件。时间积分方法采用分区隐式因式分解(AF)方法。2 内边界耦合条件211 区域分解与内边界耦合条件区域分解方法将总体区域分成几个子区,然后在每个子区域网格上对Euler方程分别求解。各个子区域的解通过内边界耦图1(a)重叠网格 (b)对接网格合条件光滑过渡。Fig.1(a)overlappinggrid(b)patchedgrid分区方式有重叠和对接两种。本文讨论对接方法,各子区域无重叠。362计 算 物 理第17卷一般设A,B为重叠分区,重叠区域为O,
12、如图1(a
