跨音速非定常流动的Euler方程隐式解

在随时间变化的贴本坐标系中给出求解非定常Euler方程的连续通量分裂法。在此基础上建立了可用于跨音速非定常流动的Euler方程隐式求解法。采用特征向量变换,可在保证原方程组离散化精度的条件下使计算大为简化。针对振动翼绕流特点建立了固连于物体的动坐标与固定坐标间的关系。数值计算在动坐标中进行,既简化网格生成又保证在物面上满足边晃条件。对NACA64A-10冀型绕1/4弦点做简谐俯仰振动的非定常气动力进行了计算,给出了与实验结果基本相符的计算结果。此外,还给出翼型做沉浮及同时进行沉浮与俯仰二自由度振动的非定常气动力的计算结果。     

翼-身组合体绕流的Euler方程数值模拟

用代数方法生成翼-身组合体H-O型网格,并用有限体积法研制出翼-身组合体绕流三维Euler方程计算程序。该方法的特点是改进了机身与机翼表面网格点分布,机翼后缘有后掠时也能保证后缘与网格线一致。程序除能按常规提供横流截面展向压强分布外,还能提供弦向压强分布。对NASA TND-712的翼-身组合体模型的计算结果与实验符合很好。     

跨音速大迎角Euler方程数值分析

应用有限体积法离散三维Euler方程,选用三步显式格式进行时间推进求解;并通过当地时间步长、残值光顺、焓修正等方法加速收敛。对一维波动方程的von Neumann稳定性分析说明这种三步格式的最大Courant数是2。对三角翼跨音速大迎角流动的数值分析表明文中的三步格式及焓修正方法优于人们经常采用的四步、五步格式以及焓阻尼技术。     

飞机电网零序阻抗的数值计算

从麦克斯韦方程出发,推导了飞机电网零序阻抗参数的计算公式。对几种不同敷设方式载流导线的零序阻抗进行了计算,其结果与标准值较为接近。这种数值计算方法比传统的解析近似估算方法更为合理和准确。     

再入弹抛壳气动特性研究

本文用有限差分方法求解再入母弹头与其所抛壳体之间的无粘干扰流场,获得了体—壳的干扰气动力,并在准定常假设的基础上近似模拟了弹体与壳体的分离过程。数值研究的结果表明,体—壳间存在的干扰会导致壳的压心位置急剧变化,危及母弹安全;抛壳速度是体—壳安全分离的关键参数。     

高精度差分求解气动方程的几个问题

本文探讨了发展高精度格式的必要性，研究了高精度格式与网格以及熵增条件的关系，并发展一种半离散化的空间为三阶精度的格式。从此半离散化的格式出发，可建立多步显式格式及隐式格式。模型问题的计算表明，该三阶精度的格式具有好的精度，且激波附近基本上没有波动。     

用欧拉方程计算机翼与翼身组合体绕流

本文用保角转绘加剪切变换生成三维机翼与翼身组合体Ｃ－Ｈ型网格，并在这类网格拓扑上，用有限体积法研制出可供分析机翼与翼身组合体绕流的三维欧拉方程计算程序。本方法的特色是改进了机翼表面网格点的分布，使机翼后缘与网格线一致；采用了当量机身，使机岙表面到对称面光滑过渡；嵌入了边界层粘性修正。改善了计算结果。数值计算表明，对粘性影响较大的超临界绕流，边界层修正效果显著。     

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

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

燃烧室污染物生成大涡模拟模型的研究

采用亚网格动力学模型和联合概率密度函数(PDF)亚网格模型对污染物NO的生成进行大涡模拟.计算中采用的数学模型有:k方程亚网格尺度模型来估算亚网格湍流粘性;亚网格EBU燃烧模型估算化学反应速率;在交错网格体系下, 采用SIMPLE算法和混合差分格式来求解差分方程.数值分析两种不同亚网格污染物生成模型和两种进口气流温度对NO排放的影响, 计算结果与实验数据比较表明:两种亚网格污染物生成模型都能较好的预估燃烧室污染物NO的生成, 但亚网格动力学模型的计算工作量要比联合概率密度函数亚网格模型少得多, 更适合工程应用.