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

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

双锥Bump压缩面设计及气动特性

采用反设计方法研究双锥Bump压缩面设计技术,该方法实质是融合运动间断边界在运动网格条件下对轴对称欧拉方程组求解.为避免激波捕获法对激波型面位置求解误差,应用计算域分块的方法预估第2道激波,在此基础上采用流线追踪法生成Bump压缩面;结合平面机身设计一双锥Bump实例,运用计算流体力学仿真手段对其进行黏性数值模拟.研究结果表明:①该模型流场结构仍保持较强的附面层扫掠能力;②在来流马赫数为2.0条件下,相比传统的正圆锥乘波体Bump压缩面设计,新型设计方法可使Bump外压缩系统总压恢复系数提高0.04左右,为Bump进气道性能提高奠定基础.      

动态间断装配法模拟斜激波壁面反射

基于非结构动网格技术和边界装配思想提出了动态间断装配法,该方法能够应用于求解含有间断的流动问题。无论入射激波还是反射激波都是作为边界进行处理,激波运动速度由兰金-许贡纽(Rankine-Hugoniot)关系确定。激波作为动网格的一部分,其运动由动网格技术实现。采用该方法模拟了超声速二维流场中激波与壁面相交问题,并且与捕捉法进行比较,二者的流场结构符合良好,但是在细节上还是存在明显差异。通过对流动结构的分析,得出采用装配方法得到的流场要优于捕捉方法的结论。激波壁面反射的问题模拟,也说明了边界激波装配方法对于复杂的激波相交问题是具有处理能力的。     

管道中复杂马赫反射流场的高精度数值模拟

本文分析了一类高精度迎风TVD欧拉方程差分格式的精度,并用240×80个网格点进行马赫反射的计算,获得了目前实验和计算不易得到的两个激波三叉点的复杂波系结构。本文还研究了激波管上壁面对马赫反射的影响,并计算了轴对称扩张管道中激波的运动及激波与半球形障碍物的相互作用。计算结果分析指出,本文所采用的高精度迎风TVD格式与计算方法可以用来模拟具有复杂波系的工程问题研究。     

一种改进的非定常激波装配算法

基于非结构动网格技术和格心型有限体积方法,提出一种改进的非定常激波装配算法,进一步拓展了其在包含有运动激波的非定常流场的应用范围。首先,针对激波在直/曲壁面传播这类问题,分别建立了壁面间断节点的运动模型;其次,为保证激波在大范围运动时装配阵面不产生失真,基于Bézier曲线拟合方法实现了间断节点分布的自动重构;接着,通过嵌入局部网格自动重构模块,提高了算法的计算效率和自动化程度;最后,对于激波相交点的运动,设计了一种根据位移推算速度的方法进行装配。数值算例表明,所提算法能够有效地处理激波传播问题,相比激波捕捉方法可以提取更多的流场信息,同时可以获得流场间断更加直观清晰的图谱。     

高效欧拉方程解法求解轴对称引射喷管三股流问题

推导了三维流动的广义欧拉方程黎曼增量梯度形式的欧拉方程及其轴对称形式表达，针对引射喷管三股流问题的特点在运动网格下作了H变换。在此基础上将高效欧拉方程解法应用于求解轴对称引射喷管三股流问题，为解决次流通道狭窄带来的问题，采用了弱耦合法处理滑流边界。     

超声速平面钝体流动的快速解法

本文介绍了一种平面流动的快速欧拉方程解法。该方法将原参数非定常欧拉方程组重新组合成以广义称曼变量表示的欧拉方程组，再用二点二步迎风格式离散解。针对钝体流动，本文先建立了动网格下的方程，构造了动网格下的算法。提出了一种简单的激波处理方法。计算结果表明，该方法速度快，稳定性好，对初场不敏感。     

燃气自由射流的正格式数值模拟

采用二阶正格式方法对复杂燃气自由射流进行了数值模拟。将二维守恒方程的正格式方法发展到轴对称Euler方程组的求解,并对不同欠膨胀压力比下的燃气射流进行了数值计算。计算结果表明,该方法能够较好捕捉到波胞、滑移线、射流边界以及三波点等复杂射流流场波系结构,与实验照片反映的流动特征及高精度、高分辨率的三阶间断有限元方法的计算结果吻合较好,马赫盘位置与理论预估值和实验测量结果误差较小。证明该方法对激波具有较强的捕捉能力,为此类流动的数值研究提供了一种新手段。     

间断探测器在间断Galerkin方法中的应用

在二维非结构网格上,对高阶间断Galerkin方法求解定常二维欧拉方程进行数值研究。为了很好地抑制解在流场间断附近处产生的伪振荡,采用了间断探测器和限制器相结合的方法,并分别对Krivodonova间断探测器和基于激波物理特征的激波探测器进行了比较和研究。数值结果表明：如果流场中只存在激波时,两者探测效果相似,且后者更适合应用于求解激波问题。如果流场更加复杂,即存在多种间断时,前者依然可以比较准确地用来探测间断单元,能够很好抑制流场间断附近处产生的伪振荡,而后者失效。     

一种用于间断装配的流场结构辨识算法

使用激波装配法时,初始激波是否准确将会对计算过程产生影响。为了确定初始激波的位置,提出了一种新的流场结构辨识算法。该算法以捕捉法计算得到的流场作为系统观测数据,根据密度、压力等参数从该数据中获取激波和接触间断等流动特征周围的网格节点作为离散点集。通过将该离散点集分割成若干子区域,在各子区域内进行分片拟合,最终将离散点集拟合成连续光滑的实体模型,并将此作为初始激波面。在二维方法的基础上,通过引入单位球模型成功将该辨识算法拓展到三维应用。结果表明,采用该方法获得的间断曲面(激波和接触间断)与捕捉法流场中的间断分布吻合较好,作为初始间断面用于装配法可快速得到收敛解。该方法解决了应用激波装配法时确定初始间断面的难题。此外,该方法还可用于网格自适应方法。选择不同流动参数,可以获得相应流场特征结构的空间曲面,在此曲面的基础上可进行网格局部加密或重剖分。该流场结构辨识算法用于网格自适应具有网格尺度自由设置的优势。