高超声速再入飞船气体热辐射计算边界虚网格方法

针对再入飞行器激波层辐射加热,应用有限体积法计算辐射传输方程,引入了一种边界处理方式——虚网格方法,根据化学非平衡流场数据,计算出飞行器表面辐射加热热流。通过分析发现:采用虚网格方法处理辐射传输边界,能够确保真正的边界条件得到应用,满足真实的物理意义,且辐射传输计算的空间网格、立体角网格、初始温度值的收敛性很好,计算结果精度误差小于2.5%。利用虚网格方法计算飞船再入返回时的辐射加热,结果表明计算所得壁面辐射热流与其他研究者辐射热流计算结果符合良好。      

跨声速欧拉方程的一种隐式有限元解法

本文用有限体积Ｇａｌｅｒｋｉｎ法求解了跨声速二维欧拉方程。求通量运用了Ｏｓｈｅｒ格式，通过外插把原来只有一阶精度的逆风格式改进成二阶精度。并采用基于一性化流量的隐式格式和非结构化网格，使程序具有较高的收敛效率。     

基于非结构网格的多段翼型流场数值模拟

将非结构网格生成技术应用于多段翼型流场数值模拟，针对多连通域可以快捷、迅速地生成非结构网格。通过控制背景网格节点的网格生成参数，进行网格的局部加密，提高多段翼型流场分析的精度，应用有限体积法求解欧拉方程，应用各种加速收敛措施提高流场求解速度，并将数值计算结果与风洞试验结果进行了对比分析，证明该方法具有型号设计应用价值。     

基于非结构网格的不可压N-S方程多矩有限体积法

提出了一种基于三角形及四面体非结构网格的有限体积法(FVM),用以鲁棒且精确地求解不可压粘性流动问题.与传统的FVM方法仅将体积分平均值(VIA)作为计算变量的做法不同,本文提出的方法将VIA及点值(PV)同时作为计算变量并在每个迭代步进行计算更新.VIA以通量形式进行计算以确保数值守恒,PV可以通过控制方程的不同形式进行求解更新,无需守恒,因此可以采用非常高效的方法进行求解.将PV作为增加的变量使得紧致网格模板得以实现更高阶精度的重构,而且由此获得的数值模型对于非结构网格变得更鲁棒.本文针对二维/三维的三角形/四面体非结构网格提出了数值格式,给出了几个基准测试算例,验证了本文提出的数值方法在采用非结构网格求解不可压粘性流动问题时的精确性和鲁棒性.     

基于Lattice Bohzmann方法的翼型绕流数值模拟

将通过物理平面的结构化C型贴体网格变换到计算平面的均匀正交网格,引入非均匀网格插值方法,采用基于Lattice-Boltzmann方程的D2Q9模型,对低雷诺数下的NACA0012翼型绕流进行了数值模拟,对五个不同位置的边界层速度型与基于有限体积法的CFL3D程序的结果进行了比较.结果表明,两种方法的计算结果一致.改进后的Lattice-Boltzmann方法适用于曲边边界,非均匀网格,同时对计算速度影响不大.Lattice-Boltzmann方法计算较简单,因而用LB方法来模拟翼型绕流是一种新的选择,并可应用于更复杂的低雷诺数流动中.     

非结构网格上N—S方程解法研究及在多段翼型绕流中的应用

在非结构网格上求解N－S方程时，靠近物面的网格伸展比相当大，用有限体积法进行显式求解时不易收敛并且雷诺数有一定的限制。本文采用的人工耗散项中含有网络的方向尺度，因而在不同区域选取不同：在粘性区内该项比物理粘性项小得多，在粘性区外与Jameson的人工耗散项相同。生成非结构网格采用的是层推进技术与传统的阵面推进法相结合的方法，这样能保证在粘中性区内网格有一定的方向规则，而在粘性区外网格完全非结构化。     

基于浸入边界法的低雷诺数流固耦合数值模拟(英文)

提出一种基于SIMPLE算法的非定常流固耦合计算方法。流体Navier-Stokes方程空间采用非结构化网格有限体积法离散,时间项采用了欧拉隐式方法。利用浸入式边界方法模拟静止或者运动固体区域,流固界面作用力通过流体体积(VOF)方法进行处理。从而可以用固定网格求解任意复杂区域中的流固耦合作用。本文模拟了低雷诺数静止及振荡圆柱绕流,所得结果与文献中贴体网格计算结果吻合,从而验证了本文方法的合理性和正确性。     

