二元喷管跨音速流场的数值计算

本文采用有限体积法进行空间离散,Runge-Kutta法进行时间推进的办法,通过求解Euler方程来模拟二元收敛—扩散喷管跨音速流场。运用隐式残差平均技术和当地时间步长法加速计算收敛。为了消除气流参数的波动和激波前后的振荡,在方程中添加了自适应耗散项。在100个时间步以内便可获得稳态解。计算结果与试验数据和其他数值方法吻合良好。      

三维非结构自适应多重网格技术

发展了一种基于几何外形的非结构自适应多重网格技术。根据流场参数的变化梯度确定加密边,由加密准则对粗网格进行自适应剖分得到分布合理的较密网格。通过预先生成的初始极密表面网格来将边界网格的加密点投影到边界上,使得边界保持初始外形。为了提高流场求解速度,将自适应加密后的一系列网格作为多重网格,采用了多重网格算法进行计算。M6机翼的流场计算表明该方法可以准确捕捉激波位置,并且表面的压力系数分布与实验值相吻合。     

A mesh optimization method using machine learning technique and variational mesh adaptation

Computational mesh is an important ingredient that affects the accuracy and efficiency of CFD numerical simulation. In light of the introduced large amount of computational costs for many adaptive mesh methods, moving mesh methods keep the number of nodes and topology of a mesh unchanged and do not increase CFD computational expense. As the state-of-the-art moving mesh method, the variational mesh adaptation approach has been introduced to CFD calculation. However, quickly estimating the flow fi...     

基于控制理论的Euler方程翼型减阻优化设计

在Jameson的控制论气动优化设计思想下，应用Euler方程研究了翼型的反设计和减阻问题。设计过程中的梯度是通过求解一个伴随偏微分方程而得到的，每个设计迭代只需求解一次Euler方程和一次伴随方,一与设计变量数无关。在具体数值执行中：①利用分部积分将目标函数变分的计算转化成了类似于计算Euler方程残值的形式，节省了机器时间；②根据控制理论的要求，将减阻问题转化成了相应地修改伴随方程的物面边界条件和目标函数的变分表达，使问题得到了简化；③利用特征不变量分析法，处理了伴承方程的物面和远场边界条件。设计算例证明了本文方法可靠性好、收敛快、大大节约设计时间。     

编队卫星队形重构防碰撞最优轨迹规划

 针对编队卫星队形重构的轨迹规划问题,提出了直接配置混合整数线性规划(DCMILP)方法。首先将卫星编队飞行问题进行简化,整个过程使用线性化描述;继而将三阶Simpson方法扩展至编队卫星的队形重构过程中,将各卫星的状态量和控制量在各节点处离散化;然后根据目标函数及碰撞规避问题等各种约束条件,将整个过程转化为混合整数线性规划问题,从而可以找到该非凸问题的全局最优解。最后,通过对三维和二维两组编队卫星队形重构进行仿真,由结果可以看出,与传统方法相比该方法快速有效,能够满足实时性的要求,使得卫星编队的自主运行成为可能。     

基于自适应混合网格的脱体涡模拟

基于混合网格和CGNS（CFD General Notation System）数据结构,建立了一种各向同性加密/稀疏的网格自适应方法。在悬空点的后处理中,让含有悬空点的单元转化为任意多面体,从而简化了自适应单元剖分模版,同时自适应网格单元之间可完全相容,自适应生成的网格能够直接用于可处理任意多面体的流场求解器。将该自适应方法与脱体涡模拟（DES）算法相结合,开展了65°后掠三角翼大迎角流动的数值模拟应用,并与初始网格的模拟结果进行了详细比较。对比表明：采用网格自适应方法适当增加局部网格量,能够以较小的成本迅速提高三角翼背风区的空间分辨率,增强数值模拟对小尺度涡系结构的解析能力,从而弥补了基于混合网格的脱体涡模拟中常用二阶格式计算的空间分辨率相对偏低、不利于湍流多尺度结构精细模拟的不足。     

基于混合优化的运载器大气层内上升段轨迹快速规划方法

针对运载器大气层内的最优轨迹快速规划问题,提出一种将求解最优控制问题的间接法与直接法相结合的混合优化方法。首先,基于最优控制问题的一阶必要条件,将运载器大气层内的三维最优上升问题转化为Hamiltonian两点边值问题;然后,采用直接法中能以较少的节点获得较高求解精度的Gauss伪谱法进行求解,提高算法的求解效率;最后,采用真空解析解初值及密度同伦技术,解决初值猜测与算法收敛困难的问题。仿真结果表明,混合优化算法能够准确、快速地对运载器大气层内的最优上升轨迹问题进行求解,并在计算精度与效率上均优于间接法,可应用于运载器的轨迹在线规划与闭环制导。     

