超声速化学反应流场迎风格式数值模拟及并行计算

基于有限体积迎风格式对超声速燃烧流场进行了的数值模拟.由于超声速燃烧流场绕流的复杂性,要求对多组分Euler/N-S方程求解的数值模拟方法应具有较高的计算精度及效率.本文引用辅助点方法建立了具有空间二阶精度的van Leer迎风矢通量分裂格式,并应用于超声速燃烧流场绕流的数值模拟.化学反应为氢气/空气十反应模型,采用考虑了化学反应特征时间的当地时间步长显式Runge-Kutta时间推进格式.对钝头体模型爆轰现象、后向台阶氢气喷射及二维内外流超声速燃烧流场模型进行了区域分裂技术的并行计算.计算结果与参考文献作了对比,得到了满意的结果.     

飞机座舱围护结构瞬态温度分布的数值计算

针对客机座舱围护结构的瞬态传热问题,论文以ARJ21-700型飞机为例,通过对飞行包线内的环境条件、座舱围护结构布置等内外部因素的分析,建立了相应的数学计算模型,并采用有限差分法进行温度场求解,获得了飞行包线内的飞机座舱围护结构各层温度随时间及飞行高度的变化情况.在此基础上,考察了材料层布置方式对座舱围护结构温度分布的影响,研究工作对于飞机座舱瞬态热载荷的计算具有参考价值.     

基于混合网格上的高温非平衡流的数值模拟

新一代再入飞行器及空间传输系统需要了解其飞行时的气动问题及热力学问题。数值模拟高空、低密度、高熵及非平衡流动是一具有挑战性的问题。本文针对带座舱飞船高超声速再入大气过程中存在的严重气动加热现象,利用混合网格及Osher逆风格式,数值求解了三维化学非平衡Navier-Stokes方程,其中化学模型是5组分1 7个化学反应的空气化学模型,对带座舱飞船再入高度为40 km和马赫数为2 0 ,1 0的化学非平衡流场进行了数值模拟,给出了迎角为0°和2 0°情况下的各个组分的密度分布、压力等参数,并与量完全气体的计算结果进行了比较。     

空气导管对飞机密闭空间传热流动的影响

在分析高温高压空气导管作用下飞机货舱三角区内的空气传热和流动机制的基础上,采用数值计算方法对舱内空气的导热-辐射-对流耦合换热问题进行研究,考察了空气导管布置位置、导管管径和辐射对传热流动特性的影响。研究结果表明:舱内空气传热和流动受空气导管在垂直方向上的位置变化影响较大;当导管布置位置在水平方向上变化时,仅对温度场影响较大,而对流动强度的影响较小;随着导管直径的增加,空气导管与舱内空气的总换热量增加,舱内各处空气的温度也随之升高,相应的自然对流强度减弱;辐射传热会减小空气导管与货舱三角区壁面之间的换热温差,增加总的换热量,并且削弱自然对流的强度。     

一种高温瞬时灭菌民用飞机空调系统（英文）

飞机座舱是一个狭小的封闭空间环境,空调系统的供气品质关系着舱室乘员的健康。为了保障新型冠状病毒(2019-nCoV)疫情期间飞机座舱内乘员的健康安全,本文提出了一种新型的高温瞬时灭菌空调系统(Ultrahigh-temperature instanteous sterilization air conditioning system,UHT-ACS)。在此基础上,对不同飞行状态下的UHT-ACS系统和传统飞机空调系统(Air conditioning system,ACS)进行了仿真计算。通过分析可知,UHTACS系统中的座舱再循环空气与高温高压的发动机引气相混合后,混合空气温度可达148.8℃,从而实现高温瞬时灭菌的目的。当热天UHT-ACS系统进行制冷时,系统送风温度为12℃左右,当冷天UHT-ACS系统进行加热时,系统送风温度为42℃左右,满足飞机空调系统的要求。通过与传统飞机空调系统的对比可知,当系统新风量为60%,再循环空气为40%时,本文提出的新型UHT-ACS系统供气温度与传统飞机空调系统供气温度一致,并且UHT-ACS系统的涡轮出口空气温度要高于传统飞机空调系统,这有助于降低涡轮出口结冰的风险。因此,本文提出的新型UHT-ACS系统可满足飞机不同飞行状态下的空气调节需求,能够为飞机座舱提供安全供气。     

三维典型超声速燃烧室流场数值模拟

超声速燃烧是实现以吸气式冲压发动机为动力进行高超声速飞行的关键技术之一.本文对三维超声速燃烧室流场进行了基于有限体积方法的并行数值模拟,还对JPNal三维典型超声速燃烧室内的激波诱导燃烧流场进行了数值模拟.控制方程为三维多组元带燃烧反应模型的Euler方程,空间离散采用二阶精度VanLeer迎风通量分裂格式,时间推进为考虑了化学反应特征时间的5步Runge-Kutta方法,燃烧反应为氢气/空气十反应模型.为提高计算效率,采用了基于区域分裂技术的并行计算.最后对比分析了单、双气态氢气喷流对燃烧室内流场结构及燃烧特性的影响.     

