高速运动弹丸诱导的瞬态流场数值研究

采用基于任意拉格朗日欧拉(ALE)方程的动网格及嵌入网格法,结合二阶精度Roe格式和沉浸边界法,对弹丸由高压气体驱动从静止状态加速至超声速,射出膛口到完全飞离初始流场的整个过程进行了数值模拟.根据数值结果,详细讨论了初始流场、火药燃气流场的形成与发展以及与弹丸的耦合和相互作用过程.      

柱面气相散心爆轰波胞格演化规律研究

本文采用二阶精度NND格式,应用改进的二阶段化学反应模型,通过求解二维Euler方程对柱面气相散心爆轰波胞格演化过程进行了数值模拟.计算结果表明在传播过程中,空间尺度的扩张导致了散心爆轰波后气流的自然膨胀,使得多波结构的爆轰阵面呈现出显著的胞格自组织特性.根据计算结果与理论分析,本文归纳了五种胞格演化模式,分别命名为内凹波阵面会聚、波阵面扭结、褶皱波面失稳、胞格自合并和三波点滑移,并定义了各种模式的物理特征,分析了其相关的演化机制和规律.     

基于二阶滑模的导弹控制系统设计

过载控制是近年来导弹控制律设计领域出现的一种新的非线性控制方法,引起越来越多的重视,但是在采用过载控制方法对侧滑转弯(STT, Skid-To-Turn)超音速巡航导弹纵向回路控制系统设计过程中,由于过载输出与舵偏输入之间存在非最小相位现象,导致系统零动态不稳定.提出用输出重定义将系统状态进行转换的思路,利用内部状态在平衡态时应具有的响应作为其期望参考轨迹,将输出跟踪转化为状态跟踪,以保证系统的零动态稳定.采用二阶滑模设计控制律,消除了一阶滑模所固有的颤振现象,用Lyapunov函数确定出控制器参数,保证了对系统输出的跟踪.鉴于攻角不易精确测量而导致了系统状态不完全可测,控制律设计过程中设计了攻角观测器解决了上述问题.系统仿真研究表明,根据上述方法设计的系统具有比较理想的稳定性和鲁棒性,同时具有更好的动态性能.     

基于GPU和隐式格式的CFD并行计算方法

从图形处理器(GPU)架构特点出发,提出了基于数据并行的隐式计算流体力学(CFD)求解方法,空间离散格式采用迎风Roe格式,计算网格适用于结构和非结构网格。采用统一计算设备架构(CUDA)技术实现了GPU上的隐式CFD并行计算。分别在Intel Core2 Quad3.0GHzCPU和NVIDIAGTX280 GPU上进行了计算,结果表明隐式格式计算速度是显式格式6倍以上,采用显式格式的计算加速比达到28倍,采用隐式格式计算加速比达到了28.7倍,同时计算加速比随计算规模的增加而增加。计算结果和实验结果较为吻合。     

带谐波抑制的新型异步电机无功补偿控制策略

针对由电容器和静态无功发生器组合而成的异步电机无功补偿装置仍然存在的谐波放大的问题,研究了相应的控制策略。首先,研究了基于多个二阶广义积分器的负载谐波电流准确提取算法;其次,在传统的d轴电网电压定向矢量控制的基础上,加入谐振控制器来控制脉宽调制并网变换器向电网输出的谐波电流,从而实现对电网谐波电流的抑制。构建系统的仿真模型,进行了仿真研究,仿真结果验证了谐波提取方法和控制策略的有效性。     

基于自适应二阶终端滑模的SRM直接转矩控制

传统的直接转矩控制（DTC）可以在一定程度上减小开关磁阻电机（SRM）的转矩脉动,但是在换相区的转矩脉动抑制效果较差,并且传统的PI控制存在超调量大、鲁棒性差、抗负载扰动能力有限的缺点。因此,设计了基于负载转矩变化率自适应估计的二阶终端滑模转速控制器,在有限时间内使得系统状态迅速收敛到平衡点,实现无需已知扰动上界快速输出给定转矩。此外,对传统的扇区规则进行改进以优化电压矢量的选取,减少由于换相产生的转矩脉动。仿真和实验结果表明:改进后的DTC系统拥有良好的调速控制性能,抗干扰能力强,鲁棒性好,能进一步减少转矩脉动。      

跨声速操纵面嗡鸣的数值研究

我国的一些高速歼击机，在飞行中曾多次出现操纵面嗡鸣问题，它严重影响了飞机的性能和飞行安全，本文采用数值计算方法，对此问题做较深入的分析和研究。气动力计算采用Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，并采用ＬＵ－ＴＶＤ混合差分格式求解由此气动力解操纵面的弹性振动方程，求出现嗡鸣的Ｍ数及振动的频率、振幅等参数。计算中采用Ｂ－Ｌ代数湍流模型，网格生成采用代数方法及保角变换方法，文中讨论了隐式算法在非定常气动力计算     

非等间距显式算法在重力波非线性传播数值模拟中的初步应用

采用3阶TVD型显式时间积分方法,建立了2维可压缩大气中重力波非线性传播的非等间距显式数值模式.对小振幅重力波传播过程的模拟结果表明,该模式能够很好地再现小振幅重力波的传播过程,并能保持能量守恒关系和各扰动量之间的相位关系,与线性重力波理论预测的结果吻合很好.对有限振幅重力波的模拟结果表明,与格点数相同的等间距网格中的模拟结果相比,采用非等间距网格能够以更高的分辨率模拟重力波的不稳定传播直到破碎的整个过程;在非等间距网格中的模拟结果与在加密一倍的等间距网格中的模拟结果相同,说明了用非等间距网格模拟重力波的饱和与破碎过程是正确的,并且采用非等间距网格能够大大减少计算量.      

一类新型自适应反扩散近似Riemann求解器及其应用

对于包含激波、剪切层等复杂结构的流动问题，为了精确模拟剪切层等精细结构，且保证激波计算的稳定性，必须采用低耗散且强鲁棒的数值通量方法。传统的HLL近似Riemann求解器的耗散性较大，Roe、HLLEM和HLLC等近似Riemann求解器在计算某些含有强激波的物理问题时会出现非物理解，容易导致不稳定。针对这一问题，本文在Riemann求解器中通过合理设计反扩散矩阵，发展了一类具有自适应反扩散的新型Riemann求解器，并将其应用到高阶加权紧致格式，实现了高阶精度求解。通过典型数值算例验证了新型方法的计算精度和稳定性，结果表明本文提出的新型自适应反扩散Riemann求解器克服了传统Riemann求解器的缺陷，既能准确识别剪切层等精细结构，又能保证激波解的稳定性。     

计及横向剪切效应的环形谐振子理论研究CSCD

经典环形谐振子理论是一种基于忽略横向剪切变形假设建立的环理论,因而无法准确求解高径比较大的环问题。为了克服上述不足,基于Timoshenko理论和哈密顿原理,建立考虑横向剪切效应的环形谐振子新理论,包括广义本构关系、平衡方程和周期性条件等。利用经典的布勃诺夫-伽辽金法,并结合有限元法,获得谐振子在缓慢、匀速转动过程中二阶弯曲角频率和进动系数的理论解。研究表明,横向剪切变形会在一定程度上影响环形结构二阶弯曲角频率和进动系数的计算结果;环形谐振子的进动系数不是恒定的,它会随着高径比h/r的增大而缓慢地减小。