测量三维湍流场的热线响应方程及热线校正方程

从Ｊｏｒｇｅｎｓｅｎ等效冷却速度公式出发推出了高阶精度的热线响应方程,它适用于密度不变低速条件下求解湍流场中时均速度和雷诺应力张量。应用该热线响应方程时不严格要求热线探头轴线与流场速度矢量完全一致,因此可简化热线测量的支撑和传动机构。导出的显式热线校正方程可保证校正精度。     

Interacting acceleration compensation algorithm for trackingmaneuvering targets

The two-stage Kalman estimator has been studied for state estimation in the presence of random bias and applied to the tracking of maneuvering targets by treating the target acceleration as a bias vector. Since the target acceleration is considered a bias, the first stage contains a constant velocity motion model and estimates the target position and velocity, while the second stage estimates the target acceleration when a maneuver is detected, the acceleration estimate is used to correct the estimates of the first stage. The interacting acceleration compensation (IAC) algorithm is proposed to overcome the requirement of explicit maneuver detection of the two-stage estimator. The IAC algorithm is viewed as a two-stage estimator having two acceleration models: the zero acceleration of the constant velocity model and a constant acceleration model. The interacting multiple model (IMM) algorithm is used to compute the acceleration estimates that compensate the estimate of the constant velocity filter. Simulation results indicate the tracking performance of the IAC algorithm approaches that of a comparative IMM algorithm while requiring approximately 50% of the computations     

喷丸强化残余应力场三维数值分析

建立了喷丸强化三维有限元模型,实现了喷丸强化处理残余应力场的数值模拟.根据显式动力分析各种能量变化过程,确定了显式分析求解时间选择方法;研究了喷丸速度、弹丸直径、入射角度及摩擦力等因素对喷丸残余应力场影响的一般规律,并与现有相关文献试验结论进行对比;分析了单次和多次冲击下材料内部的残余应力场分布及塑性应变分布特点;建立了不同数目的多丸粒喷丸强化三维有限元模型,研究了不同覆盖率对残余应力场分布的影响规律.      

双层模型机匣和叶片的包容性问题数值模拟

采用显式动力分析有限元技术,研究了飞断叶片模型与机匣模型的撞击过程,并且计算了不同类型机匣模型包容情况下的临界转速。计算结果表明,双层机匣模型的包容能力明显地优于单层机匣模型,间隙情况对双层机匣模型的包容能力有一定的影响。研究结果对多层航空发动机机匣的包容性设计有一定的参考价值。     

爆震波点火器管路传播过程数值模拟

为了进一步研究液体火箭发动机爆震波点火器的工作过程，结合实际试验发动机用爆震波点火器的工作条件，在混合比为3，初始氢氧混气压力为0．1MPa情况下，采用常用的两步显式MacCormak差分格式和有限速率化学反应动力模型，对爆震波点火器管路中爆震波的传播过程进行了数值模拟，研究了爆震波形成后在一维管路中传播特性及性能参数。从数值模拟结果来看，稳定传播的爆震波温度可达3500K，压力达到2MPa，其传播速度可达3430m／s，从爆震波形成到其传播到2m距离的出口，用时不到0．6ms。并与试验结果进行了对比分析，数值模拟结果和试验结果符合很好，为爆震波点火器的工程设计提供了理论依据。     

湍流模型对三维翼梢涡流场数值模拟的影响

 梢涡流场数值模拟对减少梢涡的不良影响意义深远。应用雷诺平均Navier-Stokes(RANS)方法建立了三维翼梢涡流场数值模型,选用的湍流模型包括两种经典两方程模型和一种代数雷诺应力模型。计算内容包括翼梢涡流场速度分布、压力分布及涡核位置,为了准确模拟梢涡涡核内系统旋转和流线曲率影响,将旋转和曲率修正方法应用于选用的湍流模型中。计算时对翼梢涡涡核处网格进行了加密处理。将计算结果与实验结果对比后发现,修正后的湍流模型较标准湍流模型明显提高了预报精度,与实验结果吻合较好,具有工程应用价值。     

一拍跟踪的直升机显模型跟踪解耦控制系统设计

针对直升机动力学表现出的高阶,强耦合的特点,本文基于一拍跟踪的原理,设计了显模型跟踪解耦控制系统,开发了控制阵,为提高系统动态和稳态性能,对系统参数进行了优化,数字仿真表明,模型跟踪解耦性能良好,同时可以迅速改变显模型,以达到改变操纵特性的目的,文章最后对如何选取采样周期及显模型带宽进行了讨论。     

