非结构混合网格鲁棒自适应技术

数值格式、湍流模型和计算网格是影响CFD数值模拟精度的3个主要因素。结合流场信息的网格自适应技术具备动态优化计算网格的能力，被NASA列为未来CFD发展的一项关键技术。本文针对非结构混合网格，发展了网格单元分布优化、表面网格几何投影和空间网格协调匹配3项关键技术，建立了高鲁棒性几何保真的网格自适应系统。首先，为了提高自适应方法的鲁棒性和通用性，发展了基于标准面网格的多面体网格单元分布优化方法。其次，发展了仅依赖表面网格信息的局部曲面重构技术，采用参数点映射方法实现了新增表面网格点的几何投影，消除了自适应系统对几何CAD系统的依赖。再次，采用改进的距离函数方法实现了空间网格与投影后表面网格的快速匹配。最后，结合基于流场特征的自适应探测器，采用二阶格式的有限体积方法，开展了30P30N三段翼绕流和三角翼大迎角绕流的网格自适应数值模拟。结果表明，通过网格自适应对网格单元的分布进行优化后，流场求解的收敛性和模拟精度都得到了显著提高。     

压气机转子三维流场计算中叶尖间隙网格的影响

对一带叶尖间隙的低速轴流压气机转子三维流场进行了数值模拟,通过与实验结果对比,分析了叶尖网格匹配、网格质量对压气机转子三维流场计算的影响,发现网格正交性和网格疏密过渡会影响二阶计算精度的收敛,研究还显示压升特性线不能准确的反映网格数目对计算结果的影响。      

基于二维结构化网格的间断Galerkin方法研究

在二维结构化网格上,对求解Euler方程的间断Galerkin方法进行了研究。应用二阶、三阶精度的问断Galerkin方法对绕NACA0012翼型的跨音速无粘流动进行了数值模拟。在整体逼近精度基本保持不变的前提下,采用moment限制器以消除数值解在激波附近的伪振荡。数值结果表明,间断Galerkin方法能够很好地模拟流场,准确捕捉激波;moment限制器成功抑制了激波附近的伪振荡;高阶格式具有比低阶格式更强的激波捕捉能力。     

高速运动弹丸诱导的瞬态流场数值研究

采用基于任意拉格朗日欧拉(ALE)方程的动网格及嵌入网格法,结合二阶精度Roe格式和沉浸边界法,对弹丸由高压气体驱动从静止状态加速至超声速,射出膛口到完全飞离初始流场的整个过程进行了数值模拟.根据数值结果,详细讨论了初始流场、火药燃气流场的形成与发展以及与弹丸的耦合和相互作用过程.      

基于非结构网格的多段翼型流场数值模拟

将非结构网格生成技术应用于多段翼型流场数值模拟，针对多连通域可以快捷、迅速地生成非结构网格。通过控制背景网格节点的网格生成参数，进行网格的局部加密，提高多段翼型流场分析的精度，应用有限体积法求解欧拉方程，应用各种加速收敛措施提高流场求解速度，并将数值计算结果与风洞试验结果进行了对比分析，证明该方法具有型号设计应用价值。     

基于超网格的重叠网格守恒插值方法

保证重叠网格边界数据插值的守恒性是计算流体力学面临的一大难题。基于超网格方法与格心格式有限体积法,发展了一种新型混合重叠网格守恒插值方法。详细研究了用网格切割技术在重叠网格边界构造局部超网格,用网格求交算法合理地扩大贡献单元模板,建立了一种适用于任意网格类型的隐式并行重叠网格守恒插值方法,以超网格为媒介可实现重叠网格二阶精度的守恒插值。数值结果表明,本文方法对二阶分布的流场变量具有严格的守恒性,相比三线性插值方法和逆向距离权插值方法,本文方法减小了数值误差,提高了重叠网格边界的插值精度,加快了计算收敛速度,改善了重叠区域网格尺度相差较大时流场的光滑性和连续性。     

入口温度剖面对喷管流场结构的影响

应用质量平均的Navier-Stokes方程和B－L代数湍流模型，对超燃冲压发动机尾喷管的流场进行了数值模拟研究，在计算过程中，对方程中的对流项采用了空间为二阶精度的TVD格式，扩散项则采用了二阶中心差分离散，通过数值模拟，对比研究了温度非均匀性，三维效应对尾喷管的流场结构的影响。     

S型进气道的数值模拟研究

 本工作改进了Denton的有限体积的时间推进法,将这种对时间一阶精度对空间二阶精度的方法改进为对时间二阶精度对空间二阶精度的计算方法,并应用于S型进气道的内流及具有进口外流的进气道流场的数值模拟,进行了实例计算,取得了合理的计算结果。     

