非结构网格中LU-SGS隐式算法的非平衡性影响

在非结构网格中用LU-SGS隐式算法求解欧拉方程时,两个近似分解因子的项数可能相等,为平衡;也可能不相等,为非平衡.采用对网格重新编号的方法很难达到我们希望的平衡.本文对非平衡性的影响进行了探讨.对二维问题,四边形单元的非结构网格,设计出平衡与非平衡的编号方式.先对标量模型方程分析LU-SGS隐式算法的增长因子,然后通过数值试验来验证这种非平衡性的影响.结果表明,尽管非平衡时也能达到收敛,但平衡却远优于非平衡的情况.     

Java实现Fibonacci数列算法

通过分析Fibonacci数列的特性，利用Java语言面向对象的编程思想实现了三种求解Fibonacci数的方法，即交互变换与累加法、递归解法和生成函数解法；还说明了递归问题转换为非递归问题的一般方法。     

遗传算法在航空发动机非线性数学模型中的应用

收敛性是航空发动机非线性数学模型的重要指标。现有发动机平衡方程迭代解法还不能保证模型大范围收敛。以某涡扇发动机为对象,采用遗传算法求解发动机非线性数学模型,将模型中的发动机平衡方程求解转换为极小值优化问题,建立了遗传算法计算模型,重点分析了采用遗传算法求解的适应度函数设计方法。数值仿真结果表明,与牛顿-拉夫逊解法相比,采用遗传算法方法可实现模型大范围收敛。     

非稳态流体网络方法在发动机空气冷却系统中的应用

将流体网络法与非稳态模拟方法相结合应用到发动机空气冷却系统沿程温度、压力、流量的非稳态特性预测中.本方法将空气系统简化为节点和元件组成的网络,节点和元件处建立守恒方程,针对传统解法的计算稳定性问题引入了压力修正和积分方法相结合的求解方法.程序应用到可参考对照的模型中得到验证.验证后的程序被应用到发动机轴部冷却通道空气系统非稳态特性预测中.      

自适应变搜索域遗传算法及其在发动机模型中的应用

针对发动机数学模型中非线性方程插值解法的不足,提出非线性方程自适应变搜索域遗传算法解法。论文详细分析了方程求解转化为遗传算法优化的数学描述、依据方程解临域特征的自适应变搜索域机理及算法的具体实现技术。实例计算与发动机仿真结果说明:所设计的变搜索域遗传算法可作为通用的非线性方程解法,相比较于简单遗传算法,能成倍提高计算效率;替代插值解法应用于发动机模型,可有效提高模型的收敛性。     

航天器再入大气层热力分析

以OREX（orbital reentry experiment vehicle）的飞行试验数据和相关的CFD数值模拟结果为基础,采用传热理论及相关公式,分析计算了OREX再入大气层过程中的轨迹,驻点处热流密度非平衡假设和平衡假设下的换热问题.所计算的轨迹、热流密度非平衡假设下计算的驻点温度和热流密度值同试验数据及相关的CFD数值结果取得了很好的一致,相应的计算方法可作为航天器驻点热力分析的通式.然而在平衡假设条件下,尽管计算所得到的驻点热流密度与之前的CFD数值结果差别不大,但这种情况下计算得到的温度与试验数据不符,这应归结于计算的热流密度结果对驻点处温度变化的不敏感.比较非平衡假设和平衡假设下的换热计算结果表明,对于航天器再入过程中的热力探讨不能仅仅满足于热流密度分析,对温度的考察或许更重要.      

基于谐波平衡法的含Iwan模型干摩擦振子非线性振动

通过谐波平衡法,研究了由Iwan模型和1个质量块构成的干摩擦振子系统的自由振动和受迫振动问题.对于自由振动,推导了系统卸载和加载过程中各自对应的非线性运动方程;通过谐波平衡法,求得了各方程的1阶谐波解;为了得到方程正确的数值解,提出了确定零速度时刻的策略;数值算例结果表明:谐波平衡法与数值方法得到结果吻合较好,从而验证了谐波平衡法求解此类问题的有效性.对于系统微滑移下的受迫振动,采用谐波平衡法求得了系统的幅频关系.最后,对Iwan模型阻尼特性的研究发现,模型等效黏性阻尼比与其振幅存在非线性关系.      

VIE-ODDM在电磁散射特性分析中的应用

在重叠型区域分解法的基础上,对体积分方程方法进行了改进,实现了基于区域分解法的体积分方程方法.通过几个算例,验证了这种算法占用内存小、迭代速度快,计算结果准确,具有分析电大尺寸目标电磁散射问题的能力.     

梁在短时受热受力下的弹塑性弯曲

本文考虑梁在短时受热受力情况下的弹塑性弯曲,应用平面假定,考虑非均匀温度引起的热应力和非均匀材料性能。对于不同温度的材料应力应变曲线,提出了用温度函数和应变函数的乘积表达应力的近似关系式。并用多项式和冪函数作为应变函数,分段近似材料曲线。近似表达式对RS-120的材料曲线给出了很好的结果。本文对弯矩较小、热应力不能忽略的弹塑性弯曲问题给出了同多项式近似解法, 结果准确, 计算相当简单。对于弯矩较大、热应力可以忽略的情况给出了同冪函数近似解法, 计算非常简单, 结果亦好。文中并从分析和计算结果看出: (1)热应力对弯矩曲率关系的影响较小; (2)对于热弹塑性弯曲梁变形问题较热应力问题更为重要。     

高超声速喷管三维热化学非平衡流场的数值模拟

