基于气动力降阶的弹性飞机阵风响应仿真分析及验证

低速飞机在阵风作用下容易产生非线性气动力,从而引发非线性气动弹性效应,对飞行安全造成威胁。针对此类问题的分析,经典面元法无法满足计算精度要求,计算流体力学(CFD)/计算结构动力学(CSD)全阶耦合分析效率低下,因此需建立满足工程应用的高精度、高效率的飞行动力学仿真分析模型。针对以上问题提出了一种适用于工程的非线性气动力降阶模型(ROM)用以实现弹性飞机飞行动力学仿真,特别是低速飞机在遭遇大幅值阵风情况下的阵风响应仿真。以风洞试验飞翼飞机模型为对象,利用CFD方法获得了该模型的气动力数据,利用自回归移动平均(ARMA)方法和径向基函数(RBF)神经网络方法分别建立了该模型的线性气动力ROM和非线性修正气动力ROM。结合模型的刚弹耦合飞行动力学方程对模型遭遇阵风情况下的响应进行仿真分析,并将仿真结果和风洞试验结果及CFD/CSD计算结果进行对比。结果表明建立的基于非线性气动力ROM的弹性飞机仿真模型在气动力预测、稳定性分析及阵风响应分析方面的表现都优于基于线性气动力ROM的仿真模型,和试验结果及CFD/CSD分析结果一致性较好,且所建模型在相同工况下的仿真时间远低于CFD/CSD分析方法...     

经典跨声速翼型RAE2822数据分析

经典跨声速翼型RAE2822风洞试验数据长久以来被广泛用于CFD计算方法和软件的验证与确认，但是数据的正确使用或者说合理使用仍存在一些需要研究和注意的问题，包括计算网格、风洞试验数据修正、中弧线修正、翼型几何定义和建模，以及摩擦阻力系数和边界层速度剖面的原始定义等。在开展CFD研究之前，必须首先对计算方法进行验证，尤其是要先尽可能消除计算结果对计算网格的依赖性；经过对目前可开放使用的计算网格的不足之处进行分析，研制了一套高品质的计算网格并获得了预期的一阶网格收敛性；通过计算软件交叉验证，进一步确保所用计算软件的可信度。在将CFD计算结果与翼型风洞试验数据进行比对时，通常需要对马赫数和攻角进行修正，且如何修正是一个需要持续研究的问题；翼型中弧线修正是一种有效的方法，但需要考虑流动参数的影响。原始翼型几何构型采用有限离散点定义，计算网格生成过程中需要采用插值方法布置型面网格点，不同插值方法对于翼型前缘附近流动的数值模拟有细微影响。多数相关研究工作只比对压力分布；少数研究工作会比对摩擦阻力系数、边界层及尾迹速度剖面，但在比对时，需要注意原始风洞试验相关参数定义与CFD计算常用定义的区别。     

改进的综合CFD与畸变合成法预测动态畸变(英文)

对综合CFD与畸变合成方法进行了改进。通过CFD求解,获得AIP上各稳态气动参数。采用BP人工神经网络直接求解高度非线性的湍流相关方程,该方程将CFD求解的稳态气动参数与试验测得湍流度进行关联。经过训练的网络可根据稳态CFD结果预测湍流度。之后根据CFD计算的稳态总压和神经网络预测的湍流度合成动态压力。最后根据合成的动态压力寻找最大瞬时畸变。文中使用了六个工况的飞行数据验证该方法的准确度并且与原方法的计算结果进行了比较。结果表明该方法预测的畸变值与试验畸变值吻合较好,且预测结果较原方法与试验值更接近,该方法能够在试验前用来预测动态畸变,具有很强的工程价值。     

尖前缘平板的高超声速粘性干扰效应引起的压力分布规律研究

采用CFD方法计算了多种来流条件下高超声速平板绕流的粘性干扰现象,计算与实验结果吻合良好。针对理论分析在平板前缘失1效的情况,结合理论分析、CFD结果和试验结果对原有沿流向壁面压力分布关系式进行了修正,使得修正后关系式在整个平板上与相关数据吻合,显著增加了压力分布关系式的应用范围。同时根据CFD计算结果,分析了平板附面层内压力等流动参数的分布规律,发现附面层内法向压力分布存在峰值,壁面压力分布大于来流静压,并给出了峰值压力及其位置的拟合式。     

伴随气泡和气穴低压管路瞬态的建模与分析