航天器在轨全过程表面辐射热计算数值仿真研究

对在轨航天器表面辐射热计算进行了全过程数值仿真研究。航天器结构较复杂,针对不同结构进行区域分解,对几何模型进行相应的规则化,同时采用结构化网格和非结构化网格建立通用的计算网格自生成技术。仿真过程重点考虑了任意曲面的网格自动划分和任意形状交界面的数据传递,兼顾几何结构、物理过程、计算精度和计算速度。将有限元法和能束均匀分布法相结合计算角系数和辐射传递系数。将积分法和能束均匀分布法相结合计算外热流。由于在轨航天器表面多用多层隔热组件包裹,针对这部分结构采用节点网络法和控制容积法计算其表面温度,而未被包裹的结构采用有限元法计算其表面温度。对具有辐射换热关系的非连通区域温度场的有限元计算进行了分析和公式推导。最后,用Microsoft VC++6.0编程设计开发了近地轨道航天器表面辐射热计算仿真软件。     

基于CFD动网格技术的飞艇动导数计算方法

由于飞艇和飞机这两种飞行器在飞行原理、气动外形等方面都存在较大的差别,因此现有的飞机动导数工程估算公式和经验数据并不适用于飞艇动导数求解。为了解决上述问题,提出了一种新的基于CFD动网格技术的飞艇动导数计算方法。通过非结构化动网格技术,建立飞艇在流场中的强迫振动模型,采用有限体积法对N—S方程求解,研究飞艇在强迫振动中的气动特性变化规律,并基于数值计算结果求解飞艇的动导数。研究发现,得到的计算结果与风洞试验结果基本一致,表明基于CFD动网格技术进行飞艇的动导数求解是有效可行的。     

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

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

多级轴流压气机设计及其特性预估

基于大量的二元平面叶栅吹风实验数据关联结果,发展了一种适用于多级轴流压气机设计的计算程序,以具有实验数据的某跨音速两级轴流风扇为研究对象,采用在子午面上叶片排间隙中设置计算站的流线曲率法进行数值模拟,计算结果与实验数据进行对比,结果表明:模拟的结果能够满足工程应用的需要,验证了该程序在设计多级轴流压气机的有效性。     

神经网络在高空长航时无人机翼型多点优化中的应用

NSGAII算法在翼型多点设计中有着广泛的应用价值,然而其巨大的计算资源和计算成本限制了它的使用。为了解决这个问题,本文引入具有较强非线性映射能力的神经网络代理模型,采用实验设计结合BP法训练神经网络响应面来代替N-S方程求解翼型的性能。在实验点的数值模拟中,为了进一步节省计算资源,提高计算效率,采用网格的变形代替网格的重新划分,使得计算网格的更新速度提高了约50%。在翼型的参数化过程中,采用改进的PARSEC方法,用较少的参数实现了翼型的精确控制。为了增强神经网络的泛化能力,采用12-7-4-3-1的隐层结构。对NLF1015翼型的多点优化算例表明,此方法不仅显著降低了整个优化过程的计算量,而且对翼型的气动性能预测也具有较高的可信度,在高空长航时无人机的翼型设计中具有一定的潜力。     

一种基于高阶剪切变形板模型的有限元算法

为了较好地模拟复合材料层合板的变形并得到较为精确的应力分布,本文基于Reddy提出的高阶剪切变形板模型建立了一种新型有限元。该模型保证横向剪应力沿板厚度方向呈抛物分布,在板的上下表面为零,从而很好地满足力边界条件。有限元算法构造的出发点是修正的虚功泛函,引入Lagrange乘子放松了挠度在单元边界之间的连续性,使得挠度的构造变得简单;由于通过常应力分片检查,所以有限元算法的收敛性得到了保证。此外,本算法的数值结果具有对单元网格形状不敏感、精度好和计算效率高的性质。     

近距-1协议在星间测距中的应用