基于无网格算法的化学非平衡流数值模拟

为了验证AUSMPW+格式在无网格算法的精度等以及无网格算法能否用于复杂化学非平衡流数值模拟中,首先将AUSMPW+格式推广到无网格算法中,给出了AUSMPW+在无网格算法中的具体形式。其次,将无网格算法应用于带化学反应的多组分气体Euler方程的求解,给出无网格条件下求解化学非平衡流的具体过程。最后,对NACA0012跨声速绕流场、高超声速弹丸绕流场、氢气/氧气诱导燃烧流场、球锥高超声速绕流场、楔体斜爆轰流场进行了数值模拟。数值模拟结果与相关文献计算结果和实验结果吻合较好,表明了AUSMPW+格式在无网格算法中表现较高精度。同时,采用无网格算法能较好地模拟复杂化学非平衡流场,正确分辨复杂的物理现象,为化学非平衡流数值模拟提供一种新的算法,拓宽了无网格算法的应用范围。     

冲压增程弹丸进气道特性分析

采用块结构网格与二阶精度流场分区求解技术，对固体火箭冲压发动机增程弹丸超声速进气道特性进行了深入研究。通过数值模拟得到了对应于不同来流马赫数和攻角情况下，临界工况时，超声速进气道内外粘性流场复杂的波系结构，详细分析了来流马赫数和攻角对进气道性能的影响。结果显示，随着来流马赫数的增大，总压恢复系数显著降低，流量系数增大，同时随着来流攻角的增大，总压恢复系数及流量系数逐渐降低，而流场畸变指数则明显增大。     

一种包含运动边界的高精度流场数值计算方法

为了准确快速地模拟运动边界的流场,提出一种反馈力源形式的包含运动边界的非定常流场数值计算方法。该方法采用完全正交的网格,以反馈力源作用点的运动来模拟边界的运动。采用物理量及其各阶导数在边界两侧的突跃修正中心差分格式,使之达到二阶精度,以此离散求解二维不可压Navier-Stokes方程。并且提出了与运动边界相适应的反馈力源构造方法及对边界上速度进行插值的方法。基于此数值计算方法,对低雷诺数的圆柱绕流、静止流体中的振荡圆柱以及昆虫振翅运动的二维非定常流场进行了数值计算,计算结果与以往的数值及实验结果非常吻合,表明本文方法与Peskin的浸入式边界方法在处理运动边界问题时具有同样的高效率,且精度高于浸入式边界方法。     

基于激波拟合与柔性壁的反设计数值方法研究

将激波和待求激波压缩面作为运动间断边界拟合,研究了超声速来流下的型面反设计数值计算方法.算法秉承拟合法思想,采用二阶格式时间推进欧拉方程组,在运动边界和动网格条件下,按预估的波系结构划分网格块,反求可产生指定压力分布的三维型面,并获得精确的激波型面.通过包括凸包进气道在内的三个反设计算例验证了本文算法的合理性及其对网格质量的适应性,证明了在型面反设计中,采用激波拟合法可大大提高型面精度.本算法可用于超声速进气道压缩面和唇口完全三维反设计.     

非紧致气动噪声半空间传播基本解的边界积分方法

为了获得低马赫数流动诱发的非紧致气动噪声在半空间内传播的基本解,结合复等效源方法和边界元方法建立了半空间精确格林函数的边界积分方程,当半空间边界为刚度型阻抗边界时可避免奇异性积分。同时基于等效源方法提出一种半空间二维非紧致圆柱声散射模型,推导了静止介质中声散射基本解的理论表达式。对静止介质中的二维圆柱声散射,数值解在研究的频率与观察点处与理论解一致。采用数值方法计算二维圆柱绕流诱发的半空间声场基本解,结果显示半空间边界强化了声源的声辐射能力,小于马赫数0.2的介质运动对声传播的影响可以忽略。     

球柱三维分离及涡流的数值模拟

对三维可压非定常Ｎ－Ｓ方程进行有限体积中心离散,并用Ｒｕｎｇｅ－Ｋｕｔａ显式格式进行数值求解。数值模拟了球柱分离及涡流的形态结构,给出了具有激波、边界层、大分离区、涡流及其间相互干扰的复杂流型     

数值求解可压缩N—S方程的增量迭代法

本文在Beam-Warming格式的基础上,引入通量分裂的概念,应用行高斯迭代技术提出了一个二阶隐式的求解BFC坐标变换下的N-S方程的增量迭代法。它消除了因式分解带来的寄生误差,合理地处理了控制方程中的粘性项,具有较好的收敛性和稳定性。二维Poissuill流动和Couette流动的计算结果和精确解吻合良好。S形进气道的数值解反映出的现象和实验结果是一致的。     

