示踪气体法求客舱通风量和空气龄的方法

 传统的用浓度变化求通风量的公式是建立在单室通风模型上的。通过建立示踪气体质量平衡方程式可得到用单种示踪气体和多种示踪气体浓度变化求客舱（多室）通风量和渗风量的方法。为了提高客舱内空气品质,引入了空气龄的概念并给出了由客舱内典型点的示踪气体的浓度求空气龄及通风效率的方法。     

亚/超声速楔状流层流边界层速度与温度相似解及拟合解

利用相似变换获得了楔状流层流边界层无量纲流函数的3阶非线性常微分方程,用Runge-Kutta法求解微分方程获得了不同楔形角楔状流层流边界层无量纲速度随相似变量的变化曲线;推导了亚声速和超声速楔状流层流边界层无量纲温度关于相似变量的2阶线性齐次和非齐次微分方程,获得了温度分布的通解,恒壁温条件下亚声速楔状流和绝热壁面条件下超声速楔状流层流边界层无量纲温度解析解及指数函数形式的拟合解.以楔形角为0为例利用相似变换研究了超声速条件下气体压缩性及黏度随温度变化等因素对层流边界层速度与温度的影响,得出不可压缩常物性与可压缩变物性条件下无量纲速度相对误差绝对值小于9.8%的结论.研究表明:Pr越大贴近壁面处无量纲温度变化越剧烈;超声速条件下壁温低于绝热壁温时黏性耗散作用可以使层流边界层气体温度从壁面到主流间出现先升高后降低的变化.      

发动机故障诊断的MONTE CARLO通解算法

给出了发动机故障诊断的 Monte Carlo通解算法,该算法可以有效地解决发动机故障诊断算法中由于故障方程存在多重共线性所引起的误诊、漏诊与多解问题。首先给出求解故障方程的 Monte Carlo算法,该算法能够保证得到满足故障方程的全部合理解;能够简单地用于各种故障相关性准则 (例如各种残差向量范数准则 )以及亚定故障方程的求解问题;并且算法简单易行,无需为每个特定情况专门编制计算程序。其次提出了基本解与通解的概念,并且给出了利用主成分分析与利用主因子模型求基本解的两种有效算法。利用基本解与通解算法可以将由于多重共线性引起的复杂故障诊断结果表示为简单明了的形式,有利于对诊断结果作出正确决策。用计算机模拟方法对算法的有效性进行了分析,对于 JT9D发动机气路方面的 2 4个实际故障样本,所给出的算法的确诊率为 86%～ 92 %     

Simulation of low-pressure water flow through the Laval nozzle taking into account evaporation on the basis of a two-temperature model of two-phase flow

A water flow through the Laval nozzle is studied using the model of the one-velocity two-temperature two-phase medium motion written on the basis of the Navier-Stokes equations. As an equation of state for the carrying phase, Tait??s equation is used. The vapor phase dynamics is described by the conservation equation for the concentration of bubbles having the critical radius, the Rayleigh-Lamb equation, and the equation of vapor pressure inside a bubble. The equations for the one-velocity two-temperature medium dynamics are solved by the Beam-Warming method. The explicit contraflow scheme is used to solve the Rayleigh-Lamb equation and the conservation equation for the bubble concentration. The solutions of the quasistationary type obtained for the plane Laval nozzle make it possible to estimate the local flow parameters and the rate of evaporation at with different values of counterpressure. We compare the solutions for the models of independent phase motion and those taking into account the interphase exchange of mass, impulse, and energy.     

燃气排导箱三维场的数值计算

对燃气排导箱内的典型三维稳态流动过程进行了粘性不可压数值计算,采用较为成熟的ＳＩＭＰＬＥ方法。所发展的计算程序可求解质量、动量和能量守恒方程,ｋ－ε双方程湍流模型方程,以及热辐射通量方程,并进行了程序的算例验证。三维燃气排导箱内的计算结果也得到了一些实验测量的数据的验证,得出一些有工程实用价值的结果和结论。     

A Comparison of Numerical Methods for Solving Optimal Control Problems

Linear systems with quadratic performance criteria typically lead to two-point boundary value problems. The numerical solution of such problems is often difficult because of roundoff and truncation errors. In this paper numerical methods for digital computer solutions are considered for a system described by a given second-order linear differential equation and the results are compared. The results show that where accurate analytic solutions cannot be obtained, the invariant imbedding technique yields more accurate results, for a given integration grid interval, than the other methods considered.     

考虑低温增压真实气体效应的运输机气动特性数值模拟研究