运动激波自适应网格算法中权函数问题的研究

变分原理的自适应网格技术被应用到运动激波问题的求解上,在解的大梯度区自动加密网格,从而非常成功地算出了激波。通过分析发现权函数选择和取值对自适应网格技术至关重要     

基于非结构动网格的非定常流数值模拟方法

给出了基于非结构动网格的非定常流动问题的数值模拟方法,重点说明两种动态非结构网格方法的原理以及在非定常流模拟中的算法.一种是基于非结构网格的重叠网格方法,该方法结合了非结构网格方法和重叠网格方法的优点,以物面距为准则参数对网格节点进行分类,进而建立重叠网格的网格问插值关系,通过采用动态非结构重叠网格算法,可以处理多体间...     

基于POD-RBFN降阶模型的串列叶栅流场预测

通过传统实验或CFD 手段获取流场信息的方法往往需要耗费大量资源或时间,这在需要快速获取大量流场信息时产生的成本是无法接受的,发展比传统CFD 更快速的流场预测方法具有重要意义。采用本征正交分解（POD）方法对样本流场进行模态分解,提取流场的主导模态;而后采用径向基函数神经网络（RBFN）响应POD 基函数的系数,实现流场降阶预测模型的构建,并在模型中采用基于函数响应偏差的自适应抽样方法;通过某串列叶栅非定常流场数据对预测模型进行验证。结果表明：本文构建的POD-RBFN 混合模型可以快速准确地预测出串列叶栅的流场参数分布;与静态采样相比,本文采用的自适应采样方法在采样效率上表现出明显优势,同样重构精度所需的样本数降低了25% 左右。     

Aerodynamic benefits of adaptive wing technology

Especially in the case of large transonic transport aircraft, flight conditions change considerably during a typical mission. This is accounted for by multiple but fixed design points which compromise aircraft performance. Employing adaptive wing technology where the wing geometry, or other means of flow control, adjusts the flow development to the changing freestream and load conditions allows us to explore fully the flow potential at each point of the flight envelope. Various means of flow control by geometric adaptation and by direct boundary layer control have been investigated within the German national program ADIF and the EU-project EUROSHOCK II and their potential explored. Here, corresponding results are presented and discussed, indicating the applicability and benefits of the adaptive control methods considered. It is also demonstrated that, generally, flow control must be adaptive to work in real aeronautical conditions since these conditions change within the mission flight envelope.     

自适应机翼技术的分类和实现途径

自适应机翼技术研究可分为通过机翼结构较小尺度变形的流动控制设计和较大尺度改变机翼几何构型的自适应结构设计两个范畴。改变机翼构型的自适应结构又包括可变前后缘结构、扭转机翼盒段结构、可变展弦比机翼结构这三种实现方式。根据目前自适应机翼技术的研究现状,归纳出了实现机翼自适应功能的两种途径,其中,采用智能材料结构进行驱动控制的研究代表了自适应机翼技术的发展趋势,而基于传统材料结构的自适应机翼技术则在现阶段更具有工程应用价值。     

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

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

气体运动论BGK格式的翼型绕流数值模拟

以气体运动论BGK(Bhatnagaar-Gross-Krook)格式为基础,通过改进非均匀网格下的通量计算方法,把气体运动论BGK方法应用到二维贴体网格上模拟了可压缩的无黏和黏性流动。在重构阶段,以交界面周围4个网格单元为一组,先把每个网格单元的宏观量(密度、动量和总能量)转换在相对于交界面的坐标系下,然后使用van Leer限制器插值求出交界面的宏观量,在求解阶段再使用BGK理论求出流过交界面的通量。通过上述方法使贴体网格和均匀网格,在计算格式上一致,与其他求解二维贴体网格的气体运动论BGK方法相比,该方法相对简单,计算量减小。结合贴体网格通量计算和气体运动论BGK格式方法理论,分别对亚、跨、超声速下的NACA0012翼型绕流进行了数值模拟和对比研究,所得可压缩无黏流动的数值结果与北大西洋公约组织航天研究与发展咨询组(AGARD)的结果吻合良好。本文还模拟了低雷诺数状态下的可压缩黏性流动,所得结果与文献中的数值计算结果吻合良好,证明了该方法在数值模拟及激波捕捉方面具有鲁棒性,为在复杂流动中的应用提供了一种简单的途径。     