An SOR Implicit Time-accurate Scheme in Calculating Unsteady Flows

A new time-accurate marching scheme for unsteady flow calculations is proposed in the present work. This method is the combination of classical Successive Over-Relaxation (SOR) iteration method and Jacobian matrix diagonally dominant splitting method of LUSGS. One advantage of this algorithm is the second-order accuracy because of no factorization error. Another advantage is the low computational cost because the Jacobian matrices and fluxes are only calculated once in each physical time step. And, the SOR algorithm has better convergence property than Gauss-Seidel. To investigate its accuracy and convergency, several unsteady flow computa- tional tests are carried out by using the proposed SOR algorithm. Roe’s FDS scheme is used to discritize the inviscid flux terms. Un- steady computational results of SOR are compared with the experiment results and those of Gauss-Seidel. Results reveal that the numerical results agree well with the experimental data and the second-order accuracy can be obtained as the Gauss-Seidel for unsteady flow computations. The impact of SOR factor is investigated for unsteady computations by using different SOR factors in this algorithm to simulate each computational test. Different numbers of inner iterations are needed to converge to the same criterion for different SOR factors and optimal choice of SOR factor can improve the computational efficiency greatly.     

二维热化学非平衡流动的非结构网格DSMC方法及其应用

研究了二维热化学非平衡流动非结构网格DSMC方法实现的过程。提出了一种新型的高效搜索算法,该算法不仅可以跟踪模拟分子在网格之间的迁移,而且在搜索过程中可以准确判别分子与物面是否相互作用,避免了原有算法中分子表面反射非确定论判据。为加快流场的时间发展历程,设计了适合DSMC方法的动态局部时间步长技术,使其可以应用到定常和非定常流场的计算。利用Fortran90的动态分配内存技术编制了计算程序。最后对过渡流域高超声速圆柱绕流进行了数值试验,计算结果初步验证了该算法的可行性。     

一种适合迭代求解的反馈力浸入边界法

对Goldstein提出的反馈力浸入边界法进行了新的思考，改进了其对力源项的计算，拓展了该浸入边界法的使用范围。传统的反馈力浸入边界法在进行力源项的计算时，含有对速度误差的时间积分项，只能用于含时间项的Navier-Stokes （N-S）方程的求解，且在显式时间推进时有严格的时间步长限制。本文改进的方法则直接通过迭代过程中的速度误差求和来计算力源项，避免了时间相关的参数，使其不仅能适合非定常隐式时间推进，还能与不含时间项的定常N-S方程求解方法结合。为了验证改进方法的可靠性，对二维静止圆柱绕流、静止流体中的振荡圆柱、运动椭圆翼以及三维静止圆球的流场进行了计算，计算结果均与文献结果符合较好，表明本文改进的方法是有效的。得出的结论为：可以直接基于迭代次数进行反馈力源项的计算，改进的反馈力浸入边界法不仅可与非定常N-S方程结合，进行隐式求解，还可以与定常N-S方程结合用于定常流动的模拟，可将改进的方法运用到更多的流动问题当中。     

Kirchhoff方法在旋翼前飞噪声预测中的应用研究

将CFD(计算流体力学)技术与噪声预测的Kirchhoff方法相结合,发展了前飞状态下直升机旋翼远场气动噪声的定量预测方法,得到了旋翼噪声的时域解。先用具有解析解的“点源”流动代替旋翼流动,通过比较计算值与理论值,验证了该预测方法的可行性和可靠性;尔后针对AH—1／OLS旋翼模型的两种前飞状态进行计算,通过比较计算值与实验值,进一步验证了本文方法及其所发展的程序的正确性。     

任意运动变形边界流场测量PIV算法的研究：薄膜涡激振动

提出了一种精确测量任意运动变形边界流场的粒子图像速度场PIV算法,能智能识别强流固耦合问题流场测量中任意运动变形的边界信息,然后生成贴体自适应的图像互相关计算窗口,计算获取运动变形边界附近的速度场数据。这一算法可大大提高任意运动边界附近流场的测量精度,为流固耦合问题的理论分析和数值计算验证提供可靠的高质量实验数据。针对数字合成给定流场的粒子图像序列,采用所发展的PIV算法对粒子图像区域中的运动变形边界进行了精确识别,高质量地复现了原始流场信息。最后,对低速闭式循环水洞中的钝体尾流柔性薄膜涡激振动现象进行了PIV实验测量,获得了柔性薄膜大变形运动状态下的瞬态流场特征。