为了给低温风洞试验数据修正提供参考,本文利用数值模拟手段研究了低温真实气体效应相比较于完全气体对飞行器气动特性的影响,以及该影响与雷诺数影响相比所占比例的大小等问题。文章应用 Aungier-Redlich-Kwongz 方程,发展了适用于模拟氮气低温高压真实气体效应的 RANS 求解软件。与 NIST 数据的对比表明,该状态方程在5倍大气压下,比热比等参数误差在0.3%以下。同时,标模测试结果表明本文软件计算精度与国外软件相当。应用本文方法研究低温增压风洞中氮气真实气体效应：以典型运输机构型 DLRF6为模型,分别计算了高速、低速状态下各种不同温度和压力工况下的流动。计算结果表明,在低温增压情况下,真实气体效应引起的气动力差异很小,升力、阻力、力矩最大相对误差均在0.3%以下,与雷诺数效应引起的偏差相比可以忽略不计。因此,以氮气为介质的低温风洞试验研究可以采用完全气体假设。     

基于二维厚板模型的波箔片轴承静特性

建立平箔片的二维厚板有限元模型,运用有限单元法和有限差分法耦合求解Reynolds方程和气膜厚度方程,研究了在两个工作转速下气体波箔片轴承在中截面和边缘处最小气膜厚度随轴承承载力变化规律.通过数值仿真对该模型、一维梁模型、二维薄壳模型和文献实验数据进行对比分析,结果表明:在轴承中截面处,3个模型的最小气膜厚度仿真结果都与实验结果符合得很好,但在轴承边缘处,由于二维厚板模型考虑了平箔片的剪切效应,因此其最小气膜厚度比二维薄壳模型的结果更接近实验值,而一维梁模型只考虑轴承圆周方向,因此不能体现气膜厚度沿轴承长度方向的变化规律.通过研究,为分析箔片轴承动力学特性奠定了理论基础.      

悬臂式挠性薄片气体动压径向轴承分析研究

本文建立了悬臂式挠性薄片气体动压径向轴承的特性分析计算方法,用差分方法耦合求解径向轴承的雷诺方程和薄片弹性变形方程。给出了径向轴承刚度、阻尼的计算方法和稳定性计算方法。计算分析了设计的悬臂式薄片气体动压径向轴承的静特性和动特性。计算结果表明,对于所有的偏心率范围和可能的工作运行范围内是恒定稳定的。      

有限长度的过渡领域微槽道流动(-分析预测和试验的比较)

在连续介质和滑流区有限长度槽道内的稳定气流,通过整体质量守恒规律得到了沿槽道轴向的压力分布和质量流率的分析表达式.在过渡领域区,也提供了槽道流动压力分布的严格的动理学理论解.在Poiseuille流动的质量流率动理学理论解给定的情况下,压力分布的控制方程被证实是被前一作者从轴承润滑问题移植到解微槽道问题的退化Reynolds方程.同时得到了质量流率分析表达式,它们和现有试验数据的比较相符很好.     

求解输流管道动态响应问题的一种有效方法

研究了两端铰支输流管道在脉动内流作用下的参数共振问题。根据数值分析中发现的参数共振响应特性,给出了第一和第二振型运动相互解耦的控制方程,并用谐波平衡法导出了参数共振响应曲线的近似解析解。算例表明,第一振型近似解与数值解吻合得很好,而第二振型解在共振点附近与数值解略有差距,但二者的整体走向完全相同。     

亚、超音速小侧滑角机翼升力面理论和数值计算方法

本文推导出侧缘与自由流不平行的常值面涡半无限角翼基本解,其特征是从面涡小块背风侧缘拖出的涡不是顺着自由流,而是沿着侧缘拖出,从而克服了用Woodward方法计算侧滑翼绕流中使用对称分划小块的困难,提出了适用于从亚音速直到超音速范围的推广的Woodward方法。本方法将问题分解为若干个对称与反对称问题,大大地节省了计算机内存量或计算时间。本文还从基本方程,或者直接从数值计算方法出发,推导出准相似律。计算结果与精确解以及和其他方法的结果符合得很好。     

热空腔-喷气流的组合辐射

讨论了热空腔(喷气发动机的加力筒)和喷气流的组合辐射问题, 推演了组合辐射的辐射传输微分方程, 导出了有边界辐射条件下的解。讨论了由多个温度不同的壁面组成且腔内充满高温燃气的热空腔辐射传输方程解的具体形式, 导出了它向后半球的光谱辐射强度的表达式;同时还推导了热空腔-喷气流组合辐射的光谱辐射强度的表达式。最后给出算例, 算例结果表明, 热空腔-喷气流组合辐射的空间分布具有与实测分布较一致的梨形包络, 其光谱分布也与实测的光谱分布的一致性较好。      