用数值模拟法研究战斗部侵彻多层间隔靶

利用显示动力学有限元程序ANSYS/LS-DYNA模拟着靶速度为2.0Ma和2.5Ma的战斗部对多层间隔靶板侵彻的过程。结果表明：钛合金抗弹性能良好可作为战斗部的壳体材料;战斗部壳体结构强度满足侵彻需求。     

喷丸强化残余应力场三维数值仿真分析

建立了喷丸强化三维有限元模型,实现了喷丸强化处理残余应力场的数值仿真。根据显式动力分析各种能量变化过程,确定了显式分析求解时间选择方法;研究了喷丸速度、弹丸直径、入射角度及摩擦力等因素对喷丸残余应力场影响的一般规律,并与现有相关文献试验结论进行对比;分析了单次和多次冲击下材料内部的残余应力场分布及塑性应变分布特点;建立了不同数目的多丸粒喷丸强化三维有限元模型,研究不同覆盖率对残余应力场分布的影响规律。数值分析结果表明,喷丸速度、弹丸大小、入射角度对残余应力场分布有显著影响。冲击次数对残余应力幅值影响小,但等效塑性应变增加明显。喷丸覆盖率的增加会显著改变残余应力场,并使残余应力场分布更加均匀,但残余压应力幅值与喷丸覆盖率没有必然相关性。     

一种适合迭代求解的反馈力浸入边界法

对Goldstein提出的反馈力浸入边界法进行了新的思考，改进了其对力源项的计算，拓展了该浸入边界法的使用范围。传统的反馈力浸入边界法在进行力源项的计算时，含有对速度误差的时间积分项，只能用于含时间项的Navier-Stokes （N-S）方程的求解，且在显式时间推进时有严格的时间步长限制。本文改进的方法则直接通过迭代过程中的速度误差求和来计算力源项，避免了时间相关的参数，使其不仅能适合非定常隐式时间推进，还能与不含时间项的定常N-S方程求解方法结合。为了验证改进方法的可靠性，对二维静止圆柱绕流、静止流体中的振荡圆柱、运动椭圆翼以及三维静止圆球的流场进行了计算，计算结果均与文献结果符合较好，表明本文改进的方法是有效的。得出的结论为：可以直接基于迭代次数进行反馈力源项的计算，改进的反馈力浸入边界法不仅可与非定常N-S方程结合，进行隐式求解，还可以与定常N-S方程结合用于定常流动的模拟，可将改进的方法运用到更多的流动问题当中。     

Modelling of Solar Wind, CME Initiation and CME Propagation

Simulations of coronal mass ejections (CMEs) evolving in the interplanetary (IP) space from the Sun up to 1 AU are performed in the framework of ideal magnetohydrodynamics (MHD) by the means of a finite-volume, explicit solver. The aim is to quantify the effect of the background solar wind and of the CME initiation parameters, such as the initial magnetic polarity, on the evolution and on the geo-effectiveness of CMEs. First, three different solar wind models are reconstructed using the same numerical grid and the same numerical scheme. Then, different CME initiation models are considered: Magnetic foot point shearing and magnetic flux emergence. For the fast CME evolution studies, a very simple CME model is considered: A high-density and high-pressure magnetized plasma blob is superposed on a background steady state solar wind model with an initial velocity and launch direction. The simulations show that the initial magnetic polarity substantially affects the IP evolution of the CMEs influencing the propagation velocity, the shape, the trajectory (and thus, the geo-effectiveness).     

超声速有翼飞行器的无粘分区数值模拟

本文采用ＭａｃＣｏｒｍａｃｋ两步显式格式，用空间推进法求解Ｅｕｌｅｒ方程来对有翼飞行器的超声速绕流流场进行数值模拟，应用特征边条件处理场面边界及装配激波的计算，在翼缘及表面斜率不连续处，采用局部解析处理，提高了算法的可靠性，复杂外形的计算采用分区处理技术，网格由代数形式生成，在大攻角的计算中，通过对横流速度型的修正改善了背风面压力分布特性，解决了推进中断问题，同时模拟出背风面涡。     

叶轮机械中间隙流与通道二次流相互作用的数值研究

采用数值解法对一压气机静叶栅中叶顶泄漏流与通道二次流的相互作用进行研究 ,详细揭示了在顶部区通道涡、泄漏涡及分离涡的生成和发展。结果表明 ,由于泄漏流动与通道二次流动的方向正好相反 ,适量的泄漏流动可以在一定程度上起到削弱通道涡强度、改善顶部区流动状况的作用。     

基于标准弹道点的最优制导方法

