激波钝体相互作用的数值模拟1.PDFVIP

激波钝体相互作用的数值模拟1 许 东 洋，吴 锤 结2 （工业装备结构分析国家重点实验室，大连理工大学航空航天学院，大连 116024 ） 摘要：本文采用五阶加权本质无振荡（Weighted Essentially Non-oscillatory, WENO ）格式，通过数值求解二维非定 常 Naivier-Stokes 方程，模拟了马赫数为 1.2，雷诺数为1.0 ×104 下绕钝体的可压缩流动，系统地分析了流动中流 场的演化过程中的若干基本现象，如复杂的多激波系结构、激波与边界层及激波与旋涡的相互作用、尾流中的发 卡状相干结构等，以及激波、旋涡、边界层相互作用产生声场的机理。 关键词：超声速流动；激波；尾流；声场 引 言 激波与钝体的相互作用是一个基本的空气动力学问题，它不仅与一些基础性研究问题相关，如激波绕 尖角后旋涡的形成、激波-旋涡的相互作用、激波与边界层干扰问题等，而且广泛应用于飞行器的设计与防 护工程。通常与激波相互作用问题的现场实验（或飞行实验）是相当昂贵的、复杂的，并且很难得到实验 数据，有时仅能有效分析模型上一些测量点的数据。另一种研究激波相互作用现象的方法就是通过数值模 拟与激波管实验，这种方法在相对较低的成本下即可给出流场的的完整描述，同时为工程设计提供了有价 值的数据。 本文采用五阶 WENO 格式，数值求解了二维非定常 Naiver-Stokes 方程，对激波与钝体相互作用问题 进行了数值模拟。在本文中，计算域为 ［－30≦x≦20，－25≦y≦25］，绕长为 3、宽为 1.4 的矩形的超声速 流动。 1 数值方法 1.1 控制方程 采用三维非定常可压流的 Naiver-Stokes 方程作为控制方程： ∂U ∂F ∂G 1 ⎛∂F ∂G ⎞ + + ⎜⎜ v + v ⎟⎟ (1) ∂t ∂x ∂y Re ⎝ ∂x ∂y ⎠ 其中， ρ ρu ρv ⎡ 0 ⎤ ⎡ 0 ⎤ ⎡ ⎤ ⎡ ⎤ ⎡ ⎤ ⎢ ⎥ ⎢ 2 ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ρu p +ρu ρuv ⎢ τxx ⎥ ⎢ τxy ⎥ U ⎢ ⎥，U ⎢ ⎥，U ⎢ ⎥，F ，G . ⎢ ⎥ ⎢ ⎥ ⎢ 2 ⎥ v ⎢ τ ⎥ v ⎢ τ ⎥ ρv ρuv p +ρv xy yy ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ + + + + u e p v e p uτ vτ q uτ vτ q e ( + ) ( + ) ⎢ xx xy x ⎥ ⎢ xy yy y ⎥ ⎣ ⎦ ⎣ ⎦ ⎣ ⎦ ⎣ ⎦ ⎣ ⎦ 1 国家自然科学基金资助项目（） 2 cjwudu