在航空航天液压系统中,为了合理地预测伴随气泡和气穴低压管路压力瞬态脉动,给出了用来描述管路流动瞬态特性的数学模型.采用有限差分法实现了偏微分方程中空间域的偏导,并建立了计算有效体积弹性模量、摩擦阻力项、气泡和气穴体积的数学模型.利用遗传算法对模型中三个参数即油中的初始气泡体积、气体析出时间常数和气体溶解时间常数进行参数辨识,得到了参数辨识以后的低压管路压力瞬态脉动模型.通过对仿真结果和实验结果的分析比较,表明该低压管路瞬态模型用于计算伴随气泡和气穴的压力瞬态脉动是可行的.      

充气式蓄能器参数选型方法研究

容积和预充气压力是充气式蓄能器的关键参数，对其动态性能（工作状态下，其压力与流量的输出性能）有较大影响。以某型飞机舵面应急型充气式蓄能器计算选型为例，基于Hypneu 仿真软件中的液压元件库，结合该型飞机液压系统的主要设计参数以及应急型蓄能器的工作原理，建立蓄能器作动仿真模型；在Hypneu中设置测试点记录仿真曲线，通过分析仿真曲线探究充气式蓄能器的关键参数对其动态性能的影响；结合理论计算与仿真分析，提出一套蓄能器参数选型方法，并通过试验验证该选型方法。结果表明：蓄能器的预充气压力越大，初始蓄能器输出的有效容积越小，减少预充气压力增加有效容积的同时会降低蓄能器的输出压力，试验结果与仿真结果接近，使用该方法选型的蓄能器可以满足作动要求。     

基于不同相变模型的微重力液氢贮箱内压力与温度分布仿真预示

受空间热流的作用,相变是低温推进剂在轨压力控制中需要被考虑到的影响因素。为研究液氢贮箱内的流体行为特性,建立了低温流体CFD仿真模型,对于相变过程,基于不同的相变传质理论,建立了四种相变仿真模型。根据NASA开展的AS-203液氢贮箱压力上升试验数据,对封闭贮箱内压力上升和温度分布开展仿真预示,分析了不同相变仿真模型对压力上升和温度分布预示的结果。结果表明,相变模型1和相变模型3得到的压力上升速率和温度场结果与试验结果趋势较为一致。受到算法和适用性的影响,相变模型2和相变模型4对AS-203液氢贮箱的温度预示偏差较大,相变模型4对压力上升的预示偏差较大。     

防热试验用高温超声速燃气流场热环境分析

对用于防热系统试验考核的高温超声速燃气流场的流场结构和热环境特性进行了计算和分析.将CFD模拟方法与工程计算方法相结合,对两类锥形喷管形成的燃气流场进行了数值模拟计算和流场结构分析,并与试验过程中高速拍摄的图片进行了比对.利用模拟得到的流场气动参数,计算了沿流场轴向位置球头驻点模型的热流和压力,并与试验测量数据进行了对比,结果吻合较好.锥形喷管形成压缩波和膨胀波交替的流场结构,导致流场气动参数与模型热流和压力形成正弦变化的特点.     

两相湍流的非线性k-ε-kp模型

由两相湍流的代数应力模型，提出一种两相湍流的非线性k-ε-kp模型，可以合理地模拟各向异性较强的旋流，浮力流等两相流动，具备二阶矩模型的优点，又比二阶矩模型简单。文中推导出气相、颗粒相的雷诺应力和两相脉动速度关联的非线性应力应变关系式。这些代数式和两相各自的湍动能k，kp以及两相脉动关联湍动能kpg的方程联立，就构成非线性k-ε-kp模型。本文将该模型用于模拟突扩湍流两相流动，给出两相时均速度场及雷诺应力各分量，并且将模拟结果和PDPA实验数据以及二阶矩模型的模拟结果对照。结果表明，该模型预报各向异性两相湍流的能力和二阶矩模型的能力接近，但是计算量比二阶矩模型的小得多。     

两相湍流的非线性k-ε-kp模型

由两相湍流的代数应力模型，提出一种两相湍流的非线性k-ε-kp模型，可以合理地模拟各向异性较强的旋流，浮力流等两相流动，具备二阶矩模型的优点，又比二阶矩模型简单。文中推导出气相、颗粒相的雷诺应力和两相脉动速度关联的非线性应力应变关系式。这些代数式和两相各自的湍动能k，kp以及两相脉动关联湍动能kpg的方程联立，就构成非线性k-ε-kp模型。本文将该模型用于模拟突扩湍流两相流动，给出两相时均速度场及雷诺应力各分量，并且将模拟结果和PDPA实验数据以及二阶矩模型的模拟结果对照。结果表明，该模型预报各向异性两相湍流的能力和二阶矩模型的能力接近，但是计算量比二阶矩模型的小得多。     

