三维瞬态导热温度场边界元计算及一阶数学奇异的处理

推导了三维瞬态导热温度场边界元解法的积分方程和离散代数方程，并编制了相应的计算机程序，采用坐标变换的解析消去法，处理边界元的一阶数学奇异点，从而提高了算法的数值精度。用本程序计算了正方体和圆柱环瞬态导热温度场，数值解与解析解的吻合程度令人满意。     

发动机燃烧室主动冷却管道的热-力耦合分析

为了准确高效地进行发动机冷却管道结构的热-力耦合数值分析,基于边界单元法提出一种冷却介质中结构对流换热过程的新型计算方法,计算过程中采用径向积分法将对流换热边界积分方程中的域积分转换为等效的边界积分,从而显著降低计算难度。采用改进后的边界单元法和有限元方法分别进行发动机燃烧室主动冷却管道处的热-力耦合分析,计算并获得了该处的温度场、位移场和应力场,发现了不同物理量随管道轴向的变化规律。通过比较基于两种不同数值方法的计算结果可以发现本文数值方法在显著降低计算模型复杂度的同时,取得了合理的计算结果。因此利用本文方法可以简便有效地进行发动机燃烧室主动冷却管道结构的热-力耦合分析。      

MTPS蜂窝夹芯结构传热性能及热应力分析

对于金属防热结构的蜂窝夹芯结构的稳态情况,基于材料的全灰体假设,同时考虑热传导和热辐射两种传热形式对温度场的耦合作用,利用热流量守恒建立了蜂窝芯层温度场的非线性积分方程。离散化后利用数值方法得到方程组的数值解。对于美国兰利研究中心的实验结果,与本文方法的对比计算结果基本吻合。进一步,利用计算结果讨论了给定面板温度边界情况下,下面板、柱体层的灰度、蜂窝结构长径比对夹芯温度场的影响。并根据温度场和近似的应力分析模型,用半解析结果讨论了稳态情况下蜂窝芯层上的热应力。     

直升机非线性运动方程及其数值分析

基于牛顿法建立了单旋翼带尾桨直升机的非线性运动方程，提出了一种求解直升机非线性运动方程的数值方法。建立方程时，用欧拉角描述直升机在空中的姿态，旋翼采用有挥舞铰且在铰上带有弹性约束的模型，桨盘处的入流采用线性模型，通过对一算例直升机的操纵响应进行数值仿真，将非线性与小扰动线性化的直升机运动方程进行了对比分析。结果表明：从操纵输入到随后的2s左右，非线性和小扰动线性化的直升机运动方程的计算结果非常接近，但在此之后，非线性运动方程的计算结果更能反映直升机的运动规律。     

跨音速流中的一种边界元技术

边界元的名称源于C.A.Brebbian,边界元技术的特点是节点未知量均位于物面边界上,从而大大简化数值求解过程和降低对计算机容量的要求。本文将边界元技术推广应用于三维机翼的跨音速绕流的非线性问题,建立了一种适用于跨音速范围的带权余数表示式,并将其直接应用于全速位方程。为了保证冲波的不可逆特性,并使解在超临界流区域容易收敛,我们在全速位方程中增加有人工粘性项,而且将流场区域划分成许多小子域,并对每个子域应用格林定理。如果对每个子域我们所选用的插值函数对速位φ具有C~1连续性,则相邻元素交界面的面积分将相互抵消,使剩余的元素面积分(因而节点未知量)均位于机翼边界面上,从而导得边界元的积分方程。     

涡轮盘腔内流动和传热的数值模拟

将实际涡轮盘腔简化为转子系统, 并用混合长度模型、N-S方程和能量方程模拟了该系统的流场和温度场, 采用SIMPLE法数值计算了自由转盘周围的温度场;封闭转子-静子腔内的流场和温度场;有入流和出流的转子-静子腔内的流场。计算结果与实验值符合较好, 因而为数值模拟涡轮盘腔内的流场和换热打下了基础。      

数值解法在燃烧室出口燃气温度计算中的应用

为了求解燃气温度,利用双变量迭代法和赫夫法求解根据物质守恒、化学平衡、压力和能量守恒方程建立的非线性方程组,通过算例验证了2种数值解法求解化学平衡成分和温度的计算结果基本一致,在燃烧室出口温度场试验中,2种数值解法都能很好地满足测试需求,但双变量迭代法的计算收敛更快。试验结果表明:燃气分析的温度计算结果与热电偶的测温结果沿着燃烧室周向变化的趋势基本一致,燃气分析计算的温度基本大于热电偶测试结果。     