利用数值求解三维热化学非平衡Navier-Stokes方程的方法和多块网格技术,数值模拟了典型高焓激波风洞锥型喷管的三维非平衡流场,分析了热化学非平衡效应对喷管流动的影响,给出了平动温度、组分N2和O2振动温度、马赫数和组分质量分数在轴对称线上的分布,也给出了喷管出口附近的速度和温度沿径向分布等结果,获得了喷管流场的详细信息,并与已经存在的采用轴对称Navier-Stokes方程计算的结果,进行了比较,二者吻合得很好.计算结果表明:目前的代码是可以模拟多块结构网格下的热化学非平衡流动的.     

用随机选取法（RCM）计算高温平衡气体的激波管流动

本文用RCM方法计算了高温平衡气体的激波管流动。在求解Riemann问题时,激波关系式采用精确解,而对于稀疏波关系式,本文建立了等效比热比的局部理想气体近似关系式,从而使RCM方法大为简化,并具有足够精度。     

Closed-form solution for the Kalman filter gains of time-varyingsystems

A method for derivation of closed-form solutions for the differential Riccati matrix equation for specific time-varying systems is presented. It allows more insight into the nature of the solution. It reduces the on-line computation requirements, since it does not require on-line solution of a differential equation. Sufficient conditions for the existence of the closed-form solution are given. The method is applied to a target tracking problem     

数值求解二维定常跨音速流场多重解现象探讨

本文用时间推进法求解欧拉方程组计算二维平面叶栅绕流问题数值实验是否有多重解现象。计算结果表明,求解欧拉方程组进行流场数值模拟最终收敛与初场无关,不存在解不唯一问题,进而计算初场熵值大小,也发现初场熵值大小不影响终场结果。     

Analysis of generalized guidance laws for homing missiles

It is shown that the optimal guidance problem can be handled as a Bolza-type problem which can be solved by straightforward calculation. The solution proves to be different from that given by Yang, et al. (1989) and is verified to be a minimum performance index solution. Also presented is an alternate unified approach to the solution of the nonlinear system of differential equation of the generalized proportional navigation law discussed earlier by Yang, et al. This unified equation has been solved for all the different cases of proportional navigation for which closed-form solutions have been published by earlier authors     

Gaseous plume flows in space propulsion

This paper presents a gaskinetic study on high-speed, highly rarefied jets expanding into a vacuum from a cluster of planar or annular exits. Based on the corresponding exact expressions for a planar or annular jet, it is convenient to derive the combined multiple jet flowfield solutions of density and velocity components. For the combined temperature and pressure solutions, extra attention is needed. Several direct simulation Monte Carlo simulation results are provided to validate these analytical solutions. The analytical and numerical solutions are essentially identical for these high Knudsen number jet flows.     

关于跨声速流动问题的多重解

本文讨论了跨声速流小扰动方程Riemann问题的多重解及物理上可实现的解的判别。     

A quasi-one-dimensional model for hypersonic reactive flow along the stagnation streamline

This study proposes a quasi-one-dimensional model to predict the chemical non-equilibrium flow along the stagnation streamline of hypersonic flow past a blunt body. The model solves reduced equations along the stagnation streamline and predicts nearly identical results as the numerical solution of the full-field Navier-Stokes equations. The high efficiency of this model makes it useful to investigate the overall quantitative behavior of related physical-chemical phenomena. In this paper two important properties of hypersonic flow, shock stand-off distance and oxygen disso-ciation, are studied using the quasi-one-dimensional model with the ideal dissociating gas model. It is found that the shock stand-off distance is affected by both chemical and thermal non-equilibrium. The shock stand-off distance will increase when the flow conditions are changed from equilibrium to non-equilibrium, because the average density of the shock-compressed gas will decrease as a result of the increase in translational energy. For oxygen dissociation, the maximum value of its dis-sociation degree along the stagnation line varies with the flight altitude. It is increased at first and decreased thereafter with the altitude, which is due to the combination effect of the equilibrium shift and chemical non-equilibrium relaxation. The overall variation of the maximum dissociation is then plotted in the speed and altitude coordinates as a reference for engineering application.     

变比热平行混合近似解

平行混合变比热改进解析解，合理、简明、快捷，定比热解误差大。在准确的改进解与定比热解之间，探索了几种简化的近似解，从中看出，冷、热、混三气各简化按其自身静温比热比ｋ＿Ｔ作计算的近似解，结果总的误差较小，可供参照取用。算例包括不等总压及等总压混合，并对等总压混合后总压降低作出初步分析。     

助推-滑翔式弹道中段的近似解

研究了助推-滑翔式弹道中段的近似解。首先根据不同假设条件,给出以下四种情况的近似解:(1)自由飞行段的解;(2)只考虑气动力常迎角飞行段的近似解;(3)考虑引力和气动力的合成解;(4)平衡滑翔解。再由上述近似解合成助推-滑翔弹道中段的分段近似解。通过与数值解的对比,验证了近似解的合理性。研究结果可为新型助推-滑翔式飞行器在方案论证阶段时的弹道特性分析以及射程、峰值热流和过载等特征参数的估算提供依据。     

高阶面元法计算多段翼型之位流

本文给出了一种求解多段翼型位流的高阶面元法。面元和奇点的分布分别采用了二阶和三阶的形式。由于高阶面元法较为精确地反映了绕流体的几何形状和奇点表面分布,本文给出了比低阶面元法更为精确的解,并且消除了低阶面元法计算中经常出现的翼型后缘处解的波动。