直通道和弯曲通道中超声速气膜冷却研究

对直通道和弯曲通道两种不同几何形状中的二维平行缝槽形式的超声速气膜冷却进行了数值模拟,分析了在有、无斜激波入射冷却层时的冷却效果.计算结果表明,斜激波的入射使壁面的冷却效率较无斜激波入射时要低.在该计算模型中,无斜激波入射时,弯曲通道壁面的冷却效果好于直通道壁面,而有斜激波入射时,弯曲通道壁面的冷却效果不如直通道壁面.      

圆形出口内转式进气道流动特征

采用数值仿真的方法研究了内转式进气道的流动特征。研究表明:设计状态在近壁面唇罩激波诱发了二次流,进而发展形成流向涡,造成低能流堆积,隔离段出口流场分布不均,消弱了进气道的抗反压能力。有攻角条件下,口面激波偏离唇罩前缘,激波形态发生改变,激波波面中部展向具有准二维特性,压缩面两侧气流压缩变弱,激波层变薄,出现局部膨胀区;有攻角条件下的无黏流场,在进气道压缩段形成三维流向涡,该流向涡促进高能高速气流向壁面迁移,改善了黏性条件下隔离段出口流场的均匀度。     

跨音速喷管流动的数值模拟

数值方法采用的差分格式是具有高分辨率和快速稳定收敛性质的全变差衰减(TVD)格式。喷管壁面形状和反压可任意给定,对于典型轴对称与二维喷管(即壁面由收缩角为45°,扩张角为15°的两条直线用圆弧光滑连接而成。圆弧的无量纲半径等于0.625,其参考长度为喉道半宽),进行了有激波与无激波流动的数值模拟。计算结果与实验结果以及其他计算结果符合良好,此法可推广到非定常跨音速喷管流动的计算,并可用于工程中喷管设计。     

凹面腔内二维激波会聚特性研究

为了研究凹面腔内二维激波会聚的特性,设计、建立了二维激波会聚冷态实验系统,进行了二维激波会聚冷态实验,结果表明:当射流入口宽度增大时,相同射流入射压力下,喷入凹面腔内的射流强度增大,导致凹面腔内的激波会聚强度增强,从而压力脉动幅值增大,频率也随之增大.数值模拟了二维冷态激波会聚过程,从模拟结果中的压力和频率上看,模拟结果与实验结果基本上是吻合的.     

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

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

二维高超声速进气道加速启动过程数值研究

为深入了解大规模分离区在进气道启动过程中的变化规律、影响因素及自持机理,针对简化的二维高超声速进气道加速启动过程进行了数值研究,对比了不同唇口角构型的启动性能、启动过程和分离区变化规律,分析了其中的流动机理。结果表明:(1)上壁面单侧压缩的二维进气道启动性能受唇口角影响显著,随着唇口角从0°增加到12°,启动马赫数呈现先减小后增加的趋势,4°时启动马赫数最小。(2)不同唇口角构型不启动状态都存在大规模分离区,分离区的前半部分接近一致,后半部分差异明显。(3)在加速启动过程中,分离区主要依靠分离激波在上壁面的反射激波维持自身的存在,不同的唇口角构型在相同来流马赫数下分离激波在上壁面的反射激波强度不同,这影响了分离区的自持能力,从而影响了启动性能。     

运用GAO-YONG湍流模型对扩压器内跨声速流动的数值模拟

运用GAO-YONG可压缩湍流方程组,采用同位网格SIMPLE算法,对扩压器跨声速流动中的二维激波/湍流边界层干扰现象进行了数值模拟。将计算得到的流场的时均参数与实验值进行比较,数值模拟结果在激波强度、壁面压力分布以及分离点和再附点位置等方面,与实验值吻合较好,表明GAO-YONG可压缩湍流方程组能够比较准确的模拟较强激波/湍流边界层干扰流动,从而进一步为GAO-YONG湍流模型的正确性及其在可压缩流场模拟方面的适用性提供了佐证。     

TSTO并联分离激波/边界层干扰流动特性分析