有限元法在光学侧窗热防护设计中的应用

光学侧窗热环境极为严酷,可采用液体对流内冷却技术对光学侧窗进行热防护。以有限元法对光学侧窗温度场进行了分析研究,研究结果表明:该方法计算精度可靠;液体对流内冷却技术能有效起到降温的作用。同时还定量的给出了冷却液用量对结构温度场的影响,为设计的经济可行性提供了保障。     

基于狭缝湍射流多尺度特征研究的机舱条缝送风口内结构的优化设计

针对 MD-82真实客机座舱内现有的壁面条缝型送风口流场平均流速沿座舱轴向分布的不均匀性,提出了一种送风口条缝形芯内结构的设计,使条缝型送风口处气流混合更加均匀,湍流发展更充分,平均流速沿座舱轴向分布更加均匀。并利用热线测速技术,精细测量条缝型送风口流场进行验证。从吹风感舒适性要求出发,利用子波分析,对条缝型送风口流场不同位置的瞬时速度时间序列信号作多尺度湍涡成分分析。结果表明:加装了条缝形芯内结构的条缝送风口流场的气流脉动特征频率始终维持在21.83Hz 左右,而导致人最不舒适的空气速度脉动频率为0.2~0.6Hz,从而证实加装了条缝形芯内结构的条缝送风口流场的气流提高了人体感觉的舒适性。     

飞机空调系统对座舱微环境空气质量和热舒适性的影响综述

对近20年来的飞机座舱空气调节方法进行了综述.首先,对飞机座舱微环境下的空气质量、污染物扩散、热舒适性等进行了分析评估.然后,阐述了飞机座舱空气环境研究的实验和数值方法.最后,进一步分析了相关研究方法的优缺点,并对不同研究方法的前瞻性研究对象和新一代环控系统进行了展望.     

基于CFD的直升机旋翼流场及气动力计算

建立了一个基于结构运动嵌套网格的流场求解器,用来精确模拟复杂的旋翼流场,为更好地预测旋翼气动载荷提供一套计算方法.在该求解器中,控制方程为惯性坐标系下的三维非定常Navier-Stokes方程,空间方向上采用低数值耗散的Roe格式结合三阶逆风格式(Monotonic upwind scheme for conservation law,MUSCL),湍流模式采用了一方程的Spalart-Allmaras模型.应用所建立的方法,分别针对Caradonna模型旋翼、UH-60A直升机旋翼的悬停流场和7A旋翼、SA349/2直升机旋翼的前飞流场以及旋翼气动载荷进行了数值模拟.计算结果表明,本文的方法在一定程度上提高了对旋翼流场的模拟能力,进而提高了旋翼气动载荷的计算精度.     

大迎角三角翼旋涡运动及其破碎特性的数值研究

从流体力学的基本方程出发，利用Ｈａｌｌ的涡核准柱假设，导出反映涡核运动的Ｎ－Ｓ方程。采用差分方法计算旋涡流场，进而分析三角翼上前缘分离涡的运动特点及其破碎机下。从计算结果可以看出，旋涡的轴向速度向下游逐渐下降，且涡心处于降较快，外缘下降较慢，反映了粘性作用自涡心外缘逐渐下降的特点；涡核外缘的径向速度开始为负，说明开始阶段有流体流入涡核，随着旋涡向下游运动，径向速度有所增加，到一定位置后增加迅速，说     

NUMERICAL SIMULATION OF THREE DIMENSIONAL INTERACTING TURBULENCE FLOW FIELD OVER PROJECTILE WITH LATERAL JETS

INTRODUCTIONIn recent years,the study of lateral jets hasaroused more interest of aerodynamicist all overthe world.On the one hand,the interacting flowfield caused by lateral jet includes complicatedshock wave patterns and has important values tostuding flow phenomena.It has been not veryclear to this day and it is very difficult for nume-rical simulation.On the other hand,lateral jethas important applications for the field of aero-nautics,astronautics and weapons:for example,lateral jet …     

槽道湍流相干结构的研究

采用高阶差分方法求解非定常不可压Navier-Stokes方程，进行槽道流相干结构的数值研究，与常用的谱方法相比，不需限制边界条件，可用于更一般的粘性流体的流动问题，文中所用时间分裂法和三维耦合高阶差分格式适用于包括邻近边界点在内的整个计算区域，克服了三维各自用高阶中心差分格式不适用邻近边界点的困难。通过算例模拟了槽道流动中单个相干结构的演化发展规律，结果与实验一致。     

基于格子Boltzmann方法的可压缩翼型绕流模拟

