改进的电热除冰系统仿真

建立了电热除冰系统矩形微段的二维简化物理模型,采用焓法模型建立相变导热问题的微分控制方程,用控制容积法对该模型的控制方程及其边界、初始条件进行离散,使用块修正技术对原有的离散方程组的求解方法进行了改进,得到各个时刻的温度分布及冰的融解情况.计算结果与文献实验结果吻合,验证了改进的算法和程序的正确性.并计算分析了时间步长的选择、网格的划分对计算结果的影响.      

三维数值模拟再生冷却喷管的换热

为了解液体火箭发动机喷管再生冷却的换热特点，采用数值模拟的方法，对内喷管燃气、壁面和冷却液建立不同的三维控制方程，进行流动和传热的耦合计算。在计算中，假定喷管流动为冻结流动，考虑燃气向壁面的对流换热和辐射换热；采用二阶迎风格式离散控制方程，采用DO模型离散求解辐射换热方程，水蒸气的吸收系数根据Leckner公式计算。计算模型采用缩比热试车发动机，数值计算结果与实验结果吻合较好，较准确地模拟出了喷管的壁面热流密度，得到了喷管燃气和冷却液的流场和温度场，对高压再生冷却喷管的设计具有指导意义。     

旋转离心叶轮内湍流的非结构化网格数值计算

给出了同位非结构化网格数值求解旋转离心叶轮内湍流流场的SIMPLE算法。网格划分采用改进的阵面推进法。控制方程采用通用的对流扩散方程并利用有限控制容积法离散。湍流模型采用标准的k-ε模型。形成的代数方程组用带有预条件器的共轭梯度平方法（CGS）求解。最后利用所给出的算法对一实验用的旋转离心叶轮进行计算。计算结果和实验数据基本相符,说明该算法合理,编程正确,可作为流场预测的理论工具。     

脉冲式燃烧风洞起动特性数值研究

利用计算流体力学软件Fluent,采用隐式双时间步方法模拟非定常过程,控制方程的离散采用二阶迎风格式,湍流模型选择剪切应力输运(SST)k-ω模型,并求解组分扩散方程,对脉冲式燃烧风洞的起动过程进行了数值模拟,分析了喷管起动时流场结构的变化情况和自由射流区超声速流场的建立过程.计算结果为改进风洞起动特性提供了一定的参考依据.      

过冷水滴撞击结冰表面的数值模拟

为了计算过冷水滴撞击结冰表面,本文应用雷诺平均Navier-Stokes方程和k-ε两方程紊流模型计算结冰表面外的空气流场。用欧拉法建立了空气中过冷水滴的运动方程,并采用有限控制容积法对方程进行了数值计算。为了提高计算精度,方程中的对流—扩散项采用MUSCL格式进行了离散。计算结果表明过冷水滴的大小和机翼攻角的变化对结冰区的大小和水收集系数的影响很大。      

二元矢量喷管流场的N－S方程计算

对二元矢量喷管在不同工况下的流场进行了数值模拟。求解的控制方程为三维非定常Ｎ－Ｓ方程。空间方向用有限体积离散，时间方向用Ｒｕｎｇｅ－Ｋｕｔｔａ法推进。粘性项的计算上采用了一种有限体积框架下新的离散化方法。湍流模型采用了广泛使用的Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数模型。计算表明，本文方法能准确地计算二元矢量喷管、非矢量喷管在设计与大设计落压比下的流场。     

CFD技术在动态压力校准中的应用研究

采用计算流体力学(CFD)流场可视化的手段,用数值模拟的方法研究了动态压力校准装置.研究工作主要集中于激波管动态压力校准装置和正弦动态压力校准装置的数值模拟,控制方程采用雷诺平均的N-S方程描述,方程的封闭采用精度相对较高的高雷诺数湍流模型,用有限体积法对控制方程进行空间离散和时间离散.并通过实验验证了控制方程和计算方法的正确性.     

基于多重网格方法的跨声速颤振数值模拟研究

使用多重网格的方法对二维翼型和三维机翼的跨声速颤振进行了数值求解.流场控制方程为N-S方程,湍流模型采用SA模型网格.计算网格采用结构网格,空间离散采用中心格式,使用双时间法进行时间推进.流场与结构之间的数据通过径向基函数(RBF)插值方法进行交换.通过耦合求解流场控制方程和结构运动方程得到Isogai(case A)机翼模型和Agard445.6机翼的跨声速颤振边界.通过与使用单重网格时得到的结果相比较得出,多重网格方法能够节省计算时间,提高颤振分析计算效率.     

电弧喷射推力器流动区域电磁场数值解法

为了有效地进行电弧喷射推力器流动区域全场(流场、电场、磁场和化学反应)耦合数值模拟并揭示其电磁场特征，建立了适当的电磁场模型，对电磁场控制方程的数值解法进行了详细研究。模型基于推力器稳态工作特征，控制方程由麦克斯韦方程组简化得到，采用有限控制容积积分方法离散椭圆型控制方程，采用9种不同的迭代方法求解离散方程。给出了不同迭代方法的收敛速度、数值稳定性和最终精度。研究表明，Gauss—Seidel逐线超松弛迭代法是推力器流动区域电磁场离散方程的一种快速有效的数值解法。     

超燃发动机燃烧室三维化学反应流场的数值模拟

