强涡-激波干扰及其声场的数值模拟.pdfVIP

2001 第六届全国流体力学会议论文集 上海 强涡一激波干扰及声场的数值模拟 、————一 一 …～—、’h，一 马延文‘ 傅德薰” ’LHD．“LNM，中国科学院力学研究所．北京100080) 摘要气动声学是气动力学与声学之问的交叉学科，着重研究流体流动及其与物体作用导致 噪声的产生机理．对它的研究有重要应用背景．本文由两个部分组成．为了更好地捕捉激波 和多尺度物理量发展了带群速度控制的超紧致格式，利用这一方法数值模拟了强涡+撒泼干扰 问题，绐出了四极于声场的发展及细致的流动结构． 关键词 超紧致格式，涡一激波干扰．(声场的数值模蛾 1}、 1 带群速度控制的超紧致格式 计算气动声学所面对的是求解非定常多尺度问题。其特点是一些流动参数的量级和尺度 差异大．声能和主流能量差异大、频带宽、传播距离远等。由于条件限制在对区域进行有限 截断时给边界条件的处理也带来困难。这些特点要求计算方法有着小的色散和耗散误差，高 精度方法是可取的方法之一。 利用文献[1】中发展的超紧致格式咀提高对小振幅高频波的模拟能力，对这一方法加以群 速度控制咀提高对激波的捕捉能力。利用完全可压Navier-Stokes方程来模拟涡一激波干扰问 题。下面对数值方法作一简单描述，考虑咀下之模型方程及对应的半离散格式 ㈨ 乱 融 詈+鼍=上Re窘缸2 丝+生一；上竺一 (2) 西△x RP△r。 这里g≯’／Ax‘为导数船‘，缸‘的差分逼近式．【2)式中的f“’可由如下之差分方程中求得‘‘1 一亡—L乃一I+(工+一)一一÷．L日“=占：乃el (3) 这里F．和e、为N维矢量．一和工为NxN矩阵 F=[f‘“，f咖，．，厂。““】7 土 !… ! l!’3l’ ’【2N一11【 m A=。，击，…，丽1 扣 … ‘一，‘“，’”， 一㈨h||㈨一 m，仉 习 ooI．．-， 去 在计算中取N：3，对一阶导数逼近精度为六阶。这一方法有利于捕捉小振幅高频波物理量· 经线性分析可知，如此构造的格式属慢型(sL聊】，此时数值解中的高频波的群速度将慢 378 于对应的物理波群速度。用其计算激波时将在激波后产生非物理振荡。为消除非物理振荡， 需对群速度加以控制f2J，使得格式在渡后为快型(FST)，波前为SLW型。可以证明．耗散比 拟修正Ⅲ可以控制群速度。文中对(2)式的右端加了经过控制的7阶导数项，使得耗散比拟系 数为正，且不降低原格式的精度。7阶导数的计算利用了由(3)式中求得的，“’。即在计算中 所采用的半离散格式为 2 J、’ V’㈣ 却 鲁+竖Ax=上Re譬A,r+喇 式中D，0)为D(u)的离散形式，而 D(“)=出6昙一罢