超音速、小高超音速三元薄翼非定常二次理论的准确度及其适用范围

 本文处理了超音速三元薄翼非定常问题,通过PLK法使二次解均匀有效。首先考虑零攻角或初始攻角时,已知基本定常绕流叠加高-量级的非定常小扰动流,把它线性化。本方法从健全的基本方程出发,使用高马赫数近似,将非定常二次方程化简,其形式与定常二次方程类似,因而有可能利用定常二次理论的方法求解。特解是求解的关键。鉴于精确特解的复杂性,本报告采用了一种近似特解。 本方法适于一般超音速和完全高超音速之间的马赫数区域（约3～8）,折合频率可达至1左右。可较精确地估计厚度,初始攻角对非定常气动力,力矩的影响。 目前据我们所知,还没有有关实验数据,只能和一些理论结果进行比较。为此对低频有初始攻角的超音速前缘平板三角翼进行了计算,在马赫数3～8,与D.D.Liu^[6]比较吻合。计算结果表明,三元薄翼二次理论可用到高超音速相似参数Mδ=1.0。     

一阶微分方程终值问题拟解的存在性

利用拟上下解方法,研究了一阶微分方程终值问题u′(t)=f()t,u(t),u(T)=xT拟解的存在性,并通过构造单调迭代序列得到该终值问题的最大最小拟解。     

气体静压止推轴承静动态性能及振动抑制

为探讨抑制气体轴承自激振动的新手段,研究了圆盘形单个小孔节流的气体静压止推轴承振动问题。通过两平行圆盘间的气体流动方程,推导建立了气体轴承静动态特性分析模型。通过数值分析的方法,研究并总结了供气压力和气膜厚度对气体轴承静态承载力、静态气体质量流量、静动态刚度以及动态阻尼的影响,计算得到了一系列气体轴承静动态特性变化曲线。然后针对系统出现负阻尼情况,通过引入非线性能量阱(NES)来抑制自激振动的问题。研究结果表明,提高供气压力可以有效地提高气体轴承的静态性能;气体轴承的自激振动主要源于气膜挤压的负阻尼特性;同时,当NES的阻尼超过临界阻尼时,可以使系统的振动得到有效的抑制。     

分区和有限容积法计算长尾喷管中的湍流流动

为了给长尾喷管热防护结构设计提供热边界条件，从轴对称二维雷诺平均N－S方程出发，采用分区算法和有限容积法计算了长尾喷管中的湍流流动，分析了流场的结构特性，提出了长尾喷管中燃气温度，压强和马赫数分布的某些规律以及燃气温度沿壁面的分布，计算和分析结果提出了长尾喷管热防护结构设计应注意的事项。     

Supersonic Two-Dimensional Minimum Length Nozzle Design at High Temperature. Application for Air

When the stagnation temperature of a perfect gas increases, the specific heat ratio does not remain constant any more, and start to vary with this temperature. The gas remains perfect, its state equation remains always valid, except it will name in more calorically im-perfect gas or gas at High Temperature. The goal of this work is to trace the profiles of the supersonic Minimum Length Nozzle with centered expansion when the stagnation temperature is taken into account, lower than the threshold of dissociation of the molecules and to have for each exit Mach number several nozzles shapes by changing the value of the temperature. The method of characteristics is used with a new form of the Prandtl Meyer function at high temperature. The resolution of the obtained equations is done by the second order of finite differences method by using the predictor corrector algorithm. A study on the error given by the perfect gas model com-pared to our model is presented. The comparison is made with a calorically perfect gas for goal to give a limit of application of this model. The application is for the air.     

考虑稀薄气体效应的止推箔片轴承静特性分析

以波箔型气体止推箔片轴承为研究对象,基于1阶滑移速度边界条件,建立了考虑稀薄气体系数的修正雷诺方程.结合Newton-Raphson迭代法和有限差分法,耦合求解Reynolds方程和润滑膜厚度方程,仿真获得了止推箔片轴承的轴向承载力、起飞转速等静态性能,并研究了稀薄气体效应和轴承结构参数对止推箔片轴承静态特性的影响.结果表明:稀薄气体效应将使止推箔片轴承的轴向承载力减小,起飞转速增大,且稀薄气体效应的影响随轴承工作转速、箔片变形柔度系数的提高有所减弱;止推箔片轴承的安装间隙小于10μm时,起飞转速随安装间隙的减小急剧增大,因而在止推箔片轴承装配时须严格控制轴向间隙.      

粒子群混合算法在变导向器面积计算中的应用

求解非线性方程组经典方法具有严格的局部收敛性;粒子群等进化算法解决了全局收敛的问题,但计算效率偏低,存在最优解不稳定的问题。结合经典Newton-Raphson法的超线性收敛速度和粒子群算法全局收敛能力的粒子群混合算法具备2类算法的优点。在迭代初期采用粒子群算法获得的近似全局解作为Newton-Raphson算法的初始值,以确定高精度的解。利用粒子群混合算法在发动机变导向器面积的大偏离计算中获得了较好的收敛效果,解决了常规Newton-Raphson法不收敛的问题。