舰载燃气轮机间冷器传热与流动的数值模拟

针对平直型翅片板翅式间冷器矩形通道的结构特点,建立了流动换热分析的耦合计算模型,应用计算流体动力学方法对间冷器的通道流场进行了数值模拟.分析了计算区域内多个截面的温度、压力、局部传热系数等参数的分布图形和变化趋势,并考察了不同工况下间冷器的换热性能和压力损失.所得结论可为舰载燃气轮机间冷器的设计与优化提供有益的参考.      

低温推进剂贮箱压力变化的CFD仿真

为预示低温推进剂贮箱在地面停放阶段的压力变化并研究贮箱内物理过程的相互作用关系,建立了包含液体推进剂和混合气体两相的二维轴对称volume of fluid(VOF)计算流体力学(CFD)模型,并引入了基于热力学平衡假设的推进剂相变模型.对实验液氢贮箱进行仿真得到的压力上升速率与实验结果相差9.1%.通过对地面加压停放阶段下的液氢和液氧贮箱的仿真发现:造成液氢贮箱压力上升的主要因素是壁面漏热对气枕的加热作用,而液氢蒸发影响更小,液氧贮箱在加压停放阶段初期明显受到液氧相变的影响.两个贮箱中液面附近的对流运动在不同的气液传热过程作用下有不同的变化趋势,对流运动会影响推进剂的相变进而影响贮箱的压力变化.      

Pressure wave propagation characteristics in a two-phase flow pipeline for liquid-propellant rocket

The aim of this paper is to investigate the pressure wave propagation behavior in LOX/kerosene rocket engine pump pipeline. A pressure wave propagation model for gas–liquid two-phase flow with phase change taken into consideration is first proposed by using ensemble-averaging techniques. Then condensation of gas oxygen in subcooled liquid oxygen and the corresponding mixing process in pump pipeline are numerically simulated with the application of thermal phase change model in Computational Fluid Dynamics code CFX. Finally, based on the established model and the predicted flow parameters, pressure wave in gas–liquid pipeline flow, characterized by the propagation speed and the attenuation coefficient, is studied by personally compiled code on MATLAB. The calculation results indicate that the propagation speed of two-phase pressure wave first slightly falls and then increases with the decreasing void fraction along the flow direction, and the increase of angular frequency causes the propagation speed and the attenuation coefficient of two-phase pressure wave to increase. The vapor bubble diameter, mass flow rate and inlet temperature of gas oxygen are shown to have strong effects on the mixing flow condensation process of gas oxygen and liquid oxygen. Bubble with larger diameter can weaken the condensation process, and leads to greater cross-sectional mean void fraction of gas oxygen. This phenomenon results in an overall decrease in pressure wave propagation speed and an increase in attenuation coefficient except for the short slug-flow region located at the pipeline inlet. Similar effects can also be achieved by increasing the mass flow rate or the inlet temperature of gas oxygen.     

混合式动压气体轴承结构优化与可靠性分析

通过将径向、止推螺旋槽动压气体轴承相结合,建立了混合式动压气体轴承的润滑分析模型。阐述了其结构特点与润滑机理,建立轴承无量纲稳态Reynolds控制方程。提出混合式动压气膜压力耦合计算方法,推导气膜压力差分表达式,定义边界条件,构建气膜压力分布的数值计算方法。以最大径向承载力为目标优化结构参数。基于最优结构参数建立轴承气膜有限元模型,运用CFD分析轴承转子系统受不同冲击载荷时径向稳定性变化规律,研究混合式动压气体轴承动态特性与可靠性。搭建混合式动压气体轴承试验台,验证数值计算方法和有限元仿真分析的正确性。结果表明:提出的压力耦合计算方法可以准确地计算分析稳态气膜压力分布、承载力和承载性能,有限元仿真能更好地模拟动态流场变化,计算分析动态承载力、动态特性系数和稳定性。高转速下混合式气体轴承承载力、稳定性较好,对单向阶跃力、单向矩形力的抗冲击能力强,可靠性强。混合式动压气体轴承在优化承载性能与刚度的同时,应考虑抗冲击特性和稳定性以提高轴承的综合性能。      

硝酸羟胺基单组元液体火箭发动机起动过程的模拟（投稿 空间推进专刊）