针对两级入轨飞行器的缩比模型，通过试验与数值模拟相结合的方式，在马赫数6条件下开展典型级间距状态的激波/边界层干扰流场研究，详细分析干扰区壁面及空间的流动结构与特性。结果表明：试验中模型壁面边界层在激波入射之前为层流状态，在强激波干扰后迅速转捩为湍流状态，因此试验结果在第一道激波作用结束之前与层流计算结果吻合，而在第一道激波作用结束之后与湍流计算结果一致；激波/边界层干扰呈现复杂的三维流动特征和明显的开放结构，强激波在壁面形成的高压区呈弧状向下游展开，轨道级头部产生的入射激波在级间来回反射，强度依次递减；同时，干扰区内存在展向弯曲的主分离线与再附线、沿流向排列的二次分离线与再附线、流动剪切形成的旋涡结构以及包括鞍点、结点、焦点在内的临界点；层流边界层受到激波作用形成的分离区明显大于湍流边界层，同时开放特征更为显著。     

抽吸腔中心隔板对弯曲壁面边界层抽吸性能影响研究

为了研究抽吸腔内隔板对边界层抽吸性能的影响,在弯曲壁面上均匀分布22个抽吸槽,并在抽吸腔内布置中心隔板,改变抽吸腔出口大小以及入射激波的位置。数值计算采用基于有限体积法的二阶迎风格式来离散二维可压N-S方程,湍流模型采用标准k-ε模型,从流场结构、抽吸流量、弯曲壁面的表面摩擦阻力系数、平均总压恢复系数以及平均马赫数等方面对流场进行了分析。结果表明在弯曲壁面抽吸腔内布置隔板会对抽吸效果产生影响,并且影响程度随抽吸腔出口大小以及激波位置的改变而改变。抽吸腔出口非节流时:布置隔板后,激波位于弯曲壁面中段、后段,抽吸效果均无明显变化,进气道性能亦无明显改变;激波位于弯曲壁面前段时,平均总压恢复系数增加8.33%,质量流量增加3.27%,抽吸效果有所改善,进气道性能也有所增强。抽吸腔出口节流时:布置隔板后,与非节流时相反,当激波位于弯曲壁面前段时,抽吸效果无明显改变,进气道性能不变;当激波位于弯曲壁面中段时,弯曲壁面分离泡的长度减小近2倍,质量流量增加66.15%,抽吸效果显著增强,进气道性能变好;当激波位于弯曲壁面后段时,分离泡长度增加了5倍,分离泡高度增加了近2倍,平均总压恢复系数降低8.33%,同时质量流量也减少1.83%,抽吸效果变差,进气道性能恶化。     

斜激波与内凹半圆柱面无粘相互作用

针对三维内转式进气道流动中激波入射内凹壁面主导的复杂波系干扰问题，采用斜激波入射内凹半圆柱面的简化构型，在来流马赫数为6的条件下，通过无粘数值模拟结合理论分析，研究了激波角βi=14°～29°的斜激波在内凹半圆柱面反射形成的三维流场。结果表明，流场对称面均会出现显著高于二维情况的逆压梯度。当βi≤25°时，从侧壁到对称面，斜激波经历了从马赫反射（MR）到规则反射（RR）的转变，形成了MR-RR型流场，转变点处产生的桥激波向对称面延伸，桥激波在对称面反射后产生的压力峰值高于二维斜激波反射；当βi≥25°时，斜激波在侧壁和对称面均发生马赫反射，形成了MR-MR型流场，两种马赫反射分界点处产生的桥激波向侧壁发展，侧壁气流在对称面相撞后产生的压力峰值高于正激波后的压力；当βi=25°时，流场存在MR-RR型和MR-MR型双解现象。通过降维分析理论，揭示了两类流场中转变点和分界点的形成机制，并厘清了桥激波的产生原因和初期演化特征。当βi≥18°时，无粘激波干扰所主导的侧壁气流加剧向对称面汇聚，并在对称面附近产生流向涡对。无粘分析获得的认识，有助于揭示内转式进气道中流动汇聚和流向涡对等现象的形成机...     

谱方法的基本问题与有限谱法

本文概述了谱方法的特点，讨论了谱方法的四个基本问题及解决方法，由此得出谱方法局域化的重要性，重点介绍了作者提出的有限谱法两个发展阶段，特别是最近的基于截断傅立叶积分的有限谱法，列举了有限谱法的一系列计算格式，给出了计算高超声速流，各向同性湍流，地震波的传播，分区算法，复杂形状物体，激波管问题，KdV方程，湍流和激波相干扰等问题的计算结果。     

