基于Euler／N—S方程的跨音速非线性静气动弹性问题研究

在C-H网格的基础上,采用Jameson的中心差分有限体积法求解Euler/N-S方程,采用结构影响系数法计算结构的弹性变形,用三角元面积加权法和常体积转换法(CVT)实现流固耦合,计算了弹性后掠机翼在跨音速状态的静气动弹性问题     

弹性三角翼的静气动弹性方法研究

静气动弹性问题的研究关键是要解决流固数据耦合的问题，即把气动力从气动网格点插值到结构点，由柔度影响系教法求得结构变形后，再把结构变形插值到气动网格上。本文在O-H网格的基础上，采用Jameson的中心差分有限体积法求解Euler／N—S方程，采用结构影响系数法计算结构的弹性变形，用三角元面积加权法和常体积转换法（CVT）实现流固力的耦合，采用无限平板样条法（IPS）实现位移的耦合。     

跨声速机翼操纵面非定常气动力Euler方程计算

采用代数法与椭圆型方程优化相结合的结构生成方法为带操纵面F5机翼生成面对接块结构网格;弹性变形技术生成动态网格;有限体积法、双时间推进法求解非定常Euler方程.给出一种适于面对接块结构网格计算的守恒量插值方法.数值计算了F5机翼操纵面作简谐运动时的非定常气动力,计算值与试验值吻合较好.     

基于Euler方程的三维自适应笛卡尔网格应用研究

本文采用非结构自适应笛卡尔网格和有限体积法求解了三维Euler方程。在网格生成过程中，采用了基于模型几何数据、模型表面曲率和流场特征为基础的网格自适应生成及加密方法，采用模型特征恢复技术保证了计算模型与原始外形的一致。在计算中，将以中心差分为基础的Jameson有限体积法在三维笛卡尔网格中进行了推广。本文中对导弹、歼击机等复杂问题进行了数值实验，计算结果与风洞实验结果符合良好，并在工程实际中成功应用。     

跨声速风洞机翼半模型试验洞壁边界层影响的数值模拟研究

采用数值模拟方法，模拟垫块法风洞试验二维垫块厚度对半模型试验结果的影响，计算模型采用M6机翼半模型。控制方程采用全场三维可压、雷诺平均的Navier—Stokes方程；采用Jameson中心有限体积、多步Runge—Kutta时间步长格式进行求解；湍流模型采用Baldwin—Lomax模型。通过本文的研究，可以为跨声速风洞半模型试验垫块厚度的选取提供理论计算方法和设计依据。     

飞行器复杂外形流场NS方程数值模拟

本文采用分块对接和“O”型网格技术，首先在飞行器表面长出“一层”分块结构网格来，再将计算区域用若干分块结构网格填充，这样既保证了网格的质量和模拟复杂外形的能力，又大大提高了计算效率。在计算中，采用有限体积法，MuScL型迎风差分格式求解Ns方程，根据多块网格的排序，研究了计算方法。最后，对某飞行器外形进行了亚跨超流场数值计算，计算结果与试验符合较好。     

基于粘性非结构网格的三维增升装置N-S方程数值模拟

应用基于粘性非结构网格求解N—S方程方法数值模拟了某型飞机三维增升装置的绕流。在物面附近粘性主导作用区域内，采用三棱柱型非结构网格，在其他流场区域，采用传统四面体非结构网格的混合非结构网格生成技术。湍流模型采用Spalart—Allmaras一方程模型，运用格心格式有限体积法，对控制方程和湍流模型方程进行数值离散。用四步龙格-库塔方法作显式时间推进，采用当地时间步长、隐式残值光顺等加速收敛措施．求解某型飞机全机带增升装置构型的粘性绕流流场，其数值计算与风洞试验结果显示了良好的吻合度。     

轴对称推力矢量喷管的静态内部性能分析

采用Jamson提出的有限体积法对4种结构的轴对称推力矢量喷管的内流场进行了数值分析，预测了其内流场性能参数，计算结果与参考文献符合的很好，表明该方法对建立轴对称推力矢量喷管的数学模型具有很高的模拟精度和较快的收敛速度。     

二维Euler方程的自适应非结构直角网格算法

本文用四叉树型非结构直角网格和有限体积法求解了二维Ｅｕｌｅｒ方程,在基本网格生成过程中,采用了以几何外形为基础的网格加细方法,包括物面边界网格加细和曲率加细;并提出了一个用自适应的方法在激波处及其附近局部加细基本网格过程,提高了激波分辨率,在计算中将结构网格经常采用的以中心差分为基础的Ｊａｍｅｓｏｎ的有限体积法做了推广,减少了计算工作量,节省了机时,为今后推广到三维创造条件,本文对单个翼型,多段翼     