针对近距航天器测距的需要,将CCSDS(空间数据系统咨询委员会)近距-1协议的定时业务与微波双向单程测距法相结合进行测距。该法以直序扩频技术来获得高精度传输帧ASM(附属同步标识)最后1位下降沿的时刻,并以时间码形式,通过帧数据域在测距星之间相互对传,综合考虑协议给出的信号物理特性,对所给出的方法进行了误差分析。分析表明,航天器之间距离不远,相对速度不大时,该法具有较高精度,且与数据通信共用1套设备,降低了成本。     

含分层损伤复合材料层合板压缩剩余强度计算

本文针对含分层损伤复合材料层合板,基于三维复合材料累积损伤分析方法,应用了非线性有限元技术,考虑了子层屈曲对层合板剩余强度的影响,提出了一种含初始分层损伤的复合材料层合板压缩剩余强度计算方法。在ANSYS软件平台上利用APDL语言编写相应的计算程序,对五块含不同初始分层损伤的复合材料层合板进行了压缩剩余强度计算和结果分析,计算结果与试验结果相比较,吻合良好。     

某型飞机发动机的安装推力计算

本文主要研究了一种计算涡扇发动机安装推力的方法。将INSTAL标准进气道/尾喷管设计参数与实际飞机进气道/尾喷管的设计参数相比较,根据INSTAL资料库中介绍的特性数据转换方法,将标准进气道/尾喷管特性参数转换为实际进气道/尾喷管特性参数。再结合飞机发动机的非安装推力计算程序,计算出飞机发动机的安装阻力和安装推力。经过与实际安装推力的资料数据对比,说明本文所用方法误差小,并且充分考虑了不同飞机之间的差异,比一般利用经验公式的方法更可靠。     

轴对称喷气式飞机后体减阻优化设计

本文通过数值计算二维轴对称雷诺平均N-S方程,采用-ωSST模型对喷气飞机后体—尾喷管的内外场进行冷态数值模拟,运用遗传算法以后体阻力系数为设计目标,在飞行马赫数Ma∞=0.6的情况下,对后体外形进行优化设计。将计算结果与实验结果进行对比,验证了在特定后体外形条件下,喷流参数变化对飞机气动性能的影响。同时,计算结果也表明优化外形能够保证后体流动不发生分离,从而降低了后体阻力。     

运输机机翼外形和吊舱位置一体化优化方法

针对运输机翼吊布局中机翼-吊舱之间有可能存在较大干扰阻力的问题,本文提出一种机翼扭转角和吊舱位置的同时优化的设计方法。气动分析模型采用高阶面元法程序PANAIR,应用参数化方法和CATIA二次开发技术实现飞机三维外形的自动生成;应用网格生成工具Gridgen的脚本语言实现气动网格的自动生成。为了减少调用气动分析模型的次数,在优化中应用了代理模型方法。算例的优化结果表明:优化后机翼的展向升力分布更加合理,显著地降低机翼和吊舱间的干扰阻力。     

航空连接结构的参数化有限元分析系统设计与开发

本文应用弹塑性有限元数值计算及参数化设计方法,对飞行器典型连接件结构的有限元计算进行了参数化设计,并基于ANSYS平台开发了参数化的连接件结构有限元仿真系统,借助通用有限元软件ANSYS的参数化设计语言APDL与界面设计语言UIDL,完成参数化建模的宏文件设计与参数输入的界面设计及参数传递,初步实现了有限元建模、计算和后处理的全程参数化;实现了对ANSYS的二次开发;利用有限元的子模型技术,在计算开销相对较小的情况下保证了计算精度,从而形成专业化、参数化、高效率的连接件有限元仿真设计系统。文中给出了采用该系统与传统方法比较的算例,可以看出结果有较好的吻合。     

一种高效的壁面距离计算方法

采用CFD(Computational Fluid Dynamic)技术进行湍流效应数值模拟时,经常需要计算流场网格点到最近壁面的距离。当网格规模很大,特别是针对一些动网格问题,壁面距离的计算量很大且费时,影响整个流场求解效率。本文对壁面距离计算的直接算法和基于二叉树技术的ADT(Alternating Digital Tree)快速算法进行了对比分析和研究,提出了一种高效的、基于方盒切割技术的快速计算方法,并采用不同外形的CFD计算网格对该方法进行了验证。计算结果表明,新方法的壁面距离计算效率明显高于直接算法和ADT算法,并具有较好的鲁棒性和通用性,可以便捷地移植到现有的CFD计算程序中,从而提高整个流场的数值计算效率。