采用ＬＵ隐式方法强耦合求解有限体积法离散的完全N-S方程及组份方程,数值模拟超燃发动机燃烧室内的化学反应流场。计算中采用了AUSM通量分裂格式,两方程及八方程的化学反应模型。计算得到了流场的结构与组份分布,并对两种化学反应模型的结果进行了比较与分析。     

Structural Nonlinear Flutter Characteristics Analysis for an Actuator-fin System with Dynamic Stiffness

The flutter characteristics of an actuator-fin system are investigated with structural nonlinearity and dynamic stiffness of the electric motor. The component mode substitution method is used to establish the nonlinear governing equations in time domain and frequency domain based on the fundamental dynamic equations of the electric motor and decelerator. The existing describing function method and a proposed iterative method are used to obtain the flutter characteristics containing preload freeplay nonlinea...     

粘弹壁板颤振的非线性动力特性

研究粘弹材料壁板在超声速气流作用下颤振时的分岔及混沌等复杂动力学特性。采用von Karman大变形理论及Kelvin粘弹阻尼模型建立壁板的动力学方程,通过线性活塞理论建立气动力模型。利用迦辽金法将壁板颤振模型转化为常微分方程组,并使用Gear的BDF方法进行数值求解。通过数值模拟研究了该系统在粘弹阻尼作用下的动力学行为以及粘弹阻尼的影响。计算结果表明,粘弹壁板颤振系统表现出丰富的动力学行为,其二次分岔特性很复杂。随着粘弹性阻尼的增大,系统的稳定解区域在减小,而静态屈曲解几乎不受影响,同时发现混沌运动区域也随着粘弹阻尼的增大而减小。     

等温不可压流体系统稳态数值计算方法研究

航空工业领域广泛存在着流体系统,任何流体系统均可化为只有节点和分支组成的通用流体网络,本文研究了这种通用流体网络的数值计算方法.将连续方程施加到节点,将动量方程施加到分支,建立起等温不可压流体系统的非线性控制方程组.根据Newton-Raphson迭代思想,研究了基于高斯-约当消去法的流体系统控制方程组数值求解方法.在此基础上,采用压缩矩阵存储技术开发了流体系统仿真平台FOCUSS-FS.通过算例进行验证,结果表明:本文计算方法能够准确获得任意流体网络的压强和流量分布,收敛速度快且效率较高,内存消耗小.     

气流激励下叶片动力响应分析方法

利用在物理过程上的耦合代替数学方程上的耦合,完善了双向顺序耦合求解的理论,进一步给出了双向顺序耦合的数值求解方法,并描述了其在物理上对应的过渡过程.结合具体工程算例,分析得到某小型发动机二级静叶在气流激励下的流场分布形式和叶片动力响应的特点,指出了双向顺序耦合数值求解方法的优越性.      

多点敷设支撑层的粘弹性约束阻尼梁振动分析

针对夹层约束阻尼梁结构的减振特性进行深入研究,在其基础上扩展多处敷设增设支撑层的约束阻尼,建立了多点铺设增设支撑层约束阻尼梁的振动特性模型,利用粘弹性材料的Maxwell模型、假设模态法、有限元理论和力学原理,结合功能特性方程和Lagrange方程,导出了增设支撑层约束阻尼梁的振动位移响应方程.通过对敷设双支撑约束阻尼简支梁在支撑层厚度变化的参数研究,可以看到增设适当厚度的支撑层能够达到更好的减振效果.     

导弹垂直发射系统三维超声速射流的数值模拟

对导弹垂直发射系统(VLS,Vertical launching system)内的三维湍流流场进行了数值模拟.采用了单元中心有限体积法和三阶迎风TVD差分格式离散Navier-Stokes方程,采用双时间步和预处理过程提高计算收敛速度,引入对角化的双时间ADI算法求解控制方程.仿真结果准确地模拟了发射系统内的激波形成过程,真实地反映出VLS中的复杂流动.研究表明,压力室底部射流中心的波后压强以大约1 200Hz的频率振荡.      

求解可压缩流动的低耗散E-CUSP迎风格式性能分析

采用对比分析的方法,对求解可压缩流动的低耗散E-CUSP(LDE)迎风格式的性能进行测评,目的是进一步更好地将其用于分离涡数值模拟。对比格式采用广泛使用的Roe格式。控制方程为Euler和雷诺平均Navier-Stokes方程耦合Spalart-Allmaras一方程湍流模型,采用全隐式无分裂Gauss-Seidel线松驰时间推进方法。结果表明,LDE格式具有强的捕捉激波能力,计算结果与Roe格式一致,和实验数据吻合良好。相对于Roe格式,LDE格式具有较快的收敛速度。     

单列叶栅紊流流场有限分析数值研究

本文将有限分析方法用于曲线座标系上紊流 N-S方程数值计算,研究了高雷诺数时叶栅粘性紊流流场。有限分析方法在网格单元上对非线性偏微分程进行线化处理,在解析边界条件下求出线化偏微分方程的解析解,以此解构造离散方程。有限分析方法能够根据对流的方向和大小自动改变格式系数,并具有数值扩散小、精度高和稳定性好等优点。本文以 k-ε紊流模型模化紊流,以壁面函数方法处理近壁区流动参数。      

高精度的机翼静气动弹性数值模拟方法研究

气动结构耦合是研究静气动弹性问题的关键所在。采用三维Euler方程为控制方程,以有限元方法计算结构弹性变形,用细平板样条法(TPS)实现气动结构数据间的转换,进行机翼静气动弹性数值的模拟。算例结果表明,本方法是一种具有很高计算精度和工程实用性的气动弹性模拟方法。