复杂外形飞行器热流的NND有限元数值计算方法

有限元方法是一种加权积分的方法,其特点和权函数的特性使得有限元方法对于复杂外形的边界处理和计算诸如热流(与温度的一阶导数相关)等物理量时有其内在的优越性.本文采用张涵信的NND思想发展了复杂流场的NND有限元计算方法,算例表明,在相同的网格雷诺数下,采用积分的方法提高了物体表面热流的计算精度,对高超声速复杂外形飞行器的绕流获得了良好的数值模拟结果.     

高超声速复杂外形热流的NND有限元数值计算方法研究

有限元法是一种加权积分的方法,这个特点和权函数的特性使得有限元法在计算诸如热流(与温度的一阶导数相关)等物理量时有其内在的优越性.本文研究了高超声速复杂外形热流的NND有限元数值计算方法,并对高超声速球钝锥外形和带翼飞行器外形的绕流进行了数值模拟,算例结果令人满意.     

复杂外形飞行器热流的NND有限元数值计算方法研究与应用

与微分法相比,有限元方法作为一种加权积分的方法,其特点和权函数的特性使得有限元方法对于复杂外形的边界处理和计算诸如热流(与温度的一阶导数相关)等物理量时有其内在的优越性.本文采用张涵信的NND思想发展了复杂流场的NND有限元计算方法,采用积分的方法提高了物体表面热流的计算精度,对高超声速复杂外形飞行器的绕流获得了良好的数值模拟结果.     

复杂组合体亚音速气动力计算改进的面元法

发展一种低阶面元法计算复杂组合体的亚音速气动力。本方法采用两种形式的速度势的内部第一类边界条件，通过对速度势的数值微分得到速度，进而计算气动力参数。本方法克服了早期低阶面元法对面元剖分敏感和高阶面元法的运算量大的缺点，兼有高阶面元法的精度高和早期低阶面元法的运算量小的优点。用本文给出的方法，计算了多个模型（包括型号）的气动特性，计算结果和实验结果有很好的一致性，表明本方法灵活、稳定、准确和有效。     

复杂外形网格生成技术

复杂外形网格生成是进行流场数值模拟的一个关键技术及前提。本文论述了采用保角变换方法,求解偏微分方程方法以及它们在代数方法结合生成二维及三维具有多块网格结构的贴体且与边界正交网格的技术。     

复杂全机带多个外挂物气动力数值计算

把镜像理论用于计算机翼与机身,机翼与外挂物的气动干扰,并开发了复杂全机带多个挂物的面元生成程序。     

战术导弹零升阻力雷诺数效应及修正方法

叙述了雷诺数对战术导弹零升阻力的影响,并给出了适合战术导弹的零升阻力系数雷诺数效应修正方法,即变雷诺数试验外推修正方法及工程计算方法。修正结果表明,修正方法是可行的,变雷诺数试验外推法得到的修正量比工程计算得到的修正量更为合理,对于外形简单的战术导弹工程计算仍有较好的精确度。     

航空发动机试车台附加阻力修正方法

附加阻力(进气冲量阻力和外部冲量阻力)是航空发动机试车台架推力的重要组成部分,准确地确定附加阻力对提高发动机试车台推力测量的准确性有十分重要的意义。本文对发动机总推力和试车台的对比标定进行了介绍,对发动机在露天试车台、室内地面台和高空台上附加阻力的确定方法和修正方法进行了详细地研究,给出了高空台与高空台、高空台与地面台的对比标定试验结果。      

基于阻力分解策略的翼型阻力计算方法研究

研究并应用了基于翼型阻力分解策略的中场积分计算阻力的方法。中场积分方法是近年来提出的一种新的计算阻力的方法,方法将流场分为三个不同的区域：波阻区、型阻区、数值耗散区,通过对不同区域进行积分将阻力分解为波阻、型阻、数值耗散阻力。计算结果表明中场积分方法与传统的表面积分方法和远场积分方法相比,不仅可以将阻力分解成不同部分以便于进一步分析阻力产生机理,而且在计算中对网格量的要求较低,即使在网格量较小的情况下也能得到较精确的结果。     

增升状态下的全机位流数值模拟

研究了势流理论下的全机增升装置复杂外形气动力数值计算方法,该方法以低阶面元法为基础,具有计算量小、速度快的优点.复杂流动的粘性效应主要通过附面层粘性修正技术加以体现.采用该方法对某增升装置翼身组合体外形压力分布特性进行计算,并和风洞实验数据相比较,验证了方法的正确性和实用性.在此基础上对某型飞机真实复杂外形的压力分布特性进行了计算,为复杂外形飞行器增升装置的初步设计、方案选型与评估提供了非常理想的工程技术手段.     

阻力精确分解的方法研究与数值模拟