柔性决策

提出了柔性决策的思想,并对决策科学的发展进行了回顾。为解决柔性决策问题,提出了极大－极小方法、人机交互方法和集包含方法。     

连续体结构拓扑优化方法评述

连续体结构拓扑优化在优化中能产生新的构型,对实现自动化智能结构设计具有重要意义。目前,连续体结构拓扑优化方法主要有：均匀化方法、变厚度法、变密度法、渐进结构优化方法、水平集法、独立连续映射方法。本文首先系统回顾了以上方法的发展历程,介绍了它们的研究现状。其次,通过对比以上拓扑优化方法对若干典型算例的优化结果,表明以上方法都有较好的减重效果。最后,对以上方法进行了总结,列出了它们的优缺点和发展方向。     

基于水平集方法的连续体结构拓扑优化方法研究

考虑到现有各种连续体结构拓扑优化方法的优缺点,提出了一种新的连续体结构拓扑优化新策略。该方法将无网格法和水平集方法结合起来,充分发挥各自的优点,使得该拓扑优化新策略具有良好的收敛性和计算精度。     

一种改进型洞映射法

针对传统洞映射法存在挖洞曲面必须封闭的局限性和消耗大量内存的缺点，提出了改进型洞映射法。主要做了两点改进：1）采用“反向标记法”判别洞映射笛卡尔网格单元属性，使得挖洞曲面不再局限于封闭的情况；2）运用“矢量射线法”进一步判断洞映射笛卡尔网格边缘单元内的计算网格点属性，使得笛卡尔网格分辨率低时仍能精确挖洞。研究表明：相比于传统洞映射法，改进型洞映射法保持了很高的计算效率和精准挖洞的能力，还极大地降低了内存消耗，并且具有较强的可靠性。     

含裂纹铆接加筋板应力强度因子的计算

本文提出了选择复连通域边界(多环路)计算J积分的一种新方法,成功地解决了单连通域内有铆接点时计算J积分的问题。文中还提出了一种“虚设铆钉法”,采用这种方法,可利用同一网格有限元模型计算不同裂纹长度下裂纹尖端的应力强度因子,这大大简化了模型准备工作并减少了计算结构总刚度矩阵的工作量。本文采用刚度导数法、J积分法及能量释放率法,计算了含裂纹铆接加筋板的应力强度因子,通过计算结果对比,对各种方法作了粗浅的评价。     

基于混合 MOM-PO的复合电磁散射问题分析

海面(或地面)环境中目标电磁散射特性,对实际战场环境下军用目标的雷达探测、雷达特性形成以及雷达目标识别等技术具有基础支撑作用。因此,快速准确地分析目标和粗糙面复合电磁散射特性具有重要理论意义和应用价值。文章采用基于矩量法(MOM)和物理光学法(PO)的混合方法MOM-PO分析了目标和粗糙面复合电磁散射问题。与传统的矩量法相比,该方法在保证精确度的前提下大大加快了计算速度。     

典型翼型参数化方法的翼型外形控制能力评估

现代飞行器的精细化设计需求要求翼型参数化方法具有精准的翼型外形控制能力。依据标准的风洞模型公差要求,以五个典型翼型作为测试翼型,分别对Hicks-Henne参数化方法、样条参数化方法、PARSEC参数化方法和CST翼型参数化方法的翼型外形控制能力进行评估,依据每个翼型参数化方法的评估结果并结合其固有特点进行适应性分析。结果表明:四种翼型参数化方法对翼型外形的控制能力存在较大差别,具有不同的适用范围。评估典型翼型参数化方法的所得结果可指导设计人员依据具体的设计需求选择适当的翼型参数化方法。     

关于电磁场的自然单元法

有限元法和无网格法在电磁场数值计算中已经得到了广泛的应用,然而有限元法存在前处理网格剖分问题,无网格法存在计算时间长、边界条件和不连续面难处理等问题。针对以上问题,本文提出了利用自然单元法求解电磁场的方法,该方法根据场域中的离散点的信息构造自然邻点插值函数,最后算出场值。通过算例证明了该方法在电磁场计算中的可行性。