由于基于需要速度的显式制导方法要确定虚拟目标点,虚拟目标点的确定要对地球扁率和再入阻力的影响进行修正,存在着百米以上的修正误差,为此提出了基于标准弹道点的最优制导方法,即将标准弹道的某点作为目标点进行需要椭圆轨道计算及最优入轨控制。提供了采用弹上迭代制导方法实现最优控制,同时给出了预测关机参数的需要椭圆轨道终端约束校正方法。该制导方法无需进行地球扁率和再入阻力的影响修正,可以减少方法误差,仿真计算证明了这一点。     

航空发动机宽弦风扇叶片鸟撞损伤模型标定

为建立航空发动机风扇叶片抗鸟撞载荷能力的量化预测方法,针对特定的航空发动机宽弦风扇叶片设计,依据发动机在典型工作状态下的吸鸟速度、角度等撞击参数开展叶片鸟撞试验,采用显式动力学数值仿真方法,建立叶片鸟撞试验仿真分析模型,并通过对模型中叶片材料参数、鸟体本构模型参数、鸟与叶片耦合接触参数进行敏感度分析,对模型进行标定.结果表明,标定后的鸟撞分析模型所预测的叶片损伤模式与试验结果一致,预测的损伤位置与试验测量结果误差小于10%.      

二维横向喷流的N-S方程数值解

让一股射流横向注入超音速流产生所需要的干扰效应是一种飞行器飞行控制的常用手段,如航天飞机和机动导弹上的反作用控制系统(RCS).所以,对射流与主流相互干扰的某些流动特性应该有个基本的了解.本文通过应用MacCormack显格式和Baldwin-Lomax修正的代数湍流模型求解二维RANS方程对带横向射流的绕后台阶的M_∞为2.19和6.00的超音速外干扰流场进行了数值计算.与无喷射流的情形相比,回流涡旋区由2个增加到4个,且底部压力有大的回升,喷口宽度越大,回升效应越显著.另外,数值计算表明,在压力梯度变化比较平缓的区域附加一四阶人工粘性能有效地抑制数值振荡,加快收敛.     

直升机起落架抗坠毁性能的有限元仿真评估

 建立了直升机起落架抗坠毁有限元模型,模拟硬着陆和坠毁过程,以考察起落架在坠毁事故发生时的吸能能力。有限元建模使用真实的几何模型,用一个虚拟的框架代替机身把起落架连接起来,全机质量折算为减缩质量并分布到该框上。缓冲器用弹簧—阻尼器单元替代,其参数由性能曲线给出。使用非线性瞬态动力学的显式解法求解冲击问题。对两种情况进行了仿真:分别以6.0m/ s 和10.2m/ s 的垂直速度撞击地面。坠毁实验表明仿真结果与实验结果具有较好的一致性。     

火箭发动机塞式喷管流场的数值研究

用时间相关法和张涵信的无振荡、无自由参数ＮＮＤ显式格式,对火箭发动机塞式喷管二维轴对称定常流场进行了数值模拟。计算得到了有粘和无粘中两种情况下,喷管流场的马赫数压强、温度、流线以及速度矢量分布图,计算网络采用代数方法生成,并通过求解Ｌａｐａｃｅ方程对网格进行修正。计算结果表明：ＮＮＤ格式用于计算塞式喷管流场具有较高的效率和良好的分辨率;无粘和有粘情况得到的流场结构有较大的差别;基础喷管的出口马赫数     

超音速气流中横向喷射氢气流场数值模拟

本文通过应用Mac Cormack显式时间分裂法和Baldwin-Lomax修正的代数湍流模型求解二维可压缩的N-S方程和混气组分方程,对超音速气流中横向喷入氢气后的流场进行了数值模拟。计算结果示出了流场中速度矢量分布,静压、静温和氢气质量分数的等值线分布,反映出了喷嘴附近流场回流区的大小及可能的点火位置,从而为超燃冲压发动机的设计工作提供理论指导。     

一种求解透平叶栅三维流场的高精度TVD格式

本文提出了一种应用三阶精度Chakravarthy-OsherTVD格式求解环形透平叶栅三维定常流场的隐式近似因子分解方法。控制方程为柱座标速度分量表示的任意曲线座标系下的三维欧拉方程。采用Beam-Warming近似分解, 对隐式项进行迎风差分, 得到了时间精度为一阶、空间精度为5点三阶的隐式格式。算例对5点三阶和5点二阶精度的格式进行了比较。表明空间三阶精度隐格式对计算结果没有明显改进, 但收敛残差有较大降低。