干线客机机翼气动/结构综合设计研究

为了妥善处理机翼设计中气动和结构之间的矛盾,提高飞机的整体性能,对干线客机机翼气动/结构综合设计方法进行了研究.将非等熵全位势方程CFD方法与基于工程梁理论的机翼结构设计方法相结合进行分析计算,用序列二次规划法对全局敏度方程所构造近似系统进行优化,通过对一系列近似系统的优化,逼近综合设计的最优解.最后以一种干线飞机的机翼气动/结构综合设计为算例验证了这一设计方法的可行性.     

求解跨音位流问题的改进SIP方法

本文研究Sankar提出的,用全位势方程求解跨音流场的强隐式方法(Strongly Implicit Procedure)。针对强隐式方法在求解跨音流场中的问题作改进,即采用了新的人工压缩密度计算方法,计算绕孤立机翼、叶栅通道的跨音流场,降低了在跨音位流计算中激波前后的数值振荡,并提高了算法的稳定性,     

二维翼型绕流的CE/SE方法

采用时空守恒元和解元CE/SE(space-time Conservation Element and Solution Element method)法,求解二维Euler方程,开展了翼型绕流的无粘数值模拟研究.用非敏感克朗数计算格式消除克朗数过小引起的数值耗散对解的污染,结合当地时间步长法,解决网格不均匀引起的当地克朗数变化跨度大的问题.对NACA0012翼型的无激波流场进行了二维数值模拟,并与AGARD算例做了对比.结果表明:CE/SE方法的计算结果与AGARD结果吻合得很好,为该数值计算方法对翼型绕流数值模拟的进一步应用奠定了基础.      

高超音速化学非平衡流的隐式格式迭代解法

针对高超音速化学非平衡流计算中的隐式TVD格式给出了一种迭代解法.该方法是在原有的Harten-Yee格式基础上,在每一时间步对无粘部分进行迭代运算,目的是为了加速收敛.借助对模型方程的收敛性分析,证明该迭代算法是收敛的.二维数值实验表明该方法对于高马赫数时的化学非平衡流能明显提高隐式格式的计算效率,其计算结果是满意的.     

用GAO-YONG湍流模式对翼型跨音粘流的数值模拟

应用GAO-YONG可压缩湍流模式,数值模拟了NACA0012,RAE2822翼型的定常跨音速粘流算例.对流项采用三阶ROE格式,扩散项采用二阶中心格式,用多步Runge-Kutta显式时间推进法求解空间离散后的控制方程.计算结果很好地预测了翼型表面的压力系数的分布、激波的位置、马赫数等值线的分布等情况,并且对翼型表面激波与边界层相互干扰以及层流向湍流的转捩问题进行了分析计算.计算结果与实验值符合很好,表明GAO-YONG可压缩湍流模式应用合适的计算方法能够高精度模拟翼型跨音粘性流动问题,并且基于GAO-YONG可压缩湍流模式各向异性湍流粘性的机理,提供了一种预测转捩起始位置的判别方法.      

翼身组合体的变可信度优化设计

基于信赖域思想,发展了一种变可信度优化设计方法组织和管理不同可信度模型,有效利用低可信度模型建立近似模型,进行优化设计,高可信度模型仅修正优化近似模型,使最终优化结果收敛到高可信度模型上.采用全位势边界层有粘无粘迭代方法作为低可信度模型,多重网格Navier-Stokes方程方法作为高可信度模型,对翼身组合体进行优化设计,优化结果表明:发展的变可信度优化方法优化效率高,且能保证优化精度,具有实际的工程应用价值.     

数值模拟高超音速化学非平衡流的隐式格式

针对高超音速化学非平衡流提出了两种隐式格式.一种是近似因子分解算法对角化形式在隐式TVD(全变差递减)格式上的推广,另一种可视为隐式TVD格式的LU分解形式.数值试验表明,本文给出的算法保持了TVD格式的高分辨率特性,同时计算量大大减少,是高超音速化学非平衡流数值模拟的有效方法.     

非热电子束流在耀斑环中能量的沉积

本文计算了非热电子束流通过库仑碰撞和反向电流过程在均匀截面的耀斑环中的能量沉积,讨论了非热电子束流对耀斑软X射线等离子体的加热和耀斑非热模型的能量收支平衡问题。      

一种利用正交流线坐标系求解跨音速管流的新方法

本文以正交流线坐标系为工具,采用位流模型,导出了正交流线坐标系下平面或轴对称流动的控制方程组,利用流函数与位函数的存在性,进而导出了两个首次积分。设计了一套求解匀熵管流的分凝迭代方案,模拟了管内无激波的平面或轴对称喷管的跨音速流动。数值实验表明该方法能够处理小喉道曲率半径平面或轴对称喷管跨音速流动,与已有的数值计算结果及实验结果符合较好。     

带副翼的翼身组合体绕流的 Euler方程解

提出了一种求解带副翼偏转的翼身组合体绕流的Euler 方程计算方法.将对接分区网格与分区求解算法相结合,有效地求解了绕此外形的复杂流动.提出了一种满足通量守恒的内边界耦合条件.数值方法中选用Van Leer 分裂格式离散无粘通量项,并构造了一种限制器(Limiter)函数以保证TVD性质.数值算例表明本文方法是求解带操纵面偏转的翼身组合体绕流的有效方法.     

毛细管束中层流粘性振荡流强化传热的若干问题

