基于欧拉方程的二维振荡机翼非定常气动设计反命题方法（Ⅱ）

建立了二维振荡机翼的非定常反问题的有限差分解法,并通过若干算例对求解方法进行了验证.首先根据基于欧拉方程的非定常气动反问题数学模型确定了具体求解方案,接着重点讨论了壁面非定常反问题边界条件和远场非定常无反射边界条件的处理方法,编写了相应的计算程序,并以绕1/4弦长点做小幅俯仰周期振动的NACA-65-0012翼型为初始机翼,在不同目标压强的条件下进行了气动反问题计算得到新翼型,它们的压强分布与预设目标相吻合,这表明本文的反问题求解方案与计算方法都是可行的.     

基于欧拉方程的二维振荡机翼非定常气动设计反命题方法

建立了基于欧拉方程组的二维振荡机翼非定常气动设计反命题方法的数学模型.通过一系列变换,将物理时空中的求解域转换成映象坐标系中的规范区域;并导出了映象坐标下的欧拉方程,结合非定常反命题的边界条件,便可用有限差分方法求解.本方法引入映象坐标系解决反命题边界形状的不确定性,并能利用欧拉方程现有的各种差分格式于反命题求解.     

超声速机翼－尾翼组合体的非定常气动力边界元法

本文以Ｇｒｅｅｎ函数为基础的边界元，应用时间历程法，计算机翼－尾翼组合体的非定常气动性能。这为飞行器的设计及其非定常气动特性分析提供了一种有效和完善的计算方法。本文介绍的方法，特别适用于研究、分析各类飞行器非定常动态响应过程的数值计算，例如飞机对突风动态响应等。本文提供的三维机翼的超声速启动过程非定常气动力数值结果，与Ｌｏｍａｘ理论解作了比较，可以看到吻合得很好，说明本文所使用的方法是可行的。对于机翼－尾翼组合体的数值结果，则进一步提供了尾翼受到机翼的洗流影响过程。     

适用于纯物质流体工质S2流面半反命题和A型杂交命题的变分原理族

文献[2]建立了轴流式叶轮机S_2流面半反命题、A型杂交命题的变分原理族。文献[1]在此基础上提出混流式叶轮机S_2流面半反命题和A型杂交命题的变分原理族,避免混流式中可能出现解的不确定性及分区处理的麻烦。但这些文章中所提出的变分原理均只适用于完全气体工质,使这些变分原理的应用受到一定的限制。例如对水蒸汽及氟里昂等工质,在某些极限工况下,如果再以完全气体的模型求解,会与实际情况相差较大。本文提出可将以上各变分原理推广至纯物质流体工质,从而使这些变分原理具有较强的适应性。     

用N—S方程计算翼型非定常粘性大攻角绕流

本文给出了用ＮＳ方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二维时均可压缩完全ＮＳ方程；湍流模型采用双层代数涡粘性模型。使用近似因式分解ＡＤＩ差分格式离散求解，网格是用保角变换方法生成的相对翼型固定的Ｃ型网格。本文给出两类典型非定常绕流数值模拟结果：翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较，表明了本方法准确、高效的特点。     

一维变截面非定常气流杂交命题的广义变分原理

应用坐标变换,把可变动的物理平面,在映像平面上转化成的十分规则的矩形区,这给有限元计算带来了极大的方便。     

超声速机翼—尾翼组合体的非定常气动力边界元法

本文以Ｇｒｅｅｎ函数为基础的边界元，应用时间历程法，计算机翼－尾翼组合体的非定常气动性能。这为飞行器的设计及其非定常气动我分析提供了一种有效和完善的计算方法，本文介绍的方法，特别适用于研究、分析各类飞行器非定常动态响应过程的数值计算，例如飞机对突风动态响应等对突风动态响应等。本文提供的三维机翼的超声速启运过程非定常气动力数值结果，与Ｌｏｍａｘ理论解作了比较，可以看到吻合得很好，说明本文所使用的方法     

牛顿流理论及其在谐振机翼上的应用