联合使用激波捕获—拟合法求解激波间断的跨声速流场

吸取激波捕获和激波拟合两种方法各自的优点,发展了一种求解有激波的跨声速流场的方法,并对回转叶栅中定常跨声速流场做了计算。计算中,通过捕获法确定初始激波位置,然后经反复修正,拟合出确定的流场通道激波。计算表明,这种方法可自动获得清晰的激波,流场中气流各参数分布合理。这一方法可适用于复杂边界和不同进口M数,计算时间仅比势函数方法多一倍左右。     

超燃冲压发动机燃烧流场纹影图像无监督分割方法

超燃冲压发动机燃烧室流场纹影图像常存在大量噪声信号,如何高效、准确提取燃烧流场图像的主要波系结构成为当前亟需探索的问题。以超燃冲压发动机燃烧室内冷流到氢燃料点火阶段流场为研究对象,基于深度神经网络方法,发展一种燃烧室内流场的关键波系结构快速识别方法。首先,采用基于图论的超像素分割方法对纹影图像进行聚类分割,为语义信息明显相同区域分配伪标签;其次,设计了一种全卷积特征提取神经网络,并使用残差结构对各个通道进行加权,提取纹影图像高级语义特征;最后,使用交叉熵目标函数优化网络模型,并通过阈值滤波操作去除噪声像素点,提升语义分割效果。结果表明：与K-means及自适应高斯阈值方法相比,本文提出方法在准确率、召回率、F1分数和交并比指标性能明显是最优的,能够准确完成燃烧流场纹影图像语义分割任务,可以更加清晰地反应流场内的主要波系和剪切层结构     

发动机短舱溢流阻力的数值模拟

 为分析发动机溢流阻力的产生原因并为进气道-发动机匹配研究提供一定的技术支持,以典型的发动机短舱为研究对象,分析了短舱上的进气道阻力,并详细介绍了相应溢流阻力的数值计算方法。在对发动机短舱外流模拟时,通过引入两个和发动机参数相关的数值边界条件,即一个是定流量的发动机进气道入流边界,另一个是定总温总压的发动机喷流边界,避免了对发动机复杂内流的模拟。利用NACA-1-81-100发动机进气道作为标准算例,验证了数值方法的可行性。通过模拟发动机短舱在不同工作流量下的流动,认为发动机短舱的溢流阻力产生归因为溢流造成的发动机外表面的阻力变化,以及捕获流管变细导致的附加阻力增加。     

基于QSS的自适应多步校正算法及其在柔性航天器动力学中的应用

量化状态系统（QSS）算法是一种基于状态变量离散化的数值积分方法,该方法与基于时间离散的传统方法显著不同,QSS通过计算状态变量每次跃迁所需的时间来推进下一步的积分。在求解非刚性常微分方程时,QSS算法比传统算法更具优势,但它不适合求解刚性问题。为此提出一种基于QSS的自适应多步校正算法（AMCQSS）,该算法以QSS为基础、结合隐式多步法思想,在计算过程中可以自适应选择二步法或三步法,以有效提高求解刚性问题的精度和效率。通过对柔性航天器动力学的仿真求解,验证了算法的可行性。将该算法与ODE23tb、ODE15s、ODE45以及QSS等算法进行对比,结果表明AMCQSS算法既能保证求解的效率及精度,又具有较好的收敛性和稳定性。     

Fluid dynamics of unsteady separated flow. Part I. Bodies of revolution

Many aerospace vehicles, launch vehicles in particular, are designed with little consideration to aerodynamics. As a consequence, the aerodynamic characteristics are often dominated by the effects of separated flow. The vehicle designer must, therefore, be able to assess what these effects will be especially in regard to the vehicle dynamics, as those effects are often adverse and always of large magnitude relative to the attached flow characteristics. In spite of rapid development of computational means, purely theoretical methods for prediction of the effect of separated flow on rigid and elastic body dynamics are not presently available and will not be for some time. To compound the problem experimental dynamic data by which one can determine the effect of flow separation on full scale vehicle dynamics are usually not available either. The present paper shows a practical solution to this dilemma. An analytic method is described that uses static experimental results to predict the separated flow effect on rigid and elastic vehicle dynamics. A key parameter in the analytic relationship between steady and non-steady aerodynamics is the time lag occurring before a change of flow conditions can affect the separation-induced aerodynamic loads. The time lag causes a statically stabilizing load to produce negative aerodynamic damping, and vice versa. This time-lagged quasi-steady theory provides predictions that are in good agreement whith dynamic experimental results for a great variety of vehicle geometrics and associated types of flow separation.