本文是作者从事振荡管流换热研究的总结。研究结果表明:流体在管内作轴向振荡运动时,可以极大地强化轴向的传热。本文给出了计及管壁效应的振荡管流换热的精确解析解,公布了变Pr数的实验结果。在此基础上,讨论了影响振荡管流轴向换热的6个无因次参数的物理含义,并指出其中的某些无因次参数在物理本质上就是傅里叶数的倒数的平方根,这将有利于理解振荡管流换热这种非定常换热过程。由本文的精确解析解,经过简化,可以得到国内外已公布的近似解。实验数据和解析解预言的基本规律符合很好。     

军航飞机流穿越民航航线冲突探测与解脱问题

针对军航飞机穿越民航航线的飞行冲突问题,对飞机流汇聚飞行场景进行建模分析,提出了基于滑动窗口的汇聚飞行冲突探测方法和基于合作博弈的多机冲突解脱方法。当飞机流进入预先划设的汇聚控制区,通过滑动窗口判断向前看时间内飞行冲突。基于此,潜在冲突机间组成一个联盟,以解脱边界条件为安全约束,以合作博弈理论中联盟福利最优解均衡各机效益,在保证安全的前提下实现效益均衡。根据最优机动方向特点,利用免疫粒子群优化算法快速求解策略。仿真结果表明,该场景下提出的方法能有效解算出满足要求的解脱策略,并均衡军民航飞机解脱效益。      

关于非结构网格生成中的背景函数

背景信息对于用阵面推进法生成非结构网格是一个关键的环节.文章基于众多文献中提到的背景信息提供方法,提出一种更为实用的背景信息提供技术.此方法不仅避免了不必要的机时,内存的消耗,又使网格具有光滑性,网格的生成具有灵活性,方便性的特点.基于生成的网格采用有限体积法研究无粘高超音速化学非平衡流的流场分布,生成的网格能够成功的运用于高超音速化学非平衡流的数值计算。     

TVD格式与高超音速气动加热数值模拟

为了在化学非平衡流动中获得准确的流场解以及表面热流分布,将总变差减小TVD(Total Variation Diminishing)格式中的熵修正函数,由各向同性分布改为各向异性分布,同时让熵修正函数中的参数与流场中的压力梯度分布相关.将改进后的熵修正函数运用到高超音速化学非平衡绕流流场的数值模拟中,获得了较使用原有熵修正函数更为准确的流场参数和表面热流分布.采用改进的熵修正函数,可以提高壁面附近的粘性分辨率,降低热流计算结果对壁面附近法向网格尺度的敏感性.      

一种新型双环控制并联均流稳压电路

航天器多路输出DC DC变换器非主控输出电路常采用线性稳压器进行二次稳压调整，以确保非主控输出电压的稳定度和调整率满足负载使用要求.当变换器非主控输出电路电流大于线性稳压器额定工作电流时，常规下垂法并联均流扩容方式往往导致输出电压稳定度和调整率较差.为了解决该问题，本文在传统下垂法的基础上提出一种新型双环控制并联均流稳压电路，采用电压和电流双环控制方法，在保证了均流度的前提下，提升了变换器非主控输出电压的稳定度以及调整率.     

气液两相流中气泡速度的电磁互相关测量

气液两相流作为一种工业领域中普遍存在的流态,在测量时需考虑气相和液相间的速度差。相较于流动规律平稳的气相而言,研究不稳定的气泡流动状态具有更高的技术难度。本文针对空气-水为工质的两相流中一项重要参数——气泡速度进行测量,引入基于互相关运算的电磁检测方法,通过提高系统激励频率得到较大的检测信号。实验设计了高频(1 MHz)电磁检测系统,选用垂直上升管道,并在管道的2个平行截面分别安装电磁传感器,每组传感器均包含一个激励线圈及一个接收线圈,对2组接收线圈上的相位信号进行互相关运算,从而求得气泡的速度。实验中,对3种不同速度气泡采集到的信号进行了对比分析,相对误差控制在10%以内,提供了一种完全非接触非侵入的测量气泡速度的方法。此方法可以进行后续补充完善,用于其他工业场合,如金属液中气泡参数的测量等。     

气体跨流非圆管羣时散热与压差的关系

本文目的,是研究跨流非圆管羣时的散热及压差,从而找出非圆管羣散热与压力损失之间的关系,作为设计参考。 研究方法是由实验方法测量圆形、椭圆形及流线型三种管形在义排及顺排中有空气以不同流速跨流时的散热及压差。 由此得出结论:在现用的一种管距情况下,流线型及椭圆形管羣在同样散热时流阻可较圆形管羣小得多。     

球锥表面压力分布和静态气动力特性的工程计算方法

本文沿用内伏牛顿流基本思想,改进爆炸波理论在钝头附近的逼近精度,利用球锥无粘流数值解的零攻角表面压力数据(M_∞=4～8,θ_c=2.5°～20°)和轴对称比拟方法,拟合了一套既适用于短钝锥,又适合于细长锥的气动力工程方法。计算表明,本文提出的拟合公式计算精度优于Ericsson的内伏牛顿流理论和Blick的分段公式的精度,球锥几倍头部半径范围之内零攻角和有攻角表面压力分布计算误差分别为5%和10%左右,在稍远的地方逼近内伏牛顿流理论;气动力系数(C_A、C_N、C_M和X_(CP)/L)的计算值与数值解和风洞测力试验的测量结果吻合良好。