根据牛顿流理论,文中讨论了高马赫数、任意迎角、一般剖面在均匀来流中作沉浮和俯仰谐振的翼型和机翼。导得翼型和机翼表面的牛顿撞击压强和安定性导数公式。给出楔、平板三角翼、双抛物线剖面三角翼俯仰刚性导数(-C_(mθ)和阻尼导数(-C_(m))的数值例子。数值结果与现有的牛顿-布兹曼流理论、薄激波层理论、欧拉方程摄动理论的结果作了比较,算例表明本文的结果令人满意。     

二维振荡机翼跨声速非定常绕流的变域变分有限元解

本文首先对文献（１－３）提出的二维振荡机翼含激波跨声速非定常绕流的变分原理进行了改造，使之能适用于时间推进法；构造了三维时空的可变节点有限元来捕获自由尾涡面，而跨声流中的激波采用人工密度法捕获。在远场边界上使用简化的无反射条件和新型非定常Ｋｕｔｔａ条件。用本文方法求解了作俯爷振荡的ＮＡＣＡ６４Ａ０１０翼型的跨声速绕流，计算结果令人满意。本文方法可以推广到三维机翼及二维和三维叶栅的同类气动力问题上去     

机翼和叶栅非定常流的Kutta条件

本文在现有的实验研究和数值计算的基础上，对经典Ｋｕｔｔａ条件在非定常流动中的适用性作了进一步探讨，并对文献〔１，２〕中提出的一种新的Ｋｕｔｔａ条件－Ｇｉｅｓｉｎｇ－Ｍａｓｋｅｌｌ（Ｇ－Ｍ）模型进行了分析与比较，本文给出了二维振荡机翼和叶栅绕流的变域变分有限元解法中Ｋｕｔｔａ条件的处理方法。本文旨在为机翼和叶栅非定常绕流的数值计算提供更合理可靠的尾缘边界条件。     

多目标自然层流翼型反设计方法

 进行了基于扰动放大N-因子的目标压力分布设计方法的多目标自然层流（NLF）翼型反设计方法研究。流场分析和转捩位置计算用XFOIL程序,大大减少了计算花费。用N-因子设计方法进行有NLF范围要求和满足气动约束的目标压力计算,压力恢复段的压力分布用Stratford分离准则来进行设计。用基于响应面方法的优化方法来进行反设计计算,使用不含二阶交叉项的二次多项式模型的响应面模型,大大减少了构造模型所需的试验次数;设计空间内试验点的选取满足D-优化准则。根据设计目标的设计状态,进行了多目标翼型反设计。计算结果表明,设计结果的层流范围和设计目标基本吻合,该方法可以用于NLF翼型的多目标反设计中。     

跨声速无粘流绕振动翼型流动的数值模拟

本文从Euler方程出发,利用LU-ADI格式对绕二元振动翼型跨声速流场进行了数值模拟。文章中给出了绕NACA0012和NACA 64A 010翼型不同状态的几组数值结果,在同某些实验和其它文献结果基本吻合的前提下,研究了振动对翼型表面压力、激波位置的非定常效应以及其它振动参数对流场非定常效应的影响。     

考虑地面效应的翼型动态特性数值模拟

为了研究考虑地面干扰的飞行器非定常气动特性,基于滑移网格技术,通过求解二维非定常Euler方程,对不同近地高度下的NACA0012翼型的动态振动特性进行了数值模拟。分析了存在地面干扰时的升力系数以及俯仰力矩系数迟滞环变化,并结合强迫振动法,计算了翼型各个近地高度下的纵向组合动导数。结果表明,地面效应不仅对定常流场产生影响,更显著地影响了非定常气动力及力矩,近地高度越小,这种效应越明显,使得升力系数的迟滞环面积变小,幅值减小,而俯仰力矩系数的迟滞环变化不规律,力矩系数的动导数随着近地高度减小而增大(绝对值减小),纵向的动态稳定性产生损失。地面效应干扰在飞行器非定常气动研究中应该引起重视。     

翼型舵面偏转非定常流动数值模拟

针对NACA0012翼型舵面偏转问题,数值模拟了不同参数对翼型气动特性的影响。基于非结构动网格技术,采用ALE有限体积描述下的二维可压缩非定常N-S方程,计算通量采用Vanleer格式、时空二阶格式,利用Venkatakrishnan限制器抑制数值振荡。非定常计算结果表明,NACA0012翼型绕1/4弦点作周期性俯仰振动的升力系数和俯仰力矩系数结果与实验数据吻合良好,验证了数值方法的准确性;在翼型舵面表面有分离区产生,升力系数和俯仰力矩系数形成滞回环,在亚声速情况下,滞回环幅值较小,进入超声速阶段以后,幅值增大,随着翼型间缝隙宽度逐渐增加,翼型升力系数和俯仰力矩系数与无缝翼型相比逐渐降低。     

基于时间谱方法的振荡翼型和机翼非定常黏性绕流数值模拟

 传统的双时间方法在非定常计算中长时间的过渡迭代推进求解导致其计算效率相对较低,针对周期性非定常流动问题的流动特征,发展了一种基于离散傅里叶变换的高效时间谱方法,用于求解振荡翼型和机翼的非定常黏性绕流。在时空耦合的雷诺平均Navier-Stokes(RANS)方程的求解中,对流项的离散应用了Roe的通量差分格式,物理时间项的离散方法为时间谱方法,伪时间推进采用了隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)格式。考虑到湍流的时空耦合效应,时空耦合的Spalart-Allmaras一方程湍流模型的物理时间项同样采用时间谱方法进行离散。为了进一步提高计算效率,当地时间步长和多重网格技术等加速收敛的措施均被采用。算例对俯仰振荡NACA0012翼型和Lann机翼的周期性非定常流场进行了数值计算。结果表明:对于周期性非定常流场的数值模拟,相比于传统的双时间方法,用时间谱方法近似物理时间项,不仅能够提高流场的计算精度,而且更能够大幅度提高计算效率。     

返回舱自由振动跨声速非定常流场数值模拟

本文利用NND差分格式对返回舱俯仰运动和跨声速绕流流场进行了耦合求解,结果表明来流M