栅格冀导弹流场混合网格N-S方程数值计算

采用结构和非结构混合网格技术，对栅格翼导弹黏性流场进行数值模拟，预测其气动特性。计算的马赫数为0．7～2．5。用有限体积法和LU-SGS算法求解N-S方程。计算与实验数据以及栅格翼附近流场结构进行了比较。     

一种MWENO谱体积（SVMWEN05）格式

构造一种新的谱体积（spectral volume）格式来求解双曲型守恒律，记为SVMWEN05（Spectral Volume method by Multi—Weighted Essentially Non-Oscillatory）格式。其主要思想是：第一步，将空间计算区域划分为一系列单元，称为谱体积，等分每个谱体积为一些子单元，称为控制体积（controlvolume）。第二步，在谱体积内部采用类似MWEN05（Multi-Weighted Essentially Non—Oscillatory）的格式进行重构，而在谱体积的边界处采用传统WEN05格式进行处理。第三步，利用Runge—KutmTVD离散方法对半离散格式进行时间离散，得到时空一致高精度全离散格式。最后．在文中给出几个经典数值算例用以验证本格式的计算能力。     

贴体与边界正交网格及大迎角涡流计算

本文应用空间理想非均匀线性弹性力学平衡关系生成三维贴体与边界正交网格。通过弹线性的刚度来控制网格的密度，并且自动调整边界节点上外力分布使网格线与边界正交。文中还用欧拉方程有限体积法三步显式时间推进及焓修正方法求解了绕三角翼的大迎角流动。计算结果表明这种方法能够较好模拟跨声速胶体涡流动。     

二元喷管无粘流的有限体积法分析

本介绍了一种计算收敛－扩散喷管内流特性的有限体积法，它将求解域分成许多有限体积单元，将欧拉方程在小体积单元上积分，从而得到半离散方程，然后采用Mac Cormack格式进行时间推进，求得定常解。计算结果与试验值吻合良好，是一种有效的计算具有喉道小曲率半径不连续的收敛－扩散喷管特性的方法。     

微射流作动器三维外流场数值模拟

采用有限体积法求解Faver平均N—s方程，对三维、粘性、非定常、不可压微作动器外流场进行数值模拟。计算结果揭示了微射流产生、发展与耗散的过程，和实验中观察到的现象比较吻合。旋涡对沿展向的不稳定性是涡环破裂、旋涡对的相干性消失的主要原因。二次涡的发现为以后研究微射流旋涡对的细微结构打下良好的基础。     

基于Navier Stokes数值计算的风洞收缩段设计

运用Navier—Stokes数值模拟对不同收缩曲线和设计方案的风洞收缩段和试验段的流动进行模拟，并对收缩段的分离特性、出口速度均匀性以及试验段的速度均匀性和附面层厚度作出评价。计算格式在空间上采用中心有限体积离散，在时间上采用多步Runge-Kutta时间步长格式进行积分。     

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

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

非结构多重网格流场数值模拟研究

本文在非结构网格上应用多重网格技术来加速Euler方程的收敛过程。此方法在一系列粗网格和细网格之间进行值的传递，并且粗网格和细网格之间结构上无关。在每一重网格上用有限体积法进行计算。采用本文的方法数值求解了NACA0012翼型的绕流。     

Jameson有限体积法对非结构网格推广的改进

本文将Ｄｅｌａｕｎａｙ法用于复杂物形的非线性结构网格生成,并利用黄明格提出的布点原理避免边界穿透。绕流的计算采用Ｅｕｌｅｒ方程和Ｊａｍｅｓｏｎ的有限体积法。本文对三维问题提出新的激波敏感器,对方法做了改进,并用于ＡＧＡＲＤ Ｍ６机冀、导弹模型和飞机模型的绕流计算。     

纹影流动显示技术在超声速射流数值模拟验证中的应用

为了研究总压探针测量方法在欠膨胀超声速射流中的可行性，论文采用非结构网格有限体积法对含探针和不含探针的射流流场进行数值模拟，并通过对两个计算结果的比较，给出这种测量方法的参考性评价。为了证明数值方法及计算结果的可靠性，给出了口径6mm的喷嘴、粗1．3mm总压探针在压比4的流动条件下纹影流场显示和总压探针测量等实验结果。     

结构/直角切割网格下时域有限体积法在计算电磁中的应用研究

 采用有限体积法直接求解时域麦克斯韦方程组，研究了二维圆柱、翼型和三维球体、立方体、双球体等典型完全导电目标的电磁散射问题。空间离散使用近似黎曼解构建网格通量，并针对二维结构网格、三维自适应直角切割网格分别给出重构方法，时间方向采用二步龙格库塔法。计算与理论、实验或文献结果符合良好。