根据阻力产生的物理机理将阻力分解为近场阻力和中场阻力,建立实现了阻力精确分解方法:近场阻力法和中场阻力法,前者对飞行器表面压强和切应力积分将阻力分解为压差阻力和摩擦阻力,后者通过区域划分式积分将阻力分解为波阻、粘阻、诱导阻力和数值耗散阻力.选取RAE2822翼型、M6机翼和某宽体标模进行数值模拟,验证阻力分解方法的正确性,对比两种方法进行阻力分解辨识的能力与不足.研究结果表明,中场阻力法可以给出详细的阻力构成,更有利于进行阻力产生分析和减阻优化设计;近场阻力和中场阻力计算可以很好地达到阻力平衡,但由于消除了数值耗散阻力,中场阻力法的计算结果要比近场阻力法更接近实验值.     

翼型—扰流片的分离气动特性计算

本文介绍了用涡面元法模拟带扰流片的翼型低速无粘分离绕流。在翼型和扰流片的面元上分布线性变化涡。在扰流片后的上下分离流线的面元上分布等强度的涡。上分离流线始自扰流片的梢部，下分离流线自翼型的后缘引出。分离所泡由两离散涡结尾。气泡内总压为常值，它与涡强大小一同求解。分离气泡的形状在迭代求解过程中确定。压强分布和升力系数的计算值与现存文献的数值结果和实验数据是一致的。     

粘性不可压流动压力修正的有限元计算

将ＳＩＭＰＬＥ方法压力修正思想应用于有限元方法，对粘性不可压缩流动的ＳＵＰＧ有限元方程，建立将压力与速度分离，进行压力修正迭代的计算格式；对典型算例进行多个Ｒｅ数的层流，湍流计算，显示此方法具有良好的收敛性和精度。     

基于非定常面元/黏性涡粒子混合法的旋翼/平尾非定常气动干扰

旋翼/平尾非定常气动干扰是导致直升机纵向“抬头（Pitch-up）”现象的主要原因。为在直升机设计阶段准确分析旋翼/平尾非定常气动干扰以及由此引起的低速纵向操纵特性变化,通过涡量等效原则和Neumann物面边界条件建立了适用于旋翼/平尾气动干扰分析的非定常面元/黏性涡粒子混合法。该方法耦合了考虑尾迹时变效应的非定常面元法、黏性涡粒子法及涡量镜面法,以准确模拟旋翼和平尾的非定常气动载荷、旋翼尾迹的非定常特性以及旋翼尾迹对平尾的气动干扰效应。首先通过计算NASA ROBIN（Rotor Body Interaction）旋翼尾迹几何和诱导速度分布,并与实验测量值、时间精确自由尾迹及CFD计算结果对比验证方法的准确性。相比于时间精确自由尾迹,本文方法计算精度更高。随后分析了旋翼/平尾非定常气动干扰对平尾向下气动载荷和气动导数的影响,并分析了平尾构型对旋翼/平尾非定常气动干扰的影响规律。分析表明：旋翼尾迹与平尾干扰导致低速状态的平尾载荷突增,气动导数反号;低平尾气动载荷突增较大,高平尾较小,但高速气动导数反号;前置平尾载荷突增量减小,但对应速度范围较宽;右旋直升机右平尾载荷突增量较小,但气动导数特性基本不变。     

某运输机全机构型雷诺数效应阻力修正研究

采用CFD方法分析了雷诺数效应对某运输机全机构型阻力的影响。基于对某运输机的翼/身/架/舱组合体绕流进行的计算分析,归纳总结出了一种适应于复杂构型的低雷诺数到高雷诺数的阻力修正方法。通过验证可知,该修正方法修正效果良好,为工程中的阻力预测提供了一种简洁有效的方法。     

基于黏性伴随方法的跨声速机翼气动优化设计

 进行了基于黏性伴随方法和Navier Stokes方程的跨声速机翼气动优化设计研究。分别推导了适用于三维跨声速机翼气动反设计和减阻设计的黏性伴随方程、边界条件和梯度求解表达式,并研究了伴随方程的数值求解方法。通过将网格生成、流场计算、黏性伴随方程数值求解、梯度求解和拟牛顿优化算法等几方面的有效结合,发展了一种跨声速机翼气动优化设计方法。为了提高计算效率,将多重网格方法应用到方程的数值求解中来加速收敛。跨声速机翼反设计和减阻设计算例验证了本文所发展的方法的正确性。采用本文的方法进行优化设计,一般通过20~30次迭代就能得到满意的结果。     

1种发动机自适应模型的修正方法

为解决由于工艺落后和性能衰退所导致的发动机模型与实际发动机失配的问题,从而达到视情维修的目的,建立了1种基于遗传算法的自适应模型修正方法.利用遗传算法高效率的全局搜索能力,以发动机性能参数为目标,优化部件参数,实现了模型的修正,并通过分析发动机性能参数对部件参数的敏感性,研究了选取不同组合的待修正部件参数对模型修正精度的影响.结果表明:采用该修正方法所建立的发动机模型准确性高、鲁棒性强,且具有良好的通用性,为发动机的健康管理和故障诊断提供了重要依据,具有很好的军事和工程应用价值.