无毒单组元液体火箭发动机是空间推进的发展方向之一。本文针对硝酸羟胺基单组元液体火箭发动机的起动过程,建立了零维模型、不考虑相变过程的仿真模型,和考虑相变过程的仿真模型,并采用这些模型对60 N发动机进行了起动过程的模拟和对比分析。计算结果显示,推力室升压过程历经快速的气体充填和较慢的催化室升温两个阶段。零维计算模型和不考虑相变的CFD方法计算得到的两个阶段时长基本一致,但第一阶段时间显著低于试验结果,第二阶段时长与试验符合。考虑推进剂相变过程的模型计算结果与试验结果符合较好。一维计算得到HAN基推进剂在进入催化床约4mm的长度内完全分解。升压第一阶段受催化分解影响很大,HAN基推进剂的分解反应速率低于肼的分解反应速率。     

污垢沉积影响叶片表面换热的研究

选用了两个实际透平叶片MARKⅡ静叶和某高压透平静叶,前者为光滑叶片,后者为表面覆盖有圆台的粗糙表面叶片.利用CFD数值模拟技术对叶栅流道特性和叶片表面换热分布进行数值模拟计算,数值计算结果与实验结果对比验证了数值方法的可靠性.研究表明,准确预测转捩位置和选取合适的y+值是计算叶片表面换热特性的关键,剪切应力输运(SST)转捩模型预测的结果最令人满意.在粗糙叶片表面数值计算中,对表面粗糙度模型进行验证,然后详细分析了污垢沉积对叶片表面换热的影响.结果表明:随着叶片表面粗糙度的增大,叶片近壁面湍流增强,叶片的表面传热系数增大.      

导弹气动特性在亚跨音速下的风洞试验研究

在导弹的设计过程中,导弹的气动特性作为重要因素直接影响导弹飞行的动态品质。在亚跨音速段气动特性呈现剧烈非线性的情况下,工程估算以及CFD数值计算方法所能提供的气动计算精度有限,导致对舵效特性的辨识精度较低,需要进一步采用风洞试验的方法精确计算气动参数,进而确定导弹的舵效。本文应用风洞试验方法研究导弹飞行马赫数在亚跨音速段对导弹气动特性的影响。研究结果表明：亚音速时导弹的气动特性基本一致,跨音速时发生剧烈的非线性变化;导弹的俯仰舵效先增加后减小,滚转舵效先减小后增大。结论对导弹控制律的设计以及后续的工程型号研制有参考价值。     

涡轮叶片内部通道流动换热算法改进

在使用工程计算方法对涡轮叶片温度场进行计算时,往往将叶片内流通道简化成光滑或带肋的换热管元件,容易忽略 各内流管段之间的影响,造成计算得到的叶片3维温度场与真实温度场存在较大差异。针对上述问题,为了提高对涡轮叶片3维 温度场模拟的准确度,对涡轮叶片内流通道的换热流动算法进行改进。考虑涡轮内部蜿蜒通道中弯转区和弯转效应2种因素对 涡轮内部流动换热的影响,使用试验得到的数据对2种因素影响区域的换热情况进行修正,利用修正后的算法对某工作叶片进行 温度场计算,并对修正前后叶片温度场进行了对比分析。结果表明：采用修正后算法得到的蜿蜒通道内的气体温度相较于修正前 算法得到的沿程升高更多,修正后算法求得的叶片整体平均温度降低,最大温差增大。     

民机空调组件舱通风仿真计算研究

民机的空调组件舱位于非增压区,其内部包含有换热器、压缩机等高温环控组件。为避免高温对舱体结构产生危害,需要采用冷却措施。以某机型的空调组件舱通风冷却方案为参考,建立了空调组件舱通风系统的几何模型,并以此模型为基础进行了CFD仿真计算。对计算结果进行分析后,提出了方案的缺陷,并为后续的方案设计提出了关注点。     

飞行器阵风响应的CFD-6DOF仿真分析

飞行器在大气中飞行,不可避免地受到阵风的影响。阵风所附加的气动载荷引发飞行器飞行状态的改变,过大幅值的阵风影响飞行的性能与安全。针对这种状况,首先采用改进的Lamb-Ossen涡模型,建立尾涡形式的阵风场;然后采用基于CFD技术的非定常N-S方程求解,并在计算网格中引入"网格速度"来模拟阵风,对SWIM(Subsonic Wall Interference Model)尾涡中的定常气动特性进行验证;最后通过CFD-6DOF的耦合,对SWIM俯冲穿越尾涡场的飞行轨迹进行研究。结果表明:计算结果与实验值符合较好;SWIM在尾涡中飞行时出现抖动、下沉、改变飞行状态、剧烈翻转的现象,与实际飞行器进入尾涡中的轨迹特性类似。