采用能够恢复可压缩Navier-Stokes方程的耦合双分布(DDF)格子Boltzmann方法(LBM)模拟了贴体网格下的NACA0012翼型绕流.首先在Ma＝0.5,α＝0.0,Re＝5 000状态下通过模拟低雷诺数流动检验了该方法;然后分别在Ma＝0.5,0.85和1.2,α＝-0.05,1.0和0.0状态下进行了模拟,并在Ma＝0.85时通过加密网格得到了更好的结果.为了避免转捩问题采用了无粘边界条件进行模拟,但得到的压力系数与实验结果吻合良好,证明了该方法的可靠性并显示了其模拟可压缩的潜力.该方法在翼型绕流上的应用为湍流流动的数值模拟提供了基础.     

微分求积单元法分析变曲率井内受径向约束管柱的非线性屈曲

基于Love的空间弯扭杆平衡方程,通过引入井壁约束条件导出了平面变曲率井内受径向约束管柱的平衡方程。采用微分求积单元和增量迭代法求解了曲率线性变化的井内管柱的非线性屈曲问题,通过与有限元计算结果的对比验证了所构建方法和编写程序的正确性,而且还表明微分求积单元法有方法简单、易于实施,计算量少、精度较高等优点。计算结果表明,等曲率井内管柱屈曲的临界载荷明显大于曲率线性变化井内管柱屈曲的临界载荷,另外,变曲率井眼的曲率变化对管柱弯矩、井壁约束力有显著的影响。     

旋翼悬停性能的Re数影响研究

为了解决由旋翼模型风洞试验数据预估全尺寸旋翼性能的问题,开展了旋翼模型悬停性能试验数据的Re数影响及修正方法研究.采用3种旋翼的不同尺寸的模型进行了对比试验,通过试验结果的分析,研究了胁数对旋翼悬停效率的影响规律;按照旋翼理论的有关公式和定义,参考粘性流体力学的研究结果,推导了对悬停效率进行Re数效应修正的公式;利用推导的公式,对模型试验结果的悬停效率进行了修正计算和比较.对比试验结果表明,Re数较小的旋翼模型的悬停效率相对较低;修正计算结果表明,利用笔者推导的修正公式,可以对旋翼模型试验的悬停效率进行合理的Re数影响修正.     

7°尖锥高超声速边界层转捩红外测量实验

为了推动高超声速边界层转捩研究的深入开展,给边界层转捩机理研究、物理模型验证、转捩数据库构建和转捩天地相关性的建立等提供基础风洞实验数据,在中国空气动力研究与发展中心的Φ1 m高超声速风洞开展了边界层转捩规律红外热图实验。针对半锥角7°尖锥模型,研究了不同单位雷诺数、迎角和马赫数对尖锥边界层转捩位置的影响规律。实验单位雷诺数（0.49~2.45）×107/m,迎角范围-10°~10°,马赫数5~7,模型头部半径0.05 mm。通过红外热图技术测量模型表面温度分布,获得了较为详细的转捩位置和转捩参数影响规律。实验结果表明:在马赫数5~7范围内,马赫数增大,尖锥转捩位置提前,分析认为是高马赫数时的雷诺数较大、自由流噪声水平较高引起;随着单位雷诺数的增大,边界层转捩位置前移,转捩雷诺数保持不变,约为3.0×106;小迎角时,随着迎角的增大,迎风面边界层转捩推迟,背风面边界层转捩前移,在10°大迎角时,迎风区中心线转捩前移,出现迎角"转捩逆转"现象,背风区出现了流动分离导致的低热流条带。     

用预处理方法及可压Navier-Stokes方程计算低马赫数流动

用可压方法及预处理方法求解了低马赫数紊流流动,其中空间离散格式采用了Roe格式,紊流模型采用了DES模型.以绕NACA0012的紊流流动为例,分析了该方法对提高流动速度和精度的特性,并将其用到多段翼型及绕圆柱的非定常流动中.数值模拟结果表明,该方法能有效求解低速流动及混合流场问题,且具有较高的数值精度.     

低雷诺数效应对0.6m连续式风洞性能影响

风洞试验中通常采用降低运行总压的方法来扩大风洞雷诺数模拟的下边界。同常压试验相比,在低雷诺数条件下,风洞流场是否存在明显变化,风洞流场品质是否满足指标要求,直接影响风洞试验数据的精准度。为了研究低雷诺数效应对0.6m连续式跨声速风洞性能（包括轴流式压缩机性能、总压及马赫数控制精度、流场均匀性）的影响,调试人员在0.6m连续式跨声速风洞中开展了大量相关试验,本文在对试验结果进行整理、分析的基础上,给出了Re对风洞性能的影响。结果表明:（1）Re对压缩机性能、总压控制精度、马赫数控制精度、流场均匀性都有明显影响,当Re_c<5×105（c=0.1√A）时,雷诺数效应明显,且Re越小,影响越大。（2）0.6m连续式跨声速风洞能够真实、准确反映Re对风洞流场性能和测力试验数据的影响规律,是开展高空低雷诺数飞行器、翼型、发动机等性能研究的理想地面模拟试验平台。