复合喷管热结构耦合计算的一种策略

基于流固耦合完成复合喷管流动换热、传热过程温度场的数值模拟,编写了计算流体动力学(CFD)与有限单元法(FEM)的接口程序,快速实现温度载荷和机械载荷数据的传递,提出了复合喷管热结构耦合计算一种策略.利用直接约束法考虑了复合喷管部件界面问高温碳化带来的边界非线性因素.计算结果表明耦合方法得到的复合喷管温度场是合理的,高...     

卫星姿控发动机混合物羽流场分区耦合计算研究

研究求解喷管内流场N-S方程数值计算方法,发展基于N-S方程物面边界滑移流理论计算技术。提出求解羽流核心区轴对称DSMC模拟方法与远场三维DSMC仿真方案,发展多组元混合物羽流DSMC仿真方法。研究求解卫星姿控发动机内外近场、远场、倒流区和物面相互作用影响区多流域流场分区耦合计算技术,建立了一套用于求解混合物燃气羽流及对太阳电池帆板与卫星体表面撞击污染影响数值模拟方法。通过对分别安装于某在轨卫星不同位置两个典型姿控发动机燃气五组元混合物羽流计算研究及相关结果对比分析,证实本文数值方法可靠性。     

强迫耗散非线性发展方程与完全平方守恒格式

从描述大气和海洋运动的强迫耗散非线性发展方程出发，对强迫耗散非线性大气和海洋方程组显式差分格式的计算稳定性进行了分析，构造了一类强迫耗散非线性发展方程的显式准安全平方守恒差分格式。理论分析和数值试验证明，这类显式准完全平方守恒差分格式是计算稳定的，值得推广应用。     

气冷叶片表面局部对流换热系数的数值计算

本文提出了联立求解二维可压缩边界层方程组和稳定导热方程的迭代法,从而修正了文献[1]假定物体表面温度等于来流总温的局限;同时提供了用这种方法计算气冷叶片表面局部对流换热系数的完整的FORTRAN程序。文中用近似积分法求解边界层方程组,用有限元素法计算叶片的稳定温度场。算例的结果与有关的理论及实验结果相接近。     

涡扇发动机高压涡轮前温度组合估计方法

考虑到目前暂无法实现机载条件下高压涡轮前温度直接、可靠的测量,提出一种用于涡扇发动机高压涡轮前温度估计的方法.基于涡扇发动机的能量守恒原理,建立高压涡轮前温度与气路参数的热力学关系,进而推导出高压涡轮前温度的6个估计模型.将各温度模型中不易测量的参数以整体的形式作为温度模型系数,并利用某涡扇发动机性能仿真模型建立温度模型系数与可测状态参数的多项式关系,最终确立高压涡轮前温度的组合估计模型.验证结果表明:组合估计方法在发动机健康及性能衰退状态下都具有较高的精度,其性能最好模型的方均根误差不超过1%.与已有线性拟合、神经网络等方法的对比也表明组合估计方法不论在精度还是性能稳定性方面都具有明显优势.      

LES of film cooling for different jet fluids

The present paper investigates the impact of the velocity and density ratio on the turbulent mixing process in gas turbine blade film cooling.A cooling fluid is injected from an inclined pipe at α=30° into a turbulent boundary layer profile at a freestream Reynolds number of Re∞=400000.This jet-in-a-crossflow(JICF) problem is investigated using large-eddy simulations(LES).The governing equations comprise the Navier-Stokes equations plus additional transport equations for several species to simulate a non-reacting gas mixture.A variation of the density ratio is simulated by the heat-mass transfer analogy,i.e.,gases of different density are effused into an an air crossflow at a constant temperature.An efficient large-eddy simulation method for low subsonic flows based on an implicit dual time-stepping scheme combined with low Mach number preconditioning is applied.The numerical results and experimental velocity data measured using two-component particle-image velocimetry (PIV) are in excellent agreement.The results show the dynamics of the flow field in the vicinity of the jet hole,i.e.,the recirculation region and the inclination of the shear layers,to be mainly determined by the velocity ratio.However,evaluating the cooling efficiency downstream of the jet hole the mass flux ratio proves to be the dominant similarity parameter,i.e.,the density ratio between the fluids and the velocity ratio have to be considered.      

热风洞中涡轮叶片温度场红外热像测量方法

提出了在热风洞中利用红外热像来测量涡轮叶片表面温度场的方法.针对热风洞特有的干扰因素,即石英玻璃窗口的透射比随着叶片表面温度变化以及燃气中的二氧化碳和水蒸气等组分参与热辐射所带来的干扰,在取得红外热像后按照燃气工况对温度场测量结果进行综合修正.考虑到叶片表面曲率的变化,通过几何上的变换重现了实际叶片表面上的温度场.结果表明:在热风洞的叶片温度场的红外热像测量中存在着110~140K的修正量.高温燃气环境中的红外热像测量结果必须按工况进行修正.      

压气机和涡轮转子三维温度场计算研究

为了对压气机盘和涡轮盘低循环疲劳寿命分析提供温度场数据,采用ANSYS计算软件的热分析模块对某型航空发动机高压压气机转子瞬态温度场和低压涡轮转子三维稳态温度场进行计算研究。重点分析旋转盘腔、旋转轴系、封严篦齿、榫头装配间隙等部位的换热规律;计算结果与相关文献进行对比,验证转子温度场计算方法的可行性。计算结果表明：从慢车到最大状态过程中,高压压气机盘最大径向温差先增加后减小直至稳定,中心孔附近较厚区域温度梯度最大;低压涡轮盘中下部沿轴向盘面温差很小,整个轮盘高温区域集中在轮盘盘缘。     

舱外航天服手套加热方案探讨

 针对舱外作业中航天员常常因手冷而不得不提前回舱的情况,提出了主动加热的方法。利用边界元法以边界点上未知函数满足积分方程为基础进行近似计算的特点,建立了手部和舱外航天服手套温度场三维计算模型,通过模拟试验件验证了模型。模型针对航天员处于太阳和地球辐射阴面时提出了舱外航天服手套主动加热的优化方案,并通过模拟试验件验证加热方法合理。     

Hoff型夹层板振动分析的边界积分方程法

 本文利用静力问题的基本解建立了Hoff型夹层板振动分析的边界积分方程。在数值实现过程中,边界上和区域内的未知函数均釆用线性插值以提高计算精度。数值结果表明本文方法具有未知量少和计算精度高等优点。对于低阶频率,本文结果与精确解的误差不大于1％。     

基于有限元/边界元耦合方法的管道进口声传播及声辐射模型

发展了一种基于有限元和边界元耦合方法的管道进口声传播及声辐射计算模型.该模型将整个声场分为内部有限域和外部无界域,分别用有限元和边界元方法求解控制方程,在两者之间的界面上使用具有物理意义的声阻抗参数进行匹配,并通过一种快速迭代方法实现全声场求解.这种迭代方法可以保证有限元刚度矩阵等带宽以及对称的特性不被破坏,有助于提高计算效率.该模型先得到了Levine-Schwinger标准解的检验,进而在无流动情况下对于简化的航发短舱进口管道模型进行了噪声辐射现象的数值模拟,最后基于计算结果分析了声衬对远场声辐射的影响.      

船舶燃气轮机进气系统中过滤装置数值模型化处理方法研究

为了研究船舶燃气轮机进气系统中,过滤装置数值模型化处理方法的模拟精度及其适用性。针对一种双通道船舶燃气轮机进气系统,分别采用Fan边界和多孔介质模型方法对过滤装置（百叶窗、滤清器）的复杂结构进行处理,开展进气系统内流场特性的详细数值计算研究。结果表明：进气系统大部分总压损失源于燃烧空气滤清器产生的总压损失和竖井中的流动损失;与Fan边界相比,多孔介质模型方法能够获得更贴合实际的进气系统流场,对系统总压损失的预测更具可靠性;两种方法计算出的进气系统流动损失相差不超过10%;使用Fan边界需保证过滤装置不产生明显的整流作用,且压力损失实验数据完备,而选用多孔介质模型需要已知过滤装置阻力系数、孔隙率等物理参数。     

压气机轮盘孔边应力的边界元分析

 本文根据压气机轮盘的轴对称问题,取出轮盘包含孔的一个扇形部分,以有限元法计算的结果作为扇形部分的边界条件,采用边界元法计算轮盘孔边的应力。这与直接用有限元法计算结果对比表明,边界元法具有数据准备工作量小、方便,计算点选择灵活,计算精度高等优点,用于计算孔边的应力集中问题是有效的。文中还推导丁常单元的二阶差分逼近式,用以求解边界上的应力。